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

книги из ГПНТБ / Москаленко Г.М. Инженерные методы проектирования в ракетодинамике

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

ставной частью (чередующейся фазой) рикошетирующего полета.

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

гаться по круговой орбите

на различных высотах в широком

диапазоне скоростей

Ѵ к р

gR.

Полет со скоростью, превы­

 

 

шающей первую космическую скорость, может происходить при отрицательной подъемной силе, направленной к центру Земли (см. рис. 192). В главе приведены некоторые характеристики суборбитального полета ракетоплана по тяге, расходу топлива, дальности и времени полета.

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

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

13.1. ДВИЖЕНИЕ РАКЕТОПЛАНА В РЕЖИМЕ РИКОШЕТИРОВАНИЯ

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

Рис. 180. Типовая траектория полета с отражением

Чередующаяся последовательность баллистических и рико­ шетных фаз определяет нестационарность режима рикошетиро-

31Т

вания. Минимальная высота рикошетирования всегда остается меньше высоты квазистационарного планирования.

Уравнения движения ракетоплана в рикошетной фазе могут быть записаны в следующем виде

 

dt

 

X — G sin Ѳ,

 

 

 

M V — - Y — Q

cos I

gR

(13.1)

 

d H

 

dt

 

 

 

dt

 

■ V

 

 

 

 

 

d Z.

=

1 г sin Ѳ. ,

 

 

 

 

----

V

cos Ѳ,

 

 

 

где Ѳ — угол

dt

 

 

 

траектории

к местному

горизонту;

наклона

М

 

 

 

 

(постоянная в полете).

 

— масса ракетоплана

 

Расстояния, проходимые ракетопланом в рикошетной фазе

сравнительно

невелики и составляют порядка 4- ІО2

км. По этой

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

имеет место (тэіпѲ-СУ и G cos0<cK

В режиме

рикошетирова­

ния обычно Ѳ<10°. При этом можно

положить

также sin 0=0

и cos 0 = 1.

 

 

 

 

При сделанных допущениях уравнения движения (13. 1) упро­

щаются и принимаютd вид

— g

V \2

e -ß tf,

 

у

 

 

 

dt

 

v :

 

 

V db dt

;i3.2)

* * L = V B , dt

dL

■ V,

dt

где

2G/S 11/2y ß

Г 2 G/S 11/2

CxQo 1 •

vT ] ■

Как показывает интегрирование уравнений (13.1) численным методом, изменение скорости на участке рикошетирования не превышает ~ 8% . Это означает, что интегрирование уравнений (13.2) можно выполнить при допущении K=const. Поскольку при больших гиперзвуковых скоростях полета зависимость аэро-

3 1 8

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

V

невелика, полагаем

также

сх =

const

и

сѵ —

 

имеется

в виду, что вели­

 

 

 

const. Везде

чина

g

 

отнесена

к условиям

на

средних

 

высотах отражения.

Однако,

принимаяg

= во внимание

малые

 

высоты отражения

(^min<100 км), при выполнении практических расчетов целесо­

образно положить

 

go-

 

 

 

 

Найдем распределение углов в рикошетной фазе траекто­

рии. На

основании

предпоследнего уравнения системы (13.2)

имеем

 

 

 

 

 

d2H

у М

 

 

(13.3)

 

 

 

 

 

 

 

d fi

dt

 

 

 

После подстановки во второе уравнение той же системы получим

 

 

 

 

 

(РН _

 

 

і Ѵ е " .

 

 

 

(13.4)

 

 

 

 

 

dft

 

 

 

 

 

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

при выполнении начальных условий

Н = Н 0

и замены

dH/dt=VQ

о дает

 

 

 

е'-ряо) 1/2

 

 

 

 

 

 

VI

 

 

 

 

(13. 5)

Если положить в

 

р

 

(13.5)

Н = Н К = Н 0,

то

получим

формуле

 

 

Ѳк = —Ѳк. Это

означает,

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

входа и выхода в рикошетной

фазе

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

абсолютной величине и противоположны по знаку.

 

 

 

Минимальную

высоту отражения

Н — Нтіп

найдем,

положив

в формуле (13.5)

Ѳ =2

0.£

При/ J

этом_

Н„

-РЯо),

 

 

 

 

откуда

 

Ѳо =

 

Р

Уу

 

 

 

 

 

 

И min

 

 

 

■ J<»o

