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

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

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

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

Добавлен: 09.07.2024

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

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

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

числа М (ß = l,4 ), представлены на рис. 30 в сравнении с экспе­ риментальными данными [78].

 

 

Таблица 26

 

 

Таблица 27

 

 

а/а*

k

(0°

Р2 1

м

 

 

k = l ,4

А’= 1 ,25

 

 

 

 

 

 

 

 

 

 

1,2

14

2,36

0

1

1

1.4

13

2,5

1

0,92

0,94

1,67

12

2,7

2

0,78

0,83

 

 

 

3

0,68

0,73

 

 

 

5

0,58

0,62

 

 

 

При фиксированном значении числа М отношение давлений в косом скачке увеличивается при увеличении k :

 

 

 

 

Р-2

^

I _д_

fr>M2

 

 

 

 

 

Pl

~

1 / М 2 — 1 '

 

С другой стороны, угол при вершине турбулентной области

при увеличении к уменьшается

(см. табл.

26). Вследствие этого

 

 

 

 

 

 

отношение рг/р\ в скачке, отходящем

 

 

 

 

 

 

от точки отрыва, мало зависит от k:

 

 

 

 

 

 

отношение рг/рі слабо растет с уве­

 

 

 

 

 

 

личением k (табл. 27 при М = 3).

 

 

 

 

 

 

Таким образом, из условия рав­

 

 

 

 

 

 

новесия системы, состоящей из косо­

 

 

 

 

 

 

го скачка и зоны турбулентного сме­

 

 

 

 

 

 

шения, приближенно вычисляется пе­

 

 

 

 

 

 

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

 

 

 

 

 

 

ходящем

от

точки турбулентного

 

 

 

 

 

 

отрыва сверхзвукового потока.

Рис. 30. Отношение давле­

 

На основании полученного ре­

 

зультата

рассчитываются место от­

ний в скачке, отходящем от

 

рыва потока от стенок dT/d Kp и тяга

точки

турбулентного

от­

 

 

 

рыва:

 

 

 

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

/ —опыт

[78];

2—расчет

при

по­

 

Rr. При этом необходимо учесть, что

вороте

в скачке на

12°

 

 

в зоне свободной турбулентности,

чение

дозвуковое

и

 

 

примыкающей к стенкам сопла, те­

градиенты

давления

относительно малы.

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

ва) до ра происходит полностью в косом скачке

(см. рис. 29):

=

^ ---- = - ^ .

(103)

Pt

РФ (Мт) Р\

 

92


Величина числа Мт в месте отрыва определяется как корень трансцендентного уравнения (103) по известному отношению pjpo, табулированной газодинамической функции л(М) =

= [1+(А — 1)М3/2]~*/А_1и зависимости перепада давления в косом скачке от числа Маха р2/Рі = /(М). По числу Мт определяется площадь сечения, проходящего через точку отрыва,

q (Мт) ’

где

*4-1

_____1______ 2 ( * - 1)

q (М) —М

I + М2(/ѵ>— 1)/2

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

Расчеты показывают, что зависимость Fy/FKP=f(pJpu) может

быть приближенно аппроксимирована степенной Функцией

(£= 1,15-М ,4)

~

1 + 0,33й (fo- — — )0,G.

(104)

Ф<р

\ р н Якрі

 

Условие безотрывного течения на протяжении всего сопла имеет вид с/аМф^^тМф для заданного ро/рп-

В соответствии с предполагаемым ступенчатым распределе­ нием давления (см. рис. 29) тяга двигателя при работе сопла в режиме отрыва потока определяется формулой

Ят = /крфс^'кр2 W Ро—

( 105)

где М — приведенная скорость перед точкой отрыва.

В табл. 28 приведены значения коэффициента тяги идеально­ го сопла СН= Д / (роЛф) при d j d 1<v = 3; /г=1,4; рсрс = 1 в зависи­ мости от давления в двигателе ра и при рп=1,01 МПа с учетом

[выражение

(105)] и без учета отрыва [выражение (87)].

 

 

 

 

Таблица 28

 

Ро

rfT

 

СR

 

с учетом отрыва

без учета отрыва

 

Рп

da

 

(105)

(87)

 

30

2,6

1,4

1,34

 

20

2,2

1,3

1 , 1 9

:

10

1,7

1,2

0,74

 

5

1,4

1,1

—0,16

 

В действительности, возрастание давления от рт до рп про­ исходит не скачкообразно, а постепенно, на участке длиной

5

3734

93


« 1 0 6. Для относительных площадей FT/Fl<p и До.оэ/Дкр (относи­ тельная площадь Fojo/Fkp соответствует сечению, в котором дав­ ление достигает 0,95 рн) при отрыве потока в конических соплах

получены следующие

зависимости

[47]: F0,gs/FKpДт/Дкр =

= 0,416 с/7Т/Дкр — 1) •

при Дт/Дкр <

0,625 (Д0/ДКр) + 0,38;

F q,9s/FkP—F J F kp= Q,69 { F jF KP—1) при Дт/Дкр> 0,625 Да/Дкр+0,38;

(т. е. точка отрыва

находится на близком — меньшем

10 6 —

расстоянии от среза).

от Д0,95 до Fa давление изменяется

 

На участке сопла

мало:

от 0,95 рп до ри, и при расчете тяги можно использовать среднее значение давления: 0,975/7н(Да — До,95). Тяга конического сопла при этом определяется формулой

