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

книги из ГПНТБ / Волков Е.Б. Ракетные двигатели на комбинированном топливе

.pdf
Скачиваний:
27
Добавлен:
24.10.2023
Размер:
8.82 Mб
Скачать

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

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

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

В качестве основной характеристики заряжания можно использовать применяемый при проектировании РДТТ коэффи­ циент заполнения сечения камеры зарядом (см. работу [42]):

є = S J F j

(4.25)

где ST — текущее значение площади поперечного сечения заряда; Fк.с — площадь поперечного сечения камеры.

Масса заряда твердого компонента определяется как

 

 

(4.26)

Подбор начального значения параметра 8o=STo/FK.c,

диамет­

ра камеры DK.C, числа камер п и длины заряда L

при

принятой

его форме и составляют содержание поставленной

задачи.

Рассмотрим связь между параметром заряжания ео и пол­ ным временем работы двигателя т при заданном диаметре ка­ меры сгорания DK.C. Для установления этой связи используем элементарную зависимость: ^

-FKCdz

= U(i)udt,

(4.27)

где П(е) — текущее значение периметра горения, рассматривае­ мое для заряда принятой формы как однозначная функция теку­ щего значения г; и — текущее значение скорости газификации.

Отсюда полное время работы определится так:

о

(4. 28)

При решении данной задачи используется осредненная по длине заряда скорость газификации, определяемая зависи­ мостью (4.9), которую можно представить в виде

и = иг

а ж.г

АХ(М).

(4. 29)

 

 

 

 

по

Подставляя

уравнение

(4.29)

в (4.28),

получим

 

 

 

о

 

 

 

 

r =

_ ^ L

Г

( 1 - О р * »

,

.

( 4 . 3 0 )

Входящие в подынтегральное

выражение

 

величины

Gm.T и р н

следует в общем случае рассматривать как переменные во вре­ мени, связанные с изменением е. Зависимость П(е) определяется формой заряда.

В простейшем случае для заряда с цилиндрическим каналом, горящим изнутри, П = яс?. Поскольку для такого заряда

 

 

 

 

Dl.

а"*

 

 

 

 

 

 

 

к -с

 

 

-

(4.31)

 

 

 

 

 

D2

 

 

 

 

 

 

 

 

 

к.с

 

 

 

 

Т О

 

tf = D K C f l - e ;

 

 

 

 

 

П ( е ) = я / > к . е

і Л - є ,

 

(4.32)

где d — текущий диаметр

канала.

 

 

 

 

Следовательно, для заряда с цилиндрическим каналом

 

 

 

 

 

о

 

 

 

 

 

т =

-

 

к -с

 

 

— .

(4.33)

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

ствующие типичным режимам работы ГРД.

 

 

Случай I . Работа

двигателя

при постоянном

расходе

жидкого

компонента

(Gm = const),

подаваемого

только

через

форсуноч­

ную головку,

при v = 0 (чистое диффузионное

горение).

 

Если принять, что для процесса

K(fi,

гр)—const, из уравнения

(4.33) получим

 

 

 

 

 

 

 

 

 

 

 

 

Fl+*

 

i _ ( i _ e o ) P + o . 5

 

 

 

 

 

к.с

 

 

 

 

(4. 34)

 

Р +

0,5

 

nD^Gl

 

 

*(р,<Ю

 

 

 

 

 

 

 

Поскольку н х =

F9

 

=

ик

— скорость газификациитвер-

 

 

к.с

 

 

 

 

 

 

 

 

'

 

 

 

 

 

 

дого компонента, соответствующая концу работы двигателя, уравнению (4.34) можно придать следующий вид:

х=

\

^ Е £ . [ і _ ( і _ е 0 ) Р + о , 5 ] .

(4.34а)

 

Р + 0,5

4 « к

 

При р = 0 , 5 т =

D k - c S ° .

 

(4.346)

 

к

 

 

Согласно уравнению (4.29) скорости газификации твердого компонента в начале и в конце работы двигателя связаны соот­ ношением

 

 

 

и °

-

1

 

(4.35)

 

 

 

«к

 

( 1 - £ о )

 

 

 

а для промежуточного значения є

 

 

 

 

 

 

 

 

 

 

(4.

36)

Относительное

изменение

расхода

твердого компонента

в процессе работы двигателя

определится

так:

 

 

G T

_

ad

 

I 1 — е 0

,5

 

 

G t

0

Utflfy

 

 

 

 

 

Максимальное

отклонение

расхода

твердого компонента

от

его начального значения, имеющее место в конце работы двига­ теля, составит:

а.

о т 0

Для топливных пар, у которых показатель степени в законе горения р>0,5, расход твердого компонента в процессе горения уменьшается (рис. 4.7), следствием чего является изменение температуры горения и давления в двигателе, что в свою оче­ редь приводит к изменению удельного импульса.

Максимальное относительное изменение параметра г|), харак­ теризующего соотношение расходов жидкого и твердого компо­ нентов, составит:

Фо_

1 — фр ч- ф 0 ( 1 — Ё 0 ) ^ - ° - 5

, 4 3 9 )

