
книги из ГПНТБ / Казакевич, В. В. Автоколебания (помпаж) в компрессорах
.pdfПри перемене числа оборотов и неизменном положении дрос селя может меняться левая часть этого неравенства, а в правой части будет принимать другие значения только величина к (так как полагаем 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 + kCaF— 2k 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