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

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

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

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

 

X

 

 

 

а О У — ~ ^

sin Bät,

(10 .8)

где Ѳ

о

к местному

горизонту; т —

угол наклона траектории

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

Рис. 132. Зависимость основных параметров движения ра­ кеты от дальности при полете по баллистической траекто­

рии (без учета дальности активного участка)

Таблица 10.1

Параметры баллистической траектории

2?

К

L, км

Кк, км/с

т

Я , км

1667

3772

10'47,4"

495

15°

41015'

30°

37°30'

3333

5033

16' 0 9,0"

827

45°

33°45'

5000

5839

2 0 '1 8 ,6 "

1093

60°

30°00'

 

6408

2 5 '2 9 ,8 "

1286

75°

26°15'

8333

6828

2 9 '2 7 ,5 "

1402

°

15°00'

10000

7144

32' 5 9,7"

1441

90°

22°30'

6666

105°

18°45'

11667

7383

3 6 '1 1 ,8 "

1402

120

7°30'

13333

7562

38' 37,4"

1286

135°

11°15'

15000

7692

4 0 '3 8 ,2 "

1093

150°

 

16667

7780

4 2 '0 4 ,6 "

827

165°

3°45'

18333

7832

4 2 '5 6 ,1 "

495

180°

°00'

20000

7849

4 3 '1 3 ,1 "

100

 

0

 

 

 

П р и м е ч а н и е . — центральный угол, Ѳк—конечный угол (в момент вык­ лючения двигателя), Ѵ к—конечная скорость, L, Т, Н —соответственно даль­

ность, время и высота баллистического полета.

237

Другие траекторные данные, связанные с конечным углом Ѳк различными параметрическими зависимостями, приведены в таб­

лице 10.1 и на рис. 131 и 132.

 

Ѳ к =

( Ѳ к ) о пт [4 ].

Таблица составлена для оптимальных углов

Эти данные необходимы также

для отыскания

значений

<т(Ѳк) при вычислениях оптимальных

перегрузок

и других па­

раметров ракет, в том числе и многоступенчатых (см. гл. X II).

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

Изучение оптимальности параметра по при показательном законе изменения массы произведем как и в предыдущем пара­ графе, воспользовавшись весовой зависимостью (6.9). После подстановки этой зависимости в формулу (2.47) получим

Ѵ „ =

( 1 0 . 9 )

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

связывает

б(в„)

между собой основные пара­

 

метры движения и параметры,

 

характеризующие влияние на

 

конечную скорость

весовых

 

факторов.

 

 

о

 

2

 

 

8 п£

 

 

 

 

 

Рис.

133.

Влияние

перегрузки

Рис.

134.

Зависимость оптимальных

 

4

6

 

п0п на конечную скорость по­

значений стартовой перегрузки от ве­

лета в пустоте при показатель­

совой

отдачи

для случая

полета в

ном

законе

изменения массы

пустоте

при

показательном

законе

 

 

 

 

 

 

 

 

изменения массы

 

238

На рис. 133 представлены кривые Ѵ к = Ѵ к(по),

построен­

ные по уравнению (10.9) для различных значений

параметра

гравитационных потерь а(9к). Как и в случае линейного закона

изменения массы, при

показательном

законе

функция

Ѵ к — Ѵ к

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

влиянием

весовых

(яо)имеет максимум,

факторов. Уменьшение гравитационных потерь также уменьша­ ет оптимальную перегрузку, которая, как и при линейном зако­

 

 

 

 

 

Ѵк =

Ѵ к (гц)

 

 

не изменения массы, становится равной нулю при а(0к) =0.

 

Экстремальное значение функции

 

 

dVJdtio.

 

 

по уравне­

нию (10.9)

найдем, приравняв нулю производную

( .

)

При этом,

пп

 

 

3 (9к)

разд.

 

поступая как и ранее (см.

10.2), получим 10 10

 

0

