Файл: Козобков, А. А. Электрическое моделирование вибраций трубопроводов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.10.2024

Просмотров: 80

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Опыт эксплуатации трубопроводных, гидравлических, газо­ вых и воздушных систем летательных аппаратов, машин и раз­ личных устройств показывает, что наибольший процент всех встречающихся разрушений трубопроводов составляют уста­ лостные разрушения, происходящие под действием вибраций [3]. В связи с этим борьба с вибрациями трубопроводов приобрета­ ет важное народнохозяйственное значение.

Динамические режимы по своей природе сложнее статиче­ ских, так как они характеризуются не только величиной усилия, точкой его приложения и направлением, но законом его измене­ ния во времени. Естественно, что расчет динамики трубопровод­ ных систем оказывается более сложным. Основной целью такого расчета является определение параметров системы, амплитуда вибрации которой не превышала бы допустимую величину.

Трудность задач, связанных с проектированием трубопровод­ ных систем, испытывающих воздействие пульсирующего потока, обусловлена тем, что теоретические основы их решения бази­ руются на выводах комплекса закономерностей из различных областей науки.

При проектировании в большинстве случаев удается добиться создания трубопроводной системы с допустимым уровнем виб­ рации. Однако при этом на отдельных элементах трубопровод­ ной системы уровень вибрации может достичь опасной величины.

Таким образом, существующие методы проектирования тру­ бопроводных систем не исключают ошибок в их построении и: размещении различных элементов, например, технологических аппаратов и запорной арматуры. Степень влияния тех или иных ошибок разрабатываемого проекта различна, но все они в боль­ шей или меньшей степени приводят к необходимости «доводки» объекта до работоспособного состояния в процессе его освоения, т. е. к переделке трубопроводной системы.

В общем случае решение задачи о вибрации трубопроводов необходимо искать в исследовании двух систем: газогидравличе­ ской и механической. Причем, воспользовавшись методом газогидродинамики, необходимо получить уравнения для расчета характеристик пульсирующего потока, определить колебания давления и скорости в функции времени для любого сечения трубопровода и найти значения собственных частот. Затем необ­ ходимо найти величины усилий, возникающих в элементах тру­ бопровода под действием пульсирующего потока, их направления и точки приложения и, воспользовавшись методами динамики сооружений, определить величины, характеризующие трубопро­ водную систему как колеблющийся объект, т. е. определить амплитуды колебаний трубопроводов и аппаратов, формы коле­ баний и значения собственных частот.

Для отыскания наилучшего проектного решения такой расчет необходимо выполнить многократно, изменяя схему расположе-

8


ния элементов трубопроводной системы, арматуры и опорных устройств.

Однако решение дифференциальных уравнений, описываю­ щих пульсирующий поток жидкости или газа и колебания трубо­ проводных систем, при помощи известных аналитических мето­ дов в большинстве случаев оказывается затруднительным из-за сложности трубопроводных систем и сложного характера взаи­ модействия потока жидкости или газа и трубопровода.

Сложность задачи и необходимость выполнения при ее реше­ нии большого количества вычислительных операций вызывает естественное стремление проектировщиков использовать вычис­ лительную технику.

В книге сформулированы основные требования, которым должно удовлетворять вычислительное устройство, предназна­ ченное для расчета параметров вибрации трубопроводных систем. Анализ этих требований, а также характеристик различных ви­ дов вычислительных устройств показывает, что наиболее мощ­ ным и удобным средством для достижения поставленной цели является электрическое моделирование, которое за последние годы получило большое развитие как в нашей стране, так и за рубежом, и которое, как показывает опыт, успешно может быть использовано при проектировании трубопроводных систем и ис­ следовании их динамики. Следует отметить, что особенно боль­ шой вклад в развитие методов электрического моделирования сделан советскими учеными и, ib первую очередь, проф. И. М. Тетельбаумом.


Глава 1

СРЕДСТВА РАСЧЕТА ДИНАМИКИ ТРУБОПРОВОДНЫХ СИСТЕМ

Расчет динамики трубопроводных систем может быть выпол­ нен при помощи различных средств и методов. Прежде чем оце­ нивать их сравнительную эффективность проанализируем основ­ ные уравнения колебаний трубопроводных систем, причем пос­ ледние будем рассматривать как стержневые системы.

§ I. ОСНОВНЫЕ УРАВНЕНИЯ. ГРАНИЧНЫЕ УСЛОВИЯ

Как указывалось выше, трубопроводные системы совершают изгибные, крутильные и продольные колебания. Рассмотрим основные уравнения, описывающие их.

Изгибные колебания

Пусть трубопровод длиной I (рис. 1), С ПОСТОЯННОЙ ПОГОННОЙ

массой р0 и с постоянной жесткостью EJ находится под давле­ нием р. Если при этом скорость протекающего по трубопроводу

потока v, а погонная масса про­ текающего по трубопроводу про­ дукта р, то кинетическая энергия системы составит

У

