О ВЛИЯНИИ ПРЕЦЕССИИ ОРБИТЫ ЛУНЫ НА ЭВОЛЮЦИЮ И ВРЕМЯ СУЩЕСТВОВАНИЯ ВЫСОКОАПОГЕЙНЫХ ОРБИТ ИСЗ Виктория И. ПРОХОРЕНКО Институт Космических Исследований РАН ИКИ РАН Семинар «Механика, Управление, Информатика, 24 февраля 2005

АННОТАЦИЯ Проблема выбора орбит ИСЗ с длительным временем существования не является новой. Она всегда была актуальной, однако в последние годы эта проблема приобрела большую актуальность из-за того, что существенно повысился ресурс бортовых приборов и бортовой энергетики. Теперь необходимое время существования орбит измеряется не годами, а десятилетиями. Это заставило произвести переоценку ценностей и вернуться к рассмотрению факторов, влияющих на длительность времени существования. 2

Введение П олученное в классической работе М.Л. Лидова [1961] решение для спутникового варианта двукратно- осредненной ограниченной задачи трех тел содержат необходимые, но не достаточные рекомендации для решения рассматриваемой задачи. Численные эксперименты в рамках ретроспективного анализа времени существования спутников серии «Прогноз», запущенных на высокоапогейные орбиты в , показали существенное влияние некоторых факторов, которые лежат вне решения осредненной задачи трех тел

Введение (2) Для орбиты «Прогноз-6» с начальным значением высоты апогея (перигея) км (498 км) и датой старта время существования составляло около 40 лет. Численный расчет эволюции под влиянием гравитационных возмущений от Луны и Солнца для гипотетического варианта орбиты «Прогноз-6» с измененной ровно на один год датой старта (1978) и фиксированными значениями всех остальных начальных условий дает время существования более 500 лет. А гипотетический перенос запуска на 1988 год привел бы ко времени баллистического существования не более 7 лет.

Введение (3) Здесь дело в прецессии орбиты Луны (которая, сохраняя постоянное наклонение к плоскости эклиптики около 5.5, прецессирует с периодом 18.6 года). Упомянутые орбиты типа ПРОГНОЗ-6 с разными датами старта отличаются начальными значениями углового расстояния между линиями узлов орбиты спутника и орбиты Луны на плоскости эклиптики. Это различие никак не влияет на процесс эволюции, описываемый решением осредненной задачи.

Введение (4) Это приводит к необходимости создания инструментов для анализа эволюции орбит с учетом новых факторов. Использование для этих исследований численных методов решения полной системы дифференциальных уравнений требует очень больших затрат машинного времени на той технике, которой мы располагаем на сегодняшний день. Расчет эволюции орбиты высокоапогейного спутника с учетом влияния Луны и Солнца на 10 лет требует около 50 мин машинного времени.

Введение (5) Мы предлагаем алгоритм, основанный на использовании эволюции по Лидову в качестве невозмущенного движения, а прецессию орбиты Луны – в качестве возмущающего фактора. В докладе обсуждаются различные аспекты разработки этого алгоритма и первые результаты его использования.

Полученные М.Л. Лидовым [1961] аналитические решения двукратно осредненной ограниченной задачи трех тел в хилловском приближении a - большая полуось, = 1 - e 2, e – эксцентриситет; i,, и - наклонение, аргумент перицентра и прямое восхождение восходящего узла орбиты ИСЗ, отнесенные к плоскости орбиты возмущающего тела; N – номер витка; 1 – параметр орбиты возмущающего тела; M, M 1 – масса центрального и возмущающего тел С1С1 С2С2 Область возможных значений интегральных констант с 1, с 2 8

Примечание: Здесь приведены выражения для параметров k,, для c 2 > 0 при c 2 < 0 в выражениях (5), (6) нужно поменять местами и Зависимость орбитальных элементов от времени 9 = 1-5/2c 2 ;, корни уравнения: Использование замены переменных T u = 2 K(k) - период изменения (u) ; T = 2 K(k) - период изменения ( ) В работе Ю.Ф. Гордеевой [1968] получено выражение зависимости эволюции орбитальных элементов от времени через эллиптический интеграл первого рода. приводит к выражению через эллиптический интеграл первого рода где где u = / (5)(5) (6)(6)

Безразмерный период либрационных составляющих эволюции орбитальных элементов Выразим период либрации T через безразмерный период T C, который в свою очередь выражается через безразмерные параметры L С (c 1, c 2 ) = 2 K(k); 10 K (k)- полный эллиптический интеграл первого рода. a * - большая полуось спутника; L D - параметр подобия возмущений; (7) (8)

Используемые варианты безразмерного времени и соответствующие им значения либрационного периода t* = t/ (либрационный период T C ) = t* /B (либрационный период Lc ) u = / (либрационный период 2K) = u/K(k)+ /2 (либрационный период 2 ) 11

Выражение для в функции параметра 12 Далее возможны два пути : Использование полученного М.А Вашковьяком [1999] решения, основанного на представлении через sn эллиптический, что приводит к выражению ( ) через эллиптические интегралы третьего рода Использование приближенной апроксимации с помощью sin( ), что позволяет выразить ( ) через элементарные функции

Применение гармонической апроксимации 13 Табличный интеграл I( ) Замена переменной = 2 / Lc + /2 (9) (10)

Представление в функции параметра в виде суммы ротационной и либрационной составляющих 14 (11) (12) (13) (14)

Период ротационной составляющей эволюции, при аппроксимации синусом 15 (15) (16)

На следующих слайдах представлены полученные двумя способами результаты расчетов эволюции параметров,, в функции на двух периодах либрации при различных значениях параметров с 1,с 2 : Сплошной линией показаны результаты расчетов с использованием обращения эллиптических интегралов первого рода [Ю.Ф. Гордеева, 1968], Штриховой линией – результаты расчетов, основанные на использовании аппроксимации = + sin. Расчеты эволюции параметра выполнены только одним способом (с применением аппроксимации) Линии уровня функции k 2 (c 1,c 2 ) показаны для значений k 2 = 0.1, 0.2, 0.4, 0.6, 0.8. Значения модуля эллиптического интеграла k 2 в функции с 1, с 2 16

c 1 = c 2min = c 2max =.3996 c 2 < 0:.1c 2min =-.06;.5 c 2min =-.28;.99 c 2min =-.55 c 2 > 0: -.1 c 2min =.06;.5 c 2max = c 2min =.28;.99 c 2max = 0.40

c 1 = c 2min = c 2max = c 2 < 0:.1c 2min =-.002;.5 c 2min =-.011;.99 c 2min =-.21 c 2 > 0: -.1 c 2min =.002; -.5 c 2min =.011;.5 c 2max =.180;.99 c 2max =.360

c 1 = c 2 < 0:.1c 2min =-.005;.5 c 2min =-.026;.9999 c 2min =-.051 c 2 > 0: -.1 c 2min =.005; -.5 c 2min =.026;.5 c 2max =.140;.99 c 2max =.280 c 2min = c 2max = 0.280

c 1 = c 2min = c 2max = c 2 < 0:.1c 2min =-.0005;.5 c 2min =-.0023;.9999 c 2min =-.005 c 2 > 0: -.1 c 2min =.0005; -.5 c 2min =.0023;.5 c 2max =.100;.99 c 2max =.200

c 1 = 0.6 c 2 > 0, c 2max =.16.1 c 2max =.016;.5 c 2max =.008;.9999 c 2max =.16 c 1 = 0.6, c 1 = 0.8 c 1 = 0.8 c 2 > 0, c 2max =.08.1 c 2max =.008;.5 c 2max =.04;.9999 c 2max =.08 21

Эволюция орбитальных элементов гипотетической версии орбиты ПРОГНОЗ-6 с датой старта = Время существования более 500 лет RpRp i 1, 1m 2, 2m 22 Численное интегриование Время расчета – 48 часов

Гипотетическая версия орбиты ПРОГНОЗ-6 с датой старта Эволюция параметров 1, 2, 1m, 2m на интервале времени 1978 –

Полуаналитический метод прогноза эволюции орбит под влиянием гравитационных возмущений внешнего тела с учетом прецессии орбиты возмущающего тела 24 В качестве основной плоскости отсчета будем использовать плоскость, перпендикулярную к оси прецессии орбиты возмущающего тела и проходящую через центральную точку. В нашей задаче возмущающее тела это Луна, а основная плоскость совпадает с плоскостью эклиптики. На следующе слайдах представлены тригонометрические формулы используемые для перехода от одной плоскости отсчета к другой.

Связь между угловыми элементами i,, отсчитанными относительно плоскости эклиптики, и соответствующими элементами i m, m, отсчитанными относительно плоскости орбиты Луны Зависимость значений i m, m от i, и углового расстояния между восходящими узлами орбит спутника и Луны на плоскости эклиптики определяется следующими соотношениями: cos i m = cos sin i sin + cos i cos ; i = i m - i; cos = (sin i cos - cos cos i sin ) / sin i m ; sin = - sin sin / sin i m ; m = + cos m = (sin i m cos i - cos cos i m sin i)/sin ; cos m = (-sin cos i + cos cos sin i) / sin i m ; sin m = sin i sin / sin i m ; = m – где = наклонение плоскости орбиты Луны к плоскости эклиптики, m - угловое расстояние линии узлов орбиты спутника на плоскости орбиты Луны от линии пересечения орбиты Луны с плоскостью эклиптики 25 (17) (18)

Зависимость i и от при разных значениях i i m = i + i(i, ) m = + (i, ) Использованы значения i от 20 до 80 с шагом 20. Красным цветом отмечена линия, соответствующая максимальному значению i 26 (19)

