книги из ГПНТБ / Казакевич, В. В. Автоколебания (помпаж) в компрессорах
.pdfПреобразуем это выражение, подставив в него значения ве личин Ац, 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
о
2л
а,= ——Г F^acos u)cos udu;
яJ
|
2л |
|
|
Ok |
_l_ |
a cos u)coskudu\ |
|
я |
|||
(4.158) |
|||
|
|
||
|
2л |
cos и)sin udu\ |
|
a\ |
я |
||
|
|
||
' |
2л |
|
|
l |
cos u)sin ku du, |
||
ak= ----- |
|||
|
я |
|
|
где и = cat.' |
|
симметрична относительно нача |
|
Если характеристика Fi(^) |
|||
ла, то ао = 0; если симметрия не соблюдается, то а0 Ф 0, что как |
бы эквивалентно смещению центра колебаний. Будем полагать,
для упрощения расчетов, что ао = 0. |
Нетрудно видеть, |
что а \ = |
|||||
= |
0. Действительно |
|
|
|
|
|
|
|
2л |
|
|
|
|
2л |
|
|
а.\ = -----( 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 |
|
|
|
|
о |
|
2Я |
|
|
|
|
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С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