книги из ГПНТБ / Козобков, А. А. Электрическое моделирование вибраций трубопроводов
.pdfu0dx\
о
( 1 3 8 )
т
Легко показать, что электрический эквивалент граничных условий, описываемый уравнениями (138), может быть получен путем включения сопротивлений R/ls и R/1l параллельно с ин дуктивностями Lei и LCl в схеме, указанной в четвертой строке табл. 2. При этом величины Ал. и А?/,, находятся из соотношений
(139)
т /И
Таким образом, если коэффициенты h\ и /гг, определяющие трение в опоре, известны, то величины, моделирующие трение Rn1 и Rhi, находятся из соотношений (139) и включаются в мо дель по указанной выше схеме.
Однако при расчете колебаний действующих трубопроводных систем диссипативные характеристики опор (h\ и /г2) не всегда могут быть известны, так как это связано со сложностью их экс периментального определения.
В этом случае можно воспользоваться некоторой общей дис сипативной характеристикой трубопровода, а именно его доброт ностью, которая определяется экспериментально по кривой зату хания свободных колебаний трубопровода. Тогда сопротивление, моделирующее неизвестное трение в опорах, наиболее целесо образно равномерно распределить по звеньям электрической мо дели. При этом величины сопротивлений, включаемых в каждое звено модели, подбираются исходя из равенства добротности объекта и модели на первых резонансных частотах.
Выясним правомерность такой замены и получим зависимо сти между коэффициентами сосредоточенного в опорах и распре деленного по длине трубопровода трения, сравнивая для того и другого случая одни и те же динамические характеристики для одних и тех же точек трубопровода. Выберем в качестве таких характеристик механические импедансы трубопровода (отноше ние силы или момента силы к скорости или, соответственно, к угловой скорости).
Выражение импедансов в любой точке трубопровода через импедансы его концов, а также частотные уравнения легко полу чить из соотношения, связывающего импедансы начала и конца трубопровода. Подобное соотношение для динамических жест костей было получено в работе [14].
Запишем его в импедансной форме:
60
|
ll°ll=P> (—/ЧХ^> |
i / wM^z ;/z o. (-•/4X0)11 X |
|
||||||||||
|
а |
Ь |
Ь |
с |
с |
d |
d |
a |
|
1 |
|
||
|
л |
В |
|
В С |
С D |
D Л |
|
|
|||||
|
л В |
|
В С |
С D\ |
D A |
|
1 |
|
|||||
X |
е |
/ |
|
f |
g |
g |
h\ |
h |
e |
X |
1 |
(140) |
|
е |
/ |
|
/ |
g |
g |
h |
/1 |
e |
|||||
|
|
(/“м)"Zy н2д |
|
||||||||||
|
Ei В |
|
И G |
G H |
H E1 |
|
|||||||
|
Е1 F F G G H\ H E , |
|
1 |
|
|||||||||
|
а |
b |
b |
с |
с |
d |
d |
a |
|
/'“MZn |
|
||
|
0(0) |
|
|
|
|
|
|
|
|
|
|||
где г уи- |
-импеданс левого |
|
|
|
|
||||||||
|
i/(0) |
|
динате; |
|
|
|
|
|
|
|
|||
Z„и» = |
M(0) |
— импеданс левого |
конца трубопровода по углу; |
||||||||||
-.6(0) |
|||||||||||||
Zy = |
|
---- импеданс правого |
конца трубопровода по |
ко- |
|||||||||
|
у(П |
|
ординате; |
|
|
|
|
|
|
||||
Z B= |
X -—----импеданс правого конца трубопровода по углу; |
||||||||||||
|
0(0 |
|
|
|
|
с = —EJk?T(kl); |
|
|
|
||||
a=S{kl)\ b = - k V { k l ) \ |
d = |
- E J k 2U{kl)\ |
|||||||||||
A = |
T(kl) . |
Л=5(А/); |
C= EJkHJ[kl); |
D = EJkV(kl); |
|
||||||||
|
|
; |
|
||||||||||
|
y m |
|
. |
/ |
= ^ |
; g= S (& ); |
A = |
r(W) |
|
||||
|
' |
ЛЗEJ |
' |
|
|
№E! |
|
|
|
|
|
||
£ 1 = |
U(kl) |
|
F _ |
Т Щ ) . |
G= kV (kly, H — S[kl), |
|
|||||||
|
|
№EJ |
|
|
kEJ |
’ |
|
|
|
|
|||
где 5, Г, |
(У, К — функции Крылова, k = k(aM). |
|
|
Сначала рассмотрим задачу о изгибных колебаниях трубо провода с одной шероховатой опорой. Будем считать, что трение подчинено линейному закону — момент трения пропорционален первой-степени угловой скорости. Вторую опору для простоты анализа возьмем шарнирной. Граничные условия в этом случае будут
1/(0) =0; у{1) =0;
|
(141) |
М (0) =0; |
М (/) =Ш |
или в импедансной форме |
|
Z y a — 00, Z y — со ; |
|
Z 0h= °; |
( 1 4 2 ) |
Z 0= yo)M/z, |
где h — коэффициент сопротивления.
Л
61
Решая совместно уравнения (140) и (141), получим частотноерешение задачи
I/2 (V ) - |
Т2 (К1) + |
[Г (Лс/) U (V ) - 5 (£C/)IW )] = °- (143) |
|
|
k\EJ |
|
|
Для того чтобы выполнялось соотношение (143), |
необходимо |
||
|
V 2{k j . ) - T 2{kzl) = 0; |
(144) |
|
|
T{kcl)U (kJ .)-S(kzl)V{k,l)=0. |
(145) |
|
Как известно, уравнение (144) справедливо для |
граничных |
||
условий, определяющих шарнирное закрепление, |
уравнение |
||
(145) — для |
случая, когда |
один конец трубопровода |
закреплен, |
шарнирно, а другой — жестко. Из этого следует, что при приня тых нами граничных условиях возможны свободные колебания с двумя группами частот:
—■первая группа частот определяется уравнением
корни которого |
sin |
kcl= Q, |
(146) |
|
|
|
|
*c= |
- f |
, (s=l,2,....); |
|
— вторая группа частот определяется уравнением |
|||
tg k cl = thkcl, |
(147) |
||
его наименьшие корни |
/ес2/ = 7,07; /гс8/= |
10,21. |
|
kcil= 3,93; |
Этим группам частот соответствуют две группы собственных форм. В дальнейшем нам понадобится форма, соответствующая, частотному уравнению (144). Она может быть записана как
ус{х)= s i n - ^ - . |
(148) |
Найденные по уравнениям (146) и (147) |
частоты являются соб |
ственными. |
|
Получим выражения собственных частот и форм трубопрово да в случае, когда трение распределено по его длине. Запишем, силу трения в виде, наиболее удобном для моделирования
F TP=2b |
д5</ |
|
дх*dt |
где 2b — коэффициент сопротивления. |
|
Рассмотрим наиболее простой |
случай, выбирая трубопровод |
сшарнирными опорами.
Вимпедансной форме граничные условия имеют вид
2 ^ = 0 , Z 0—0, ZyH^=co, Zу—.со. |
( 149) |
62
Подставляя граничные условия (149) в (140), получим частот ное уравнение
V2{kpl)—T2(kpl)=Q. (150)
Форма колебаний может быть записана как
k0x |
, _ |
z/„=sin — . |
(151) |
Можно показать, что если сопротивление имеет принятый нами вид, то
ил НОым(£ / — 72К .) .
р(£/)2 +
|
н Ч |
« ч |
1_ |
к \ |
2 |
||
EJ |
( £ / ) 2 |
|
При Ь —0 получим
Н-0мм
£7
Вдальнейшем будем считать, что
«ч
(£7)2 « 1 .
т. е. влиянием сопротивления на собственные частоты прене брегаем.
Сравнение формул (148) и (151) показывает, что отрезок трубопровода с одной шероховатой опорой при колебаниях на частотах из группы (146) имеет такую же собственную форму, как и отрезок трубопровода, в котором трение распределено.
Одинаковые формы выражений (148) и (151) позволяют сравнивать одни и те же точки трубопровода для определения связи между коэффициентами сосредоточенного и распределен ного трения, причем этой точкой может служить любая точка, принадлежащая интервалу (0, /). В качестве такой точки выбе рем точку с абсциссой Z/2.
Из выражения (140) найдем импедансы для сосредоточен ного и распределенного трения Zoy(i/2) и Zpy(;/2). Для этого необ ходимо в выражение (140) подставить Z0(i/2) и Zy^/2) вместо Z с„ и Zyu, а вместо Z„ и Zj —их значения из выражения (142) и, соот ветственно, из выражения (149). Тогда, считая b малой величи
ной, т. е. пренебрегая в искомых |
выражениях членами, содер |
жащими b в степени два и выше, получим выражения |
|
|^Ci(Z/2) I И |
\Z Jy (l / 2 ) \- |
63
Приравнивая их, найдем связь между коэффициентами рас пределенного и сосредоточенного трения
1 + 28,38' |
А |
|
0,96 1 + 4 8 ,6 |
62 |
1 — 607 - |
62 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
EJy^l*I |
(152) |
||||
|
|
|
|
|
|
|
|
62 |
|
|
|
|
|||
+ |
15,43 |
|
|
|
|
1 — 1216 |
|
|
|
|
|
|
|||
/3Но |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, мы установили правомерность замены сосре |
|||||||||||||||
доточенного в опоре трения h = М/д распределенным |
и |
опреде |
|||||||||||||
|
|
|
|
|
|
лили связь |
между |
коэффици |
|||||||
|
|
|
|
|
|
ентами h u b . |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
Для |
|
экспериментальной |
||||||
|
|
|
|
|
|
оценки |
погрешностей, |
вызыва |
|||||||
|
|
|
|
|
|
емых подобной заменой, на |
|||||||||
|
|
|
|
|
|
электрической |
модели |
первого |
|||||||
|
|
|
|
|
|
приближения |
однородного от |
||||||||
|
|
|
|
|
|
резка трубопровода |
были |
оп |
|||||||
|
|
|
|
|
|
ределены три первые формы и |
|||||||||
|
|
|
|
|
|
частотные |
характеристики |
од |
|||||||
|
|
|
|
|
|
нородного участка трубопрово |
|||||||||
|
|
|
|
|
|
да |
при |
различных |
граничных |
||||||
|
|
|
|
|
|
условиях и при сопротивлении, |
|||||||||
|
|
|
|
|
|
сосредоточенном |
в |
опорах |
и |
||||||
|
|
|
|
|
|
распределенном |
по длине тру |
||||||||
|
|
|
|
|
|
бопровода. |
частотных |
характе |
|||||||
|
|
|
|
|
|
|
Снятие |
||||||||
|
|
|
|
|
|
ристик |
и |
форм |
колебаний на |
||||||
|
|
|
|
|
|
модели |
производилось |
|
по ме |
||||||
|
|
|
|
|
|
тодике, |
изложенной |
|
в преды |
||||||
|
|
|
|
|
|
дущем параграфе. |
|
|
|
од |
|||||
Рис. 29. Три собственные |
формы |
|
Электрическая модель |
||||||||||||
нородного отрезка трубопрово |
|||||||||||||||
участка трубопровода с распреде |
да с учетом внутреннего сопро |
||||||||||||||
ленным трением (точки) |
и |
тре |
|||||||||||||
нием, сосредоточенным в |
опорах |
тивления состояла из |
1 1 |
звень |
|||||||||||
|
(кружки): |
|
|
|
ев с параметрами |
|
|
|
|
|
|||||
а—первая |
форма; |
б—вторая |
форма; |
Ci = 10- 9 |
Ф; L,= 1,5-10- |
3 Г; |
|||||||||
|
в—третья |
форма |
|
|
|||||||||||
|
|
|
|
|
|
|
£i = 4; У?| = 0,5 |
Ом. |
|
||||||
На рис. 29, а, б и в и рис. 30 |
приведены в безразмерном виде |
||||||||||||||
три первые собственные |
формы и |
частотная |
характеристика |
||||||||||||
(указаны |
кружками) |
трубопровода |
с |
шарнирными |
|
опорами |
|||||||||
с сосредоточенным |
в них |
трением. |
Выбранному |
трубопроводу |
|||||||||||
с трением в опорах |
(рис. |
31, а) |
соответствует электрическая мо |
дель, показанная на рис. 31,6. Пунктирная кривая на рис. 30 соответствует частотной характеристике трубопровода с учетом трения только внутри металла.
64
Частотная характеристика трубопровода с трением в опорах (обозначена кружками на рис. 30) позволила вычислить его добротность
|
|
|
|
Q |
|
/ 1 |
5. |
|
|
|
|
|
|
|
|
|
|
2 |
Д/ |
|
|
|
|
|
|
Следовательно, |
величина |
|
полного |
сопротивления |
одного |
|||||||
звена будет |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ ? = |
^ |
= 4,5 |
Ом, |
|
|
|
||
|
|
|
|
Л, |
|
|
|
|
|
|
|
|
где юэ — первая резонансная частота модели. |
|
|
|
|||||||||
7,0 |
т |
|
|
II |
|
|
|
|
Л |
|
|
|
И3 |
|
|
|| |
|
|
|
|
II |
|
|
|
|
"Н" |
|
И |
|
|
|
|
и |
|
|
|
||
Ц8 |
и |
|
|
|
|
|
|
11 |
|
|
|
|
0,1 |
и |
|
|
тг |
|
|
|
|
И |
|
|
|
и |
|
|
]Г |
|
|
|
|
II |
|
|
|
|
0,6 |
It1 |
|
|
|
|
|
II |
|
|
|
||
0,5 |
■| |
|
|
|| |
|
|
|
|
11 |
|
|
|
1 |
|
|
п |
|
|
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
0,4 |
н- |
|
|
г* |
|
|
|
|
д |
|
|
|
|
1\ |
|
|
1 |
|
|
|
|
|
|
|
|
0,3 |
1 . |
|
!\Г |
|
|
|
|
1\ |
|
|
|
|
О,г |
~ 1 |
V |
^7 |
Л |
) |
tsk |
|
Л |
1 |
|
|
|
0,1 |
Г - |
|
|
|
Г |
|
|
|
||||
О |
1 I f , 3 |
4 |
5 6 7 f j |
9 |
10 |
11 12 |
13 |
n f J 5 |
16 |
17 18 19/( к Г ц ) |
|
|
|
'/ |
|
|
Z |
|
|
|
|
J |
|
|
|
Рис. 30. Частотная характеристика |
участка |
трубопровода |
|
|||||||||
с распределенным |
трением |
(точки) |
и трением, |
сосредото |
|
|||||||
|
|
|
ченным в опорах (кружки) |
|
|
|
||||||
Электрическая модель трубопровода с распределенным по |
||||||||||||
длине сопротивлением показана на рис. 31,6. На рис. 29, |
а , б , в |
|||||||||||
и рис. 30 приведены |
частотные |
характеристики и собственные |
||||||||||
формы трубопровода |
(указаны |
точками) |
с |
распределенным |
||||||||
трением. |
|
экспериментальные |
данные |
показывают, что |
||||||||
Приведенные |
при моделировании возможна замена трубопровода с сосредото ченным трением трубопроводом с распределенным трением. При этом амплитудная погрешность и погрешность моделирования по форме не превысит 10—15%. Исследование подобной замены на электрической модели трубопровода со смешанными и жест кими опорами на концах трубопровода и с промежуточными опорами с трением также показали, что погрешность моделиро вания по форме не превысит 10—15%.
Следует отметить, что заменой сосредоточенного трения рас пределенным следует пользоваться только в случае, если не мо гут быть заданы характеристики трения опор, так как этот метод не исключает амплитудной погрешности и погрешности формы.
3 |
366 |
65 |
Выше мы предполагали, что трение, сосредоточенное в опоре, пропорционально первой степени угловой скорости. Это предпо ложение справедливо в тех случаях, когда ширина опоры доста точно велика. В том случае, когда ширина опоры незначительна, и в опоре наблюдается проскальзывание, определяющим будет трение, вызываемое продольными колебаниями трубопровода. Рассмотрим это трение, определим, как оно отражается на попе-
0 — |
4 ,5 Ом 4,50 м 4 ,5 0 м .— 0 |
В)
Рис. 31. Схема модели трубопровода с трением:
а—трение сосредоточено; б—трение распределено
речных колебаниях трубопровода и сравним его с распределен ным трением.
Пусть трубопровод закреплен как показано на рис. 32. В опо рах возникает сила трения /чр, которая вызовет дополнитель ную нагрузку на трубопровод дтр. Уравнение свободных колеба ний можно записать в виде
h> |
й - у |
£ У 0 |
+ ,„ = О . |
(153) |
|
d t °-
Найдем нагрузку <7 тр, предполагая силу трения в виде
•ЕтР = аД, |
(154) |
где а — коэффициент трения, А — скорость относительного удли нения трубопровода.
' Тр
4 -
э-
щг,
Рис. 32. К расчету декремента колебаний
Величина А может быть вычислена как
|
|
|
й 5 |
|/^ +U9 d x ~ i |
|
|
|
д = - |
д |
I |
(155) |
|
|
|
d t |
|
|
Если |
д у |
“С 1> |
что допустимо в рассматриваемом случае, |
||
|
то рыражеиие Ц55). можно переписать в виде
66
|
д |
i |
(156) |
|
д |
dx . |
|||
dt |
||||
|
|
|
Тогда, учитывая соотношение (156), можно записать выражение для qтр:
?-гр |
д2у |
_ д _ |
|
|
|
(157) |
дх2 |
dt |
|
|
|
||
|
|
|
|
|||
Подставим выражение |
(157) |
для силы |
qTp в |
уравнение |
(153) |
|
и получим уравнение |
для |
свободных |
колебаний трубопровода |
|||
в окончательном виде |
|
|
|
|
|
|
<9x4 1 дх2 |
— — |
Г ( - ^ - )2dx-\-\>-0^ - = 0 . |
(158) |
|||
d t |
2 |
J \ d x } |
1 0 |
dt* |
|
Запишем граничные условия в следующей форме:
0( О ) = |
£ ( / ) = О ; |
|
|
= |
0 . |
= 0 . |
|
(159) |
|
|
|
|
|
0ЛГ2 .Г-О |
дх2\х=1 |
|
|
||
Решение уравнения |
(159) |
будем, |
как и ранее, искать в форме, |
||||||
предложенной Фурье |
|
|
со |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
У = |
У s i |
|
|
|
|
||
где ys=Xs(x)Ts(t). |
|
|
|
j=i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Считая нелинейность в уравнении |
(158) |
слабой и не влияю |
|||||||
щей на собственные значения уравнения |
(158) |
с граничными |
|||||||
условиями (159), придем к решению уравнения |
(126),, |
которое |
|||||||
для первой моды запишем как |
|
|
|
|
_ . |
||||
|
y = T(t) sinkx. |
|
|
|
(160) |
||||
Подставим решение (160) |
в уравнение (158), получим1 |
|
|||||||
|
|
I |
|
|
|
|
|
- = 0. |
(161) |
Р Е У Т - Р - ^ - Т |
j |
-^-(Г cos k x f dx-\- y-0f |
|||||||
Выполним в уравнении (161) |
квадратуру |
|
|
|
|||||
i |
|
|
|
|
|
|
|
|
(162) |
— (Т cos k x f d x = — T%— |
— sin 2kx |
||||||||
dt v |
' |
|
|
dt |
2 |
1 |
2k |
|
|
Подставляя пределы |
интегрирования в выражение |
(162) и, |
|||||||
учитывая, что sin 2 £ /= 0 , получим |
|
|
|
|
|
3* |
67 |
|
l |
|
|
|
Г — (T cos k x f dx = — 2TT. |
(163) |
|
|
.) dt |
2 |
|
|
0 |
|
|
Подставим значение интеграла (163) в уравнение (161) |
|
||
|
k^EJT — Id -Щ- ТгТ -|- \>-0Т = 0 |
(164) |
|
д |
akU |
№Е1 |
|
и обозначим |
а —----- ; |
г = ------- . |
|
|
2fj<) |
Ди |
|
Тогда уравнение для определения Т (t) будет иметь вид |
|
||
|
f + аТ*Г + И1Т = 0. |
(165) |
|
Решение нелинейного уравнения (165) будем искать в виде |
|
||
|
Г = |
А sin [o(t)\ е - Е<'\ |
|
где А — начальная амплитуда. |
|
||
Подставим это решение в уравнение (165) |
|
||
— (ср) 2 |
А е_Е sin cp-j-cpA е_Е cos ср — 2гееА е_Е cos ср — |
|
—А е_Е [—( е )2 -(- е]sin tр -f-а.А3 е~ЗЕср sin2 ср cos ср —
—аА3г е—3еsin3 со— А е—Е[—s2-)- s] sin ср -j- £2А е_Е sin ср = 0.
(166)
В уравнении (166), как и раньше,предполагая
siп3 ср; — sines; |
cos^cp^ — cos ср, |
4 |
4 |
т. е. пренебрегая высшими гармониками процесса, сгруппируем члены с одноименными тригонометрическими функциями
(?)2 — яА2 ее_2 Е-\-Ь2 — (s) 2 + ej sin cp-J— |
|
-f j^cp—2ecs-1—^-aA2e-2Ecpj coscp= 0. |
(167) |
Приравнивая выражения в квадратных скобках в уравнении
(167)нулю, получим два уравнения для определения функций
Фи е
— (ер)2 — - i- a А2 е -2б£ + & 2 — (е) 2 + Е= 0; |
(168) |
. tp'_2Ecp + i a A 2 e -2 EcP= 0 . |
(169) |
68
Предположив затухание слабым, т. е. пренебрегая (е) 2 |
и е в урав |
||||||||||||
нении (168), получим |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
_ (ср) 2 |
---—а Л2 e -2Es-f Ь -= О, |
|
(170) |
||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
аА^е—2s |
|
)= |
|
, , |
/ |
Г . |
1 |
Л2 |
2 |
■ь 1 |
I |
|
|||
|
о 1 |
|
1 |
------ |
а ---- |
- е . |
8 |
b2 |
)■ |
||||
|
|
У |
|
|
4 |
62 |
|
||||||
Ограничиваясь нулевым приближением |
в |
определении частоты |
|||||||||||
колебаний, т. е. считая, что ср = 6 , |
|
|
|
|
|
||||||||
получим |
|
|
|
|
|
|
f |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
___ |
<?= f bdt=wa0t-]r <o0, |
|
|
|
||||||
|
|
|
|
6 |
|
|
|
|
|
|
|||
где о)МО" |
\ |
f |
kAE-~ ; tp0 |
—начальная фаза |
колебаний. |
|
|||||||
|
I |
|
|
но |
|
|
|
|
|
|
|
|
|
Подставим нулевое приближение для ф в уравнение (169): |
|||||||||||||
|
|
|
|
|
|
_ 2 =6 + — Л2 е -2^ = 0 , |
|
|
|
||||
|
|
|
|
|
|
|
1 |
4 |
|
|
|
|
|
откуда е = |
|
1 |
|
|
|
е~-26. |
|
|
|
|
|
|
1171) |
|
— аЛ2 |
|
|
|
|
|
|
||||||
|
|
8 |
|
|
решим, разделяя переменные |
|
|
||||||
Уравнение (171) |
|
|
откуда Ле~Е
\ / ~ at + Л2
Следовательно,
Т = A sin («м0^ + 9 с
V a t + #
Подберем такую функцию Ле~6(, чтобы функция
V
отличалась от нее на небольшую величину для любого значения t. Для этого запишем
|
Ле—5/. |
4 |
|
|
|
|
a t+ A2 |
|
|
ИЛИ |
л— 5/— . |
1 |
(172) |
|
|
||||
1 f |
«Л2 |
|||
|
|
|||
|
V |
1 + — ' |
|
69