Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы качественного анализа в динамике твердого тела

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
7.68 Mб
Скачать

4. Разложение возмущающей функции

51

Значит и>1 = K V W , где

 

 

 

 

к - 1 (р; A ,B ,C )

=

- ±

- i

/ И

1

 

 

 

 

 

7

й 2Z_| cos21 )]

 

 

 

 

 

 

sin2l

>

0.

 

 

 

 

 

cos2l '

 

Функция 7 зависит

лишь от

р, Л,

В, С;

следовательно,

с?7/дЗ' = 0. Поэтому

 

 

 

 

 

 

 

7)

-^(р)

c?a;i

к (р )

д у

 

det

2- /^

 

 

 

2yfW др'

 

д(& , Р)

О

с? 7

 

 

 

 

др

 

 

 

 

 

 

 

 

 

 

Учитывая результаты предыдущих вычислений, получаем окончательно

г =

1 d j

d l2

/|

Чтобы решить вопрос о знаке гессиана, осталось применить лемму 4, которая утверждает, что производная д'у/др > 0 (< 0), когда 1< р < 1/В (1/В < р < 1/С).

Замечание. Теорема 3 является уточнением соответствую­ щего результата, полученного В. И. Арнольдом [9, дополнение], ко­

торый доказал, что в общем случае Г ^0.

§4. Разложение возмущающей функции

Спомощью формул сферической тригонометрии можно показать [26], что в специальных канонических переменных L, G, Н, I, g, h возмущающая функция °U имеет вид

а), _ х

, z

(4.1)

— f h

+ 772

+ 77з»

52 Глава 2

Tl = f\/1 “ (§)2sinZ+§\/1 “ ( f ) 2sin/cos^+

+ 0 . - ( | т ) COS /s i n g ,

T2 = f\/1“(§) cosZ+§^/1 _(f)aCOslcOSg~

~ \ 1 ~ ( § ) sin/sin^’

Тз = ^ У ! - ( § ) 2 ]/ l - ( § ) 2cosg,

где (ж, у, z) — координаты центра тяжести в главных осях инерции, г = у/х2 + у2 + z2 — расстояние от центра тяжести до точки подвеса.

В переменных

действие-угол /| / 2/3 <pi <р2<рз функция

I213 Ц>1 № ‘Рз)

не зависит от ipz и 2тг-периодична по уч