3 (ѲЭ

I

 

 

 

 

 

 

~1“

+

 

п^—п^

((\), построенные

по

Кривые оптимальных значений

 

формуле (10.10) для различных значений параметра о(Ѳк), при­

ведены на рис. 134.

для=случая показательного расхода

Как видно из графика,

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

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

ra(J ra{J (рк)

имеет максимум. Как

и ранее (см. разд. 10.2), этот максимум характеризует влияние на оптимальные перегрузки весовых факторов. По сравнению с линейным законом, максимумы кривых rajj—/г" (|хко (Ѳк)

при показательном законе сдвинуты в сторону малых значений коэффициента цк. Это свидетельствует о несколько меньшем вли­

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

яоІ;==/гопт

и объясняется

тем, что конструкция ракеты работает не на максимальные (ко­

нечные) перегрузки, а на минимальные,

равные начальным осе­

вым перегрузкам rag.

 

по сравне­

Увеличение абсолютных значений параметра

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

Значение коэффициента цк, при котором имеет место макси­

мум функции rag = rag ((А«) Для оптимальных значений rag,

как

и ранее,

найдено по выражению

(10.10), приравняв нулю произ­

водную

dn" d\^K.

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

цк= 1/е~0,368. Таким

обра­

 

езомт ( ,/ г онеобходимость уменьшения гравитационных потерь сдвига-

п т ) т а х

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

нению с линейным законом на величину Ацк~0,132.

По формуле (10.10) имеем также rag— 0 при цк = 0 и rag =

= а(Ѳк) при рк=1.

239

Как и при линейном законе изменения массы, график (см. рис. 134) не нуждается в коррекции крайних областей по пара­ метру рк, поскольку такое построение функции п%= п"[ІАК)

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

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

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

ным законом имеют место на всем

участке 0 < р к < 1 ,

причем

расхождение в значениях

Пд

тем больше, чем меньше величи­

на параметра рк. Указанное явление

имеет место для

любых

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

Пользуясь формулами (2.36) и (2.38), найдем относитель­ ное увеличение времени полета при показательном законе изме­ нения массы по сравнению с линейным законом. Оно составит

 

T

= т =

- I s(«’°г!)S

--ІпЦ= - -

,

- - -

ѵ -

°

^( 1 и н0 . 1 1

 

 

Тлин

( я д ) гі 0к (1

(Чс)

 

 

 

 

где индексы «пок» и «лин» означают параметры, относящиеся к

показательному и линейному законам изменения массы.

выра­

Для случая

одинаковых

перегрузок

(л.0)лин =

(я[5)

 

жение (10.11) принимает вид

ИД

 

 

 

 

 

 

( 10. 12)

 

 

 

1— Нк

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

значений

Кривые т= т(р к) для одинаковых

перегрузок я о ,

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

(10.11)

и

(10.12),

приведены на рис. 136. График показывает, что увеличение вре­

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

 

законе

изменения

массы особенно заметно при малых значениях параметра рк. Из

графика видно

также,

что оптимальные

перегрузки

заметно

уменьшают это время, что и способствует уменьшению гравита­

ционных потерь.

Для большей иллюстрации сказанного на рис.

240

137 представлены кривые Гк = Гк(Рк)> построенные по формулам

(10.2) и (10.9) при одинаковых значениях перегрузок п о . Гра­ фик показывает, что ракеты с расходом массы, подчиняющимся показательному закону, значительно проигрывают в конечной скорости ракетам с расходом массы, подчиняющимся линейно­ му закону, если значения перегрузок в одном и другом случаях одинаковы. На другом рис. 138 приведены аналогичные кривые Й к=Гк(|ік), но построенные с учетом оптимальных перегрузок. Как видно из графика, оптимальные перегрузки значительно

Рис. 135. Сравнение оптимальных перегрузок п'д для линейного и по­

казательного законов изменения массы:

/ —линейный закон ; 2—показательный закон

Ѵк

16

 

\

f

і

 

1

Ь2

 

V

і

 

 

 

 

 

 

 

 

!

 

 

\

і

 

 

 

_______ А ___А __I

 

 

0,8

X

1 \

 

!

 

ОЛ

 

І

>А ч

і

 

0

0 , 2

0 ,1 t

0 , 6

0 , 8 р к

Рис. 137. Сравнение безразмерных

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

Пд = const:

I —линейный закон; 2 —показательный закон

Рис. 136. Зависимость относительного увеличения времени полета ракеты на активном участке траектории для по­ казательного закона изменения массы по сравнению с линейным законом:

П

1—при одинаковых л о; 2 —при оптималь-

п

0 0,2 0 ,6 0,6 0,8 JU

Рис. 138. Сравнение_безразмерных ко­

нечных

скоростей

Ѵк для

линейного

и показательного

законов

изменения

массы

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

по ~ (п0^о.іт'■

1—линейный закон; 2 —показательный закон

241

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

10.5. СТАРТОВАЯ П ЕРЕГРУЗК А

СУЧЕТОМ ВЛИЯНИЯ АТМ ОСФЕРЫ

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

Если в пустотных условиях оптимальность стартовой пере­ грузки определяется двумя противоположно действующими фак­ торами (гравитационными и весовыми), то в сопротивляющейся среде к этим факторам добавляется еще один фактор, учитыва­ ющий потери скорости от аэродинамических сил лобового сопро­ тивления АѴх - Величина этих потерь может быть оценена изу­ чением кривых АѴХ = АѴХ (п0, Ѳк), построенных интегрированием выражения вида

 

 

 

 

 

Т

 

 

 

 

 

 

 

(10.13)

 

 

 

 

ді/ — [ JiM L d t,

 

 

 

 

где

X(t)

 

 

 

о

 

 

 

 

 

M (t)

—теку­

 

— текущая сила лобового сопротивления;

 

 

щая масса ракеты.

Ѳк) показывают,

что потери скорости

 

Кривые

АѴх = АѴх (п0,

пропорциональны стартовой перегрузке

п0

и обратно пропорцио­

нальны углу Ѳк-

Такой характер этих

кривых объясняется тем,

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

 

скорость

полета, а

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

X(t)

на участке атмосферы; при

 

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

тери скорости

АѴХ .

 

 

 

 

Ѵѵ = Ѵк(пй)

 

 

Изучение экстремума

функции

произведем при­

менительно к линейному закону изменения массы.

 

Представляя полет в атмосфере с пустотными характеристи­

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

Циолков­

ского, запишем

Wn

1п

± -

о (0J

gx

-

Ѵ й -

Ѵ х ,

(10.14)

 

1/к =

 

 

 

Д

Д

где в правой части первый член определяет идеальную скорость в пустоте по Циолковскому; второй член — потери в скорости на

242

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

деляемых увеличение тяги от стартовой Р 0 до пустотной Ро ; четвертый член — потери в скорости на аэродинамические силы

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

 

уравнений

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

движения, можно путем аппроксимации получить

(10.15)

ДНп = 0 ;-1 ) ( а Г о + ѵ/го);

Ь.Ѵх = 'Ц(п0-\-Ь)-\-с,

(10.16)

где а, V, г), Ь, с — коэффициенты

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

Зависимости (10.15) и (10.16)

с удовлетворительной для па­

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

Подстановка выражений (10.15) и (10.16)

в равенство (10.14)

и затем его преобразование дает

 

 

 

 

 

 

 

 

1/.

Ѵк

 

 

 

g (Вк)

- t o

- £

h +

v

( e - i ) ] -

(10.17)

Wn

 

 

 

 

 

 

 

t\b + c

 

 

 

 

 

 

П0

 

 

 

W q

 

 

 

 

 

 

 

 

 

Wo

m {l— 1).

имеющие

 

Коэффициенты аппроксимации b,

с и m,

место в

уравнении

(10.17), из дальнейшего рассмотрения выпадают.

Учет весовых факторов произведем, как и ранее, по зависи­

мости

(6.9). Для решения уравнения (10.17)

на экстремум необ­

ходимо произвести дифференцирование

по параметру

п0

и при­

равнять нулю производную. Выполнив эти действия, получим

 

пп

0,5а (Ѳк)

 

+

0,5а (Ѳк)

 

 

 

 

1/2

 

(10. 18)

 

k

k

 

°

(вк) <?т

 

 

 

 

 

 

 

 

 

 

 

 

■По

Ж

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

Ңс

 

 

 

 

 

 

,

е*

 

 

где

 

 

 

 

 

 

 

 

 

|Чс

 

 

= W о [Л + ѵСе— 1)].

 

 

 

 

 

 

 

 

 

Безразмерный параметр т]о учитывает влияние атмосферы на величину оптимальной перегрузки п0. Характер кривых цо = = тіо (Wo, 0К) представлен на рис. 139. Как видно из графика, ве­ личина цо в значительной мере зависит от конечного угла 0К) что в свою очередь является определяющим и для перегрузки п0. Из графика видно также, что увеличение скорости истечения пони­ жает величину г)о> что сказывается на некотором увеличении оп­ тимальной перегрузки.

243

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

различных значений а(0к), построенная

по

формуле

(10.18),

представлена на рис. 140.

Как

видно

из

графика,

характер

кривых Яо = «о[рк. <т(Ѳк)] остается таким же,

как

и для случая

полета в пустоте. При этом сами значения перегрузок

тіа всем

диапазоне весовых отдачп

цк = 0-1-1 уменьшены,

что

является

следствием влияния атмосферы.

Как и в случае полета

в пусто­

те, величина перегрузки 0пт тем больше, чем больше гравитаци­ онные потери. Области на графике, где «опт^К практического смысла не имеют, поскольку при вертикальном старте с поверх­ ности земли всегда имеет место н0> 1.

0

20

60

60

80

0

 

 

Рис. 139.

Характер

зависимости

Рис. 140. Изменение начальной

осе-

иараметра

тіо

от конечного

угла

вой перегрузки п0 в зависимости от

Ѳк и скорости истечения W0

весовой отдачи с учетом влияния ат­

При т]о = 0 и |= 1

формула

мосферы

(10.6).

(10.18) дает выражение

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

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

10.6.ОПТИМАЛЬНЫЙ ДИАМЕТР КОРПУСА РАКЕТЫ

СУЧЕТОМ ВЛИЯНИЯ АТМОСФЕРЫ

Получение простой зависимости рк = Рк(-Оо) чисто аналитиче­ ским путем не представляется возможным. Наиболее целесооб-

244

разным здесь является метод, основанный на аппроксимации кривых весового анализа, представленных на рис. 141. Такие

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

где значения ро и

Do

 

(10. 19)

 

соответствуют экстремальным точкам кри­

вых рк= Рк(Пм)

при

dyiK/dDM= 0.

AVX(DM),

 

Воспользовавшись формулой (2. 44) и учитывая потери в ско­

рости

за счет

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

 

за­

пишем

-lnК(ö„)l - s m . [ 1 _ ^(Du)) _ y,x D l ( J £ f ,

 

 

 

 

 

«о

\ Do 1

 

где параметр ѵх определяется статистически на базе весового анализа по аналогии определения параметра АѴх (10.16).

Рис. 141. Влияние диаметра и гру­

Рис.

142.

Изменение отношения

зоподъемности на весовую отда­

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

чу Рк

А > _

 

метра D

 

для

 

чи G т

различных значений диа­

Приравнивая нулю производную dVK/d(Du/D q),

получим

2

IN, ( %

) 4+ < а д + 1 ) ( f )

-

 

- ( D J D 0) +

- ^ r 0=0,

(10.20)

 

 

ZriQ

 

245

Решение этого уравнения для случая переменной нагрузки на мидель <7m= <7m(Gt) (по данным весового анализа или зависи­ мостей разд. 5.11) представлено на рис. 142. Как видно из гра­ фика, влияние сопротивления воздуха на уменьшение диаметра

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

малых

 

и боль­

ших значениях весовых отдач

GT. При малых

значениях пара­

метра G T

э т о

происходит вследствие уменьшения потерь

/±.ѴХ =

=

Ѵ х

(GT)

(активный участок

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

атмосфе­

ру),

при

больших значениях GT — вследствие

увеличения на­

грузки на мидель ^m= <7m(Gt)-

И з

графика видно также, что ве­

 

 

D M/D

0) пропорциональна параметру

qM.

 

 

личина (

 

 

 

 

 

Рис. 143. Зависимость соотно­ шения диаметров D M/D0 от па­

раметра От при различных значениях нагрузки на ми­ дель ды

Рис. 144. Влияние параметра qM

на

соотношение

диаметров

D M/Do при

различных значени­

 

ях

перегрузки

п0

Представление о влиянии параметра qM на величину (Au/Аз) дает график (рис. 143), построенный по тому же урав­

нению при различных значениях параметра <7М= const.

 

Из графика видно, что влияние сопротивления

воздуха на

диаметр ракеты особенно сказывается при малых

нагрузках на

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

D M/D0

на участ­

ке весовых отдач Gt= 0,6-h 1 остается

практически постоянной.

Это позволяет положить в уравнении

(10.20) ро = 0 и после ре-

246

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