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

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

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

§ 4. Несуществование дополнительного интеграла

71

Следствие 4. Если А > В > С и хотя бы один из интег­ ралов (4.4) отличен от нуля, то уравнения Эйлера-Пуассона задачи о вращении твердого тела вокруг неподвижной точки в ньютоновском поле сил не имеют дополнительного анали­ тического интеграла, независимого от классических.

Доказательство теоремы 2.

 

Предположим,

 

что

существует

новый интеграл

&(p,q,r, 71, 72, 73)5

аналитический в

Е

С

R 6.

Введем

в уравнения Эйлера-Пуассона ма­

лый параметр р,

заменяя л

на рл-

Тогда &(р, q, г, р л , р л , Р1з)

можно

представить в виде

степенного ряда

по р.

Этот ряд будет сходиться для достаточно малых значений р, ес­ ли (р, q, г) £ G (G — некоторая малая окрестность нуля в R 3{p, q, г}), а л принимают значения из сферы Пуас­ сона

S2 = {TI + т ! + 7з : !}•

После введения малого параметра ц уравнения (4.1) будут иметь тот же вид, а потенциал станет равным

У = рУ1 + р2У2 + ■- •-

Таким образом, новые уравнения (4.1) будут описывать дви­ жения твердого тела в слабом потенциальном поле. Заметим, что форма рУ1имеет смысл потенциальной функции силы тя­ жести, а ее коэффициенты можно интерпретировать как про­ изведения веса тела на координаты центра тяжести в главных осях инерции.

Переменные 71, 72, 7з выражаются через специальные ка­ нонические переменные L, G, Н, /, g, h по формулам (4.1)

гл. II, а переменные р, q, 1-----по формулам:

 

р = ± V G 2 - £ 2sinI,

q = ± -V G 2 - L 2cosl, г =

(4.5)

A

i>

О

72 Глава 3

Обозначим через D область на кольце К , содержащую при каждом значении G G (ад, аг), 0 < ai < аг ((ад, аг) — малый интервал R) обе сепаратрисы неподвижных точек 3 и 4 (см. § 3).

Согласно формулам (4.1) (гл. II) и (4.5) переменные Эй­ лера-Пуассона являются аналитическими функциями специ­ альных канонических переменных, если

G G (ai, a 2), \Н\ < а±, (L ,l)E D , g G Т 1{ж mod 27г}.

Подставляя выражения (4.1) (гл. II) и (4.5) в функцию &, получим при фиксированном значении переменной Н = Н° (|if°| < ai) интеграл канонических уравнений с функцией Гамильтона (3.5), аналитический в области

D х (ад, аг) х mod 27г} х (—е, е).

Так как Ж и 9 как функции от переменных L, I, G, g и параметров Н, ц независимы, то при некотором фиксирован­ ном значении Н = Н° из интервала (—а±, ai) новый интеграл уравнений с гамильтонианом (3.5) не зависит от интеграла энергии. Но это противоречит заключению теоремы 1 (с уче­ том замечания 3 предыдущего параграфа). ■

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

Задача о существовании дополнительного интеграла урав­ нений вращения тяжелого твердого тела вокруг неподвиж­ ной точки, аналитического по каноническим переменным и параметру //, впервые поставлена А. Пуанкаре в п. 86 его «Новых методов небесной механики». Анализируя разложение возмущающей функции, А. Пуанкаре показал, что (в нашей терминологии) вековое множество не является всюду плот­ ным, и, следовательно, его общая теорема об отсутствии новых аналитических интегралов не применима: <г. . . ничто не препятствует существованию третьего однозначного ин­ теграла, если только якобиан трех интегралов обращается в нуль, как только п' (у нас и>2, В. К.) становится кратным п (у нас LJi, В. К.); отсюда следует, что этот третий интеграл не может в общем случае быть алгебраическим.

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

73