Описание алгоритма (1) Исходя из начальных условий a 0, 0, i 0, 0, 0, 0 = 0 - 0M, заданных в момент времени t 0 относительно плоскости эклиптики, рассчитываем значения 10, 20 по формулам, правые части которых совпадают с правыми частями первых интегралов c 1, c 2 (4) Лидовского решения осредненной задачи. Переходим к плоскости орбиты Луны, рассчитываем угловые элементы i 0m, 0m относительно этой плоскости по формулам (17) и значения 10m, 20m по формулам (4). 27

Описание алгоритма (2) Делаем шаг t, используя Лидовское решение (6 - 16) учитывающее в качестве орбита возмущающего тела орбиту Луны замороженную на момент времени t 0. Получаем значения 1, i 1m, 1m, 1m, соответствующие моменту времени t 1 = t 0 + t. Возвращаемся к плоскости эклиптики, преобразуем угловые элементы по формулам (17), (18), получаем i 1, 1, 1, 1 = 0M + 1. Первый шаг завершен. Для нового шага используем орбиту Луны, замороженную на момент времени t 1 28

Сопоставление результатов полуаналитического расчета эволюции орбитальных элементов и результатов численного интегрирования На следующих слайдах на примере орбиты типа ПРОГНОЗ-6 с гипотетическим стартом в 1978 году приводится сопоставление результатов полуаналитического расчета эволюции орбитальных элементов и результатов численного интегрирования на интервале времени 150 лет. И в том и в другом расчете учитываются гравитационные возмущений от Луны и Солнца и прецессия орбиты Луны Для сопоставления используются следующих параметры: три угловых элемента орбиты, i, ; угловое расстояния между линиями узлов орбиты спутника и орбиты Луны на плоскости эклиптики; параметры 1, 2, правые части которых соответствуют двум первым интегралам двукратно-осредненной задачи трех тел. Различие состоит параметре, характеризующем форму орбиты: в первом случае используется параметр = 1-e 2, а во втором случае – R p, радиус перицентра (измеряемый в радиусах земли). Эти параметры связаны между собой взаимно однозначным соответствием при известном значении большой полуоси a и имеют аналогичный характер эволюции. 29

Результаты полу- аналитического расчета с шагом 12 суток (а) и результаты численного интегрирования (б) на интервале времени 150 лет i 1 2 R p (R E ) (а) (б) 30 Время расчета а – 5 сек б – 14 часов

Эволюция положения точки 1 (t), 2 (t) в области возможных значений этих параметров на интервале времени 150 лет (а) б) Сопоставление результатов полуаналитического расчета эволюции орбитальных элементов с шагом 12 суток ( а) и результатов численного интегрирования (б). Начальная точка показана светлым кружком, а конечная точка – темным кружком. 31

i 1 2 б) (а) Сравнение результаты полу- аналитического расчета с шагом 12 суток (а) и с шагом 24 сутки (б) на интервале времени 150 лет 32

i 1 2 б) (а) Сравнение результато полу- аналитического расчета с шагом 12 суток (а) и с шагом 6 суток (б) на интервале времени 150 лет 33

Эволюция орбиты ИСЗ типа ПРОГНОЗ-6 в отсутствие прецессии орбиты Луны 34

Заключение Разработан эффективный инструмент исследования эволюции и времени существования орбит ИСЗ подверженных гравитационных возмущений от Луны и Солнца. 35 Полученный алгоритм отличается очень большим быстродействием. Расчет эволюции орбиты на 150 лет занимает не более 5 секунд против 14 часов при численном интегрировании. Это открывает широкие возможности для его применения как для выбора орбит, так и для исследования закономерностей эволюции различных классов орбит.

Список литературы Лидов М.Л. Эволюция орбит искусственных спутников планет под действием гравитационных возмущений внешних тел. // Искусственные спутники Земли С. 5. Моисеев Н.Д. О некоторых основных упрощенных схемах небесной механики, получаемых при помощи осреднения ограниченной круговой проблемы трех точек Труды ГАИШ, т.15, ч.1, с.100, Гордеева Ю.Ф. Зависимость элементов от времени в долгопериодических колебаниях в ограниченной задаче трех тел // Космич. исслед Т С Вашковьяк М.А. Об эволюции орбит далеких спутников Урана // Письма в "Астрон. журн." Т С Прохоренко В.И. Геометрическое исследование решений ограниченной круговой двукратно осредненной задачи трех тел // Космич. исслед Т С Прохоренко В.И. Исследование периодов эволюции эллиптических орбит в двукратно осредненной задаче Хилла // Космич. исслед Т С. 22. Назиров Р.Р., В.И. Прохоренко, А.И. Шейхет Ретроспективный геометрический анализ долгопериодической эволюции орбит и времени баллистического существования ИСЗ серии ПРОГНОЗ // Космич. исслед Т С Вашковьяк М.А. Тесленко Н.М. Построение периодически эволюционирующих орбит спутника сжатой планеты в осредненной задаче Хилла с учетом прецессии орбиты возмущающей точки // Письма в Астрон. журн Т , С