Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Гликман Б.Ф. Автоматическое регулирование жидкостных ракетных двигателей

.pdf
Скачиваний:
56
Добавлен:
25.10.2023
Размер:
13.98 Mб
Скачать

(1 — Р) (РыІРгг)

г

" г

(7.27)

р +

(1 —

р) (Ры!Ргѵ)(%г 1>/<г

 

,

,,

р

ѵ(І'р—1)/-

 

W ,

 

Р -г (1 —

р) (РыІРгг)

г

 

р

 

 

 

 

 

рс — давление за соплами турбины.

приведе­

Графики для коэффициентов расхода турбины ат и аТр

ны на рис. 4.1, для коэффициента скорости истечения газа из со­ пел турбины ßT — на рис. 5.1. При расчете динамики газогенера­ тора, подающего газ в чисто активную турбину, коэффициент рас­ хода ат находится по формуле (7.22), в которой в качестве про­ тиводавления используется давление в газовом тракте рм, а при расчетах-для турбины с заметной степенью реактивности р — по формуле (7.6), в которой в качестве противодавления использу­ ется величина давления за соплами турбины рс.

В коэффициенты ги г2 и гз, описывающие влияние степени реактивности р, входят параметры b и с из эмпирической зависи­ мости для степени реактивности турбины:

где и — окружная окоростъ лопаток турбины; с\ — скорость истечения таза из сопел турбины.

Подставив в уравнение газового тракта (7.24) вместо амплитуды вариации расхода на выходе 6GT из соотношения (7.23) амплиту­ ды вариации расхода из зависимости (7.25), находим уравнение газогенератора с учетом реактивности турбины на выходе:

(7.28)

278

Кроме уравнения 'баланса массы в газогенераторе (7.21) или (7.28), для расчета динамики двигателя необходимо второе урав­ нение— для температуры газа на входе в газогенератор (7.20).

7.3.2.Учет гидравлического сопротивления

вгазовом тракте газогенератора

Анализ результатов расчетов динамических характеристик газового тракта с учетом и без учета гидравлического сопротив­

ления столба газа в тракте показал

(§ 4.5), что при относи­

тельно небольших перепадах давления

вдоль тракта (порядка

процента от общего давления) гидравлическое сопротивление не оказывает заметного влияния на динамические характеристи­ ки тракта. При больших перепадах давления (приближающихся или выше 10% от среднего давления) влияние этого перепада на АФХ газового тракта уже достаточно существенно (см. рис. 4.11), причем особенно резко отличаются кривые АФХ для вариаций давления на входе и на выходе газового тракта.

Таким образом, расчет по некоторому среднему давлению в газогенераторе приводит к ошибкам в расчетах вариаций расхо­ дов жидких компонентов на входе и газа на выходе из газогене­ ратора, прямо связанным е колебаниями давления в газогенера­ торе. Кроме того, на динамические характеристики газогенера­ тора достаточно существенно влияют перепады давления на форсунках (см. рис. 4.7) и на турбине (см. рис. 4.9). При расче­ тах по среднему давлению в газогенераторе вносится ошибка как в величину амплитуды перепада давления на форсунках, так и в величину амплитуды перепада давления на турбине, что приводит к дополнительным погрешностям в расчетах АФХ газогенератора.

Перепад давления по газовому тракту из-за гидравлического сопротивления определим, используя обычное уравнение пидравлики:

А Ar г Ргг Рт ^гг

>

Ргг

где рт— давление у соплового аппарата турбины; р гг— давление у головки газогенератора; Сгг— коэффициент сопротивления; ргг— средняя плотность газа в газогенераторе.

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

р

Аіт + Рт

. ГГ

Агг(7'гг + 7'т) ’

(7.29)

(7.30)

где Тт — температура газа на выходе из газогенератора

перед

соплами турбины.

 

 

После линеаризации и приведения вариаций к

безразмерно­

му виду уравнения гидравличеакого сопротивления

(7.29)

с уче-

279

том зависимости (7.30) (Принимает вид

S^Prr==-^I- Sfirc--- ~ ~ Ърх= 2SGT—

 

 

Лргг

Д/?гг

 

- ^ h - t P r r -------5

8pT+ ^ - 8 f rr + -i-87\.

(7.31)

Ргг + Рт

Ргг + Рх

 

 

Амплитуда вариации температуры перед турбиной определяется

решением уравнения для транспорта энтропийных волн

(4.1.5):

8ГТ= е- '“иЬТ„ -(- —— 1(1 _ е- ,(« 0 8Лг-

(7-32)

»г

 

Учет гидравлических потерь приводит «к появлению в уравнении газогенератора дополнительной переменной — амплитуды вариа­ ции давления перед сопловым аппаратом турбины брт. Эта ва­ риация войдет также в уравнение для расхода газа через сопла турбины (без учета реактивности):

85т= - - L 8ГТ- И 1- схт) ѣрг+ ат8Я,.

(7.33)

Из уравнений (7.31) и (7.33) можно исключить амплитуду вариа­ ции брт и после преобразований найти соотношение для ампли­ туд вариации расхода газа через турбину при учете гидравличе­ ского сопротивления по газовому тракту газогенератора:

80т=

Dx

 

 

e - ,mx-}8prr+

 

 

 

 

+ — 8Ам----

- « ,[ ( ■ Pvr

Н-

(1 — ат) X

1 Dx

AD-

P -у

 

 

 

X

 

 

(7.34)

где DT= l - f ( l —ат

Ртт

 

 

 

 

Рт

 

 

При рт=ргг, т. е. при отсутствии ігидравличеокого сопротивления, несколько ироімоздкое соотношение (7.34) переходит в обычное выражение (7.33). В уравнении (7.34) пренебрегли из-за их отно­ сительной малости членами

-'-г — 1 Y _ P r r \ 2 _ l'

4хг ІА Рх )

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

280

 

 

 

■80т=

[! + «■

 

 

 

 

 

 

 

 

 

 

 

 

=

I Г3

РтГ 1

 

 

 

 

 

 

 

 

 

тр

 

 

 

 

 

-j-aTp

2 +

~y

ßTrjj 8pu--- —(1 +

a-rpr i) °TT-\-atpr.fiii,

и тогда 'получаем

другой

вариант

уравнения

для

амплитуды

вариации

 

расхода

через

турбину,

пренебрегая

членами

■Х-Р ---- 1

 

 

 

 

 

 

 

 

 

 

 

 

 

4%г

L ß f M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 a

 

Р ГГ

/ 1

'*т '

 

 

 

 

 

 

 

 

t

i l t

Ру

 

 

2%г— h arp) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

г

1 ( 1 4 -

a Tp r x)

 

Ьргг - ( -

a

Tp

( r 2 + ^

-

ß T r

r ) 8 / > M +

 

 

2xr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'Tp

)

-

l] +

 

 

 

 

 

 

 

 

 

 

^

 

+

{2+ (l +

a.rp) [ ( - ^ ) 2-

l] +

 

2aT„r1} e - '- . \ 8 7 ’rr} , (7.35)

где

Drp= 1 —|—(1 —)—oTp)

Prr

 

 

 

 

 

 

 

 

 

P t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

y.p —

1

 

I

 

 

 

 

 

 

 

CtTp— CCTp I г3

 

 

 

- y ß x ' - l - l ) -

 

Подставив 6 G t и з зависимости (7.35) в уравнение (7.24), найдем уравнение газогенератора с учетом гидравлического сопротивле­ ния газового тракта и реактивности турбины на выходе из газо­ генератора:

ІшTj

 

%г—1

l _ [ Pt l Y (

 

 

°У ? \

Рт

 

1

/11

\

1

 

D.

(1 + атргі) —------

 

тр

 

 

 

 

 

тр

 

 

-

I

1

Д гг

 

4D.тр

Р у

 

 

 

 

- { ! — ^W. ^ ( ^ “W + C l+ ^ p )

! _ JhL^L + aTp)_|_

2-х.р

е - '" т‘

8 л г +

тр

 

х.г--1 ,

г

 

1 ]]}е- '^ } } х

