CC BY
52
'■Ш
МАТЕМАТИКА
УДК 517.928 ББК 22.161.6
ГАМИЛЬТОНОВЫ УКОРОЧЕНИЯ ОГРАНИЧЕННОЙ ЗАДАЧИ ТРЕХ ТЕЛ
Н.В. Батхина, И.А. Тимофеева
Целью данной работы является определение методами степенной геометрии предельных случаев для пространственной ограниченной задачи трех тел и пространственной задачи Хилла, а также сравнение полученных результатов с результатами работы [3].
Ключевые слова: степенная геометрия, системы Гамильтона, теория возмуще-
Введение
Данная работа посвящена исследованию гамильтоновых укорочений пространственной ограниченной задачи трех тел (далее — ОЗТТ) и пространственной задачи Хилла методами степенной геометрии.
Основная концепция степенной геометрии заключается в изучении нелинейных задач не в исходных координатах, а в логарифмах этих координат. Тогда многим свойствам и соотношениям, нелинейным в исходных координатах, можно поставить в соответствие определенные линейные соотношения. Алгоритмы степенной геометрии основаны на этих линейных соотношениях: они позволяют упрощать уравнения, разрешать их особенности, выделять их первые приближения и находить либо их решения, либо асимптотики этих решений. После одного шага таких упрощающих преобразований можно получить степенные решения или степенные асимптотики решений. После кратного применения этих алгоритмов можно получить решения или их асимптотики, содержащие кратные логарифмы и экспоненты. Этот подход позволяет также вычислять асимптотические и локальные разложения решений. Алгоритмы степенной геометрии применимы к уравнениям разных типов: алгебраическим, обыкновенным дифференциальным и в частных производных, а также — к системам таких уравнений.
Над разработкой методов степенной геометрии работали И. Ньютон, М.М. Вайн-берг, В.А. Треногин, С. Лефшец, П.Г. Айзенгендлер, А.Д. Брюно (см. обзор [9], книги [4; 5], статью [2] и указанную в этих работах библиографию). Применение методов
степенной геометрии к небесно-механическим задачам выполнено в работах А.Д. Брю-но [3;4; 11].
В 2007 г. в работе А.Д. Брюно и В.П. Варина [3] с помощью методов степенной геометрии найдены все предельные случаи плоской ограниченной задачи трех тел при стремлении массового параметра ¡л к нулю: задача двух тел, задача Хилла, промежуточная задача Энона и основная предельная задача. В каждой из них выделяются решения, которые являются пределами периодических решений ограниченной задачи трех тел при ¡л —> 0, и пределы семейств периодических решений (называемые порождающими семействами).
Методы степенной геометрии могут быть эффективно применены к исследованию гамильтоновых систем, при этом достаточно выполнять основные операции степенной геометрии над функцией Гамильтона, а не над системой связанных с ней канонических уравнений. Напомним основные идеи степенной геометрии, следуя [4, гл. I, II, IV].
Пусть z = (х, у) G М2га — вектор фазового пространства, х G М™ — координаты, у G R" — канонически сопряженные импульсы, а Я(z) — функция Гамильтона, аналитическая в точке z0. Тогда она разлагается в ряд Тейлора
Я(2) = ^ Як(2 - z0)R (1)
Res'
где R G 1?™ — векторный показатель степени, S — некоторое дискретное множество в 1?™ — носитель суммы (1), Яр G М, сходящийся в некоторой окрестности U точки z = z0. Носителю S разложения (1) соответствует внешняя выпуклая оболочка
Г (Я) =f conv/S, которую называют многогранником Ньютона. Граница <9Г многогранника Ньютона Г(Я) состоит из граней Г^ размерности d, а каждой грани Г^ соответствует граничное множество Sjd\ касательный конус Т^ в М2га и нормальный конус U^ в сопряженном пространстве М2'\ Каждой грани Г^ соответствует свое укорочение Н^ гамильтониана Я, определяемое следующим образом:
Hf = noRGSf.
R
Оказывается (см. [4, гл. IV]), что не каждое укорочение Н^ функции гамильтона Я будет гамильтоновым. В работе [11] было предложено находить укорочения функции Гамильтона Я так, чтобы они все были гамильтоновы. А именно, каждому векторному показателю R = (P,Q) G Р = (pi,... ,р„), Q = (qi,---,qn) ставится в соответствие точка — проекция 7г(Р, Q) на двумерную решетку Z2, определяемую по правилу:
7г(Р, Q) =f (p,q) G Z2, где p = Yn=iPi> Q = Yn=iQi- На плоскости R2 для проекции 7tS носителя S строится выпуклая оболочка 7гГ и для каждой грани 7 ее границы <97гГ находится соответствующее укорочение Я7 функции Гамильтона Я.
Если гамильтониан Я зависит дополнительно от вектора малых параметров М = (¿¿1,... ,/j,i), то разложение (1) производится не только по фазовым координатам z, но и по параметрам М:
Я = ^ЯРдтхр yQMT
где носитель S = {(P,Q,T) G Z2ra+/ : ЯРдТ ф 0}. В этом случае проектирование 7г носителя S выполняется в трехмерную целочисленную решетку Z3 по правилу
7г(Р, Q, Т) =f (р, q, t), где р и q определены выше, a t = ^¿=1U-
Гамильтоновы укорочения Н^ — это различные первые приближения гамильтониана Н, справедливые в различных областях фазового пространства М2га или фазового пространства М2,г+/, расширенного пространством параметров. Выбор этих областей определяется соответствующим конусом задачи К € М2га+/, представляющим собой открытый выпуклый конус.
Гамильтониан Н можно рассматривать как возмущение некоторой более просто устроенной гамильтоновои системы, определяемой соответствующим укорочением Н - . Эти укорочения могут быть устроены существенно проще исходного гамильтониана, и, возможно, удастся проинтегрировать соответствующие этим укорочениям системы канонических уравнений и получить начальные приближенные решения полной задачи.
1. Десингуляризация пространственной ограниченной задачи трех тел
Рассмотрим гамильтониан пространственной ОЗТТ, канонические уравнения которого описывают движение частицы пренебрежимо малой массы под действием сил гравитации двух тел, которые равномерно вращаются вокруг общего центра масс по круговым орбитам. Система координат выбрана так, что тяготеющие тела массами 1 — р и р (р < 1/2) расположены в точках (—¿/,,0,0) и (1 — ¿/,,0,0) соответственно. Тогда гамильтониан Н в безразмерных координатах имеет вид
я(х>У) = \{У\ + Уъ + Уз) + Х'2У1 ~ хгу2 ~ -—- - —, (2)
2 р 1 р2
где
Р1 = \/+ И)2 + х'22 + хз-> Р-2 = \/- (1 - н))2 + Х2 + хз-
Для того чтобы при малых р, найти все первые приближения ОЗТТ вблизи меньшего из двух тяготеющих тел, нужно ввести локальные координаты
& = X! + р - 1, $2 = Х2, 6 = Хз, т = Уъ 42 = У2 +Р~1, Щ = Уз
и разложить функцию Гамильтона (2) по этим координатам. После разложения функции 1 /V(и + I)2 + 42 + £з в РЯД Маклорена функция Гамильтона (2) примет вид
2 1 2 м 2- ^а+ж+ш
- (1 - б+£ -1 -1 - & е2)) + ^ (1 - б+£ -1 -1 - & е2)) •
После упрощения получим следующую функцию Гамильтона:
н = + 2р - ^ +1 № + п2 + %2) + Ьт - Ът - & + | + |
+/ & , & е2) +»(е? -1 -1 - & ?з2) -
1
(3)
л/О+Ж+1:
2
где / — сходящийся степенной ряд, не содержащий членов порядка меньше трех. Носитель (¿>1) ряда в правой части равенства (3) состоит из точек
Я = (Г1, Г 2, Гз, г4, г5, г в, г7) = (огс1 , огс1 £2, ОГС1 , огс11]1, огс11}2, огс1 ?7з, огс1 ц) : (0,0,0,2,0,0,0), (0,0,0,0,2,0,0), (0,0,0,0,0,2,0), (-1,0,0,0,0,0,1), (0,-1,0,0,0,0,1), (0,0,-1,0,0,0,1), (0,1,0,1,0,0,0), (1,0,0,0,1,0,0), (2,0,0,0,0,0,0), (0,2,0,0,0,0,0), (0,0,2,0,0,0,0), (2,0,0,0,0,0,1), (0,2,0,0,0,0,1), (0,0,2,0,0,0,1), (к, 21, 2 т, 0, 0, 0, 0), (к, 21, 2т, 0, 0, 0,1),
где к, I ^ 0, к + 2/ + 2т ) 5 и плоскости, соединяющей точки (—1,0,0,0,0,0,1), (О, —1, 0, 0, 0, 0,1), (0, 0, —1, 0, 0, 0,1). Эта плоскость — носитель члена ц,/\/х'{ + +
Конус задачи будем выбирать так, чтобы укорочения гамильтониана (3) соответствовали малым значениям параметра ¡л и описывали движение тела в окрестности начала координат. Таким образом, конус задачи будет
К = {]¥ е М* : и>1 < 0,и>2 < 0,и>3 < 0,ги7 < 0},
так как ->■ 0.
Построим проекцию
тгЯ = В" = (р, д, в) е
где р = Г1 + г2 + г3, д = г4 + г5 + г6 и в = г7.
Множество Б" этих точек состоит из точек (0,2,0), (2,0,0), (А:, 0,0), (—1,0,1), (1,1,0), (2,0,1), (А:, 0,1), А: = 5, 6, 7, ...
Замыкание выпуклой оболочки множества Б" — это многогранник Г С М3, состоящий из граней ребер Г^ и вершин Проекция исходного конуса задачи К есть конус задачи
К" = {\¥" е М* : «'1 < 0,ги3 < 0} .
Многогранник Г представляет собой полубесконечную трехгранную призму с косым основанием, которая имеет четыре грани, шесть ребер и три вершины (см. рис 1).
Грань служит косым основанием призмы Г, она содержит вершины (0,2,0), (2,0,0), (-1,0,1) и точку (1,1,0). Ее вектор нормали = -(1,1,3) € К". Ей соответствует укороченная функция Гамильтона
я!2) = \ (■ч1 + п! + %2) + Ьт - Ът - £ +1 +1 - 2,
2 2 2 + Й
которая при ¡л = 1 описывает пространственную задачу Хилла.
Рис. 1. Проекция многогранника Ньютона гамильтониана (3) ограниченной задачи трех тел
Для малых р > 0 вблизи меньшего тела имеется сингулярное возмущение случая р = 0. При р —> 0 объекты ограниченной задачи в координатах стремятся к пре-
дельным положениям. При этом для некоторых объектов предельные положения имеют меньшую размерность, чем допредельные. В этом и проявляются сингулярные возмущения. Чтобы устранить сингулярность возмущения при р > 0, надо вводить новые (зависящие от р) координаты, в которых размерность предельного объекта равна размерности допредельного. При этом в качестве предельной получается задача, соответствующая первому приближению сингулярного возмущения. Введение таких координат в [3] названо десингуляризацией.
Степенное преобразование
& = У
1/3 г.
Чг
г = 1,2,3
приводит пространственную задачу Хилла к системе с функцией Гамильтона
1 й
Ннт = - {гЦ + /)22 + /)32) + Ып ~ - + ^
1
I? + й + 42
(4)
Задача Хилла является неинтегрируемой, но имеет многочисленные применения (см., например, [1]), а ее гамильтониан исследован методами степенной геометрии в п. 2.
Грань Ц содержит вершины (0, 2, 0), (2, 0, 0), (1,1, 0) и точки (к, 0, 0), к = 5, 6,____
Ее вектор нормали N¡1 = (0,0,-1) € К". Ей соответствует укороченная функция Гамильтона
Н.
которая описывает задачу, занимающую промежуточное место между интегрируемои задачей Энона (6) и неинтегрируемой задачей Хилла (4).
Остальные две грани имеют нормальные векторы (0, —1,0) и (0,1, 2), лежащие вне конуса задачи К", поэтому соответствующие укорочения гамильтониана не рассматриваются.
Рассмотрим ребра. Из шести ребер одно несобственное. Оно проходит через точку (0, 2, 0) параллельно вектору (1, 0, 0). На трех ребрах q = 0, то есть для них укороченная функция Гамильтона не зависит от импульсов г = 1,2,3, и у решений соответствующей системы Гамильтона = const. Остаются два ребра Г^ и Г^.
Ребро Т[1] содержит точки (0,2,0) и (—1,0,1) множества S". Соответствующая укороченная функция Гамильтона
которая описывает пространственную задачу двух тел в неподвижной системе координат. Эта задача интегрируема и все ее финитные движения есть периодические орбиты (см., например, [8]).
(О)
Пусть Гд — грань, проходящая через точки (0, 2, 0), (2, 0,1), (—1,0,1). Ее внешняя нормаль Nq = (0,1,2), а вектор нормали к грани есть N" = —(1,1,3). Поскольку ребро Г^ — пересечение граней Гд2^ и то его нормальный конус состоит из векторов N = (п\, п-2, /гз) = cxNq + ¡3N", где а, /3 > 0. Положим а = /3 = 1, тогда получаем вектор N = —(1,0,1), лежащий в нормальном конусе ребра Ему соответствует степенное преобразование
6 = Пг = Н'Пг, t = pi г = 1,2,3, которое приводит к системе с функцией Гамильтона
n = \ ivi + + " , \
2 V^i + й + ^
Эта функция отличается от гамильтониана (5) только коэффициентом р = 1.
Ребро Г^ содержит точки (0,2,0), (2,0,0), (1,1,0) множества S". Ему соответствует укороченная функция Гамильтона
Щ1] = \ (vl + vt + ni) + bm - Cm - £ +1 + f, (6)
которая описывает промежуточную задачу Энона. Название задачи было предложено в работе [3], поскольку ее плоский вариант был подробно исследован М. Эноном (М. Нёпоп) в работе [12]. Эта задача рассматривалась многими исследователями, в частности, она применялась для построения и исследования квазиспутниковых орбит (см. [6; 7; 10]).
Укорочения, соответствующие вершинам, не являются содержательными, поскольку содержат либо члены, зависящие только от координат, либо члены, зависящие только от импульсов.
Итак, как и в случае плоской ОЗТТ [3], в окрестности меньшего из двух тел
имеются три различные десингуляризации, соответствующие грани Г^ и ребрам Г^ и r(i)
2. Десингуляризация пространственной задачи Хилла
Рассмотрим гамильтониан пространственной задачи Хилла, который имеет одну особенность в начале координат.
Hum = \ (у2 + У 2 + Уз)--/ 2 1 2 I + Х'2У1 ~ Xiy'2 ~ \ (2;Г1 ~ - 4) ■ (?)
Перейдем в пространство показателей степеней и построим носитель Si правой части равенства (7). Он состоит из точек
R= (ri, Г2, Г3, Г4, Г5, г в) = (ord^i,ord^2,ord^3,ordyi,ordy2)Ordy3) : (0,0,0,2,0,0), (0,0,0,0,2,0), (0,0,0,0,0,2), (-1,0,0,0,0,0), (0,-1,0,0,0,0), (0,0,-1,0,0,0), (0,1,0,1,0,0), (1,0,0,0,1,0), (2,0,0,0,0,0), (0,2,0,0,0,0), (0,0,2,0,0,0)
и из плоскости, соединяющей точки (—1,0,0,0,0,0), (0,-1,0,0,0,0), (0,0,-1,0,0,0). Эта плоскость — носитель члена \/\Jx\ + х22 + #3. Рассмотрим два конуса задачи. Первый конус задачи Ki определим так, чтобы укорочения гамильтониана (7) описывали движение тела в окрестности начала координат:
Ki = {W G М® : wi < 0, w2 < 0,гу3 < 0},
так как ^1,^2,^3 —> 0. Второй конус задачи Кг определим так, чтобы укорочения гамильтониана (7) описывали движение тела вдали от начала координат:
К2 = {W G М® : wi > 0, w2 > 0,гу3 > 0},
так как . х,- хл —)■ оо. Построим проекцию
7Г R = R" = (р, q) еж2,
где р = ri + г2 + г3 и q = r4 + r5 + r6.
Множество S" этих точек R" состоит из (0,2), (—1,0), (1,1), (2,0). Замыкание выпуклой оболочки множества S" — это многоугольник Г С М2. Граница <9Г многоугольника Г состоит из ребер Г^ и вершин Г^ (см. рис. 2). Каждо-
-nid) - Tj{d)
му такому элементу соответствует укороченный гамильтониан Н - , являющийся суммой тех членов (7), точки которых R" G Г^. Укороченные функции гамильтониана Н^ — это различные первые приближения функции (7), справедливые в разных областях фазового пространства М6. Конусы задачи есть К'/ = {W" G М2 : wi < 0} и Щ = {W" G М2 : Wi > 0}.
Рис. 2. Проекция многогранника Ньютона гамильтониана (7) пространственной задачи
Хилла
Ребро Г^ содержит вершины (2,0), (0,2) и точку (1,1). Его вектор нормали N'1 = (1,1) € К'з. Ему соответствует укороченная функция Гамильтона
Я,
Нёпоп
а= н[1] = ^ (у\ + у\ + у\) + х2уг - Хгу2 - ^ (2х{ х'1)
которая описывает промежуточную задачу Энона. Эта задача занимает промежуточное место между задачей двух тел и задачей Хилла. Степенное преобразование
Хг — у |С| Аг, Уг — у |С|1'г, С 0, — 1,2,3, приводит гамильтониан (7) к виду
£
ЯНШ(Х, У) = Янёпоп(Х, У) - +
где е = |С|"3/2. Следовательно, при больших значениях \С\ задачу Хилла можно рассматривать как возмущение интегрируемой промежуточной задачи Энона. Такой прием в случае плоской задачи Хилла был использован Л. Перко (Ь. Регко) для доказательства существования семейств симметричных периодических решений задачи Хилла (см. [13; 14]).
Ребро Г^ содержит вершины (—1,0) и (2,0). Его вектор нормали Щ = (0,1) лежит вне конусов задачи К7/ 2, поэтому соответствующее укорочение гамильтониана (7) не рассматривается.
Ребро Гд^ содержит вершины (—1, 0) и (0, 2). Его вектор нормали Л^' = (—2,1) е К'/. Ему соответствует укороченная функция Гамильтона
ЯкеР1еГ = Щ1] = I (у! + у! + у!) - 1
2 \/ х'1 + Х22 + х'1 '
которая описывает задачу двух тел (задачу Кеплера) в неподвижной системе координат. Степенное преобразование
хг = еХг, Уг = £~1/2¥г, Ь = £"3/25, е < 1, г = 1, 2, 3, приводит гамильтониан (7) к виду
#нш(Х, У) = #Кер1ег(Х, У) + ^(ВД - ВД) - ^ (2Х2 - XI - XI) .
Следовательно, при малых е задачу Хилла можно рассматривать как возмущение интегрируемой задачи двух тел.
Укорочения, соответствующие вершинам не рассматриваются в силу тех же причин, что указаны в п. 1.
Заключение
Методами степенной геометрии [4] были получены четыре предельных случая пространственной ограниченной задачи трех тел и два предельных случая пространственной задачи Хилла. Показано, как с помощью степенных преобразований представить неин-тегрируемую задачу в виде возмущения интегрируемой.
СПИСОК ЛИТЕРАТУРЫ
1. Батхин, А. Б. Задача Хилла / А. Б. Батхин, Н. В. Батхина. — Волгоград : Волгогр. науч. изд-во, 2009. - 200 с.
2. Брюно, А. Д. Асимптотики и разложения решений обыкновенного дифференциального уравнения / А. Д. Брюно // УМН. - 2004. - Т. 59. - Вып. 3. - С. 31-80.
3. Брюно, А. Д. Периодические решения ограниченной задачи трех тел при малом отношении масс / А. Д. Брюно, В. П. Варин // Прикладная математика и механика. — 2007. - Т. 71. - Вып. 6. - С. 1034-1066.
4. Брюно, А. Д. Степенная геометрия в алгебраических и дифференциальных уравнениях / А. Д. Брюно. - М. : Наука, 1998. - 288 с.
5. Вайнберг, М. М. Теория ветвления решений нелинейных уравнений / М. М. Вайн-берг, В. А. Треногин. - М. : Наука, 1969. - 528 с.
6. Коган, А. Ю. Далекие спутниковые орбиты в ограниченной круговой задаче трех тел / А. Ю. Коган // Космические исследования. — 1988. — Т. XXVI. — Вып. 6. — С. 813-818.
7. Тучин, А. Г. Проектирование квазисинхронных орбит космического аппарата вокруг Фобоса для решения задачи посадки на его поверхность / А. Г. Тучин // Космические исследования. — 2008. — Т. 46. — № 6. — С. 536-546.
8. Уинтнер, А. Аналитические основы небесной механики / А. Уинтнер. — М. : Наука, 1967. - 532 с.
9. Чеботарев, Н. Г. Многоугольник Ньютона и его роль в современном развитии математики / Н. Г. Чеботарев // Исаак Ньютон. — М. ; JI. : АН СССР, 1943. — С. 99-126.
10. Benest, D. Libration effects for retrograde satellites in the restricted three-body problem / D. Benest // Cel. Mech. - 1976. - V. 13. - № 2. - P. 203-215.
11. Bruno, A. D. On periodic flybys of the Moon / A. D. Bruno // Celestial Mechanics. — 1981. - Vol. 24. - P. 255-268.
12. Henon, M. Numerical Exploration of the Restricted Problem. V. Hill's Case: Periodic Orbits and Their Stability / M. Henon // Astron. & Astrophys. - 1969. - V. 1. -P. 223-238.
13. Perko, L. Families of symmetric periodic solutions of Hill's problem I: First species periodic solutions for С < -1 / L. Perko 11 Amer. J. Math. — 1981. — V. 104. — № 2. - P. 321-352.
14. Perko, L. Families of symmetric periodic solutions of Hill's problem I: Second species periodic solutions for С < -1 / L. Perko 11 Amer. J. Math. — 1981. — V. 104. — № 2. - P. 353-397.
HAMILTON IAN RESTRICTIONS OF THE THREE-BODY PROBLEM
N.V. Batkhina, I. A. Timofeeva
The purpose of this study is to determine of limiting cases for the spatial restricted three-body problem and the spatial Hill's problem by Power Geometry methods, and compare our results with those of [3].
Key words: power geometry, Hamiltonian systems, perturbation theory.
