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

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

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

При перемене числа оборотов и неизменном положении дрос­ селя может меняться левая часть этого неравенства, а в правой части будет принимать другие значения только величина к (так как полагаем F <С 1 и рк « ро), Но, как известно, величина на­ пора пропорциональна л2, а расхода — пропорциональна л, т. е.

„2

 

 

 

 

£ ki_ _

_ ^ i _ .

 

 

 

(1.28)

 

 

 

 

Рк2

2

*

 

 

 

 

 

 

 

п2

 

 

 

 

 

 

 

 

 

Qki _

п\

 

 

 

 

(1.29)

 

 

 

 

QК2

Лг

 

 

 

 

 

где рк -

ролк ро = F{QK) — р0 = F(QK).

при изменении

числа

 

Поэтому характеристики вентилятора

оборотов будут меняться так, как показано на рис.

1.1, а.

Будем

приписывать функции F(QK) и ее аргументу индекс, совпадаю­

щий с индексом числа оборотов.

Например,

числу

оборотов Л|

будет соответствовать характеристика /•'i(Qki).

 

 

 

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

 

 

 

 

 

 

 

По

 

- Ц - F , ( Qk2 - Ь - ) =

F2(Qk2).

 

Рк2 Рк1

=

- Ц -F x ( Q k i ) =

( 1 . 3 0 )

 

 

ft?

ft?

ft?

\

*2

>

 

 

 

Далее

F2(Qk2) — <P(Qk2)-

30

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

о

 

- J - / 7. ( Q k2 - ^ - )= < P (Q k2).

 

Пусть это значение будет равно Q^,.

 

Подсчитаем величину

^Ркг

в точке Q*2

 

 

 

dQK2

 

 

 

dp**

4

 

 

^

)

< 2

«I

d ( q*2

) rfQ*.,

 

 

 

 

dF, (Qki) Ч\

(1.31)

или

 

 

 

 

 

 

dF2 (Q k2) _

n2

d F i (Q ki)

 

 

< 2

 

"•

dQ‘Kl

 

Выясним теперь,, как будет изменяться величина k. 1 Пусть ре = aQ l , тогда

k = _d£§_ _ 2aQK. dQK

Если п = п\, то

Qk— Qki» k — — 2aQKi .

При п = П.2

Qk2 = Qki — ;

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

k = k2 = 2aQKl^ - , «I

t . e. с возрастанием n возрастает и k.

Это обстоятельство, как легко видеть, имеет место в общем

случае характеристики дросселя, если только производная dpe dQк

возрастает с увеличением QK (поскольку увеличение числа обо­ ротов влечет за собой возрастание расхода QK).

Итак, если рабочая точка лежит на восходящем участке ха­ рактеристики, то с повышением числа оборотов и при неизмен­ ном положении дросселя левая часть неравенства (1.27) моно­ тонно возрастает, а правая, вследствие увеличения в знаменате-

31

 

ле величины k, убывает, стремясь к нулю, когда п-*-оо. Вслед­

 

ствие этого область устойчивости с увеличением п будет

 

уменьшаться. При этом всегда можно указать такое достаточно

 

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

 

начнется помпаж.

 

 

 

 

область

устойчивости

 

Наоборот, при уменьшении оборотов

 

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

I

указать такое достаточно

малое число оборотов,

при

котором

помпаж невозможен.

 

 

 

Предположим, что рабочая

■"

Влияние изменения температуры.

 

точка характеристики соответствует мягкому режиму возбужде­

 

ния. Как указывалось выше, в этом случае область устойчивости

 

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

 

рдс2/ -

 

 

 

 

 

ар

 

 

 

(1.32)

 

dQ*

 

-— -.

 

 

 

 

ksls2l2

 

 

 

 

Следовательно, граница устойчивости

будет соответствовать

 

рабочей точке, для которой

 

 

 

 

 

 

 

dF

__ dFrр

__ Росо li

 

 

 

dQ*

dQ*

ks\s^

'

 

 

 

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

 

границе устойчивости и пусть температура наружного

воздуха

 

уменьшилась от То до T0i.

 

 

 

 

 

 

 

 

Обозначая параметры р, с и А, соответствующие температуре

 

Той индексом 1, получаем

 

 

 

 

 

 

 

 

Т0

 

— Pi

Т0

 

 

 

р! = ро— ; рк = рк- ^ р к-г~.

 

 

 

101

 

 

Ро

*01

 

 

 

