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