Научный журнал
Международный журнал прикладных и фундаментальных исследований
ISSN 1996-3955
ИФ РИНЦ = 0,593

РЕШЕНИЕ ЗАДАЧИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ ОДНОСТУПЕНЧАТОЙ РАКЕТОЙ ГОРИЗОНТАЛЬНОГО СТАРТА ПРИ ВЫВОДЕ ЕЕ НА НИЗКУЮ ОРБИТУ ЗЕМЛИ

Мозжорина Т.Ю. 1 Рожков А.А. 1
1 ФГБОУ ВО «Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет)»
В данной работе решается задача оптимального управления одноступенчатой ракетой горизонтального старта при выводе ее на низкую орбиту Земли. Оптимизация основана на применении принципа максимума Понтрягина и численных методов. Предполагалась возможность изменения направления вектора тяги и наличия аэродинамических рулей, поддерживающих постоянство нулевого угла атаки. Ограничений на допустимый угол вектора тяги по отношению к оси ракеты, а также на допустимый уровень перегрузки не накладывалось. Ставилась задача исследования возможности вывода на орбиту одноступенчатой ракеты горизонтального старта с дозвукового носителя при оптимизации угла действия тяги двигателя. Минимизировался потребный расход топлива на выполнение задачи или время вывода на орбиту, что является эквивалентным при постоянном расходе топлива жидкостного реактивного двигателя. Численное решение с использованием принципа максимума Понтрягина было найдено методом пристрелки. В качестве исходных данных были выбраны характеристики первой ступени трехступенчатой американской ракеты легкого класса «Пегас». Численный эксперимент показал, что при оптимальном управлении вектором тяги и применении топлив с высоким удельным импульсом возможен вывод одноступенчатой ракеты на низкую орбиту Земли.
оптимальное управление
принцип максимума Понтрягина
метод пристрелки решения краевых задач
одноступенчатая ракета горизонтального старта
управляемый вектор тяги
1. Пегас (ракета-носитель). Материал из Википедии. [Электронный ресурс]. URL: https://ru.wikipedia.org/wiki/%D0%9F%D0%B5%D0%B3%D0%B0%D1%81_(%D1%80%D0%B0%D0%BA%D0%B5%D1%82%D0%B0-%D0%BD%D0%BE%D1%81%D0%B8%D1%82%D0%B5%D0%BB%D1%8C) (дата обращения: 07.11.22).
2. Пегов В.И., Киселев В.И. Аналитическое представление аэродинамических характеристик летательных аппаратов сложной формы // Наука ЮУрГУ: материалы 66-й научной конференции секции технических наук. 2014. С. 1728–1739.
3. ГОСТ 4401-81 Группа Т27 межгосударственный стандарт атмосфера стандартная. М.: Стандартгиз, 1949.
4. Федоренко Р.П. Приближенное решение задач оптимального управления. М.: Наука, 1978. 486 с.
5. Мозжорина Т.Ю. Численное решение задач оптимального управления с переключением методом пристрелки // Математическое моделирование и численные методы. 2017. № 2 (14). С. 94–106.
6. Мозжорина Т.Ю., Чуванова Л.О. Моделирование и оптимизация перелета спутников малой массы с земной орбиты на орбиту Марса с помощью ионных двигателей // Математическое моделирование и численные методы. 2021. № 2. С. 54–67. URL: https://mmcm.bmstu.ru/articles/250/ (дата обращения: 09.12.2022).
7. Синицын А.А. Расчет траектории межпланетного перелета Земля – Марс с малой тягой без использования метода грависфер. [Электронный ресурс]. URL: http://trudymai.ru/published.php?ID = 80987 (дата обращения: 29.05.2021).

