Файл: Шишков А.А. Газодинамика пороховых ракетных двигателей. Инженерные методы расчета.pdf

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

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

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

Добавлен: 09.07.2024

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

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

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

потока в выходном сечении); £тр— потери из-за трения;

потери из-за химической иеравновесности сечения;

потери из-за многофазности;

£пр — прочие неучтенные потери.

На рис. 24 показаны профиль 1 дозвуковой части сопла Витошннского, рассчитываемый по формуле

г

( і + ^ О а

где Бд. ч — длина дозвуковой части, и профиль 2 радиусного до­ звукового сопла, состоящего из грех элементов: а) области входа

Рис. 25. Сверхзвуковая часть сопла:

 

/—промежуточная

линия тока

с радиусом

б) конического участка с углом 0ВХ и в) области

критического сечения с радиусом г2При этом должно быть

 

П О .« : 0„х<45°;

гкр < г2< 2гкѵ.

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

Коэффициенты расхода сопел Витошинского и радиусных со­ пел при г2>0,5 /'кр близки к единице.

Минимальную длину сверхзвуковой части ( L c . ч) имеют сопла с угловой точкой К в критическом сечении и равномерным по­ током на выходе (рис. 25) [44].

Быстрый разгон потока до скорости, соответствующей числу Мц, от плоской звуковой поверхности ОК осуществляется в цен­ трированной волне разрежения ОКВ, вызванной обтеканием уг­ ловой точки К. Выравнивающий участок КС переводит течение из свободнорасширяющегося на характеристике КВ в однородное и параллельное на характеристике ВС. Расход газов через кри­ тическое сечение ОК и характеристику ВС одинаков.

73

Газодинамические параметры в области отхода от звуковой поверхности рассчитываются аналитическими методами. Даль­ нейший расчет как в области ОКБ, так и в области ВКС про­ водится методом характеристик [44, 45]. Веер волны разрежения ОКБ остается одним и тем же для всех контуров при заданных свойствах газа (изменяется только положение характеристики КВ соответственно значениям числа М„). Таблицы осесимметрич­ ного сверхзвукового течения свободнорасширяющегося газа с плоской звуковой поверхностью для /г= 1,14; 1,33; 1,4 и 1,66667, рассчитанные на быстродействующей электронной счетной маши­ не (БЭСМ) АН СССР, приведены в работе [45]. Исходное одно­ родное звуковое течение не всегда можно состыковать с под­ ходящим дозвуковым потоком, но это приводит к незначитель­ ным потерям тяги, если дозвуковая часть выполнена достаточно гладко.

Однако идеальные сопла с угловой точкой получаются слиш­ ком длинными и тяжелыми, например: Lc. ,,= (5-f-6)r„ для газа А’=1,25 при Ma= 3-h5 [63]. Вместе с тем прирост тяги на конце­ вых участках аС настолько мал, что он может погаситься мест­ ной силой поверхностного трения. Следовательно, укорочение идеального сопла не только уменьшает длину и площадь поверх­ ности сопла (и, следовательно, массу), но и увеличивает тягу.

Таким образом, оптимальное сопло может быть построено на основе идеальных сопел, укороченных до выходного сечения а. Этот метод профилирования сверхзвуковых частей является до­ статочно простым, охватывает различные задачи на экстремум (минимальные длина, поверхность или диаметр выходного сече­ ния) и дает результаты, близкие к результатам расчета опти­ мальных сопел методами вариационного исчисления (метод Гудерлея-—Шмыглевского) [5, 63]. В плоском случае укорочен­ ные идеальные сопла п сопла минимальной длины совпадают.

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

ного идеального

сопла осуществляется для постоянного k =

= ln ^2- /

ln — 1

равного k в выходном сечении сопла [63]. По

Po I

Qn

 

технологическим и эксплуатационным причинам иногда целесо­ образно вершину угла скруглять. В этом случае профиль может быть выбран по промежуточной линии тока 1. Примеры укоро­ ченных идеальных сопел для /г=1,2 приведены в работе [5].

Коэффициент потерь на непараллельность истечения в про­ филированных соплах приблизительно равен [58]

где Ѳ„ и Gi — полууглы расширения сопла в выходном сечении и в начале раструба.

74


