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

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

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

или

 

 

V II —

 

 

(9. 58)

где индексы

I и II

к режимам

 

относятся

планирования при

Vi<VgR

и

ѵп> Y

соответственно.

Ѵі

 

 

(9.57)

и

 

 

Формулы

 

(9.58)

идентичны.

Они устанавливают

связь (соотношение)

между

скоростями

 

и Ѵп,при которых

высоты планирования остаются одинаковыми, несмотря на не­ равенство на этих высотах скоростей полета Ѵ іФ Ѵц.

Рис. 127. Параметры квазистационарного битраекторного планирования

Существует единственное решение, при котором, в соответст-

ствии с формулами (9.57) и (9.58), имеет место Ѵ 1= Ѵ п Y g R . Это означает, что при входе в атмосферу со скоростью, превы­ шающей первую космическую скорость, максимальная высота будет тогда, когда скорость планирования упадет до величины, равной первой комической скорости. Эта точка в свою очередь является общей для начальной высоты планирования при поло­ жительной подъемной силе.

Если в формуле (9.58) положить равным нулю знаменатель,

то получим

= — ^ 0 ,7 1 .

(9.59)

 

Екр

/2

 

8*

227

 

Это означает, что при скорости входа, равной бесконечности (1/0ІІ= оо), существует предельно малая высота входа (началь­ ного планирования) Я 0 = Н 0тіп, при которой обеспечивается за­ хват ЛА атмосферой планеты и требуемый режим планирования с отрицательной подъемной силой. Этой минимальной высоте входа соответствует скорость на нисходящем участке планиро­ вания (при положительной подъемной силе), составляющая по

величине 71% от первой космической скорости.

 

 

Общий характер кривых

Н = Н(Ѵ)

для квазистационарного

битраекторного планирования,

построенных по формулам (9.44)

и (9.54) для различных значений удельной .нагрузки

G/S,

пред­

ставлен на рис. 127. Как видно из графика, влияние

удельной

нагрузки

G/S

на высоту планирования одинаков для двух уча­

 

 

 

 

 

 

 

стков при V gR. Очевидно, что поворот вектора аэродина­ мической подъемной силы Y (V) при переходе с восходящего на нисходящий участок планирования происходит в точкб, где име­

ет место V = У gR.

 

Пример 9. 1. Определить

дальность

активно-инерционного планирования

если

дано: скорость истечения 117=4000

м/с, весовая

отдача р,к=0,27, циклов

я =1

= 4 , аэродинамическое

качество

(су/сх) = 4,5,

стартовая перегрузка

«0=1,5, коэффициент гравитационных потерь о(Ѳи)= 1. Силами сопротивле­ ния воздуха на участке выведения пренебречь.

Решение. Весовая

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

риодам определяются

выражениями

 

 

 

=

V K = W Mn Ц

а (Ѳк) — J ,

L ИЛ

RißyjCx) ,

\

«о /

------------ ln

 

 

Вычисления по этим формулам сводим в таблицу.

Количество периодов п

1

2

3

4

Весовая от.~ача пк

0,27

0,52

0,645

0,72

Число Ц

3,7

1,925

1,55

1,39

Величина In Ц

1,31

0,65

0,44

0,33

Величина От/п0

0,49

0,32

0,24

0,19

Скорость выведения Ѵк, м/с

3280

1320

800

520

Величина (17к/Ккр)2

0,16

0,0174

0,01

0,0042

Величина 1— (Ук/Икр)2

0,84

0,983

0,99

0,9958

Величина In ----------- ------------

0,18

0,03

0,01

0,004

1 - (Ѵк/Ѵкр)2

 

 

 

 

Дальность одного периода Lt, км

2600

430

144

58

Суммарная дальность, L^—nLi, км

2600

860

430

230

228

Как видно из таблицы, дальность одного периода больше дальности че­

тырех периодов в 2600/230«11 раз.

 

 

