Периодические и условно-периодические движения спутника-гиростата под действием гравитационного момента на круговой орбите Текст научной статьи по специальности «Физика»

УДК 531.352

А. А. Панкратов ПЕРИОДИЧЕСКИЕ

И УСЛОВНО-ПЕРИОДИЧЕСКИЕ ДВИЖЕНИЯ СПУТНИКА-ГИРОСТАТА ПОД ДЕЙСТВИЕМ ГРАВИТАЦИОННОГО МОМЕНТА НА КРУГОВОЙ ОРБИТЕ

Проведено исследование периодических движений осесимметричного спутника-гиростата, центр масс которого движется по кеплеровой круговой орбите. В качестве переменных была использована модификация переменных Андуайе, в уравнениях движения спутника относительно центра масс учитывали только гравитационный момент. Система дифференциальных уравнений рассматриваемой задачи содержит малый параметр, характеризующий близость осевого и экваториального моментов инерции спутника-гиростата, а также малость проекции собственного кинетического момента гиростата на экваториальную плоскость. Уравнения движения задачи, допускающие интеграл Якоби, с помощью которого осуществляется понижение порядка системы на две единицы, приведены к так называемой форме Уиттекера. Методом Пуанкаре доказано существование многопараметрических семейств периодических решений приведенной системы, найдены порождающие периодические решения.

E-mail: [email protected]

Ключевые слова: спутник-гиростат, переменные Андуайе, гравитационный момент, метод Пуанкаре, периодические и условно-периодические движения.

Твердое тело с полостями, заполненными идеальной несжимаемой жидкостью, совершающей потенциальное движение, и тело-носитель, в котором жестко закреплена ось вращения симметричного ротора, являются наиболее известными, но не единственными примерами механической системы, называемой гиростатом. Понятие гиростата неоднократно обобщалось и уточнялось. Работа Н. Е. Жуковского [1], опубликованная в 1885 г., сыграла важную роль в формировании представлений об основных характеристиках движения гиростата. Жуковский показал, что потенциальное движение жидкости в полости определяется движением тела, само же движение тела происходит так, как если бы жидкость была заменена эквивалентным твердым телом. Несмотря на то что многие классические задачи динамики твердого тела удалось обобщить на случай гиростата, специфические свойства движения гиростата остаются еще недостаточно изученными.

Постановка задачи. Под спутником-гиростатом будем понимать твердое тело (движущееся по орбите) с расположенными внутри него симметричными роторами (маховиками). Оси вращения роторов жестко связаны с несущим телом и совпадают с осью симметрии каждого ротора. Предположим, что угловые скорости роторов относительно несущего тела постоянны.

В настоящей работе проведено исследование различных семейств периодических решений в задаче о вращательном движении спутника-гиростата, центр масс которого движется по круговой орбите, собственный кинетический момент роторов мал, а эллипсоид инерции спутника-гиростата близок к сфере. Для нахождения периодических решений были использованы достаточные условия Пуанкаре [2, 3] и условия, полученные в работах [4, 5]. Для этого уравнения движения приведены к форме, необходимой для применения метода малого параметра Пуанкаре. Для получения требуемых уравнений сначала необходимо вывести дифференциальные уравнения вращательного движения спутника-гиростата в углах Эйлера и сопряженных канонических переменных, а затем с помощью канонических преобразований перейти к безразмерным переменным, типа переменных Анду-айе. Анализ найденных периодических решений ограничен порождающими решениями.

Уравнения движения. Рассмотрим движение динамически симметричного спутника-гиростата вокруг центра масс £ под действием гравитационной силы тяготения центрального тела О (однородного шара), относительно которого центр масс спутника описывает кепле-ровскую орбиту радиусом К. Для определения движения спутника-гиростата введем следующие системы координат:

инерциальную Охуг с началом координат в точке О и осями Ох, Оу, расположенными в плоскости орбиты;

декартову £хуг с началом координат в точке £ — центре масс спутника-гиростата и осями £х, £у, £г, параллельными осям Ох, Оу, О2 соответственно;

подвижную £%т]С, оси которой направлены по главным центральным осям инерции спутника-гиростата.

