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

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

.pdf
Скачиваний:
20
Добавлен:
23.10.2023
Размер:
5.5 Mб
Скачать

Сравнение эффективности различных органов регулирования вектора тяги можно провести по их качеству, т. е. по отношению управляющей силы к потерям тяги Ry/AR, и по результатам ве­ сового анализа всей ракеты. В табл. 31 приведены соответству­ ющие данные по ряду органов управления для мощной транс­ портно-космической ракеты при одинаковом уровне потребных управляющих сил [23].

3. 10. ОБТЕКАНИЕ ВЫДВИЖНОГО ЩИТКА

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

Возьмем часть стенки плоского сверхзвукового сопла, на ко­ торой установлен щиток высотой /г. Поворот сверхзвукового по­ тока осуществляется в косом скачке, отходящем от линии отры­ ва. Рассмотрим течение в окрестности отрыва (рис. 42). Погра-

в)

Рис. 42. Обтекание сверхзвуковым потоком препятствия на стейке:

о —схема течения при обтекании

выдвижного щитка;

б—распределение давления:

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

потока со струйным

препятствием

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

Из условия равновесия системы, состоящей из косого скачка и зоны турбулентного смешения, вычисляется перепад давления в скачке, отходящем от точки турбулентного отрыва сверхзвуко­ вого потока p2.lpi = f{Mi) ä; 1+0,5 Mi (см. § 3.5), где Mi — число Маха перед плоскостью отрыва хт.

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

109

ступенькой (различного вида), противодавлением, падающей ударной волной, поперечным вдувом газа или впрыском жидко­ сти.

Экспериментальные данные показывают, что давление при отклонении сверхзвукового потока перед ступенькой повышает­ ся не скачком, а достаточно плавно: на участке длиной -—2,5 ö

перед

плоскостью

отрыва

давление

увеличивается

до

» 0 ,6 (р2Рі) +Рі [36].

После

плоскости

отрыва, находящейся

на расстоянии ~4,5 Іг от ступеньки, давление продолжает плав­ но расти и достигает значения р2 в плоскости, расположенной иа расстоянии — 2/г от ступеньки. Непосредственно перед ступень­ кой давление снова возрастает и на ее лицевой стороне дости­ гает величины » 1 ,6 р 2- Это объясняется натеканием на ступень­ ку высокоскоростной части клиновидной области турбулентного размыва тангенциального разрыва и образованием у верхней кромки отошедшей ударной волны у. Граница области размыва составляет со стенкой угол ~ 16,5°. Над этой линией располагает­ ся более или менее однородная область А сверхзвукового тече­ ния. Линия, соединяющая точку отрыва с верхней кромкой сту­ пеньки, составляет со стенкой ~13°. Эта граничная линия отде­ ляет поток, проходящий через ступеньку, от того потока, который попадает в зону обратных токов и образует клиновидную об­ ласть, поворачивающую поток в отходящем от ее вершины ко­

сом скачке

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

p2jpі (р11с. 42, е).

Боковая сила,

вызванная

турбулентным

отрывом плоского

сверхзвукового

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

интегралов

прироста давления по области

перед

отрывом

(~ 2,5 б ) и по области отрыва ( — 4,5 /г).

 

 

 

Приближенно

 

 

 

Ry_

4,5 (рч— рi) hr (.МО

 

 

 

R

Р\Н

п

 

 

где R = piH/r(lAi) =РіН (1+kMr) — пустотная

тяга

двигателя

 

при

выведенном щитке;

 

Н — высота поперечного сечения

 

плоского сопла шириной=1;

 

Кѵ= 2,3 Мігі) — коэффициент относительной

 

боковой силы;

 

 

Рщ = Іі/Н-— относительная площадь по­

 

гружения щитка.

 

Потери тяги вследствие погружения щитка в поток (или от­

носительное

сопротивление щитка) без

учета

давления на его

тыльной стороне равны

дR __

1 ,%p2hr (М|) __ —

R

Р\Н

где /Сі = 1,6р2/[/Оіг(М1)]= 1,6 г(Мі) (1+ 0,5 M i)— коэффициент относительных потерь.