и ip2. Согласно результатам работ [22, 18], разложение этой функции в двойной ряд Фурье по переменным (р± и ц>2 можно представить следующим образом:

/ х/г

cos S = h

°U = (sin^sin^, sin<Scos<^>2, cosS)S I y/r

V /r.

V

(4.2)

 

Здесь элементы квадратной матрицы 5 = ЦауЦ третьего по­ рядка не зависят от ср2. Выпишем разложения в ряды Фурье по угловой переменной tp\\

S 1 1

2тг

 

«»+ i/2( l - g 2n+i) eh J

— — Г

у/х2 + Л2 ^

sin(2n+l)^>i,

 

к

1—2g2n+1 ch.2cr-\-q4n+2

S12 — —

1 + x 2

^ g«+i/2(l + д2п+1) ch<5

Л2 + x 2

cos(2n+l)^>i,

 

 

1+2q2n+1 ch 2o-+g4n+2

«13 = 2тг

 

(1 - q2n) cho-

X

An sin2ny>i,

 

К

V x 2 + А2 h

1 - 2q2n ch2a + q

 

 

 

 

 

 

 

 

 

Исторический очерк

53

,

_

2тг

 

 

 

 

 

 

„n+l/2

 

д2и+1^

 

 

 

 

 

-

_____

 

 

cos(2n +

« 21

____________ V

 

 

 

 

к

 

\/х2 + Л2

1 -

2q2n+1 ch2(T + qin+2

 

 

 

27Г

/ 1 +

X

 

 

^

«П+

 

 

2™+1) sh<r

 

 

 

 

 

«

 

1/2(1 _____________ sin(2n +

22 -

K

t l A

V

n=0

i _| 2q2n+1 ch 2(7+ q4n+2

,

_

2TT

 

 

 

 

 

 

 

 

 

 

qn{\ + q2n) sh(7

 

 

 

 

 

 

 

 

 

 

 

An C O S 2i n<^ij,

«23 -

IF

Vx2 2

 

 

 

 

n=l

 

 

 

к

 

{

 

 

 

1 2g2™ch2(T+g

 

 

 

 

 

 

 

 

4sh(5^" ^

 

 

 

2TT

 

 

 

 

 

 

n + l / 2

 

 

 

«31

 

 

 

 

 

E

4

 

— cos(2n + l)y>i,

 

К

 

 

 

 

 

 

„2n+l

 

 

 

 

 

yjx2 + A2 n=0

1 + 2q

 

 

 

 

 

 

27Г

/ 1 + X 2

 

 

,n + l / 2

sin(2n + l)y>i,

« 3 2 = - К

A2+ x2 n = 0n 1 u

 

n

 

 

 

 

2TT

 

X

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

r f ^ cos2n^ );

« 3 3

=

^

 

V x 2 + Л2 (i + E

~FT

 

 

 

 

 

 

 

 

 

 

 

 

x 2 = C(A

B)

 

Д2 = ^

2С^

 

J2

exp

 

 

Л ( В - С ) ’

 

 

 

 

I2 -

2

 

К ' = К ( ^ 1 - Л 2),

 

 

a = 2K ^ (a rctg | ,

V l - Л 2) .

Здесь К(Л) — полный эллиптический интеграл первого рода с модулем Л, F — эллиптический интеграл первого рода.

Исторический очерк

Задача о вращении тяжелого твердого тела с неподвиж­

ной точкой как возмущение интегрируемого

случая

Эйле­

ра Пуансо впервые поставлена А. Пуанкаре в

пятой

главе

«Новых методов небесной механики».

 

 

Пуанкаре принадлежит важное замечание

о том,

что

в некоторых канонических переменных I, <р гамильтониан сво­ бодного вращения твердого тела имеет вид .'У (/i, I2). Им же введена функция а(2/У //|; А, В, С ); отношения а/2тг суть числа вращения (опять-таки определенные впервые Пуанкаре) на двумерных торах интегрируемого случая Эйлера Пуансо. Пуанкаре первым указал вид разложения возмущающей функ­ ции в кратный ряд Фурье по угловым переменным ipi, q>2. Ссы-

54

Глава 2

лаясь при этом на некоторые промежуточные формулы как на известные, А. Пуанкаре, наверное, прежде всего имел в виду ре­ зультаты работы Якоби [22], в которой указаны разложения направляющих косинусов свободно вращающегося тела в три­ гонометрические ряды по величинам и>Д и из-Я, где u>i и и>2 зависят только от постоянных первых интегралов. Так как возмущающая функция сводится к линейной комбинации на­ правляющих косинусов, то разложения Якоби после замены ui[t и u>2t на углы (fix и g>2 дают как раз ряд Фурье этой функции

впеременных ipi и tp2-

Внекоторых работах (например, [26, 33, 34]) специальные канонические переменные L, G, Н , /, g, h несправедливо назы­ вают переменными Депри. Это связано, вероятно, с тем, что

водной из работ Депри [15], где вводятся эти переменные, от­ сутствуют ссылки на другие источники. Однако специальные канонические переменные давно применялись в небесной меха­ нике при анализе вращательного движения небесных тел (см., например, трактат А. Андуайе [16]).

Используя разделение специальных канонических перемен­ ных в функции Гамильтона задачи Эйлера-Пуансо, Ю.А.Садов получил явные выражения для переменных действие-угол [18]. Отметим, что формулы, определяющие переменные дейст­ вие, были найдены иным способом в квантовой механике уже

вначале XX в., в связи с исследованием спектров многоатом­ ных молекул [19]. Дело в том, что свободно вращающееся твердое тело является в «классической» квантовой механике простейшей моделью «невозбужденной» молекулы. Как извест­ но, переменные действие играют определяющую роль в услови­ ях квантования Бора - Зоммерфельда.

Отметим в заключение, что вычисления Ю. А. Садова ко­ эффициентов Фурье разложения возмущающей функции повто­ ряют в современных обозначениях вычисления Якоби из уже упоминавшейся работы [22].

Глава III

Неинтегрируемость задачи о вращении несимметричного тяжелого твердого тела вокруг неподвижной точки

В соответствии со стандартными обозначениями класси­ ческой теории возмущений функцию Гамильтона рассматри­ ваемой задачи обозначим

Ж Ж$ + fi2?t1

(здесь Ж0 = 5", Ж\ =°U}. Канонические уравнения с гамиль­ тонианом Ж имеют циклический интеграл — интеграл пло­ щадей (в переменных действие-угол I , он равен / 3). Фикси­ руя его постоянную (13 = Г®), сведем задачу о вращении тя­ желого твердого тела с неподвижной точкой к системе с двумя степенями свободы, которую будем называть приведенной. На протяжении всей главы твердое тело предполагается несим­ метричным, т. е. А > В > С.

§ 1. С тр уктур а векового м нож ества

Из формулы (4.2) главы II следует, что разложение воз­ мущающей функции Ж\ в двойной ряд Фурье имеет вид

00

00

00

—оо

Н т ,

+ Y , Н п ,0eimvi.

—оо

—оо

 

 

( 1.1)

Коэффициенты этого разложения зависят от 11: / 2, / 3. Их лег­ ко вычислить, использовав формулы § 4 главы II.

56

Глава 3

Функция Гамильтона Ж определена и аналитична по пе­ ременным I\, I2 в области

Д0а = Д а П {(Д ,/2):| /з°| < /2}.

В частности, коэффициенты (1.1) аналитичны в Д°. Согласно определению 1 § 1 главы 1 вековым множеством

системы с гамильтонианом Ж§ + цЖ\ называется множество точек I = (Ii, I2) £ Д°, удовлетворяющих условиям:

1)

mwi(I) ±ш 2{1) = 0; Шг = дЖ0/д1{ (г = 1, 2), т £ Z;

2)

Нт,±1{1)ф 0.

Теорема 1. В каждой из четырех связных подобластей области Д° множество 3% состоит из бесконечного числа пря­ мых, проходящих через начало координат. Эти прямые накап­ ливаются у одной из двух прямых линий

2Ж0(11, I2) = Ii/B,

лежащих на границе Д°.

Доказательство.

Из результатов § 2 гл. II следует, что множество точек из области Д° = Д П { |Zg |< / 2}, удовлетворяющих уравнению

ти>1 ± и>2 = 0

(1.2)

при достаточно больших фиксированных т и выборе знака пе­ ред частотой и>2 (когда 1\ > 0, то при т > 0 (< 0) выбирается знак — (+); когда 1\ < 0, то наоборот), суть четыре прямые линии (по одной в каждой связной компоненте области Д°), проходящие через начало координат Ii = 12 = 0. Ниже будет показано, что для бесконечного числа значений т G Z функ­ ции Hm,±i(-0 отличны от нуля на прямых (1.2). Использовав теорему 1 гл. II, легко найти множество предельных точек для £$. Для этого устремим в уравнении (1.2) т к бесконеч­ ности. Тогда

0>2tJT = Ы р ш, А, в , <7)1 -> (X).

Согласно теореме 1 (гл. II) величина р = 2Ж^)1\ стремит­ ся к 1/В (В — средний момент инерции). Итак, множество

§1. Структура векового множества

57

предельных точек для % — две прямые линии

2Ж0(1и h ) =

которые лежат на границе Да и пересекаются в начале коор­ динат.

Осталось показать, что в разложении (1.1) бесконечно много коэффициентов вида Нт>±i отличны от нуля. Рассмот­ рим три случая:

1) X2 + у2 = О, г ф 0; 2) ж2 + у2 ф 0, г = 0; 3) х 2 + у2 ф 0, г ф 0.

Пусть сначала центр тяжести лежит на оси Oz (случай 1). Согласно формулам § 4 гл. II в этом случае т — четные це­ лые числа. Пусть т = 2п, пЕ N — натуральное число. Случай т = —2п, пЕ N рассматривается аналогично. Используя фор­ мулу (4.2) (гл. II), разложения Фурье для S13, S23>5зз> а также формулы

sin^ = ^(е*> -

cos ip = ± (eiv + е " * ),

 

нетрудно установить, что

 

 

Н 2 п , - 1 — Н - 2 п , 1

z 2ж

х

 

4Г К

+ Л2

 

 

 

 

«™(1 + q2n) sh<7qn(1 q2n) ch<r

(1.3)

x sino------------------- ------------------------------

 

1 -

2q2n ch 2a + q4n

 

Так как sin<$ ф 0 (в противном случае |/j| = I2), то коэффи­ циент H2n, - i равен нулю только в том случае, если

(1 —q2n) ch <7—(1 + q2n) sh <7= 0.

(1.4)

Когда п —>• сю, то согласно теореме 1 (гл. II) функция р =

=2Ж0Щ стремится к 1/В. При этом, очевидно,

Л= х,

58 Глава 3

Так как Л — непрерывная функция от р в некоторой окрест­ ности точки 1/В, то существует номер Ni(A, В, С), такой, что при tn2/wi = 2п > N± справедливы неравенства

О < Ai < А < Л2 < 1; Ai, Л2 = const.

При этом, очевидно, |сг| < сг0 и |g| < q0 < 1 (числа Ai, Л2, Сто и до в конечном счете зависят только от А, В, С). Следо­ вательно, существует номер N(A, В, С) (N > Ni) такой, что при 2п > N равенство (1.4) не может выполняться.

Рассмотрим теперь случай 2, когда х2 + у2 ф 0, г = 0. В этом случае числа m нечетные. Согласно формуле (4.2) (гл. II)

Ж\ = |)(sin(5sim£>2Sn + sin(5cos^>2S2i + cos<Ss3i) +

+ | (sinSSin <£>2S12 + sin(5cOS^>2S22 + CO<SS32).

Используя разложения Фурье функций Sij, нетрудно устано­ вить, что для натуральных п

Н;

— Н - 2 п - 1 , 1

х 2тг sin ё-

 

 

2 п + 1 ,

4г К

фус2 + А2

 

 

g™+i/2[(l + q2n+1) sha — (1 q2n+1) cha]

 

 

1 2q2n+1 ch 2(7 + qin+2

+

 

 

1

+ x2

 

 

 

A2 + x 2

 

 

qn+1/2[(i - q2n+1) sha - (1 + q2n+1)ch a ]

 

X

1 + 2g2n+1 ch 2(7+ g4n+2

'

Аналогичная формула

справедлива

для коэффициентов

# 2 п + 1 ,1 = Н - 2 п - 1 , - 1 ~

 

=

0 только в том случае,

Так как sin <5 ф 0, то ff2n + i,-i

если обращаются в нуль его действительная и мнимая части. Значит, когда х ф 0, то

(1 + g2n+1) sho- — (1 — q2n+1) cha = 0,

(1.5)

а когда уф 0, то

 

(1 - q2n+1) sh <т — (I T q2n+1) cha = 0.

( 1.6)

§1. Структура векового множества

59

Эти уравнения имеют вид равенства (1.4). Совершенно анало­ гично можно доказать существование такого N (A, В, С), что соотношения (1.5) и (1.6) не имеют места при 2 п + 1 > N. Следовательно, при достаточно больших п коэффициенты Я 2я+1,1 ф 0.

Общий случай, когда х2+ у 2 ф 0,2 ф 0, очевидно, сводится к двум рассмотренным. ■

Напомним (§1, гл. I), что вековым множеством мы на­ зываем также множество резонансных торов в фазовом про­ странстве невозмущенной задачи, отвечающих значениям пе­ ременных действие I € £$. Опишем это множество, используя специальные канонические переменные L, G, I, g (значение интеграла площадей Н = const зафиксировано).

Рассмотрим двумерное кольцо

К = {(/, L) <ЕR 2 : 0 < I ^ 2тг, \L\^G ,G = G°, g = 0},

секущее трехмерный уровень интеграла модуля момента G = const задачи Эйлера Пуансо (ср. с § 1 гл. II). Траектории невозмущенной системы трансверсальны к кольцу К всюду, кроме границ, которые представляют собой периодические ре­ шения — постоянные вращения вокруг большей оси инерции в противоположных направлениях. Любая точка к кольца К через некоторое время снова вернется на К .

 

Действительно, в силу уравнения Гамильтона

 

к =

^

 

ё

8G

где

определяется формулой

Ж0 = |(asin2I + bcos2l)(G2 - L2)^L 2,

будем иметь

g = G(a sin21+ bcos21).

Так как

G = G° ф 0 (безразличное равновесие тела,

ког­

да G° =

0, мы исключаем из рассмотрения), то g > е

> 0.

Следовательно, за конечное время переменная g станет рав­ ной 27г , и точка к снова вернется на кольцо К , подойдя к ней с «противоположной стороны».

60

Глава 3

Таким образом, получаем естественное отображение S внутренности кольца К на себя. Используя интегральный ин­ вариант Пуанкаре-Картана, нетрудно показать, что S сохра­ няет меру

v «»-// dL dl

D

(см., например, [23, 24]). Преобразование S по непрерывности вращает границы К в противоположных направлениях.

Задача Эйлера-Пуансо интегрируема, поэтому фазовое пространство расслоено на замкнутые двумерные поверхности

{G = G0, Ж0 = ft},

(1.7)

которые являются в общем случае двумерными торами. Оче­ видно, что поверхности (1.7) пересекаются с кольцом К по замкнутым инвариантным кривым отображения S, которые совпадают с линиями уровня функции fflo(lgLG°). Значит, кольцо К расслоено на замкну­ тые инвариантные кривые отобра­ жения S. Эти кривые показаны на рис. 9 (чтобы получить эту картин­ ку, надо отождествить на рис. 5 точ­ ки, Z-координаты которых отлича­

ются на 27г).

Неподвижные точки 1 и 2 соот­ ветствуют периодическим решени­ ям — вращениям вокруг меньшей оси инерции в противоположных на­ правлениях. Точки 3 и 4 тоже явля­

ются неподвижными точками отображения S. Им отвечают постоянные вращения вокруг средней оси инерции, которые имеют гиперболический тип. В невозмущенной интегрируе­ мой задаче сепаратрисы неподвижных точек 3 и 4 сдвоены: они идут из седла в седло. Эти сепаратрисы разбивают коль­ цо К на четыре связные подобласти Ki (i = 1, 2, 3, 4).

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