Поскольку сформулированные в этой главе условия явля­ ются необходимыми, но не достаточными, ничто не доказы­ вает, что этот третий интеграл существует. . . » [1, с. 226].

Фактически в этих строках содержится постановка задачи о существовании новых алгебраических интегралов

иутверждение, которое сейчас принято называть теоре­ мой Пуанкаре-Гюссона: если твердое тело несимметрично,

ицентр тяжести не совпадает с точкой подвеса, то ново­ го алгебраического интеграла нет. Общая задача А. Пуанкаре

одополнительных алгебраических интегралах решена в рабо­ тах Лиувилля, Гюссона, Бургатти. Исследования этих ав­ торов аналогичны работам Брунса и Пенлеве, посвящен­ ным задаче трех тел. Прекрасное изложение истории во­ проса об алгебраических интегралах можно найти в статье П. Я. Полубариновой-Кочиной [44]. Следует сослаться также на недавние работы А. И. Докшевича [88, 89], в которых ис­ правлены и упрощены некоторые существенные моменты до­ казательства теоремы об отсутствии новых алгебраических интегралов.

Г л а в а IV

Динамические эффекты, препятствующие интегрируемости уравнений движения несимметричного тела

Теория рождения периодических решений в каноничес­ ких системах дифференциальных уравнений, близких к ин­ тегрируемым, была разработана А. Пуанкаре для целей не­ бесной механики. В данной главе устанавливается примени­ мость этих результатов к классической задаче о движении тяжелого твердого тела с закрепленной точкой. Тем самым удается существенно расширить класс известных периодичес­ ких решений. В этой же главе исследовано воздействие воз­ мущения на сепаратрисы неустойчивых периодических реше­ ний задачи Эйлера-Пуансо — постоянных вращений вокруг средней оси эллипсоида инерции. Сложное поведение траек­ торий уравнений движения несимметричного тяжелого твер­ дого тела вокруг неподвижной точки (в частности, рожде­ ние многочисленных изолированных периодических решений и расщепление сепаратрис) несовместимо с существованием нового независимого аналитического интеграла.

§ 1. Характеристические показатели. Теорема Пуанкаре о периодических

решениях

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

§ 1. Характеристические показатели

75

Рассмотрим аналитическую систему дифференциальных

уравнений, определенную в R n:

 

 

Х\

/1

 

X = /(ж ), X =

, / =

 

Хп

fn

 

Эта система имеет единственное аналитическое решение x(t, £) такое, что

6

Ф , £ ) = £, С

£»7

(см., например, [3]). Предположим, что при £ = £* решение системы (1.1) периодическое с периодом Т*. Частные произ­ водные

дхк

дЬ

удовлетворяют линейной системе

 

dfk

%kl

дх, •Erl)

 

r = l

где в выражении для d fk/дхг вместо х подставлена периоди­ ческая функция x(t, £*). Введем матрицу n-го порядка

F = dfk дхг

Линейные дифференциальные уравнения с периодическими коэффициентами

2/1

 

y = F y , У =

(1.2)

Уп

 

называются уравнениями в вариациях периодического реше­ ния ж(£, £*).

76

Глава 4

Матрица

т =II*WWII

является, очевидно, фундаментальной матрицей линейной системы (1.2), причем Х(0) = Е, где Е — единичная мат­ рица п-го порядка. Матрица Х (Т*) называется матрицей монодромии системы линейных уравнений (1.2). Собственные числа Ai, ... , Ап матрицы Х (Т *) называются мультиплика­ торами, а величины ад, ... , а п: А* = ехр(а*Т*) — характе­ ристическими показателями системы (1.2). Числа ад, ... , ап называются также характеристическими показателями пери­ одического решения x(t, £*). Система (1.1) предполагается ав­ тономной, поэтому хотя бы один из характеристических по­ казателей равен нулю. Действительно, в этом случае

\\Х(Т*)-Е\\Г =