Для профилированных сопел коэффициент потерь импульса на рассеяние составляет сг ~ 0,54-2%; коэффициент потерь им­

пульса на трение £тр« 4 6 ГA4 —1-4-2% — толщина потери импульса в пограничном слое на стенке сопла в выходном сече­ нии).

В пороховых ракетных снарядах по компоновочным, техноло­ гическим и др. причинам применяют также конические сверхзву­ ковые сопла, в которых 2Ѳ„= 204-30°. Потери на рассеяние в ко­ нических раструбах с точностью до 0,01 оцениваются на основе модели радиального течения, из которой следует:

При Ѳп=10°,

15° и 20° значения

соответственно равны 0,008;

0,017 и 0,030.

 

 

Конический раструб для сопел с большой степенью расши­ рения (с/яМіф>2 ) длиннее и тяжелее профилированного при оди­ наковом уровне потерь.

Современные высокопмпульсные твердые топлива содержат добавки металлов, вводимых в виде тонкодисперсного порошка. Наиболее широко применяются добавки алюминия (5-1-15%); изучаются топлива с добавками бериллия [4, 68]. Продукты сго­ рания таких топлив содержат значительное количество конденси­

рованных окислов металлов: А120 3, ВеО,

Ве20 3.

Окислы находят­

ся в

камере РДТТ (7 = 30004-3500 К)

в виде

жидких капель

либо

в виде твердой пыли.

 

 

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

низированных топлив

найдено эмпирически: £с= 0,02 + 0,16 е

[68], где е — процентное

содержание конденсированной фазы в

потоке.

При профилировании сопел РДТТ следует учитывать особен­ ности гетерогенных потоков: существенную роль сверхзвуковой части сопла, траектории частиц в до- и в сверхзвуковой частях (возможность выпадения конденсата на стенку), непараллель­ ность и неоднородность полей скоростей обеих фаз потока в кон­ трольных сечениях сопла (■§ 3.4). Например, при применении сопел, частично вдвинутых в камеру ракетного двигателя, рабо­ тающего на смесевом твердом топливе, возникают дополнитель­ ные потери тяги (табл. 16) [58, 81].

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

75


 

 

 

 

Таблица I

Степень погружения в %

 

0

3

25

50

Содержание АЬ03 в

21,5

0

1

1

1,2

 

ДО

%

0,4

0,4

0,5

 

5

0

сеяние (£р)кв=1— (sin 0а/Ѳа ) 2

и появления дополнительных

по­

терь на неоднородность течения в двугранных углах £Пр~0,5% . При несимметричном выходном срезе сопла распределение давления на боковых стенках не изменяется, если линии Маха,

 

отходящие от выходных кромок, по­

 

падают на свободные границы исте­

 

кающей струи, и оно используется

 

для

оценки боковой силы, воздейст­

 

вующей на неуравновешенную часть

 

раструба.

Площадь боковой

проек­

 

ции косого среза сопла, наклонен­

 

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

 

речного

сечения

(рис.

26),

равна

 

FV = F sin a = Fatg a (1—tg20 atg2 a )-3/2,

Рис. 26. Сопло с косым

а

боковая

сила

на этом участке

срезом:

Rvm p aFy

(сечение

а а проходит

м 0 = 3 ; * = 1.25; р „ = 0

через точку

пересечения

плоскости

 

косого среза с осью сопла).

Таким

образом, отношение боковой

силы

Ry^ p aFv

к осевой

Roa =

