Файл: Козобков, А. А. Электрическое моделирование вибраций трубопроводов.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), получим
Н = ^ 2■ f [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