УДК 629.78
ОЦЕНКА ВЕСОВОЙ ЭФФЕКТИВНОСТИ РАКЕТ-НОСИТЕЛЕЙ ПРИ СТАРТАХ С ДИРИЖАБЛЯ И САМОЛЁТОВ
© 2009 В. И. Куренков Самарский государственный аэрокосмический университет
Предложена методика сравнительного анализа весовой эффективности стартов ракет-носителей с поверхности Земли, дирижаблей, дозвуковых, сверхзвуковых и гиперзвуковых самолётов. Методика основана на оценке экономии потребной характеристической скорости ракет-носителей, стартующих с дирижабля и самолётов, по сравнению с наземным стартом и определении минимальной стартовой массы ракеты с помощью оптимального перераспределения массы между ракетными блоками. Методика может быть использована на начальных этапах проектирования.
Ракета-носитель, полезная нагрузка, стартовая масса, дирижабль, самолёт, характеристическая скорость, оптимизация массы, ракетные блоки.
Стартовая масса ракеты-носителя (РН) и масса её полезной нагрузки зависят от многих взаимовлияющих факторов: грузоподъёмности дирижабля или самолёта, высоты и скорости их полёта, угла наклона траектории ракеты в момент старта, количества ступеней ракеты и схемы соединения и включения ракетных блоков, используемых компонентов топлива и двигателей, программы изменения угла наклона траектории [ 1].
Поэтому минимизация начальной массы РН при заданной массе полезной нагрузки (ПН) или максимизация массы ПН при заданной стартовой массе РН является сложной задачей, так как связана с учётом взаимовлияния элементов технической системы «дирижабль - ракета-носитель» или «самолёт - ракета-носитель». Задача минимизации стартовой массы РН осложняется ещё и тем, что на начальных этапах проектирования многие характеристики, которые необходимо использовать в расчетах, окончательно не определены. Поэтому целесообразно иметь методику хотя бы приближённой минимизации стартовой массы РН, которая требовала бы небольшое количество исходных данных.
В дальнейшем в качестве критерия весовой эффективности рассматриваемых видов стартов будем использовать критерий начальной массы РН.
Общая постановка задачи. Необходимо определить минимальное значение стар-
товой массы РН т0 при заданной (фиксированной) массе выводимой на орбиту полезной нагрузки тПН и ограничениях по характеристической скорости Ух:
употр урасп £ 0*
тПН = соті,
где у"отр - ухрасп £ 0 - функция ограничения, выражающая условие непревышения потребной характеристической скорости уХ“отр над
располагаемой Ухрасп.
Данная постановка задачи эквивалентна следующей постановке: найти минимальное значение отношения стартовой массы РН
к массе полезной нагрузки р0 = т0/тПН при
тех же ограничениях:
Р0 =■
т
-» mln
т
ПН
употр ^ расп £ 0*
тПН = соті.
(1)
Суть методики заключается в определении потребной характеристической скорости ¥х для РН, стартующей с Земли, дирижабля или самолёта и расчёте массы т0 с
учётом оптимального распределения массы РН по ступеням.
Определение массы РН при наземном старте. Прежде всего следует отметить, что минимизация стартовой массы РН с использованием аналитических зависимостей возможна только для частных случаев, а именно для схемы с последовательным соединением ракетных блоков и одинаковыми удельными импульсами топлива и двигателей для всех ступеней РН [2].
В общем случае минимизация стартовой массы РН с учётом влияния множества факторов: различных схем соединения ракетных блоков (последовательное, пакетное, смешанное), конструктивных характеристик ракетных блоков, удельных импульсов, перелива компонентов топлива, соотношения сил тяги двигателей - возможна только численными методами. Методика минимизации стартовой массы РН и оптимизации массы ракетных блоков с учётом упомянутых факторов представлена в работе [3].
Для оценки начальной массы т0 при старте с поверхности Земли должны быть определены или заданы:
- масса полезной нагрузки mПН ;
- потребная характеристическая ско -рость ¥™тр;
- количество ступеней ракеты-носителя N
- схема соединения ракетных блоков
(РБ);
- компоненты топлива для каждого РБ;
- удельные импульсы топлива и двига-
теля Jуд і ( і = 1, N );
- конструктивные характеристики РБ
£ (і = 1, N ).
Конструктивная характеристика і-го РБ
есть отношение его массы тБ к массе блоБ і
симости:
тгпотр тгпотр , дтг
x x ид п
где У'пот - идеальная потребная характеристическая скорость; АУп - суммарные потери характеристической скорости от действия гравитационных, аэродинамических сил и сил противодавления на срезе сопла двигателя, действующих на РН в полёте.
Идеальная потребная характеристическая скорость для вывода полезной нагрузки
на опорную круговую орбиту радиусом гоо
вычисляется по следующей зависимости [2]:
потр
х ид
(
2Ю>
RЗ
л
(3)
oo \ З 0
где тЗ = 3,986 -105 км / с - гравитационная
постоянная Земли; RЗ = 6371,4 км - средний
радиус Земли.
Суммарные потери А¥п по статистике
составляют от 1350 до1650 м/с [3].
Рассчитаем минимальную стартовую массу РН.
Сначала на основе общей постановки задачи (1) должна быть сформулирована задача на поиск условного экстремума с конкретной целевой функцией и конкретной функцией ограничений.
Целевая функция p0, как показано в [3],
выглядит следующим образом:
ка без топлива mТ і :
где хі = m Т^т
ПН
(4)
(5)
m
£ =-
Б і
тБ- тТ і
Б і 1 і
(2)
Расчёт потребной характеристической скорости производится по следующей зави-
отношение массы топлива г-й ступени к массе полезной нагрузки.
Вид функции ограничений зависит от количества ступеней РН и схемы соединения РБ. В частности, для последовательного соединения РБ она имеет вид [3]:
УГ„ - Jу. |п
і=1
(* -1)
N
і=1
(£ -1)
. хі - х
Jуд2 ' |й
1+2
(£ -1)
1+2
(* -1)
... £ 0. (6)
Для других видов соединений РБ функции ограничений представлены в работе [3].
Далее решается задача на поиск условного экстремума функции (4) с учетом ограничений (6). В результате решения отыскиваются значения аргументов х°рр (г = 1, N) ,
которые соответствуют минимальному значению целевой функции (4).
Затем рассчитываются массы составных частей РН в последовательности:
- по (5) определяется масса топлива г-го РБ
тт і = тпн • х ;
(7)
- по (2) определяется масса і-го РБ
т =------------— т ■
ГПБ і / л\ Т і >
( £ -1)
(8)
РБ
тКі = тБ і - тті;
(9)
- определяется стартовая масса РН
N
т01 = тПН +2 тБ і .
(10)
её экономии в связи с воздушным стартом, которая возникает, во-первых, из-за старта РН с некоторой высоты и, во-вторых, из-за того, что дирижабль или самолет летят с определённой скоростью.
В свою очередь, экономия Ух из-за старта РН с некоторой высоты имеет две составляющие:
- экономия из-за изменения потенциальной энергии РН;
- экономия из-за уменьшения потерь от гравитационных, аэродинамических сил и сил противодавления на срезе сопла двигателя.
Экономия Ух из-за изменения потенциальной энергии РН, которую обозначим через АУ1, рассчитывается как разность идеальных потребных характеристических ско -ростей при наземном и воздушном стартах
Ґ
2ю>
к,
л
г
оо \ З 0
2 го
Л
оо \ З
К + Н
(11)
с 0
- определяется масса конструкции г-го
Полученная начальная масса РН будет минимальна, а распределение массы РН по ступеням - оптимальным.
Определение потребной характеристической скорости для РН при воздушном
старте. Проведём определение Ух с учётом
где Ис - высота старта РН.
В правой части этого выражения первый член представляет собой формулу (3), а второй член - модифицированную формулу
(3) с учётом высоты воздушного старта.
Расчёты показывают, что экономия ско -
рости А V мала. Например, для высоты старта 12 км она составляет всего лишь 14,7 м/с, а для высоты 17 км - 20,8 м/с. Поэтому в первом приближении её можно не учитывать.
Экономию потребной характеристической скорости Ух, возникающей из-за уменьшения потерь от гравитационных, аэродинамических сил и сил противодавления на срезе сопла двигателя, обозначим через АУ2. Эту
экономию можно определить, зная потребную характеристическую скорость РН для достижения высоты воздушного старта. По данным [4] для достижения высот 10.. .12 км и для различного количества ступеней и схем соединения ракетных блоков РН эта скорость составляет
ЛУ2 = 500...600 м / с . (12)
Экономию потребной характеристической скорости ЛУ2 можно также приближённо оценить, если учесть, что основные потери скорости РН связаны с действием сил гравитации, а угол наклона траектории меняется от 90° на небольших высотах до 40°.. .50° градусов на высотах 10.. .12 км. По статистике [4] для различных РН (исключая конверсионные) время достижения таких высот составляет 60.70 секунд (РН «Союз» - 64 с, РН «Зенит» - 63 с, РН «Энергия» -65 с). Среднее значение угла наклона траектории составляет 67°...70°. Следовательно, потери скорости от действия гравитации составят:
ЛУС = £0 ■ і■ зіп(65о..70о) »(520...610)м/с .
Более точно экономию потребной характеристической скорости от действия сил гравитации и других сил можно получить, интегрируя уравнения движения конкретной РН. Для проектных расчетов достаточно интегрирования упрощённых уравнений движения [3].
Определим экономию Ух в связи с полётом дирижабля или самолёта с определённой скоростью. Обозначим эту экономию
через ЛУ3 .
Рассмотрим пуск РН с дирижабля, который может двигаться при старте ракеты только горизонтально со скоростью примерно до 200 км/час или 55,5 м/с. Реальная экономия потребной характеристической скорости ЛУ3 будет несколько меньше скорости дирижабля, так РН должна сменить направление полета с горизонтального на наклонный.
Рассмотрим пуск РН с самолёта. В работе [1] отмечается, что перед пуском РН самолёт должен сделать горку с максимально возможным углом наклона траектории, который составляет 20°...30°. Время интенсивного разворота РН составляет примерно 20 секунд при изменении угла наклона траектории РН от стартового значения до опти-
мального программного 50°.. .60° [1]. Реальная экономия ЛУ3 будет также несколько
меньше скорости самолёта. Однако в процентном отношении её следует ожидать большей из-за начального угла наклона траектории по сравнению с горизонтальным пуском РН с дирижабля.
Для учёта влияния пуска РН на экономию характеристической скорости ЛУ3 введём эмпирический коэффициент кввс:
ЛУ3 = кввсУнос, (13)
где Унос - скорость дирижабля или самолета.
Суммарная экономия потребной характеристической скорости при воздушном старте определится по следующей зависимости:
ЛУ = ЛУ1 +ЛУ2 +ЛУ3. (14)
Определение потребной массы РН при воздушном старте. Будем считать, что
УХпотр для РН при воздушном старте определена с учётом высоты и скорости дирижабля или самолёта.
Последовательность определения массы РН при воздушном старте совпадает с последовательностью, используемой при наземном старте. Целевая функция имеет вид
(4), а при воздушном старте изменится лишь функция ограничений, в которой должна учитываться экономия характеристической ско -рости ЛУ:
употр -ЛУ - Урасп £ 0. (15)
Определение масс составных частей ракеты-носителя и стартовой массы ракеты производится по тем же зависимостям и в той же последовательности, что и в случае наземного старта.
Результаты расчётов. На рис. 1 представлены результаты расчётов стартовой массы РН с последовательным соединением РБ при наземном и воздушном стартах. Масса полезной нагрузки во всех случаях одинако -ва и составляла 2,5 т.
Рис. 1. Влияние вида старта и топлива на начальную массу РН при заданной массе полезной нагрузки
Рассмотрены две группы РН: трёхступенчатые с компонентами топлива «керосин
- жидкий кислород» и двухступенчатые с компонентами топлива «жидкий водород - жидкий кислород». Цифры над столбиками диаграммы означают стартовую массу ракет-носителей.
Расчёты и оптимизация проводились по вышеизложенной методике с использованием системы Mathcad. Составляющая экономии потребной характеристической скорости РН, возникающая при воздушных стартах вследствие снижения потерь от сил гравитации, аэродинамических сил и сил противодавления на срезе сопла двигателя, рассчитывалась с помощью специального программного обеспечения, построенного в среде Delphi на основе интегрирования упрощённых уравнений движения.
Предварительно масса полезной нагрузки варьировалась и была подобрана таким образом, чтобы начальная масса трехступенчатой РН при наземном старте составляла примерно 100 тонн, что облегчает сравнение начальных масс РН в процентах.
Скорость полёта самолётов составляла М = 0,185; 0,83; 2,0 ; 4,0 , где М - число Маха. Высота полёта для самолётов со ско -ростью М = 0,185; 0,83; 2,0 составляла 12 км, для самолёта со скоростью М = 4,0 принималась равной 17 км. Эмпирический коэффициент кввс принимался для дирижабля 0,8,
а для самолётов - 0,9. Конструктивные характеристики РБ для трехступенчатой РН первой, второй и третьей ступеней составляли соответственно: ^ = 12, s2 = 10, £3 = 8,
а для двухступенчатой РН - ^ = 10, s2 = 8 .
Анализ результатов показывает, что чем больше скорость и высота полёта дирижабля или самолёта, тем меньшая масса РН требуется для вывода фиксированной полезной нагрузки. Наиболее эффективным в весовом отношении является старт с гиперзвукового самолёта при М=4 с компонентами топлива «жидкий водород - жидкий кислород».
Однако следует отметить, что использование водорода в качестве горючего РН влечёт за собой увеличение ее габаритов по сравнению с РН, использующими в качестве
горючего керосин. Это увеличение происхо -дит несмотря на уменьшение стартовой массы РН. Это объясняется тем, что средняя плотность топлива с учётом соотношения компонентов у пары «керосин - жидкий кислород » выше, чем у пары «жидкий водород
- жидкий кислород», примерно в пять раз (1000 кг/м3 и 200 кг/м3 соответственно).
На рис. 2 представлены результаты расчётов объемов РН для исходных данных, использовавшихся при расчете масс.
Анализ результатов показывает, что потребный объём для реализации РН, использующих водород, значительно превышает объём РН, использующих керосин. Поскольку при использовании конкретных типов самолётов габаритные размеры РН могут быть критичными с учётом её расположения по отношению фюзеляжа (над, внутри, под), то при выборе компонентов топлива с учётом удельного импульса и средней плотности топлива необходимо пользовать комплексный критерий эффективности, предложенный в [3].
Кроме того, при проектировании необходимо учитывать, что самолёты должны
включать дополнительное оборудование для подпитки постоянно испаряющихся низко-кипящих компонентов топлива. Это оборудование сложнее и тяжелее в случае, когда в качестве горючего используется жидкий водород.
В заключение отметим, что на практике воздушные старты реализованы только по программе «звёздных войн», а именно, для перехвата низкоорбитальных космических аппаратов. Проблема практической реализации проектов воздушных стартов для «мирного космоса» связана с экономическим аспектом создания сложной технической системы, в которой появляются новые элементы, в том числе дирижабли или самолёты специального назначения. Для своего создания и эксплуатации они требуют не меньших ресурсов (финансовых, материальных и людских) и времени, чем создание РН с наземным стартом. Поэтому при принятии окончательного решения о реализации тех или иных проектов воздушного старта помимо оптимизации стартовой массы РН необходимо проводить подробные исследования по оценке экономической эффективности.
Рис. 2. Влияние вида старта и топлива на объем РН при заданной массе полезной нагрузки
References
1. Ballistic missiles and boosters/ O.M.Alifanov, A.N.Andreev, V.N.Gushyn et al -edited by O.M.Alifanov.-M.: Publishing house Drofa, 2004.-512 p.
2. Kurenkov, V.I. Selection of fundamental ratings and design configuration of boosters/ V.I.Kurenkov, L.P.Yumashev - edited by corresponding member of Russian Academy of sciences D.I.Kozlov. Samara: Publishing house of SSAU, 2005.-239 p.
3. Belokonov, V.M. Computer-aided evaluation of space vehicles characteristics/ V.M.Belokonov, V.A.Vyujanyn.-Kuybishev, KuAI, 1987.-68 p.
4. International reference guide to space launch systems. Third edition. Steven J. Isakowitz, Joseph P. Hopkins Jr., Joshua B. Hopkins. Corporate Sponsors Lockheed Martin Corporation The Boeing Company. Published and distributed by American Institute of Aeronautics and Astronautics (AIAA) 1801. Alexander Bell Drive, Suite 500 Reston, Virginia 20191-4344.
ASSESSING WEIGHT EFFICIENCY OF CARRIER ROCKETS LAUNCHED FROM DIRIGIBLES AND AIRPLANES
© 2009 V. I. Kurenkov
Samara State Aerospace University
A method of comparative analysis of weight efficiency of carrier rocket launches from the Earth’s surface, dirigibles, subsonic, supersonic and hypersonic airplanes is proposed. The method is based on assessing the savings in the required characteristic velocity of carrier rockets launched from dirigibles and airplanes as compared to ground launches and determining the minimal launching mass of the rocket by optimal mass redistribution among the rocket units. The method can be used at the initial stages of design.
Carrier rocket, payload, launching mass, dirigible, airplane, characteristic velocity, mass optimization, rocket
units.
Информация об авторе Куренков Владимир Иванович, доктор технических наук, профессор кафедры летательных аппаратов, Самарский государственный аэрокосмический университет. Область научных интересов: методы проектирования космических аппаратов и ракет-носителей. E-mail: [email protected].
Kurenkov, Vladimir Ivanovitch, doctor of technical sciences, professor, professor of aircraft design department, Samara State Aerospace University. Area of research: methods of designing space vehicles and carrier rockets. E-mail: [email protected].
71
Библиографический список
1. Баллистические ракеты и ракеты-носители: пособие для студентов вузов / О. М. Алифанов, А. Н. Андреев, В. Н. Гущин и др. / Под ред. О. М. Алифанова. - М.: Дрофа, 2004. - 512 с.
2. Белоконов, В. М. Расчет характеристик летательных аппартатов с применением ЭВМ: учеб. Пособие /В. М. Белоконов, В. А. Вьюжанин. - Куйбышев, КУАИ, 1987. - 68 с.
3. Куренков, В. И. Выбор основных проектных характеристик и конструктивного облика ракет-носителей: учеб. пособие / В. И Куренков, Л. П. Юмашев. / под ред. чл.-корр. РАН Д. И. Козлова. - Самара: Изд-во Самар. гос. аэрокосм. ун-та, 2005. - 239 с.
4. International reference guide to space launch systems. Third edition. Steven J. Isakowitz, Joseph P. Hopkins Jr., Joshua B. Hopkins. Corporate Sponsors Lockheed Martin Corporation The Boeing Company. Published and distributed by American Institute of Aeronautics and Astronautics (AIAA) 1801. Alexander Bell Drive, Suite 500 Reston, Virginia 20191-4344.