Х Щ Г

Gpl

e -imx^Q

^оі_ е_/“тг8(7о І= 0.

(7.36)

 

П

281

Вторым уравнением и в этом случае служит уравнение для ва­ риации температуры (7.20) *).

7.3.3. Двухзонный газогенератор

Для двухзонного газогенератора уравнения динамики были выведены в гл. IV. Однако форма этих уравнений удобна для анализа особенностей динамических характеристик газогенера­ тора и неудобна для расчета динамических характеристик двига­ теля в целом, так как приводит к появлению в одном уравнении членов с разными запаздывающими аргументами для одной пе­ ременной (вариации параметра). С этих позиций лучше использо­ вать в качестве внешних возмущающих воздействий не колеба­ ния давления перед форсунками, а колебания расходов компо­ нентов.

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

стирующее устройство

(решетку и т. д.)

сср, вариация

давления

в емкости

Ьргг — на Ьр'п , а

вариация

 

давления

за

емкостью

бри — вариацией давления во второй зоне öpr/ J-

 

 

Учтя все это, запишем уравнение для

колебания давления в

первой зоне двухзонного газогенератора в следующем виде:

/СОТ,

 

1

 

* г-1

е

 

-1 3^rr +

apS/7rr —

 

 

+

1 - а п

Ъ.

 

 