Пусть А, В, С — моменты инерции спутника-гиростата (вместе с роторами) относительно осей £%, £п, £С соответственно. Обозначим К%, Кп, К£ проекции собственного кинетического момента роторов на оси £%, £п, ££". Тогда К = (К%, Кп, К^)— кинетический момент роторов относительно системы координат £%г/С. Отметим, что в рассматриваемой задаче К%, Кп, К^ — постоянные величины.

Ориентацию системы координат £%т]С относительно системы £хуг зададим углами Эйлера р, у, в. Кинетическая энергия вращательного движения спутника-гиростата с точностью до аддитивной постоянной имеет вид

Т = 1 (Аа( + Ва>П + С|) + К(% + Кпю71 + К^ю^, (1)

где ап, — проекции угловой скорости спутника на оси £%, £п, £С- Зависимости проекций ап, а^ от углов р, у, в и их производных заданы кинематическими уравнениями Эйлера. Углы Эй-

лера можно рассматривать как обобщенные координаты, которые вместе с обобщенными импульсами pp, py, pe определяются формулами

дТ дТ дТ

pp=dp • py=w< • pQ=dQ> <2)

и образуют в рассматриваемой задаче систему канонических переменных.

Опуская промежуточные выкладки, запишем гамильтониан задачи Н* в канонических переменных p, у, Q, pp, py, pe:

1 ( pp Л2 (sin2 p cos2 pi

H* = 2Uq"p(ctgQ,J brJ +pQsinpcospx

xl-py-ppctge^i 1 - !) +1 pQ(cos^p+sinpp -

lsinQ Fp 5 ЛA b) 2FQl A B ) 2C

py ctgQY K^sinp + KncosP + Knsinp К4 cosp sinQ Pp * A A B B A

- ppf- U (p, y, Q). (3)

Далее выполним преобразование (Q, p, y, pQ, pP, py) ^ (l', g', h', L', G', H'), где G' — модуль вектора кинетического момента вращательного движения спутника-гиростата, т. е. вектора G' = (Am¿¡ + K^)i + (Bav + Kn)j + (Ca^ + K^)k; H' — проекция вектора G ' на ось Sz; L' — проекция вектора G на ось S£. Пусть £ — плоскость, проходящая через центр масс спутника-гиростата и перпендикулярная вектору G. Тогда l' — угол между осью и линией пересечения плоскости Z с плоскостью S^rj; g' — угол между линиями пересечения плоскости Z с плоскостями S^r/ и Sxy; h' — угол между осью Sx и линией пересечения плоскостей Z и Sxy. Переменные l', g', h' являются углами, изменяющимися по модулю 2п. Отметим, что связь между импульсами pQ, pp, py и L', G', H'

можно получить с помощью формул (1), (2), учитывая определения переменных L' , G', H' , l' , g', h':

Pe

= V G '2 - L'2 sin(/'-p), pp = L' , pw = H(4)

Используя формулы сферической тригонометрии [3], получаем следующее дифференциальное соотношение:

dg' = cos6d(<-1') + cospd(y- h') - sin(<-1')sm6d6, (5)

L' H' G'2 -La . G'2 -L'2

cosO =—, cos p = —, sinO = J-, sin p = J-.

G' G' V G' V G '

Умножив соотношение (5) на величину G', на основании формул (4) имеем

G 'dg' + L'dl' + Hdh' = p(dp + pydy + podO.

Из этого соотношения следует, что преобразование (в, р, у, po, Рр, Ру) ^ (l', g', h', L', G', H') является однородным каноническим

преобразованием [3].

Используя определения канонических переменных L', G', H', l', g', h', соотношения (4) и формулы сферической тригонометрии, получаем

Aa^ = G' sin Osinl' - ^^ = >/G '2 - L'2 sinl' - K#,

„ = G ' sin O cos l' - Kv = VG2-L2 cos l' - Kn, Ca^ = G ' cosO-K^ = L ' - K^.

С помощью формул сферической тригонометрии можно получить формулы, связывающие углы Эйлера с переменными Андуайе [3]. Здесь приведем необходимую в дальнейшем запись направляющих косинусов a13, a23, a33 через переменные L', G', H', l', g', h ':

«13 =

f 1 -(f I cosg + f 1-Í£T

Л

sin h +

+ A^ -f f J sin g COs h^

Г

«23 =

f '-ff Jcosg + f 1-Íf)2

Л

sin h +

(6)

+ AI1 -f f J sin g 'sin h^

«33 = fr -i -ff) J1 -(§} cos g'

Таким образом, кинетическую энергию и силовую функцию U, а следовательно, и гамильтониан задачи (3) можно записать в переменных L', G', H', l', g', h'. С точностью до аддитивной константы имеем

' = (G'2 - L'2 )( sin21' cos21' 1 L* K VG'2 - L'2

H ' =

5

+-+ Kp-sin/-

2C p A

-^G'2 - L'2 cosl' - ^L ' - U. (7)

5 С

Ограничимся следующим общепринятым [6] приближенным значением силовой функции:

U = f [(Л - В)у'2 + (A - С )у"2 ], (8)

где f — гравитационная постоянная; m — масса спутника-гиростата; R = R - модуль радиус-вектора центра масс спутника-гиростата; у' — косинус угла между радиус-вектором R и осью спутника Srf; у" — косинус угла между радиус-вектором центра масс спутника R и осью SZ. В рассматриваемом случае R = const (центр масс спутника-гиростата движется по круговой орбите) и моменты инерции А и В равны (спутник осесиметричный). Тогда силовую функцию (8) можно представить в виде

U = |о2(A - С)у"2. (9)

Здесь о — среднее орбитальное движение спутника. Так как

X y z R • ez

Y = a13 —+ a23 —+ a33 — =--,

R R R R

где ez — орт оси Sz; x, y, z — координаты центра масс спутника-гиростата,

x = R cosot' , y = R sinot' , z = 0,

то

y" = a13 cosot' + a23 sinot'. (10)

Формулы (6), (9) и (10) позволяют представить силовую функцию U явной функцией элементов Андуайе и времени t'. Опуская выкладки, окончательно представим силовую функцию в тригонометрическом виде:

U (L', G' , H' , g' , h', t' ) =

(A-C) 2 Ukl,b(ß,p)cos[kig' + k2(h'-01')].

ki, k2

Здесь суммирование проводится по индексам к = 0, 1, 2; = 0, ± 2, а коэффициенты разложения Uk1,к2(в, р) определяются следующими формулами:

3 3 ( 3 ^ 3

U00 = — sin2 в +—\ 1 — sin2 в I sin2 р, U20 = — sin2 $sin2 р, 0 4 41 2 J 2 8

3 3 ( 3 ^

U1.0 = ^sin2psin2e, U0.2 =-4^ 1 -2sin2 eJ sin2 р, (11)

3 3 U2 ±2 =--sin2 в(1 ± cosр)2, U1 ±2 = +—sin2$sinр(1 ±cosр).

16 8

Выполним каноническую замену переменных И" = H', G'' = G', L " = L', h'' = h' - mt', g'' = g', l'' = l', что соответствует переходу к вращающейся вокруг оси Oz с угловой скоростью m системе координат Oxyz (в этой системе центр масс спутника-гиростата покоится). В результате такой замены уравнения движения принимают автономный вид. Для упрощения последующих выкладок перейдем от переменных L'', G'', И", l", g'', h " к новым (безразмерным) каноническим переменным L, G, И, l, g, h и безразмерному времени t:

L' G'' И''

L =-, G = —, И =-, l = l'', g = g '', h = h'', t = cot',

л ' л ' л ' " О О " " "

Ac Ac Ac

H пп = '

пр " Аф2 '

где Н пр — преобразованный гамильтониан; Н ''— гамильтониан (7),

записанный в соответствующих переменных. С учетом изложенного запишем гамильтониан (7) для случая А = В:

Нпр = — + А—— 1} -Н -(K^sin/ + K„cos/)-пр 2 2С Аф к 4 п ' (12)

Cm

Здесь

KzL - U(L, G, H, g, h).

и - — G, H, g, h) = 2 Utbk2(p,d)cos(g + k2h). A- A ki, k2

Введем вспомогательные обозначения ж0, а, в, определяемые формулами

^K— ^K ^K—

——= ж0 sin^sina, —-= ж0 sinPcosa, ——= ж0 cos^, A- A- C-

тогда гамильтониан (12) можно представить в виде

G 2 A - C

Нпр =~2" + 2C L2 -H - x0G[sinesin#cos(/-а) + cos^cos^]-

-Ü(L, G, H, g, h), (14)

или

Hпр = "Y + A2CcL2 - H - жо {{2 - L2 sin ecos(/ - a) + L cos^} - U

(15)

Запишем дифференциальные уравнения вращательного движения спутника-гиростата в виде

d/ = д Н dg = д Н dh = д Н d Н пр, d Н пр, d ~~т Н пр,

(16)

dL=_дН dG=АН H=АН dt д/ пр' dt dg пр' dt dh пр

Поскольку при значениях G = 0, |H| = G или |L| = G уравнения (11),

(13), (15), (16) имеют особенность, рассмотрим их на множестве D х T3:

D = {G, р,в: 0 <е1 < G; 0 <е2 <р<п-е2; 0 <е3 <6<n-s3},

(17)

T3 = {/, g, h mod 2n}.

Без ограничения общности рассматриваемой задачи в выражениях

(14) и (15) можно положить а = 0 (т. е. K— = 0) . Выберем в качестве малого параметра динамическое сжатие спутника-гиростата / = (A - C)/C и модуль собственного кинетического момента гиростата Ж0 = /ж. Тем самым приведем гамильтониан (14) (или (15)) к стандартному виду:

Нет — Нет о + ДНсг! + ДНсг 2 + •••

(18)

Нст0 — —--Н

Нст1 — L- - x ((—2 -L2 sin ßcosl + L cos ß)

- 2 Uki,fc(P,0)cos(kig + k2h);

ki, k2

L2

Hст2 —--xLcosß,...

(19)

Периодические решения приведенной системы. Система (16) с гамильтонианом (18), (19) допускает интеграл Якоби

G- - H + /Нст1 (L, G, H, l, g, h) + /2Нст 2 +... = Ci = const, (20)

который используем для понижения порядка системы на две единицы, путем приведения уравнений (16) к так называемой форме уравнений Уиттекера [7]. Для этого решим уравнение (20) относительно переменной Н и получим

Н = Нст — Но + л Hi (L, G, l, g, h) + ¡л2 H2 +...

(21)

где

Н — —2 Г Но — —:— Г1,

2

Н1 — Нст1 (L, G, Н, l, g, h)|

I G2 ;

H — Нст0 — —— Г '

Н2 —

д

Нст 1 (L, G, Н, l, g, h) — Нст1 + Нст 2

дН

G2

Н — Нст о ——— Г

(22)

(n-1 1 йп-k-1 ,

Z1 _д_/Н \n-k

к k—1 (п - k )!дНп-k-1 (Н ст k+11

Н п —

G2

Н — Нст о ———Г

Так как дНсто/дН ^ 0, ряд (21) является сходящимся при достаточно малых значениях параметра ¡л.

Приняв угловую переменную к за новую независимую переменную, уравнения вращательного движения спутника-гиростата запишем в виде

ёО _дН _ дН ёк д е' ёк д1 '

ф __дН ^ _ дн аи дО аи дь'

где гамильтониан Н определяется формулами (21), (22).

В случае если будет найдено какое-либо решение системы (23), то соответствующее ему решение исходных уравнений (16) можно представить формулами

H = H(L, G, /, g, h) = H(h), j

hQl dH CT

dH

dh

-= t (24)

L=L(h),/=/(h), G=G(h), g=g(h), H = H (h)

При значении л_ 0 уравнения (21)—(23) несложно проинтегрировать, т. е.

О _ Оо, Ь _ Ьо, g __Оок + go, I _ 1о, (25)

где Оо, Ьо, go, ¡о — постоянные интегрирования. При рациональном значении Оо решение (25) будем называть периодическим. Рациональность величины Оо означает соизмеримость между невозмущенным значением угловой скорости вращательного движения спутника-гиростата п(о) и средним орбитальным движением с, т. е. Оо _ п(о): с _ £1 : к2 (£1 и £2 — взаимно простые целые числа). Период такого решения (по параметру И) То _ 2пк2.

Для исследования вопроса о существовании периодических решений системы (23) при значении л ^ о воспользуемся стандартным достаточным условием Пуанкаре, которое сформулируем в виде следующей теоремы [2— 4].

Теорема 1. Система (23) допускает при малых значениях л ^ О голоморфные по малому параметру л То-периодические решения, близкие к порождающему решению (решение (25) при Оо _ к2), если параметры порождающего решения удовлетворяют условиям

д 2H, dG2

0 * 0; (26)

dHi = dHi = — = 0, (27)

д/о dgo dLo

Hes (Hi)

lo, go, L

* o,

(28)

Hi = -1 f Hi Lo, Go, - = h + go, lo, h

л ■> ir^,

To. f

To

o

dh;

(29)

Hes( F) \x — определитель Гессе, вычисленный от функции F по

величинам, стоящим справа от вертикальной черты.

Вычислив интеграл (29), получим:

а) при Go = ±1 (n(0): а = 1: ±l)

Hi = L - жо У Go2 - LO sin в cos lo + Lo cos в} - Uo.o - u2. ± 2 cos 2 go;

(30)

б) при G0 = ±2 (n(0): а = 2 : ±l)

Hi = L2 - Ж0 {V G2 - L0 sin ecos I0 + L0Cos в} - wo.o - и.± 2 cos go;

(31)

в) при G0 = кJk2, исключая пп. а), б) и к1 = 0 (n(0) = 0)