0, /*

=

/ ( П ^ 0 .

(1.3)

Мы будем, в основном, иметь дело с периодическими ре­

шениями автономных канонических систем

 

dx =

дЖ^

 

dy _

 

{)■/(!

 

dt

ду

dt

 

Ox

 

 

 

 

XI

 

2/1

 

Ж = Ж (х, у),

х =

>

У =

 

 

 

 

Хп

 

Уп

 

Для таких систем справедлива

Теорема Пуанкаре —Ляпунова. Характеристичес­ кий многочлен р(А) матрицы монодромии периодического решения гамильтоновой системы (1.4) возвратный: (р)А = = А2> (1 /А ).

Доказательство можно найти, например, в [4, 20]. Значит, если А — корень характеристического уравне­

ния р(х) = 0, то 1/А — корень того же уравнения. Так как многочлен р(х) имеет действительные коэффициенты, то его корнями являются также числа А и 1/А.

§ 1. Характеристические показатели

77

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

В случае п = 2 оставшиеся характеристические показате­ ли либо оба действительны, либо оба чисто мнимы. Если они отличны от нуля, то периодическое решение называется невы­ рожденным. Невырожденное решение с действительными ха­ рактеристическими показателями называется гиперболичес­ ким, а с чисто мнимыми — эллиптическим. Эллиптическое решение устойчиво в первом приближении, а гиперболичес­ кое неустойчиво.

Предположим теперь, что система (1.1) имеет т (0 ^ т ^ ^ п 1) первых интегралов ... , &т(х), аналитических в R ". Докажем, что если в точках траектории периодического решения x(t, £*) ранг матрицы

д&1 dxi

д&т

1—1

дхп

д&т

дх„

равен то, то по меньшей мере т + 1 характеристических по­ казателей этого решения равны нулю.

Действительно, так как функция ^ (ж ) — первые интег­ ралы, то

@ k {x{tA *))= & k { О ,

к = 1, . . . , то .

 

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

 

 

 

 

 

д&кХ Ц + . . . +

- j r - ^ X n i

=

(/ “ 1 ; ' ' '

’ п)'

(1'5)

дхг ................

дохп ~'ш

_ дх{

Положим t = Т*. Тогда ж(Т*, £*) =

и равенство (1.5) можно

представить в следующей матричной форме:

 

 

\\Х(Т*) - Е\\щ =

0, щ

= g rad ^ fc, к = 1,

... , то.

(1.6)

С другой стороны,

 

 

 

 

 

( щ , /) = 0.

78

Глава 4

Так как /(£*) = /* ф 0, то из этих равенств вытекает линей­ ная независимость векторов щ, ... , i)n, f* . Из формул (1.3) и (1.6) следует, что щ, ... , т]п, /* — собственные векторы матрицы

Y = ||Х(Т*) - Е\\

с собственными значениями А = 1. Следовательно, по крайней мере ш + 1 ее характеристических чисел (мультипликаторов) равны единице. Это и требовалось доказать.

В качестве следствия получаем следующее утверждение: если п — m характеристических показателей решения x(t, £*) отличны от нуля, то на его траектории интегралы S i, ... , S m зависимы. В частности, на траекториях невырожденных пе­ риодических решений двумерной гамильтоновой системы ин­ теграл энергии Ж и любой первый интеграл 3- зависимы.

Рассмотрим теперь систему дифференциальных уравне­

ний

х = /(ж, а),

(1.7)

правая часть которой аналитична в прямом

произведении

R " х (—е, е; х G R ", a G (—£, е). При всех значениях парамет­ ра a G (—е, г) существует аналитическое решение x(t, £, а) системы (1.7) такое, что ж(0, £, а) для всех а. Предположим, что решение x(t, £*, 0) — периодическое с периодом Т*. Бу­ дет ли «возмущенная» система (1.7) допускать периодические решения, если а ф 0, но мало?

Для рассматриваемых ниже приложений особый интерес представляет случай, когда система (1.7) имеет первый ин­ теграл g(x, а), аналитический в R™ х (—е, е).

Без ущерба общности можно считать, что /„(£*, 0) ф 0. Так как

dg

dxi fi +

то одна из производных

 

 

Sr =

dg

(Г = 1,

П —1)

дхг «=е,«=о

 

 

отлична от нуля. Пусть, например, gi ф 0.

§ 1. Характеристические показатели

79

Введем в рассмотрение матрицу V, которая получается из матрицы Y «порождающего» решения x(t, £*, 0) вычерки­ ванием последнего столбца и первой строки.

Теорема (А . Пуанкаре). Если F / 0, то при доста­ точно малых а существует аналитическая функция £(а), т = £*, такая, что решение x(t, £(а), а) системы (1.7) — периодическое с периодом Т * «порождающего» решения x(t, £*, 0).

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

Условие периодичности решения x(t, £, а) имеет следую­ щий вид:

а) = хк(Т*, £ , а ) - £ к = 0, к = 1, ... ,п .

Полагая £п = £*, найдем п -

1 функций ^I (QI),

... , £n_i(a),

удовлетворяющих уравнениям

 

 

a ) = 0 ,

к = 2, ... , п.

(1.8)

Эти уравнения, очевидно, удовлетворяются, если положить а = 0, = ££ = 2, ... , п). Соответствующая функцио­ нальная матрица

&Z2

 

&Z2

dti

'

e t„ -i

dzn

'

dzn

Oil

d£n-i

совпадает с матрицей V. Так как |V| ф 0, то по теореме о не­ явных функциях существуют аналитические решения £fe(a), £&(0) = 1, ... , п —1) системы (1.8). Осталось показать, что при этом выполнено уравнение zi(£, a) = 0. Действительно, по теореме о среднем значении

g(x(T*, £, a), a) -g (£ , а) =

dg

а) = 0,

(1.9)

 

dxi

х = х

 

где х к = £к (к ф 1) и xi лежит между и х\(Т *, £, а). Так как gx = dg/dxi ф 0 при £ = £* и а = 0, то при малых а будем иметь dg/dxx ф 0. Но тогда из (1.9) вытекает, что zi(£, a) = 0.

80

Глава 4

Для некоторых целей выгодно искать периодические ре­ шения «возмущенной» задачи, расположенные на фиксирован­ ном уровне g(x, а) = g* = g(£*, 0). Тогда к уравнениям пери­ одичности

Zk(T, i, а) = х к(Т, £, а) - £*. = 0 (к = 1, ... , га)

(1.10)

надо добавить уравнение

2n+i(£> а) = g(£, а) - g* = 0.

(1.11)

При этом в уравнениях (1.10) период Т является неизвестной функцией а.

Полагая £„ = £*, ищем решение £ i(a ,), ... , ^n_ i(a ), Т(а) уравнений (1.10), когда к ф 1, и (1.11). Этим уравнениям при а = 0 удовлетворяют величины , £*_l5 Т*. Соот­ ветствующая функциональная матрица

 

 

дгг

dZ2

д ь

din-i

dT

dzn-\-i

dZn+1

dz„+i

dii

din-1

dT

совпадает при а = 0 с матрицей W , которая получается из

матрицы

 

 

 

 

 

Y

Г

 

( 1.12)

 

г]

0

 

 

 

 

вычеркиванием n-го столбца

и первой

строки. В матри­

це (1.12) /* = /(£*, 0), а г) =

(dg/dxx, ...

, dg/dxn), где вмес­

то х подставлено i* . Если |W| ф 0, то по теореме о неявных функциях при малых а существуют периодические решения системы (1.7), расположенные на интегральной гиперповерх­ ности g(x, а) = g*.

§ 2. Возмущение равномерных движений

Когда параметр Пуанкаре /х равен нулю, имеем случай Эйлера-Пуансо. В этой невозмущенной интегрируемой задаче