
книги из ГПНТБ / Москаленко Г.М. Инженерные методы проектирования в ракетодинамике
.pdfфик построен по данным численного интегрирования уравнений движения как
|
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 |
|
|
|
П р и м е ч а н и е . 2ß— центральный угол, Ѳк—конечный угол (в момент вык лючения двигателя), Ѵ к—конечная скорость, 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 |
|
° |
(вк) <?т |
|
||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
■По |
Ж |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
4о |
+ |
— |
|
||
|
|
|
Ңс |
|
|
|
|
|
|
, |
е* |
|
|
|
где |
|
|
|
|
|
|
|
|
|
|Чс |
|
|
||
= 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