Математическая модель движения спутников Марса Текст научной статьи по специальности «Математика»

Банъщикова M. А., Титаренко E. Ю.

13

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДВИЖЕНИЯ СПУТНИКОВ МАРСА

*

М. А. Баньщикова, Е. Ю. Титаренко

Научно-исследовательский институт прикладной математики и механики, Томского государственного университета, 634050, Томск, Россия

* Институт кибернетики Национального исследовательского Томского политехнического университета, 634050, Томск, Россия

УДК 521.1

Представлены результаты численного моделирования движения спутников Mapсa (Фобос, Деймос) с использованием дифференциальных уравнений движения в элементах Роя. Проведено сравнение с моделью движения в прямоугольной системе координат.

Ключевые слова: спутники Mapсa, численное моделирование, орбитальные параметры.

The results of numerical modeling of motion of Mars' satellites (Phobos and Deimos) by using the differential equations in Roy's orbital elements are presented. Comparison with motion model in rectangular system of coordinates is carried out.

Key words: Mars' satellites, numerical modelling, orbital parametres.

Введение. Исследование спутников Mapсa Фобоса и Деймоса представляет большой интерес для астрономов. Это малые объекты, которые быстро движутся вблизи плоскости экватора планеты по почти круговым орбитам с периодом 0,3 и 1,2 сут на расстоянии приблизительно 9000 и 23 000 км от Mapсa соответственно. Разработано большое количество динамических моделей движения спутников Mapсa [1-4]. Первые численные модели движения спутников Mapсa появились в конце 1990-х - начале 2000-х гг. [5, 6].

Поскольку спутники быстро вращаются вблизи планеты, численное моделирование динамики этих объектов, основанное на пошаговом интегрировании дифференциальных уравнений в прямоугольных координатах, сопряжено с известными трудностями. Правые части уравнений представляют собой быстроизменяющиеся функции времени, поэтому численное интегрирование таких уравнений необходимо выполнять с достаточно малым шагом, что приводит к большому объему вычислений, а также к интенсивному накоплению ошибок округления в численном решении. Кроме того, дифференциальные уравнения в прямоугольных координатах сильно неустойчивы по Ляпунову, что является причиной быстрого (сверхлинейного) увеличения количества методических ошибок в процессе численного интегрирования. Так как уравнения Роя (со средней долготой) регулярны и обладают стабилизирующим эффектом [7], можно предположить, что численное моделирование на их основе будет значительно более эффективным по сравнению с моделированием на основе классических уравнений движения.

В данной работе представлена численная модель движения спутников Mapсa Фобоса и Деймоса, построенная на основе дифференциальных уравнений в элементах Роя, а также прове-

ден сравнительный анализ эффективности численного моделирования движения спутников Марса в прямоугольных координатах и элементах Роя.

1. Формальное представление модели движения спутников Марса. Формально численную модель спутникового движения рс в пространстве угловых координат р = (р1,р2) относительно стандартного земного экватора J2000.0 можно представить в виде [8]

рс = рс (/, ) = Та, х(/, qD ), дт ). (1)

Здесь t - эфемеридное время; цот — (цD, ) - вектор всех параметров модели; Т - преобразование перехода от ареоцентрической системы координат в пространство угловых координат на небесной сфере, отнесенной к топоцентру; цв — (х0, х0, t0, д8,...) , - параметрические векторы, связанные соответственно с движением спутника и с координатным преобразованием; х0, х0 - векторы динамического состояния небесного тела в начальный момент времени t0.

2. Уравнения движения спутников в элементах Роя. В модели (1) положения х и скорости х определяются через элементы Роя е(И, А, X) [9] следующим образом:

х3к1 . _ 1

, X — 2

к + к г

X — гсовX + 3 , Х1 — —(х3к2 - х2кз + гГх^),

х2 — гвшX + Хзк2 , хХ2 — —х1к3 -х3к + гГх2),

к + к3 г

г1

х3 — — (к1 сов X + к2 вт X), х3 — — (х2к1 - х1к2 + ггх3).

к г

Здесь

к $ л /2 2\1/2

г —-, г — -, А — (с2 + $2),

ц + с к

с — у сов X + 5 вШ X, $ — у вт X-5 С0в X,

у — А —А3к1_, § — А2 - М

к + к3 к + к3

Элементы Роя вычисляются путем численного интегрирования дифференциальных уравнений, которые можно представить в виде [9]