=/,<,,Tel*Д1іР~ (К) Ро+ 0,55 (Рт-і- 0,95/7,,) (До,95 - Fr) —

— /7,, (0 , 9 7 5 /^0,95 + 0 ,0 2 5 z 7 ,) .

Заметим, что на рис. 29 отрыв предполагается осесимметрич­ ным, однако он может также быть несимметричным, в частно­ сти, при больших углах расширения сопла или большом проти­ водавлении; в этом случае струя отклоняется от оси сопла, при­ соединяясь к стенке раструба.

Из изложенного видно, что в расширяющейся части сопла происходит частичное восстановление давления. Следовательно, в случае сверхзвукового сопла в критическом сечении устанав­ ливается течение со скоростью звука при давлении окружающей среды /7*, большем критического давления ркр= Пі<рро. Другими

словами, наибольшая (критическая) плотность потока наступает

при Рп= Ря^> ркр «0,53/70;

значение рнзависит от степени вос­

становления

давления

в

расширяющейся

части

сопла. При

лн = РнІРо ^

лн = Рн Ро

расход не зависит

от рп.

В области

я* <1 pJPo •< 1расход газа приближенно определяется с помощью формулы, аналогичной выражению (36):

Таким образом, тяга сопла, работающего в режиме отрыва потока, в первом приближении определяется формулой (105), структура которой совпадает со структурой обычной формулы тяги (90). Но в формулу (105) вместо параметров на срезе соп­ ла [Да, z(%a)] входят параметры в сечении, проходящем через точку отрыва [Дт, z (\T)I Положение точки отрыва dr/dKp прибли­ женно определяется аппроксимационным соотношением (104), полученным из условия равновесия косого скачка и турбулентной области, отходящих от точки отрыва. Восстановление давления в системе косых скачков внутри раструба при отрыве потока при­ водит к тому, что критическое течение в горле сверхзвукового

сопла устанавливается приян= р«/ро^>пкр-

94


3.6. НЕРАСЧЕТНАЯ СВЕРХЗВУКОВАЯ СТРУЯ

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

Построить качественную картину нерасчетной сверхзвуковой струн и рассчитать ее начальный участок можно методом одно-

Рпс, 31. Схема сверхзвуковой струн при ра> р п

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

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

Faq{la)= Fmq{lm)\

(106)

*ß») = *(XJ + ( ^

) ____ !___ ,

 

1 /крпУ О^а)

где п = ра/рц — степень нерасчетное™ сопла. При выводе урав­ нений (106) предполагались постоянными расход, температура торможения, полное давление на участке а т и сила от внешнего давления pn{Fm— Fa). При п=ра/рп— будет z(km)— кг(А,а), Т. е. импульс струи сохраняется (этого следовало

5*

95


ожидать, так как при рп— ИЗ исчезает сила внешнего давления на границе струи).

Система двух уравнений (106) позволяет найти две неизве­ стных переменных Х,„ и Fm по известным характеристикам сопла Ха, Fa и степени его нерасчетное™ п (для недорасширенной струи /г> 1). Уравнения (-106) являются трансцендентными, и их решение удобно проводить графически, откладывая по оси абсцисс отношение Fm/F„, а по оси ординат — функции z и q. Абсцисса точки пересечения линий z(Fm/Fa) и q(Fm/Fa) является искомым значением площади максимального сечения первой «бочки» Fm/Fa. Затем по первому уравнению системы (106) и найденному значению Fm/Fп находятся значения q(X,n), Хт и остальные осредненные параметры для этого сечения. В работе [1] приведены решения системы (106) при /г = 1,4 для различных чисел Ma=l-f-2,5 и различных степенен нерасчетное™ сопла п = = 2,5-ь50. При этом показано, что давление в максимальном се­ чении первой бочки независимо от числа Ма=1ч-2,5 определяет­ ся следующей линейной формулой:

Рт =

1

-f-0,07.

(107)

Ра

П

 

 

Формулу (107) можно использовать также для расчета всех осреднеиных параметров в сечении т пѵ. по формуле (107) на­ ходится рт, затем вычисляется ят—РтІРа и с помощью таблиц газодинамических функций определяются Хт, qm—FKVIFm, Fm и т. д. Таким образом, средние характеристики потока в сечении т т слабо зависят от числа Ма на срезе сопла.

Из формулы (107) видно, что давление в максимальном сече­ нии первой бочки ниже атмосферного рт/Ри< 1- Дальнейшее па­ раллельное течение газа невозможно, так как разность давлений

приводит

к искривлению линий

тока. Вблизи

кромок

сопла

граница

струи

отклоняется

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

угол,

определяемый

теорией веера

разрежения Прандтля— Майера

(гл. I).

ввести в рассмотрение

дополнительное

переменное —

Если

средний угол наклона скорости потока к оси, то в рамках одно­ мерной теории можно не только найти характеристики струи в се­ чении т т, но и построить профиль струи на участке а т. Этот метод расчета начального участка струи подробно изложен в работе [1]. Если степень нерасчетности струи близка к едини­ це, то сверхзвуковая струя вначале будет иметь почти периоди­ ческую структуру. Экспериментально установлено, что длина вол­ ны L этой почти периодической структуры осесимметричной сверхзвуковой струи определяется следующей формулой [8]:

96