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

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

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

Подставляя в уравнение (3.61), найдем

d W

dx =

■оил плп

dp' =

(

dW \

y R T

дх

y g R T

 

 

-----= ----I “

'

 

 

 

dt

dt

 

\

дх )

 

 

 

Допускаем,

что

d W

 

W 2 — «7,

 

 

д х

 

l

 

 

 

 

 

 

 

 

 

Это приводит к выражению

 

 

 

 

 

- ^ . =

(«7

Ц72) Ж

 

(3.62)

 

 

dt

 

1

21 у

 

 

где V — объем, выделяемый сечениями 1 п 2. Учитывая, что W = Qp, имеем

рс2

Переходим к полным давлениям и температурам:

 

 

у

 

Pt — Р ^1 +

М2

V - 1 7\ =

]•

 

 

= Г [ 1 +

Будем полагать, что изменение скорости потока, текущего че­

рез ступень, мало. Тогда можно приближенно считать М =

const.

Дифференцируя при этих условиях pt по времени, получаем

dp(

 

1 y -1 M2

v

 

dp

v-1

 

dt

 

dt

 

 

 

 

Подставляя выражение (3.62), имеем

у

 

 

 

 

 

 

— ---- (И?! —

1 + - ^

M 2 у—1

 

dt

 

 

 

 

 

с учетом зависимости T(Tt)

 

 

1

 

 

 

 

Y-1

 

dt

 

1MEL

М2 у—1

 

 

 

 

 

 

Если число Маха не очень высоко, порядка 0,3—0,4, то

 

 

 

1

 

 

 

 

1 +■

М2 у— 1

 

 

 

 

^ - = (U71_U72) -lME_Tt

(3.63)

или

dt

1

V

 

 

 

 

 

 

 

 

Ca-^ j- = Ql- Q

2,

 

 

 

 

d t

 

 

 

120

где

Как следовало ожидать, такая идеализация оказалась не­ удовлетворительной. Тогда наряду с уравнением (3.63) было рассмотрено уравнение вида (1.7) в форме

= (P/i — P t i ) + Ар/2Ь

где Др/21 — приращение полного давления, создаваемое лопастью ступени в соответствии с ее характеристикой.

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

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

Аналогичным путем была исследована устойчивость турбо­ вентиляторного двигателя; результаты расчета оказались доста­ точно близки к экспериментальным.

Было изучено влияние формы характеристики компрессора и геометрии ступеней на область устойчивости.

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

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

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

121

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

лопастям может возникать нелинейный резонанс.

помех — в

При тех же условиях исследовался

другой вид

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

Эти помехи

фиксиро­

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

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

Некоторые вопросы, связанные с мягким или жестким помпажом, рассматриваются в литературе [15]. Здесь в основу иссле­ дования неустойчивости берется «вариационный принцип макси­ мума потока механической энергии» [14], согласно которому на­ блюдаемое (устойчивое) движение вязкой среды отвечает уело-. вию максимума потока механической энергии через определяю­ щее сечение. Под определяющим сечением в практических рас­ четах принимается поверхность, на которой концентрируются процессы энергообмена и диссипации.

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

В предположении р = const рассматривается характеристика

Н = f(Q ), где Н — полный напор.

В качестве

определяющего

сечения берется плоскость рабочего

колеса.

Отмечается, что

жесткий режим происходит в области, близкой к Нтах. Мягкий и жесткий помпаж возникает на восходящей ветви характерис­ тики самопроизвольно.

Условия неустойчивости жесткого режима, согласно «вариа­ ционному принципу максимума», принимают вид

дЦНО) < 0 ;

дН > 0, где

Q — Р^я*

dQ2

dQ

 

Тогда на характеристике компрессора определяются две точки:

одна — соответствующая

изменению

знака

дН

,

,

------

(точка 1 на

 

 

 

 

dQ

 

 

рис. 3.11, где Н = # Юах)

и точка 2',

где

д

 

изменяет знак

(точка 2 — такая, в которой dH > 0

и

достигается

максимум

д (HQh *

dQ

что правая ветвь характе-

 

величины —5— 0 . Автор утверждает,

dQ

122

ристики,

где

дН

Ои

дЦНО)

 

< 0 , И,

дНСа

dQ

dQ2

 

 

 

 

 

’ дСа

обеспечивает абсолютную устойчивость

 

потока.

 

 

 

 

 

 

 

В точке 1 возникают условия для

 

возникновения жесткого помпажа, а даль­

 

нейшее дросселирование вызывает пере­

 

ход к режиму, начинающемуся

с

точ­

 

ки 2, где наступает абсолютная неустой­

 

чивость движения, т. е. возникает

мяг­

 

кий помпаж.

 

рассматривает

вопрос

 

Далее

автор

 

об образовании зон срыва.

 

преды­

 

Из сравнения с материалами

 

дущих глав очевидна ошибочность полу­

 

ченных в работе [15] результатов.

 

 

 

В работе [35] несколько страниц по­

 

священо вопросам помпажа, которые

 

рассматриваются на неправильных

предпосылках, поэтому да­

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

РАЗДЕЛ II

ПОМПАЖ В СИСТЕМАХ С РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ.

ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ; ВОЗДЕЙСТВИЕ НА ПОМПАЖ

ГЛАВА 4

ИССЛЕДОВАНИЕ ПОМПАЖА В РАСПРЕДЕЛЕННЫХ СИСТЕМАХ

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

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

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

ми результатами показывает, что пока основная

предпосылка

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

системы систе­

мой с одной степенью свободы — оправдывается,

результаты

упрощенной теории

подтверждаются экспериментом

с хорошей

точностью. В то же

время, как видно из работ [39]

и, частично,

[1], в ряде случаев наблюдаются скачкообразные изменения час­ тоты и амплитуды колебаний, форма колебаний часто далека от гармонической, встречаются постепенные переходы от одной частоты к другой.

Все эти существенные обстоятельства не получают должного объяснения в рамках теории систем с сосредоточенными посто­

янными, между тем как они имеют важное значение

с научной

и прикладной точек зрения.

можно, лишь

Выяснить причины указанных особенностей

рассматривая систему, включающую компрессор

с

сетью, как

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

Исследование распределенной системы, включающей ком­

1 2 4

прессор с коммуникациями в рамках линеаризованного прибли­ жения, выполнено в работах (5 и 6], по которым и ведется даль­ нейшее изложение с необходимыми дополнениями и изменениями. Кроме того, развивается приближенный метод исследо­ вания нелинейной распределенной системы, учитывающий нелинейные свойства характеристики компрессора. В результате выводятся соотношения, определяющие амплитуды и частоты периодических колебаний, находятся условия существования мягкого или жесткого режима возбуждения автоколебаний, оп­ ределяется характер устойчивости периодических движений и даются аналитические условия, накладываемые на характеристи­ ку компрессора, при выполнении которых возможно возбужде­ ние колебаний на нисходящих ветвях характеристики компрес­ сора.

4.1. УРАВНЕНИЯ ДВИЖЕНИЯ СИСТЕМЫ

При выводе уравнений сделаем некоторые упрощающие пред­

положения.

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

На рис. 4.1 представлена схема рассматриваемой системы. Ко входу и выходу компрессора присоединены трубопроводы со­ ответственно длин /i и /2 и поперечных сечений Si и S2-

Поведение системы может быть описано посредством уравне­ ний гидродинамики в форме Эйлера и уравнений состояния

(4.3)

где U, р и р — скорость, плотность и давление, являющиеся функциями времени t и координаты х;

у— показатель диабаты;

аи т — коэффициенты внутреннего и турбулентного

трения;

d — диаметр трубопровода.

125

 

 

St

St

Рис. 4.1

 

 

 

 

 

V

 

 

 

 

A

 

X

 

 

 

 

T

( t

LA

 

*

 

 

x

 

 

Уравнения (4.1) — (4.3)

должны быть взяты для каждого из

трубопроводов системы.

Начало координат поместим в начале

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

другого компрессором и на их концах помещены акустические сопротивления Z\ и Z2, то для граничных условий задачи получа­ ем выражения

 

Pi =Poi — 4\Z\ при * = 0;

(4.4)

 

P2 = PiF(Q2)

 

(4.5)

 

v - ^ - = QiPi- Q 2P2 }при

Х = 1и

(4.6)

 

р2 = Р02 + Я2%2 ПРИ X= /, + /2,

(4.7)

где р1

и р2— полные давления во входном и выходном трубо­

Poi

проводах;

 

 

и ро2 — постоянные составляющие тех же давлений;

Qi

и Q2 — объемные расходы;

объемные

скорости

Я\ и <7г — избыточные колебательные

 

(расходы).

и граничные условия

Уравнения гидродинамики (4.1) — (4.3)

(4.4) — (4.7) нелинейны. Однако при рассмотрении устойчивости

и автоколебаний уравнения движения (4.1) — (4.3)

можно взять

в линеаризованном виде.

 

 

 

 

 

 

Полагая, как обычно,

 

 

 

 

 

 

 

 

 

 

 

 

Pi = Poi + Pil

 

 

 

 

 

 

 

 

 

 

Ut= U0l + U',;

 

 

 

(4.8)

 

 

 

 

 

 

Pi =

Poi + pi;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qi = Qoi + <7i

 

 

 

 

(здесь p'. , p /,

UI и q[

— малые приращения соответствующих

величин)

 

и подставляя эти значения в уравнения

(4.1)— (4.3),

после очевидных преобразований получаем

дифференциальные

уравнения в приращениях:

 

 

 

 

 

 

д%_

 

dqi

,

ОП

62qi ,

О„11

6qi

 

П2 \ д2Я\

(4.9)

■26,

dt

+

2 UQi-

dtdx

2a^ 0i

дх

= ( c l - u U

dx2

dt2

 

 

 

 

 

 

dqi

+ 26,91-

U2 —с2

9qi

 

dPi

(i =

1.2), (4.10)

u 0i

ci

 

dt

 

 

 

U,i

 

dx

PoiU0i dt

 

 

 

126

где <7t и pi — приращения объемной скорости .(<?» = Sif/») и дав­ ления (штрихи у pi и U{ отброшены); U0i и ро,-— постоянные со­

ставляющие

скорости

и плотности; 26,= а, + -^ -.

Индекс i

 

 

 

 

 

<%i

 

указывает номер трубопровода (1 — ввода, 2 — выхода).

Введем в уравнениях (4.9) и (4.10) безразмерное время и

координату, полагая

 

 

 

 

 

т =

yi = -2-\

/ = /, +

/»; c2 = V c l - U h .

(4.11)

Производя преобразования координат в уравнениях (4.9) и

(4.10), получаем

 

 

 

 

 

d2q.

d2q.

 

d2q,

dt>i ,

de>t

 

 

+ — 11----2a,— -4i-

(4Л2)

dy2

dx2

 

dy dx "

Ц| {~5Г+а'~£г) ;

dqt

i

d4(

Si

dp{

( i = l , 2),

(4.13)

- r 1— M

i -----

dy

PoAii

dx

dx

a

 

 

где

 

 

 

 

 

 

 

 

ci

 

ai— Ui

 

 

 

 

c.

 

 

Из уравнения (4.12) определяется объемная скорость <?,• как функция времени т и координаты у , т. е. qi — qi{x, у ) . Давление P i ( т, у ) найдется из уравнения (4.13) по известному qi(x, у). В самом деле, из уравнения (4.13) следует

о

Преобразуем также граничные условия (4.4) — (4.7), выразив входящие в них переменные через приращения. Не следует, од­ нако, думать, что при таком преобразовании граничные условия во всех случаях будут рассматриваться как линейные. При ис­ следовании устойчивости и условий самовозбуждения граничные условия будут приниматься как линейные с постоянными коэф­ фициентами. При исследовании автоколебаний необходимо учи­ тывать существенные нелинейности в характеристике компрес­ сора. Этот учет можно осуществить тем, что коэффициенты линеаризованных уравнений (граничных условий) будем прини­ мать зависящими от основных переменных.

• Выражения (4.4) и (4.7) могут быть записаны в виде

P i --P o i = Z’ \ q \ , /?2

Р02 = ^2?2>

причем в дальнейшем разности р\— /?ш и р2— р<я, представляю­ щие собой избыточные давления, будем обозначать соответст-

127

венно р\ и р2. При линеаризации уравнений компрессора (4.5)

и (4.6) воспользуемся соотношениями (4.8). Находим

Рог + Рг = (ро\ + pi) F(Qq2 + 92) —

 

= (ро\ + pi) ^ (Q o2) +

Я2 +

••• ;

(4.15)

У ——(рог + Р2) + (Q01 + <7 i)(P o i +

P i ) — (Q02 + ?г)(Ро2 + р г ).

a t

 

 

 

 

Производя вычисления и ограничиваясь первыми степенями

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

 

 

 

Р02 + Р2 = Poi^(Qo2) + Р01

dQ2Я2 + F { Q

02) p i- ,

dH

 

 

 

(4.16)

= Q oiPoi---Q02P02 +

P01<7l + Q o iP l--- Q02P2--- P02^2-

dt

 

 

 

 

Поскольку для установившегося режима справедливы соот­

ношения

 

 

 

 

Р02 = PoiF(Q 02y,

(4.17)

 

Q01P01 — Q02P02 =

О,

то вместо выражений (4.16) находим (штрихи у р и р опускаем)

 

 

 

AF

(4-18)

 

 

 

Р2 = Poi-rp— Ч2 + F(Q02)pi',

 

 

 

0Q2

 

 

 

У ~ ~

= Poi<7l + Q o iP l---Q02P2--- Р02<7 2 -

(4 . 1 9 )

 

 

a t

 

 

Исключим из уравнения (4.19) величину р, воспользовавшись

уравнением состояния (4.3), которое в линеаризованном

виде

будет

 

 

 

(4.20)

 

 

 

Pi = CiPh

где с

 

Рос

—скорость звука в соответствующей

среде.

у —- -

 

 

рог

 

 

Следовательно, уравнение (4.19) с учетом выражения (4.20)

принимает вид

Q01

 

V

dpi

(4.21)

 

2

dt

Р01Р1 + “ ~ P l — P02P2

с

Cl

 

2

 

 

Представим уравнения (4.18) и (4.21) в виде уравнений че­ тырехполюсника. Для этого1решим их относительно переменных

Р \ и q2. Находим

 

_ 1

Р 01

Р01

 

 

Pi = (

Я1 +

 

Ai

Р 02

Р о 2

 

+ 4 -

•—

Г 1 + f ^

+ У Р о г

/) рпр/]w

(4 •22)

Ai

Рог

L

\

 

J

 

128

 

 

Pol

<7i

£0. Sp2,

 

(4.23)

где

 

Po2

 

*1

 

 

 

 

 

 

 

 

 

 

Co-

A, = l

Q02P01

PoiF'\

F'-

dF

-

d

------;

s

P02c2

 

VP02P02

 

 

dQ2

 

dx

Следовательно, граничные условия вместо выражений (4.4) — (4.7) примут вид

Pi — — Zi4i при t/ = 0;

Pi = h\\q\ + hi2p2

q2= h2iq1+ h22p2 J при t/ = -y -; p2 = Z2q2 при y = 1,

где

Лц

Л12 = —— • P o l

Ai Рог

1

п2\----

Ai

1- . l o L ' - P o L p 0lF'.t

 

А]1

so2

Рог

 

 

Г 1 +

f CoS -Ь

PO!/7' ] ;

L

\

 

VP02 /

J

Pol

.

U _

c 0

~

P02

»

n22

A

 

 

 

Л|

 

(4.24)

(4.25)

(4.26)

(4.27)

i

(4.28)

Уравнения (4.25) и (4.26) являются эквивалентными лине­ аризованными уравнениями (4.5) и (4.6) компрессора. При этом компрессор рассматривается как четырехполюсник, свойства ко­ торого характеризуются четырьмя коэффициентами Ли, h\2, Л21, Л22. При исследовании устойчивости системы будем полагать, как указано выше, эти коэффициенты постоянными,

дF

что соответствует постоянству производной F ' = —— . Посколь- dQi

ку функция F' зависит от расхода Q2, то при нарушении усло­ вий устойчивости уже нельзя полагать F'(Q2) = const. Установ­ ление автоколебаний в системе как раз и обусловлено зависи­ мостью производной F' от расхода Q2. Следовательно, предпо­ ложение о постоянстве коэффициентов четырехполюсника Лц, Л12, Л21, Л22 справедливо до тех пор, пока система устойчива.

Таким образом, поведение системы описывается линеаризо­ ванными уравнениями в частных производных (4.12) и (4.13) при линейных (постоянные Лц, h\2, h2U Л2г) или'нелинейных (пере­ менные Лц, Л12, Л21, Л2г) граничных условиях (4.24) — (4.26).

4.2. ПОСТРОЕНИЕ РЕШЕНИЯ

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

9 Заказ 1516

129

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