- •Фролов, А.Д.
- •ОГЛАВЛЕНИЕ
- •ВВЕДЕНИЕ
- •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.
6.5. Параметрическое проектирование ракет с РДТТ из различных материалов
Параметрическое проектирование ракет минимальной стартовой массы с применением и без применения композитных материалов
Известно, что условия работы отдельных элементов и ракеты в целом характеризуются высокими механическими и температурными нагрузками, возникающими как в процессе эксплуатации, так и в полете, а также коррозионным воздействием окружающей среды в процессе длительного процесса эксплуатации. К тому же необходимо обеспечить минимальную массу конструкции с учетом наличия сырьевой и промышленной базы, стоимость производства и т.п. Отсюда требуются высокая механическая прочность и жесткость при минимальной плотности, сохранение в течение всего периода эксплуатации характеристик в требуемом диапазоне, высокая коррозионная стойкость, технологичность и недефицитность исходных материалов. Определяющими критериями при выборе материалов являются, в частности, предел прочности ав, удельная масса материала у и их отношение - удельная прочность Ов/у. Характеристики применяемых конструкционных материалов [29, 84] приведены в табл. 6.5:
|
|
|
|
Таблица 6.5 |
|
|
Характеристики применяемых конструкционных материалов |
||||
№ |
|
Предел |
Удельная |
Удельная |
|
Материал |
прочности о„ |
масса у, |
|||
варианта |
прочность, см |
||||
|
кгс/см2 |
кг/см* |
|||
|
|
|
|||
1 |
Алюминиево-магниевый сплав |
3250 |
0,0027 |
1 203 703 |
|
2 |
Алюминиево-бериллиево- |
4500 |
0,0024 |
1 875 000 |
|
|
магниевый сплав |
|
0,00785 |
2 165 605 |
|
3 |
Малоуглеродистая сталь |
17 000 |
|||
4 |
Титановый сплав |
10 500 |
0,0046 |
2 282 608 |
|
5 |
Сталь |
19 000 |
0,0078 |
2 435 897 |
|
6 |
Высокопрочная сталь |
21 000 |
0,008 |
2 625 000 |
|
7 |
Стеклопластик на основе |
10 500 |
0,00205 |
5 121 951 |
|
|
стекловолокон |
|
|
5 500 000 |
|
8 |
Стеклопластик |
11 000 |
0,002 |
||
9 |
Углепластик |
13 000 |
0,0015 |
8 666 666 |
|
10 |
Органопластик |
14 000 |
0,0014 |
10 000 000 |
|
11 |
Композиционный материал |
16 000 |
0,0017 |
9 411765 |
Частные исследования для одной и той же ракеты (L ^ = 10 000 км, тпг - 1700 кг) с оптимизацией по критерию дальности полета позволили провести сравнительную оценку полученных вариантов 3-ступенчатой ракеты минимальной стартовой массы. Основные результаты расчетов приведены в табл. 6.6:
|
|
|
Таблица 6.6 |
|
Результаты параметрического проектирования вариантов ракет |
||||
из металла и композиционного материала |
Отношение |
|||
Параметр |
Материал корпуса, днищ, соплового раструба |
|||
Металл |
Композит |
(композит/металл) |
||
|
||||
Характеристики материала: |
|
0,0017 |
|
|
- удельная масса у, кг/см3 |
0,0080 |
0,2125 |
||
- удельная прочность а„ кгс/см2 |
25000 |
16000 |
0,6400 |
|
Масса то, кг: |
60512 |
47801 |
0,7900 |
|
- 1-й субракеты (ракеты на |
||||
старте) |
23467 |
16076 |
0,6851 |
|
- 2-й субракеты |
6235 |
5419 |
0,8691 |
|
- 3-й субракеты |
36902 |
31682 |
0,8585 |
|
- стартовой ДУ |
17140 |
10631 |
0,6202 |
|
-2-Й ДУ |
4445 |
3662 |
0,8238 |
|
- 3-й ДУ |
|
|
|
|
|
|
Окончание табл. 6.6 |
|
Параметр |
Материал корпуса, днищ, соплового раструба |
Отношение |
||
Металл |
Композит |
(композит/металл) |
||
Коэффициент весового |
||||
|
|
|
||
совершенства аду: |
|
|
|
|
- 1-й ДУ (стартовой ДУ) |
0,0821 |
0,1123 |
1,3678 |
|
-2-йДУ |
0,0687 |
0,0641 |
0,9330 |
|
- 3-й ДУ |
0,0951 |
0,0450 |
0,4732 |
|
Характеристики конца АУТ: |
|
|
|
|
- скорость VKi м/с |
6890 |
6408 |
0,9300 |
|
- высота Як, м |
245128 |
382458 |
1,5602 |
|
- координата х к, м |
436 983 |
623918 |
1,4278 |
|
- угол бросания 0 К, град |
19,23 |
20,86 |
1,0847 |
|
Длина ракеты на старте /0, м |
28,17 |
26,04 |
0,9244 |
|
Потери скорости ДКк, м/с |
-1085 |
-1540 |
1,4194 |
|
Полное время АУТ / с |
162 |
234 |
1,4444 |
Главным результатом применения композиционных материалов можно считать снижение стартовой массы ракеты на ~ 21 %.
Влияние удельной прочности конструкционных материалов на основные характеристики твердотопливной баллистическойракеты
Приведены результаты параметрического проектирования оптимальных ракет с РДТТ, доставляющих полезный груз тпг = 1700 кг на дальность 1 = 10 000 км в условиях применения различных конструкционных материалов обечаек и днищ двигательных установок. Диапазон исследуемых значений удельной прочности конструкционных материалов 5-60 км. Очевидно, что каждое конкретное значение удельной прочности может теоретически реализоваться бесконечным числом сочетаний пределов прочности и удельных масс конструкционных материалов (на практике реально ограниченных). Результаты параметрического проектирования ракет представлены на рис. 6.76-6.78:
Удельная прочность, км
Рис. 6.76. Изменение массы ракет (ось ординат, кг)
в зависимости от удельной прочности конструкционных материалов (ось абсцисс, км)
120000
. 100000
. 80000
3 60000
g40000
20000
in l f > l 0 c N < 0 0 ^ c 0 c g < 0 ° * 00 см со о
ф з з ' ч- т- с м с м с м с о с о ^ ^ ^ - ю ю с о
Удельная прочность, км
Рис. 6.77. Изменение запасов топлива ракет (ось ординат, кг) в зависимости от удельной прочности конструкционных материалов (ось абсцисс, км)
Рис. 6.78. Изменение давления в ДУ (ось ординат, атм) стартовой, средней и верхней ступеней ракет (сверху вниз) в зависимости от удельной прочности конструкционных материалов (ось абсцисс, км)
Из расчетов следует, что с увеличением удельной прочности конструкционна/ материалов уменьшаются начальные массы субракет и запасы твердого топлива в ДУ ? условиях перераспределения этих запасов между ступенями. Давления (а с ними г удельные пустотные импульсы ДУ) растут с увеличением удельной прочное-^ конструкционных материалов.
Параметрическое проектирование ракет с РДТТ из нетрадиционных материалов
В теоретическом плане представляет интерес изучение возможности создан^ ракеты с РДТТ из нетрадиционных материалов. В результате расчетов были получен?* оптимальные трехступенчатые ракеты с минимальными стартовыми массам**» доставляющие полезный груз массой 1700 кг на одну и ту же дальность 10 000
Результаты расчетов сведены в табл. 6.7:
Таблица бJ Результаты параметрического проектирования ракет из нетрадиционных
конструкционных материалов
|
Характери- |
|
Численные значения характеристик для конструкционных материалов типа: |
|
|||||||
стики ракеты |
Чугун |
Медь |
Сталь |
Бронза |
Латунь |
Сосна |
Алюм. |
Титан |
Ст/вол |
Fe-HK |
|
|
ст^МПа |
387 |
510 |
510 |
612 |
816 |
85 |
480 |
2548 |
1631 |
13400 |
|
Y, кг/см3 |
0,007 |
0,009 |
0,008 |
0,007 |
0,009 |
0,001 |
0,003 |
0,008 |
0,002 |
0,008 |
|
т 0> Т |
143,2 |
128,5 |
127,4 |
104,1 |
99,6 |
86,6 |
70,3 |
61,8 |
48,3 |
47.7 |
|
ГПг.Т |
127,3 |
113,6 |
113,7 |
91,9 |
88,5 |
78,0 |
62,9 |
55,4 |
43,2 |
43,0 |
|
6щ1Х 1 |
1,09 |
0,88 |
0,82 |
1,00 |
0,71 |
6,40 |
1,76 |
0,44 |
1,20 |
0,14 |
|
^max * М |
2,88 |
2,76 |
2,76 |
2,40 |
2,48 |
2,28 |
2,04 |
1,88 |
1,76 |
1,80 |
|
/о,м |
34,0 |
32,5 |
31,4 |
30,9 |
30,0 |
29,3 |
28,7 |
27,7 |
25,2 |
24,2 |
|
* |,с |
67,6 |
55,3 |
57,7 |
54,0 |
67,3 |
51,4 |
60,3 |
50,1 |
63,2 |
57,9 |
|
/иТ1 кг/с |
69,6 |
70,8 |
73,3 |
68,0 |
72,8 |
59,5 |
67,7 |
61,9 |
60,2 |
55,8 |
|
/VI. МПа |
1,32 |
1,42 |
1,63 |
1,88 |
1,73 |
2,17 |
2,62 |
3,78 |
5,92 |
6,50 |
|
Л |, МПа |
0,03 |
0,03 |
0,03 |
0,03 |
0,03 |
0,02 |
0,03 |
0.02 |
0,02 |
0,02 |
_ |
т2,с |
70,0 |
70,0 |
45,8 |
62,8 |
55,6 |
58,6 |
54,6 |
57,5 |
47,2 |
38,7 |
|
mj2 f кг/с |
262,1 |
332,3 |
282,3 |
313,4 |
222,5 |
233,6 |
245,9 |
212,8 |
158,6 |
149,8 |
|
Р а , МПа |
1,25 |
1,27 |
1,63 |
1,64 |
2,73 |
2,17 |
3.14 |
4,54 |
10,2 |
9,34 |
|
р ш2 , МПа |
0,07 |
0,06 |
0,06 |
0,06 |
0,05 |
0,06 |
0,06 |
0,05 |
0,03 |
0,03 |
|
* з,с |
69,0 |
49,5 |
56,1 |
47,7 |
50,4 |
56,0 |
48,6 |
50,5 |
54,2 |
56,8 |
|
т тз> кг/с |
1511 |
1744 |
1720 |
1437 |
1413 |
1094 |
933 |
793 |
603 |
599 |
|
Ркз, МПа |
1,59 |
1,73 |
1.-60 |
2,78 |
2,54 |
2,60 |
4,50 |
6,50 |
12,2 |
11,2 |
|
Ра, МПа |
0,10 |
0,13 |
|
0,13 |
0,12 |
0,07 |
0,10 |
0,08 |
0,06 |
0,06 |
|
Не составляет большого |
труда определить |
положение |
точек |
материалов с их |
удельными прочностями из табл. 6.7 на рисунках 6.76—6.78 и сделать выводы.
Результаты расчетов, приведенные в табл. 6.7, позволяют сделать заключение о том, что теоретически возможно создание ракет, решающих одну и ту же задачу, из различных
материалов. Однако диапазоны изменения характеристик полученных 3—ступенчатых ракет варьируются в весьма широких пределах (стартовая масса и суммарный запас
топлива |
- в 3 раза, «сухая» масса ракеты - в 3,4 раза, длина ракеты на старте |
" |
1,4 |
раза, |
|
максимальный диаметр ракеты - |
в 1,6 раза, максимальная толщина обечайки |
' |
в 46 |
раз). |
|
В |
целом непреодолимыми |
для нетрадиционных материалов, используемых в |
качестве конструкционных, могут стать как технические, так и экономические проблемы^ Отсюда безусловными способами совершенствования конструкции ракет с Р Д Т Т остаются традиционные уменьшение удельной массы у и увеличение предела прочности а в конструкционных материалов.
6.6. Влияние количества ступеней ракеты с РДТТ на ее основные характеристики
Компьютерная программа KAMFAD [48, 49, 50] параметрического проектирования ракет с РДТТ в детерминированной постановке позволяет получать оптимальные, например по критерию дальности, проектные параметры. В условиях одного и того же заданного полезного груза тпг= 1700 кг, одной и той же заданной дальности 1зад = W 000 км и одинаковых пределов начальных и конечных перегрузок были получены оптимальные варианты составных ракет тандемной схемы с числом ступеней р = 2, 3 и 4 (одноступенчатую ракету в условиях принятой математической модели получить невозможно, а ракеты с числом ступеней больше 4 пока не представляют практического интереса). При этом заданная дальность полета ракеты 1зад обеспечивалась в условиях оптимизации проектных параметров подбором исходной для проектирования дальности 1исх. Результаты расчетов приведены в табл. 6.8:
Таблица 6.8
Результаты параметрического проектирования ракет с различным числом ступеней
Характеристики |
|
Вариант ракеты (количество и номера ступеней и субракет) |
|
||||||
двигательных |
_____ Р г 2______ |
|
Р шз |
|
_______________ Р = А _______________ |
||||
установок и субракет |
1 |
2 |
1 |
2 |
3 |
1 |
2 |
3 |
4 |
/лоь кг |
65 411 |
10 370 |
60 960 |
21 790 |
5786 |
58 650 |
39 453 |
12 528 |
5798 |
mrh кг |
51 858 |
7992 |
36 158 |
14 899 |
3636 |
16 746 |
25 039 |
6066 |
3623 |
|
0,2005 |
0,2173 |
0,4002 |
0,3071 |
0,3589 |
0,7123 |
0,3594 |
0,5095 |
0,3619 |
VI |
0,0584 |
0,0718 |
0,0793 |
0,0681 |
0,0992 |
0,1374 |
0,0709 |
0,0986 |
0,1058 |
Т/,с |
59,2 |
71,8 |
45,9 |
63,1 |
55,6 |
11,5 |
43,0 |
30,0 |
47,9 |
mTjt кг/с |
876,2 |
111,3 |
787,7 |
236,2 |
65,4 |
1463,6 |
582,8 |
202,6 |
75,6 |
атм |
44,3 |
32,8 |
63,5 |
41,7 |
36,9 |
108,8 |
52,9 |
63,3 |
44,1 |
А . атм |
0,85 |
0,30 |
0,.85 |
0,51 |
0,22 |
1,88 |
0,92 |
0,76 |
0,25 |
Д , см |
208 |
112 |
192 |
132 |
100 |
180 |
160 |
100 |
100 |
ЮтЕ, кг |
59 850 |
|
54 693 |
|
|
51474 |
|
||
ЮсухХ, КГ |
5561 |
|
6267 |
|
|
|
7176 |
|
|
/о, см |
2765 |
|
2270 |
|
|
|
3175 |
|
|
Vxy м/с |
6964 |
|
6893 |
|
|
|
6932 |
|
|
*«,м |
353 664 |
|
448 621 |
|
|
384 608 |
|
||
Я „ м |
188 278 |
|
251 263 |
|
|
211440 |
|
||
КМ |
8500 |
|
8500 |
|
|
|
9160 |
|
|
^всал. КМ |
10 060 |
|
10 030 |
|
|
10 040 |
|
Результаты расчетов позволяют сделать выводы о том, что с увеличением количества ступеней для ракет с оптимальными проектными параметрами, доставляющих полезный груз на одну и ту же заданную дальность, уменьшаются как стартовая масса,
так и суммарные запасы топлива; одновременно увеличивается «сухая» масса ракеты. Длина ракеты на старте обусловлена ограничениями как по предельно допустимому
уГ ^ |
, ? У; 1“ , " ю™ |
п и |
и ком„о»о.кой р » ™ J цело». Необходшо |
||
|
|
|
“""'"Г |
«* <tfOOS лл» 2- |
|
6.7. |
Баллнстнчеркое проектировав, ракеты с мггамазаинао |
||||
|
|
™««ты вопадавня а ™ |
|
■“ » |
|
10 000 |
ш . Д е т е р м |
я ^ о ^ " ^ ^ " ? ^ |
|
>500кт надаланоет, |
|
ракеты (количество итераций равно 50) в самках nnrrvupn к |
рию конечной скорости |
||||
дальность до 12 000 км. |
Р |
^ |
Ной стартовой массы увеличила |
Далее полученная оптимальная по детерминированному критерию ракета подверглась оптимизации по двум вариантам стохастического критерия J - вероятности и частоте попадания в цель (заданный диапазон скоростей VKmin-VKmax). Оптимизации подлежали проектные параметры РДТТ - времена работы TJ, секундные расходы топлива mTj, давления в ДУ рку и на срезах сопловых блоков p ^ ,j = 1,2, 3.
Результаты расчетов на ЭВМ (приводятся результаты оптимизации только для давлений рк у и рй j, j = 1, 2, 3, дающих более резкие экстремумы по сравнению с распределением топлива по ступеням ракеты, и их производных - коэффициентов влияния вида Jp i = dJ/дрк1, атм""1), отражены в табл. 6.9:
Таблица 6.9 Результаты оптимизации по критериям вероятности или частоты попадания в цель
|
|
|
Количество статиспытаний |
z = ;180 |
||||
Параметр |
z = |
12 |
z - |
60 |
* = 1 2 0 |
|||
|
Вероятн. |
Частота |
Вероятн. |
Частота |
Вероятн. |
Частота |
Вероятн. |
Частота |
№ опт. ит.—► |
19 |
21 |
21 |
20 |
20 |
11 |
6 |
16 |
Ржи атм |
37,83 |
38,53 |
38,78 |
38,06 |
39,70 |
38,31 |
36,46 |
37,14 |
У ^ .а т м '1 |
-2,67 |
0,00 |
0,00 |
-0,02 |
0,02 |
-0,02 |
0,00 |
-0,15 |
р * 2 , атм |
45,56 |
44,39 |
44,28 |
44,25 |
44,02 |
47,03 |
45,13 |
45,00 |
Jpa> а™ '' |
-0,07 |
0,01 |
0,04 |
-0,17 |
0,06 |
0,00 |
-0,04 |
0,14 |
|
|
|
|
84,27 |
78,09 |
79,76 |
86,29 |
|
0*з, атм |
90,74 |
78,43 |
82,39 |
88,09 |
||||
У ^ .а т м '1 |
0,00 |
0,06 |
0,08 |
-0,02 |
0,05 |
0,00 |
0,02 |
0,07 |
0аь атм |
0,27 |
0,21 |
0,26 |
0,27 |
0,27 |
0,24 |
0,25 |
0,27 |
^ 1 ,а т м -’ |
-3,06 |
7,50 |
-2,97 |
-3,33 |
1,46 |
-5,85 |
-6,66 |
2,29 |
0 ,2, атм |
0,56 |
0,56 |
0,40 |
0,54 |
0,58 |
0,48 |
0,52 |
0,58 |
J ^ , атм’1 |
-1,65 |
0,06 |
-4,00 |
-0,71 |
0,36 |
-3,49 |
-5,93 |
1,52 |
0*3, атм |
0,68 |
0,95 |
0,72 |
0,75 |
0,72 |
0,88 |
0,74 |
0,77 |
У ^ .а т м '1 |
-0,73 |
0,52 |
89,12 |
-3,74 |
0,65 |
0,00 |
-1,67 |
89,50 |
|
11676 |
12 108 |
И 775 |
11617 |
11 749 |
12213 |
11 527 |
11 504 |
mi укм |
11 531 |
11416 |
11431 |
11 523 |
11495 |
11483 |
11 256 |
11 533 |
c L , км |
199 |
185 |
180 |
182 |
220 |
223 |
218 |
221 |
К ном , м/с |
7 170 |
7 194 |
7 162 |
7 157 |
7 185 |
7 201 |
7 128 |
7 145 |
туК1м/с |
7 153 |
7 141 |
7 143 |
7 152 |
7 149 |
7 148 |
7 125 |
7 153 |
СТик, м/с |
20 |
20 |
18 |
18 |
22 |
22 |
23 |
22 |
|
0,4370 |
0,5833 |
0,4399 |
0,5000 |
0,3706 |
0,3917 |
0,3576 |
0,3944 |
_________ i m ------------- |
0,9258 |
0,9167 |
0,9246 |
0,9667 |
0,8938 |
0,9000 |
0,6496 |
0,9000 |
|
|
|
|
|
|
|
|
Эти данные позволяют сделать вывод о слабой зависимости результатов стохастической оптимизации как от числа статиспытаний z, так и от вариантов критериев оптимизации.
Из расчетов следует, что вполне удовлетворительные результаты имеют место при относительно небольших объемах статиспытаний в рамках обоих критериев. При этом необходимо заметить, что результаты оптимизации по критерию частоты попадания в цель являются более универсальными, чем вероятность попадания в цель, поскольку в этом случае отпадает необходимость формирования композиции теоретической формы критерия J в функции исходных законов распределения случайных характеристик параметров, формирующих объект исследования - ракету с РДТТ.
Эти выводы справедливы для стохастической оптимизации параметров как термостатированной, так и нетермостатированной ракеты.
6.8. О теоретической возможности создания унифицированной системы ракетного вооружения на базе трехступенчатой ракеты с РДТТ
В рамках идеи об унифицированной системе ракетного вооружения (частично реализованной в двух ракетах средней дальности на базе ракеты 8К98: 8К96-РТ15 из 2-й и 3-й ступеней на 2500 км и 8К97-РТ25 из 1-й и 3-й ступеней на 6000 км) по компьютерной программе КАМА02 [46] для всех возможных сочетаний ступеней исходной номинальной и оптимальной ракет, полученных в результате параметрического проектирования ракеты KAMFAD [48, 49, 50] с массой полезного груза тпг = 1800 кг на дальность 9300 км, были проведены компьютерные исследования. Результаты этих расчетов отображены в табл. 6.10 (1-я ступень - стартовая, 2-я - средняя, 3-я - верхняя; вариант «1+2+3» означает исходную 3-ступенчатую ракету):
Таблица 6.10
Результаты исследования теоретической возможности создания унифицированной системы ракетного вооружения
Варианты исходных |
Начальная масса, |
|
Вариант сочетания с-пдтеней ракет |
|
|||||
ракет |
дальность |
1+2+3 |
2+3 |
1+3 |
1+2 |
1 |
2 |
3 |
|
|
той кг |
74 825 |
22 169 |
59 093 |
70 063 |
54 250 |
17 325 |
6503 |
|
Номинальный |
т02, кг |
22169 |
6503 |
6503 |
17 559 |
- |
- |
- |
|
/Поз, кг |
6503 |
- |
- |
- |
- |
- |
- |
||
|
|||||||||
|
L, км |
9300 |
1840 |
5730 |
6160 |
3200 |
1220 |
190 |
|
|
/Л01, кг |
74 825 |
31 777 |
49 485 |
70 066 |
44 628 |
26 920 1 |
6499 |
|
Оптимальный |
/П02, кг |
31 777 |
6499 |
6499 |
27 188 |
- |
- |
- |
|
/Поз, кг |
6499 |
- |
- |
- |
- |
- |
- |
||
|
|||||||||
|
L, км |
11350 |
4100 |
5630 |
6960 |
2733 |
2700 |
420 |
Очевидно, что в результате унификации полученными 6 вариантами ракет [68] из исходной номинальной ракеты перекрывается диапазон дальностей 190-6160 км (т.е. около 6000 км); из исходной оптимальной ракеты - диапазон 420-6960 км (т.е. около 6500 км). Разумеется, практическое применение из рассмотренных вариантов может найти только их небольшая часть, например, упомянутая по ракете 8К98.
6.9.Частная оптимизация РДТТ в составе баллистической ракеты
Впроцессе исследования проектных параметров ракет с РДТТ возможна ситуция, когда оптимизировать необходимо параметры отдельной двигательной установки, что характерно для процесса создания ракеты в кооперативном режиме различными организациями. Алгоритмы и соответствующая компьютерная программа KAMFAD [48, 49, 50] позволяют осуществлять такие исследования при наличии информации по
«замороженным» (обозначены индексом «*») характеристикам РДТТ.
В рамках неизменной стартовой массы то = 68 600 кг 3-ступенчатой ракеты, доставляющей полезный груз массой 1700 кг на дальность 10 000 км, была проведена как общая, так и раздельная оптимизация проектных параметров для каждой из трех РДТТ по критерию дальности. При этом «свободные» параметры, обеспечивающие неизменность стартовой массы ракеты, - соответствующие времена работы РДТТ ту, j = 1, 2, 3). Результаты расчетов приведены в табл. 6.11:
Таблица 6.11 Результаты частной параметрической оптимизации параметров ракет_________
Характеристика |
|
Численные значения характеристик для вариантов: |
|
|||
Исходный |
Оптимальный |
Оптимальный |
Оптимальный |
Оптимальный |
||
РДТТ |
||||||
|
|
для верхней ДУ |
для средней ДУ |
для нижней ДУ |
||
|
|
|
||||
Т|.С |
98,33 |
54,46 |
67,20 |
98,33* |
98,33 + |
|
/йТ1 кг/с |
39,32 |
55,47 |
60,11 |
39,32* |
39,32* |
|
Рхи атм |
60,00 |
46,82 |
33,92 |
60,00* |
60,00* |
|
Pal. атм |
0,200 |
0,171 |
0,203 |
0,200* |
0,200* |
|
TiC |
78,69 |
56,99 |
78,69+ |
57,28 |
78,69* |
|
т Т2> кг/с |
166,56 |
203,51 |
166,56* |
232,64 |
166,56+ |
|
Аа, атм |
60,00 |
54,13 |
60,00* |
48,49 |
60,00* |
|
Рв2 , атм |
0,328 |
0,403 |
0,328+ |
0,495 |
0,328* |
|
тз.с |
69,81 |
53,23 |
69,81* |
69,81* |
53,15 |
|
/wT3f кг/с |
636,35 |
874,98 |
636,35* |
636,35+ |
824,09 |
|
р хз, атм |
60,00 |
72,61 |
60,00* |
60,00* |
77,70 |
|
Раз» атм |
0,633 |
0,857 |
0,633* |
0,633* |
0,829 |
|
I , км |
10 058 |
11 772 |
10 974 |
10 389 |
10 444 |
Данные табл. 6.11 позволяют в рамках одинаковой математической модели для всех двигателей ракеты заметить, что:
-с одной стороны, полная оптимизация ракеты дает прирост дальности 1714 км; оптимизация параметров только верхней РДТТ - 916 км, средней - 331 км, стартовой - 386 км; отсюда с точностью до 5% результаты раздельной оптимизации (916 + 331 + 386 = 1633) совпадают с результатами полной оптимизации;
-с другой стороны, значения оптимальных проектных параметров РДТТ для вариантов полной и частной оптимизации отличаются до 28% для верхнего РДТТ, до 23% - для среднего РДТТ и до 7% - для стартового РДТТ.
Очевидно, что максимальный эффект дает полная одновременная оптимизация всех проектных параметров ракеты. При раздельной оптимизации неизбежны потери возможной максимальной дальности.
Влияние отдельных проектных параметров РДТТ нарезультаты детерминированной оптимизации ракеты
В условиях неизменного стартового веса 3-ступенчатой баллистической ракеты $ РДТТ с массой полезного груза тПГ = 1700 кг и номинальной дальностью 10 000 ю* (неизменность обеспечивается «свободной» характеристикой — продолжительностью времени работы стартовой ДУ), детерминированная оптимизация всех проектных параметров ракеты по критерию дальности L обеспечивает максимально возможны^ общий прирост дальности AZ, = 1713 км.
Раздельная оптимизация групп проектных параметров (приj = 1, 2, 3):
-продолжительностей работы ДУ т/. А1Т= +201 км (11,73 %);
-секундных расходов топлива ДУ mTJ АЬт= +512 км (29,89 %);
-давлений в ДУ pKj ALpK= +441 км (25,75 %);
-давлений на срезах сопл ДУ рау АЬрг = +334 км (19,50 %),
что в сумме дает *87 % от AL.
Раздельная оптимизация характеристик отдельных ДУ‘
-1-й (стартовой) ДУ: ALi = +385 км (22,48 %);
-2-й (средней) ДУ: AL2 = +223 км (13,02 %);
-3-Й (верхней) ДУ: AL3 = +778 км (45,42 %), что в сумме дает «81 % от AL.
Раздельная оптимизация приводит, как правило, К потерям максимальной дальности по сравнению с результатами комплексной оптимизации всех проектных параметров. Тем
не менее можно ориентироваться на тот факт, 4X0 в процессе оптимизации перераспределение топлива между ступенями 3-студенчатой баллистической ракеты обеспечивает около 41 % общего прироста дальности; оптимизация только давлений в ДУ - около 26 % ц оптимизация только давлений на срезах сопловых блоков Ду - около 20 %. С Другой стороны, оптимизация тех же параметр0®только для 3-й (верхней) ДУ обеспечивает около 46 %, 2-й (средней) ДУ - около I3 % и 1-й (стартовой) Ду - около 23 % общего прироста дальности.
Выводы справедливы для варианта использования конструкционных материалов обечаек, днищ и сопловых блоков ракеты в меТалличесК0М варианте (в рамках математической модели ракеты и РДТТ). Для иных, например композиционных, материалов возможно иное перераспределение приведенных величин.
6.10. Влияние удлинения РДТТ на результаты детерминированной оптимизации ракеты
В процессе параметрического проектирования [71] с заданием различных предельных значений удлинений цилиндрической части РДТТ Xjm&Xl j = 1, 2, 3, были получены оптимальные ракеты с минимальными стартовыми массами, доставляющие полезный груз массой тиг = 1700 кг на одну и ту же дальность L « 10 000 км. Результаты расчетов приведены в табл. 6.12:
Таблица 6.12 Результаты исследования влияния удлинения РДТ на оптимальные параметры ракеты
Характеристики |
|
Численные значения характеристик при |
с |
||
РДТТ и ракеты |
1 |
2 |
3 |
4 |
5 |
т0,т |
66,0 |
62,2 |
61,8 |
61,4 |
61,8 |
т т, т |
59,4 |
55,6 |
55,2 |
55,0 |
55,4 |
Я |
2,84 |
2,44 |
2,16 |
1,96 |
1,88 |
^4imi2 |
1,14 |
1,77 |
2,52 |
2,66 |
4,84 |
1о. М |
21,2 |
22,8 |
24,2 |
25,6 |
27,7 |
^иил1 |
1,14 |
1,78 |
1,99 |
2,43 |
2,28 |
T l . C |
53,0 |
60,7 |
50,0 |
54,5 |
50,0 |
m n , кг/с |
55,8 |
64,6 |
54,7 |
60,6 |
61,9 |
А ьатм |
29,7 |
37,2 |
37,2 |
40,5 |
37,2 |
Аьатм |
0,12 |
0,16 |
0,18 |
0,20 |
0,21 |
^иил2 |
М 4 |
1,77 |
2,52 |
2,66 |
4,84 |
Х 2 . С |
67,8 |
57,9 |
54,3 |
43,6 |
57,5 |
т Т 2 , кг/с |
304,1 |
214,2 |
192,7 |
206,9 |
212,8 |
p x2t атм |
23,1 |
40,5 |
42,0 |
48,5 |
44,5 |
p t2, атм |
0,21 |
0,25 |
0,32 |
0,39 |
0,52 |
^тшлЗ |
U 2 |
1,98 |
3,14 |
4,20 |
4,91 |
Т з . С |
40,0 |
47,2 |
51,9 |
52,7 |
50,5 |
т тъ, кг/с |
897,7 |
833,3 |
808,8 |
801,7 |
793,0 |
Ркз» атм |
38,4 |
48,2 |
60,3 |
69,6 |
63,8 |
р лз, атм |
0,70 |
0,63 |
0,66 |
0,72 |
0,84 |
Анализ полученных результатов показал в частноптн |
„ |
„ |
|
|
||||
1 s w S 5 |
увеличение |
сумм.™.,* |
у |
*Ж№ет* |
||||
растут стартовая масса ракеты т0 (до 7 %) |
|
с уменьшением Ядш |
||||||
диаметр Dmix (до 4 %); уменьшает™ длина ракеты нГ |
6 задасы Топлива тТ (до 8 |
%), |
||||||
полученные оптимальные |
варианты ракет д о с т ^ о ^ |
|
^® /о (до 33 %)• Поскольку все |
|||||
обеспечивая небольшие разбросы (до 4 %) |
|
|
|
На заданную дальность L |
||||
активного участка траектории, следует считало |
ы |
фазовы* |
характеристик кмпа |
|||||
указанных пределах из соображений п р ^ ^ ^ |
м |
подбор значений), |
в |
|||||
устойчивости ракеты, ее |
управляемости |
СКого |
порядка |
Значений |
W |
|||
характеристик проектируемой ракеты |
И Т'п0 |
без |
Z f |
(обеспечение |
запаса |
|||
|
Р ТЫ' |
|
|
|
существенного ухудшения |
|||
6.11. Влияние метода а шага ЧИслеп |
|
|
|
|
|
|
||
|
роцессе ее параМетричес^ х ар ак тер Иегик движеНия |
|||||||
Как правило, алгоритм численного интет |
|
|
оект«рования |
|
|
|||
моделирующих процесс движения ракеты н! ^ ИроВания Дифференпи». |
|
|
||||||
разложение искомых функций в пеги Тей™ * активном участи/ |
Циальных Уравнений, |
|||||||
разложения ряда о п р й Г е т м ^ д |
|
|
|
|
аСП° " » > " |
разложения ряда; метод Эйлера-Коши использу<^з^ц^ 1Ггод Э»*'Р« и ш о т я ^ Т ^ |
|
Купа 4-го порядка использует 5 член0, р а з д о р |
улож ения ряда; метод Р д а” |
же шага интегрирования точность результата воздетая Естественно, для одного и того
увеличением числа учитываемых членов ряда) |
с |
° услояшением алгоритма (т.е. с |
||||||
интегрирования может привести к удовлетворите |
|
ДРУГ°Й |
СТ0Р0НЫ> уменьшение шага |
|||||
точности для менее сложных методов. |
|
льным в ряде случаев результатам по |
||||||
В процессе, в частности, баллистического проектирования ракет с РДТТ |
возникает |
|||||||
необходимость |
выбора |
как метода, |
так |
и |
тяга |
Р |
Г’ возникает |
|
, , |
r |
м ’ |
ldK |
и |
шага |
интегрирования |
системы |
диффнреншальнык Уравнении движения ракеты на ее акиьиом >£ае™ траекторют (АУТ)- с применением компьютерное программы КАШ02 [45] были проведены расчеты
[75] для 3—ступенчатой ракеты с РДТТ со стартовой массой 50 072 кг, доставляющей полезный груз на дальность около 10 000 км. Полное время работы всех двигательных
установок составляет 177,84 с (74,04 с + 58,00 с + 45,80 с); |
конечный угол тангажа срк = |
= 20° Основные результаты (характеристики конца АУТ - |
скорость VK, координата хк, |
высота Нк , суммарные потери конечной скорости ДРке, дальность полета L) отражены в табл. 6.13:
Таблица 6.13 Результаты исследования метода и шага численного интегрирования
Метод |
Шаг |
|
Характеристики конца АУТ |
|
||
интегрирования |
интегрирования, с |
VK, м/с |
JCK, м |
Нк, м |
ДКкЛ, м/с |
L, км |
|
0,1 |
6808,75 |
353 328 |
181 089 |
-1472,69 |
9226,35 |
|
0,5 |
6739,66 |
351818 |
178 847 |
-1480,34 |
9007,52 |
Эйлера |
1,0 |
6621,84 |
347 447 |
175 443 |
-1489,95 |
8740,98 |
|
5,0 |
6244,09 |
315 446 |
160 406 |
-1560,94 |
7000,89 |
|
10,0 |
5909,91 |
288 931 |
147 439 |
-1621,47 |
5751,88 |
|
0,1 |
6832,65 |
356 197 |
181 776 |
-1470,79 |
9281,98 |
Эйлера - Коши |
0,5 |
6832,77 |
356 175 |
181791 |
-1470,99 |
9284,19 |
1,0 |
6831,10 |
356 102 |
181 786 |
-1471,35 |
9285,73 |
|
|
5,0 |
6843,51 |
355 154 |
181 651 |
-1477,48 |
9347,84 |
|
10,0 |
6875,56 |
354 268 |
171 497 |
-1486,48 |
9510,69 |
|
0,1 |
6832,54 |
356 197 |
181 773 |
-1470,77 |
9281,21 |
Рунге - Кутта |
0,5 |
6833,05 |
356211 |
181 787 |
-1470,82 |
9283,90 |
1,0 |
6832,99 |
356 197 |
181 788 |
-1470,87 |
9283,66 |
|
|
5,0 |
6832,86 |
356 110 |
181 835 |
-1470,98 |
9283,18 |
|
10,0 |
6832,92 |
356 113 |
181 794 |
-1471,11 |
9283,29 |
Анализ результатов расчетов показал, что в рамках изменения тага интегрирования в пределах от 0,1 с до 10 с (в 100 раз) наибольший разброс, как и следовало ожидать дает наименее точный метод - Эйлера (по дальности от -0,67 % при шаге 0,1с шаге 10 с). Учет всего лишь еше одного члена разложения для метода Эйлер^ошидает
' |
^лтшлг'ти - более, чем на 1 |
десятичный |
порядок (по дальности |
|||
о/Увелич^ни® |
45 % при шаге 10 с)- Наиболее точные результаты и261 MeTOn |
|||||
-41,02 % при шаге 0,1 с до |
прш« |
т |
ширенном |
wrajkpo.»»™ систем |
||
Рунге-Кугга, |
пшроко |
ШМММШн |
|
. умзаиныл |
||
дифферешшальных |
£ |
лиш ает ±0,007 %. Для росчто! т рат крт |
||||
пД1);Тпрк«деш .ы я расчета* - г |
1 с Для метода Руше-Ю™ яаляешя |
|||||
базовым. |
|
|
|
|
|
|
6.12. Связь оптимальных проектных |
|
сочаввип шердив- |
баллистической ракеты с конструкционными материалами
Рассмотрены все возможные варианты исполнения РДТТ баллистической ракеты,
доставляющей полезный груз массой 1700 кг на дальность ~10 000 км.
понимается изготовление корпуса или из металла (М) с характеристиками у = 0,008 кг/см и ов = 25 000 кгс/см2, или из стекловолокна (С) —у = 0,0017 кг/см3 и ов - 000 кгс/см . Исследования проводились с помощью компьютерной программы KAMFAD [48, 49, 50] для ракеты минимальной стартовой массы с оптимизацией проектных параметров. Результаты расчетов приведены в табл. 6.14, 6.15 и 6.16.
Втабл. 6.14 и 6.15 обозначены: I - «нижняя» ДУ (стартовая) и 1-я субракета (ракета на старте), II -«средняя» ДУ и 2-я субракета, III - «верхняя» ДУ и 3-я субракета.
Втаблице 6.16 обозначены: I - «нижняя» ДУ (стартовая) и 1-я субракета (ракета на старте), II - 2-я «верхняя» ДУ и 2-я субракета.
ДУразличных диаметров S-ступенчатойракеты [76]
Таблица 6.14 Результаты исследования вариантов ракет с различивши диаметрами ДУ
Характеристики ДУ |
|
|
|
Варианты ракет |
|
|
|
||
и суораке!г |
I |
c c c |
M C C |
с с м |
C M C |
M C M |
M M C |
CM M |
МММ |
Начальная |
47 634 |
50 736 |
50 935 |
50 337 |
63 168 |
55 682 |
54 632 |
61 085 |
|
масса, кг |
II |
1455 |
23 801 |
14 169 |
12 851 |
34 646 |
20 402 |
13 485 |
21 837 |
Запас |
III |
4699 |
4758 |
4630 |
6041 |
5253 |
7341 |
5750 |
5790 |
I |
30 798 |
24 564 |
33 721 |
35 182 |
25 484 |
32 282 |
38 906 |
36 229 |
|
топлива в |
II |
9381 |
17932 |
9642 |
6120 |
27898 |
12092 |
7014 |
14 941 |
ДУ, кг |
Ш |
2677 |
2715 |
2536 |
3939 |
3125 |
5176 |
3589 |
3639 |
«Сухая» |
I |
2036 |
2232 |
2198 |
2257 |
2876 |
2860 |
2193 |
2871 |
масса ДУ, кг |
П |
646 |
1076 |
667 |
611 |
1454 |
885 |
641 |
1016 |
|
ш |
263 |
285 |
306 |
344 |
339 |
406 |
370 |
360 |
Давление |
I |
109,7 |
70,2 |
109,7 |
105,9 |
90,1 |
71,8 |
93,4 |
63,5 |
в ДУ, кгс/см2 |
п |
76,5 |
70,2 |
76,5 |
51,4 |
62,8 |
45,3 |
47,2 |
41,7 |
|
ш |
53,3 |
58,6 |
37,2 |
58,0 |
36,5 |
60,0 |
39,4 |
36,9 |
Время |
I |
51,5 |
37,8 |
52,8 |
55,9 |
32,3 |
45,3 |
56,8 |
45,9 |
работы ДУ, с |
II |
55,5 |
73,5 |
56,0 |
42,4 |
81.5 |
56,5 |
46,3 |
63,1 |
|
111 |
48,3 |
53,5 |
48,8 |
63,5 |
59,5 |
69,3 |
57,9 |
55,6 |
Диаметр |
I |
172 |
176 |
180 |
180 |
192 |
184 |
184 |
192 |
ДУ, см |
II |
116 |
140 |
120 |
116 |
168 |
124 |
120 |
132 |
Длина |
111 |
100 |
100 |
100 |
100 |
100 |
100 |
100 |
100 |
I |
2544 |
2578 |
2568 |
2507 |
2685 |
2778 |
2576 |
2767 |
|
субракеты, |
II |
1461 |
1664 |
1446 |
1369 |
1813 |
1728 |
1374 |
1669 |
см |
III |
726 |
730 |
712 |
822 |
758 |
917 |
793 |
797 |
Анализ результатов расчетов для 3-ступенчатой ракеты показал, что с увеличением удельной прочности конструкционных материалов увеличиваются оптимальные значения рабочих давлений, уменьшаются запасы топлива, времена работы, «сухие» массы ДУ и стартовая масса ракеты с ДУ различных диаметров.
Начальные массы в чисто металлическом и в чисто композитном вариантах ракеты различаются на ~28 %. Влияние ограничений диаметров ДУ ракеты приводит к увеличению стартовой массы на ~11 %.
ДУ одного и того же диаметра 3-ступенчатой ракеты [77]
Таблица 6.15 Результаты исследования вариантов ракет сДУ одного и того же диаметра
Характеристики ДУ |
|
|
|
В а р и а н т ы |
р а к е т |
|
|
|
|
и субракет |
|
С С С |
М С С |
С С М |
С М С |
М С М |
М М С |
С М М |
МММ |
Начальная |
I |
51 628 |
52 846 |
61 792 |
55 181 |
59 726 |
57 178 |
62 430 |
66 184 |
масса, кг |
II |
24 898 |
37 272 |
20 157 |
21 578 |
39 369 |
34 328 |
20 077 |
37 566 |
|
III |
7821 |
7816 |
8899 |
8885 |
8622 |
10 605 |
8657 |
10 145 |
Запас |
I |
24 748 |
13 182 |
38 572 |
31 288 |
17 342 |
19 501 |
39 635 |
25 074 |
топлива |
II |
15 986 |
27 622 |
10 475 |
11442 |
28 922 |
22 130 |
10 305 |
25 411 |
в ДУ.кг |
III |
5605 |
5586 |
6470 |
6629 |
6263 |
8224 |
6239 |
7667 |
«Сухая» |
I |
1940 |
2272 |
3018 |
2272 |
2889 |
3223 |
2672 |
3411 |
масса ДУ, кг |
II |
1048 |
1792 |
737 |
1 128 |
1780 |
1468 |
986 |
1877 |
|
III |
446 |
459 |
612 |
485 |
535 |
609 |
588 |
643 |
Давление в |
I |
154,8 |
154,8 |
154,8 |
142,0 |
154,8 |
154,8 |
130,3 |
129,3 |
ДУ, |
II |
75,2 |
90,1 |
67,9 |
57,6 |
86,0 |
43,8 |
52,9 |
52,4 |
кгс/см2 |
1П |
36,5 |
36,5 |
21,3 |
33,5 |
21,3 |
36,5 |
25,7 |
25,5 |
Время |
I |
39,3 |
20,8 |
53,0 |
46,7 |
24,1 |
28,4 |
53,5 |
30,8 |
работы ДУ, с |
II |
65,6 |
78,7 |
46,8 |
50,7 |
76,5 |
66,7 |
47,4 |
67,9 |
|
III |
69,2 |
68,5 |
163,2 |
77,6 |
98,9 |
86,8 |
76,4 |
77,8 |
Диаметр ДУ, |
I |
148 |
148 |
156 |
152 |
156 |
156 |
154 |
160 |
см |
II |
148 |
148 |
156 |
152 |
156 |
156 |
154 |
160 |
|
III |
148 |
148 |
156 |
152 |
156 |
156 |
154 |
160 |
Длина |
I |
2431 |
2446 |
2565 |
2498 |
2534 |
2537 |
2596 |
2637 |
ракеты, см |
п |
1352 |
1767 |
1109 |
1222 |
1734 |
1654 |
1109 |
1632 |
|
III |
550 |
561 |
525 |
582 |
543 |
658 |
550 |
586 |
Анализ результатов расчетов для 3-ступенчатой ракеты показал, что с увеличением удельной прочности конструкционных материалов увеличиваются оптимальные значения рабочих давлений, уменьшаются запасы топлива, времена работы, «сухие» массы ДУ и стартовая масса ракеты с ДУ одного диаметра в условиях ограничения на удлинение ракеты.
Начальные массы в чисто металлическом и в чисто композитном вариантах ракеты различаются на ~28 %.
Влияние ограничений диаметров ДУ ракеты приводит к увеличению стартовой массы на -11 %.