W

О - * , ) * - 0 - 5

 

Из зависимости (4.13) для определения

рабочего давления

при v = 0, а = 0, с учетом

выражений (4.32) и

(4.39), получим

Рй А ЖР.фо)'

Таким образом, падение давления при G)K = const в процессе работы двигателя определяется как снижением суммарного рас­ хода топлива, так и понижением температуры горения (коэф­ фициента истечения А ) вследствие изменения \р.

ГРД, работающий на режиме Gm = const, может быть исполь­ зован в объектах, для которых требуемый режим полета обес­ печивается при силе тяги двигателя, меняющейся в ограничен­ ных пределах.

При

решении

вопроса о целесообразности использования

режима

бж = const в проектируемом двигателе приходится,

с одной

стороны,

руководствоваться соображениями простоты

конструкции двигателя и функционирования системы топливо-

подачи на таком режиме, а, с другой стороны,

учитывать

про­

игрыш от

снижения

/ ] , вызванного

отклонением

ip от

оптималь­

ного значения и падением давления в двигателе.

 

 

 

Преобразуем зависимость

(4.31) в вид, удобный

для расче­

тов

при баллистическом

проектировании

изделий

с ГРД.

 

 

 

GT0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

• ^ j ^ j ^ ^ l / ^ i

 

У\

 

A

 

 

 

 

0,3 ^№^^

 

-

" "

i

 

 

\/

1

/

\

-

і

 

 

 

 

 

 

 

^ > ^ - l

^ 1

1 / ^

1

 

і

 

 

 

 

OS

 

 

1

>

^

і

 

 

і

 

і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

і

>

^

і

і

 

і

 

і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

!

 

!

 

!

!

 

;

 

'

 

 

 

 

0,6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.5 I

 

 

 

1

I

 

1

1

 

1

 

і

 

 

 

 

 

 

 

1I

'>

 

I1 .

I1

 

I1

 

1 1 —

 

 

 

 

0

0,2

 

 

OA

 

0,6

 

0,8 £

 

 

 

 

 

 

 

 

 

Рис.

4.

7.

 

 

 

 

 

 

 

 

Выразим расход жидкого

компонента так:

 

 

 

 

 

 

 

q

_

ы

ж

_

(1 — Ф)Н<?0

 

 

 

 

(4 . 40 )

где

Qo — стартовая

масса

изделия;

 

 

 

 

 

 

 

 

| i K — относительный

запас топлива;

 

 

 

 

 

 

 

 

— среднее значение -ф за

время

работы

двигателя,

опре­

 

 

деляющее

соотношение

компонентов

в

бортовом за­

 

 

пасе топлива;

 

 

 

 

 

 

 

 

 

 

 

 

 

п — число камер

сгорания