по

Коэффициенты относительной боковой силы Кѵ и относитель­ ных потерь тяги К а , а также отношение RV/\AR\ =Куа д л я плоского сопла со щитком на срезе приведены в табл. 32 в за­ висимости от числа Мі потока (/г=1,2).

 

 

 

 

Таблица 32

Mj

Т'д/Т'кр

Ку

Ка

/ э д д я і

2

1.9

0,80

0,55

1,5

3

6,4

0,59

0,34

1,7

4

28

0,46

0,24

1,9

Таким образом, относительная боковая сила и потери тяги при погружении плоского щитка_ пропорциональны относитель­ ной омываемой площади щитка Fm, а коэффициенты пропорцио­ нальности Кѵ и К а зависят от числа Мі и /г. Качество щитка, т. е. отношение Ry/\&R\, равно приблизительно 1,7.

Рис. 43. Схема течения перед цилиндром, установленным на пластине и облекаемым сверхзвуковым потоком:

а—вид

сверху;

б—пространственная

картина

течения;

/—косой

скачок

уплотнения; 2—отошедшая ударная

волна;

Л—скачок уплотнения в отклоненном потоке; 4—граница застой­ ной зоны; 5—линия отрыва потока; 6—хвостовой скачок уплот­ нения; 7—застойная зона за цилиндром; 8—линия, проходящая через точку максимального давления в застойной зоне

Выдвижной щиток РДТТ имеет не только верхнюю, но и бо­ ковые кромки [41]. Турбулентный отрыв потока при натекании на такое трехмерное препятствие имеет более сложный харак­ тер. В области отрыва возникают течения, параллельные линии пересечения плоскостей щитка и сопла. Трехмерные эффекты становятся особенно существенными, если высота щитка превы­ шает размах в 5 и более раз (рис. 43). Например, в случае обте­ кания достаточно высокого цилиндра, установленного на пла­ стине, длина зоны отрыва (расстояние от передней кромки ци­

111

линдра до линии отрыва) зависит только от диаметра цилиндра: /T» 2af«4r; максимальное повышение давления в отрывной зоне

Р-2ІРі ~ 1+0,5 Мі [17].

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

ментально. Опытные характеристики интерцепторной

системы

управления вектором тяги даны в табл. 33

[41].

 

 

 

 

Таблица 33

Относительная площадь щитка в выходном се­

0,05

0,10

0,15

чении сопла Fm

 

 

 

Угол отклонения вектора тяги

Относительные потери тяги 1R/R

0,0'2

0,05

0,08

Обтекание дефлектора аналогнпчно рассмотренной картине течения в сопле около выдвижного щитка. Однако поверхность дефлектора составляет отличный от нуля угол с направлением оси у (перпендикулярной к оси сопла). Поэтому равнодейству­ ющая сил давления на дефлектор является не только сопротив­ лением, но и имеет составляющую по оси у. Кольцо дефлектора воспринимает ~2/3 полной боковой силы, а ~1/3 боковой силы приходится на стенку сопла в области отрыва потока.

3. 11. СОЗДАНИЕ БОКОВОЙ СИЛЫ ВТОРИЧНОЙ СТРУЕЙ

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

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

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

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

Рассмотрим поток в плоском сверх­

 

звуковом сопле (с выходным сече­

 

нием высотой Я), а также воздей­

 

ствие иа этот поток вторичной струм

 

бесконечно малой интенсивности сЮ

 

(рис. 44). Размеры возмущенной об­

 

ласти ограничены линией Маха,

от­

 

ходящей от начала этой области под

 

углом а и заканчивающейся на про­

 

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

 

за. Длина области равна

Рис. 44. Взаимодействие вто­

1 = Н ctg а = Н ctg arcsin — =

ричной струн бесконечно ма­

лой

интенсивности с пото­

М

ком

в плоском сверхзвуко-

= н V м -- 1 .

Повышение давления в потоке вследствие ввода массы clG под углом ß к оси сопла оценим по теории одномерных пото­ ков [80, 93]:

 

dp

_

А'М2

f dh