Пример 9. 2. Составить

аналитическую

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

разность пути, отнесенного

к величине

R

! °у \

I — I, при движении ракето­

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

яйцах

при Ц = 6 ,

 

Ѵ кр =

^кр =

 

2,3 и а (Ѳк) = 0.

 

 

 

 

Решение. Дальность гиперзвукового планирования определяется по фор­

муле

 

 

 

 

 

 

LII.I

=

 

 

 

 

 

 

 

 

где Ѵк — конечная

 

 

 

 

 

 

 

 

 

 

 

скорость выведения (или все

равно, что

начальная ско­

рость

планирования).

 

 

в режимах одного и п

 

 

 

Разность

пути,

 

пройденного

периодов,

составит

Ln=1Lп—• —

R_ Су_

 

 

 

 

 

 

 

 

2

In

 

 

 

 

 

 

 

 

 

 

 

с х

 

 

 

 

 

 

 

кр

где Кк

и Ѵк„ —

скорости в конце активного участка при одном и п периодах

 

2

 

 

 

истечения W.

 

 

 

 

 

полета,

отнесенные

 

к скорости

 

Ѵя„ —1п ln Ц \2

Поскольку

о

 

(Ѳк)

= 0 имеют

место

ln Ц \ 2

и

 

AL

R_

1

 

 

 

— ln

1—

^кр J .

+

n ln

1—

nVJ

J’

2

и затем

n

AL = ln

ln Ц \2

кр

Подстановка в полученную формулу значений параметров, указанных в условиях примера, дает следующую разность относительной дальности:

 

Д L ----- In