в пакете.

 

 

 

 

 

 

Как отмечалось ранее, линейная скорость газификации боль­

шинства

твердых

компонентов,

используемых

в

настоящее

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

При такой компоновке двигательной установки необходимо, чтобы отношение наружного диаметра отдельной камеры сгора­ ния к калибру изделия равнялось величине D, соответствующей плотной укладке камер в пакет при заданном п. Величина D определяется из элементарных геометрических соотношений

и зависит только от числа камер в пакете п. Значения D для различных п приведены в табл. 4. 1.

 

 

 

 

 

 

Таблица 4.1

п

3

4

5

7

12

14

19

ъ

0,464

0,414

0,370

0,333

0,244

0,220

0,200

пЪ

0,645

0,686

0,685

0,776

0,715

0,675

0,760

Площадь сечения одной камеры пакета связана с площадью миделя .FM зависимостью

При дальнейших преобразованиях используется обозначение

<2о/^м = По — начальная поперечная

нагрузка.

 

Подставляя

полученные соотношения в

равенство (4.34)

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

камеры DKX, получим

А , с = 4 ( Р +

0 , 5 ) т і - Ч

 

(4.41)

 

пЪв1

1 - ( 1

- £ „ )(3+0,5

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

Введем ограничительное условие:

Кс Я.ДОП-

Расход газов в конечном сечении заряда в начальный момент времени определится так:

 

Ос =Л/г к.с(1-е0 )<?(Яс )рсо,

(4.42)

где

рсо — полное давление

в газовом

потоке в этом

сечении;

 

q(kc) —газодинамическая

функция.

 

 

 

Поскольку « G c 0 = G s o , мы можем G c 0

выразить так:

 

 

О - _ 3 L _

 

(4.43)

где

Ро — начальная сила тяги двигателя;

 

 

Ло — удельный импульс при соотношении расходов компо­ нентов -фо и рабочем давлении в двигателе рк о-

Записав из формул (4.42) и (4.43) выражение для площади сечения камеры и подставив это выражение в формулу для массы заряда (4. 26), получим

1 wAq с ) рк0

1 — г0

При проектировании двигателей и летательных аппаратов (ЛА) различного назначения в качестве отправных параметров обычно задают начальную тяговооруженность r\o = P0/Qo и отно­ сительный запас топлива |лк=(й/<2о. ДЛЯ того чтобы ввести их в уравнение (4.44), разделим обе его части на стартовую массу ЛА <2оПри этом получим

 

т К = - 7 - т ^

 

 

(4-45)

 

 

А9 с ) рс0

1 — є 0

 

Введем

следующие

обозначения:

 

 

 

v _

Лон*

.

_

Ад(ксс0

 

 

If—

 

,

Иг

QT

 

Первый

.

 

(хі)

 

 

собой приведенное

из комплексов

представляет

время работы двигателя,

рассчитанное

по

значениям Ло и т|0.

Второй комплекс («/) представляет собой некоторую характери­ стическую скорость.

Решая уравнение (4.45) относительно ео, с использованием

введенных обозначений, получаем

 

* U f X f

(4.46)

ф и / Т / +

L

Комплекс tyuftf представляет собой некоторую характерную длину Lf. Заметим, что при проектировании заряда РДТТ при' заданных времени горения (толщине свода) и параметре Побе­ доносцева % формула для определения ео принимает следующий вид (см. работу [42]):

 

 

t-Єі

 

 

(4.47)

Легко заметить,

что

формула

(4.46)

по

структуре

подобна

формуле (4.47), с

тем

лишь различием,

что

в последней роль

характерной длины играет комплекс у.е\ (или хыС рт).

 

Используем далее зависимость,

определяющую возможность

размещения требуемого запаса твердого компонента

топлива

при принятой длине заряда L в летательном

аппарате

опреде­

ленного калибра:

 

 

 

 

 

 

 

 

^=FK.cns0QTL.

 

 

 

(4.48)

Разделив обе части

равенства (4.48)

на Q0 и

подставив

в него значение L из формулы

(4.46), получим

 

е 0 =

1

3°5о

.

(4. 49)

 

 

BlQrLfIwnD2

 

 

Формулы (4.41), (4.44) и (4.49) образуют систему уравнений для решения поставленной задачи. Последовательность решения представляется следующей:

1) рассчитав предварительно