(1/;)2+ е - ря°

 

 

 

(13.6)

 

 

 

 

 

 

 

. 2g

 

 

 

 

 

 

Как видно из полученной формулы, минимальная высота отра­ жения тем меньше, чем больше угол входа Ѳо и тем больше, чем больше величина коэффициента подъемной силы сѵ. Из фор­ мулы (13.6) видно также, что при увеличении нагрузки на еди­ ницу характерной площади (площади крыла G/5KP), высота отражения уменьшается.

Время рикошетирования найдем, переписав выражение (13.5) для левого участка отражения (Ѳо<Ѳ<0) в следующем виде:

dH_

T I — 2

£

V

е-ря

1/2

 

d t

ß

v t

 

(13.7)

319

где

» = w + 2 f ( j ;

е - р я 0 = const.

( 1 3 . 8 )

Введем обозначения для текущих /( те ) и начальных /(те0 ) условий отражения

 

 

у

Ш £

 

 

 

1

/ 2

 

 

 

 

 

 

*

-

2

f

 

Ѵ ж — ѳ

 

 

 

/ т -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

/ 2

/ І

 

+

ѳ

 

/

т

е

 

+

 

 

 

 

 

 

 

(л)Ѵ№

 

 

 

 

(13.9)

 

 

/

т

е т— е

2

1 / 2

 

 

 

/

(ЭЛ)о=

1 / 2

УV т ет е

Ѳ+

ѳ 0

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

0

 

 

 

/

т

е+

т е

2g

 

 

 

 

 

 

5Ш= —те

 

 

2

у

(0—л М

 

 

 

 

 

 

 

 

 

н к г

 

 

 

 

 

 

1/2

 

 

 

 

 

 

 

 

 

 

 

После интегрирования и определения постоянной интегриро­ вания при условии, что Н = Н 0 и /( те )=/(те0 ), выражение

(13.7) дает

t-.

/те

■ In

/ т )

(13. 10)

/ т 0)

Заменив знак «минус» на «плюс» в выражении (13.7) приме­ нительно к правому участку отражения (Ѳо>Ѳ>0), легко убе­ диться в справедливости выражения (13. 10) и на этом участке полета.

Полное время полета т при отражении получим, положив в формуле (13. 10) Ѳ = —Ѳо. Обозначив индексом «к» параметры конца отражения и принимая во внимание, что

получим

Шк

 

/Ѵк \2

 

 

%

U o / ’

 

 

 

v_

к

 

+

(13.11)

 

/>Kтl

е

 

ß

т / a t —

ѳ 0

Легко видеть, что время отражения тем больше, чем больше угол входа Ѳо и чем меньше конечная скорость торможения Ѵк.

Текущую высоту на участке рикошетирования, выраженную в функции времени полета # = #(/), легко определить по фор­ муле (13.10), решив ее относительно параметра Н. Тогда, с уче­ том обозначений (13.9), получим

320

 

 

l - l

1

cm

\ 2"

 

 

f m

 

н=-

■ln m

1+

/(3Ro)e <р/ эк

(13. 12)

Характер зависимости

H —H(t),

ß I /

по формуле

 

 

построенной

(13.12), представлен на рис. 181. Из графика видно, что траек­

тория отражения

симметрична

относительно минимальной вы­

соты

Н = Н т

 

как мы уже видели, определяется по

 

щ. Последняя,

формуле (13.6).

траектории

 

в функции

времени

отражения

Угол

наклона

 

найдем из выражения (13.10),

решив его

относительно пара­

метра Ѳ.

При этом найдем

'

-

/<Ю1

.

(13.13)

 

 

 

:V S

 

 

 

 

1

+

/ (Silo)

 

 

Зависимость Ѳ = Ѳ(() представлена на рис. 182. Как видно из графика, в процессе полета траекторные углы Ѳ(^) распределены симметрично относительно половинного времени отражения.

Время, при котором угол наклона траектории становится рав­ ным нулю (половинное отражение), найдем по формуле (13. 13), положив в ней Ѳ = 0. Тогда получим

1__1п У ЯК -

В0

(13.14)

ß / З К / Э К +

Ѳ0 ’

Этому же времени полета соответствует минимальная высота отражения. В этом легко убедиться по формуле (13. 12), прирав­ няв нулю производную dH/dt.

Имея выражения (13.11) и (13.14), найдем

t0=0

1

(13.15)

тѴп

,+ ГГ

Поскольку Ко~Кк, то (ѳ+о/т~1/2. Таким образом время до­ стижения минимальной высоты отражения составляет примерно половину общего времени полета в рикошетной фазе.

Найдем далее распределение скоростей на участке рикошета. Тогда сделав замену dt = dH/VQ в первом уравнении системы

(13.2), получим

(13.16)

V 04)2

321

 

Перепишем затем формулу

(13.5)

в следующем виде

 

 

Обозначим

 

 

 

2g

e-ßtf

1/2

(1 3 . 17)

 

 

 

 

 

 

 

 

 

Н К ) 2

 

 

 

 