Выведение спутников на низкую околоземную орбиту возможно как с применением вертикального, так и с горизонтального старта, когда ракета-носитель выводится на некоторую высоту с дозвуковой скоростью самолетом, после чего происходит ее сброс и включение двигателей первой ступени. Горизонтальный старт имеет ряд преимуществ, особенно в классе легких ракет, выводящих небольшие по массе грузы на невысокие орбиты Земли. Примером таких ракет может служить американская ракета «Пегас» [1]. Она имеет 3 ступени с двигателями, работающими на твердом топливе. Относительная масса топлива составляет для нее около 90 % от стартовой массы. Представляет интерес оценить возможность вывода на низкую околоземную орбиту одноступенчатой ракеты с горизонтального старта, которая будет оборудована двигателем с изменяемым вектором тяги. Проведенный численный эксперимент в данной работе представляет собой предварительную оценку такой возможности с рядом допущений и упрощений математической модели. Задача оптимизации управления вектором тяги была решена методом пристрелки с использованием принципа максимума Понтрягина. Расчет аэродинамического сопротивления проводился при допущении постоянства коэффициента аэродинамического сопротивления missing image file, где сx – коэффициент аэродинамического сопротивления, ρ – плотность воздуха, зависящая от высоты, V – скорость полета, S – характерная площадь летательного аппарата (в данной задаче площадь миделя ракеты). Это связано с тем, что достижение сверхзвуковых скоростей на высотах выше 12000 м происходит достаточно быстро, а основное изменение коэффициента аэродинамического коэффициента происходит на трансзвуковых скоростях [2]. Также предполагалось, что ракета имеет аэродинамические рули, позволяющие поддерживать нулевой угол атаки и не создающие значительной подъемной силы. Кроме этого, не рассматривались ограничения на возможный угол отклонения вектора тяги и прочностные ограничения конструкции ракеты. Плотность воздуха аппроксимировалась экспоненциальной зависимостью, дающей достаточно хорошее совпадение со значениями Международной Стандартной Атмосферы [3] на высотах от 13 км и выше.

Цель исследования – провести численный эксперимент для легкой ракеты горизонтального старта, идентичной по параметрам первой ступени ракеты «Пегас», решить задачу оптимизации управления (угол действия тяги двигателя), используя принцип максимума Понтрягина методом пристрелки, по полученным результатам оценить возможность вывода на низкую орбиту Земли одноступенчатого варианта ракеты.

Математические модели и численные методы

Критерий эффективности управления – минимум затрат топлива на выполнение миссии или минимум времени вывода ракеты на орбиту при постоянном уровне тяги, а следовательно, и расхода топлива.

missing image file.

Система уравнений движения в полярной системе координат:

missing image file

где R – расстояние от центра Земли, м; V – скорость, перпендикулярная радиусу, м/с; U – скорость вдоль радиуса, м/с; φ – угол, рад; P = JGT – тяга двигателя, Н; J – импульс ракетного двигателя, м/с; GT – расход топлива, кг/с; Xaer – сила аэродинамического сопротивления, Н; m– масса ракеты, кг; missing image file – ускорение свободного падения, м/с2, RЗ – радиус Земли, м; t – время, с, θ – управление, угол действия тяги (между осью ракеты и направлением силы тяги), рад, θtr – траекторный угол (между скоростью и перпендикуляром к местному радиусу), рад,

missing image file; missing image file;

ρ = 1,703746625(R – RЗ)exp(–0,000142422) – плотность воздуха по высоте, кг/м3, S / m0 – относительная площадь миделя ракеты.

В качестве исходных данных были приняты технические характеристики первой ступени ракеты Пегас. Момент сброса ракеты с самолета-носителя принимался при высоте 13000 м и горизонтальной скорости полета при М = 0,85. При использовании принципа максимума Понтрягина функция Понтрягина имеет вид

missing image file,

где ψR, ψV, ψU, ψφ, – сопряженные переменные.

Максимум функции Понтрягина по управлению определим из missing image file.

Отсюда для оптимального управления

missing image file.

Дополним систему уравнений движения сопряженной системой:

missing image file

missing image file

Получим П-систему с краевыми условиями:

missing image file