£ М с) pc0IIQUft

и определив по принятому значению п соответствующее ему зна­

чение пВ2

(см. табл. 4. 1), по формуле

(4.49) определяем

е0 ;

2)

по

формуле (4.41), подставляя

в нее

полученное

значе­

ние єо, находим

DK,C;

 

 

 

 

 

3)

находим длину заряда L :

 

 

 

 

 

 

 

L =

Lf±=^-;

 

 

4)

находим

массы отдельных

элементов

камеры сгорания;

5)

находим

суммарную

массу

твердого компонента сот

= псот ь

а затем и массу жидкого

компонента:

 

 

 

 

 

 

 

і —

Ї

-

 

 

 

 

 

№... =

— ш т ;

 

 

6) по массе жидкого компонента определяем массы топлив­ ного бака, системы топливоподачи и арматуры;

7)определяем стартовую массу ЛА Qo как сумму масс его основных элементов;

8)проверяем, обеспечено ли получение заданных для проек­ тирования значений |1К , По и ц0. При несовпадении полученных значений с заданными повторяем расчет сначала, задаваясь другим значением п. Из вариантов, обеспечивающих удовлетво­

рительную

сходимость

по

указанным

параметрам,

выбираем

тот, который обеспечивает

наименьшую

стартовую

 

массу ЛА.

Случай

П.

Работа

двигателя

при

постоянном

 

давлении

(р = const)

и

постоянном

соотношении

расходов

компонентов

(г|) = const),

при

двойном

 

впуске

в камеру сгорания

жидкого

компонента

(см.

рис. 1.4).

 

 

 

 

 

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

Условие постоянства расхода твердого компонента выра­ жается равенством

wd = const, которое можно записать в виде

ud=uKDK.c. (4.50)

Подставив выражение (4.50) в (4.28), после интегрирования получим

„,

М<.сЕ0

 

 

4 а к

 

Рассмотрим два варианта

решения этой

задачи: а) при

Р>0,5 и б) при р<0,5.

 

 

а) При 0>О,5 для соблюдения равенства

(4.50) требуется,

чтобы расход жидкого компонента в канале Gm.T

возрастал, обес­

печивая необходимый рост скорости газификации твердого ком­ понента.

Обозначим отношение

расхода

жидкого

компонента,

пода­

ваемого в канал заряда,

к общему

расходу

этого

компонента

символом ср:

 

 

 

 

 

Очевидно, при р>0,5

коэффициент ср в процессе

работы

дви­

гателя должен возрастать, достигая значения, близкого к еди­

нице в конце работы.

Для

простоты

выкладок примем

срк =1-

Из условия ud=const

 

следует, что

 

 

 

 

С І .

:const=0j>,

(4.51)

(1

-

£ ) Є - 0 '

5

 

 

поскольку для конца

работы

при е = 0,

Gm,T=Gm.

 

Отсюда получаем

 

 

 

 

 

 

 

 

Р—0,5

 

 

 

 

<р =

( 1 - е ) Р

 

Для начала работы

(є = єо)

 

 

 

 

 

р—0,5

 

 

 

 

< Р о = ( 1 - в 0 ) Р

 

Расход газов в конечном сечении заряда в начальный момент времени определится так:

Ое0=АГкл{1-в0)д(\е)Рс0. (4.42а)

Здесь обозначения те же, что и в уравнении (4.42). Исключение составляет коэффициент Ж, который рассчитывается для состава

газов в сечении «с», отличающегося от состава конечных про­ дуктов сгорания перед соплом.

 