[ 2 / j

I

к — 1 доз\

[___

 

 

р

М2— 1 I СрТ

 

\

 

2

j a

 

 

 

 

- [ 1 + ( / г -

1)Ма]

^

C0S^ G

N

I ,

 

 

 

1

1

 

;

'

 

vG

I

 

где /?

и р ■— энтальпия и давление потока;

 

 

ѵ = М.а — скорость потока;

 

 

 

 

 

 

 

ÜBT — скорость вторичной струи;

 

вторичной

струи;

увт cos ß — осевая

составляющая

скорости

 

N — молекулярный вес.

 

по возмущенной области

Равнодействующая

сил давления

длиной

l — H ] f

— 1 равна

dRyv = ldp,

а

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

удельный импульс определяется соотношением

 

 

 

dRy р

._

kpH М2

Gdh

=*(■

4 м Й' -

\

 

УР'

dG

G

'М2— 1

I cpTdG

 

 

_

[ l - L

( / ? _

1 ) M *

Уat COS ß

GdN I

 

 

 

 

 

 

 

 

V

 

 

NdG I '

 

Результаты расчета Iy для случая

вдува

в основной поток (k=

= 1,24;

7о=3000 К; М = 2,5)

различных газов и жидкостей при­

ведены в табл. 34 по данным работы [80].

 

 

 

ИЗ

 

 

Таблица 34

 

Тепло, поглощае­

 

 

 

мое впрыскиваемым

 

 

 

веществом при на­

 

 

Вдуваемое рабочее тело

греве от началь­

I,/, м/с

ЛЯ*.

ной температуры

 

до статической

 

кг/м3■с

 

температуры в по­

 

 

 

токе = 1700 К),

 

 

 

кДж/моль

 

 

Инертные газы

 

 

 

Н2

41,8

7300

0,06

n 2

46

2320

0,27

Продукты сгорания твердого топ­

—54

4260

7,1

лива (от= 1 , 6 6 г/см3)

 

 

 

Нейтральные жидкости

 

 

 

НоО

144

1240

1 , 2

HCI2 F2 (фреон-12)

156

1560

2 , 0

Вго

83,6

1700

5,3

Реагирующие жидкости и газы

 

 

 

2 содержится в потоке)

—426

6600

0,9

0 2=2Н оО—2Н2

—860

4700

6 , 8

N2О4 = N2“Ь4-Н2О—4-Н2

2 1 0

1470

2 , 2

N2 0 4 = 2 N0 2 ( д и с с о ц .)

 

 

 

* 6 — плотность рабочего тела. Газы хранятся при давлении 10,3 МПа

и впрыскиваются со звуковой скоростью под углом ß=90°

к оси

сопла.

Боковая сила при вдуве вторичной струи складывается из равнодействующей сил давления dRv р и собственной реактивной силы dRm = !BTdG.

 

riF

Здесь

/»т = г'втsin 3+ (Рв' — Р) du ,

где рвт— давление во вторичной струе; F-ач ■— площадь сопла вторичной струи.

В соответствии с этим удельный боковой импульс струйного органа управления равен (табл.35)

=/вТ'

Из результатов расчета Іѵ по изложенной линейной теории и сравнения с опытными данными (экстраполированными к зна­ чениям dG?xO) следует:

114

Таблица 35

 

Газовая постоян­

Удельный импульс

Коэффициент

Вдуваемый газ

ная,

усиления

Дж

вторичной струп

 

 

 

 

/цт> М/С

к

І у

 

кг ■К

К

у - j

 

 

 

1ВТ

н2

4160

2000

 

3,7

n 2

295

500

 

4,7

Продукты сгора­

314

1800

 

2,4

ния твердого топ­ лива

линейная теория дает удовлетворительное описание зави­ симости боковой силы от газотермохимических факторов;

удельный боковой импульс, возникающий при вдуве рабо­ чего тела в закритическую часть сопла, в несколько раз выше собственного удельного импульса вторичной струи /вт; коэф­

фициент усиления для газовой вторичной струи равен

(dG— М3)

=2,44-4,7;

'ПТ

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

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