Решение было получено методом пристрелки. Метод пристрелки или стрельбы дает наиболее точные результаты численного решения краевых задач [4]. Возможность решения реальных задач оптимального управления указанным методом как для задач без ограничений по управлению, так и для задач с переключением управления показана в [5–7]. Внешний цикл – модифицированный метод Ньютона обращается к численному методу решения задачи Коши – методу Рунге – Кутты 4-го порядка. Для более точного выхода из метода численного интегрирования вводился новый аргумент missing image file, где Т – время вывода ракеты на орбиту соответственно, П-система претерпевает изменения, связанные с ведением нового аргумента.

Параметры пристрелки и невязки для этого варианта расчета выбирались следующие:

ψr(0), ψV(0), ψU(0), T, – параметры пристрелки

missing image file

missing image file

missing image file

Точность вычислений определялась условием: max{|δi|}<10–5 (i ∈ [1, 4]).

Результаты расчетов

Расчеты были проведены для различных значений удельного импульса топлива (2000 м/с и 4000 м/с) и для дозвуковой скорости старта и небольшой сверхзвуковой скорости старта ракеты (соответствующим на 13 км числу Маха полета М = 0,85 и М = 2).

Результаты расчетов приведены на рис. 1–6 и в таблице. Варианты 1 и 3 расчета соответствовали дозвуковой скорости старта и удельным импульсам топлива, равным соответственно 2000 м/с и 4000 м/с. Варианты 2 и 4 соответствовали сверхзвуковой скорости старта и удельным импульсам топлива, равным соответственно 2000 м/с и 4000 м/с.

Заключение

Увеличение начальной стартовой скорости с дозвуковой до сверхзвуковой при М = 2 оказывает незначительное влияние на время вывода ракеты на орбиту, а следовательно, почти не снижает потребные затраты топлива. Вместе с тем переход через трансзвуковые скорости значительно увеличивает нагрузку на начальный носитель ракеты и усложняет конструктивное решение.

Увеличение значения импульса ЖРД сокращает потребное количество топлива почти на 20 %, что делает реализацию одноступенчатой ракеты-носителя с горизонтальным стартом при оптимальном управлении вектором тяги вполне вероятной.

missing image file

Рис. 1. Изменение радиальной скорости ракеты по времени при выводе ее на низкую орбиту Земли

missing image file

Рис. 2. Изменение окружной скорости ракеты при выводе ее на низкую орбиту Земли

missing image file

Рис. 3. Изменение радиуса (расстояния от центра Земли) при выводе ракеты на орбиту

missing image file

Рис. 4. Изменение оптимального угла действия тяги при выводе ракеты на орбиту

missing image file

Рис. 5. Изменение траекторного угла при выводе ракеты на орбиту

missing image file

Рис. 6. Изменение величины относительного аэродинамического сопротивления Xaer / m (ускорения от силы аэродинамического сопротивления) при выводе ракеты на орбиту

Потребная для вывода на орбиту относительная масса топлива

Номер варианта

V0, м/с

J, м/c

Т, c

Относительная масса топлива, %

1

236,11

2000

105,819

95,24

2

600

2000

105,269

94,74

3

236,11

4000

86,5976

77,94

4

600

4000

85,4353

76,89

Следует отметить, что такая пара окислитель – горючее, как кислород и водород, дающие импульс, близкий к 4000 м/с, может не обеспечить теоретического выигрыша, так как водород обладает низкой плотностью, что приведет к увеличению размеров топливных баков и изменению габаритов ракеты. Это может сказаться на ухудшении аэродинамических показателей и в результате на увеличении времени вывода ракеты на орбиту.


Библиографическая ссылка

Мозжорина Т.Ю., Рожков А.А. РЕШЕНИЕ ЗАДАЧИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ ОДНОСТУПЕНЧАТОЙ РАКЕТОЙ ГОРИЗОНТАЛЬНОГО СТАРТА ПРИ ВЫВОДЕ ЕЕ НА НИЗКУЮ ОРБИТУ ЗЕМЛИ // Международный журнал прикладных и фундаментальных исследований. – 2022. – № 12. – С. 103-109;
URL: https://applied-research.ru/ru/article/view?id=13492 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674