Hi = -2 - ж0 {^G2 - LL0 sin ecos l0 + L0cose} - Щ.0. (32)

Условие (26) выполняется при любом значении G0. Для случая в) условие (28) заведомо нарушается. Следовательно, с помощью теоремы 1 можно исследовать периодические решения системы (23) только для резонансов k1 : k2 = 1: ±1 и 2: ±1.

Подставим выражения (30) и (31) в уравнения (27) и с учетом неравенства (28) получим следующие решения:

А. При G0 = ±1 (n(0): а = 1: ±1)

п 3п

1) Р0, 00 и в произвольны; I0 = 0, п; go = 0, —, п, —;

ssin2#o

3 (6 -a) cos2 po - 6sacos do + 2 - 3a

16sin (0o -aß)

п 3п п

2) po и ^o произвольны; lo = o, п; go = o, —, п, —; ß = 9o = —;

3) в0 и ж0 произвольны; l0 = 0, п; g0 = 0, п;

S

cosро = -3±в=во + 0,5(1+ с)(п-2во).

Б. При G0 = ±2 (п(0): ® = 2: ±l)

4) ро, во и в произвольны; 1о = о, п; go = о, п;

ssin2e0 (l3 + 9cos2 р0) - 6cos2e0 (l + scos р0) sin росг Жо = 16 sin (во -сф) '