Картина отклонения части сверхзвукового потока вторичной струей аналогична картине обтекания препятствия, установлен­ ного на стенке (см. рис. 42, в). Перед проходящей по всему пло­ скому соплу поперечной щелью, через которую вдувается вто­ ричная газовая струя, образуется отрывная зона с возвратными течениями. Повышение давления в косом скачке, отходящем от вершины клинообразной отрывной зоны, приближенно равно РгІРі = 1+ 0,5Mi, где Mi — число Маха перед линией отрыва [77].

Длина отрывной зоны, как и прежде, пропорциональна высо­ те проникновения вторичной струи: Ita (4,2-=-4,5)/гСТр- Высоту проникновения /гстр можно определить из равенства количества движения струи в поперечном направлении ее сопротивлению основному потоку, т. е.

°вт'г'иг= [(!+ ?) Р г ~ аРЛ Астр.

где (l+ ß )p 2 — среднее давление на переднюю грань контроль­ ной поверхности, ограничивающей вторичную струю, причем

і + ( й — і)М?/4

115

api — среднее давление на заднюю

грань ( а ^ І ) .

 

Таким образом, для суммарной боковой силы при вдуве газа

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

 

Ry = Ry и

= 4,2 ( /;2 — Л ) Лстр -f А>„,

2 , 2 ДА х О ң ' У ң г _____

+ ^вт>

ß)(l 4- ÜjöMj)—а

 

 

 

где Rur — тяга вторичной струп.

Так как в пустоте (RnT)cc = k+\/k[GaKPz(K)]ßT, для коэффици­ ента усиления в пустоте KU= RUI (Rin)** имеем

д- =

______________ 2,2Л']/сЛпт_______________ I

/г„т

 

у

(А- + 1) [(1 + Р) (1 -I- 0 , 5 Д \ ,) - а] г (Х„т)

(/?,„)„

'

Расчет коэффициента усиления по этой формуле

при вдуве

азота в воздушный поток дает значения /\?у= 2,94-3,2

при Мі =

=2-^6, согласующиеся с данными специальных опытов.

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

Результаты экспериментальных исследований взаимодействия воздушного потока с поперечной газовой струей, вытекающей из отверстия, проанализированы в работе [19] с позиций теории

подобия. При этом предполагалось, что влияние кривизны сте­ нок мало, а зоны перераспределения давления, обуславливаю­ щие возникновение дополнительной боковой силы, полностью умещаются на стенке сопла. В широкой области значений па­ раметра Gmvj(pL2) удалось провести корреляцию опытных зна­ чений коэффициента усиления:

К у=

-

1,87 lg- ^ + 1 , 1 2 ,

 

 

 

p L °-

 

где V и р — скорость и давление основного потока

(в месте

вдува);

от

критического сечения до

центра от­

L — расстояние

верстия вдува.

В исследованном диапазоне параметра GnTv/(pL2) =0,06-1-0,9 коэффициент усиления изменяется от 3 до 1.

Влияние трехмерных и нелинейных эффектов приводит к то­ му, что удельный боковой импульс / у и коэффициент усиления Ку уменьшаются при увеличении относительного расхода вторич­ ной струи GBT/G (табл. 37). Относительная боковая сила Ry/R сначала растет с увеличением GnT/G, достигает максимума при некотором его значении, а затем уменьшается при дальнейшем

116

Таблица 36

 

РДТТ; G =

 

 

 

 

 

 

РДТТ ракеты

 

 

 

=•■16,7 кг/с;

ЖРД;

G =7,25

кг/с;

 

 

 

 

 

 

„Титан-ЗС“;

 

 

 

вдув из ПГГ

вдув из ЖГГ

[40]

 

 

 

 

 

 

впрыск N0O4 [64]

 

 

[40]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gbt/G

0 , 0 1 2

0,03

0,03 0,05

0,07

0,09 0,18

0 , 0 2

0,05

0 , 1 0

0,15 0 , 2 0

0,30

Ry/R

 

 

 

 

 

 

 

1

 

 

 

0,115 0,14

0,023

0,05

0,05 0,065 0,077 0,07 0,045|0,02 0,0-15 0,075 0 , 1 0

ly, м/с

4300

