CC BY
44
ЧЕБЫШЕВСКИЙ СБОРНИК
Том 14 Выпуск 1 (2013)
УДК 531.391
ДВИЖЕНИЕ СИСТЕМЫ «ПЛАНЕТА С ЯДРОМ — СПУТНИК» В ГРАВИТАЦИОННОМ ПОЛЕ СИЛ ВЗАИМНОГО ПРИТЯЖЕНИЯ
А. В. Шатина, Е. В. Шерстнев (г. Москва)
Аннотация
Исследуется движение спутника в гравитационном поле вращающейся планеты. Планета моделируется телом, состоящим из твердого ядра и вязкоупругой оболочки из материала Кельвина — Фойгта. Спутник моделируется материальной точкой. Из вариационного принципа Даламбера — Лагранжа выводится система интегро-дифференциальных уравнений движения механической системы в соответствии с линейной моделью теории упругости. С помощью асимптотического метода разделения движений строится приближенная система уравнений движения в векторном виде, описывающая динамику механической системы с учетом возмущений, вызванных упругостью и диссипацией. Выводится усредненная система дифференциальных уравнений, описывающая эволюцию параметров орбиты спутника. Для частных случаев построены фазовые траектории, найдены стационарные решения и исследована их устойчивость. В качестве примеров рассмотрены некоторые планеты Солнечной системы и их спутники.
THE EVOLUTION OF A SATELLITE MOTION IN THE GRAVITATIONAL FIELD OF A VISCOELASTIC PLANET WITH A CORE
Shatina A. V., Sherstnev E. V.
(Moscow State Technical University of Radio engineering,
Electronics and Automation)
Abstract
We investigate the motion of a satellite in the gravitational field of a massive deformable planet. Planet is modeled as body that consists of a solid core and a viscoelastic shell of a Kelvin-Voigt material. The satellite is modeled as a point mass. The system of integro-differential equations for a motion of a mechanical system is got out from the variational principle of the
d'Alembert-Lagrange according to the linear theory of elasticity. Approximate equations of motion in vector are constructed with asymptotic method of motions separation. This system of equations describes the dynamics of the "planet-satellite"with regard to the perturbations caused by elasticity and dissipation. To describe the evolution of the orbital parameters of a satellite, averaged differential equations were derived.
Phase trajectories were constructed for particular cases, their stationary solutions were found and investigated on stability. In the case of the existence of two stationary orbits stationary solution that corresponding to the motion along the orbit of larger radius is asymptotically stable, and the orbit of smaller radius is unstable. Some of the planets in the solar system and their satellites are considered as examples. This problem is a model for the study of the tidal theory of planetary motion
Исследование приливной эволюции системы <планета — спутник> проводилось многими авторами [1, 2, 3]. В работе используются методы аналитической механики систем с бесконечным числом степеней свободы [4]. Ранее указанный подход был применен, в частности, к ряду задач о поступательно-вращательном движении вязкоупругого шара [5, 6, 7, 8].
1 Постановка задачи. Уравнения движения
Рассмотрим задачу о поступательно-вращательном движении системы <планета — спутник> в гравитационном поле сил взаимного притяжения. Спутник будем моделировать материальной точкой P с массой ц. Планету будем моделировать телом, состоящим из твердого ядра и вязкоупругой оболочки, занимающим область V = V0 U V1 в трехмерном евклидовом пространстве при отсутствии деформаций. Здесь V0 = { r G E3 : |r| ^ Го}, V1 = { r G E3 : r0 < |r| ^ ri}. Пусть p0, pi — плотности ядра и вязкоупругой оболочки соответственно, а m0, m1 — их массы. Предполагается, что материал оболочки планеты является однородным и изотропным.
Введем инерциальную систему координат OXYZ с началом в центре масс системы <планета-спутник>. Для описания вращательного движения планеты введем подвижную систему координат Cx1x2x3 жестко связанную с ядром и систему осей Кенига C^1^2^3, где C — центр масс планеты в естественном недеформированном состоянии (рис. 1).
Положение точки M планеты в инерциальной системе координат OXY Z определяется векторным полем
Rm (r,t) = OC + Г(г + u(r,t)), (1)
где Г — оператор перехода от подвижной системы координат Cx1x2x3 к системе осей Кенига C^1^2^3, u(r, t) — вектор упругого смещения, равный тождественно
Рис. 1: Иллюстрация задачи
нулю для точек твердого ядра У0. Так как О — центр масс рассматриваемой
механической системы, то
J Им (г, ^ рйу + /1 ■ ОР
V
Здесь р =
!
Ро, г е У0
рг, г е У1
Введем в рассмотрение вектор И = СР. Тогда из (1) и (2) получим:
ОС = -
ОР
И
1
Тир1йь1,
VI
1
Гир1йу1
VI
Здесь т — масса планеты, т = т0 + т1.
Потенциальная энергия гравитационного поля определяется функционалом
рйу
V
—И + Г (г + и)
где f — универсальная гравитационная постоянная.
Функционал потенциальной энергии упругих деформаций зададим в соответствии с линейной моделью теории упругости:
Е = j Е [и] йу1, Е [и] = а1 (IЕ — а211Е)
VI
0
где
Е (1 - V) 2(1 - 2v)
а1 — 77ТЇ—,—Тлі—> а2 — —;-> а1 > °> 0 < ак < 3,
2(1 + V) (1 — 2v) 1 — V
— еіі’ пе — ^ Ыей — е2к1) , еы — 2 ^ дк + , и = (щ} ^ ^
j=1 к<1
где Е — модуль упругости Юнга, V — коэффициент Пуассона вязкоупругой оболочки планеты, 1е, 11е — инварианты тензора малых деформаций.
Диссипативные свойства вязкоупругой оболочки опишем диссипативным функционалом
$ = J $ [и] дю1, $ [и] = хЕ [и] ,
VI
соответствующим модели Кельвина-Фойгта (здесь х > 0 — коэффициент внутреннего вязкого трения).
Положим Ир = ОР. Уравнения движения системы <планета-спутник> получим из вариационного принципа Даламбера-Лагранжа:
^Им, 8Им^ Р^у + і ^Ир, 8Ир^ + 8П +
м, «Им
V
+ J (VuE [и] + V*$ [и] ,8и) ^1 = 0 (4)
VI
Согласно равенствам (1) и (3) имеем:
Им =--------1— И-1—Г I [о х [о х и] + 2о х и + о х и + и} р1йю1
т + 1 т + 1 У
VI
+Г [о х [о х (г + и)] + 2о х и + (о х (г + и) + и} ,
8Им —-----8И —------Г I [8а х и + 8и} р\$У\ +
т + ц т + і ]
Уі
+Г [8а х (г + и) + 8и} , (5)
т1
Ир —----И--------Г [о х [о х и] + 2о х и + о х и + и} р\йь\,
т + ц т + і ]
Уі
т1
8ИР —-----8И------Г [8а х и + 8и} р^ь\.
т + і т + і У
Уі
Здесь о — вектор угловой скорости планеты, 8а — вектор, возникающий при варьировании ортогонального оператора Г:
о х (•) — Г-^О , 8Г (•) — Г [8а х (•)].
Подставляя в равенство (4) выражения (5) для Им, 8Им, Ир, 8Ир и приравнивая коэффициенты при независимых вариациях 8И, 8а, 8и, получим уравнения движения системы <планета-спутник> в виде:
і [ ^ , іт - „ Г И — Г (г + и)
И м рЛю +--------IIР + /і --------------------------з рЛю — 0, (6)
т + і У т + і У |И — Г (г + и
УУ
г + и —
1
У
т +1
ирі Лю1
Уі
х Г 1ИмрЛю —
і
ир1Лю1 х Г 1ИР —
т +1
Уі
— Iі
У
(г + и) х (Г 1И — (г + и)) |И — Г (г + и)|3
рЛю — 0, (7)
Р- < Г-‘И“ — т+11 Г-‘И"Р‘к — т+1 Г-‘ИР — ^ |И —Г~(г(++и)|>з+
+ VuE [и + хи] = 0. (8)
2 Возмущенная система уравнений движения
Будем полагать, что жесткость вязкоупругой оболочки планеты велика, т.е. мал безразмерный параметр е = рц^Е-1, где ш0 — величина модуля начальной угловой скорости планеты. Выбрав соответствующим образом масштабы размерных единиц, можно ввести малый параметр е = Е-1. При е = 0 вектор упругого смещения и полагается равным нулю. В этом случае получим задачу о движении механической системы, состоящей из абсолютно твердого тела сферической формы и материальной точки, в поле сил взаимного притяжения. Невозмущенная система уравнений движения имеет вид:
И+ И = 0, Ао = 0,
л3
где А = 11 [рог0 + р1 (г5 — гц)] — момент инерции планеты в недеформирован-ном состоянии относительно диаметра.
Используя метод разделения движений [4] получим векторную систему дифференциальных уравнений, описывающую поступательно-вращательное движение системы <планета-спутник> с учетом возмущений, вызванных упругостью
и диссипацией. Будем искать вектор-функцию и из (8), описывающую квази-статические деформации планеты под действием внешних сил и сил инерции в виде: и — єи1 + є2и2 + • • •.
Было найдено первое приближение и — еи1. После подстановки этого решения и преобразований получим векторную систему дифференциальных уравне-
ний:
И +
/(т + , 3/(т + і)р\єО
Я3
И +
б/і
Я3
тЯ
1 + 3ХЯ^ | + 6*№
Г {—5 І (і, о) + |о2 + 2о(|, о)+
0,
Я )' ' Я3 *}
І х I + (І , о)[о х і] |
:ю)
0
п 4пг7 { \
В —--------Ф(х, V),
105 ^
ф(х, V) — { —1б(9к + 14)х17 — 200(3к + 8)х14 + б72(4к + 9)х12—
Ао
— (210к2 + 3044к + 5824)х10 + (525к2 + 125бк + 1576)х7 + 84(17к + 12)х5 —
— 25(21к2 + 92к + 5б)х3 + 210к2 + 71бк + 416} .
Здесь 1 — вектор кинетического момента планеты относительно центра масс:
Ь — [ Г(г + и) х Л [Г(г + и)] рЛю. (11)
3 аЬ
V
Система уравнений (9)—(10) имеет первый интеграл — закон сохранения момента количеств движения системы <планета-спутник> относительно общего центра масс:
т,гИ х И + Ь — Со,
где тг — , С0 — постоянный вектор.
3 Эволюция орбитального движения спутника
Так как с точностью до членов порядка е вектор кинетического момента планеты имеет вид Ь = АГ^ , то Г^ = А 1Ь. Тогда при учете уравнений (9) и
(11) получим векторное дифференциальное уравнение орбитального движения спутника:
рИ = Ео + еЕ1 + е\^2,
где
Е ____ /|0Р И
Ео = — Я?И
Го0 — тгИ х И) 2 (И, С0 ) Го0 — тгИ х И
Е1 = —СА^--------—-----------------------}— И + 4 1
А2 Я5
А2 Я5
5И(И, С0 )2 6/р
А2 Я7
—С2 * Я+Я И—
+
Я8
И
}
Сю — т,г И х И
И
АЯ8
С1 = З^рр^т , С2 = 6/рСь /0 = /(т + р).
Для получения эволюционной системы уравнений движения спутника перейдем к каноническим переменным Делоне Ь, О, Н, д, к [10, 12]. После усреднения уравнений по быстрой угловой переменной I — средней аномалии — получим замкнутую систему обыкновенных дифференциальных уравнений относительно переменных <действие> Ь, О, Н и медленных угловых переменных д, к.
Выпишем эволюционную систему уравнений орбитального движения спутника в безразмерных переменных п0, е, г, д, к, где п0 = пАО-1 ( п — среднее движение по орбите, е — эксцентриситет, г — наклонение орбиты и д — долгота перигелия от восходящего узла). Ось OZ инерциальной системы координат OXYZ направим по вектору С0 = (0, 0,00). Указанная система ОДУ имеет вид:
3Д1п16/3
п 0
(1 — е2
Д1п0^3 е (1 — е2)13/2
)15/2 |
Р (л 2\1/2
сое г--------1/3 (1 — е ) '
Р 2 1/2 сое г-----------1/3 (1 — е ) '
1/3
■ F2(e) ■ (1 — е2)3/2 — п0 ■ *3(е)
п,
1/3
■ ^5(е) ■ (1 — е2)3/2 — п0 ■ ^4(е)
}
&
Д1п03/3 в1п
(1 — е2)
яп г Г1 /9 3 . 2 \ 2 /5 1 . 2 N Л
—— < —+--------------вт д е +------------------------вт д е > ,
2)5 \ 2 V4 2 у V16 4 ) )
}
Д2п0
Зсов г--------1/3 (1 — е2)1/2 \ +
(1 — е2)3/2
п.
1/3
ЫпТ в1п г Г 5 2 1
+ Ц—ё?\2С082 г — 2 +
15рп(
(т + р)(1 — е2)3
0
н = -
Д1п“/3 . 2 | Зе2 е4\
(1 - е2)5 ' 31П 9 \ ~4 + “8 } + (1 - е2)3/2 (1 - е2)2 ■
Д2Щ Д3п0/3 соя і
Здесь
Ді —
ШХ^р2 о /с<Л 16/3
т(т + м)/о V А
о2/Ч А) ’
Д2 —
Р —
Аі/3/
Зє^р2 0с3
(т + ^)А4
2/3
Д3 —
ЗєріБ (Со
т/?3 У
Со \
А
і3/3
тг
С
4/3
^ / ч , 15 2 45 4 5 6 , 31 2 255 4 185 6 25 8
г(е) - 1 + уе9 + -8е‘ + 16е‘, гзм . 1 + -е9 + -+ -е» + 64е8,
^ / ч ~ 135 2 135 4 45 6 ^ , ч 11 33 2 11 4
г,(е)- 9 + —е2 + -+ 64е=, я(е) — у + 7^ й.*.
Рассмотрим два частных случая движения спутника: а) г = 0, б) е = 0. В
обоих случаях п0 — п0* является корнем уравнения:
по +
Р
п
і/3
— 1,
которое в зависимости от значения параметра р либо не имеет решений, либо имеет одно решение, либо имеет два решения П01 и П02. При этом стационарное решение п01 асимптотически устойчиво, а п02 неустойчиво.
Фазовые портреты для частных случаев систем уравнений при значении параметра р — 0.375 изображены на рис. 2.
а) б)
Рис. 2: Фазовые портреты для случаев: а) г = 0, б) е = 0
о
о
В таблице 1 представлены численные значения параметра р, стационарные значения п01, п02 и значение величины п0 — п0 (0) в настоящее время для различных систем <планета-спутник>. Безразмерная переменная пропорциональна среднему движению спутника по орбите и связана с большой полуосью орбиты спутника а соотношением: а — /01/3/ (С0А-1 п0)2/3.
Для всех приведенных примеров, за исключением системы Марс-Фобос, имеет место двойное неравенство: п01 < п0(0) < п02. Значение переменной п0 во время движения уменьшается. Это означает, что большие полуоси орбит спутников увеличиваются, стремясь к асимптотически устойчивым стационарным значениям. При этом для спутников Юпитера и спутника Марса Деймоса текущее значение по(0) ближе к неустойчивому стационарному значению п02, а для системы Земля-Луна по(0) ближе к асимптотически устойчивому значению п01. Для системы Марс-Фобос п01 < п02 < по(0). Значение переменной п0 увеличивается. Это означает, что Фобос приближается к Марсу.
Таблица 1: Параметры систем <планета — спутник>
Планета — спутник Р п>(0) п01 п02
Земля — Луна 0.15535 7. 2 сл 00 со о 1 со 3.7922 • 10-3 0.83503
Марс — Фобос 1.5235 • 10-6 3.2171 3.5358 • 10-18 1 —1.5235^10-6
Марс — Деймос 2.5391 • 10-7 8.1268•10”1 1.6369 • 10-20 1 — 2.5391 • 10-7
Юпитер — Ио 5. 8 СО 00 0 1 о 1 0.2335 2.0517 • 10-10 0.9994
Юпитер — Европа 3.1685•10“4 0.1163 3.1811 • 10-11 0.9997
Юпитер — Ганимед 9.7643•10-4 5.7650•10-2 9.3093 • 10-10 0.9990
Юпитер — Каллисто 7.0902 • 10-4 2.4717 • 10-2 3.5644 • 10-10 0.9993
Заключение
В рамках исследования модели системы <планета — спутник> получены уравнения движения указанной системы. Эти уравнения учитывают возмуща-ютттие эффекты, вызванные упругостью и диссипацией, которые влияют на динамику эволюции элементов орбиты согласно приливной теории.
Для двух частных случаев движения, когда наклонение орбиты равно нулю и когда эксцентриситет орбиты так же равен нулю, найдены стационарные решения и исследована их устойчивость. Показано, что в случае существования двух решений, одно из них является асимптотически устойчивым, а второе неустойчивым. Для пар <планета — спутник> Солнечной системы вычислены соответствующие значения и определён характер их движения.
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
[1] Приливы и резонансы в Солнечной системе. Сб. статей / под ред.
В.Н. Жаркова М.: Мир, 1975. 287 с.
[2] Белецкий В.В. Движение спутника относительно центра масс в гравитационном поле. М.: Изд-во МГУ, 1975. 308 с.
[3] Марков Ю.Г., Миняев И.С. Роль приливной диссипации в движении планет и их спутников // Астрон. вестн. 1994. Т. 28. №2. С. 59-72.
[4] Вильке В.Г. Аналитическая механика систем с бесконечным числом степеней свободы. Ч 1,2. М.: Изд-во мех.-мат. факультета МГУ, 1997. Ч.1. 216 с.; Ч.2. 160 с.
[5] Вильке В.Г. Движение вязкоупругого шара в центральном ньютоновском поле сил // ПММ. 1980. Т.44. Вып. 3. С.395-402.
[6] Шатина А.В. . Эволюция движения вязкоупругого шара в центральном ньютоновском поле сил // Космич. исследования. 2001. Т. 39. №3. С. 303315.
[7] Вильке В.Г., Шатина А.В., Шатина Л.С. Эволюция движения двух вязкоупругих планет в поле сил взаимного притяжения // Космич. исследования. 2011. Т. 49. №4. С. 355-362.
[8] Шатина А.В., Шерстнев Е.В. Движение спутника в гравитационном поле вязкоупругой планеты // ПММ. 2012. Т. 76. Вып. 6. С. 913-922.
[9] Лейбензон Л.С. Краткий курс теории упругости. М.;Л.: Гостехиздат, 1942. 304 с.
[10] Мюррей К., Дермотт С. Динамика Солнечной системы. М.: ФИЗМАТЛИТ, 2010. 588 с.
[11] Куликовский П.Г. Справочник любителя астрономии. М.: Кн. дом, 2009.
[12] Дубошин Г.Н. Небесная механика. Основные задачи и методы. М.: Наука, 1975. 799 с.
Московский государственный технический университет радиотехники, электроники и автоматики МГТУ МИРЭА
Поступило 03.04.2013
