Файл: Казакевич, В. В. Автоколебания (помпаж) в компрессорах.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 19.10.2024

Просмотров: 150

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

но по потоку и против потока. Видно, что трубопроводы при на­ личии распределенного сопротивления и потока газов являются несимметричными четырехполюсниками. Если распределенное сопротивление отсутствует (щ = р,2 = 0), то волновые сопротив­ ления принимают значения Z\ и Z22, при этом они одинаковы в обоих направлениях.

4.3.ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПЕРИОДИЧЕСКИХ ДВИЖЕНИЙ

ВСИСТЕМЕ

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

Из выражений (4.36) — (4.39) следует, что характер зависи­ мости функций q и р от времени определяется знаком декремен­ та затухания б.

Если б > 0, то с течением времени функции q и р убывают, что соответствует устойчивому движению системы. При этом ус­ ловии в системе будет установившееся течение воздуха, а вся­ кие начальные возмущения будут убывать.

Если б = 0, то в системе могут установиться незатухающие периодические колебания. Наконец, если б < 0, то приращения р и q с течением времени будут неограниченно нарастать, т. е. движение в системе неустойчиво.

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

Для определения комплексной частоты е воспользуемся, как указано выше, уравнениями (4.40), переписав их в виде

ЕпАх+ ЕХ2ВХ+ 0 + 0 = 0;

0 + 0 + Е2зА2+ E2iB2= 0;

:Еп Ах+ Е^2ВХ+ Е22А2+ E2iB2= 0; ■£41^1 + Е±2В\ + Е^А2+ Е^В2= 0.

Для того чтобы эта система имела решения, отличные от ну­ левых, необходимо и достаточно приравнять нулю ее определи­ тель, т. е.

я »

е х2

0

0

 

Д (е ) = 0

0

е £24

(4.43)

£31

£32

£33

£34

 

£41

£42

£43

£44

 

133