При изменении Т0 и р0 объемный расход на выходе, отнесен­

 

ный к значениям Тк и рк,

меняется,

так как зависит только от

 

числа оборотов вентилятора.

 

 

 

 

 

 

 

Следовательно, характеристика вентилятора примет вид

 

P*i = ^ F ( Q k) = Fh(Qk).

 

 

 

_

'01

 

 

 

 

 

 

 

Тогда

Т0

dF(QK)

_ ,df„

 

 

 

dpKl

 

 

 

dQK

Т0l

dQK

dQK

 

 

 

Уравнение сети (рис. 1.1,6)

примет вид

 

 

 

Pei = Р6

 

= ^г- Ф(<2к)-

 

 

 

 

Рк

Т

о.

 

 

 

Для определения равновесного режима приравниваем значе­ ния ркI и pel и получаем

!^ ( 0 к ) = ^ ф ( < ? к ) 101 *01

32

или

^ ( Q k) = < P ( Q k)-

Отсюда следует, что значение Q* , соответствующее равно­

весному режиму, не меняется при изменении Го. Далее имеем

d-Рб _ T0d<?(QK)

dQK T0idQK

Наконец, так как скорость звука c — YkgRT, то

Подставляя измененные значения в условие, определяющее границу устойчивости, находим

dFrp

 

 

 

 

 

dQK

*0

Tо

SiSil2

k^SiSzli

dQк

 

_

 

 

 

 

1 01

 

 

 

Отсюда следует, что при изменении Г0 граничное значение тангенса угла наклона касательной к характеристике вентиля­ тора не меняется. Однако при понижении температуры крутизна характеристики вентилятора, как показано на рис. 1.1, б, растет, так как

dFH _ dF T0

dQK dQK T01

Учитывая, что величина Q* при этом сохраняется неизмен­

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

Влияние изменения давления наружного воздуха. Пусть дав­ ление наружного воздуха изменилось от р0 до poi. Тогда при не­ изменном положении дросселя будем иметь

Pta = - £^ F ( Q Kl) = Fp(Q«ti

(1.33)

Po

 

 

Р1б = Рб—

= — <P(QKl)-

(1.34)

Po

Po

 

3 З а к а з 1516

33

Так как для равновесного режима рк = рei,

получаем

Ро F ( Q k 1 ) =

Ра

ф ( Q k i ) -

 

Следовательно, равновесное значение Q*,

определяется вы­

ражением

 

 

 

^(Qki) = <p(QkO-

 

В то же время равновесное

значение Q*

при неизменном

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

т) = ф(0*к).

Так как уравнения одинаковы, то и корни их должны быть одинаковы. Поэтому Q^i — Qk

Следовательно, объемный расход при перемене давления на

входе и неизменном

положении

дросселя

не меняется

(см.

рис. 1.1, в).

 

 

 

 

=

с0 и

 

Так как Т0 = const, то скорость звука с,

 

dpK, _ Pot

dF(QK)

dTp(Qk)

 

 

dQx

Po

dQк

 

dQK

 

 

Далее имеем

 

 

 

 

 

 

 

£ _

dpiл _

Pqi d<p(Qk) _

Poi

£

 

 

dQK

Po

dQK

 

po

 

 

Подставляя найденные выражения в формулу (1.32), опреде­

ляющую область устойчивости, получаем

 

 

 

^ Р№к) ^

Ро( Ро ) 0

1

Ррс0 *1

Ро,

 

<

k

 

 

ks's2‘*

ро

 

 

Ро

 

 

 

 

 

 

Отсюда следует,

что при уменьшении

давления граничное

значение тангенса угла наклона касательной

к характеристике

вентилятора уменьшается пропорционально величине

Но

 

 

 

 

 

 

Ро

 

наклон характеристики вентилятора убывает, как показано вы­ ше, по такому же закону:

dFp

Ро,

dF

dQK

Ро

dQK

Следовательно, значение неприведенного объемного расхода, соответствующее границе устойчивости, не меняется, и область устойчивости сохраняется неизменной.

34

Разумеется, если перейти к приведенным расходам, то гра­ ница устойчивости будет сдвигаться влево, т. е. область устой­ чивости будет возрастать. Отсюда следует, что с подъемом на высоту область устойчивости, отнесенная к параметрам на выхо­ де, будет меняться только при изменении температуры.

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

