Файл: Гликман Б.Ф. Автоматическое регулирование жидкостных ракетных двигателей.pdf

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

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

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

Добавлен: 11.04.2024

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

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

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

мо изменить исходную систему уравнений гидромеханики. Рас­ смотрим вначале только входной 'коллектор турбины, пред­ полагая, что сечение его постоянно, а расход газа" по длине уменьшается вследствие оттока через сопла, равномерно рас­ пределенные по длине коллектора. Для течения газа с перемен­ ным расходом изменяется первое уравнение системы (4.42) — уравнение неразрывности *, в котором появляется дополнитель­ ный член [24]

dp

I д(рта>)

pi« dG

(4.49)

dt

дх

G

дх

 

 

 

Перед линеаризацией

уравнения

(4.49)

 

удобно произвести сле­

дующую замену:

 

 

 

 

 

 

dG __ _q

 

^nx

 

(4.50)

 

 

 

ЙѵОЛ

 

 

 

 

 

 

где Дол — длина коллектора турбины;

аппарат турбины на еди­

q — расход газа через сопловой

ницу длины коллектора;

согласно принятому предпо­

ложению в статике 9 = const.

Для определения величины q можно воспользоваться обыч­

ным уравнением (4.9),

которое в линеаризованном виде с учетом

зависимости

(4.7) запишется так:

 

 

 

 

 

 

 

 

b q = — — S s - ff l —ат---- j 8/?-|-ат8/?н.

 

(4.51)

Подставив

в уравнение (4.49)

величину âGjdx

из

соотношения

(4.50),

линеаризовав уравнение неразрывности

и

учтя

зависи­

мость (4.51), находим

 

 

 

 

 

 

 

<ЭЗр_,

до (рэд)

1— а„

X — 1

8/J —у о х + а то/7и = 0, (4.52)

dt '

вх

дх

Р^ко

2 х

 

 

 

 

 

 

где Fкол — площадь проходного сечения коллектора. При приве­

дении вариации к ‘безразмерному виду вариацию

скорости

от­

несли к скорости газа

во входном сечении коллектора

швх,

так

как по длине коллектора скорость потока изменяется. Обратная величина коэффициента ql(pFK01l) = GBX/ (lKoapFK0„) =w DJ l li0n=

= 1 /ткол — имеет размерность времени.

Полученное линеаризованное уравнение неразрывности (4.52) заменяет первое уравнение исходной системы (4.42), остальные два уравнения остаются неизменными. Если сделать предполо­

жение, что акустическими эффектами можно

пренебречь,

т. е.

* Сохранение неизменными уравнений импульса и

энергии для

тракта

с отводом массы по длине в нашем случае объясняется тем, что отводимая масса имеет такие же параметры, как и параметры среды в данном сечении тракта.

о;


dop

dbw

 

.. ,n,

исчезает.

что — = ---- = О, то второе уравнение системы

(4.42)

дх

дх

системы

(4.42)]

скорость

