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

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

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

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

Добавлен: 19.10.2024

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

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

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

проводе начинает уско­ ренно двигаться в ту или иную сторону.

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

при ср < 0 система устойчива при любых положительных зна­ чениях х.

На рис. 4.3 в тех же координатах нанесены линии постоян­ ных значений декремента затухания (нарастания) системы. От­ метим наиболее характерные линии на этой фигуре. Линия 6 = 0

состоит из отрезка прямой ф = х и дуги АВ гиперболы ф = — .

Линия 6 = сю параллельна оси ординат и отстоит от нее на еди­ ницу. Очевидно, что на линии 6 = оо система абсолютно устойчи­ ва. Эта линия, как легко видеть, соответствует нагружению трубы на волновое сопротивление. Линия 6 = —°° параллельна оси абсцисс и отстоит от нее на единицу. На этой линии система аб­ солютно неустойчива.

Следует отметить

особую роль линий ф =

1 (6 = — оо) и

х = 1 (6 = оо), первая

из которых соответствует

равенству ве­

личины Po\F' волновому сопротивлению Z22, а вторая — равенст­ ву сопротивления нагрузки R2 волновому сопротивлению Z22Эти линии совместно с линиями ф = х и фх = 1 определяют то­ пологическую структуру разбиения области параметров системы.

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

е~16Ра_ £2 = COnSt.

143

Решая уравнение (4.82) относительно <р

(предполагаем, что

ф = 0), получаем

 

±[fe(i+ * ) - ( ! - * ) ]

 

Ф ±[А(1 + х) + (1-к)]

 

или

 

ft(l +х) — (I — н)

 

й(1 +х)+(I—х)

(4.91)

fe(l +х) + (1 — х)

 

ft(l +х) — (1 —У.)

 

Придавая величине k (или, что все равно, величинам декре­ мента б) различные численные значения, получим семейство кри­ вых (см. рис. 4.3), проходящих через точку <р = 1, х = 1.

Рассмотрим более детально влияние активного и реактивно­ го сопротивлений выхода трубопровода на поведение системы. Для этого обратимся к выражению (4.78). Если принять напор и скорость потока в компрессоре малыми и считать, что труба на конце снабжена достаточно широким фланцем, сравнимым по величине с длиной волны, то можно приближенно принять излучение трубы за излучение круглого отверстия в бесконечной стенке [13]. При этом

 

Яг = — 0О;

(4.92)

 

 

■’ др

 

 

 

S 0c

 

 

 

 

Хо,

 

 

 

 

4др

 

 

где

 

 

 

 

 

Л (2/сг) .

Н ! ( 2 к г ) .

>

 

К Г

* АО —

 

 

КГ

 

г — радиус выходного сечения трубы;

\

к — волновое число

/

Шп

с

1к = —£.;

o)p = cl>— J;

с — скорость звука;

 

 

 

 

7i(2/cr) — функция Бесселя первого рода;

 

Hi (2кг) — функция Струве первого рода.

разложения вида

Для функций J(2кг)

и Н\ (2кг)

известны

У, (2кг) = кг-

 

 

 

 

2-6

2* К4г4 +

(4.93)

 

 

 

 

К6Г6— .

32-5

32 -52-7

 

На рис. 4.4 даны графики зависимости функций 0 и хо от ар­

гумента кг. Из этих графиков,

как впрочем и из выражений

144


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

противлению -2^-, а реактивная, sap

&0,+O

 

 

 

 

 

 

/

 

 

 

 

 

 

f r

 

 

---- - Xo__

Оt

V

j

2

 

5

6

8

7r

Рис. 4.4

пройдя

через максимум при кг =

стремится

= 1,25,

затем сильно уменьшается и колебательно

к нулю.

Поэтому,

если кг достаточно велико (кг > а), то мож­

но считать 0о = 1

и хо = 0- Тогда б-> оо. Если sap =

s2, , т. е.

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

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

Р. А. Шиповым [36] получено точное решение для случая круглой полубесконечной трубы с бесконечно тонкими стенками при наличии во всем пространстве равномерного однородного потока газа в направлении оси трубы, с числом Маха равным М.

На графиках, построенных Р. А. Шиповым (рис. 4.5 и 4.6), даны найденные зависимости для различных значений числа Ма­ ха (сплошные кривые) и для случая трубы с фланцем на выхо­ де (пунктирные линии).

Аналитические выражения для трубы без фланца имеют вид

при х

1

 

 

0 о = О,22к2г2;

 

 

Х о = 0,61кг.

(4.94)

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

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

kr < 3,83 / 1 — М2.

После возникновения двумерных (асимметричных, а затем и неасимметричных) волн полученные выражения оказываются не­ справедливыми.

Рассмотрим вопрос о возбуждении колебаний. На границе устойчивости декремент бг равен нулю и

с-дя. / gM+p0,F '\a

(Z22-*2)2 + *2

l

\ z **— PnF'/