После публикации наших исследований появились статьи, в которых утверждалось, что диапазон устойчивой работы не за­ висит от атмосферных условий. В работе [22], например, говорит­ ся: «В ряде работ, посвященных анализу устойчивости компрес­ сорных систем, материал излагается таким образом, что соз­ дается ложное представление, что диапазон устойчивой работы компрессора должен зависеть от параметров воздуха перед входом в систему, например от температуры Тн.

Вместе с тем, из соображений размерности ясно, что при представлении характеристик компрессора в параметрах подо­ бия n * = n * ( q ) или я* = л * (р ,п р и в ) точка на характеристике, соответствующая режиму самовозбуждения, не должна зависеть от величин давления, температуры, плотности газа в системе».

Рассмотрим идею приводимого доказательства и установим, в чем заключается допущенная там ошибка.

Автор развивает приближенный энергетический метод иссле­ дования, основывающийся на предположении, что среднее за период колебаний изменение потока механической энергии в. компрессоре должно быть равно изменению потока механичес­ кой энергии на дросселе. Он пишет: «Рассмотрим теперь работу системы при различных значениях рн и Тн, однако при физиче­ ских оборотах ротора компрессора, устанавливаемых в каждом

случае таким образом,

чтобы величина п„р =

пф

остава-

лась неизменной».

 

 

 

I/ - 5 -

 

Характеристика компрессора берется в виде

 

 

 

я« — ~р^~ •Р(Мпр),

 

 

где м „Р= м

Рн норм

 

' н

массовый расход,

приве­

д я

у

1 Н норм

 

 

|

/

 

 

 

денный к нормальным атмосферным условиям.

 

Характеристика дросселя представлена

в виде

 

■Мпр.др = ф ( л др)>

где Л4пр.др — величина расхода через дроссель, приведенная по

параметрам окружающего систему воздуха, а ядр — Р д р

Рн

3*

35

В результате исследования автор получает для границы ди­ намической устойчивости выражение

L

Ф'Спр ’

где F' — крутизна характеристики компрессора; L — акустическая масса;

Ф' — крутизна характеристики дросселя; Сп'р — приведенная акустическая гибкость.

Это условие совпадает с ранее найденным нами, но записано в приведенных параметрах подобия.

Далее автор говорит: «Итак, граница устойчивости может

быть определена в параметрах подобия по величине

=

ф Сар

вне зависимости от параметров газа в окружающей среде». Теперь ясно, что ошибка автора заключается в требовании

постоянства приведенных параметров, в том числе параметра

^пр— ^физ УЪ- Но при изменении Тн и постоянном Пф„3 ме­

няется Лдр, а это, как мы видели, приводит к изменению границы области устойчивости.

Пусть температура Тн стала иной, тогда иной будет ппр. Мы видели, что с уменьшением числа оборотов область устойчивости возрастает и наоборот —- с увеличением параметра пф„3 — умень­ шается. Если Тн уменьшается, то для сохранения ппр = const должно уменьшаться Лфиз. Итак, для того чтобы при уменьше­ нии температуры граница устойчивости сохранялась, нужно уменьшать физическое число оборотов. Значит, если не умень­ шать Лфиз, то устойчивость потеряется. Следовательно, уменьше­ ние Та ведет к снижению устойчивости, поскольку для сохране­ ния рабочей точки на границе устойчивости пришлось уменьшать число оборотов, что привело к повышению устойчивости, ском­ пенсировавшему ее снижение, вызванное понижением Тн- Анало­ гичным образом можно показать, что увеличение Тн ведет к увеличению области устойчивости.

Далее в работе [22] рассматривается граница статической устойчивости, определяющаяся соотношением F' = Ф', и делает­ ся вывод о независимости ее от рн, Тн, рн- Этот вывод является справедливым и легко вытекает из проведенного выше анализа.

1.5.ПОВЕДЕНИЕ СИСТЕМЫ В ОБЛАСТИ САМОВОЗБУЖДЕНИЯ

ИСТАТИЧЕСКОЙ УСТОЙЧИВОСТИ («В МАЛОМ»)

Рассмотрим поведение колебательной системы и характер самовозбуждения в ней при изменении ее параметров. Для этого

построим плоскость параметров k и dp линеаризованной сиdQ*K

36

стемы, получаемой из уравнения (1.15), если положить dF dQ

dF

dQ*

. Характеристическое уравнение такой линеаризованной

системы, описывающей поведение колебательной системы «в ма­ лом», будет иметь вид

v2 +Y - kCj ra

Lg^ V + —

) - о .

(1.35)

\

и )