= paFJr(Ka) =paFa(l+kNi„'1) в пустоте

 

 

 

 

 

__________ tga_________

(89)

( 1 +АМI)

(l-tg->0o tg->a)3/'2

 

Результаты расчета Ry/R«, по этой формуле и при интегри­ ровании давления р— пЩро на неуравновешенной части сопла практически совпадают (табл. 17).

Величину R,j/Roc можно также получить, исходя из уравнения количества движения [62].

Площадь критического сечения сопла Fup, геометрическая сте­ пень расширения (da/dKр) 2 и тип контура сопла используются в расчете основных характеристик РДТТ:

массового расхода продуктов сгорания в единицу вре­ мени G;

тяги в пустоте Rac',

эксцентриситета реактивной силы (см. '§ 3. 8 ).

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

76


 

 

Q _Лж'р^кр

 

 

 

 

1

_

с*

 

 

R ™ =

G d Kpt^ ( ^ a ) =

f 'K

f P

\ ) F ( \ j ) =

Ро^ a ff t'" I'

 

 

k_______________

 

 

 

где С* = 1у / RT0j k

{ j ^ [ f k+mk

l) = v R T 0l{mKV) = A~1 — xap ак-

терпстическая

скорость;

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

17

 

 

 

 

 

R y / R

 

 

а 0

“в

П о ф о р м уле

(8 9 ) при

П р и численном

 

 

р а

 

 

М 0= 3 ,

/ г = 1 ,25,

р „ = 0

и н тегри рован и и

 

0

0 ,0 0 7 1 4

 

0 ,0 0 7 1 4

 

 

15

0 ,0 0 7 1 5

 

0 ,0 0 7 1 2

 

20

0

0 ,0 2 9 7

 

 

0 ,0 2 9 7

 

15

0 ,0 3 0 1

 

 

0 ,0 2 9 3

 

 

 

 

 

30

0

0 ,0 4 7 1

 

 

0 ,0 4 7 1

 

15

0 ,0 4 8 9

 

 

0 ,0 4 6 1

 

 

 

 

 

Действительный расход газа через сопло отличается от рас­ четного, соответствующего равновесному одномерному потоку при тех же значениях (RT0)KP, р 0 кр и АКР и плоской поверхности перехода от дозвуковых скоростей к сверхзвуковым. Это откло­ нение характеризуется коэффициентом расхода сопла р. В соп­ лах с резким переходом от дозвуковой части к «горлу» (малый радиус скругления /'2< 0,5 гІф, большой угол Овх конического до­ звукового сопла) звуковая поверхность искривляется, может произойти местный отрыв потока от стенок сопла. Эти наруше­ ния одномерности потока приводят к наиболее существенным отклонениям р от единицы; при этом в рамках модели кусочно однородного течения газа коэффициент расхода совпадает с ко­ эффициентом сжатия струи (см. § 3. 2 ).

При известных р1Ърк, ц и р расход газа через сопло и реак­ тивная сила определяются отношениями (установившееся исте­ чение)

■'ІР-АДкр

с

R = K Tr\ppKF KP — p nF a.

77


Здесь

Кт — коэффициент

тяги

(приближенно

/Ст»

~срс/кРг (1 я), где ф0 — коэффициент сопла);

ц = р 0іф/Рк — коэф­

фициент восстановления давления в РДТТ

(см. гл. II).

 

Влияние характеристик внешней среды на тягу соплового ап­

парата не

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

pnFa в формуле

(90).

При увеличении ри/ра> 1 возможно изменение режима газового течения в раструбе — отрыв потока (см. § 3.5). В условиях по­ лета ракеты при уменьшении p j p a< 1 возникают сложные явле­ ния взаимодействия реактивных струй с внешним потоком и между собой (в случае многосоплового блока). В результате таких взаимодействий возможно изменение донного и лобового сопротивлений ракеты [9].

3.2. КОЭФФИЦИЕНТ РАСХОДА СОПЕЛ

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

с коэффициентом расхода степенной зависимостью /і~!і~'І/<1“ ѵ) , где показатель степени 1/(1— ѵ )> 1 . Поэтому относительная ошибка в оценке и приводит к большей относительной ошибке в определении давления: Ар/р = —[1/(1 — ѵ)]Д

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

Рассмотрим докритическое истечение газа через внезапное сужение (рис. 27). Уравнение количества движения для сечений 1—1 и 3—3 можно записать в виде

А п / I х /F i =

Роз/

F 3“Г{А п — АР) (F i — Г 31

(91)

Импульс сил давления на торцовую стенку равен

 

(А)і — А р ) ( К 1F 3) = Ап —О —Сі)

 

А п

(1

^ i) A n

U-

 

где коэффициент смягчения входа у

из-за уменьшения площади

проходного сечения принимается таким же, что и при течении

несжимаемой жидкости:

]/'^1 = 0,637.

Подставляя это выражение

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

движения

(91), используя приближен­

ное уравнение неразрывности %iFitak3F3 при оценке малых чле­ нов, получим

Роз/ (h) F 3~ Po1 (l + ^ _ Xi) F3- ^ - J L - PoiltF3(1 - ^ - ) •

78