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

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

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

 

Длина

отдельных

элементов

 

берется

такоц,

чтобы

внутри

 

каждого

 

из

них

величина

 

F(Qk) — Рб

или

QK— фДРб) из­

 

менялась незначительно.

Нужно

 

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

 

жением (2 .6), при помощи зара­

 

нее найденного предельного цик­

 

ла легко построить кривые изме­

 

нения рк или QKпо времени.

 

2.3.

ФИЗИЧЕСКИЙ СМЫСЛ

 

ФАЗОВОЙ ПЛОСКОСТИ

 

 

Пусть фазовая диаграмма по­

 

строена. Это значит, что в коор­

динатной системе QK, Рб нанесены характеристики вентилятора

Рк = -F(Qk)

и сети QK= фДРб) и фазовые траектории системы

уравнений

(1.7) и (1.8), в том числе

предельные

циклы

(рис. 2.4).

 

 

 

 

 

 

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

На рис. 2.4 видно, что расход QK через вентилятор меняется при помпаже от Qc, до QCj на величину AQ„ = Qc, — Qc,- В то же время разница расхода через дроссель менее значитель­ на, она характеризуется величиной AQH= Qo> — Qo, .

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

2.4. ИССЛЕДОВАНИЕ ПОМПАЖ А В КОНКРЕТНОЙ СИСТЕМЕ. ВЛИЯНИЕ НА ХАРАКТЕР. ПОМПАЖ А ИЗМЕНЕНИЯ ПОЛОЖЕНИЯ

РАБОЧЕЙ ТОЧКИ НА ХАРАКТЕРИСТИКЕ ВЕНТИЛЯТОРА

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

70

рических

размеров

сети

(параметра

 

-^ -)

на

структуру

 

 

Са

пло­

фазовой

скости

и проследим,

как

сказывается на ней изме­

нение

положения дрос­

селя.

 

 

вентилятор,

Выберем

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

которо­

го дана на рис. 2.5: про­

изводительность

равна

примерно 1 м3/с, а макси­

мальный

напор — сколо

100 мм вод. ст.

 

 

Будем вести расчет в

системе CGS. На рис. 2.5

видно,

что

масштаб

по

оси

ординат

взят

в

100 раз большим, чем по

 

 

 

 

 

оси абсцисс (цена деле­

мм вод. ст.== 1,0 г/см2 = 103 дин/см2;

ния по оси ординат —

100

по оси абсцисс —

0,1

м3/с

=

105 см3/с, следовательно, масштаб-

 

 

 

1Q5

 

 

 

Как было показано ранее,

ный коэффициент т = —— = 100.

£а = Ро(— + - ^ М

; с а

о15/2^2

Рос2

 

 

 

 

 

h

) '

° а

Примем площади входного и напорного трубопроводов

одинаковыми и равными

 

 

 

 

 

S1 = s2 = s = 0,5

м2 = 0,5. 104 см2

Тогда

 

 

 

 

 

 

 

 

г

—*

Ро

'(/.! + 0„5/2);

с а =

0,5/feS

 

 

5

 

 

 

 

 

Рос2

 

и

_

1<?

, [t

+ 0,5V

V

 

Са

 

 

s2

 

\

Имеем

 

 

 

г

v

 

 

 

 

 

 

 

 

ро = 1,22 • 10_3 г/см3;

с 340 м/с = 3,4 • 104

Если положить = 25 см, k — 100 см, то

,

1 ,22-ю- 3

(25 + 50)^ 0,182-Ю -4;

 

0 ,5 -104

 

 

 

 

 

0,5-100.0,5-104

0,18;

1,22,10_ 3 .3,42.109

 

71

к

2-1,222-10~6-3

,4 2-108

,

q 5 \

с,

0 ,2 5 -108

 

\ 12

, Ь ) ~

= 1,34.10-“

----0,5) = 1,34-10-“ ( - ^ - + 0,5) ^ 10-“.

То же самое получаем, естественно, поделив непосредственно подсчитанные значения La и Са.

Подсчитаем период малых колебаний в окрестностях равно- df

весного режима в случае, к о г д а ------ .При этом по формуле

(1.26) имеем

dQK

 

 

 

1

---------:--------т ~ з - 105

1/с2;

со2

CnZ-я

0,18-0,18-10-“

 

со = 5,5-102

1/с; Г = — = 1,ЫО-

2 с.

 

 

Таким образом, период помпажных колебаний в момент воз­ никновения помпажа равен « 0,01 с.

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

положениях дросселя. Подсчитаем параметру =

mPL

------- , который

 

____

показывает, во сколько раз необходимо изменять отрезок М\А\,

или М\А%. Подставляя значения т и получаем Са

v = m2-^ - = 1002 - 10—“ = 1 .

Са

Следовательно, в данном случае длину отрезков М\А\ и М\А% изменять не следует.

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

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

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

72

Рис. 2.6 Рис. 2.7

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

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

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

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

Легко видеть, что стабильность помпажа зависит от расстоя­ ния между устойчивым и неустойчивым предельными циклами: чем это расстояние меньше, тем меньше стабильность помпажа. В то же время расстояние между циклами тем меньше, чем бли­ же положение дросселя к тому, которое соответствует полуустойчивому предельному циклу. Этим объясняется известный из экс­