LaCa \

k )

 

Здесь мы обозначили

Г - ( - * * - )

..

\ dQK JQK-QK

Корни характеристического уравнения (1.35) будут

Условие существования комплексных корней имеет вид

(La— kCaF')2— 4CaLa(A2— kF)<. 0.

(1.37)

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

(La + kCaF ' f - { 2 k V c X f <

0

 

или

 

 

(La + kCaF' + 2k V C aLa) (La + kCaF2k V C aLa) < 0 .

(1.38)

Границе области комплексных корней

соответствуют

два

уравнения 2-й степени:

 

 

La + kCaF' + 2kV cJ7a = 0-

(1.39)

La+ kCaF'— 2 k V c aLa = 0.

 

(1.40)

Каждое из этих уравнений определяет гиперболу, отнесенную

к асимптотам: одной из асимптот для обеих

кривых является

ось F'\ другой асимптотой для первой кривой является прямая

а для второй — прямая

 

 

На рис. 1.2 эти ги­

 

перболы

представле­

 

ны

кривыми 1

и

2,

 

асимптоты — прямыми

 

3

и

4.

Пространство

 

между ними

представ­

 

ляет

собой

область

 

распространения

ком­

 

плексных корней.

 

 

 

Область

динамиче­

 

ской

устойчивости

«в

 

малом»

соответствует

 

положительному

коэф­

 

фициенту

при

 

первой

 

производной, т. е. оп­

 

ределяется

условием

 

F'-----^ - < 0 .

 

(1.43)

 

 

Границей

области

 

динамической

 

устой­

 

чивости

является

ги­

 

пербола

 

 

 

 

 

Рис. 1.2

F' =

С*а

•4 “ >

 

(1-44)

 

 

 

Л

 

 

 

имеющая асимптотами оси F' и А. На рис.

1.2 она показана кри­

вой 5.

 

 

 

 

 

следует,

Из условия динамической устойчивости «в малом»

что в областях распространения узлов и фокусов,

расположен­

ных ниже этой гиперболы,

режим будет устойчивым

и

 

в этих

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

Область статической устойчивости соответствует положитель­ ному свободному члену характеристического уравнения, т. е. оп­ ределяется условием (1.17) A > F'. Граница области статичес­ кой устойчивости определяется условием F' — А, представляю­ щим собой уравнение прямой 6. Эта прямая касается гиперболы

(1.39) в точке, где ___

Над прямой 6 расположена область распространения седел. Она соответствует статически неустойчивым точкам пересечения Л и С характеристики вентилятора и сети (см. рис. 0.4).

Таким образом, рис. 1.2 дает полную картину поведения сис­ темы «в малом» в зависимости от изменения параметров F' и А.

38

1.6.АНАЛИТИЧЕСКОЕ ОПРЕДЕЛЕНИЕ

АМПЛИТУДЫ И ПЕРИОДА ПОМПАЖНЫХ АВТОКОЛЕБАНИЙ

Произведем в уравнении (1.16) замену независимой перемен­ ной.

Пусть

I

 

dF_

 

 

 

Vm

 

где kx= k

 

 

 

 

 

 

dQ

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ

_ rw

1

d*Q r,„

1

 

 

dt

4

Vm

dP

 

m

 

где штрихами обозначено дифференцирование по т.

(1.15),

Подставляя выражения для Q' и Q"

в уравнение

имеем

 

1

 

 

 

 

 

 

Q"

 

dF

L

Q '+Q = 0 .

(1.45)

 

 

 

dQ

 

V'kkfizU

 

 

 

 

Мы получили уравнение вида Ван дер Поля.

Для его исследования применим разработанный нами при­ ближенный метод [21], позволяющий найти ряд величин при лю­ бом значении нелинейности (см. приложение III).

Прецизируем вид функции F. Предположим, что она аппрок­ симируется кубической параболой

F(Q* + Q) = F(Q*) + b Q - c Q \

тогда

+= b -3cQ2.

dQ

Следовательно, уравнение примет вид

« " - / ^ [ ( 6- ^ ) - 3cQ2] Q' + Q - 0' <‘ -46)

Из изложенного в приложении III метода следует, что вели­ чина периода колебания уравнения

х " — р(а + у*2 — Ьх4)х' + х — 0

(1-47)

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

TH = 2 ^ f я2 + -^ -^ 2 а + уа2---- ^-6а4^2,

(1.48)

где а — амплитуда автоколебаний расхода Q.

39

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