SELECTING THE OPTIMUM CONTROL OF A SPACECRAFT WITH THE LOW-THRUST ENGINE FULFILLING TRANSPORT TO THE INCONVERTIBLE POINT L1 OF EARTH-MOON SYSTEM

Full Text

ВВЕДЕНИЕ Луна - это естественный спутник Земли, исследование и освоение ресурсов которого может приносить не только новые научные знания, но и практическую выгоду. Одной из главных проблем освоения Луны является высокая стоимость доставки грузов на селеноцентрические орбиты. Применение перспективных электроракетных двигателей в качестве базы для маршевой двигательной установки на космических аппаратах, функционирующих в системе Земля-Луна, позволит повысить эффективность транспортных операций. Высокая скорость истечения реактивной струи, характерная для двигателей малой тяги [1], обеспечивает значительно меньший расход рабочего тела по сравнению с двигателями на химическом топливе. Это позволяет уменьшить массу рабочего тела и объем баков для его хранения а, следовательно, стартовую массу космического аппарата (КА) и затраты на осуществление перелетов. Некоторые технические особенности двигателей малой тяги и теоретические аспекты поиска оптимальных законов управления описаны в статьях [2-6]. Однако, низкий уровень тяги, создаваемый двигательными установками такого типа, приводит к значительным длительностям перелётов. Например, перелет КА SMART-1 Европейского космического агентства с низкой околоземной на низкую окололунную орбиту длился более 450 суток. При выборе типа двигательной установки для осуществления миссий, необходимо знать какова будет минимально возможная длительность манёвра, не превысит ли она заданные ограничения? Решение задачи об оптимальном по быстродействию перелёте позволит ответить на эти вопросы. В статье рассматривается методика оптимизации перелётов c заданной околоземной орбиты к точке либрации L1, в окрестности которой планируется размещение посещаемой международной станции. На движение космического аппарата в окрестности точек либрации гравитационные поля Земли и Луны оказывают практически равнозначное воздействие, поэтому законы программного управления необходимо получать в рамках ограниченной проблемы трех тел. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДВИЖЕНИЯ КА И ОПТИМАЛЬНОЕ ПО БЫСТРОДЕЙСТВИЮ УПРАВЛЕНИЕ Движение аппарата описывается в инерциальной полярной барицентрической системе координат (рис. 1). Основная плоскость совмещена со средней плоскостью движения Луны, центр координат находится в барицентре системы Земли-Луна, полярная ось направлена в сторону точки весеннего равноденствия. Считается, что КА движется под действием гравитационных полей Земли и Луны, которые вращаются вокруг общего барицентра по плоским круговым орбитам с постоянной угловой скоростью. Влияние нецентральности Земли и Луны, а также атмосферы Земли и гравитации других небесных тел не учитывается. Рис. 1. Инерциальная полярная барицентрическая система координат Вектор реактивного ускорения от тяги двигателей лежит в основной плоскости, его величина зависит только от массы КА. Не учитываются изменения величины тяги при изменении светотеневой обстановки на орбите, деградации энергоустановки, а также особенностей углового движения КА. Направление реактивного ускорения задается углом , который отсчитывается от барицентрического радиус-вектора против часовой стрелки. Положение Земли и Луны определяется постоянными по величине радиус-векторами и и углом , имеющим постоянную скорость изменения , равную средней угловой скорости движения Луны относительно барицентра. Положение центра масс КА определяется радиус-вектором и полярным углом φ. Все фазовые координаты системы приводятся к безразмерному виду. Для этого все линейные расстояния относятся к среднему расстоянию между центрами Земли и Луны км, скорости - к средней скорости Луны относительно барицентра, ускорения - к гравитационному ускорению на орбите Луны, время - к среднему периоду обращения Луны, угловые величины переводятся в радианы. В этом случае уравнения движения имеют следующий вид: (1) где - отношение массы Луны к суммарной массе системы Земля-Луна; - номинальное безразмерное ускорение от тяги двигателей; с - безразмерная скорость истечения рабочего тела; - функция включения-выключения двигателей; - относительный расход топлива. Заменим в (1) угол φ на угол Луна - барицентр - КА по формуле . Дифференциальное уравнение, описывающее изменение этого угла примет вид: . В этих обозначениях расстояния от КА до Земли и Луны равны соответственно: . (2) Выбор программного управления проводится с использованием в качестве критерия оптимальности минимального времени перелёта: . (3) Введем в рассмотрение вектор фазовых координат КА , который подчиняется граничным условиям, соответствующим цели перелёта и возможным ограничениям; вектор функций управления, где U - область возможных значений функций управлений и - вектор проектных параметров КА. Формально задача оптимизации описывается следующим образом: требуется определить вектор функций управления, удовлетворяющий граничным условиям, уравнениям движения (1-2) и доставляющий минимум критерию оптимальности (3): . (4) В соответствии с формализмом принципа максимума Понтрягина вводим вектор сопряженных переменных и формируем функцию Гамильтониана . Из условия максимума Гамильтониана определим оптимальное направление вектора ускорения и функции включения- выключения двигателей: , (5) , . (6) Система сопряжённых уравнений, согласно принципу максимума определяется соотношением и имеет вид: (7) Значение относительного расхода рабочего тела в момент старта равно нулю (рабочее тело еще не израсходовано). Для задачи на оптимальное быстродействие расход рабочего тела в конечный момент времени не фиксирован и согласно условию трансверсальности . Согласно уравнениям (5-7) , при всех , поэтому при всех , . Поэтому, согласно (6), при всех функция включения-выключения двигателей , то есть двигатель работает без выключений, кроме, может быть, конечной точки траектории, в которой возможно , что соответствует выключению двигателя в момент завершения движения. Кроме того, в силу монотонного убывания функции , удобно использовать эту сопряженную переменную для выхода из интегрирования. Остальные сопряженные переменные и фазовые координаты могут изменяться немонотонно, а минимальная длительность перелета для заданных проектных параметров заранее неизвестна. Таким образом, задача об оптимальном по быстродействию плоском движении в системе Земля-Луна сводится к следующей двухточечной четырехпараметрической краевой задаче: требуется найти такие начальные значения сопряженных переменных, чтобы на концах оптимальной траектории выполнялись начальное и конечное условия, . Для задачи достижений точки либрации L1 конечные значения вектора фазовых координат задаются следующим образом: , здесь - безмерная скорость точки либрации L1, - безразмерное расстояние между точкой либрации L1 и барицентром. Для решения поставленной краевой задачи использовалась модификация метода Ньютона, описанная в работе [5], и метод продолжения по параметру [6]. В качестве начального приближения для решения краевых задач использовались значения, вычисленные по аппроксимационным зависимостям, полученным в [5] в рамках задачи двух тел (без учета гравитации Луны) при не фиксированной суммарной угловой дальности перелета. Метод продолжения по параметру, подробно описанный в [6], был использован для получения требуемых решений в рамках ограниченной задачи трех тел. В качестве изменяемого параметра использовался параметр - отношение массы меньшего притягивающего центра к суммарной массе системы. РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ И ОПТИМИЗАЦИИ ДЛЯ ЗАДАЧИ ДОСТИЖЕНИЯ ТОЧКИ ЛИБРАЦИИ L1 Проектные параметры аппарата выбраны близкими к параметрам исследовательского КА Европейского космического агентства SMART-1: кг, Н, c = 15 км/с. Рассматривается перелет КА с круговой барицентрической орбиты радиусом км в точку либрации L1 системы Земля-Луна. Система уравнений (1-2, 5-7) интегрировалась методом Рунге Кутты четвертого порядка. Краевая задача решалась до достижения суммарной погрешности по краевым условиям 0,001%. При решении задачи использовалась следующая методика. 1. В качестве начальных приближений для сопряженных множителей использовались значения, полученные в результате решения задачи об оптимальном по быстродействию перелёте между круговыми, компланарными орбитами без учета гравитации Луны [5]. В этом случае, согласно условию трансверсальности для угловой дальности , что при в любой момент времени, согласно второму уравнению (7) приводит к и порядок краевой задачи снижается до трёх. 2. Используя метод перемещения по параметру от до , определяется изменённый вектор начальных значений сопряжённых множителей, соответствующее управление и траектория движения с учетом гравитации Луны с нефиксированной конечной угловой дальностью КА-Луна (достигается значение ). 3. Условие трансверсальности в конечный момент времени заменяется на условие достижения требуемой угловой дальности КА-Луна , осуществляется процедура метода продолжения по параметру от к и получаются оптимальное управление и соответствующая траектория движения при требуемых граничных условиях. На рис. 2 -5 приведены полученные с использованием этой методики оптимальная программа управления, изменение компонент скорости КА и соответствующая траектория движения. Рис. 2. Оптимальная программа управления углом направления вектора тяги Рис. 3. Изменение радиальной составляющей скорости Vr Рис. 4. Изменение трансверсальной составляющей скорости Vφ Длительность оптимального перелёта зависит от относительного углового положения между Луной и КА в момент старта, на рисунке 6 приведена эта зависимость. Полученные данные показывают, что наилучший результат для КА с такими проектными параметрами 41,95 сут достигается, когда начальный угол между КА и Луной равен 147,7 град. Рис. 6. Зависимость оптимального времени перелета от начального углового положения КА Найденное оптимальное относительное угловое положение КА и Луны позволяет определить наилучшую дату начала движения, которая в рамках рассматриваемой круговой ограниченной задачи трех тел будет повторяться один раз за лунный месяц. ЗАКЛЮЧЕНИЕ Предложенная авторами методика, основанная на формализме принципа максима Понтрягина и методе продолжения по параметру, продемонстрировала свою эффективность при решении задач об оптимальных по быстродействию перелётах КА с электроракетными двигателями к точке либрации L1 системы Земля-Луна для различных граничных условий и проектных параметров КА. Кроме того, данная методика легла в основу разработанного программного комплекса, позволяющего автоматизировать процесс выбора оптимального по быстродействию управления и соответствующей траектории движения, которые могут быть использованы для проектно-баллистического анализа транспортных операций КА с ЭРДУ в системе Земля-Луна. БЛАГОДАРНОСТИ Результаты, представленные в статье получены в рамках выполнения государственного задания Минобрнауки России 9.5453.2017/8.9. Рис. 5. Оптимальная траектория достижения точки либрации L1

About the authors

Yu Dong Hu

Samara National Research University named after academician S.P. Korolyov

Email: huyudong6688@163.com
Master

Olga Leonardovna Starinova

Samara National Research University named after academician S.P. Korolyov

Email: solleo@mail.ru
Doctor of Science, Professor at the Space Engineering Department

References

Supplementary files