УДК 521.1
Возмущения первого порядка во вращении Земли, обусловленные гравитационными моментами со стороны Луны при высокоточном описании ее орбитального движения
© Ю.В. Баркин1, М.Ю. Баркин2, 3
1 ГАИШ МГУ им. М.В. Ломоносова, Москва, 119991, Россия
2 МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
3 Московский авиационный институт, Москва, 12599, Россия
На основе уравнений вращательного движения в переменных Андуайе для несферической Земли разрабатывается аналитическая теория вращательного движения (теория прецессии, нутации и вынужденных колебаний полюсов). Построены возмущения первого порядка во вращении планеты под действием гравитационных моментов со стороны Луны в условиях ее реальной орбиты, задаваемой высокоточной теорией ЕЬ 421. Вычисления, выполненные в эклиптической системе координат даты, представлены в форме, удобной для анализа.
Ключевые слова: переменные действие-угол, невозмущенное движение Эйлера — Пуансо, задача Лиувилля, ряды Фурье, эллиптический интеграл.
Введение. В предыдущих работах авторов изучалось вращательное движение Земли как изолированной планеты с изменяемой геометрией масс на основе специального подхода, учитывающего специальные формы канонических уравнений в переменных Андуайе и действие-угол (гамильтонов формализм) [1, 2]. Гравитационные моменты, действующие на несферическую Землю со стороны Луны, Солнца и планет при этом не учитывались. В рамках данного подхода можно обобщить выполненное исследование и в возмущающий гамильтониан задачи добавить силовую функцию гравитационного взаимодействия несферической Земли с Луной и другими небесными телами. Это позволяет построить новую теорию прецессии и нутации, поскольку в качестве невозмущенного вращательного движения Земли принимается не осевое вращение (вращение вокруг полярной оси инерции планеты), а коническое движение, подробно описанное в работе [1]. Следует отметить, что возмущения первого порядка, обусловленные наблюдаемыми вариациями тензора инерции Земли [1], и возмущения, обусловленные силовой функцией ньютоновского взаимодействия несферической Земли и Луны, в качестве которой мы принимаем вторую гармонику, имеют одинаковый порядок и харак-
_7
теризуются малым параметром ц ^ 10 . А это означает, что возмущения первого порядка для указанных возмущающих факторов могут быть вычислены отдельно в соответствии с принципом суперпозиции.
Разложение силовой функции в задаче о вращении Земли.
Силовая функция в теории вращательного движения Земли определяется формулой для второй гармоники гравитационного потенциала [3]:
U = 3n¡(A - C)[a^ (2 + а3), (1)
где n 2 = fmSj a3; r — расстояние между центрами масс Земли и Луны; а — невозмущенное значение большой полуоси орбиты Луны; 5 — динамический параметр,
8 = (B - A)I(C - A) > 0; A, B, C — главные
центральные моменты инерции Земли, соответствующие ее главным центральным осям инерции CE,, Сц, CC,; а2, а3 —осинусы углов, образованных осями инерции Cц и CC, с линией центров Земля — Луна. В переменных Андуайе эти направляющие косинусы определяются выражениями
al = \jc1l cos (K-h) + c2l cospsin (K-h) - c3i sinpsin (K-h^^
+ (c2i sin p + C3i cos p) sin ф, i = 1,2,3. (2)
В выражениях (2)
c11 = cos g cos l - sin l sin g cos 9;
c21 = sin g cos l + sin l cos g cos 9;
c12 = - sin l cos g - cos l sin g cos 9;
c22 =- sin l sin g + cos l cos g cos 9; (3)
c13 = sin g sin 9;
c23 = - cos g sin 9;
c31 = sin l sin 9;
c32 = cos l sin 9;
c33 = cos 9
являются направляющими косинусами осей инерции планеты CE, ^ и CC, в промежуточной системе координат Андуайе CG1G2G3, связанной с вектором кинетического момента вращательного движения Земли G. Ось CG3 направлена вдоль вектора G, оси CG1 и CG2 расположены в промежуточной плоскости Q, ортогональной
вектору кинетического момента Земли, ось С01 направлена к восходящему узлу плоскости плоскости эклиптики даты на основной координатной плоскости эклиптики эпохи (какпринято в известной теории вращения твердой Земли Киношита [4]. Переменные р и к — наклон и долгота восходящего узла плоскости Q в основной эклиптической системе координат, 9 — угол между полярной осью инерции Земли и вектором кинетического момента ее вращательного движения. Невозмущенное значение этого угла является малым и в приведеных исследованиях в соответствии с данными наблюдений принимается 9 = 0,24".
В выражение силовой функции (3) сомножителями входят функции (а / г)3 а2 и (а / г)3 а2. Для этих функций получены следующие представления [3]:
+
2
+ ( sin2 р + с21еЪ1 sin 2р + с3у р) эш2 ф | +
2(ей - 4 р - 4 sin2 р + с21сЪ1 sin со§2 ф cos 2 (Х - к))
+ (С2, cos р - СцСЪ1 sin р) -II а I cos2 Ф sin 2 (Х - к) и
+2((С^ sinр + С1;С3;cosр)| [ — 1 з1пфcosфcos(Х-к) |> +
"[(г - Сзг )sin2р + 2с21С31 COS2р Л"-) БШ Ф cos ф вт (Х- к)
(4)
Выражения в фигурных скобках зависят от сферических эклиптических координат гелиоцентрических координат Луны, которые в свою очередь представляются известными тригонометрическими рядами по кратным средних долгот планет. Эти ряды построены с учетом планетных возмущений в орбитальном движении Луны с высокой точностью [3].
Полученные данные представляют собой новые разложения функций (5.9), (5.10а), (5.10Ь), (5.11а) и (5.11Ь) (от сферических ко-
ординат орбитального движения Луны r , ф и Xm ) из работы Киношита [4]. Разложения выполнены на интервале времени 1000AD-3000AD на основе долгосрочной численной эфемериды Луны LE-406 (система координат ICRF) и представлены в виде рядов Пуассона [3], [5]. Эти разложения, имеющие более общий вид, определяются более точными формулами по сравнению с рядами Киношита и учитывают новые эффекты в орбитальном движении Луны:
-2 Г -1 (1 - 3sin2 ф) = £ A^ cos 0V + at0) sin 0V; (5.1) (5)
a| cos2 фcos2(b-h) = Y A¡,2) cos0V + aV2) sin0V; (5.2)
r J V
cos2 фsin2(b-h) = 242) sin0V+ ¿V2) cos0V; (5.3)
a
43
— I sinфcosфsinh) = ^a|,1) cos0V + aV1) sin0V; (5.4)
r r
a I sinфcosфcosh) = ^bV1) sin0V + bV1) cos0V. (5.5)
В выражениях (5) ©у — любые линейные комбинации классических аргументов теории орбитального движения Луны,
0у=у1/ + v2/' + v3F + у4 D + у5 (О- h), (6)
V
V
V
V
Здесьу = (V!,у2, ..., у5) — целочисленные коэффициенты; /, /', ^ D, О — аргументы Делоне и средняя долгота восходящего узла орбиты Луны, определяемые, соответственно, выражениями (3.5.Ь) и (3.4.Ь.3) из работ [3], [5].
Коэффициенты рядов (5), (6):
А/), 47) и 47), ^) (7)
с высокой точностью представлены квадратичными функциями времени и учитывают вековые планетные возмущения в орбитальном движении Луны [3], [5]:
д( у) = д( у) , д( У Ь , д( У) 12
где
А = (А, В, а, Ъ), у = (0,1,2). (8)
9 3 2
Фундаментальная постоянная /шЕ = 398 600,5 • 10 м /с . Значение большой полуоси а = 383 397 772,5 м было принято при построении рядов Пуассона (5), (6). Для фундаментальной частоты п0 = V №Еа-32 и для значений частот орбитального движения Луны пп, пр, входящих сомножителями при времени в выражениях аргументов теории орбитального движения Луны, имеем следующие численные значения:
п0 = 17 311 058"6 854 464 ед./год;
па = -69 679" 193 631 ед./год;
пр = 17 395 272"628 478 ед./год,
Как и в теории физической либрации Луны [3] за начальный момент времени принимаем юлианскую дату 12000.0 (ШБ2451545.0). В рядах Пуассона (5), (6) коэффициенты
АЙ, а^ (у = 0,1,2) и В£, Ъ^П (у = 1,2) (9)
представляют собой численные коэффициенты (амплитуды) (у = 0,1,2; п = 0,1,2).
Ниже приведены [3], данные коэффициенты, умноженные на соответствующую степень, для каждой из пяти шаровых функций (5.1) -(5.5):
А^0, а^ (у = 0,1,2), , Ъ$, у = 1,2 (107);
А^, ай (у = 0,1,2), , , у = 1,2 (108); (10)
а$, ((=0,1,2), у=1,2 сю9>.
В первом столбце указан номер наборов коэффициентов, в последующих пяти столбцах по строкам приведены наборы пяти целочисленных индексов у = (уь у2, у3, у4, у5) — коэффициентов в линейных
комбинациях классических аргументов теории орбитального движения Луны (6). В седьмом столбце указаны значения периодов (в сутках) соответствующих возмущений или тригонометрических слагаемых разложений (при условии, что пренебрегаются вековые квадратичные изменения аргументов со временем). Здесь не приводятся таблицы указанных значений коэффициентов из разложений (9), (10),
(11). Отметим лишь, что для этих разложений получены таблицы значений коэффициентов (5.1) для 286 наборов индексов аргументов ©v = v1/ + v2/' + v3F + v4D + v5 Q, в разложении (5.2) — для 336 наборов, в разложении (5.3) — для 332 наборов аргументов, в разложении (5.4) — для 525 наборов индексов аргументов и, наконец, в разложении (5.5) — для 526 аргументов ©v.
После громоздких преобразований на основе формул (1)-(7) получаем окончательное тригонометрическое разложение второй гармоники силовой функции задачи:
U = 2nl(A -C) jNoZ[R),v (P, t)cos0V + r0,v (p, t)sin 0v] +
+ZNo,2^Z[Ro,v(P,t)cos(0v -2м/)-ro,v(P,t)sin(0v -2м/)] +
м v
+ N1,0 Z Z [ Rv (P, t )cos (g - 8 0v ) + /^ (P, t )sin (g-80v)^^
8 v
+ N2,0 Z Z [R2Ev (P, t) cos (2g - 80v ) + r^v (P, t) sin (2g - 80v )
8 v
+Z N1, м Z Z [Rv(P, t)cos (g - 80v+ 2м/) +
M 8 v
+r^v t) sin (g -80v+ 2м/)] +
Z N2, м Z Z [ REÍ (P, t) cos (2 g - 80v + 2м/) -
8 v
+/g (P, t) sin (2 g - 80v + 2м/).
+
м
(11)
Коэффициенты Nn,m (9) и функции наклона Rn m (p), rn,m (p) при
тригонометрических функциях в разложении (11 ) являются углов наклона 9 и P. Они определяются последовательностью формул:
N0 = (2 + 5) (2 - 3sin2 9); N10 = i(2-5)sin29;
1 2 3 2
N2o = -(5-2)sin29; NQ2 = --5sin29; (12)
4 4
1 12 N1,м=-45sin9(м + cos9); N2^ = -^м5cos9(1 + мcos9) .
Для функций наклона запишем следующие представления:
R° v(p,0 = -1 (зcos2 р-1)aV0) - 1sin2p^í/1) - jsinW2); 6 v ' 2 4
r° V(p,t) = -1 (3cos2p-1)0) - 1sin2pa(1) - 1sin2paV); 6 ^ '2 4
RV = sin2p|aV0) -2aV2)|-2cos2pA^ + 2scospB^ - esinpB^); r/V) = 2cosp¿V-1) - sinp¿V2) - e sin2p|aV0) - 2a(,2) |
+2ecos2pa
(1).
R2eV = 42) + sin2 p | aV0) -142) | - sin 2pA^ + 2e sin pBV1) + e cospBV2);
r2(ev) = 2 sinp¿V1) + cospé|,2) - eaV2) - e sin2 p | aV°) -1 a(2) | + e si
sin 2pa
(1)
(e = ±1).
(13)
Полученное разложение силовой функции задачи (11)—(13) обобщает аналогичное разложение из работы Киношита [4].
Возмущения первого порядка во вращательном движении Земли в переменных Андуайе.
Уравнения вращательного движения твердой Земли под действием гравитационного притяжения Луны имеют вид:
dG = SU;
dt Sg '
de 1 _ . j 1 11 . 1 nsu 1 nsu
— = — G sin el---I sin2/ н— ctg e---cosec e-;
dt 2 l A B J G S g G S /
dp 1 S U 1 S U
— = — ctg p---cosec p-;
dt G S g G S h
dt
dg dt
í
= -G cos e
2 Л
1 sin / cos /
C
B
1 ÜS U н— cosec e-:
g Se
(14)
= G
2
sin / cos / -+-
AB
J
1 „ SU 1 SU
- Gctg Gctg ;
dk 1 д и
— = — cosec р-,
dt G др
где силовая функция задачи определяется разложением (11)—(13).
Уравнения (14) легко приводятся к стандартному виду уравнений многочастотных колебательных систем, содержащих малый параметр (или несколько малых параметров). В качестве таких малых параметров принимаются обычно малые динамические сжатия эллипсоида инерции Земли или коэффициенты второй гармоники геопотенциала J2, С22 .
В отличие от теории Киношита в данной работе принимается более общее невозмущенное вращательное движение Земли, а именно равномерное осевое вращение планеты и коническое чандлеровское движение вектора кинетического момента относительно полярной оси инерции (конус с углом полураствора 0 = 0,24"). Период движения полюса твердой осесимметричной Земли составляет около 305 сут. Однако это движение (и уравнения враща-тельного движения) легко обобщается на случай Земли с упругой мантией и период блуждания полюса относительно полярной оси инерции уже будет составлять 432 сут. [2], как и для чандлеровского движения полюса.
Таким образом, невозмущенное вращательное движение рассматриваемой модели Земли определяется простыми формулами:
/ = щг + /о; g = V + &о; Ь = % 0 = 0,24"; р = Ро = -23°45. (15)
Отметим, что указанные переменные Андуайе определяют движение вектора кинетического момента по отношению к главным центральным осям инерции Земли (средним осям); ng, щ — невозмущенные (постоянные) значения частот эйлеровского — чандлеровского движения [2].
Для основных параметров задачи можно принять следующие значения [2]. Экваториальный радиус Земли 637 8140 м, I = с/тг02 — безразмерный момент инерции. Коэффициенты второй гармоники геопотенциала: 32 = 1082,6265 • 10-6 , С22 = 1,81537-10-6. Параметр 5 = = 4С22/(32 + 2С22 ) = 0,668479 • 10-2.
Согласно рассматриваемому методу построения приближенного решения задачи для возмущений первого порядка, имеем простые квадратуры:
•д и
ОХ = \dUdf; д g
dt;
1 „Гди . 1 „г ди
01 =— йе 0-dt--еоБее 0 -
1 О 6 } д g О ) д /
1 гди. 1 гди ,
р1=—р.!-^ -—с08есР)~дЬ
1 ди
/1 = — совес 0J-dt
(16)
О д0
ди , 1 гди
^ =- —* - —Р^др
dt;
Ъ1 = —сое ее р|
ди
О * др
dt.
Отметим, что в подынтегральных выражениях и в правых частях (уравнений) (16) переменные Андуайе принимают невозмущенные значения (15). Постоянные интегрирования определяются сответствую-щими начальными условиями задачи, которые для краткости не рассматриваются.
Возмущения первого порядка переменных Андуайе. В результате подстановки в подынтегральные выражения в (18) невозмущенных значений переменных Андуайе (17) вычисление интегралов сводится к интегрированию известных тригонометрических функций времени. Поэтому, опуская довольно громоздкие преобразования, приведем окончательные формулы для возмущений первого порядка, обусловленных гравитационным притяжением Луны при ее сложном орбитальном движении (второй гармоникой). Для простоты в выражениях коэффициентов А (8) сохраним лишь их постоянные слагаемые (из общего разложния силовой функции) и примем постоянные значения коэффициентов геопотенциала, У2 и С22, которые связаны с соответствующими разностями главных центральных моментов инерции А, В, С известными формулами:
_ 2С - А - В _ В - А
•2 --Г-2-; С22 ---
(17)
2т0г0 4да0г0
где т0, г0 — масса средний радиус Земли.
Запишем ряд формул, определяющих возмущения первого порядка переменных Андуайе (О, р, 9, /,И):)
Чо ЕЕ
О---^(•2 + 2С22К X О 4 1&к 2 22' 0
^(Р,*) ( п К Г,у)(Р,Г) • ( п \
' 008 (д -вПу ) + —-— БШ (д-вПу)
пд -
пд - вОу
+Е ^1, .ЕЕ
я&Чр, г)
пд - вОу + 2.п/
Г, ?(Р, г)
соб (д -вПу + 2|/) +
+2Е ^ЕЕ
пя - вОу + 2|п/
^(р, г)
(18)
2пд - вОу + 2|п/
г2(5(р, г)
Бт (д -вПу + 2|/)
соб (2 д -вПу +2|/) +
2п„ - вОу + 2|п/
Бт (2 д - вПу + 2|/)
Pi = _ 2 (j2+2C22) 4 IqV 2 22 '
n x
_ N0(G)csc p^v5
R0,v (P, t) © , r0,v (P, t) • ©
cos ©v + —-sin ©v
a
a
+ N0,2(G)csc P^^v5 I v
, r0,v QP t ) •/ © 2 A +—-sin (©v _ 2|I/)
R0,v (p, t )
cos (©v _ 21/) +
av _ 2 |in/
Rl(,sv(p, t ) ng _ sQv
av _ 2|n/ + N1,0 csc p^^ (cos p + sv5 ) :
cos (g _ s(v ) + r1, v (p,t) sin (g _ s©v
ng _s£2v
+N20 csc p^^(2cos p + sv5 ) x
s v
r2sv(p,t) , ч r1(v}(p,t) '
' -cos (2 g _s©v) + ^-— sin (2 g_s©v)
2ng _sQv
2ng _ sQv
+Z N1 ,I cscp^^(cosp + sv5)
I s v
Rl(,sv)(p, t )
ng _ sQv + 2|n/
cos (g _ s©v + 2|/ ) +
+-
rj v)(P, t )
ng _ sQv + 2|n/
sin (g _s©v +2|/ )
+Z N2 ,I csc ^ (2 cos p + sv5 ) :
(19)
r2 sv(p, t )
2ng _ sQv + 2|n/
cos (2 g _s©v + 2|/) +
ÍvV t )
2ng _ sQv + 2|n/
sin (2 g _s©v + 2| / )
v
9^1,0 ^^
Е
+ 2cot 9^2,0 ЕЕ
91 =- 3 ТШ (2 + 2С22 )п0 X ^(р,t) ( п )+ ^У)(р,t) • ( 0 ч
, ( - Е0у ) + ,-— Sin (-Е0у)
па - еоу
па -еоу
4еу(р, t) ч г2(е.)(р, t)
2пё - еоу
2п§ -е оу
-2csc 9N0,2 (9)Е Е ц
ц у
, г0,ц (P, )
-ОУ + 2цп
Я),у (р, t) -ОУ + 2цп/
•вт (0у - 2ц/)
cos(0v - 2ц/) +
(20)
+Е( 9- 2ц csc 9)Nl, цЕЕ
ц Е у
. 1(у)(р, о
яечр, t)
п„ - еОл, + 2цп/
п8 - еоу + 2цп/ sin ( - Е0у + 2ц/)
cos ( - Е0у + 2ц/) +
+2Е(cot 9-ц csc 9)),цЕЕ
Я ЕУ(р, t)
2п„ - еОу + 2цп/
cos (2 £-Е0у + 2ц/) +
+-г£(р/)-sin ( - Е0у + 2ц/)
2п„ - еОу + 2цп/
/1 =- (J2 + 2С22 )п0
N0 (9)Е
4 Iш sin 9 я, у (р, t)
Оу
sin 0у - Г°,у (р, t) cos 0у
у Оу у
+^1,0(9)ЕЕ
Е У
+^,0( 9)ЕЕ
Е У
+ ^,2(9)ЕЕ
^(р,t) • ( 0 ) яеу t) ( 0 )
' вШ ( - Е0у )----— COS (-Е0У)
па - еоу
па - еоу
+ (21)
ЕУ(р,о • (2 0 ) г&уt) (2 0 )
' -вт ( - Е0у ) - — :-— cos (2^-е0у)
2п„ - еОу
2п„ - еОу
ц у
^0,У (P, О • /0 2 А Г0,У (P, t) /0 2 А -вт(0У -2ц/)--:-cos(0У - 2ц/)
ОУ - 2цп/
ОУ - 2ци/
х
у
+1N, i (G)SS
r(,sàp, t )
ng -s Qv + 2|nl l( ?(P, t )
sin (g - s©v + 2|l ) -
+E N2,| (G)2S
ng -s Qv + 2|nl
R sV(p, t )
cos (g -s©v + 2|l )
(2l)
2ng - sQV + 2|nl ^(P, t )
sin (2g - s©v + 2|l) -
2ng - sQV + 2|nl
cos (2 g -s©v + 2|l )
Qv
Qv
<jcot GNo (e)£ + cot pNo(e)^
V
+cot eM,o(e)££
s V
r( s}(P, t)
gl - 3~T(J2 + 2C22)n0 x 4 I ю
^^ sin ©V- ro^ cos ©V
Qv
R0,v(p Г )sin ©V- rkM cos ©V
Qv
R?(p, t )
ng - sQV
sin (g -s©v)-
ng - sQV
cos (g - s©v )
+ cot pNl,o(e)2S
Rl(vs)(P, t) ■ ( © ) rl(vs)(p, t) ( © )
' sin (g -s©v)--!-— cos (g-s©v)
ng - sQV
ng - sQv
+cot eN2,o(e)££
s V
rg(p, t )
RsV(p, t) .
2ng -s Qv
sin (2g - s©v ) -
2ng - sQV
cos (2g - s©v )
+ cot pN2,o(G)^^
s V
rjSV t)
R2(s)(p, t )
2ng - sQV
sin (2g - s©v ) -
2ng - sQV
cos (2g - s©v )
+ 9^,2(9)ЕЕ
ц у
г0,у (P, t)
ОУ - 2цп/
Я0,у (р, t) ОУ - 2цп/
cos (0У - 2ц/)
вт (0Л, - 2ц/) -
+ cot рN0,2 (9)Е Е
ц у
Я0,у (р, t)
ОУ - 2цщ
sin (0У - 2ц/) -
г0,у (P, t) ОУ - 2цп/
cos (0У - 2ц/)
+Е ^ц (9)ЕЕ
Я(,?(р, t)
п„ - еОу + 2цп/
sin ( - Е0У + 2ц/) -
Г1(,У)(р, t)
п„ - еОу + 2цп/
cos ( - Е0У + 2ц/)
+Е ^ рNl, ц (9)ЕЕ
Я1(уЕ)(р, t)
п„ -е ОУ + 2цп/
-вт (-е0у + 2ц/ - (22)
Г"(уЕ)(р, t)
п„ - еОу + 2цп/
cos ( - Е0У + 2ц/)
+Е cot 9N2,ц (9)ЕЕ
Я ЕУ(р, t)
2п„ - еОу + 2цп/
sin (2 g - Е0У + 2ц// -
Г2(,Еу}(р, t)
2п„ - еОу + 2цп/
^ (2 g - Е0У + 2ц/)
+Е cot рN2ц (9)ЕЕ
я2(У}(р, t)
2ng - еОу + 2цп/
вт (2g -Е0 + 2ц)-
€}(р, t)
2п„ - еОл, + 2цп/
cos (2g - Е0Л, + 2ц/)
к = - Т^-^ (2 + 2С22 ) х
х| е)£
4 I ш б1п р ^ (р, 1)
а
• © 'О, V (р, 1)
Б1П ©v--:-СОБ ©V
v а
<Е)(р, 1). ( ©) '¿е)(р, 1) ( ©) ' Б1П ( -S©v)--:-— СОБ (-©V)
па -
-
я2(еЛр, 1) .
2п„ - вО,,
б1П (2 д -©V)-
2п„ - вО.,
СОБ (2 g-в©v)
+^о,2(
ц v
^ (р, 1) •/© 2 А (р, 1) /© 2 л Б1П (©v - 2ц/)----СОБ (©v - 2ц/)
Qv - 2цпл N1, ц (
ц Е V
>йЕ)(р, 1)
Qv - 2цп/
п„ - + 2цпл
б1п ( - в©„ + 2ц/) +
(в)2Е
пя - + 2цпл
R2(V)(р, 1)
СОБ ( -Е©V + 2ц/)
2п„ - + 2цпл
б1п (2 g-Е©V + 2ц/)
€}(р, 1)
СОБ (2 д -Е©V + 2ц/)
2п„ - + 2цпл
(23)
Суммирование в формулах (18)—(23) осуществляется по всем значениям индексов у = (у1,у2, ..., у5) из разложений (5), (6) и по двум вспомогательным индексам: в = ±1 и ц = ±1. Здесь все аргументы тригонометрических функций являются известными линейными функциями времени. Переменные /, д, к принимают свои невозмущенные значения (15). Коэффициенты Nп,т (в) и функции наклона Я„,т (р), а также их производные Н'„,т (в) и Я'п,т (р) по переменным Андуайе в и р при тригонометрических функциях в формулах для возмущений (18)—(23) являются постоянными. Их можно вычислить при постоянных значениях соответствующих углов
(24)
9 = 90 и р = р0 (они определяются последовательностью формул (12) и (13).
Первые производные данных коэффициентов определяют по формулам:
N0 (9) = -3(2 + 5) sin 29; Ní,0(9) = 1 (2 - 5)cos 29;
#2,о(9) = 4 (5 - 2)sin29; N0,2(9) = - 4 5 sin29;
N1, ц (9) = - 4 5( cos 9 + cos29);
N2, ц (9) =1 ц5 sin 9(1 + ц cos 9) (1 + 3ц cos 9).
8
Для их производных по наклонности р получены следующие формулы:
R0,v(р, t) = 2sin 2р fAÍ0) - 2Л^ - cos2рЛ«; r0,v (р, t) =1 sin 2р f aV0) - 2aV2) J - cos 2ра(1);
Я^е)(р, t) = 2 cos 2р fЛ^ - 2 Л2) J + 4sin 2рЛ^ - 2s sin рB^ - e cosрbV2) ; ^(р, t) = -2sinр¿V1) - cosр¿V2) - 2e^2р^0) -2a?^ - 4esin2рaV1; RfVV,t) = sin2р| AV0)-1 aV2) |-2cos2рЛV1) + 2ecosрbV1) -esinрвV2);
r2'(V}(p, t) = 2cosр- sinpb[2) - e sin2p[40) - 2«V2^ + 2e cos2paV.
(e = ±1).
В формулах для возмущений (18)-(23) Qv = v1nM + vn + v3nF + +v4njD + v5nQ — комбинации частот основных аргументов теории орбитального движения Луны, ng, n¡ — частоты невозмущенного
эйлеровского движения.
В теории вращения Земли [4] получены вспомогательные разложения определенных функций сферических геоцентрических координат Луны: формулы (5.9)-(5.11b), а квадратичные слагаемые относительно времени с коэффициентами Аг, aVj'2, (j = 0,1,2) и В('2,
bv2 (( = 1,2) вообще не рассматриваются. В настоящей работе также
не используются упрощающие предположения Киношита, касающихся полноты учета вековых орбитальных возмущений Луны [4].
Процедура вычисления и табулирования амплитуд, периодов и наборов аргументов для возмущений первого порядка всех переменных Андуайе является весьма емкой, хотя принципиальных трудностей не вызывает и будет рассмотрена в следующей работе авторов.
Заключение. Построены возмущения первого порядка во вращательном движении несферической Земли непосредственно в переменных Андуайе. В качестве базовых использованы уравнения движения в переменных Андуайе, в которых учтены вторая гармоника силовой функции при высокоточном описании орбитального движения Земли и Луны. Приближенное решение задачи о вращении Земли построено с помощью метода малого параметра.
ЛИТЕРАТУРА
[1] Баркин М.Ю. Изучение возмущенных вращательных движений небесного тела с приложением к теории вращения Земли. Дис. ... канд. физ.-мат. наук. Москва, 2014, с. 154.
[2] Barkin Yu.V. (2000) Perturbated rotational motion of weakly deformable celestial bodies. Astronomical and Astrophysical Transactions, 2000, vol. 19, issue 1, pp. 19-65. DOI: 10.1080/10556790008241350.
[3] Barkin Yu.V., Kudrjavtsev C.M., Barkin M.Yu. Perturbations of the first order of the Moon rotation. Proceedings of International Conference "Astronomy and World Heritage: across Time and Continents". Kazan, 2000, 19-24 August. KSU, pp. 161-164.
[4] Kinoshita H. Theory of Rotation of the Rigid Earth. Celest. Mech, 1977, vol. 15, pp.277-326.
[5] Kudryavtsev S.M. Long-term harmonic development of lunar ephemeris. Astronomy & Astrophysics. A&A 471, 1069-1075. DOI: 10.1051/ 0004-6361:20077568.
Статья поступила в редакцию 29.10.2014
Ссылку на эту статью просим оформлять следующим образом: Баркин Ю.В., Баркин М.Ю. Возмущения первого порядка во вращении Земли, обусловленные гравитационными моментами со стороны Луны при высокоточном описании ее орбительного движения. Инженерный журнал: наука и инновации, 2014, вып. 12.
URL: http://engjournal.ru/catalog/fundamentals/hidden/1337.html
Баркин Юрий Владимирович — д-р физ.-мат. наук, профессор, ведущий научный сотрудник ГАИШ МГУ им. М.В. Ломоносова. Область научных интересов: теоретическая механика, небесная механика. е-mail: barkin@inbox.ru
Баркин Михаил Юрьевич — ассистент кафедры «Теоретическая механика» МГТУ им. Н.Э. Баумана, ассистент кафедры «Теоретическая механика» Московского авиационного института. Область научных интересов: теоретическая механика, небесная механика. е-mail: barkin@yandex.ru
First-order perturbation in the rotation of the Earth caused by the gravitational moments from the Moon when describing the precision of its orbital motion
© Yu.V. Barkin1, M.Yu. Barkin2,3
1 Sternberg SAI, Lomonosov Moscow State University, Moscow, 119991, Russia 2
Bauman Moscow State Technical University, Moscow, 105005, Russia
3
Moscow Aviation Institute, Moscow, 12599, Russia
The article presents the analytical theory of a rotary motion (the theory of a precession, nutation and the compelled fluctuations of poles) developed for non-spherical Earth on the basis of the equations of a rotary motion in Anduaye's variables. We constructed perturbation of the first order in rotation of the planet under the influence of the gravitational moments from the Moon in the conditions of its real orbit set by the high-precision theory of EL 421. All calculations are executed in ecliptic system of coordinates of date and presented in the form convenient for the analysis. Tabulation of perturbations for all variables of Anduaye is executed.
Keywords: action-angle variable, Liouville problem, Fourier series, elliptic integral. REFERENCES
[1] Barkin M.Yu. Izuchenie vozmuschennykh vraschatelnykh dvizheniy nebesnogo tela s prilozheniem k teorii vrascheniya Zemli [Study of perturbed rotational motions of the heavenly bodies with application to the theory of the Earth's rotation]. Ph.D. Thesis (Phys.&Math.). Moscow, 2014, pp. 154.
[2] Barkin Yu.V. Pertuibated rotational motion of weakly deformable celestial bodies. Astronomical and Astrophysical Transactions, 2000, vol. 19, issue 1, pp. 19-65. DOI: 10.1080/10556790008241350.
[3] Barkin Yu.V., Kudrjavtsev C.M., Barkin M.Yu. Perturbations of the first order of the Moon rotation. Proceedings of International Conference "Astronomy and World Heritage: across Time and Continents". Kazan, 2000, 19-24 August. KSU, pp. 161-164.
[4] Kinoshita H. Theory of Rotation of the Rigid Earth. Celest. Mech, 1977, vol. 15, pp.277-326.
[5] Kudryavtsev S.M. Long-term harmonic development of lunar ephemeris. As-
tronomy & Astrophysics. A&A 471, 1069-1075. DOI: 10.1051/00046361:20077568.
Barkin Yu.V., Dr. Sci. (Phys.&Math.), professor, leading researcher in Sternberg Astronomical Institute at Lomonosov Moscow State University. Sphere of scientific interests includes theoretical mechanics, celestial mechanics. e-mail: barkin@inbox.ru
Barkin M.Yu., assistant lecturer of the Theoretical Mechanics Department in Bauman Moscow State Technical University, assistant lecturer of the Theoretical Mechanics Department in the Moscow Aviation Institute (Technical University). Sphere of scientific interests includes theoretical mechanics, celestial mechanics. e-mail: barkin@yandex.ru