В уравнении энергии [3-е уравнение

газа перед производной во втором

члене — величина

перемен­

ная по длине коллектора. При постоянной площади поперечного

сечения коллектора и равномерном

отводе газа по длине ско­

рость газа определяется простым соотношением

 

® = - T 2L^K

=

(4.53)

‘кол

Т-КОЛ

 

Величину Ткол принимаем за условное время пребывания газа в коллекторе (по скорости на входе). Подставив в 3-е уравнение

системы (4.42) соотношение (4.53) и искомое частное

периоди­

ческое решение

8s= 8se'lu*,

находим

уравнение для

амплиту­

ды вариации энтропии:

 

 

 

 

 

 

іш8Іф- Ь р і ~ .х . dbs

— о.

 

. (4.54)

 

 

Т КОЛ

d x

 

 

 

Уравнение (4.54)

имеет решение

 

 

 

 

 

87= С (/кол- х ) іштко,,

^

 

Учтя начальное условие

 

 

 

 

 

 

х = 0 \

8s=

8sKOJI,

 

 

где 8sKOjl — амплитуда вариации энтропии

на входе в коллектор«

находим окончательную форму решения для 8s:

 

8s= 8s

н

=

°5Коле

 

(4.55)

 

 

 

 

 

 

Подставив в уравнение неразрывности (4.52) искомое периодиче­

ское решение для всех переменных Ъх— Ъхр1“1 и проинтегри­ ровав его но координате х, учтя соотношения (4.44) для вариа­ ции плотности и (4.55) для вариации энтропии, после преобразо­ ваний получаем

( - ^ + 1 - « г -

(***кол + 1 ) ^ -

 

— ^/шткол-|——j 8sK0JI — (і“ТКол+ 1) 8^кол= 0»

(4.56)

где бСкол — амплитуда вариации расхода газа на входе в кол­ лектор. При интегрированииучли, что д8р/дх = 0, а вариация расхода газа в выходном сечении коллектора (х = 1) б ( р а у ) В ы х = 0, так как скорость газа у этого сечения надает до нуля. При т К о л = 0 уравнение (4.56) переходит в обычное соотношение (4.48) для вариации расхода газа через сопла турбины.

205


Уравнение (4.56) описывает динамику только одного участка общего газового тракта газогенератора — коллектора турбины. Воспользовавшись принятым условием д8р/дх = 0, объединим уравнение газового тракта самого газогенератора до коллек­ тора (4.48) с уравнением коллектора (4.56), приняв öGKon= 6GT. Подставив соотношения (4.8), (4.16) и (4.17) в объединенное уравнение, находим окончательную форму уравнения газогенера­ тора с коллектором турбины:

“ Іті + тко.-і) (1 + г'1оТкол)+ (

й»хкол) ^ 1

ат-)— — -'j

%— 1

S p = ( l + ^ х кол) (— 2 s — Ю 0 +

80Г) -

2%

 

G0 +

Gr

G0 +

Gr

 

__ е - ,вт. - ( 1 +

icütK0J

(800

— 8C?r).

(4.57)

Ha рис. 4.15 приведены АФХ газогенератора с коллектором тур­ бины при различных соотношениях между Т| и т КОл- Кривые по­ строены по безразмерной частоте Ѳ= соті, причем во всех случаях принято ті+ тК0Л; = const. При малой величине отношения тКОл/ті кривая АФХ газогенератора с коллектором 'проходит близко к кривой для газогенератора 'без коллектора. При соизмеримых ве­ личинах Ткол и ті изменяется общий характер АФХ, а именно: исчезает петля, уменьшается амплитуда. При тКОл/ті = 0,9 общий характер кривой АФХ по своему виду приближается к кривым АФХ апериодического звена 1-го порядка (см. табл. 2.1).

206

Изменение характера кривых объясняется тем, что в уравне­ нии коллектора (4.56) отсутствуют члены с запаздыванием при вариации энтропий, которые имеются в уравнении динамики газо­ вого тракта (4.47) и присутствием которых (см. § 4.4) объясня­ ется специфическая форма кривых АФХ. При ТколМ =0,9 опре-

,деляющую роль играет коллектор. Изменение типа уравнений при переходе от тракта, имеющего сосредоточенный выход, к тракту с распределенным выходом показывает, что в последнем случае роль энтропийных воли уменьшается. Это видно и из срав­ нения решений (4.46) и (4.55). Решение (4.46), описывающее движение по тракту волн энтропии с постоянной скоростью, за­ меняется решением (4.55), в котором 'скорость волны уменьша­ ется при удалении от входа и в конце коллектора все волны эн­ тропии сливаются, так как скорость стремится к нулю. При ин­ тегрировании по длине тракта член, характеризующий энтро­ пийную волну, принимает форму двучлена для апериодического звена 1-го порядка, а не члена с транспортным запаздыванием, как в уравнении (4.47).

При распределенном выходе из газогенератора изменяются также уравнения динамики газового тракта за турбиной (тазовода). В силу разного времени пребывания в коллекторе турбины различных частей одной н той же порции г'аза с одинаковой тем­ пературой волны энтропии на входе в газовод оказываются уже не плоскими, как в газогенераторе, а пространственными. Для упрощения будем считать коллектор одномерным, а газовод (как продолжение коллектора)— имеющим прямоугольное попереч­ ное сечение, причем высота его совпадает с длиной коллектора, а ось перпендикулярна оси коллектора (см. рис. 4.14).

Воспользовавшись зависимостью (4.53), определим время пребывания порции газа в коллекторе от момента поступления

до данного сечения с координатой (по оси) х:

X

dx

 

 

 

Г

: Г _Ек2 ^ _ = _ т;кол1пЛ

 

т = \

 

:)