Поскольку входящие в уравнение (4.43) коэффициенты Ец{ являются комплексными величинами, то это комплексное урав­ нение распадается на два вещественных уравнения, из которых можно определить частоты со и декремент б:

 

Д(е) = Ai (со, б) + /А2(со, б) = 0

(4.44)

 

Ai(©, б) = 0;

(4.45)

 

A2(w, б) = 0.

(4.46)

В ряде случаев, как будет видно из дальнейшего, уравнения

(4.45) и

(4.46) просто разрешаются относительно декремента и

частоты.

Если в этих уравнениях положить 6 = 0,

то получим

уравнения для определения частоты и амплитуды автоколеба­ ний. При этом предварительно необходимо усреднить в смысле гармонического баланса входящую в эти уравнения величину F', которую следует считать функцией расхода (скорости).

 

Раскроем определитель (4.43):

 

 

 

где

Ае = — En (E23M— E2iN) + Ei2(E23Q— E2iP) = 0;

(4.47)

 

 

 

 

 

M = (0 '_ Z n )e ^ 2+X2)i/i

 

 

 

 

TV= (©I—

 

 

(4.48)

 

Q1= (0 ;_ Z n )e (" ‘ +^ !'

 

 

 

 

 

 

 

Q = (©; _ Z ,' !> e(x,+x»>

 

 

 

 

01 = Ли + Z22hi2h2\~,

 

 

(4.49)

 

02 = Лц + ^22^12^21-

 

 

 

 

 

 

ai

Если предположить, что £ц = —£ 12, что

справедливо при

1 (это соответствует малости скорости течения t/oi по срав­

нению со скоростью звука Ci), то уравнению

(4.47)

можно при­

дать вид

 

 

 

E23(M + . Q ) - E 2i(N+P) = -----Ц- [E2< ( Q - N ) - E 23( Q - M ) l

(4.50)

 

z i1

 

 

 

 

Заметим, что если £ц = —£ 12, то Z'u = —ZJ'( .

 

 

 

Подставляя в уравнение (4.50) значения величин М, N, 0 и Q

из выражений (4.48) и величин Ё23 и Еи из выражений

(4.41),

получаем

 

 

 

 

(Z22—Z2) [(02— Z'Oе^2+х2) + (02— Z,',) е^ 1+^

1,1] h

- ( Z 22— Z2)[(0 ; - Z n ) e ^ 2+^ ) i'1 + ( 02- Z n ) e ^ 1+^ ) J'‘ ] Л* =

134


= ^ L-{(Z22- Z 2)[(0 2 - Z ;i)e(Xl+:4) y\_ (e' _ z ; ,) c(*2■+:ч) уi] eh _ z i i

- ( Z m — Za) [(0 l-Z ii)e (x«+^> y (0i— Zn)e(X2+x‘) 1'«]ex*}. (4.51)

Если воспользоваться выражениями (4.34) и (4.35), то мож­ но преобразовать степени различных членов выражения (4.51):

2 + А.2) У\ + ^ 1 = — Х\ +X 2 + jb[(ai + a 2) y l + a2]\

+ Яг) ух+ Я! = х х— х2+ /в [(а! + а2)ух+ а,];

^ g

(Я2 + Я|) ух + Яг = — Х\х2+ /е [(<Z| + а2)ух+ а2];

(Я[ -)- Я^ ух -+- Я2 = Х\ + х2+ /е [(а] + а2)ух+ а2],

где

' - Л -

Из этих выражений следует, что все члены уравнения (4.51) имеют множитель е,е[(а'+а»>у>+а»1, на который уравнение мож­ но сократить. Тогда получаем

 

(Z22— Z2)[(02— Zu)e дс‘^-Л!:а-1- (02— Zii)er‘ **]—

 

— (Z22— z 2) [ 0 ;- Z [ ,) e - j:‘- j;4 - ( 0 ; - Z ; i)ej;‘+i

=

 

4

- {(Z22 Z2) [02- z ; ,)ex'~ x' + (02- Z n ) e - ^ * ] -

 

z n

 

 

— (Z2 2 — Z2) [0i-^-Z|i)eJ:,+**-l- (0i— ZnJe- *1-"'4] ).

(4.53)

Если

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

ZJj = —Zj',, то окончательно

 

 

(Z22— Z2) (02 + Z|)ch X\ + ( Zi i + 02 —j—^ sh Xi

 

 

= (Z22— Z2) (0[ + Z ,)ch x ( + { Z"i 1 + 0 [ -^ -^ s h x .

e~Xl.

(4.54)

Это уравнение является достаточно сложным. Поэтому изу­ чим влияние параметров системы на величину частоты ш и дек­ ремента б для нескольких частных случаев. Рассмотрим, в част­ ности, случаи, кргда имеется один выходной трубопровод без распределенного затухания (/2 = 0, 1\= 0, р2 = 0), один входной трубопровод (А =т%0, 12 = 0, pi = 0), входной и выходной трубо­ проводы (Л Ф 0, 12ф 0, pi = р2 = 0) и распределенное затуха­ ние (pi = 0, р2 = 0).

135


В дальнейшем воспользуемся следующими

*i

= /* i —

J

x2

= jx2— b2, j

*i =®Pi«/u

j

x2= ©P2( 1

yx)\ J

6, - ( 6P,

J J J w

.о .

обозначениями:

(4.55)

(4.56)

(4.57)

При получении этих выражений мы воспользовались уравнением

(4.30).

4.4. УСТОЙЧИВОСТЬ СИСТЕМЫ С ВЫХОДНЫМ ТРУБОПРОВОДОМ

Рассмотрим систему, состоящую из компрессора и выходно­ го (нагнетающего) трубопровода, на свободном конце которого находится акустическое сопротивление Z2. На входе в компрес­ сор сопротивление Z\\ входной трубопровод отсутствует (Л = 0), распределенное сопротивление также отсутствует (p.2 = 0). При этих предположениях получаем

ха= ©р2;

б2 — бр2;

(4.58)

х х= 6, =

0.

Полагая пока Zx — 0 и принимая, что Z2 = /?2 + jX2, где R2 и Х2— активная и реактивная составляющие акустического сопро­ тивления, получаем из уравнения (4.54)

{— Z22— R2— jX2)(cos х2+ / sin х2) 62e~e* =

=(Z22— /?2— jX2)(cosx2— j sin дс2)0!ев*.

Произведем вычисления и сгруппируем члены:

'* +

r e,=

>e’ +

Iea‘,

где

01 = Ап— Z22hx2h2x\

02 = Ац + Z22A12A2I: !

(4.59)

(4.60)

(4.61)

136

Комплексное уравнение (4.60) распадается на два веществен­ ных:

(— A cosх2 + В sinx2)e~ea = (cos*2C — D sin j^ e*1;

(4.62)

(A sinх2 + В cosх2) e-e> = (С sinjc2 + D cosх2) е®1,

(4.63)

где

D = X201-

А = (Z22-\-R2)Q2, в = *202; C = (Z22

Разрешим уравнения (4.62) и (4.63) относительно декремен­ та 6 и частоты со. Для этого разделим одно уравнение на другое. Получаем

A cos хг + В sin дг2

__ С cos — O s iris

A sin х2+ В cos

C sin * 2 + D cosx 2

Отсюда

АС sinх2cosх2 + ВС sin2 х2ADcos2x2+ BD sinx2cos x2 =

=AC sin jc2cos x2— AD sin2 *2 + BC cos2 x2— BD sin x2cos x2

или

(BC + AD)(sin2 x2—cos2 x2) — 2(AC — BD)cos x2sin x2 = 0,

либо, после деления на cos2jc2,

(ВС + /4D)tg2 х2— 2(АС— BD)tg х2— (ВС + AD) = 0.

Вычисления дают

 

 

ВС + AD == 2QiQ2X2Z22\ AC— BD = Qfi2{Zl2— R\— X%)\

(4.64)

поэтому уравнение частот принимает окончательный вид

 

tg2x2— 2atgx2— l = Q,

(4.65)

где

_р2 у2

 

,

 

1 ^22—Л2 —Л2

 

2

Zt!X2

 

Для получения уравнения декремента_разделим сначала каж­ дое из уравнений (4.62) и (4J53) на cosx2, а затем решим урав­ нение (4.63) относительно tgx2. Тогда

tgx2 = Per*— В е ~ ~ °3

Л е - 6 »— С е в *

Аналогичная операция с уравнением (4.62) дает

. - Сев» + Ле-в* lg ^2 = — ---------т-

Dee* + Ве~°а

(4.66)

(4.67)

137


Приравнивая правые части выражений (4.66) и (4.67), полу­ чаем

(А2 + В2)е-2в’ = (С2 + D2)e26'

или

с- 4б, _ с 2 + D2

(4.68)

А 2+ В 2

 

Подставляя сюда значения коэффициентов А, В, С и D и про­ изводя преобразования, найдем выражение для определения де­ кремента

 

е -46’ (

01 V

(Z22

R 2)

+ Z 2 .

(4.69)

 

 

1

02 )

(Z22 + R2)2 + x2

 

с~ 4б2

/

Ли 2 22Л|2Л2|

\ 2

(Z22

^2) + Z 2

(4.70)

 

V

Ли + Z22/i|2/i21

'

(Z22 + R 2)2+ xl '

 

 

Существенно

заметить,

что частоты

определяются

активной

и реактивной составляющими акустического сопротивления на конце трубопровода.

Уравнение частот (4.65) является квадратным относительно

tg*2:

 

tg х2= а ± У Т + а 2-

(4.71)

Получаем два ряда собственных частот: один — для случая, когда перед корнем взят знак плюс, второй — когда знак минус. Для того чтобы установить, какой из этих рядов соответствует реальным частотам системы, рассмотрим частные случаи. Под­ ставляя значения величины а в уравнение (4.71), найдем

Z22- R 2- X 2 ±

У (Z222- R 2- X l ) 2 + 4X2Z;2

(4.72)

tg *2 = — •

22

Z99X0

 

 

Рассмотрим сначала случай, когда перед корнем взят знак плюс. Если |Z2|<С Z22, то выражение (4.72) принимает вид

(4.73)

Z90X 2 zg2x2

Отсюда следует, что при |Z2|->-0 tgX2 -»-°°. Известно, одна­

ко, что при

|Z2 |—>-0, т. е.

в случае открытого конца,

должно

быть tg х2— О,

22, то вместо выражения (4.72) находим

 

Если |Z21 >

Z

 

 

 

tg * 2

Z22Z2

_|__ }_

Z22

(4.74)

 

 

 

R\ + x\

4

X2Z2

 

138