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

книги из ГПНТБ / Казакевич, В. В. Автоколебания (помпаж) в компрессорах

.pdf
Скачиваний:
23
Добавлен:
21.10.2023
Размер:
9 Mб
Скачать

Преобразуем это выражение, подставив в него значения ве­ личин Ац, h\2i h2l:

е - 4 в р , = е “ 2 - f f ( Z 22 + P q \ F + V - i Z y P x f + Iх 2z 2 2 n l x

 

(Z22 p0lF

^2^22m2) + (*2^22п2

 

x [Z22(l —Ц2т 2) — RiY + {Xj + Ц2£22я 2)2

(4.153)

[Z22( I + Щ/n,) + Л2]2 + (Хг + H2Z2Sn,)2

 

В частных случаях открытой

(|Z2 |->-0) и закрытой

(|Z2|

-*■ оо) труб последний множитель правой части уравнения (4.153) становится равным единице, поэтому получаем

е_4вра_ е—2 7Г

(Z22 + P0\F' + I V ^ l )2 + ^ г 22»1

(4.154)

 

(Z22—Р0\^ + ^2^22т г) + 1А2^22п2

 

Анализ выражений

(4.153) и (4.154) позволяет сделать ряд

выводов о влиянии распределенного затухания на величину де­ кремента.

При отличном от нуля коэффициенте распределенного сопро­

тивления (р2 Ф 0) декремент затухания системы б

возрастает

по сравнению со случаем ц2 =

 

0.

Это следует, в частности, из

того, что множитель

-г-й*.

<

1 ;

.так как тi « —т2, то первый

е

множитель выражения

(4.154)

при прочих

равных

условиях

уменьшается при возрастании ц2.

 

(4.154)

ни при ка­

Если р.2 Ф 0, то знаменатель выражения

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

порциональны множителю —2^ 2 , который при и -> -0 и со-»-оо

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

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

4.8. ОПРЕДЕЛЕНИЕ АВТОКОЛЕБАТЕЛЬНЫХ РЕЖИМОВ

При исследовании автоколебаний в распределенных системах используется метод периодических решений Пуанкаре в форме, Витта [11] или метод Линдштедта — Ляпунова в форме Гвоздовера [12]. Однако эти методы сложны.

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

160

Ранее мы показали, что колебательная граница устойчивости определяется условием 6 = 0, где 6 — декремент затухания ко­ лебаний. В частности, для случая 1\ = 0 получены условия дина­

мической и статической устойчивости в виде F' <

D{ (динамиче­

ская устойчивость) и F' < D2 (статическая

устойчивость), где

D\ и D2 определяются выражениями (4.96) и (4.97); на границе

устойчивости F' = D\.

 

в окрестнос­

Рассмотрим характеристику компрессора F(Q)

тях рабочей точки Q * PoF(Q*).

 

 

Пусть характеристика аппроксимируется полиномом

F (Q )= b0 + blQ + b2Q2 + b3Q3 + ... .

(4.155)

Рассматривая характеристику в окрестностях рабочей точки

Q*, получаем

 

 

F(Q* + q) = F(Q*) + F'(Q*)q + Г ^ 4*-

F"'(Q*)q*

 

 

Здесь

 

 

F(Q*) = bo+ 6,Q* + b2Q*2+ b3Q*3+

... =

C0;

F'(Q*) = 6, + 2 b2Q* + 3b3Q*2+ ... = Cl;

F"(Q*) = 2b2+ 6b3Q* + . . .

= C 2;

(4.156)

F"'(Q*) = 663 + ... = C';

Следовательно,

F(Q* + q) = C0+ C[q + -£|-q2+ ~ q3+ ...

Рассматривается приращение давления

F{Q* + q) — F{Q*) = C[q +

q2 + - y - q3 + ...

 

Обозначая F(Q* + q) F(Q*)

= Fi(q), имеем

 

Fl(q) = Clq + C2q2+ C 3q3+ ....

(4.157)

где

Положим q = a cos at и разложим Fi (<7) в ряд Фурье:

Fi(acos(ol) = ao + aIcoso)^ + a2cos2(Di + ... + a Acos^(o^ +- . . . +

+ a\ sin at + a2sin 2 a t + . . . + a* sin k a t + . . . ;

|1 Заказ 1516

161

здесь

а0 = —— Г Fi(acosu)du;

2л J

о

а,= ——Г F^acos u)cos udu;

яJ

 

 

Ok

_l_

a cos u)coskudu\

я

(4.158)

 

 

 

cos и)sin udu\

a\

я

 

 

'

 

l

cos u)sin ku du,

ak= -----

 

я

 

где и = cat.'

 

симметрична относительно нача­

Если характеристика Fi(^)

ла, то ао = 0; если симметрия не соблюдается, то а0 Ф 0, что как

бы эквивалентно смещению центра колебаний. Будем полагать,

для упрощения расчетов, что ао = 0.

Нетрудно видеть,

что а \ =

=

0. Действительно

 

 

 

 

 

 

 

 

 

 

 

 

а.\ = -----( F1(acosu)sinudu = ------- Г F, (a cos u)dd cos и =

 

я J

 

 

 

na J

 

 

 

 

 

2it

 

 

 

 

 

 

 

—1

с?ф(а cos u),

 

 

 

 

 

яа I

 

 

 

 

 

 

 

 

где

dtb(acosu)

,

,

 

 

с перио­

----- J------- - = Fi(acos и) — периодическая функция

 

даcos а----------------------------------------------------------------------- г

дом 2я.

 

и

переменная

составляющая

функции

 

Следовательно

 

•ф(асовц) — периодическая по и с тем же периодом.

162

Поэтому

 

 

 

а[ = —— N)(acos 2n)— Tjj(acos 0)] =

0 .

 

я a

 

 

 

Пренебрегая членами высшей кратности, получаем

 

Fi (a cos и )= F1(q) = aicos и =

 

 

2n

-I

 

 

— J Fi{acos u)cos иdu\a cos и = hq,

(4.159)

о

J

 

 

где

 

 

 

2 л

 

 

 

h = -— Г Fl(acosu)cosudu = Flcpia).

 

na J

 

 

 

Иными словами, h — усредненное

значение

dFt

в смысле

 

 

dq

 

гармонического баланса. Мы предполагаем, что амплитуда ко­ лебаний мала и значение D\ при колебаниях равно значению Dx на равновесном режиме. Иными словами, считаем, что при коле­ баниях величины активного и реактивного сопротивлений не меняются. Тогда амплитуда автоколебаний найдется из урав­ нения

Ficp(a)==T>i. (4.160)

Найдем теперь выражение для коэффициента h в явном ви­ де. Для случая полинома 5-й степени аналитическое выражение для характеристики имеет вид

F\{q) — С\Ч+ ^2<72 + Сзq3+ C t f + C5q5.

(4.161)

Тогда

h — —— Г [Cya cos и + C2a2 cos2 и + C3a3 cos3и + C4a4cos4и +

na J

 

 

 

 

о

 

 

 

 

 

j^ a - rtcos 2и

+ C6cos5 u + G5a5cos5 u]cosucfu = —!—J

 

 

0

 

 

+ C2a2 cos и + -j- cos 3uj + C3a3

+ -j- cos 2u +

cos 4uj +

+ Cta4 ( — cos и + —— cos 3u + —5— cos 5*Л +

 

\ 8

16

16

/

 

+ C5as ( - ^ - + -^ -c o s 2u + —^-cos4u + —— cos5iAl du.

\ 16

32

16

32

J j

l l

163

Все члены, пропорциональные косинусам кратных дуг после интегрирования в пределах от 0 до 2 я, в силу периодичности да­ дут нуль. В результате имеем

Л = С[ + -^ -С 3а2 + -^- С5а4.

(4.162)

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

4.9.МЯГКИЙ И ЖЕСТКИЙ РЕЖИМЫ ВОЗБУЖДЕНИЯ КОЛЕБАНИЙ,

ИХ УСТОЙЧИВОСТЬ

Если взять с1 > Du то это соответствует случаю мягкого воз­ буждения колебаний, которое характеризуется самовозбуждени­ ем колебательного процесса, то есть неустойчивостью положения равновесия; если С\ < D u то возможно жесткое возбуждение ко­ лебаний.

Определим величину амплитуды автоколебаний

для случая

Z, = 0.

 

 

 

 

 

 

 

Рассмотрим вначале случай мягкого режима. При этом

 

 

^ i ( < 7 ) = C i <7— С з <73;

 

(4.163)

 

 

 

3

 

 

 

h(a) = С ,---- - С 3а2.

 

 

 

 

 

4

 

 

 

 

Подставляя выражение для h = F c'p

в уравнение (4.160),

по­

лучаем

 

 

 

 

 

 

 

 

Ficp(a) = С [-----—C$a

= D l7

 

 

откуда

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

а " - 2 /

^ д Г -

(4 Л 6 4 )

На рис. 4.8,

а графически показано условие существования

устойчивых колебаний. При Ci > jD,

квадратная

парабола /,

отображающая

кривую

F [ср

(а),

 

пересекает

прямую

D\

в точке с абсциссой аст.

Если

Ci <

D\, то пересечения

нет

и периодическое движение невозможно (случай, показанный кривой 2).

Для случая 1\ = 0

 

/?]

Za2

V (Z22+ /?2)2 + Xj + V (Z22— R2)2+

Xj

 

/7°‘

V ( Z 22 + R2)2 + x l - l / ( Z 22- / ? 2)2 +

*2

164

Поэтому

Oqt

2

 

К

(Z22 + R2f

+ х\ + /

(Z22^ 2f

+ Х\

У1 с 8

 

V

(Z22+ R2)2+ x \- V

(z22- R 2f

+ x l

 

 

 

 

 

 

 

 

 

 

22

 

 

(4.165)

В частности, при X2 = 0 D

 

 

 

 

 

------и стационарная амплитуда

определяется выражением

 

 

PoiRt

 

 

 

 

 

 

 

 

 

 

Яст

 

 

1 /

”С

Z 22

 

 

 

К*

зс3

У

1

РоЛ ■

 

 

 

 

 

 

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

F[(q)=C i— 3Czq2.

На равновесном режиме q — 0; тогда

'F[(Q) = Ci.

Вто же время амплитуда автоколебаний разыскивается из условия

.Fcp=Ci — —+■С3аст= Di.

4

165

Следовательно,

С, — D\ н—— Сга2 > D],

4

т. е.

F [(0 )> D .

Итак, условия устойчивости равновесия не выполняются, рав­ новесие неустойчиво и режим возбуждения — жесткий.

Далее, из выражения

следует, что при Cj > D 1 существует только одна стационарная амплитуда (при С\ < Dx— корень мнимый, что подтверждает невозможность периодических движений в этом случае).

Выскажем теперь качественные соображения об устойчивости периодических движений. Если а < аст, то h(a) < Du и ампли­ туда колебаний растет; если же а > аст, то А (а) > Dь и ампли­ туда убывает. Следовательно, периодическое движение устойчи­ во по амплитуде (орбитально-устойчиво). Поэтому можно считать доказанным, что рассматриваемый режим возбуждения колебаний является мягким, а в системе имеется одно устойчи­ вое периодическое движение. В системе с компрессором будут происходить помпажные колебания с амплитудой, возрастающей от нуля при прикрытии дросселя.

Сравним теперь выражение для стационарной амплитуды, полученное в частном случае Х2 = 0, с выражением, найденным в предположении, что система описывается уравнением с сосре­ доточенными параметрами и 1\ = 0. Мы получили выше

Учитывая принятые обозначения, имеем b = Ci, с == Сз,

Следовательно,

т. е. оба выражения совпадают.

Рассмотрим теперь случай жесткого режима возбуждения при наличии устойчивых колебаний. Он имеет место, как будет показано ниже, если Ci < D\ и С5 < 0.

166

Амплитуда автоколебаний в силу уравнения

(4.160)

будет

определяться уравнением

 

 

 

 

h(аст) = С, + -j- Сйа%

^ Сба*т = DX.

(4.167)

Перепишем его в виде

 

 

 

 

flcx-

4

-тг- «ст + ( D . - C ,) - ^ - = 0,

 

 

 

5

Сб

 

5Сб

 

 

откуда

 

 

 

 

 

 

 

а ст(1,2)

5

Сб

Э

С ^ - б - в ф

. - а д .

 

 

5Сб

 

 

 

Отсюда следует, что в системе возможны два периодических

движения с амплитудами aCTi и аСТ2, если

 

 

 

 

9Сз— 5-8(D, — С,)Сб > 0.

 

(4.168)

При этом

 

 

 

 

 

 

 

асri = j / 1 5 7 [ Сз + Т ^

 

 

 

°ст2= |

/ - ^ r f Cs— \-V sC \--4 0(D ,-C 1)C5.

 

Если Ci < Du то самовозбуждения колебаний не будет, так

как положение равновесия q = 0 устойчиво.

что при Д > 0

Возьмем а2 =

а 2т1

+ Д,

|Д| «С 1

и покажем,

происходит убывание амплитуды колебаний, а при А <

0 .— воз­

растание амплитуды.

 

 

 

 

 

Имеем

 

 

 

 

 

 

 

а< = Г-L . ^

+ - 1 - V 9Сз— 40(£>!— Q C j + д]* =

L 5 С5

 

оС5

 

 

J

 

= 4 - ■ 4 - + - V

Е9 С '

- C i ) C s l + Л 2 +

 

25

С\

25С\

 

 

 

 

+ 2 4

•4 - •

Vr9C23-40(Z)l-^€,)C5 +

 

5

 

Сб

5Сб

 

 

 

+ 2 4 -

Д+

К 9C|=-40(D, - С ,) С 5 Д;

 

5

 

Сб

 

оС5

 

 

 

167

Подставляем значения а4 и а2 в выражение (4.167):

 

 

h(a) = h (V а\т\+ Л) = С, + ^ -С 3

+

 

 

+ - ^ - V 9 C !- 4 0 ( D ,- С,)С51 -

 

 

 

 

5 С 5

 

J

 

5гС5 1 ^ -. 4 - + - LT

[9C l-5.8(D 1- C 1)C6] +А 2 +

 

8

\ 25

С\

2ЬС\

 

 

 

 

+ 2 ± .-^ ~

. - L - V 9C23-40(D ,-C1)C5 +

 

 

5

С5

5

 

 

 

 

+ 2 | 4 А + - ^ K 9 C | _ 4 0 (D 1- C ,) C 5a 1.

 

 

5 С5

 

0С5

 

)

После упрощений получаем

 

 

Л(1Л&,+Д) = D1

С5Д2

1- Д / 9 С 23-4 0 (Л 1- С 1)С5.

Из этого выражения видно, что характер изменения h(а) в

окрестностях астi определяется знаком Д: при Д >0 h (V а2т »+ д )<

< Ь ,.

Если же Д <

О, то при Д достаточно малых

по модулю

л ( К < £ ,— a )< D ,.

Следовательно, амплитуды, большие aCTi. убывают, а мень­

шие a0T1, но близкие к ним возрастают, откуда

вытекает, что

большая стационарная амплитуда aCTi устойчива.

амплитуду аСТ2-

Рассмотрим теперь меньшую стационарную

Возьмем а2 = а ^.2

 

+ Д, |Д| С 1.

 

Подставляя а2 в выражение для h(a), получаем

й (^ а ?т2 + д) = С, + 4 -Сз Г А ._ £ з _ _

 

 

V 9 C S -4 0(Dl _ С ,)С 6+ д]-

- Т С‘ { ^ -

|

- + -5 j| - [ 9 C ! - 4 0 (O1- C |)Csl + д *+

+

^ 9 С | — 40(D, - С,)С6+

 

25eg

+" Г - “| _ А _ Л - ^ Г >^ 9С|— 40( D, — С,) СБ} .

168

Отсюда

h { V «Its + А) = D, —4 - C5A2 + — A V 9Сз— 40(D,— C]) C5.

 

8

4

 

Рассмотрим это выражение. Если Д > 0, то при достаточно

малых А

 

 

 

h(Va!r2 + а ) = D,

g- С5Д2 + -L А V

9С^— 40(D, — С,) С5 > £>,.

Следовательно,

амплитуды,

чуть

большие аСТ2, возрастают,

стремясь к значению а = aCTi.

 

 

Если А < 0, то

 

 

 

h ( V ас2т2 + А) = D, —4 С5Д2— 4

А V

9Сз— 40(Л>! — Сх) С5< £>,.

Поэтому амплитуда будет убывать, стремясь к нулю.

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

На рис. 4.8, б показано графически условие

существования

колебаний в случае жесткого режима.

При

Сх < Dx

кривая

F'lcp (а) пересекает прямую D\ = const

в двух

точках,

опреде­

ляющих стационарные амплитуды астi и аст2.

(а)

изображается

Если Ci = C j> D | , то зависимость

FJ

кривой 2. При этом существует только одно устойчивое периоди­ ческое движение со стационарной амплитудой а ^ при мягком

режиме возбуждения.

Нетрудно установить характер устойчивости стационарных амплитуд. Если кривая F'lCp (а) пересекает при возрастании а

прямую D! снизу вверх, то стационарная амплитуда, соответст­

вующая точке пересечения,— неустойчива

(на рис. 4.8, б точ­

ка /4). Если же кривая F [ (а) пересекает

при возрастании а

прямую D\ сверху вниз, то стационарная амплитуда устойчива (точка пересечения В)‘,

Вслучае точки А соседние значения -амплитуд будут удалять­ ся от аст2, что условно показано стрелками, исходящими из точки А.

Вслучае устойчивой амплитуды aCTi соседние амплитуды бу­ дут приближаться к точке В, что также отмечено направлением стрелок.

Если кривая FJcp (а) касается прямой D\ == const (кривая 3),

то имеется полуустойчивое периодическое движение с амплиту­ дой а "т. Здесь амплитуды, большие а"т, убывают, стремясь к

а"т, а амплитуды, меньшие а "т, убывают до нуля. В этом слу­

169

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