х 2 = Ш

2g

е-ря

 

 

и

d H -

g

 

dx.

ß(Vyf

 

 

 

 

запишем

 

(13. 18)

 

 

 

 

 

 

-ß#

 

 

 

(

К

Г

 

(13.17)

и

 

 

 

Подставляя

выражения

 

 

(13.18) в уравнение (13. 16),

получим

 

 

 

 

d V

d x

 

 

 

 

 

 

 

V

К

 

 

 

 

 

 

 

 

 

 

 

 

О 10 20 30 00 t ,c

О

10 20 30

00 t,C

Рис. 181. Зависимость H = H ( t )

Рис. 182.

Зависимость

Ѳ = Ѳ(0

Интегрирование этого выражения с учетом начальных усло­ вий Ѵ = Ѵ 0 и Н = Н 0 дает

— = ехр

(13. 19)

Из

 

Скорость полета в конце отражения (ее относительное изме­ нение при Ѳ = —Ѳо) будет равна

-exp 2 -

(13.20)

Ѵо

Изменение скорости в функции времени полета при отраже­ нии показано на рис. 183. Как видно из графика, рикошетная фаза осуществляется с потерей кинетической энергии, причем величина этой потери возрастает с увеличением угла входа

322

и уменьшением аэродинамического качества. Ввиду малых углов, входа, эти потери невелики и, как мы уже говорили ранее, со­ ставляют порядка 8%.

Поперечная перегрузка, равная отношению подъемной силы к силе тяжести, запишется так:

Пи

I у *

'

п

 

у J

 

или с учетом (13. 17) и (13. 19)

»а

1/§(да-Ѳ2)1ехр1

см

=*|0

L2^o

J

 

1!

1

 

L \сх}

1

о1 _

1. -

J

(13.21)

(13.22)

Характер зависимости перегрузки от времени полета показав на рис. 184. Как видно из графика, поперечная перегрузка при

Рис. 184. Зависимость пу—

= n y (t)

Рис. 183. Изменение скорости

юпо времени полета t

ч!—

1,0

 

 

0,9

0

= j n o{d L -L

 

 

 

 

 

 

 

 

 

 

 

0,8

иО

( Ш 7 Сх ¥

 

 

 

 

 

 

 

 

 

 

 

 

О

10

2 0 3 0 9 0 t , с

 

 

 

отражении

 

О

10

20 30

90 t,C

достигает значительных величин.

Формула (13.22)

показывает, что главными параметрами,

определяющими вели­

чину

п у ,

являются скорость Ко и угол входа Ѳо-.

 

На минимальной высоте отражения,

где

Ѳ = 0,

перегрузка

равна

 

 

 

 

 

 

 

 

 

 

о) = 2 І 1/оЖе

 

 

(13.23)

Как видно из графика (см. рис. 184), перегрузка, определяе­ мая выражением (13.23), близка к максимальной п у т ах. Значе-

323

ние последней легко определить по формуле (13.22), если изве­ стен угол Ѳ = Ѳп^тахі при котором имеет место пу= путах. Для

нахождения этого угла продифференцируем по Ѳ выражение (13.22) и приравняем нулю производную (dny/dB). При этом получим

Ѳ2 +

Cjl

6 — ЭД = 0,

Сх

откуда

(13.24)

Рис. 185. Максимальная перегрузка

Рис. 186. Максимальная перегрузка в

рикошетирующего полета в зависимо-

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

дальности при входе

сти от дальности

в плотные слои

атмосферы по балли­

 

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

Формула (13.24) показывает, что поперечная перегрузка до­ стигает максимального значения несколько раньше, чем насту­ пает Ѳ=0.

Для ракетопланов, имеющих аэродинамическое качество ( — )>1, угол Опытах можно считать практически равным нулю.

При этом для определения путах с достаточной точностью можно пользоваться формулой (13.23).

Необходимо отметить, что при движении в режиме рикошетирования возможны перегрузки порядка пу = 20—30, что яв­ ляется весьма серьезным обстоятельством с точки зрения проч­ ности и веса конструкции, а также с физиологической точки зре­ ния, если ставится задача обеспечения необходимых условий для полета экипажа.

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

324

Для определения горизонтальной дальности рикошетирования найдем из последних двух уравнений системы (13.2)

і =

(13.25)

Делая подстановку при помощи выражения (13. 17) и выпол­ няя интегрирование, находим

l

=

-р ± /=1

т

- Ш І

(13. 26)

 

/ т о )

В конце отражения Ѳ = —Ѳо, Н —Н 0 и следовательно

ß

2__in

Ѵ ж +

«о

(13.27)

 

 

Ѵ ш

Ѵ ш — ѳ0

 

Условной границей, разделяющей полет ракетоплана в ква­ зипустотной среде по баллистической траектории и в режиме отражения, целесообразно принять квазистационарную высоту планирования. Поскольку в режиме квазистационарного плани­ рования имеет место dQ/dt = 0, то из второго уравнения системы (13.2) легко получить

'Ѵ_

ß//0— g

 

0,

g к

 

откуда при V —Vo

-— ln

V ’

\(—)

(13.28)

H n=

 

 

 

^KP /

 

Как видно из полученной формулы, высота условной границы начала отражения определяется главным образом параметрами

V, Су и G/S. При этом увеличение скорости входа и коэффи­ циента подъемной силы увеличивает начальную высоту отраже­ ния, а увеличение нагрузки на единицу площади крыла умень­ шает эту высоту.

Обычно рикошетирующее движение сравнивается с квазистационарным и баллистическим режимами спуска. Главными пара­

метрами сравнения являются

дальность,

температура и

пере­

грузка. На рикошетной

фазе

возникают

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

и перегрузки, чего нет

при гиперзвуковом планировании.

В то

же время дальность рикошетирующего полета превышает

сум­

марную дальность баллистических и рикошетных фаз примерно на 12%. Что касается тепла, то его интегральная величина ока­ зывается наименьшей при баллистическом спуске [5]. Таким образом использование рикошетирующих траекторий для воз­ вращения ступеней, обладающих подъемной силой и охлаждае­ мых излучением, представляется менее перспективным, чем

325

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

13.2. ПРИБЛИЖЕННЫЕ ЗАВИСИМОСТИ ДЛЯ ПОЛОГОЙ БАЛЛИСТИЧЕСКОЙ ТРАЕКТОРИИ

Рассмотрим теперь пологую баллистическую траекторию, когда возвращаемая ступень ракетной системы многократного применения, выйдя из рикошетной фазы, попадает в квазипу­ стотную среду и затем снова входит в атмосферу. Вследствие малых углов отражения, для пологой баллистической траекто­ рии справедливы допущения co sB ^ l, эіпѲ^Ѳ и Ѵ = const. Диф­ ференциальные уравнения движения для таких траекторий при­ нимают вид

 

 

 

 

 

 

d V

= 0,

Ѵ_\2

 

(13.29)

 

 

 

 

V

±

dt

 

 

 

 

 

■ g

 

(13.30)

 

 

 

 

dt

 

 

 

 

 

 

 

=ѵв

Р к р

 

 

(13.31)

 

 

 

 

 

 

dt

 

(13.30)

и (13.31) легко

 

Делением почленно уравнений

получить

 

ѲаГ0 =

Ѵ2

['

V

2

d H ,

 

 

откуда

 

JL

V,кр

 

 

 

 

 

 

 

 

 

 

 

где

L

— постоянная интегрирования.

 

 

 

условия

 

Для

восходящего

участка

траектории начальные

равны Ѳ = Ѳ0 и

Н = Н 0.

При этом

 

 

 

 

(13.32)

 

 

 

 

 

где Н 0— начальная высота баллистического полета.

 

На максимальной высоте полета

(в апогее)

имеет место Ѳ = 0.

При этом формула (13.32)

дает

 

 

 

 

где Н щах

 

Н .

2g

(Wo)2

 

 

(13.33)

 

 

 

 

 

 

 

высота апогея.

 

 

 

 

 

 

326

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