п

5) ро, g0 и в произвольны; 1о = о, п; во = —;

3sin ро (1 + scos р0) cos g0

жо = ;

0 8cose

6) жо и Ро произвольны; 1о = 0, п; g0 = -Л, ^ в = во =

7) ж0 и р0 произвольны; l0 = 0, п; g0 = 0, п;

6(s + cosp0)sinрос _ _ л ч/ ч

1в2во = \ /о;2 ; в=во + 0,5(1 + с )(п-2во). 13 + 9cos2 р0

В пп. 1—7 использованы следующие обозначения: s = 1 при G0 = 1 и 2; s = -1 при G0 = -1 и -2, ст1 = cos l0 = ±1, с = ±1 (при п 3п

G0 = ±1 и g0 = 0, —, п, — значения с = cos2g0 = ±1; при G0 = ±2

и g0 = 0, п значения с = cosg0 = ±1). В соответствии с неравенством (28) некоторые из произвольных порождающих значений в решениях 1— 7 исключаются.

Итак, доказано существование и найдены порождающие семи различного типа семейств периодических решений системы (23). В соответствии с формулами (24) периодические решения (25), вообще говоря, есть условно-периодические решения исходной системы (16). В решениях исходной системы, соответствующих найденным периодическим решениями системы (23), произвольно выбираются значения: t0, ко, ро, 6о, во — для решений типа 1), 4); Í0, ко, ро, 6о, во — для решений типа 3), 5); t0, ко, р0, ж0 — для решений типа 2), 6), 7). Таким образом, найдены три пятипараметрических (при фиксированных значениях жо и в — трехпараметрических) и четыре четырех-параметрических семейства условно-периодических решений системы (16), соответствующих резонансам п(0): а> = 1: ±1 и 2: ±1.