О с о = О т + ? 0 0 « = «

 

(1 + < P o - L ^ 1 ) •

(4 - 5 2 )

Поскольку

Gso=^o//io",

из

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

(4.26),

(4.42а)

и (4.52) получаем

 

 

 

 

 

 

 

 

 

 

М _ _ 0 Л ^ .

. ^ о _ ф / 1

+ ? о _ 1 ^ \ .

( 4 . 53)

По

аналогии с выводом

зависимости

(4.46)

находим

 

 

 

 

 

 

4

+

^

 

 

 

(4.

54)

 

 

 

 

 

 

 

 

 

где

Lh=

АК^Р^юн

 

 

J

_

 

( 4 > 5 5 )

Зависимость

для

&о получается

такой

же, как

и в предыду­

щем случае [см. уравнение

(4.46)].

 

 

 

 

 

 

Последовательность решения

задачи остается

такой

же,

как

и в предыдущем

случае.

 

 

 

 

 

 

 

 

Различие в расчете основных параметров определяется раз­ личием формул для определения £>к.с и характерных длин —

Lf и Ljtf.

Рабочая зависимость для определения DK.C получается из уравнений (4.30) и (4.50) посредством преобразований, рас­ смотренных в предыдущем случае:

'к.с

К(Ш

(4.56)

 

 

б) При р<0,5 для соблюдения равенства (4.50) требуется, чтобы по мере увеличения диаметра канала заряда расход жид­ кого компонента в канале уменьшался, т. е. коэффициент ф дол­ жен снижаться от единицы в начале работы до минимального его значения фк в конце. Из условия ud—const следует:

 

Ж.П

 

( 1 — с / " 0 ' 5

 

поскольку при е = е0 , Gm.T=Gm;

при е = 0,

Gm.r=G»

отсюда

 

 

 

Р—0,5

 

<Р« = ( 1 - * о )

" •

Расход газов в конечном сечении заряда в начальный момент времени будет равен:

Gco~ GTo + GmT— Gi,

т. е. такой же, как в случае I .

Следовательно, зависимости (4.42), (4.44) и (4.49) для слу­ чая I оказываются пригодными и для расчета данного варианта.

При этом длина Lf и коэффициент А должны определяться

так

же, как и в случае I .

 

 

Рабочая формула для определения диаметра камеры

DK,C

выводится из зависимости (4.29), в которую подставляется

зна­

чение ик , равное:

 

 

v ^ к . с ;

fic

 

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

 

/ Л . с = 4 т * - Ч

(1-ф)м-к П0

1

 

"

/ ,

ч0.5-3

 

 

 

Л-(р.ф)

( 1 ~ е о )

 

.

(4.57)

 

 

пЪв1

 

 

 

 

 

 

 

 

 

 

«о

 

 

 

Случай

I I I . Работа

двигателя

 

при

переменном

давлении

( p K = v a r ) , регулируемом

 

по определенной

 

программе,

и

постоян­

ном соотношении

расходов

компонентов

(\|)=const), при

двойном

впуске

в камеру

сгорания

жидкого

 

компонента.

 

 

 

При

проектировании

ГРД могут

быть

выбраны

различные

законы изменения давления и силы

тяги — в соответствии с за­

данной програмой полета летательного аппарата.

 

 

При любом варианте программного регулирования давления,

если i|) = const,

сохраняется прямая

пропорциональность

между

давлением

и

расходом

продуктов

сгорания, а

следовательно,

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

 

РкО

 

 

РкО

 

 

Второе уравнение можно переписать в виде

 

 

°т _

П(Е)Д _

 

рк

(4.58)

Ото

П (є0) и0

 

рк0

 

 

 

Подставляя уравнение

(4.58)

в (4.28), получим

 

 

 

Гк.с

 

(

ds

(4.

59)

 

 

 

 

(AW/AO)

яДссО -

£ о)0 , 5

"iCJL

)

 

 

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