книги из ГПНТБ / Гликман Б.Ф. Автоматическое регулирование жидкостных ракетных двигателей
.pdfbk' — вариация соотношения компонентов в этой же порции газа в «момент ее образования у -головки 'первой емко
сти;
I
Ті — время пребывания газа в первой зоне; t
тг — время преобразования жидких компонентов в газ в первой зоне.
Из определения соотношения компонентов в газе в начале вто рой зоны
|
|
|
|
o;.p + g; |
, |
|
|
|
|
находим линеаризованную зависимость |
|
|
|
|
|
||||
|
8£"=(1 —Оо)800.р+Оо80о — (1 —Gr)8Gr.p— Gr8Gr, |
(4.33) |
|||||||
где |
|
800.Р, 80г.р— относительные |
вариации |
расхода |
|||||
|
|
|
|
окислителя и горючего, приходяще |
|||||
|
|
|
|
го с порцией |
газа к |
вводу |
балла |
||
|
|
|
|
стирующего |
компонента |
(на |
выхо |
||
|
|
|
|
де из зоны /); |
|
|
|
|
|
|
|
|
8G0, |
8Gr— относительные |
вариации |
расхода |
|||
|
|
|
|
окислителя |
и горючего, |
поступаю |
|||
|
|
|
|
щего в газ |
через балластирующий |
||||
|
|
|
|
ввод (в конкретном газогенераторе |
|||||
|
|
|
|
имеется только одна из этих вариа |
|||||
г? |
Gl |
^ |
|
ций); |
|
|
|
|
|
G” |
|
доли расходов окис |
|||||||
Go= —;----- ;; |
и г= |
-------- -— относительные |
|||||||
|
ло + G 0 |
|
G0 + |
G0 |
|
|
|
|
|
|
|
|
|
лителя или горючего, поступающего |
|||||
|
|
|
|
во вторую зону (для конкретного |
|||||
|
|
|
|
газогенератора |
один |
из этих |
коэф |
||
|
|
|
|
фициентов равен нулю). |
|
|
Расход газа через балластирующую решетку между первой и вто рой зонами Gp связан балансовыми соотношениями с расходами соответствующих компонентов, входящих в состав газа:
GО . р - |
Gr.p |
*р + |
1 |
или в линеаризованном виде —
8Go.p=80p -j---- 5^" ;
і + *Р
Gp
*р + 1
(4.34)
BGr.p=8Gp |
k’ |
8Ap, |
(4.35) |
£----- |
l+Ap
188
где öGp — относительная вариация расхода газа из первой зо ны во вторую. Для вариации расхода через местное гидравличе ское сопротивление между первой и второй зоной можно вос пользоваться зависимостью типа соотношения (4.10) и, кроме того, уравнением (4.15):
80р = — і- ЬТ'Р+ ар (Ьр" - 8р') + Ьр' =
== — 1— Ф' [SGo {t —xI — тг)] — Юг (t — х[ — Tr)] —
[ 8 р ' - 8 / , ' ( * - ^ ) ] + а р (8/ / ' - 8 / ,') + 8 р', (4.36)
где 8р', Ьр" — относительные вариации давления в первой и второй зоне газогенератора;
ар — коэффициент расхода для гидравлического со
противления между первой и второй зоной, вы числяемый по формуле (4.11);
,,k ’ д Т ' .
ф----- безразмерный наклон зависимости температуры
газа в первой зоне от соотношения компонентов;
87^ — относительная вариация |
температуры газа в |
конце первой зоны. |
’ • |
Вариации расхода жидких компонентов через форсунки в первую
и вторую зоны газогенератора связаны |
с вариациями |
давления |
||||
линеаризованными |
зависимостями |
типа |
соотношения |
(4.26): |
||
Ь О г = ± - ^ Ь Рт |
2 рг—р7 8 /; |
|
||||
|
2 рг—р |
|
||||
ЬОо- |
1 |
Р о |
1 |
|
|
|
2 Р о — Р |
2 |
P o — Р |
|
|||
|
(4.37) |
|||||
ЬѲг |
1 |
Р г |
1 |
Р |
|
|
ѵ ьр ; |
|
|||||
7 ---- °Pr |
|
|
|
|||
|
2 рг—р' |
2 рг—р |
|
|||
ь о : = ^ г ~ ^ , ь р о |
1 |
Р |
8/Л |
|
||
|
2 |
Р о — Р " |
2 Р о — Р " |
|
где рг, Ро, 6рг, бр0 — давление и относительные вариации давле ния в точках развилки трактов компонен тов, подающих их в зоны I и II.
Уравнение первой емкости двухзонного газогенератора получа ется подстановкой двух первых соотношений (4.37) в уравнения газового тракта типа (4.24) и (4.15):
/coTj X—2%1( 1 - е
199
|
|
1 |
|
|
|
s y + OpSjO" — |
|
2 \ k ' + |
I р 0 — р ' |
f t ' + |
1 Р г — p ' j |
|
|
||
|
|
|
|||||
1 ___У о - '“-1! \ ST'____ *_______ £о_ |
|
e 'ШХг8л,+ |
|
||||
---- L e - ' ^ |
s r « |
— |
2 (Po— p ') |
|
|||
2 |
/ |
|
(ft' н+- J) |
|
|
||
|
1 |
p T |
— / cut, |
|
|
(4.38) |
|
|
f t ' + |
1 2 ( р т — р ' ) |
г8а ; |
|
|||
|
|
|
|
|
|||
s |
r I+ Ciü |
^ |
Р ’ |
е ,a"«V = |
|
||
|
2 \Ро— р ’ Р г— Р' |
|
|
|
|||
£ |
- * 4 |
«"'“ '»Л - f |
Рг—Р |
е-'-МЯ |
(4.39) |
||
2 До—Д |
|
2 |
|
|
|
По аналогии с уравнением газового тракта за турбиной (4.25),
использовав соотношения (4.32), (4.33), (4.34), (4.36), (4.37),
получаем уравнения для Еторой емкости двухзонного окислитель ного газогенератора:
|
20)Т- |
|
1— ат - |
ft' + |
1 |
, |
X |
— 1 |
,, |
— |
Ыт , . |
, |
||
|
— |
|
|
|
ßn-l--------( 1 - е |
|
») + |
|||||||
|
X |
|
|
|
f t " |
+ |
1 |
Р |
2 х |
|
|
|
|
|
|
f t " — |
f t ' |
р " |
|
|
|
|
|
8//' + |
^ - ^ ( а р- 1 ) + |
||||
|
f t ' + l |
2 ( р 0 — р " ) |
|
|
|
У I |
k„+ |
j V. Р |
' I |
|||||
|
|
|
|
|
|
|
|
|
|
|||||
|
I X — |
1■(1 — е |
|
!) |
So' |
2 |
ft" + |
е ;“ті8Г — |
||||||
|
2 |
х |
|
|
|
|
|
1 |
1 |
|
|
|
||
- ( 1 |
---- -- е |
/ |
Ь Т " - \ - а |
Ь р |
= |
|
—— - |
|
/>0 |
|
е“ ,в*'18/70; (4.40) |
|||
\ |
2 |
' |
|
т " ы |
|
ft" + |
1 |
2(р 0 — р") |
|
|
||||
8^ - f |
|
|
|
|
1 |
|
|
|
1/1 —/шт.. |
|
||||
|
|
|
+ ■ 1 - S - — d - e |
>) V - |
||||||||||
|
—/сот |
|
|
' |
P' |
|
|
P' |
|
|
|
Ь р ' - |
||
|
e |
ЧТ' |
|
,Рг—Р' Po— p ' |
|
|
|
|||||||
|
2 |
|
|
|
|
|
|
|
||||||
|
|
|
1 + f |
t ' |
x — |
1 |
(1 |
|
|
|
|
|
|
|
|
|
|
1 + |
ft" |
X |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
V |
Pr> |
U J k " + |
1 с - / < о ( ч ( + т ' ) |
2 k " p 0 — p " _ |
f t ' + |
1 |
(Ä^.— Ä') e—/wTr |
8Am (4.41) |
где ft"— соотношение компонентов .в зоне I I газогенератора; ат — коэффициент расхода газа через турбину;
х "— время пребывания газа в зоне //;
ЬТ "— амплитуда относительной вариации температуры за
2üÜ
балластирующим вводом (в момент испарения жид ких компонентов, поступивших в зону II) ;
т г— время преобразования |
жидких компонентов в газ в |
доне II) . |
|
На рис. 4.13 приведены кривые АФХ Д'вухзонного газогенера |
|
тора, рассчитанные по формулам |
(4.38) — (4.41) и для сравне- |
Рис. 4.13. Сравнение АФХ двухзонного газогенератора с АФХ
однозонного газогенератора
ния;—кривые АФХ однозонного газогенератора, имеющего оди
наковые с двухзонным время пребывания |
( и другие |
характеристики: ф=<|/'; тг = Т г= тг = 0. |
Принято, что т"=т", |
а Ѳ= шт—сіДт'+т")- Для простоты перепады давления на форсун ках приняты одинаковыми, т. е. ро—р' = рг—р'=Ро—р". Возму щение гвносится по давлению окислителя перед форсунками.
Сопоставление кривых показывает, что наличие зоны II в га зогенераторе достаточно существенно влияет на его динамиче ские характеристики. При одинаковом суммарном времени пре бывания у однозонното газогенератора фазовый сдвиг больше, чем у двухзонного. Характер кривых при частоте Ѳ>3-М резко отличается. Таким образом, наличие второй зоны в газогенерато ре необходимо учитывать при расчетах динамических характе ристик двигателя.
201
4.7.ГАЗОГЕНЕРАТОР
СРАСПРЕДЕЛЕННЫМ ВЫХОДОМ
(КОЛЛЕКТОРОМ ТУРБИНЫ)
Газ из газогенератора поступает в 'турбину ТНА обычно через коллектор 1 перед сопловым аппаратом 2 турбины 3 (рис. 4.14). Если газ поступает в коллектор турбины в одном сечении, а за тем по коллектору подается к отдельным соплам турбины, то та
кой коллектор представляет из |
себя распределенный (по длине |
|||||
газового тракта) выход из |
газогенератора. Ранее, |
при |
выводе |
|||
; |
уравнений |
газового |
тракта |
в § 4.2, |
||
было сделано предположение, что |
||||||
|
молено пренебречь временем |
пробе |
||||
|
га |
акустических н |
|
энтропийных |
||
|
волн по длине сопла (соплового ап |
|||||
|
парата). В рассматриваемом случае |
|||||
|
такого предположения |
уже |
сделать |
|||
|
нельзя, так как необходимо |
учиты |
||||
|
вать |
время |
пробега |
|
энтропийных |
|
|
волн по длине коллектора. |
|
||||
|
Вернемся к выводу |
упрощенных |
||||
Рис. 4.14. Схема газогенерато |
уравнений |
газового |
тракта |
(4.21), |
||
ра с распределенным выходом |
(4.24) и др. Эти уравнения были по |
|||||
|
лучены путем разложения в ряд не |
|||||
которых членов более точного, но громоздкого решения |
(4.13) — |
(4.14). Преимущество такого метода состоит в том, что, с одной стороны, получается и точное, и приближенное решение, а с другой, — легко оценить отбрасываемые члены.
Однако эти же упрощенные уравнения можно получить непо средственно из исходной системы уравнений гидромеханики. Для этого исходную систему уравнений гидромеханики удобно запи
сать в первоначальной ферме, не |
исключая |
плотности |
(сразу в |
|||
линеаризованном виде): |
|
|
|
|
|
|
asp , |
|
( р » )= 0 . |
|
|
|
|
dt |
дх |
’ |
|
|
|
|
dbw |
, |
о dbw , |
аЗ |
дЬр |
„ |
(4.42) |
W ------ |
\- W--------------- - = |
0; |
||||
dt |
|
дх |
X |
дх |
|
|
dbs I |
дЬs |
п |
|
|
|
|
-----\-W---= 0 . |
|
|
|
|
||
dt |
дх |
|
|
|
|
|
Здесь все вариации относительные. Для упрощения введем предположение, что при относительно .низких частотах, когда со//а<СІ, акустическими эффектами можно пренебречь, т. е. что dbpldx = d6wldx—0. В этом случае 2-е уравнение системы (4.42) исчезает.
202
Проинтегрируем по х ('координата вдоль тракта) первое урав нение системы (4.42):
I |
|
8ро'л:4-™[8(рта)вЫХ-В(рда)вХ]= 0 , |
(4.43) |
о |
|
где 6 ( р г « ) в х , б ( р г и ) в ы х — относительные вариации расхода газа в начале и конце тракта. Решение третьего уравнения системы (4.42) остается тем же, что и решение для 3-го уравнения систе мы (4.3), т. е. соотношение'(4.15), Учтя термодинамическую за висимость
|
s = c 0ln р — ср 1п р -j- const, |
|
|
||
которая после линеаризации принимает вид |
|
|
|||
|
|
8 р = — Ьр— 8 5 , |
|
(4.44) |
|
|
|
|
% |
|
|
преобразуем уравнение (4.43): |
|
|
|||
|
— |
[ |
S5 — = 80вх- 8 0 т, |
|
(4.45) |
|
dt J |
|
|||
|
к. |
w |
|
|
|
|
|
о |
|
|
|
где 8/»—-производная по времени. |
|
|
|||
При |
преобразованиях учли принятое условие | ^ = 0, |
а |
вариа |
||
цию |
8(ри>)вЫХ заменили |
вариацией расхода через турбину |
80т, а |
||
M p®Obx=&öbX. Проведем в уравнении (4.45) замену |
dxjw = dx, |
где %— текущее значение времени пребывания газа в тракте до
сечения с координатой х, и, используя решения 3-го |
уравнения |
системы (4.42) в форме |
|
8s= 8sBX(^ — т), |
(4.46) |
после преобразований находим |
|
I l S/)+ 8sbX( ^ - t1) - 8 SbX=8G bX- 80t. |
(4.47) |
X |
|
Здесь 6sbx — относительная вариация энтропии на входе в газо вый тракт. Если заменить 6sBXс помощью соотношения (4.8), то получим другую форму уравнения тазового тракта:
Z |
8^ + 87'«^—t 1)- 8 7 'BX+ ^ i [ 8 ^ - 8 /, ( ^ - t 1)]= |
80ie-8 G T. (4.48) |
|
|
X |
|
|
Используя соотношение |
(4.10), это уравнение |
можно свести к |
|
уравнению (4.22). |
|
|
|
|
Вернемся к задаче о газогенераторе с распределенным выхо |
дом.' Для потока газа в тракте с переменным расходом необходи-
203
мо изменить исходную систему уравнений гидромеханики. Рас смотрим вначале только входной 'коллектор турбины, пред полагая, что сечение его постоянно, а расход газа" по длине уменьшается вследствие оттока через сопла, равномерно рас пределенные по длине коллектора. Для течения газа с перемен ным расходом изменяется первое уравнение системы (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