Для исследования периодических решений, соответствующих другим резонансам (см. формулу (32)) воспользуемся одним из достаточных условий существования, полученных в работах [4, 5], которое сформулируем в виде следующей теоремы.

Теорема 2. Система (21)—(23) в особом случае (32) допускает при малых значениях ¡ф о голоморфные по параметру л То-перио-дические решения, близкие к порождающему решению (25), если параметры этого решения удовлетворяют условиям

д Но

dGo2

* o; (33)

=dl!= o; (34)

dio dLo

Л = Hes(H1) i l * o; (35)

v 'lo, Lo

Фр= o, (36)

dgo

дф(-1) , d2H¿

+—т-^ * o. (37)

dgo dgo

Алгоритм построения функций фД 1) приведен в работе [5], индекс i (i е N) — порядковый номер приближения, в котором появляются «критические» для исследуемого резонанса члены. Учитывая структуру функций Hi и рекуррентную формулу для вычисления функций

Фд 1), можно доказать, что достаточные условия существования периодических решений (36), (37) для резонансов n(0) : а = г : p (r е Z, p е N, rp Ф 0) принимают вид

F(p,r) (po, 0o)sinkpgo = 0; (38)

kpF(p,r) (po, 0o) cos kpgo Ф 0, (39)

где при четном rk — некоторое натуральное число, при нечетном rk — четное.

