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

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

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

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

Добавлен: 19.10.2024

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

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

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

Далее автор пишет: «Условия существования автоколебаний,

как легко видеть, сводятся к неравенству

> 0. Другими сло-

 

6Q

вами, если рабочая точка на характеристике компрессора выб­

рана так, что

> 0 , то в системе возникнут

автоколебания,

 

3Q

 

о свойствах которых можно судить по предельным циклам».

Переходя к рассмотрению этой интересной

работы, следует

отметить, что в ней, к сожалению, допущены некоторые погреш­ ности. Рассмотрим уравнение (0.16). Из него следует, что если

/'

дР

 

то в системе будут самовозбуждаться коле-

----------------- > 0,

 

dQ

d.Q

 

 

из уравнения (0.1 б) усло­

бания, т. е. возникнет помпаж. Итак,

вие помпажа можно записать в виде

 

 

 

 

dF

3R

 

 

 

 

dQ

3Q

'

 

Однако из рис.

0,4 видно,

что это

условие соответствует не

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

3F

баний является условие ----- > 0, не вытекает из его уравнений.

3Q

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

В настоящей монографии излагается разработанная автором7 теория помпажа, учитывающая нелинейные особенности харак- ( теристик вентилятора и свойства присоединенной к нему сети.Э Эта теория позволяет объяснить явления, происходящие при помпаже, указать причины различий в характере помпажных колебаний у разных типов вентиляторов и оценить влияние раз­ личных факторов на область помпажных режимов. Она доказы­ вает также, что помпаж может происходить не только на режи­ мах, соответствующих восходящим ветвям характеристик, но и на спадающих ветвях, на которых режимы раньше считались ус-,-, тойчивыми. При этом показано, что возникающие на нисходящих ' ветвях колебания обычно более опасньи lrep^bie н^кл^ку-течт^/ю

2 Заказ 1516 I мц V :>о- •'•ХНИЧ Э С К а Я 1 1 7

о -- ' о

С С С Р \


амплитуду и появляются неожиданно (так называемый жесткий режим возникновения колебаний), в то время как вторые обычно имеют мягкий режим возбуждения — возникают постепенно, начиная с очень малой амплитуды. Теоретически определена 31висимость областей устойчивости от внешних условий. Указэд также метод построения кривых, определяющих помпажные ко­ лебания непосредственно по характеристикам компрессора и сети. Разработанная теория позволила найти способ подавления помпажа путем введения системы регулирования с обратными связями, воздействующими на положение дросселя или напразляющего аппарата. Результаты расчетов по этой теории хорошо подтверждаются экспериментами.

РАЗДЕЛ I

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

ГЛАВА 1

УРАВНЕНИЯ ДВИЖЕНИЯ МАЛОНАПОРНОГО КОМПРЕССОРА И ИХ ИССЛЕДОВАНИЕ

1.1. ВЫВОД УРАВНЕНИЙ ДВИЖЕНИЯ

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

обозначим через р„а, рк, ск.

вентилятора дается

Пусть относительная характеристика

зависимостью

 

п ( Q k)>

(1. 1)

где ри — абсолютное полное давление перед вентилятором; <3к — объемный расход за вентилятором.

При установившемся режиме, если пренебречь потерями на трение в напорном трубопроводе, давление Рб перед дросселем равно давлению рк, а объемный расход Q„ за вентилятором ра­ вен объемному расходу QR через дроссель. При неустановившемся режиме рб не равно рк\также не равны расходы QKи QH.

Будем считать, что характеристика сети дается зависимостью

Рба Ро ~ 4*(Qk)‘

( 1. 2)

При составлении уравнений в этой главе сделаны следующие

основные допущения:

система, образо­

1. Сложная распределенная акустическая

ванная вентилятором с присоединенными

трубопроводами и

дросселем, может быть заменена системой с одной степенью сво­

боды. Это — единственное допущение,

имеющее существенное

значение при уточнении излагаемой теории.

2. Степень повышения давления за

вентилятором, а также

изменения давления при колебательном

режиме невелики, так

что можно считать рк р0, Тк Т0, /к ~ /о- Такие компрессоры будем называть малонапорными.

Значения этих параметров берутся соответствующими равно­ весному режиму. При данном предположении объемные расходы до и за вентилятором равны.

3. Скорости потока малы.

2?

19



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

Считая перепад давления Ар — ро рia малым, пренебрегая сжатием воздуха во всасывающем трубопроводе и полагая по­ этому ро = const, получаем по теореме об изменении количества

движения

 

(mo) = F,

dt

у '

где т = poSiA — масса воздуха во всасывающем трубопроводе;

о = Qo — скорость движения воздуха; s'

F — сила, равная произведению разности давлений на площадь сечения трубы: F = Si(po— Ри)-

Отсюда

Р(М — =(Ро— Р1а)«1

s‘

или

Q0-H°^i- = p0—р1а.

«I

Обозначим

 

 

 

РоЛ

г

 

 

 

 

 

------ — Ьа1.

 

 

 

 

*i

 

 

 

 

Величина La = -^ -

называется

акустической массой

трубо-

 

провода.

5

 

 

 

 

 

 

QK, получаем

 

 

Тогда, полагая Qo =

 

 

 

 

Pla = РО— АпОк-

 

 

Рка = Plan{Qn) = (Ро

^\Qidn(Qtd‘

 

1

Рассмотрим нагнетательную часть тракта. Схематизируя яв­

 

ление, полагаем, что на конце трубопровода имеется воздушная

 

емкость 4 (см. рис.

0.2). Тогда в процессе колебания воздушный

 

объем в нагнетательном трубопроводе 3 будет играть роль акус­

 

тической массы La2Воздушный объем в емкости 4 будет обла­

 

дать некоторой акустической гибкостью Са, от которой зависит

 

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

 

лением в ней или расходованием из нее воздуха.

 

 

Подсчитаем давление рба в емкости 4. Оно будет равно дав­

^

лению за компрессором, измененному на величину, требуемую

для ускорения потока воздуха:

 

 

 

 

Рба ” Рка

J-'aiQ.K= (Ро

^al Qk)^(Qk) ^-агОк

(1.3)

20


или

Рба = Роя (Фк) — Qк I^a2 + Laiя(<2к)]-

Обозначим

Po^(Qк) = F, (QK), La2 + LaI:t(QK) = La.

Тогда уравнение (1.3) примет вид

L a Q a ~ Л (< Э к ) Рба

(1.4)

или

 

^а0к= [^1(Ск)—Ро]~ 1Рба~Рб]-

 

Обозначим

 

F(Qh) Po= / ?(Qk); Рба—Ро= Рб-

 

Следовательно, получим

 

LaQK= F(QK)— p6-

(1.5)

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

Полагаем процесс в объеме V адиабатическим. Связь между плотностью и давлением при этом выражается соотношением

-5— = const.

Рб

Дифференцируя это выражение, получаем

dp __ dpa

РРб

Масса dm воздуха, втекающего в объем V за время dt, про­ порциональна разности секундных расходов:

dm = p(QK— QH)dt.

Плотность

или

Следовательно,

V

21

Учитывая уравнение состояния

-*- = gRT,

р

получим

yp{QK—Qfj) & _

dpts

рК

pgRT

Наконец, имея в виду,

что скорость звука

t= V v g R T ,

найдем

dp6

 

V

 

рс2

Qк— Qr-

 

dt

 

Множитель Са = ——-

называется акустической гибкостью,

рс2 .

 

 

Итак, второе уравнение будет иметь вид

 

С - ^ - = QK- Q « .

(i.e i

 

at

 

Третьим будет уравнение характеристики сети (1.2):

p6 = <v(Qr)-

Исключая из уравнений (1.6) и (1.2) расход QK, получаем систему двух нелинейных дифференциальных уравнений 1-го по­ рядка, описывающую движения в рассматриваемой системе:

L^ = F(QK) -P 6 , (1.7)

И)

С а - ^ = 0 к -ф 1 Ы ,

( 1.8)

at

 

где ф1 — обращение функции (1.2).

Разделив уравнения (1.6) и (1.5) почленно одно на другое, получаем дифференциальное уравнение интегральных кривых

dpe __

Qk Qr

/ 1 g\

dQK

f (Qk) Рб '

Ca '

Рассмотрим полученные выражения.

Сделав некоторые допущения, приведем систему (А) двух уравнений 1-го порядка к одному уравнению 2-го порядка. Его анализ дает ряд приближенных выводов о мягком и жестком возбуждении колебаний, об областях устойчивости и другие ре­ зультаты. В главе 2 дан геометрический метод точного интегри­ рования выведенных уравнений движения и при помощи этого метода исследованы некоторые системы *.

* Под точным понимается метод интегрирования уравнений (1.5)— (1.8) без каких-либо упрощений. Разумеется, любое численное или графическое ре­ шение обладает погрешностями расчета.

22