(Z22 + R2f + x l

 

JO З а к а з 1516

145


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

zn

V ( z 22- R 2)2+ x 2+ \ f ( z 22 + R2)2 + xl

 

(4.95)

POi

V (Z 22 + ^г)2+ X2 ~ V (Z22~ ^г)2 + x 2

— определяет границу области динамической устойчивости сис­ темы;

j? / ( Q u) _ Z32 ^ ( ^ 2 2 + ^ г ) 2 +

{Z22~^2) + Х2 ^ g g y

Р°1 1/(^22 + Я2)2 + *2 + ^ { Z22~^2f + Х2

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

146

Введем обозначения

Z22

V (Z22~K 2)2+ XI + V (Z22+ F2) 2 + X2

_

p .

Рп

V(z 22+ F 2)2 + x 2- V

 

 

(4 .9 7 )

( z 22- R 2)2 + x l

 

 

Z22

V {Z22 + F ^ f + X\ — Y

{Z22— F ^ f + X\

__ p

 

V( z 22 + *2)2+ х\ + V(z 22“ *2)2+

 

(4 .9 8 )

Р°‘

 

 

С учетом этих обозначений неравенство F' < D

определяет

область динамической устойчивости линеаризованной системы,

т. е. область, где самовозбуждение колебаний невозможно. Нера­

венство же F' < D2 определяет область статической

устойчи­

вости.

и D2 =

Построим теперь графики зависимости D\ — D\ (кг)

= й 2(кг) по формулам Морза.

 

1 0 *

147


При 5др $2

 

 

_ Z2i

V

О + вр)2 +^о + У' (1~^р)2 + -уо

Ро1

Y

(i + e0)2 + Jfo—1^(1—0о)2 + ^о

Пользуясь выражениями для 0о и хо, получаем графики за­

висимости Di=Di(Kr)

и D2 — D2(tcr), изображенные на рис.

4.7, а, б линиями 1.

 

 

Рассмотрим характер возбуждения колебаний при возраста­ нии р' от р ' 0 (при сохранении остальных параметров). Оче­ видно, что (при выполнении условия F' < D) при возрастании F' от F' 0 первой будет возбуждаться та частота [из определя­ емых выражением (4.77)], при которой

F' = Di(Kr)= ^ ~ ~ ) -

Отсюда следует, что первыми возбуждаются колебания более высоких частот, а затем, по мере возрастания F' — и низших по­ рядков при тех значениях F', которые определяют частоты, сле­ дуемые из выражения (4.77). При этом частоты более высоких порядков сохраняются.

Мы пока не рассматриваем вопроса об амплитуде этих коле­ баний, который может быть решен при анализе нелинейных свойств системы. Данный случай — изменение F' при сохране­ нии остальных параметров — практически неосуществим.

Если же возбуждение помпажа вызывается изменением поло­ жения дросселя, например его прикрытием, то одновременно ме­ няются Хо, а0, эквивалентная величина sap и форма кривой D\ = = D\{nr). В связи с этим ранее возникшие колебания различных частот могут становиться неустойчивыми и исчезать по мере при­ крытия дросселя.

л Ро1

148

Здесь необходимо отметить следующее. Из графика Di = = D 2(Kr) следует, что при малых F' не обеспечивается условие статической устойчивости F' < D 2. Указанное обстоятельство вы­ звано тем, что мы не учли потерь на дросселе, вызванных сно­ сящим потоком воздуха. Если охарактеризовать эти потери не­ которым коэффициентом Ro, то примерный вид кривых Ь\{кг) и D 2{nr) будет таким, как показано на рис. 4.7 кривыми 2. При этом указанная трудность отпадает.

Осциллограммы Пфляйдерера и Вейнриха, зафиксировавшие развитие помпажа, позволяют проследить за его характером. Из рис. 0,15 для 12 = 1,5 м видно, что в процессе прикрытия дроссе­ ля в момент начала колебаний при расходе V « 0,59 м3/с начи­ наются высокочастотные колебания. При дальнейшем прикры­ тии дросселя (при возрастающем F') появляются колебания при­ мерно вдвое низкой частоты и большей амплитуды. При продол­ жающемся прикрытии дросселя амплитуда низкочастотных ко­ лебаний уменьшается и они переходят в высокочастотные коле­ бания.

На рис. 0.16, а отчетливо видна форма полигармонических ко­ лебаний с большей амплитудой низкочастотной составляющей.

Интересны осциллограммы, приведенные на рис. 0.16,6, по­ казывающие протекание помпажа при различных, но фиксиро­ ванных числах оборотов и при постепенном прикрытии дроссе­ ля. Для п = 14 300 об/мин, ясно видны высокочастотные колеба­ ния малой амплитуды, которые при закрывании дросселя допол­ няются низкочастотными с постепенно возрастающей амплиту­ дой. При большом прикрытии дросселя остаются низкочастот­ ные колебания, в малой степени искажаемые высокочастотной составляющей.

При п — 16 300 об/мин влияние высших гармоник проявляет­ ся более сильно, а при п = 17220 об/мин они играют основную роль в колебательном процессе.

Предположим теперь, что на входе в компрессор подключено

сопротивление Z b Полагая

в уравнении

(4.5) х\ = 0

(так как

li = 0) и пренебрегая распределенным сопротивлением

(рг = 0),

найдем

 

 

 

(-Z ^ -Z zM O z + Z

, ) ^ (Z22- Z

2)(01+ Z ,)e-*’.

(4.99)

Если положить здесь Z\ = R\ + jXi; Z2 = R2 + jX2 и учесть,

что x2 — jx2— 602, то после преобразований, аналогичных изло­ женным выше, получим выражения:

для декремента

е-4вЭ,

(е,+ я,)2+ х? (z22- r 2)2 + x 2

(4.100)

149