—ёА — (вх1Л ,(. х ^ ^ к + кк2 - к2к

— х X Р, -— (р X И) + (х X11),

& ' ^ у ; 1 '' & г1 к (к + к3)'

— Г, (2)

ёе

&

с начальными условиями, составляющими начальный вектор состояния динамической системы в пространстве элементов е0 — е(^). В (2) И - оскулирующий момент количества движения ( к — | И |); А - вектор Лапласа ( А —| А |); X - истинная долгота исследуемого тела (сумма долготы восходящего узла, аргумента перицентра и истинной аномалии). Элементы Роя выражаются через прямоугольные координаты по следующим формулам:

Банъщикова М. А., Титаренко Е. Ю.

15

а

а г. км А г. км

10 ю2

10 1

ю-1 ю-2 ю-3 ю-4

ю-5

10б ю-7 10"г

Е Q 01

1 \ □ 2

= \

1 Lü

=

III

I

I

= Зп

1

........ ........ ........ ........

10 102 10 1

Ю-1 ю-2 ю-3 ю-4 ю-5 10"® ю-7 10"8

= 171

i

1 9 \

i \ \

1 b \

I V \

! Л H

1 CD ^ П

i \

1 \ 0a

...... ........ ........

дг^ jq-9 Е_i.........i.........i...........

10" 103 104 Ю5 10" 102 ю3 ю4 ю5

Характеристики точность - быстродействие для спутников Марса: а - Фобос, б - Деймос; 1 - уравнения Роя, 2 - уравнения в прямоугольных координатах

h — X3 X3 , — X3 Xi — X3 , h^ — X X2 X2 Xi,

A1 — Dx1 - Kx1, A2 — Dx2 - Kx2, A3 — Dx3 - KX3, (3)

cos X — A. - sin X— - K — |x | • | x |, D —|x |2 --^.

| x | | x | h + h3 | x | | x | h + h3 | x |

Функция P вычисляется как сумма различных возмущающих сил, действующих на движение спутников Марса:

P — Рм + P, + Pp (4)

( PM - гравитационное поле Марса; PS - притяжение Солнца; Pp - притяжения планет Солнечной системы и Луны).

3. Оценка эффективности уравнений движения. Для сравнения представленной модели с существующими моделями выполнено численное исследование движения спутников. Сравнение проводилось с моделью в прямоугольной системе координат. Дифференциальные уравнения (2)-(4) интегрировались методом Эверхарта 12-го порядка с переменным шагом в компьютерной арифметике с двойной точностью. Движение спутников рассматривалось на временном интервале, составляющем один год. В течение этого промежутка времени Фобос совершает приблизительно 1000 оборотов, Деймос - 300 оборотов.

Проведена оценка эффективности численных моделей, характеризуемой точностью вычисляемых координат спутников и быстродействием их получения. В качестве меры быстродействия численного моделирования было принято число шагов интегрирования NS . Ошибка интегрирования оценивалась как разность численного и точного решений Дг — | Ax |.

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

жения точности Дг от 1 до 1G-7 км требуется приблизительно в 1G раз меньше шагов интегрирования NS, чем при решении дифференциальных уравнений в прямоугольных координатах.

Следует отметить, что с увеличением объема вычислений (с увеличением NS ) возрастает ошибка округления (несистематическое поведение характеристик), которая фактически определяет уровень наивысшей точности, достижимой в арифметике с двойной точностью. Таким образом, поскольку уравнения Роя позволяют получать заданную точность при меньшем числе шагов, для них уровень наивысшей точности оказывается почти на порядок выше, чем для классических уравнений в прямоугольных координатах. На рисунке видно, что для Фобоса уравнения Роя позволяют получить точность порядка Ю-7 км, тогда как уравнения в прямоугольных координатах - не более 1G-6 км.

Заключение. Таким образом, в работе получены оценки эффективности использования уравнений в элементах Роя и в прямоугольных координатах при численном моделировании орбитального движения спутников Марса. Результаты показали, что уравнения в элементах Роя позволяют выполнять численное моделирование значительно более эффективно по сравнению с классическими уравнениями в прямоугольных координатах. В частности, одна и та же методическая точность в представлении вектора положения для уравнений Роя достигается на порядок быстрее.

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

1. SINCLAIR A. T. The orbits of the satellites of Mars // Vistas Astronom. 197S. V. 22, pt. 2. P. 133-14G.

2. Shor V. A. Motion of Martian satellites // Preprint Sternberg Astronom. Inst. 19S9. V. 13. P. 61-75.

3. Chapront-Touze M. Orbits of the Martian satellites from ESAPHO and ESADE theories // Astronom. Astrophys. 199G. V. 24G. P. 159-172.

4. Emelyanov N. V., Vashkovyak S. N., Nasonova L. P. The dynamics of Martian satellites from observations // Astronom. Astrophys. 1993. V. 267, N 2. P. 634-642.

5. БОРДОВИЦЫНА Т. В., Шарковский H. А. Эффективный алгоритм численного моделирования движения Фобоса, спутника Марса // Изв. вузов. Физика. 1994. № 1G. С. S-12.

6. Lainey V., Dehant V., Pätzold M. First numerical ephemerides of the Martian moons // Astronom. Astrophys. 2GG7. V. 465, N 3. P. 1G75-1GS4.

7. Avdyushev V. Numerical stabilization of orbital motion // Celest. Mech. 2GG3. V. S7, iss. 4. P. 3S3-4G9.

S. Авдюшев В. А. Численное моделирование орбит. Томск: НТЛ, 2G1G.

9. Рой A. E. Движение по орбитам. М.: Мир, 19S1.

Банъщикова Мария Александровна — канд. физ.-мат. наук, ст. науч. сотр.

Научно-исследовательского института прикладной математики и механики Томского государственного университета;

Титаренко Екатерина Юрьевна - ст. преп. Института кибернетики Томского политехнического университета; e-mail: am@am.tpu.ru

Дата поступления - 12.G9.12 г.