Презентация на тему «Аннотация на части тела»

Об эволюции орбит исз под влиянием гравитационных возмущений от луны и

солнца и проблеме выбора долгоживущих восоко апогейных орбит

Памяти Павла Ефимовича ЭЛЬЯСБЕРГА

Институт Космических Исследований Российской Академии Наук

Виктория И. ПРОХОРЕНКО

ИКИ РАН Семинар «Механика, Управление, Информатика, 25 марта 2004

Аннотация (1из 2)

Речь идет о практической задаче выбора долгоживущих орбит ИСЗ с большим эксцентриситетом и наклонением. Орбиты ИСЗ серии ПРОГНОЗ, запущенные с 1972 по 1995 г.г. послужили экспериментальным материалом для исследований. На первой стадии исследований были использованы аналитические решения двукратно-осредненной ограниченной круговой проблемы Хилла, полученные М.Л. Лидовым [1961]. Геометрическая интерпретация этих решений позволила разработать геометрический метод анализа долгопериодической эволюции, и времени существования орбит ИСЗ. Предположение о компланарности орбиты Луны и плоскости эклиптики позволило применить упомянутые решения задачи трех тел к задаче четырех тел (Земля, Спутник, Луна Солнце).

2

Аннотация (2 of 2)

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

3

П.Е. ЭЛЬЯСБЕРГ и М.Л. ЛИДОВ

Введение

Эволюция эллиптической орбиты точки P (спутник нулевой массы) рассматривается в рамках ограниченной круговой проблемы трех тел. Точка P движется в поле притяжения центральной точки S (массы M) под влиянием гравитационных возмущений со стороны третьей точки J (массы M1), которая движется вокруг точки S по круговой орбите радиуса a1. М.Л. Лидов [1961] получил аналитическое решение двукратно-осредненной системы дифференциальных уравнений движения точки P в приближении Хилла, полагая что отношение большой полуоси a орбиты точки P удовлетворяет соотношению: ? = a/a1 << 1. Это позволило использовать первый член разложения возмущающей функции по параметру ?. Полученное аналитическое решение включает три первых интеграла и две независимых квадратуры.

5

С1

С2

Полученные М.Л. Лидовым [1961] аналитические решения двукратно осредненной ограниченной задачи трех тел в хилловском приближении

Критическое значение ?* , соответствующее соударению спутника с центральным телом радиуса R: ?* = 1- (1-R/a)2

Область возможных значений интегральных констант с1, с2

a - большая полуось, ? = 1 - e2, e – эксцентриситет; i, ?, и ? - наклонение, аргумент перицентра и прямое восхождение восходящего узла орбиты ИСЗ, отнесенные к плоскости орбиты возмущающего тела; N – номер витка; ?1 – параметр ? орбиты возмущающего тела; M, M2 – масса центрального и возмущающего тел

I. Основные закономерности эволюции высоты перицентра, использованные

в процессе проектирования орбит серии «ПРОГНОЗ»

Короткопериодическая эволюция высоты перицентра за виток

Изменение высоты перицентра орбиты спутника за виток ?hp зависит от значений большой полуоси спутника, эксцентриситета, и положения вектора возмущающего ускорения относительно орбитальной системы координат O???*)

Знак изменения высоты перицентра за виток ?hp зависит от угла ? между осью ? и проекцией вектора возмущающего ускорения на плоскости O?? : 0 < ?hp при ? ? I или III четверти ?hp < 0 при ? ? II или IV четверти

В книге П.Е. Эльясберга [1965] приведены оценки модуля максимального отклонения высоты перицентра за виток ??hp?max под влиянием гравитационных возмущений от Луны и Солнца для орбит с высотой апогея (перигея) от 2 000 до 100 000 км (от 200 до 50 000 км).

*) Правая система координат O???: начало координат совпадает с притягивающим центром, плоскость O?? совпадает с плоскостью орбиты спутника, ось ? направлена в точку перицентра, ось ? - по нормали к плоскости орбиты.

Долгопериодическая эволюция высоты перицентра

Знак долгопериодического изменения высоты перицентра зависит от значения аргумента перицентра ?, измеренного относительно линии узлов орбиты спутника на плоскости орбиты возмущающего тела: 0 < при ? ? II или IV четверти < 0 при ? ? I или III четверти

В книге П.Е. [1965] показано, что ? ?max = ? ??hp?max

P1

P8

P6

P2

P10

P3

I-1

Эволюция радиуса перицентра rp и время существования орбит ИСЗ серии «ПРОГНОЗ» (1972 –1995)

Численное интегрирование полной системы дифференциальных уравнений выполнено с учетом гравитационных возмущений от Луны и Солнца

Типичные начальные значения орбитальных элементов: 16.12 < a (RE) < 16.74; 0.930 < e0 < 0.936; ie0 = 65?(до P10); ?e0 = 290?(до P7). Угловые элементы измерены относительно плоскости земного экватора

P4,5,7

II

Геометрическое исследование первых интегралов задачи Хилла в сферической системе координат O?i?:

? = 1 - e2 (0 ? ? ? 1) - радиус; i (0? ? i ? 180?) - коширота ?; ? (0? ? i ? 360?) - широта ?

Сечения поверхностей вращения c1 =

cos2i плоскостями ? = 0?, 180?(а) и ? = 90?, 270?(b). Серым тоном здесь и далее выделена область, соответствующая значениям c2 < 0

Геометрическая интерпретация первых интегралов c1, c2

А

b

c

d

Интегральные кривые, соответствующие линиям пересечения поверхностей c2 = (1- ?) (2/5- sin2 ? sin2i) с поверхностями c1 = 0 (i = 90?) (c) и c1 = 0.2 (d)

12

? = 0?

? = 270?

? = 90?

? = 180?

Геометрическая интерпретация соударения спутника с центральным телом

конечного радиуса R

13

?* = (2a* - 1)/a*2, a* = a/R

a*

a* = 8, ?* = 0.234, c1=0.1, c2=0.1; c2= -0.1

Косой штриховкой показана область, соответствующая орбитам с конечным временем баллистического существования для a* = 8, определяемая неравенствами: c1 < 0.6 ?*2 или c2 > (1 - ?*)(c1 / ?* - 0.6) Гордеева [1968]

Соотношение между областями возможных значений начальных орбитальных

элементов ?0, i0, ?0 и интегральных констант c1, c2

c1 = ?0 cos2i0, c2 = (1 - ?0) (2/5 - sin2i0 sin2?0)

a

c

d

0 ? ?0 ? 1 superposition

?0 = 0.6

b

?0 = 0.4

14

Сферическая поверхность ?0=0.4

III

Параметрический анализ периодов долговременной эволюции элементов ?, i и мажоранты времени баллистического существования

16

Зависимость эволюции орбитальных элементов от времени

Время эволюции орбитальных элементов можно представить в виде произведения независимых параметров, используя квадратуру (2) и теорию подобия и размерностей

Ю.Ф. Гордеева [1968] выразила квадратуру L через эллиптический интеграл первого рода

T* = 4/15 a*-3/2

LC (c1, c2)?/LD,

?LC (c1, c2)? = 2L(c1, c2, ?min, ?max, ?/2)

Конфигурационный параметр подобия орбит LC (c1, c2) зависит только от c1, c2, его знак совпадает со знаком параметра c2, а абсолютное значение равно удвоенной квадратуре L, вычисленной в пределах ?min, ?max

9

8

9

8

6

6

17

Период T* долгопериодической эволюции орбитальных элементов (?, i) и безразмерный конфигурационный параметр подобия орбит LC (c1, c2)

Изолинии поверхности ?LC (c1, c2)?показаны для уровней от 6 до 13 с единичным шагом

?Lc(c1,c2)? сечение плоскостями c1 (c1<0.6)

TB* = 4/15 a*-3/2

LB (c1, c2, a*)?/LD,

18

Мажоранта TB* времени баллистического существования и безразмерный конфигурационный параметр подобия орбит LB (c1, c2, a*)

?LB (c1, c2 , a*)?= 2L(c1, c2, ?*, ?max, ?0(?*))

LB определено только для c1, c2 , при которых ?min< ?*< ?max

Изолинии для поверхностей ?LC (c1, c2)? и ?LB (c1, c2 , a*)? при a* = 16

Конфигурационный параметр LB (c1, c2 , a*) имеет тот же знак, что и c2, а абсолютное значение, равное удвоенной квадратуре L вычисленной в пределах ?*, ?max, с начальным значением ?0, определенным как функция от ?* (при sin 2?0(?*) < 0)

Линии соответствуют значениям уровня от 5 to 13 с единичным шагом

Острый пик при c2 = 0 (при c1 < 0.6) Зеркальная квазисимметрия

относительно плоскости c2 = 0 в окрестности c2 = 0 (при c1<0.6) ?LB(c1, c2, a*)? < ?LC(c1, c2)?при любых a* Выражаются через эллиптические интегралы первого рода [Гордеева, 1968]

Свойства функций ?LC(c1,c2)?и?LB(c1,c2,a*)?

19

IV

Анализ семейства орбит ИСЗ серии ПРОГНОЗ (a* = 16.6, ?* = 0.117) и метод выбора долгоживущих орбит

Для каждой орбиты значения параметров c1, c2 показаны черными точками и маркированы номером ИСЗ

Геометрический метод выбора долгоживущих орбит

Большая полуось a = 8 RE ?* = 0.234

Lc(c1,c2)

Высота перигея hp0 = 5000 km e0 = 0.777 ?0 = 0.4

LB(c1,c2, a*)

21

Область значений с1, с2, соответствующих орбитам с конечным временем существования

V. Сопоставление аналитических решений с результатами численного

интегрирования полной системы дифференциальных уравнений с учетом реальных возмущений от Луны и Солнца

· l =RE = 6371200 m,

=365 сут.; · ? = 0.39860044 1015 m3/s2 (Земля); · ?1 = 0.4902799 1013 m3/s2, a1 = 0.3844109 m, ?1 = 1 (Луна); · ?2 = 0.13271244 1021 m3/s2, a2 = 0.1495979 1012 m, ?2 = 1 (Солнце).

LD

0.00219

0.00101

0.00320

Безразмерный параметр подобия возмущений LD для системы тел: Земля, Спутник, Луна, Солнце

Система тел

Земля–ИСЗ–Луна

Земля– ИСЗ– Солнце

Земля– ИСЗ – Луна +Солнце

Значение LD в третьей колонке представляет собой сумму значений, расположенных в колонках 1 и 2

23

Использованы следующие характерные размер l, время ?, и динамические параметры центрального и возмущающих тел:

1

2

3

Луна

Солнце

Луна+ Солнце

7.60

18.0

5.20

7.72

16.8

5.29

Солнце

Сопоставление времени существования ИСЗ ИНТЕРБОЛ-1 под влиянием возмущений от Луны и Солнца вместе и отдельно, рассчитанного по аналитическим формулам и по результатам численного интегрирования с учетом реальных возмущений

Возмущающее тело

Учет реальных возмущений

Аналитическое решение

Эволюция радиуса перигея rp под влиянием Луны и Солнца отдельно и вместе

Луна

1995

2013

2000

rp = 6 RE

Луна + Солнце

Численный расчет (с учетом гравитационных возмущений от Луны и Солнца)

времени баллистического существования для гипотетических версий орбит v1-v7 типа ИНТЕРБОЛ –1 с различными значениями аргумента перицентра 314 ?? ?0e ?290?, при фиксированных значениях остальных орбитальных элементов (a =16.12 RE , e0= 0.93, i0e = 62.9?, ?0e = 260?, ??0 = 24.5?) и датой старта 03.08.1995

?0e = 298.15?

298.75?

300?

293?

292?

304?

290?

314.?

Ив-1

26

Сопоставление численных расчетов времени баллистического существования TBR с аналитическим расчетом мажоранты TB*(с1, c2, a*)

ИНТЕРБОЛ-1: a* = 16.12, c1 = 0.0179, с2 = 0.247, e0 = 0.93, ?0 = 0.123, ?0 = 338.7? и версии v1-v7 со значениям 328? ? ?0 ? 314? (0.14 ?c2?-0.036)

T*,TB*,TBR

Это позволило обнаружить «сдвиг» функции TBR относительно функции TB*

Сплошная (штриховая) линия показывает период эволюции T* (мажоранту времени баллистического существования TB*) в функции параметра c2. Расчетное время баллистического существования TBR показано в виде дискретных символов в функция значения параметра c2, определяемого начальными значениями орбитальных элементов. Светлые (темные) значки показывают расчетное время баллистического существования TBR, связанное с ротационным (либрационным) типом эволюции аргумента перицентра.

VI

Исследование влияния прецессии орбиты Луны на эволюцию орбитальных элементов ИСЗ и время их существования (Учитывается наклонение 5.15? плоскости орбиты Луны к плоскости эклиптики и прецессия орбиты Луны с периодом 18.6-года)

Вспомогательные функции

1(t), ?2(t) и ?1m(t), ?2m(t) для исследования эффекта от прецессии орбиты Луны

Для сопоставления аналитических решений с результатами численного интегрирования полной системы уравнений будем в процессе интегрирования следить за поведением функции ?1(t) и ?2(t) с начальными значениями ?1(t0) = c1 и ?2(t0) = c2 ?1(t) = ? cos2i ; ?2(t) = (1 - ?)(2/5 - sin2? sin2i). Параллельно рассмотрим другую пару функций ?1m(t), ?2m (t): ?1m (t) = ? cos2im; ?2m (t) = (1 - ?)(2/5 - sin2?m sin2im), где индекс m маркирует орбитальные элементы, измеренные относительно плоскости орбиты Луны. Из определения этих пар функций следует, что области их возможных значений совпадают с областью допустимых значений параметров c1, c2.

28

Параметр

, отвечающий за сдвиг функции TBR относительно функции TB*

29

Эволюция функций ?1m(t), ?2m (t) определяется эволюцией углового расстояния ? между восходящими узлами орбит спутника и Луны на плоскости эклиптики Для орбит с фиксированным начальным значением прямого восхождения восходящего узла ?0 начальное значение ?(t0) = ?0 зависит от даты старта, которая в свою очередь определяет позицию восходящего узла орбиты Луны. Угловая скорость эволюции параметра ? определяется как разность между угловой скоростью эволюции прямого восхождения восходящего узла орбиты спутника и постоянной угловой скоростью прецессии орбиты Луны . Эволюция параметра ? в рамках двукратно осредненной проблемы Хилла определяется квадратурой (3). М.А. Вашковьяк [1999] выразил эту квадратуру через эллиптические интегралы первого и третьего рода.

?0=40?

?0=78?

Зависимость времени баллистического существования и поведения функций ?1m(t), ?2m (t) от начального значения параметра ?0

Рассмотрены два варианта орбит с одинаковым значением c1 = 0.018 : IB1 - эквивалентен орбите ИСЗ ИНТЕРБОЛ-1 (?0 = 339?, c2 = 0.247) v4 - отличается от первого только начальным значением аргумента перицентра (?0 = 322.6?, c2 = 0.069).

T*, TB*,TR*, годы

C2 > 0

Для каждой из орбит сделан расчет времени баллистического существования TBR для набора дат старта, обеспечивающего покрытие всего интервала возможных значений параметра ?0 (0? ?0 ? 360?). Для каждой орбиты значения TBR(?0 ) отнесены к своему значению c2 и маркированы значениями ?0.

?0= -89?

Светлые (темные) значки соответствуют ротационному (либрационному) типу эволюции аргумента перицентра.

?1, ?1m

v4

v4

?2, ?2m

VII

ПРОВЕРКА ГИПОТЕЗЫ о роли параметра ?0 на примере других орбит серии ПРОГНОЗ

Цилиндрическая система координат O

??1

? = 1 - e2 (0 ? ? ? 1) -радиус; ? (0? ? ? ? 360?) - долгота; ?1 (0 ? ?1 ? 1) – координата Z

32

?1= 0

?1= 0.2

? = 270?

? = 90?

33

Эффект от начального значения параметра ?0 на примере ИСЗ ПРОГНОЗ-2

ПРОГНОЗ-2 (a*=16.7, ?*=0.116) ?0 = 284?, ?0 = 0.126 ?1(t0)=c1=0.07, ?1(t0)=c2 =-0.03

Характер эволюции параметров ?, ? и ?1, ?2 в зависимости от даты старта (определяющей значение параметра ?0)

Гипотетический запуск 29.VI.1981, ?0 = 247? Время существования ~ 60 лет

Реальный запуск 29.VI.1972, ?0 = 70? Время существования ~ 8 лет

34

Эффект от начального значения параметра ?0 на примере ИСЗ ПРОГНОЗ-6

ПРОГНОЗ-6 (a*=16.6, ?*=0.117) ?0 = 268?, ?0 = 0.126 ?1(t0)=c1=0.05, ?1(t0)=c2 =-0.19

Характер эволюции параметров ?, ? и ?1, ?2 в зависимости от даты старта (определяющей значение параметра ?0)

Реальный запуск 22.IX.1977, ?0 = 225?. Время существования ~ 40 лет

Гипотетический запуск 22.IX.1988 , ?0 = 80?. Время существования ~ 7 лет

Эволюция орбитальных элементов гипотетической версии орбиты ПРОГНОЗ-6

с датой старта 22.09.1978 ?0= 247.5?. Время существования более 500 лет

Rp

??

i?

??

??

?1, ?1m

?2, ?2m

? = 90

?1

?=1

?2

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

1978-2150

2150-2470

2150-2470

Заключение

Заключение

Сопоставление аналитических решений двукратно осредненной проблемы Хилла с решениями, учитывающими возмущения от реальных внешних тел, позволило выделить параметры, от которых зависит характер эволюции орбитальных элементов и время баллистического существования ИСЗ, обусловленное гравитационными возмущениями со стороны внешних тел (Луны и Солнца). Такими параметрами являются безразмерные константы первых интегралов двукратно-осредненной задачи c1 (0? c1 ? 1), c2 (-0.6 ? c2 ? 0.4), безразмерный параметр 1 < a*, равный отношению большой полуоси орбиты спутника к радиусу центрального тела, и параметр ?0 (0?? ?0 ? 360?) – начальное угловое расстояние между восходящими узлами орбит ИСЗ и Луны на эклиптике.

37

38

Список литературы

Лидов М.Л. Эволюция орбит искусственных спутников планет под действием гравитационных возмущений внешних тел. // Искусственные спутники Земли. 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.