2-/.„

 

 

 

 

 

 

 

-

1

— і и

 

 

 

ol '

 

 

 

1

 

о Т г г —

 

 

 

 

 

 

 

C’

o;r

 

 

 

 

 

 

 

 

 

 

 

 

(7.37)

 

 

 

п e— ,'BGri—0,

 

 

 

 

 

G'

 

 

 

 

 

 

где

Ti — время пребывания газа в первой зоне газоге­

 

 

нератора;

 

 

 

 

 

 

 

■/-г ■—показатель адиабаты газа в первой зоне;

 

Тг — время преобразования

жидких

компонентов

 

 

в газ в первой зоне;

 

 

 

 

 

Grr= Gor+ Gri

— расход газа через первую зону;

 

 

(

2

 

- к < + і ) , ч ( Р ы ' р ф ^

 

 

« р =

 

 

 

 

 

 

 

 

 

—коэффициент расхода через балластировочную решетку.

*> Изменением энтропии из-за трения газа о стенки пренебрегли.

282

Уравнение для амплитуды вариации температуры у головни в первой воне практически совпадает с уравнением (7.20):

8Ггг—фгге ш'г ( 80оі — 8Gri) = 0,

,

k'TT

дт'гг

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

где фгг^—;-------

;------

 

Т'„

дк„

 

первой зоны.

Уравнение второй зоны (4.40) выведено о многими упрощаю­ щими предположениями и в форме, неудобной для использова­ ния при расчетах динамики двигателя. Проделав такие же пре­ образования, как и при выводе уравнения (7.36), т. е. учтя гид­ равлическое сопротивление столба газа во второй зоне и реактив­ ность турбины, находим уравнения для амплитуды вариации давления во второй зоне:

/а)Tj

 

 

С?„а

 

 

 

 

 

 

!

СХтр

1

 

 

 

 

 

 

 

 

 

 

р р

I

D.тр

 

 

 

 

I

+

 

 

 

 

 

 

 

X

е

Ьргг I

 

~

Г ~

( 1 + ° т РГ j )

-

 

Xг

 

 

D

Tf,

 

 

2 х т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«тр'-І Ъп,

 

 

 

 

 

 

 

 

 

 

 

Dтр

 

 

 

Г-4

 

Ѵ -„ -

1

*тр

 

 

 

 

 

 

2ѵ.

 

 

Рт

 

 

 

 

4D.тр

 

 

 

 

1 \

 

 

 

1

2(1-}- аТргI)—(1 -)- Итр) (— Y

 

8г;г-

 

 

W тр

 

 

 

 

 

\

дт

/

 

 

 

 

и г1

—' “■'ггЯ”

 

 

е

/шТг80; + 0 р(ар-1)8/?;г II

 

— тс— е

г8(?г

 

 

 

 

0,-г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

1

-

^

= 0.

 

 

(7.38)

Здесь

Ті — время пребывания газа во второй зоне газогене­

 

 

 

ратора;

 

 

 

 

 

 

 

 

 

Ргг,

хг — показатель

адиабаты газа во второй зоне;

 

8/?гг— давление

и амплитуда

колебаний

давления за

 

 

 

балластирующим вводом, в начале второй зоны;

 

тг — время преобразования

жидкого

компонента в

 

 

 

газ во второй зоне;

 

 

 

 

 

283

On + Goi

Grr— Gn + Ooi + Gri + Оой

 

aTp = aтр

— — г,

xr - 1

 

----Ртг1— 1

 

 

2

1

/ Л р = 1 +

1

 

 

Ori, Ooi, On, Ooi— расходы

горючего и окислителя в первую

 

зону и во вторую зону;

 

ЬТр — амплитуда вариации температуры газа пе­ ред балластирующим вводом, в конце пер­ вой зоны.

Для общности в уравнение (7.37) введены вариации двух жидких компонентов, подаваемых во вторую зону. Для окисли­ тельного газогенератора надо принять амплитуду вариации рас­

хода oGri = 0 , а для восстановительного— 8С?о1= 0 . Если гид­ равлическое сопротивление газового тракта во второй зоне газо­ генератора мало, им можно пренебречь, приняв в уравнении (7.38) ргг/рт=1. В случае малой реактивности турбины р<СІ, в уравнении (7.38) можно принять

атр= —ат; г3 = г1= 0; г2= 1 .

Амплитуда вариации температуры газа в конце первой зоны, перед балластирующим устройством, определяется уравнением, аналогичным уравнению (7.32):

ЬТр— е

Ѵ гг ---- т—

t l — е

1j8^rp=0.

 

Д

 

 

 

Уравнение для вариации температуры газа за

балластирующим

устройствам, в начале второй

зоны, находим

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

(4.41), заменив вариации давления перед форсунками вариация­

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

изменив

форму записи коэффициентов:

 

+ -° ^ гг ■фггв-10^ (8Öri-

8Ön) + фгГе_ 'ШТг (Ö'oibÖh -ÖrtSÖrt) +

1 + krr

 

 

-(- фгг (Ool örl) ( 1.

Ctp) ЗрггД'фгГ (öol --On) OpS/?rr---

— J {O'oi — örO <МГр —87\r= 0 .

(7.38)

284

З д е с ь

4

/~і

^oi

.

/~л"

Оог

=

1

/~*г

 

О 0\ =

“ ---------

>

С/оі = -----;---------

 

1 — С/оЬ

 

 

Gol +

Gol

G o\

+

Gol

 

 

 

 

5 ri=

, ° rl .

; 0 A= -

° rl .

=

1

- ö ;i;

 

 

Grl +

Grl

Grl +

Grl

 

 

 

krr= 0 0ijOri — соотношение компонентов в первой зоне;

Ts= Ti-f-Tr — суммарное время задержки для первой зоны. Для окислительного газогенератора

80г1

= 0 ;

Ооі =

1,

а для восстановительного

 

 

 

8Ооі =

0;

С?гі =

1.

Вариация температуры газа перед балластирующим устройством определяется условием переноса энтропии (4.15):

8Тр= е ~ іаТі87;;г+ - --^ - ( 1 - е- ‘“Ті).

(7.39)

7.3.4. Газогенератор с распределенным выходом

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

Динамические характеристики газогенератора о распределен­ ным выходом описываются уравнением (4.57), выведенным без учета реактивности турбины. Кроме того, это уравнение второго порядка, что неудобно для использования’ его в общей системе уравнений двигателя.

Для учета реактивности турбины необходимо соотношение (4.63), определяющее вариацию расхода газа через сопла, при­ вести к виду, аналогичному виду соотношения (7.25):

— -------- (1 + а трГ l) 8scp 4" 1 + “ ip -------------ßxr i —

X, — 1

 

^ P r г "Ь а тр 8 л -

 

+ °тр

Ртг іj Ърм,

(7.40)

285

где öscp— амплитуда'Средней мгновенной вариации энтропии газа за соплами турбины.

Отличие зависимости (7.40) от (4.63) состоит в разной записи ко­ эффициентов при вариациях и появлении дополнительного члена с вариацией частоты вращения. Эти отличия не вызывают какихлибо затруднений при преобразовании уравнения газогенератора.

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

опособу осреднения аналогична (вариации средней по

колесу

адиабатической работы газа (5.18):

 

 

(4 ткол + 1) Sscp = 8зкол = е - '“н ( ЬТгг-

ъ Щ ,

(7.41)

где

SsK0JI— амплитуда вариации энтропии на входе в кол­

 

лектор соплового аппарата турбины;

 

тк01—

Ргг^ к'ол----- условное время пребывания

газа в

коллек-

 

Gf

 

 

торе.

Ргг^кол — плотность и объем газа в коллекторе. Использовав зависимости (7.40) и (7.41), получаем уравнение однозонного газогенератора с учетам распределенности выхода из негон реактивности турбины:

 

ш Ті + Ткол

-f- 1

иг— 1

(1 —2е—''“"■)Д-

 

 

 

 

2 хг

 

 

 

1

Up 1

1

 

Ь р гг -J- а т, г 2 4 - — Г,3Т Ь р ы —

— Г, — !-

 

 

^Ткаі+ - у ( И - Г1атр) 8scp- (

l - e-'-.)8 7 ’rr+

 

-)-а-грГ^я — е"

 

- ^ 8 0 о1 +

- ^ 8 0 г1 )= 0.

(7.42)

 

 

 

Gpp

 

Grr

 

Уравнения (7.20), (7.4Г) и і(7.42) .определяют дпнамичеокие ха­ рактеристики однозонного газогенератора с распределенным вы­ ходом.

7.3.5. Газовод

Газ, поступающий во входное сечение газовода (после турби­ ны), имеет температуру, отличную от температуры перед турби­ ной из-за работы, совершенной газом на турбине. Для опреде­ ления падения температуры газа в турбине воспользуемся тер­ модинамическим соотношением

С2

 

1

(7.43)

2 4 ,

286

где Діт— теплоперепад на турбине (изменение энтальпии газа);

Рт— к.'п.д. турбины;

 

 

Т ы— температура газа за турбиной.

(7.43) на­

Приняв ср= const, после

линеаризации зависимости

ходим

 

 

 

- ^

Ь Т , -

т ^ ^ Ь Т м= 2 ^ 1+ ^

(7.44)

т

*м

*т 1м

 

Воспользовавшись выведенными ранее уравнениями для вариа­ ции скорости на выходе из сопел турбины бщ и вариации клі.д. турбины ötjt (5.7) и (5.13), приведем уравнение (7.44) для ам­ плитуды вариации температуры газа за турбиной öTu к следую­ щему виду:

8 7 \

3

7Т 7М

 

 

 

Т,

т

ь -

2ТЫ

 

 

 

 

 

X

е~ - э ) ]

Тг -

Ты р / 1

 

 

 

 

 

 

А г - ^ < М *

+ ( 1

Тт Тм -Фс. )е~іап'ЪТгг.

(7.45)

 

 

 

 

 

2ТЫ

 

В уравнении

(7.45)

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

газо­

генератора без учета гидравлических потерь по газовому тракту и распределенности выхода из газогенератора. В случае двухзон­

ного газогенератора уравнение (7.45) сохраняет свой

вид, но

вместо вариации брГг необходимо подставить вариацию

давле­

ния газа во второй зоне

8ргг, а вместо вариации бГгг— соответ­

ственно вариацию 8Ггг.

При расчетах газового тракта с учетом

распределенного выхода из газогенератора необходимо изменить только один член уравнения (7.45), содержащий амплитуды ва­ риации температуры 67Ѵ Из сравнения выражений для вариаций

скорости бс, (5.7)

и (5.18)

видно, что вместо члена

8Тт= 8 f rr е - ' “ т > + - Х г ~ 1

(1 - e - , m «) f p rr

 

 

Хр

 

появляется член

—•

_2 —.

Произведя соответству-

Sscp-|---- ------- Ъргт.

Х г

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

Уточнение уравнения (7.45) путем учета гидравлического со­ противления по газовому тракту приводит к существенному ус­ ложнению этого уравнения, вряд ли’оправданному, так как раз­ ница между температурами Тши Тгг и их вариациями невелика.

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

287

Соседние файлы в папке книги из ГПНТБ