'

О ВЛИЯНИИ ПРЕЦЕССИИ ОРБИТЫ ЛУНЫ НА ЭВОЛЮЦИЮ И ВРЕМЯ СУЩЕСТВОВАНИЯ ВЫСОКОАПОГЕЙНЫХ ОРБИТ ИСЗ

Понравилась презентация – покажи это...





Слайд 0

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


Слайд 1

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


Слайд 2

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


Слайд 3

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


Слайд 4

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


Слайд 5

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


Слайд 6

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


Слайд 7

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


Слайд 8

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


Слайд 9

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


Слайд 10

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


Слайд 11

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


Слайд 12

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


Слайд 13

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


Слайд 14

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


Слайд 15

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


Слайд 16

? ? ? c1= .001 17 ? ? ? c2min = -.552 c2max= .3996 c2 < 0: .1c2min=-.06; .5 c2min=-.28; .99 c2min=-.55 c2 > 0: -.1 c2min =.06; .5 c2max = .20 -.5 c2min =.28; .99 c2max = 0.40


Слайд 17

? ? ? c1= 0.1 18 ? ? ? c2min = -0.210 c2max= 0.360 c2 < 0: .1c2min=-.002; .5 c2min=-.011; .99 c2min=-.21 c2 > 0: -.1 c2min =.002; -.5 c2min =.011; .5 c2max = .180; .99 c2max=.360


Слайд 18

? ? ? c1= 0.3 19 c2 < 0: .1c2min=-.005; .5 c2min=-.026; .9999 c2min=-.051 c2 > 0: -.1 c2min =.005; -.5 c2min =.026; .5 c2max = .140; .99 c2max=.280 ? ? ? c2min = -0.051 c2max= 0.280


Слайд 19

? ? ? c1= 0.5 20 ? ? ? c2min = -0.005 c2max= 0.200 c2 < 0: .1c2min=-.0005; .5 c2min=-.0023; .9999 c2min=-.005 c2 > 0: -.1 c2min =.0005; -.5 c2min =.0023; .5 c2max = .100; .99 c2max=.200


Слайд 20

c1= 0.6 c2 > 0, c2max= .16 .1 c2max=.016; .5 c2max=.008; .9999 c2max=.16 ? ? ? c1= 0.6, c1= 0.8 c1= 0.8 c2 > 0, c2max= .08 .1 c2max=.008; .5 c2max=.04; .9999 c2max=.08 ? ? ? 21


Слайд 21

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


Слайд 22

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


Слайд 23

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


Слайд 24

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


Слайд 25

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


Слайд 26

Описание алгоритма (1) Исходя из начальных условий a0, ?0 , i0, ?0, ?0, ?0 = ?0 - ?0M , заданных в момент времени t0 относительно плоскости эклиптики, рассчитываем значения ?10 , ?20 по формулам, правые части которых совпадают с правыми частями первых интегралов c1 , c2 (4) Лидовского решения осредненной задачи. Переходим к плоскости орбиты Луны, рассчитываем угловые элементы i0m, ?0m относительно этой плоскости по формулам (17) и значения ?10m , ?20m по формулам (4) . 27


Слайд 27

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


Слайд 28

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


Слайд 29

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


Слайд 30

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


Слайд 31

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


Слайд 32

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


Слайд 33

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


Слайд 34

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


Слайд 35

Список литературы Лидов М.Л. Эволюция орбит искусственных спутников планет под действием гравитационных возмущений внешних тел. // Искусственные спутники Земли. 1961. №. 8. С. 5. Моисеев Н.Д. О некоторых основных упрощенных схемах небесной механики, получаемых при помощи осреднения ограниченной круговой проблемы трех точек Труды ГАИШ, т.15, ч.1, с.100, 1945.  Гордеева Ю.Ф. Зависимость элементов от времени в долгопериодических колебаниях в ограниченной задаче трех тел // Космич. исслед. 1968. Т. 6. № 4. С. 536. Вашковьяк М.А. Об эволюции орбит далеких спутников Урана // Письма в "Астрон. журн." 1999. Т. 25. № 7. С. 554. Прохоренко В.И. Геометрическое исследование решений ограниченной круговой двукратно осредненной задачи трех тел // Космич. исслед. 2001. Т. 39. № 6. С. 622. Прохоренко В.И. Исследование периодов эволюции эллиптических орбит в двукратно осредненной задаче Хилла // Космич. исслед. 2002. Т. 40. № 1. С. 22. Назиров Р.Р., В.И. Прохоренко, А.И. Шейхет Ретроспективный геометрический анализ долгопериодической эволюции орбит и времени баллистического существования ИСЗ серии ПРОГНОЗ // Космич. исслед. 2002. Т. 40. № 5. С. 538. Вашковьяк М.А. Тесленко Н.М. Построение периодически эволюционирующих орбит спутника сжатой планеты в осредненной задаче Хилла с учетом прецессии орбиты возмущающей точки // Письма в Астрон. журн. 1998. Т. 24. № 6, С. 474. 36


×

HTML:





Ссылка: