Файл: Казакевич, В. В. Автоколебания (помпаж) в компрессорах.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 + |
2е |
(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