73

периментов факт, что в случае жесткого режима помпаж в мо­ мент возникновения является малоустойчивым и то возникает, то исчезает.

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

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

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

Рис. 2.8

Рис. 2.9

74

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

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

На рис. 2.10 показана фазовая плоскость в слу­ чае полного закрытия дросселя (что соответст­ вует, например, напол­ нению бака). На ней имеются два предельных цикла (из которых мень­ ший — неустойчивый) и

устойчивый фокус. Следовательно, в системе возможен жесткий помпаж.

Чтобы в случае расположения рабочей точки слева от мини­ мума характеристики вентилятора получить такой же характер фазовой плоскости, как на рис. 2.5, необходимо рабочую точку характеристики расположить еще левее, чем указано на рис. 2 .10. Для этого требуется принудительно вдувать в вентилятор воздух через напорный трубопровод. В системе тогда будет существо­ вать особая точка типа устойчивого фокуса или узла, а периоди­ ческие движения исчезнут, так как неустойчивый и устойчивый предельные циклы сольются, превратятся в полуустойчивый пре­ дельный цикл и затем исчезнут.

2 .5 . ОПРЕДЕЛЕНИЕ ПЕРИОДА КОЛЕБАНИЙ

Определим период колебаний для одного из рассмотренных случаев, например, для показанного на рис. 2.8. Точками Nu N2, N3 и N4 разбиваем предельный цикл на четыре дуги NiN2, N2N3,

N3N4, N4Ni. Дуги NiN2 и N3N4 делим на участки,

выделяемые

полосами шириной в одно деление оси абсцисс, а

дуги N2N3 и

N^Ni — Ha участки высотой в одно деление оси ординат. Концы участков отмечаем цифрами /, 2, ..., 31. От середины каждого

из участков дуг N\N2 и N3N4 проводим вертикали СД)* до пере­ сечения с характеристикой вентилятора (на рисунке отмечена вертикаль С7Д7). От середины каждого из участков дуг N2N3 и

75

NaNi проводим горизонтали С,£,- до пересечения с характеристи­ кой сети (на рисунке отмечена горизонталь С22Е2э).

Подсчитаем длины отрезков С,£г- и С,£),- и образуем выраже­

ние

 

 

 

 

 

 

 

 

 

 

 

23

 

 

 

 

 

 

 

 

 

 

 

T = L.

 

 

Dfii

 

 

Dfii

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

 

31

 

 

 

 

 

 

 

 

 

 

 

+ C.

 

V (A2B2) t

i

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ - 9

С,Et

+

Z j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ 2 4

 

 

 

 

 

 

 

Зная, что одно деление оси абсцисс соответствует 105 см3/с, а

одно

 

деление

оси ординат

103 дин/см2,I что

La = 0,18*10-4 и

Са =

 

0,18, получаем

Г

 

 

 

 

 

 

 

 

 

 

 

 

Т = 0,18-10-4

105

1

. — +

 

+ — + — + — + —

 

 

 

 

 

 

 

103 [

5,3

6

 

6,2

4

2

 

2

2,2 +

 

 

 

 

'

I

'

 

1

,

1

.

1 .

1

 

,

1

пт]

+

 

2

 

 

 

 

8,6

 

 

 

 

 

 

 

 

5,8

 

7,2

 

8,3

 

 

 

3.8

3,5

 

 

 

 

 

 

 

 

 

 

+ 0,18

I03

 

1

 

1

 

1

 

1

.

1

.

1

\-— +

 

105

 

3,5

 

4.1

 

4,8

 

5,1

 

5,2

 

5,3

 

 

 

1

1

 

 

 

 

 

'

5

 

 

 

 

 

 

1 +

1

.

1 .

1

,

!

,

0 ,7

] ‘

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3,3

 

4,5

 

5,2

5,6

 

5,8

5,6

 

5,5

 

5

 

 

 

Т = 0,18-10—2 -4,32 -4-0,18-10—2-2,93 =

1,3 •10~2,

 

 

т. е. период Т =

0,013 с.

 

 

 

 

 

по формуле

(1.26)

Ранее период колебаний был определен

без учета сил трения. При этом было получено Т = 0,011 с. Таким образом, в рассмотренном случае под влиянием сил трения зна­ чение периода увеличивается на 0,002 с, что составляет 15% от истинного значения Т.

2.6. ВЛИЯНИЕ ИЗМЕНЕНИЯ ГЕОМЕТРИЧЕСКИХ РАЗМЕРОВ СИСТЕМЫ НА ХАРАКТЕР ПОМЛАЖА

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

Из выражения

Са

Ц

/

U

следует, что —— увеличивается при уменьшении s и увеличении

I

 

— . Если s возрастает, то ——убывает неограниченно, если же

/2

са

76

убывает до нуля, то ——

12

 

 

Са

стремится к пределу

и

—— =

2

2

 

Са

~

 

Рос

 

 

= — -—

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

S ‘

 

 

вса­