w

• ) ЧІОЛ X

\

' К О

о

 

 

 

 

где т — запаздывание времени поступления данной порции

газа

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

кол­

лектор, по отношению к моменту поступления всей порции

газа

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

y = {x —xn)wt т + Ткол1п

1 - - І - ) W,

\

'кол 7.

207



где Tn — время пребывания данной порции газа в коллекторе до

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

турбины, имеющего -координату х. Приняв

т=Т2, где Т2 — время

пребывания газа в газоводе, и отнеся все

к длине газовода Д= шгТг, находим кривую профиля волны энтро­ пии в газоводе:

V -)—Пі£і- pi (1 — х), (4.58)

Іг

1 де X х}Ікая.

На рис. 4.16 приведены кривые форм волн энтропии в газоводе, рассчитанные по формуле (4.58) для различных-значений Ткол/тгВсе кривые пересекают ось х при ,т<1, т. е. волны энтропии «растягиваются» по всей длине коллектора, причем к моменту прихода головной части волны к выходу из газовода (у=і\) неко­ торая часть волны не поступила

0,5----------\ — у

 

 

еще в газовод из коллектора. Это

\ --------- естественно, так

как

[см.

(4.53)]

 

 

 

/

\

I

при X—>-1

ку-э-0

и соответственно

W

rr ^ /

 

 

запаздывание т—>-оо.

 

 

 

 

 

 

Выведем

уравнение газовода

W

Т2=1>9'

 

 

с учетом изменения формы энтро­

 

 

пийных волн. Для этого восполь­

'

0,5

 

X

зуемся общей системой линеари­

Рис. 4.16. Форма волн

зованных уравнений

(4.42),

пре­

небрегая акустическими эффекта­

энтропии в газоводе за га­

зогенератором с

 

распреде­

ми (т. е.

dbpldx = dbw/dx = 0), но

 

ленным выходом

учитывая,

что параметры

газа

в

нат: X

 

 

• тракте зависят от двух коорди­

и у. Несмотря на зависимость параметров от двух коорди­

нат, конвективные члены с производными от параметров по і

в

уравнениях неразрывности и энтропии

(4.42)

отсутствуют,

так

как в рассматриваемом случае

отсутствует

поперечная

состав­

ляющая скорости среды.

 

 

 

 

 

 

 

 

 

Подставив в первое уравнение системы

(4.42) искомое част­

ное периодическое решение, найдем

 

 

 

 

 

 

 

 

 

 

 

дойр-}-и;

= 0 .

 

 

 

 

(4.59)

 

 

 

Т г

èy

 

 

 

 

 

 

 

 

Вариация плотности связана с другими параметрами соотноше­ нием (4.44), зависимость вариации энтропии от поперечной ко­ ординаты X определяется зависимостью. (4.55). Для вариации энтропии, кроме того, необходимо учесть транспортное запазды­ вание при движении газа по газоводу. Тогда для точки газовода с координатами (х, у) изменение энтропии связано следующим

2 0 8