Условие (33) выполняется тождественно. Подставим выражение (32) в уравнения (34) и решим их вместе с уравнением (38) с учетом неравенств (35) и (39). Получим следующее решение:

lo = 0, П; go =np (s = o, ...,2kp-1); Po, 00, и в — произвольные величины;

ssin200 (9cos2 p0 + 4G0 - з)

Жо =-8G I • (0-• (40)

8 |Go|sin (По -oiß)

Естественно, что значения параметров, при которых нарушается условие (35) и F(p, r) (po, По) = 0, исключаются. Соответствующие неравенства можно представить в явном виде, например выражение (35)

А = ж2 sin2 ß ctg2 00 (l + tg3 00 ctg ß)* 0.

Таким образом, можно утверждать, что система (23) имеет бесконечное количество периодических решений (порождающие определяются решениями (40) и G0 = r/p), которые, согласно формулам (24), есть условно-периодические решения исходной системы (16). В решениях исходной системы произвольно выбираются значения fo, h),

Р0, 00, ß и G0 = r/p.

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

СПИСОК ЛИТЕРАТУРЫ

1. Жуковский Н. Е. О движении твердого тела, имеющего полости, наполненные однородной капельной жидкостью: собр. соч. в 6 т. М.; Л.: ОГИЗ, 1949. Т. 2. С. 152—309.

2. Пуанкаре А. Избранные труды. Т. 1. М.: Наука, 1971.

3. Козлов В. В. Методы качественного анализа в динамике твердого тела. М.: МГУ, 1984.

4. Баркин Ю. В., Панкратов А. А. О характеристических показателях периодических решений гамильтоновых систем // Известия АН СССР. МТТ, 1987. № 2. С. 25—33.

5. Баркин Ю. В., Панкратов А. А. Периодические решения гамильтоно-вых систем в некоторых случаях вырождения. ПММ, 1987. Т. 51. Вып. 2. С. 29—41.

6. Белецкий В. В. Движение искусственного спутника относительно центра масс. М.: Наука, 1965.

7. Беленький И. М. Введение в аналитическую механику. М.: Высш. шк., 1964.

Статья поступила в редакцию 14.09.2012