Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Козобков, А. А. Электрическое моделирование вибраций трубопроводов

.pdf
Скачиваний:
7
Добавлен:
19.10.2023
Размер:
6.3 Mб
Скачать

u0dx\

о

( 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

Соседние файлы в папке книги из ГПНТБ