ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ
«ВЫСШАЯ ШКОЛА ЭКОНОМИКИ»
Содержание
Введение
Точками либрации в ограниченной задаче трех тел, описывающей движение тела малой массы в гравитационном поле, создаваемом двумя массивными телами, вращающимися вокруг общего барицентра, называются точки равновесия тела с малой массой (космического аппарата) относительно более массивных тел (Солнца и Земли). Существуют 5 точек либрации. 3 из них называются коллинеарными и расположены на прямой, соединяющей массивные тела. 2 точки расположены на орбите одного из массивных тел (Земли). Точка L2 расположена на прямой, соединяющей Солнце и Землю, за Землей.
Рис.1. Точки либрации.
Гало-орбитами называется класс периодических орбит вокруг коллинеарных точек либрации. Они образуются при совпадении периодов обращения космического аппарата (КА) вокруг точки либрации в плоскости эклиптики и в плоскости, перпендикулярной плоскости эклиптики. Периоды обращения КА определяются амплитудами орбиты.
Гало-орбиты впервые были предложены Робертом Фаркуаром для программы Аполлон в его диссертации в 1968 г. [1]. Он предложил поместить КА на гало-орбиту вокруг точки L2 системы Земля-Луна для связи с аппаратом, находящимся на обратной стороне Луны. Хотя данная концепция не была реализована, впоследствии были запущены многие миссии, включающие в себя гало-орбиты. Наиболее известными из них являются:
· ISEE-3, ACE - запущены на гало-орбиты вокруг точки L1 системы Солнце-Земля;
· Herschel, Plank - запущены на гало-орбиты вокруг точки L2 системы Солнце-Земля.
Гало-орбиты обладают рядом преимуществ по сравнению с другими ограниченными орбитами вокруг точек либрации. В частности, КА, находящийся на гало-орбите вокруг точки L1, имеет возможность наблюдать пространство между Солнцем и Землей со стороны Солнца. При этом достаточно большая амплитуда орбиты позволит КА не терять связь с Землей. Аналогичное преимущество есть и у орбит вокруг точки L2 - при достаточно большой амплитуде аппарат не будет попадать в тень Земли, т.е. не будет возникать перебоев с энергообеспечением из-за неактивности солнечных батарей.
Движение вокруг точек либрации является неустойчивым, т.е. малое отклонение от начальных параметров приводит к тому, что аппарат сходит с гало-орбиты и впоследствии выходит из окрестности точки либрации. Для удержания аппарата на орбите используются т.н. корректирующие импульсы (импульсы коррекции). Основная задача исполнения этих импульсов - нивелирование заранее неизвестной возрастающей компоненты движения.
Комплекс корректирующих импульсов с учетом места их исполнения, периодичности и т.д. называется стратегией удержания КА на гало-орбите (station-keeping). Так же под стратегией понимается непосредственно методика расчета импульсов. Разработка стратегии удержания является одной из важнейших задач при планировании миссии.
Данная работа посвящена разработке и анализу стратегии удержания КА на гало-орбите вокруг точки L2 системы Солнце-Земля.
Целью данной работы является разработка и реализация стратегии удержания КА и моделирование движения КА на гало-орбите.
Задачами данной работы являются:
· Создание программного обеспечения (ПО) для расчета движения КА на гало-орбите;
· Моделирование движения КА на гало-орбитах с различными амплитудами с использованием разработанного ПО;
· Получение зависимости между энергетикой поддержания орбиты и направлением исполнения импульса для различных амплитуд орбит и анализ полученных зависимостей;
· Расчет направления неустойчивости и применение полученных результатов при имитационном моделировании движения КА на гало-орбите.
Результаты данной работы, в частности, разработанное ПО могут быть применены для моделирования как гало-орбит, так и других ограниченных орбит вокруг точки L2 системы Солнце-Земля.
1. Современные методики удержания космического аппарата на ограниченной орбите вокруг точки L2 системы Солнце-Земля
Стратегиям удержания КА на ограниченных орбитах (гало-орбитах, орбитах Лиссажу и прочих) посвящены многие статьи. В данном разделе приведены краткие сведения о существующих стратегиях удержания КА и методиках расчета величин корректирующих импульсов.
Первоначальным этапом планирования миссии к точке либрации является расчет т.н. номинальной орбиты КА. Как правило, для описания движения КА по номинальной орбите используются упрощенные модели. Так, в статьях [1]-[11] приводится ограниченная круговая задача трех тел. В статьях [1]-[3] приводятся уравнения движения КА в инерциальной системе координат; в статьях [1]-[7] также приводятся уравнения движения КА в системе координат, связанной с соответствующей точкой либрации. На начальном этапе исследования наибольший интерес вызывает линеаризованная система уравнений движения; такие уравнения приведены в статьях [1], [8]-[10]. Более подробно ограниченная задача трех тел и уравнения движения КА на гало-орбите будут рассмотрены в соответствующем разделе.
В [11] приводится классификация стратегий удержания КА на гало-орбите вокруг точки L1 системы Солнце-Земля. Стратегии коррекции разделяются на два класса: loose control strategy (техника свободного контроля), в которой варьируется только одна компонента импульса с целью приближения орбиты КА к номинальной и tight control technique (техника строгого контроля), подразумевающая варьирование двух и более компонент корректирующего импульса.
В качестве примера tight control technique в [11] приводится методика расчета импульсов коррекции для КА ISEE-3. На первом этапе обрабатывались данные о траектории КА за время, прошедшее с момента выполнения последнего маневра. Затем орбита аппарата численно интегрировалась на один оборот вперед (примерно 178 дней) в реальной модели сил и в 8 равномерно распределенных точках вычислялось расстояние между номинальной орбитой и полученной. На следующем этапе сумма полученных расстояний минимизировалась за счет варьирования компонент вектора импульса. Эта процедура производилась для различных дат, после чего для исполнения импульса выбиралась дата, которой соответствовали минимальные затраты топлива (хотя иногда, по техническим причинам, выбирался день с неоптимальной величиной импульса).
Всего в ходе миссии было выполнено 15 импульсов коррекции. Импульсы совершались в среднем один раз в 82 дня, а суммарный импульс составил 30 м/с. Согласно [11], суммарный импульс мог бы быть меньше, но т.к. ISEE-3 был первым КА, выведенным на гало-орбиту, ученые, управлявшие его движением, стремились сделать минимальное количество маневров с минимальным риском.
В качестве примера loose control strategy в [11] описана методика, с помощью которой рассчитывались корректирующие импульсы для КА SOHO. Использованная методика называется orbital energy balancing. Ее суть заключается в том, что если энергия аппарата слишком велика (т.е. был исполнен слишком большой импульс), он отклоняется от номинальной орбиты в одну сторону, а если слишком мала - в другую. В данном случае уравнения движения КА интегрировались до оси Солнце-Земля. В процессе управления КА в течение миссии было добавлено условие, чтобы аппарат двигался перпендикулярно данной оси при ее пересечении. Импульсы исполнялись только вдоль прямой Солнце-Земля (в отрицательном или положительном направлении). Единственным варьируемым параметром была величина импульса коррекции (с учетом знака). В результате описанная аппаратом траектория была не периодической гало-орбитой, а квазигало-орбитой. С момента выхода КА на орбиту в марте 1996 было исполнено 9 корректирующих импульсов, не считая экстренные маневры, выполненные в период с сентября 1998 по март 1999 (в тот момент аппарат сошел с гало-орбиты и экстренные маневры выполнялись для его возвращения). Суммарный импульс, без учета экстренных маневров, составил чуть более 5 м/с. Также стратегия удержания КА SOHO была описана в статье [12].
В [13] приведен аналитический подход к расчету величин корректирующих импульсов. В статье приведены 2 основные методики: target point approach (приближение целевой точкой) и Floquet mode approach (приближение модами Флоке). Несмотря на то что статья посвящена орбитам вокруг точки L2 системы Земля-Луна, данные методики применимы и для орбит вокруг точек либрации системы Солнце-Земля.
Методика target point approach, впервые описанная в [15], использует штрафную функцию, которая зависит от затрат топлива и от отклонения параметров КА от номинальных. Задача метода состоит в минимизации штрафной функции. Для этой методики вводится ряд ограничений. Во-первых, между импульсами должен быть перерыв не менее некоторого заданного времени. Во-вторых, отклонение от номинальной траектории также должно быть не меньше заданного. В-третьих, при моделировании процесса удержания КА, величины отклонений от номинальной траектории сравниваются между собой и если отклонение уменьшается, величина импульса не рассчитывается. Для успешного расчета величины импульса требуется выполнение всех трех условий.