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

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

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

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

Добавлен: 19.10.2024

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

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

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

намическим поглотителем (РАП). Во многих случаях в конечных выражениях можно считать приближенно равными мгновенные массовые расходы при входе и выходе из компрессора и поэто­ му пренебречь их разницей.

Поэтому, при достаточо малых числах М газодинамические процессы в компрессоре относительно малых отклонений полных параметров можно описать уравнением

l U ^ - a X p

W6MK= (Я* + а*)бР; _ ( 1 - а ) 8 р к,

(6.28)

at

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

дп*/дМ пр

 

_

Р - Н

, R — 1 .

 

 

 

1 + я*/<рМпр ’

Ф _

2R

'

2R

'

 

 

дМ В пр

 

 

 

 

 

 

 

 

f

МВ пр

 

 

 

 

 

 

 

 

дп В пр

я* = const

— безразмерный коэффициент,

 

 

ПВ пр

 

 

 

определяющий

смещение ха­

 

 

 

 

 

 

 

 

рактеристик

при

изменении

 

 

 

 

пщ>;

 

 

 

 

П пред —

j

 

М'1 +я*)

 

 

~Z~Z

>

эф — '

 

 

 

 

n(X,)sMnp

 

 

SB6(*«) + 5ке(^к)

 

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

акМ к

6рк + — —

о м к

акМ к6М к

 

ьР:

SK

 

(6.29)

 

е(Хк)

 

е (Х к )

 

 

Уравнение трубопровода получится из уравнения (6.28), если

положить я* = 1, а* =

0 и отбросить знак осреднения. При этом

для трубопровода, расположенного перед компрессором, будет

^ . з Ф- ^ = б р ; - б Р; .

(б.зо)

a t

 

Из трех уравнений (6.28), (6.29) и (6.30),

исключая брк и

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

дкМк

« Х р е д - О - 'О ' ®(^к)*к

6МК 8p i.

я* + а*

 

200


Имея в виду, что возбуждение компрессора начинается в окрестностях а* = 0, можно полагать, что инерционность столба воздуха в трубопроводе и компрессоре при изменении наклона напорной характеристики компрессора остается постоянной, т. е. можно принять

 

г*

,

/

,

г

 

 

''к.эф _ , •

ьк.эф

(6.31)

 

 

 

, Г — ^тр.эф ~\

—-

Обозначая

 

 

а* + я *

 

я *

 

 

 

 

 

 

 

L x— ^тр.эф

^•Кэф

 

а’ ^пред-О

t(h)sK

(6.32)

 

F' =

 

я * + а*

получаем для случая малого значения средней осевой скорости газа (Мк «С 1 ) уравнение трубопровода и компрессора в виде

L ] W}L + F'8MK= 6p\.

(6.33)

dt

 

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

длина тру­

бопровода и компрессора не превышает величины порядка чет­ верти длины волны, соответствующей частоте возбуждения. Бу­ дем считать, что дроссель работает при сверхкритическом пере­ паде давления, т. е. 6Л4др = 0.

Линеаризуя характеристику сопротивления R2, найдем

A Pr = R'Mr,

где R' = —-— , а г — коэффициент, численно равный коэффици- siPi

енту трения пористых фрикционных материалов и измеряемый в механических омах на единицу площади.

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

Схрх= — 6МК— 6М*; LXMK+ F8MK= 8px\

 

L2Mr = 8pR— 8р2, С2р2=

(6•34>

R'MR = 8p2— 8p.

 

Уравнения (6.34) с достаточным приближением описывают движение в системе лишь в тех случаях, когда максимальный линейный размер газосборника и РАП вместе взятых значитель­ но меньше четверти длины волны колебаний системы:

А

я а

V l xc x.

4

т

201


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

диненная

емкость

выполнена в

виде

кольцевого

объема (рис.

6.7),

на

котором 1 — дроссель,

2 — сетка, 3 — компрессор, 4

основной

объем

(газосборник),

5 — присоединенный объем. Полученная система четырех уравнений 1-го порядка сводит­

ся к линейному дифференциальному уравнению 4-го порядка

c,c2l,l2

+ c,c2(PL2+ RLx) -^ ~ +

 

 

dt*

at3

 

+ c2[l2+ l3(1 + A .) + c.tfp]

+

 

+ C2 | > + F '( l + - | - ) ] - ^ - + A4 = 0,

(6.35)

где Ci и C2 — параметры, характеризующие емкостные свойства соответственно газосборника и присоединенного объема

Cf = 4а2-(‘ = 1 . 2);

L2— параметр, характеризующий инерционность столба воз­ духа в горле РАП.

Согласно критерию устойчивости Льенара — Шипара, усло­ вия устойчивости компрессора в линейном приближении записы­ ваются в виде

FL + R > 0;

(6.36)

(А + С + 1) + RF > 0;

(6.37)

L(1 + C)RF3 + [R3(L + C + 1) + L 2]F2 +

 

+ R{R* + [I — (1 + C )]2 + 21} F + R2 > 0,

(6.38)

— безразмерные параметры.

Невыполнение хотя бы одного из неравенств (6.36) — (6.38) свидетельствует о неустойчивости работы компрессора. Из этих неравенств следует, что в принятых предположениях работа на

202


«*

дя*

правой ветви характеристики компрессора, где

------- < О,

 

<?МПр

обеспечивает устойчивость работы компрессора «в малом». Кро­ ме этого, устойчивой является и часть левой ветви. Диапазон устойчивой работы компрессора в сети с применением РАП ши­ ре, чем без него. При отсутствии РАП условие устойчивости за­ писалось бы в виде F > 0.

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

Условия (6.36) — (6.38) также показывают, что для практиче­

ски интересных значений параметров L и С уменьшение массо­ вого расхода через компрессор приводит к потере устойчивости при нарушении условия (6.38). Задачу можно свести к опреде­

лению такой совокупности_значений R0пт, LonT и С0Пт, при кото­

ром функция F = F(R, L, С), заданная в виде выражения (6.38) достигает минимума.

Поскольку уравнение

Ф = 1(1 + C)RF* + [R2( \ + с + 1 ) + L2] F2 + R[R2+ ( I — 1 — C)2 +

+ 2L]F + R2 = 0

(6.39)

имеет 3-ю степень относительно F, решение его

затрудне­

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

Будем считать С — параметром. Тогда

экстремум

функции

Ф(Р, R, L)

по двум независимым переменным может быть в точ­

ке, где

 

 

 

 

 

 

 

 

 

Ф = 0; ^ -

= 0; - ^

= 0,

(6.40)

 

 

d R

 

d L

 

 

 

НО

 

 

 

 

 

 

 

 

d F _ — З Ф / d L .

d F _ — д Ф / d R

 

 

d L

д Ф / d F

d R

d 0 / d F

 

Проделав необходимые выкладки, из второго уравнения сис­

темы (6.40)

получаем

 

 

 

 

 

 

F = ----- i=-[tf2 + Z 2- ( l + C

) 2],

(6.41)

 

 

R L

 

 

 

 

 

d F

а подставляя выражение (6.41) в уравнение —=- = 0, находим

(/L

KL [R2+ Г2— (1 + С)21 = R/VZ2— Я2( 1 + С—Z). (6 •42)

203