- •Фролов, А.Д.
- •ОГЛАВЛЕНИЕ
- •ВВЕДЕНИЕ
- •1. ПОСТАНОВКА ПРОБЛЕМЫ ПАРАМЕТРИЧЕСКОГО ПРОЕКТИРОВАНИЯ
- •РАКЕТ
- •1.1. Предварительные замечания
- •1.2. Сокращения, условные обозначения, индексы
- •1.3. Основные этапы процесса параметрического проектирования
- •2.1. Предварительные замечания
- •2.3. Определение массовых характеристик ракет с РДТТ
- •2.4. Определение геометрических характеристик РДТТ и ракеты
- •2.5. Определение проектно-баллистических параметров РДТТ и ракеты
- •2.6. Определение предельных секундных расходов топлива
- •2.7. Анализ и учет габаритных ограничений РДТТ и ракеты
- •2.8. Аэродинамические характеристики ракеты
- •2.9. Моменты инерции и центровочные характеристики ракеты
- •В) Расчет центровочных и моментных характеристику-й «сухой» субракеты,
- •Сtp(0 = фнавед ” 0 /
- •3.3. Назначение потребной конечной скорости и угла бросания
- •3.5. Проектирование ракеты без оптимизации параметров (Организация работы программы KAMFAD)
- •4. ДЕТЕРМИНИРОВАННАЯ ОПТИМИЗАЦИЯ ПРОЕКТНЫХ ПАРАМЕТРОВ РАКЕТ
- •4.1. Предварительные замечания
- •4.2. Адаптация метода неопределенных множителей Лагранжа
- •4.3. Метод направленного поиска оптимальных параметров
- •Вывод алгоритма решения задачи
- •Выберем X,(r),X2(r),X3(r),X4(r) из уравнений:
- •5. СТОХАСТИЧЕСКАЯ ОПТИМИЗАЦИЯ ПРОЕКТНЫХ ПАРАМЕТРОВ РАКЕТ
- •5.1. Предварительные замечания
- •5.2. Формирование случайной реализации ракеты
- •5.3. Определение основных вероятностных характеристик ракет
- •5.5. Метод направленного поиска оптимальных параметров
- •Графики изменения аэродинамических коэффициентов ракеты:
- •Графики изменения параметров движения ракеты на ПУТ:
- •6.5. Параметрическое проектирование ракет с РДТТ из различных материалов
- •6.13. Частная параметрическая оптимизация секундных расходов твердого топлива двигательными установками баллистической ракеты
- •6.16. Влияние закона распределения случайных величин на статистические параметры дальности полета ракеты
- •6.17. Связь высоты точки старта ракеты с ее эффективностью
- •6.18. Параметрическое проектирование баллистических ракет с твердотопливными двигательными установками различных диаметров
- •7. ЛАБОРАТОРНЫЙ ПРАКТИКУМ
- •7.1. Предварительные замечания
- •7.4. Лабораторная работа № 3.
- •7.5. Лабораторная работа № 4.
- •7.6. Лабораторная работа № 5.
- •7.7. Лабораторная работа № 6.
2.4.Определение геометрических характеристик РДТТ и ракеты
Врамках принятого облика ракеты (тандемная схема), с использованием полученных в результате работы соответствующих модулей значений энергетических и массовых характеристик ДУ и ракеты можно определить основные габаритные характеристики проектируемого объекта для стадии его баллистического проектирования.
Кэтим характеристикам могут быть отнесены:
•длина ПГ ракеты;
•длина ПО ракеты;
•длина ГЧ ракеты;
•длина переходного отсека субракеты;
•длина переднего днища ДУ;
•длина цилиндрической части ДУ;
•длина заднего днища ДУ;
•длина заднего днища ДУ с сопловым блоком;
•длина соплового блока;
•длина снаряженной ДУ;
•длина субракеты.
Естественно, значения этих характеристик подлежат уточнению при окончательной
конструкторской проработке, осуществляемой специалистами или соответствующими компьютерными модулями, разработанными опять же при участии этих специалистов.
Основные геометрические характеристики ракет с РДТТ определяются с учетом рекомендаций, данных специалистами [28, 31, 35, 40, 51]. Расчеты выполняются в следующей последовательности, реализующей принятую математическую модель:
Вычисление длины ГЧ ракеты:
•длина конусной части ПГ
^пг кон = ]-21ППГК КОц/(H D ПгУпг)}
•длина цилиндрической части ПГ
U\YЦИЛ“ 4ТПпг-^цил/(jiD nrYnr)j
•полная длина ПГ
/цГ ~ AirКОН ^ПГЦИЛ;
•длина ПО
/по l2mno/[n(D 1+ D\Dnr+ D2пг)];
• длина ГЧ (суммарная длина ПГ и ПО ракеты)
/гч /nr^"/noj
Вычисление длины Д У и субракет;
длина переходного отсека, соединяющего ГЧ с 1-й (верхней) ступенью
/поте 1 = /дн пер \К п отс 1>
•длинау-го переходного отсека, соединяющего ступени
/п отсJ ~ (/дн перJ /дн зад полнj)K n отсу>
у= 2, 3,
•длина переднего и заднего днищу-й ДУ
/дн перj ~ /дн зад
•длина соплау-й ДУ
/соплj = ^Kpy/[2(sy-l)ctg(pconn / 2 ) ] ;
•длина цилиндрической частиу-й ДУ
/цилу = т Тцилу/(Уту*^мy-^Т ЦИЛу)»
•полная длина заднего днища с сопловым блокому-й ДУ
/дн зад полнУ^дн заду /соплу-^с дну>
•полная длинау-й ступени
/дуу = /дн перу /цилj |
/дн зад полну» |
•длина 1-й субракеты
/субр 1 = /гч“^ /п отс 1 /ду 1 “ /дн пер 1j
•длинау-й субракеты (/' ^ 1)
/субру = (/субру-1 " /дн зад полну-0 /п отсу /дуу “ /дн перу»
Теоретический чертеж трехступенчатой ракеты, соответствующий приведенному алгоритму, отражен на рис. 2.1:
2.5. Определение проектно-баллистических параметров РДТТ и ракеты
Основные проектные параметры ракеты с сопутствующими характеристиками, полученными в результате работы алгоритмов подразделов 2.2- 2.4, полностью определяют весь облик проектируемой под ТЗ ракеты. Сам факт удовлетворения полученных характеристик ракеты предъявляемым требованиям должен быть подтвержден расчетами, одним из вариантов которых является приведенный в настоящем подразделе. Рассчитываются основные баллистические и некоторые проектные характеристики ракет с РДТТ:
•среднетраекторный угол наклона вектора скорости к стартовому горизонту на каждом АУТ; начальная и конечная перегрузка на каждом АУТ;
углы бросания в конце каждого АУТ:
-реально сформированный;
-оптимальный для достижения максимальной дальности;
начальная нагрузка на Мидель в начале каждого АУТ; относительная конечная масса ракеты на каждом АУТ; среднетраекторные потери скорости на каждом АУТ, вызванные:
-влиянием силы притяжения Земли;
-аэродинамическим сопротивлением атмосферы;
-изменением тяги с высотой;
скорость Циолковского в конце каждого АУТ; конечная скорость на каждом АУТ;
прогнозируемая дальность полета оставшейся части ракеты в зависимости от конечной скорости в конце каждого АУТ.
Основные баллистические и некоторые из основных проектных характеристик ракет с РДТТ определяются с учетом рекомендаций специалистов [2,4,14,28, 35].
Последовательность расчетов соответствует порядку перечисления зависимостей: - тягау'-й ДУ в пустоте
Rnj Jcnnj rhj g 0\
- продольная перегрузка в концеу-го АУТ
** ■
^XKj
-продольная перегрузка в началеу-го АУТ (/ * р)
К.
j =
”*0180 ’
-продольная перегрузка в началеp-то АУТ (т.е. на старте)
?^ipPhq .
m o p S o
- конечная относительная масса ракеты нау-м АУТ
- начальная нагрузка на Мидель нау-м АУТ
р . m o j8 o .
гм0у - « > Лм/
-скорость Циолковского в концеу-го АУТ
^= -£ о Л ш П п( М ;
-среднетраекторный угол наклона вектора скорости к стартовому горизонту нау-м АУТ:
•в случае, если степень параболы а > 1,
(вр - 9,)[(f, |
- ft,- t 0JY+'] |
Gc " ^о) |
у ^0у) |
• в случае, если степень параболы а < 1 |
|
(в0- е кЖ ^ |
(toj -Л ) Г 1] |
выражения получены интегрированием на у-м АУТ параболы а - й степени, задающей программу движения ракеты:
•в случае, если степень параболы а > 1,
•в случае, если степень параболы а < 1,
полагается, например, в приведенных соотношениях при р = 3:
- угол бросания в концеу-го АУТ:
•в случае, если степень параболы а > 1,
•в случае, если степень параболы а < 1,
-потери скорости на преодоление притяжения Земли (гравитационные потери) наj - м АУТ
-потери скорости на преодоление аэродинамического сопротивления атмосферы на j-u АУТ
AVCxJ = -KCxJAVCxp.
Здесь AVcxp определяется следующим образом [4]. Большое количество расчетов, проведенных в целях определения коэффициента, характеризующего сопротивление атмосферы:
позволило установить эмпирическую зависимость h от Vp и а:
V = V + AV
¥ Р Уир ^
<y = T^u^sineK-10'},
где
Т = тйр1тр,
1 |
“0р |
|
|
|
и |
= —- |
“ |
2° (^ад°я |
|
мср |
|
|
||
ЛдО/7 = Л д п р " |
|
1 171Р ’ ^ |
0332 атм. |
|
|
|
Эта зависимость интеграла h от и а приведена в табл. 2.1:
Таблица 2.1
Зависимость интеграла h от Vp и а
a \ V p |
750 |
1000 |
1500 |
2000 |
2500 |
3000 |
3500 |
4000 |
4500 |
5000 |
2 ,0 |
1180 |
1680 |
2750 |
3460 |
3860 |
4080 |
4180 |
4250 |
4310 |
4350 |
2 ,2 |
1100 |
1500 |
2430 |
3000 |
3310 |
3440 |
3500 |
3540 |
3580 |
4000 |
2,4 |
1020 |
1400 |
2200 |
.2610 |
2810 |
2930 |
3000 |
3030 |
3040 |
3060 |
2,6 |
980 |
1330 |
1950 |
2300 |
2480 |
2550 |
2600 |
2610 |
2620 |
2625 |
2,8 |
920 |
1250 |
1780 |
2030 |
2150 |
2210 |
2260 |
2300 |
2320 |
Г2340 |
3,0 |
830 |
1190 |
1610 |
1820 |
1920 |
Г980 |
2000 |
2005 |
2010 |
2015 |
3,2 |
810 |
1080 |
1480 |
1650 |
1750 |
1800 |
1800 |
1800 |
1800 |
1800 |
3,4 |
800 |
1000 |
1310 |
1500 |
1580 |
1600 |
1600 |
1600 |
1600 |
1800 |
3,6 |
760 |
930 |
1200 |
1350 |
1400 |
1420 |
1420 |
1420 |
1420 |
1420 |
4,0 |
700 |
830 |
1020 |
1120 |
1160 |
1170 |
1170 |
1170 |
1170 |
1170 |
4,4 |
610 |
730 |
900 |
990 |
1000 |
1000 |
1000 |
1000 |
1000 |
1000 |
4,8 |
580 |
680 |
800 |
860 |
880 |
880 |
880 |
880 |
880 |
880 |
5,2 |
520 |
600 |
710 |
750 |
780 |
780 |
780 |
780 |
780 |
780 |
5,6 |
1'430 |
520 |
620 |
660 |
660 |
660 |
660 |
660 |
660 |
660 |
U,V60 |
400 |
470 |
560 |
580 |
580 |
580 |
580 |
580 |
580 |
580 1 |
6 4 |
360 |
400 |
480 |
480 |
480 |
480 |
480 |
480 |
480 |
480 |
6,8 |
310 |
350 |
410 |
410 |
410 |
410 |
410 |
410 |
410 |
410 |
Полученное из табл. 2.1 значение /2позволяет определить искомые потери:
So^Qp^vOp т
Л ^ р = -
|
Р |
2 |
|
Г иОр |
|
|
° ' = *. |
|
К |
'Л дпр^р |
S a p P h b ■ |
Потери скорости на изменение тяги ДУ с высотой на j - м АУТ определяются зависимостью
у = 1,2...... |
р - 1. |
Величина AVSapможет быть найдена следующим образом [4]. Вычисляется значение интеграла 1у.
послечего определяется значение AVsap-
AV^ = -go(Janp~^<*Op^y
Значение величины r\ =fitp,vop) определяется из табл. 2.2, полученной эмпирическим путем:
Таблица 2.2
Зависимость величины т| от tp и VQP
V0p\t,, |
0 |
10 |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
90 |
100 |
0,2 |
0 |
4,5 |
8,4 |
10,0 |
11,0 |
11,7 |
12,0 |
12,5 |
12,7 |
12,9 |
13,0 |
0,3 |
0 |
3,4 |
6,7 |
8,4 |
9.6 |
10,5 |
11,1 |
11,5 |
11,7 |
11,8 |
11,9 |
0,4 |
0 |
2,8 |
5,3 |
7,1 |
8,3 |
9,2 |
9,7 |
10,3 |
10,6 |
10,8 |
11,0 |
0,5 |
0 |
2,2 |
4,3 |
6,0 |
7,3 |
8,3 |
9,1 |
9,6 |
10,0 |
10,3 |
10,4 |
0,6 |
0 |
1,9 |
3,7 |
5,3 |
6,6 |
7,6 |
8,4 |
9,0 |
9,5 |
9,8 |
10,0 |
0,7 |
0 |
1,7 |
3,2 |
4,7 |
6,0 |
7,0 |
7,9 |
8,6 |
9,2 |
9,5 |
9,7 |
0.8 |
0 |
1,4 |
2,8 |
4,1 |
5,4 |
6,5 |
7,4 |
8,2 |
8,8 |
9,3 |
9 ,6 |
0,9 |
0 |
1,1 |
2,3 |
3,5 |
4,7 |
5.9 |
6,9 |
7,8 |
8,6 |
9,1 |
9,5 |
- скорость ракеты в концеj - го АУТ
Кг = у*чг+ |
+ AV&+ AVc*j+ AVsaj> |
V*P =0’ |
J = |
- дальность полета оставшейся части ракеты в концеу'-го АУТ
определяется из табл. 3.2;
~ оптимальный угол бросания для скорости V& обеспечивающей максимум дальности полета
^копту “ /е(^к/)> J Ъ |
Р9 |
определяется также из табл. 3.2.
Необходимо заметить, что табличные значения зависимостей f L(V^) и f e(V
получены из соотношений (19.14) и (19.32) [4] при допущении RK= Яо(подразд. 3.3).