что при любой длине

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

сора, при уменьшении которого помпаж в системе невозможен.

Для полноты картины рассмотрим влияние изменения

от

 

Са

О до оо. Положим вначале, что — — >0 (рис. 2 .1 1 ). В этом слу-

Са

чае точка Л3 (рис. 2 .12 ) неограниченно приближается к точке

Mi, а точка Л4 к точке Аг. Из этого,

а также из уравнения

(1.9)

следует, что если F(QK) — ре ¥= 0,

т о - ^ — »0, т. е. все

инте-

dQK

тральные кривые будут идти горизонтально, и фазовая диаграм­ ма примет вид, показанный на рис. 2.11. Если точка пересечения характеристики сети находится на восходящем участке характе­ ристики вентилятора, то независимо от ее расположения на этом участке образуется устойчивый предельный цикл, составленный из дуг АВ и CD характеристики вентилятора и отрезков прямых

Рис. 2.12

9

AD и ВС. При этом особая точка О является неустойчивым узлом.

Рассмотрим теперь случаи, когда —— Ф 0 .

 

.

 

 

dF

 

Са

 

 

dF

Определяя

 

 

и

,

получаем

графически

------

k,

=0, 01;

k = 0,05.

 

 

 

dQ’K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

вначале случай,

при

котором т2——= 0,2, т. е.

-г—= 0,2-10-4.

Этого значения можно достигнуть,

например,

са

 

 

 

 

 

прежними, увеличить их

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

диаметр в 1,5 раза.

 

 

 

 

 

 

 

 

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

находим

 

 

 

 

 

 

kCa

 

 

 

 

 

 

 

 

0,2. 101'-4

0,0004.

 

 

 

 

kCa

 

0,05

 

 

 

Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,05

т. е.

k яг 11

 

 

 

 

 

 

V b

 

 

V 0 , 2 . 1 0 ~ 4

 

 

 

 

F0

 

0 , 0 1

 

 

 

 

 

 

=

= r l / - ^ - = 2 .2 l / -F~-

 

 

 

 

1

с а

 

У

Ca

 

 

] /

 

[Л— 4

У

 

 

 

0,2. 10'

 

 

 

 

 

 

На рис. 1.2 видно,

что точка с такими

координатами лежит

над прямой 3 и справа от прямой 6, т. е. находится

в области

неустойчивых узлов. На рис.

2.12

построена

фазовая плоскость

для рассматриваемого случая. Здесь же треугольником АВС, у

Св

масштаб для умножения отрезка

которото —— = 0,2 , показан

/1 О

 

 

 

М\А\ на величину v = m2- ^ - = 0,2 .

Построение

показывает, что

Са

 

 

типа неустойчи­

в системе действительно имеется особая точка

вого узла и один устойчивый предельный цикл.

 

Если теперь положить т2

=

1,то на плоскости параметров

Са

k и F' легко усмотреть, что особая точка будет неустойчивым фо­ кусом.

Выполним теперь аналогичное построение для случая, когда

- ^ - = 2 . 10-"

или т 2-^ - = 2 .

Са

Са

 

 

При Э Т О М

 

2-10,-4

 

 

 

kCa

5

-2 = 0,004.

 

 

- 10-

78

Получить значение -^ - =

2 - 10- “

можно, полагая,

что

Са

 

 

 

 

s = 0,5-10-“;

= 1 .

 

 

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

(см. рис. 1 .2 )

этому случаю соот­

ветствует особая точка типа

неустойчивого

фокуса.

Фазовая

плоскость для рассматриваемого случая показана на рис. 2.13, а.

Так как т2—^- — 2 > ^целесообразно уменьшать отрезок М1А2

с а

(в 2 раза).

Из сравнения рис. 2.13, а с предыдущими видно, что дальней-

U

шее увеличение ----- привело к значительному увеличению ам-

Ся

плитуды колебания давления и к большему уменьшению ампли­

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

;

при т2

3 (рис. 2.13, б) качественная картина сохраняется.

Са

 

 

При -^-/?г2= 4,5 качественная картина изменяется, состояние рав-

Са

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

(рис. 2.14).

 

Будем дальше увеличивать

При этом неустойчивый пре­

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

предельный цикл — уменьшаться. При некотором значении

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

еще больших значениях —— сделаются невозможными периоди­ ка

ческие движения в системе, т. е. помпаж исчезнет (рис. 2.15, 2.16).

Рассмотрим, например, поведение траектории, проходящей

через точку Мх при -^ -= 8 -1 0 -4 (рис. 2.15). Проводя через

Са

точку М1 фазовой плоскости вертикальную и горизонтальную

прямые, получаем отрезки М\А\ и MiA2. Уменьшая отрезок М\А% в 8 раз, замечаем, что точка Л3 оказывается близкой к точке М\, а точка А4— к точке Следовательно, интегральная кривая идет круто вверх, в точке, лежащей справа от точки Аг, пересе­ кает характеристику вентилятора, далее также круто с уклоном влево подымается вверх и вблизи характеристики сети резко изгибается влево. Далее она идет вдоль характеристики сети слева и накручивается на особую точку О.

79

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