[1 — (ln 6/4-2,З Я

 

 

 

0,77.

 

 

 

1 — (ln 6/2,3)2

 

 

При

аэродинамическом

качестве, взятом,

например

Су_ = 6 ,

 

 

 

 

с х

это означало бы, что превышение дальности одного

периода

планирования

по сравнению с суммарной дальностью шести периодов планирования состав­ ляет величину, равную

 

Д 7 —

( —

) = 0 ,7 7 (#з/2)4«

1000 км.

 

2

\ с х

J

 

 

Пример

9.3. Найти соотношение

между начальной скоростью входа

У о п > / J R

и конечной скоростью

планирования

Ѵ 1 </■ /"g R , при котором

229

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

участках.

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

лучим

 

 

 

 

 

 

 

 

^ оп1 Л = Ѵ^р.

 

 

 

 

 

 

 

(9.60)

Таким

образом,

скорость

в конце

планирования___V ,< ^ -\ fg R

указывает

на вполне

определенную скорость входа Von ]>

gRi

если

время

 

движе­

ния по восходящему и нисходящему участкам должно быть одинаково.

Если скорость входа равна второй космической скорости

Von = / 2 g R , то

формула

(9.60)

дает

V \ — -/gR /2.

 

Любопытно

заметить,

что

 

в

этом

частном случае характер

изменения

скорости

по

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

планирования

следует

закону

геометрической

прогрессии

со

знаменателем

1 Іу '2 . При

этом:

Ѵои =

/ 2 g R (при

входе), Ѵі = Ѵ\і =

у1g R

точке

восходящего и. начала нисходящего участков планирования)

и

Ѵц =

у

g R /2

(в конце

планирования).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СП И СО К ЛИТЕРАТУРЫ

 

 

 

 

 

 

 

1.

Д у б о ш и н

Г.

Н.

Небесная механика. М.,

Физматгиз,

1963,

586с.

2.

Исследования

по

динамике

полета.

Под

ред. д-ра

техн.

наук

проф.

И. В. Остославского. [Сборник статей].

М.,

«Машиностроение»,

1965,

вып. 1,

355 с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

М а р т и н

Д ж .

Вход

в атмосферу. Введение в теорию

и

практику.

Под ред. канд. техн. наук Н. А. Анфимова. М., «Мир», 11969,

320 с.

 

 

 

4.

О р л о в

Б.

В.

и др.

Основы проектирования ракетно-прямоточных дви­

гателей для беспилотных летательных аппаратов. Под ред. д-ра техн. наук

проф. Б. В. Орлова. М.,

«Машиностроение»,

1967,

424 с.

 

 

 

 

 

 

5.

Р а б и н о в и ч

Б. И.

Вариационные

режимы

полета

крылатых

лета­

тельных аппаратов. М.,

«Машиностроение»,

1966,

182 с.

 

 

 

 

 

 

6.

Р у п п е

Г.

Введение в астронавтику,

т. 1. М.,

«Наука», 1970, 612с.

 

7.

Ш к а д о в Л. М.

и др. Механика оптимального пространственного дви­

жения летательных аппаратов в атмосфере. М., «Машиностроение»,

1972,

240 с.

 

8.

Ш т е р н

Т.

Е.

Введение в небесную механику.

М., «Мир»,

1964,

242 с.

9.

Э л ь я с б е р г П.

Е. Введение в теорию полета

искусственных

спутни­

ков Земли. М., «Наука», 1965,

540 с.

 

 

 

 

 

 

 

 

 

 

10.

Э р и к е

К.

Космический

полет, т. 1.

М.,

Физматгиз,

1963,

586 с.

 

11.

Я р о ш е в с к и й

В. А. Приближенный

расчет

траекторий

входа

в атмосферу. — «Космические

исследования»,

 

1964, т.

2,

вып.

4, с.

507—531,

вып. 5, с. 579—696.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12.

B e c k e r

J.

V.

 

Re-Entry

From Space,

Scientific

American,

January,

1961,

p. 10—59.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

July,

13.

C a z 1e у C.

Atmospheric

Entry. The

RAHD

Corp.,

Rep.,

P-2052,

1960.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г л а в а X

ОПТИМАЛЬНЫЕ ПАРАМЕТРЫ ОДНОСТУПЕНЧАТЫХ РАКЕТ

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

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

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

Действительно, без учета весовых факторов и, в частности, влияния перегрузки на весовую отдачу, второй закон Циолков­ ского (2.44) запишется так:

где полагаем, что величина параметра рк не зависит от пере­ грузки п0, т. е.

Рк= Рі((«о) = const.

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

Этот чисто академический вывод не имеет практического смысла. При увеличении тяги увеличивается вес реальной конст­ рукции, что приводит к увеличению коэффициента цк = Ц к ( « о ) -

Таким образом, варьирование перегрузкой

п0

вызывает измене­

 

ние весовой отдачи по конечному весу, влияние которой на ско­ рость Ѵк требует отыскания наиболее выгодного (оптимального) значения перегрузки п0. Зависимость |Хк = Цк(яо) учитывает вли­

231

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

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

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

10.1. ОСНО ВНЫ Е ПАРАМЕТРЫ ОДНОСТУПЕН ЧАТЫ Х РАКЕТ

К основным параметрам одноступенчатых ракет относятся:

— весовая отдача по топливу G T;

весовая отдача по конечному весу цк;

число Циолковского Ц;

коэффициент грузоподъемности е;

стартовая перегрузка я0;

—'коэффициент расхода активной массы (топлива) ß;

коэффициент эффективности конструкции цЯфф;

скорость истечения активной массы (топлива) IF;

удельная тяга Р ул.

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

Функциональная взаимосвязь между основными параметра­ ми имеет следующий вид (см. гл. II, V I):

Ц — 1/Рк> ^ —:£Руд> Рк —

( 10. 1)

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

232

10.2. СТАРТОВАЯ ПЕРЕГРУЗКА В ПУСТОТЕ ПРИ ЛИНЕЙНОМ ЗАКОНЕ ИЗМЕНЕНИЯ МАССЫ

Оптимальность стартовой перегрузки для случая полета ра­ кеты в пустоте определяется двумя противоположно действую­ щими факторами. С одной стороны, желание уменьшить потери на преодоление силы тяжести -указывает на необходимость уве­ личения перегрузки п0; с другой ■— увеличение перегрузки п0 связано с понижением весовой отдачи GT, что вызывает умень­ шение конечной скорости Ѵк. Последнее обстоятельство связано с действием весовых факторов, учет которых произведем, вос­ пользовавшись основной весовой зависимостью вида

Параметр [лэл не связан функционально с перегрузкой и яв­

ляется постоянной величиной. Коэффициент k при параметре «о > как было показано ранее (см.

разд. 6. 1.), учитывает изменение ко­ нечного веса ракеты в связи с из­ менением тяги двигателя (или все равно, что перегрузки п0п). Напом­ ним также, что на величину коэф­ фициента k влияют и те элементы конструкции, вес которых зависит «от тяги. Сюда можно отнести арма­ туру двигательного отсека, несущие корпусы, частично баки и крепеж­ ные детали.

Рис. 123. Влияние перегрузки nJJ на конеч­

ную скорость полета в пустоте при линей­ ном законе изменения массы

Воспользовавшись вторым законом Циолковского, запишем

Й к = In -Нк---- а(ѲJ

nl

,

(10.2)

где V„ — —- — конечная скорость, отнесенная

к скорости ис-

к

 

 

 

течения в пустоте.

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

Параметр ст(Ѳк) в формуле (10.2)

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

висимость 10.8).

Подстановка «выражения (6.9) в уравнение (10.2) дает

(10.3)

1/к = 1п

k n 0

^ - [ І - ^ э л + Ю ]

Иэл +

 

«0

 

233

Равенство (10.3)

устанавливает

зависимость

безразмерной

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

Ѵк

от параметров, характеризуемых

учетом

динамических и весовых факторов.

 

 

 

 

 

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

V * — Ѵ к

(«о)

для различных

значе­

ний параметра а(Ѳк), построенной

по уравнению

(10.3), пред­

ставлен на рис. 128. Как видно

из графика, при значениях

о(Ѳк)> 0 , функция

\7к = 1/^(/го)

имеет

максимумы, которым

соответствуют наивыгоднейшие (оптимальные) значения

пере­

грузок

По.

 

 

 

 

 

 

 

 

Из графика видно также, что эти перегрузки тем меньше, чем

меньше гравитационные потери.

Если последние

отсутствуют,

т. е. а(Ѳк)= 0 , то функция

Ѵ к =

Ѵ к(по)

не имеет экстремума

 

 

 

и теоретически максимальная скорость достигается при значени­

ях «о =0. Этот вывод имеет смысл для космических

летатель­

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

быть пре­

дельно малы (ng<Cl).

имеет местоѴ ке== const, то оче­

Поскольку в выражении (6.9)

видно, что отыскание экстремума

функции

Ѵ к (п'о)

воз­

 

 

можно при постоянном весе полезного груза As = const и посто­

янном стартовом весе ракеты Go = G0i = const. При этом,

прирав­

нивая нулю производную

d V Jd n ",

получим

 

 

 

 

kn” =

 

 

k (л”)2

 

з (flK) (1 - і»вл)

 

0,

(10.4)

- з (Ѳк) рэл (1 - ң,л) -

 

откуда оптимальная перегрузка

по = п"пт

будет равна

(10.5)

^ = 3(ѲК)

 

 

 

 

 

Выражение (10.5) определяет оптимальное значение старто­ вой перегрузки /г"пт в функции параметров, задаваемых при

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

Выражение (10.5) может быть представлено и в другом виде. Для этого разрешим уравнение (10.4) относительно параметра По после подстановки в него выражения рэл = —kn™.

Тогда получим

 

пU

 

6

а (0К)

р,

( 10. 6)

 

 

° (2к)

 

 

 

__

 

 

— Оа

 

 

о —

 

 

 

где GT= ( 1—Цк)— весовая отдача по топливу.

 

 

График оптимальных перегрузок,

построенный по формуле

(10.6)

для различных значений параметра а = а(Ѳк),

представ­

лен на рис. 129. Характерным для кривых на графике

является

наличие максимума функции

пЦ =

Пд(р-к)

 

для оптимальных

значений параметра

^о'—^оопт-

Этот максимум указывает

на

необходимость

уменьшения

оптимальных

перегрузок

234

яо — ^оопт как на малых, так и на больших значениях пара­

метра цк. Физическое объяснение такого явления заключается в следующем.

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

чае были бы велики ввиду увеличе­

з,5

ния параметра

G T,

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

К

и увеличения времени полета т.

о

Ухудшение при этом весовой отда­

jr

чи [хк (увеличение этого параметра)

не оказывает доминирующего влия­

20

ния на конечную скорость, посколь-

ку относительный вес двигательной

25

установки в общем

весе конструк-

ции (конечном весе ракеты)

неве-

лик. При движении от средних к ма-

т

лым значениям параметра рк умень-

_______________________________________________

^

Рис. 129. Влияние на

оптимальную

пере-

грузку гаg параметров

рк и сг(Ѳк) при ли­

 

нейном законе изменения массы

 

q

шение оптимальной перегрузки связано с возрастанием относи­

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

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

довательно, и относительный вес двигательной установки) дол­ жна быть уменьшена.

При больших значениях параметра рк относительный вес двигательной установки невелик, однако при этом слишком ма­ ла скорость по Циолковскому, которую можно увеличить только вследствие уменьшения параметра рк, а следовательно, благода­ ря уменьшению перегрузки п0. График показывает также, что оптимальные перегрузки пропорциональны гравитационным по­ терям, характеризуемым параметром <т(0к).

Из формулы

(10.6) видно, что при Цк=Г оптимальная пере­

грузка

.

в пустоте равна

пЦ

= о (Ѳ к)-

Если о(Ѳк) = 1 и |лк= Г то

/г^пт= 1

Это

случай

стоянки на

стартовом столе гипотети­

ческой ракеты, у которой отсутствует топливо. При цк = 0 и всех

235

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

Крайние области по параметру рк> где äJ < 1, не имеют

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

ва от стартового стола всегда необходима

перегрузка

больше

единицы П оЖ

Однако график не

нуждается

в

перестройке,,

поскольку эти области применимы

к космическим

кораблям и

старшим ступеням составных ракет,

у которых

имеет

 

место

а(Ѳ

к )

< 1 и

по

1.

 

 

 

 

 

 

 

 

 

Выражение для параметра рк, при котором имеет место мак­

симум функции

 

^опт^^опт^к)

нетрудно

получить по вы­

ражению (10.6), приравняв нулю производную

dn^/d^.

При

этом

 

4 — о (Ѳк)

k

4 — а (Ѳк ) k

4

а (%K) k

jl/2

 

(10.7)

 

 

 

 

 

 

1

 

 

 

Рис. 130. Зависимость парамет­

Рис. 131. Зависимость дальности, вы­

ра

гравитационных потерь

соты, времени и угла Ѳк от конечно!)

а(Ѳк) от конечного угла Ѳк для

скорости ракеты

 

одноступенчатых ракет

 

Если пренебречь малыми значениями произведений а(Ѳк)& (обычно меньших 0,023 для Ж РД ), т. е. положив а(Ѳк)&= 0, То получим Цк^ 1/2, Это означает, что наибольшие значения опти­ мальных стартовых перегрузок (По)опт имеют место для ракет с половинной весовой отдачей.

Значения параметра <т(Ѳк), рассчитанные для оптимальных траекторий одноступенчатых ракет, приведены на рис. 130 Гра-

236

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