3600

3000 2400

2 1 0 0

1500

470

2600 2400

2 0 0 0

1700 1600

1 2 0 0

Ку

2 1 , 6

1,4

0,9

0,3

росте

GBT/G.

Например,

для

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

в табл.

36 ЖРД

(Ry/R)max= 0,077 при GBT/G = 0,07. Немонотонный характер зави­

симости R,j/R= f(GrsT/G) объясняется тем,

 

 

 

 

 

что при слишком большой

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

 

 

 

 

 

вдува возмущенная

область захватывает

 

 

 

 

 

газовые слои вблизи противоположной по­

 

 

 

 

 

верхности сопла, и давление повышается

 

 

 

 

 

на таких участках сопла, которые дают

 

 

 

 

 

боковое усилие в направлении вторичной

 

 

 

 

 

струн. Наибольшее значение боковой силы

 

 

 

 

 

получается

тогда,

когда

возмущенная

Рис.

45.

Схема

тече

область занимает приблизительно

поло­

мня

при

наибольшем

вину проходного сечения сопла

(рис. 45).

допустимом

интен­

На рабочем участке для зависимости

сивности

вдува:

Ry/R = f (GBT/G) может быть использована

/ —у д а р н а я

в о л н а ;

2—

линейная либо степенная

аппроксимация.

в о з м у щ е н н а я

об л а с т ь ;

 

5—м есто

в д у ва

В частности, при вдуве фреона Rv/R tt 0,35,

 

 

 

 

 

GBT/G

[90];

Ry/R=0,3 (GBT/G)0’5'

для

ЖРД

[40];

R,JR =

= 0,4

( G b t / G )

0 ' 7 2 для РДТТ ракеты «Титан-ЗС» [64].

 

 

 

Высокая эффективность системы впрыска в ракете «Титан-ЗС» обеспечи­

вается

конструкцией форсунок,

порядком и местом

расположения

их

в сопле

( 4 н р =

0,96 м,

(/„ = 2,71

м), а

также

подбором

впрыскиваемой

жидкости.

Впрыск N2O.1 осуществляется в каждом квадранте через 6

форсунок под дав­

лением 5,25 МПа, диаметр проходного сечения в месте впрыска 1,8 м. При

полете непрерывно определяется избыток N20.i, который сливается через

все

24 форсунки, не создавая боковой силы. Масса NjO/, рассматривается

как

прибавка к топливу с собственным импульсом 100=120 [64].

 

Г л а в а IV

НЕУСТАНОВИВШЕЕСЯ ТЕЧЕНИЕ ГАЗА В РДТТ

Неустановившиеся течения газа осуществляются в пороховом двигателе в периоды воспламенения заряда, отсечки тяги и по­ летного регулирования тяги. Задачи об изменении характеристик РДТТ при неустановившнхся режимах работы являются нели­ нейными даже в квазистационарном приближении.- Для их реше­ ния применяются численные методы, а также приближенные, обладающие различной степенью сложности и точности [9, 62, 65, 69, 75]. Полученные аналитические решения для переходных про­ цессов в двигателе при выходе на режим, отсечке тяги и перете­ кании газов по соединительному газопроводу сообщающихся двигателей просты и достаточно точны. Время неустаповившегося истечения в рассматриваемых случаях порядка ІО-2 времени работы РДТТ. Изменение геометрических характеристик двига­ теля (свободного объема и площади проходного сечения) за это время того же порядка малости (10~2) и в дальнейшем не учи­ тывается.

4.1. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В РДТТ ПРИ ВОСПЛАМЕНЕНИИ ЗАРЯДА

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

Т ( x ) ^ T 3Jr (Fs~ T 3)е- д'"/а,

где Ts, Т3— температура поверхности горящего топлива и на­ чальная температура заряда;

ü = X / ( c q ) — коэффициент температуропроводности; X— расстояние от горящей поверхности.

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

со

cpQ ( T - T 3) d x ^ ± ( T s- T 3).

о

Основной запас этого тепла заключен в слое толщиной б= а/и, время прогрева которого порядка 4 = б/« = а/и2 (времена тепло­

118

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