Рис. 1. Схема трубопровода, подвержен­

Рис.

2. Расчетная схема определе­

ния

нагрузки от сил внутреннего

ного изгибным колебаниям

 

давления

 

 

(1 )

10

а потенциальная энергия скоростного напора

Нагрузку от сил внутреннего давления можно вычислить как

 

 

 

р | cos yds

 

 

 

 

о______

(3)

 

 

 

d x

 

 

 

 

где

ср — координата по сечению трубы (рис. 2);

 

ds — элементарная площадка.

 

 

Из рис. 2 видно, что

 

 

 

 

 

 

ds = -y- zd<o.

 

где

 

 

 

 

 

2 = — (1 — cos w)d ( — ) .

 

 

2

\д х

)

Раскроем выражение (3):

 

 

 

Я-

 

 

 

2~

d^

(1 — cos ср) cos yd

'j dy,

 

qpdx — p J

 

откуда

 

 

 

 

 

 

- p F B L ,

(4 )

 

 

 

дхч-

 

где

ndi

 

сечения трубопровода «на просвет».

F — —-р — площадь

■Знак (—) в выражении (4) означает направление равнодейст­ вующих сил внутреннего давления в сторону выпуклости. Тогда потенциальная энергия от сил внутреннего давления составит

=

р£_ 1дУ \ 2

 

П-2

2

)

(5 )

Изгибные колебания трубопровода вызывают появление про­

дольной силы, величина которой может быть определена как

 

 

H = s.EFт,

(6)

где е — относительное удлинение трубы;

 

FT— площадь поперечного сечения тела трубы. Из рис. 1 следует

S — 1


— длина средней линии изогнутого трубопровода при условии

(<Р~) ^ Решая совместно выражения (6) и (7), получим

Н = ^ 2f [dJ - \ d x .

(8)

21 JО1дх

Учитывая выражения (1), (5), (8) и обозначая статический момент инерции трубопровода J, запишем * действие по Остро­ градскому — Гамильтону:

 

 

/,о

I

 

X

EFr

дц_

' dx-\- рт)2

dx dt.

21

J \дх

 

 

о

 

(9 )

 

 

 

 

Уравнение Эйлера—Лагранжа, минимизирующее функционал (9), получает вид

Ej * y + f o + v .)*y + * l

1™'-¥-,Я57)!" + Х Г +

 

дх4 1 v 0 Г дР ' дх2

 

4-2М.Ц <% = 0.

( 10)

 

дх dt

 

Полагая последний член в уравнении (10)

2[vo - ^ - = 0,

дх dt

т. е. пренебрегая кориолисовыми силами из-за их очевидной ма­ лости, запишем уравнение (10) в виде

EJ d*i , д*у

pv2-\-pF

21 J \дх

dx + (lA+ !io)

0-

(11)

дх$ дх2

 

 

о

 

 

 

Для оценки членов уравнения (11), стоящих в квадратных скобках, решим его. Решение будем искать в виде суммы

( 1 2 )

i=1

*При записи функционала (9) пренебрегаем, как это примято, инерцией вращения сечений.

12


Подставим решение (12) в уравнение (11) и получим

E J ^

^ '

V

U s i n

8 «+Р / W ^ +

-

A"' sin cp; (/)

2X 2

^

f

T f

/=i

 

 

i

 

о

L

 

 

 

 

X

sincp^x

’y A " " s in tp,.+ (p.+ (x0) X

 

 

 

 

 

 

/=1

 

 

 

 

 

 

 

X 2

X ?/ cos <?i —

X

^'-)2 sin ?<■

=

0.

 

(13)

./«1

 

 

/=1

 

 

 

 

 

 

Опуская в выражении

(13) знак 2

и индекс /,

а

также

полагая

 

 

1

 

 

 

3

 

 

 

 

sin3<p= — (3 sin ср— sin 3® )« — sin <р,

 

 

т. е. пренебрегая высшими гармоническими составляющими, перепишем выражение (13), сгруппировав члены, содержащие одноименные тригонометрические функции:

sin c p |£ 7 X v + X ' I*»2 + p F - ^ ^ U x ' f d x

Ы

— А' (ср )2 (|л-|- р.а)) -j- (;л-|- [aq) А'ср cos ср—0.

(14)

Приравнивая в уравнении (14) коэффициенты при

тригономет­

рических функциях нулю, получим

 

 

 

 

(t* + K - o ) X p = 0 ,

 

(15)

E J X xv +АГ"

\> tf+ pF - У р Ц ( x ') 2dx

- * ( с р ) 2 (

^ 0) = О .

 

 

 

 

(16)

Из уравнения

(15)

следует

 

 

cp=0;

ср = | v d t — Const— со; J j

срй^= ш/-|-ср,

где ш — круговая частота; ср — начальная фаза колебаний.

Учитывая последнее обстоятельство, перепишем уравнение (16), полагая в нем нелинейный член равным нулю, т. е.

( X ' f d x = 0.

13