Главная > ИНТЕГРИРУЕМЫЕ ГАМИЛЬТОНОВЫ СИСТЕМЫ (А. В. Болсинов А. Т. Фоменко)
НАПИШУ ВСЁ ЧТО ЗАДАЛИ
СЕКРЕТНЫЙ БОТ В ТЕЛЕГЕ
<< Предыдущий параграф Следующий параграф >>
Пред.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
След.
Макеты страниц

Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше

Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике

ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO

5.2.1. Топология 3-поверхности и бифуркационная диаграмма
В этом параграфе мы опишем топологический тип изоэнергетических 3 -поверхностей для интегрируемых случаев, перечисленных в параграфе 1 .

Изоэнергетическая 3 -поверхность Q3 — это совместная поверхность уровня функций f1,f2,H, заданных на евклидовом пространстве R6(S,R), где
f1=R12+R22+R32,f2=S1R1+S2R2+S3R3,

а H — это гамильтониан. Поскольку мы считаем, что f1=1, то различные 3 поверхности Q3 задаются двумя параметрами g и h, т. е. значениями функций f2 и H. Описание топологического типа Q3 будет дано нами в следующем виде. Мы рассмотрим бифуркационную диаграмму для пары интегралов f2 и H. В результате на плоскости R2(g,h) появятся кривые, разбивающие плоскость на области таким образом, что для всех точек (g,h) из одной области топологический тип соответствующих изоэнергетических поверхностей
Q3={f1=1,f2=g,H=h}

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

Исследование механических систем при помощи бифуркационных диаграмм отображения f2×H было начато С. Смейлом [178]. Итак, рассмотрим отображение
F=f2×H:S2×R3R2(g,h),

заданное обычной формулой F(P)=(f2(P),H(P))R2(g,h), где точка P принадлежит S2×R3. Отметим, что для функций f2|S2×R3 и H|S2×R3 мы сохраняем прежние обозначения f2 и H. Образом множества критических точек отображения F является бифуркационная диаграмма Σ в плоскости R2(g,h). Полный прообраз любой точки (g,h)otinΣ является неособой изоэнергетической 3 -поверхностью Q3={f1=1,f2=g,H=h}, что следует из теоремы о неявных функциях. Для гамильтонианов вида (1.9) (а только такие мы и будем рассматривать в настоящей главе) отображение F будет собственным, то есть прообразом любого компакта является компакт. Отсюда следует, что для всех точек (g,h), принадлежащих одной связной компоненте множества R2(g,h)Σ, топологический тип 3 -многообразия Q3 будет одним и тем же.

Как определить тип изоэнергетической 3 -поверхности Q для каждой компоненты связности R2(g,h)Σ ? На этот вопрос помогает дать ответ следующее утверждение. См. работу С. Смейла [178].
Предложение 5.2. Пусть H=CS,S+W,S+V, где C — постоянная симметричная положительно определенная матрица размера 3×3,W(R) — произвольная вектор-функция, и V(R) — произвольная гладкая функция. Рассмотрим проекцию π(Q3) изоэнергетической 3-поверхности Q3 на двумерную сферу Пуассона, задаваемую уравнением R12+R22+R32=1. Здесь π(R,S)=R. Пусть π(Q3) гомеоморфна сфере с т дырками. Тогда:
1) если m=0, то Q3 диффеоморфна RP3,
2) если m=1, то Q3 диффеоморфна S3,
3) если m>1, то Q3 диффеоморфна связной сумме \# S1×S2) в количестве m1 экземпляров.

Доказательство.
Рассмотрим проекцию многообразия Q3 на 2-сферу Пуассона S2. Возьмем произвольную точку R=(R1,R2,R3) из образа π(Q3) и посмотрим, как устроен слой над нею. Этот слой является пересечением 2 -плоскости R1S1+R2S2+R3S3=g с 2 -эллипсоидом H=CS,S+W,S+V= const. Следовательно, он гомеоморфен либо окружности, либо точке, либо пуст. Если точка R лежит внутри образа π(Q3), то легко видеть, что слой гомеоморфен окружности. Если точка R лежит на границе образа π(Q3), то слой является точкой. Если же точка R расположена вне образа π(Q3), то слой пуст. Следовательно, топологически 3 -многообразие Q3 может быть описано следующим образом. Если проекция π(Q3) совпадает со всей сферой Пуассона, то Q3 представляет собой S1-расслоение над сферой, топологически эквивалентное, как нетрудно увидеть, расслоению единичных касательных векторов. Следовательно, Q3=RP3.

Если проекция π(Q3) не совпадает со всей сферой, т.е. имеет дырки, то Q3 можно можно получить следующим образом. Сначала рассмотрим прямое произведение π(Q3)×S1, а затем над каждой граничной точкой проекции π(Q3) сожмем слой S1 в точку. Трехмерные многообразия, полученные таким способом, мы уже встречали ранее в главе 4 первого тома, и их топология была описана в предложении 4.5. Последние два утверждения фактически являются переформулировкой этого предложения. Предложение 5.2 доказано.

Таким образом, построив бифуркационную диаграмму отображения F и указав топологический тип Q3 для каждой из областей, на которые Σ разбивает плоскость R2(g,h), мы получаем полное описание всех топологических типов неособых изоэнергетических 3 -поверхностей исследуемой гамильтоновой системы.

В работах С.Б.Каток [78], Я. В. Татаринова [191], [192] по изложенной выше схеме исследована топология 3 -поверхностей Q3 для задач о движении твердого тела с закрепленной точкой. В этих работах построены бифуркационные диаграммы отображения f2×H для гамильтонианов достаточно общего вида, вообще говоря, неинтегрируемых. В частности, рассмотрены интегрируемые случаи Эйлера, Лагранжа и Ковалевской. Приведем здесь подробное описание этих бифуркационных диаграмм для интегрируемых случаев.

Чтобы сделать это, нам потребуется некоторый общий прием, который мы сейчас опишем. Напомним некоторые обозначения. Систему уравнений, описывающих движение тяжелого твердого тела (или, в более общем случае, гиростата) с неподвижной точкой, можно записать в виде уравнений Эйлера на коалгебре Ли e(3), где введены линейные координаты (S1,S2,S3,R1,R2,R3). Относительно них скобка Пуассона имеет вид
{Si,Sj}=εijkSk,{Si,Rj}=εijkRk,{Ri,Rj}=0.

Ядро этой скобки порождается функциями f1=R12+R22+R32 (геометрический интеграл) и f2=S1R1+S2R2+S3R3 (интеграл площадей). Все регулярные совместные поверхности уровня этих функций вида {f1=1,f2=g}

диффеоморфны касательному расслоению к двумерной сфере. Гамильтониан этой задачи имеет следующий вид:
H=(S1+λ1)22A1+(S2+λ2)22A2+(S3+λ3)22A3+U(R1,R2,R3),

где A1,A2,A3 — главные моменты инерции тела, (λ1,λ2,λ3) — гиростатический момент, U(R1,R2,R3) — потенциал, который для рассматриваемой задачи имеет вид a1R1+a2R2+a3R3. Итак, требуется описать топологию изоэнергетических поверхностей Qh,g3={f1=1,f2=g,H=h} для различных констант A1,A2,A3,λ1,λ2,λ3,g,h и функций U(R1,R2,R3).

Введем следующие обозначения S=(S1,S2,S3),R=(R1,R2,R3),λ= =(λ1,λ2,λ3),a=(a1,a2,a3),A=diag(A1,A2,A3), где S,R,λ,a трехмерные векторы, а A — диагональная матрица. Будем считать, что все эти векторы из одного и того же евклидова пространства R3, в котором задано евклидово скалярное произведение ,. Тогда f1=R,R,f2=S,R, H=12A1(S+λ),S+λ+U(R).

Рассмотрим проекцию π:(S,R)R. Поскольку f1(R)=1, то образ 3 -многообразия Qh,g3 при этой проекции есть некоторое подмножество сферы Пуассона S2={R12+R22+R32=1}R3(R1,R2,R3). Очевидно, что точка R сферы Пуассона S2 принадлежит образу π(Qh,g3) тогда и только тогда, когда существует решение системы
S,R=g,A1(S+λ),S+λ=2(hU(R)),

где неизвестными являются (S1,S2,S3)=S. Первое уравнение задает плоскость, а второе — эллипсоид в пространстве R3(S1,S2,S3), если hU(R)>0. Легко проверяется, что они пересекаются тогда и только тогда, когда
(g+λ,R)22AR,R+U(R)h.

При этом пересечение является окружностью в случае строгого неравенства или точкой в случае равенства.
Назовем приведенным потенциалом следующую функцию:
φg(R)=(g+λ,R)22AR,R+U(R),

заданную на сфере Пуассона.
Оказывается, топология изоэнергетических поверхностей Qh,g3 и их перестройки при изменении значения энергии h полностью определяются функцией φg(R). А именно, комбинируя предложение 5.2 с приведенными выше соображениями, получаем следующую теорему.
Теорема 5.1. Если значение h приведенного потенииала φg(R) не является критическим значением, то изоэнергетическая 3 -поверхность Qh,g3 является гладким трехмерным многообразием. B этом случае множество {φg(R)h} на сфере Пуассона представляет собой объединение непересекающихся двумерных подмногообразий с краем Pi1,,Pim, где Pk2-диск с k дырками, а 3-многообразие Qh,g3 гомеоморфно несвязному объединению трехмерных многообразий Ni1,,Nim, где N0 — трехмерная сфера, а Nk, при k1, — это связная сумма k экземпляров S1×S2. Если h меньше минимального значения приведенного потенциала φg(R), то Qg,h3 пусто. Если же h больше максимального значения функции φg(R), то Qg,h3=RP3, т. е. проективное пространство.

Используя эту теорему, можно предложить следующий алгоритм описания топологического типа изоэнергетических 3 -поверхностей Qg,h3.

Шаг 1. Строим граф Риба функции φg(R). Если значение h меньше минимума или больше максимума функции φg(R), то 3 -многообразие Qh,g3 либо пусто, либо гомеоморфно RP3 соответственно.

Шаг 2. Пусть h принадлежит множеству значений функции φg(R). Разрежем граф Риба по уровню h и рассмотрим его нижнюю часть. Количество связных компонент в нижней части совпадает с количеством связных компонент 3 -поверхности Qg,h3. При этом, если связная компонента графа имеет k граничных точек, т.е. следов разреза, то соответствующая связная компонента 3 -поверхности Qg,h3 гомеоморфна связной сумме k1 экземпляров S1×S2 при k1, или же является 3 -сферой S3 при k=1.
Опишем теперь метод построения бифуркационных диаграмм.
Пусть параметры гамильтониана фиксированы. Тогда топологический тип изоэнергетической 3 -поверхности Qg,h3 определяется константами g и h. Эту зависимость удобно изображать при помощи бифуркационной диаграммы, т. е. кривых на плоскости с координатами (g,h), точкам которых соответствуют особые изоэнергетические поверхности. Если такие кривые построены, то они разбивают плоскость (g,h ) таким образом, что точкам из одной области соответствуют гомеоморфные изоэнергетические 3 -поверхности.

Согласно теореме 5.1, для каждой точки ( g0,h0 ) такой кривой значение h0 является критическим значением приведенного потенциала φg(R). Разберем теперь отдельные случаи интегрируемости.
5.2.2. Случай Эйлера

Здесь λ=0, и U(R) тождественно равно нулю. Приведенный потенциал имеет здесь вид:
φg(R)=g22AR,R.

При geq0 критические точки этой функции на сфере Пуассона такие же, как и у квадрики AR,R, то есть те точки, где соответствующий эллипсоид касается сферы R,R=1. При каждом geq0 имеется шесть критических точек:
R=(±1,0,0),(0,±1,0),(0,0,±1).

Топология лиувиллевых слоений в классических интегрируемых случаях 205
Им соответствуют три критических значения: hi=g22Ai. При g=0 приведенный потенциал имеет единственное критическое значение, равное нулю. Таким образом, бифуркационная диаграмма в случае Эйлера состоит из трех парабол (рис. 5.1).

Топологический тип изоэнергетических 3 -поверхностей устанавливается при помощи графов Риба, которые изображены на рис. 5.2. В случае, когда

Рис. 5.1 среди главных моментов инерции есть одинаковые, некоторые вершины графа Риба стягиваются в одну. На диаграмме соответствующие параболы также сливаются в одну. Итак, подведем итоги.
Рис. 5.2
Бифуркационная диаграмма отображения f2×H для гамильтониана случая Эйлера
H=S122A1+S222A2+S322A3

имеет простой вид, показанный на рис.5.3. Она состоит из трех парабол h=g22Ai,i=1,2,3, разбивающих плоскость R2(g,h) на 6 областей. В каждой из этих областей на рис. 5.3 указан топологический тип 3 -многообразия Q3. Знак означает, что для всех точек (g,h) из данной области прообраз при отображеРис. 5.3 нии f2×H пуст, т.е. (f2×H)1(g,h)=. На рис. 5.3 изображен случай, когда параметры A1,A2,A3 твердого тела связаны соотношением 0<A3<A2<A1. Если же 0<A3=A2<A1, то две верхние параболы сливаются в одну. Если же 0<A3<A2=A1, то сливаются две нижние параболы. При A1=A2=A3 гамильтониан H случая Эйлера становится резонансным, поскольку для него появляются два функционально независимых интеграла, например, функции S1 и S2. Такие, т. е. резонансные, случаи мы пока рассматривать не будем.
5.2.3. Случай Лагранжа
Для классического гамильтониана случая Лагранжа
H=12(S12+S22+S32β)+R3,

где β>0, бифуркационные диаграммы качественно отличаются в следующих случаях:
a) 0<β<1,
b) β=1,
c) 1<β43,
d) β>43.
Все они изображены на рис. 5.4. Все бифуркационные диаграммы состоят из двух парабол h=g22β±1 и кривой, при βeq1, которую можно задать параметрически следующим образом:
g=t+1(β1)t3,h=t22+32(β1)t2,t21|β1|.

Эта кривая касается одной из парабол в точках
(±β|1β|,3β22|1β|)R2(g,h)

и при β>43 имеет две точки возврата с координатами
(±433β14,3β1).

Топологический тип Q3 указан на рис. 5.4. Здесь S3(S1×S2) обозначает несвязное объединение S3 и S1×S2.
Теперь мы подробно опишем, как получаются эти результаты.
Рассмотрим сначала гамильтониан более общего вида, чем для классического случая Лагранжа. Он задается следующими параметрами:
A1=A2=B,λ1=λ2=0,U(R)=U(R3).

Рис. 5.4
Приведенный потенциал имеет вид
φg(R)=(g+λ3R3)22(BR12+BR22+A3R32)+U(R3)==(g+λ3R3)22(B+(A3B)R32)+U(R3).

Он зависит только от координаты R3. Поэтому для любых значений g и h область, выделяемая на сфере Пуассона условием φg(R)h, состоит из некоторого числа колец и, возможно, одного или двух дисков. Следовательно, изоэнергетические 3 -поверхности Qg,h3 гомеоморфны несвязному объединению некоторого количества многообразий S1×S2 и, возможно, одной или двух трехмерных сфер. Рассмотрим теперь подробнее обычный случай Лагранжа, т.е. с линейным относительно координаты R3 потенциалом. Без ограничения общности можно считать, что параметры гамильтониана имеют следующий вид:
A1=A2=1,A3=β,λ=(0,0,0),U(R)=R3.

Тогда приведенный потенциал запишется так:
φg(R)=f(R3)=g22(1+(β1)R32)+R3.

Рис. 5.5
Пусть 0<β<1, т.е. эллипсоид инерции вытянут. Тогда, как легко проверить, f(1)=g22β1,f(1)=g22β+1 и f(R3)>0 на всем отрезке [1,1]. Поэтому с качественной точки зрения для вида графика возможны два варианта, показанные на рис. 5.5a и рис. 5.5b. Выбор того или иного варианта определяется знаком производной f(1). Вычисляя производную, получаем
f(1)=1g2(1β)β2.

Таким образом, случаи (a) и (b) получаются соответственно при g2<β21β и при g2>β21β. Графы Риба изображены на рис. 5.6.
Рис. 5.6
Рис. 5.7
В критических точках, соответствующих вершинам степени 1 графа Риба, функция f принимает значения ±1+g22β. Минимальное значение в случае (b) можно вычислить следующим образом. Равенство нулю производной
f(R3)=1g2(β1)R3(1+(β1)R32)2

возможно лишь при (β1)R3>0. Обозначим (β1)R3 через t2. Тогда |t|2(1β)1, поскольку |R3|1. Переписывая условие f(R3)=0 через t, а также подставляя (β1)1t2 вместо R3, получаем
g=t+1t3(β1),h=t22+32(β1)t2,|t|1|1β|.

Рассматривая эти уравнения как параметрическое задание кривой на плоскости (g,h), получаем, что бифуркационная диаграмма в рассматриваемом случае состоит из данной кривой и двух парабол h=±1+g22β. Она изображена на рис. 5.7.
Рис. 5.8
В случае β=1 приведенный потенциал линеен: f(R3)=R3+g22. Бифуркационная диаграмма состоит только из двух парабол h=±1+g2.

Пусть теперь β>1, то есть эллипсоид инерции сплюснут. Возможны четыре типа графиков функции f(R3). См. рис.5.8. Соответствующие графы Риба приведены на рис. 5.9.
Рис. 5.9
Какой из вариантов реализуется для различных значений β и g, определяется количеством нулей производной f(R3) на интервале (1,1) или, что то же самое, количеством нулей полинома P(x)=(1+(β1)x2)2g2(β1)x на интервале (1,1) и значениями функции f в этих нулях.

Исследуя все возможности, получаем ответ. Бифуркационная диаграмма состоит из двух парабол и кривой, задаваемых теми же уравнениями, что и в случае β<1. Она имеет различный вид при 1<β43 и при β>43. См. соответственно рис. 5.10a и рис. 5.10 b.

Рис. 5.10
Точки касания кривой с параболой, как и в случае β<1, имеют на плоскости (g,h) координаты
(±β|1β|,3β22|1β|).

В случае β>43 имеются две точки возврата, координаты которых таковы:
(±433β14,3β1).
5.2.4. Случай Ковалевской
Сначала рассмотрим специальный гамильтониан
H=12(S12+S22+βS32)+R1.

Рис. 5.11
Вообще говоря, он не является гамильтонианом Ковалевской. Для гамильтониана (2.1) бифуркационные диаграммы имеют в точности такой же вид, что и для случая Лагранжа. Они получаются из диаграмм, изображенных на рис. 5.4, сжатием по оси g в β раз. Однако топологический тип изоэнергетических 3 -поверхностей Q3 будет другим. Специальный гамильтониан (2.1) соответствует случаю Ковалевской при β=2. Поэтому на рис.5.11 приведена бифуркационная диаграмма лишь для β>43. При β=2 координаты точек касания кривой и параболы таковы: (±2,2). Координаты точек, в которых кривая трансверсально пересекает параболу, таковы: (±221,2(21)). Топологический тип Q3 указан на рис. 5.11. Трехмерную поверхность, обозначенную здесь через K3, можно описать следующим образом. Рассмотрим подрасслоение расслоения единичных касательных векторов к двумерной сфере, базой которого B3 является 2-сфера с тремя дырками.

Другими словами, мы выбрасываем из 2-сферы три открытых непересекающихся диска. Затем каждый слой, т.е. окружность, над границей базы B стянем в точку. Полученное трехмерное многообразие и обозначено через K3. Как мы уже знаем, в действительности K3 является связной суммой двух экземпляров S1×S2.

Напомним, что таким же способом, рассматривая базу Bm, т. е. 2-сферу S2 с m дырками, можно получить проективное пространство RP3 при m=0, трехмерную сферу S3 при m=1, прямое произведение S1×S2 при m=2. Как мы показали выше, при m2 получающиеся указанным способом 3 -многообразия гомеоморфны связной сумме m1 экземпляров S1×S2. В частности,
K3=(S1×S2)#(S1×S2).

На самом деле, такое описание изоэнергетических 3 -поверхностей для гамильтонианов вида (1.9) получается естественным образом. Рассмотрим проекцию на сферу Пуассона π:TS2S2, где
TS2={f1=1,f2=g}R6(S,R),S2={f1=1}R3(R).

При проекции π изоэнергетическая 3 -поверхность
Qh3={f1=1,f2=g,H=h}

проектируется в некоторое замкнутое подмножество P на сфере Пуассона. Причем, если z — это внутренняя точка P, то π1(z)Qh3 является окружностью. Если z лежит на границе P, то π1(z)Qh3 — точка. Выше было показано, что P={RS2R3(R)φg(R)h}, где
φg(R)=(g+λ1R1+λ2R2+λ3R3)22(A1R12+A2R22+A3R32)+U(R).

Используя явный вид функции φg(R), можно установить вид множества P в S2, то есть найти число дырок. Это позволяет определить топологический тип Q3. Как было отмечено выше, это либо RP3, либо S3, либо связная сумма #i=1mS1×S2.

Этот метод был разработан Смейлом в [178] для натуральных механических систем. Случай систем с гироскопическими силами рассмотрен М. П. Харламовым в [219].

В заключение опишем подробности вычисления бифуркационной диаграммы в случае Ковалевской. Для параметров гамильтониана случая Ковалевской
A1=A2=1,A3=12,λ=(0,0,0),U(R)=R1

приведенный потенциал имеет вид
φg(R)=g22R32+R1.

Рис. 5.12
Рассмотрим эту функцию в координатах (R1,R3 ) на полусфере {R2>0}. Отметим, что на полусфере {R2<0} картина симметрична.

Нарисуем линии уровня функции φg(R) в проекции на плоскость (R1,R2). Уравнение φB(R)=c равносильно следующему:
R1=cg22R32.

Возникают три качественно различные картины линий уровня на диске

{R12+R321}. См. рис. 5.12. Здесь R3 — горизонтальная ось, а R1 — вертикальная. При этом в случае (c) возможны два варианта в зависимости от значений функции φg(R) в точках максимума. Соответствующие графы Риба изображены на рис. 5.13.

Значения функции φg(R) в точках (±1,0,0) равны g22±1. Это задает две параболы h=g22±1, принадлежащих бифуркационной диаграмме. Вычисляя значения для остальных точек, получаем еще одну кривую (рис.5.14). Вычисления аналогичны тем, которые проводились выше для случая Лагранжа. Более того, полученная кривая отличается от кривой для случая Лагранжа, при β=2, лишь сжатием по оси g в 2 раз. Однако, как видно из построения графов Риба, топологический тип изоэнергетических 3 -поверхностей отличается от случая Лагранжа. Окончательный ответ приведен на рис.5.11.

Как отмечалось выше, описанные бифуркационные диаграммы для случаев Эйлера, Лагранжа и Ковалевской являются частными случаями диаграмм, построенных в работах С.Б.Каток [78], Я.В.Татаринова [191], [192]. Перейдем теперь к исследованию топологического типа 3 -многообразия Q3 для остальных интегрируемых случаев, перечисленных в параграфе 1.
5.2.5. Случай Жуковского
Будем предполагать, что параметры гамильтониана случая Жуковского
H=(S1+λ1)22A1+(S2+λ2)22A2+(S3+λ3)22A3

удовлетворяют условиям: 0<A1<A2<A3 и λ1λ2λ3eq0.
Здесь λeq0,U(R)=0. Приведенный потенциал имеет вид:
φg(R)=(g+R,λ)22AR,R.

Критические точки отображения
f2×H:S2×R3R2

определяются из условия:
gradH=μ1gradf1+μ2gradf2,f1=1,

где μ1 и μ2 — некоторые числа. Вводя вектора S=(S1,S2,S3),R=(R1,R2,R3), λ=(λ1,λ2,λ3) и матрицу A=diag(A1,A2,A3), запишем систему соотношений (2.3) в векторной форме:
A1(S+λ)=μ2R,2μ1R+μ2S=0,R,R=1.

Если μ2=0, то μ1=0, так как Req0. Получаем следующее решение системы (2.4):
S=λ,R,R=1,μ1=μ2=0.

Множество (2.5) является двумерной сферой в R6(S,R). В точках этой сферы значение гамильтониана равно нулю, а значение функции f2 равно g=λ,R. Следовательно, образ множества (2.5) при отображении f2×H представляет собой отрезок {h=0,|g|λ,λ12} на координатной оси h=0 в плоскости R2(g,h). Прообраз каждой внутренней точки этого отрезка является окружностью, минимальной для функции
H~=H|{f1=1,f2=g}, где |g|λ,λ.

Пусть теперь в системе (2.4) выполнено условие μ2eq0. Тогда
S=2μ1Rμ21.

Поскольку S,R=g, то из условия (2.6) получаем, что
g=S,R=2μ1Rμ21,R=2μ1μ21,

откуда S=gR. Подставляя в первое уравнение системы (2.4) соотношение S=gR, получаем
(μ2Ag)R=λ.

Учитывая (2.6) и (2.7), решение системы (2.4) можно записать в следующем виде:
Si=λigAitg,Ri=λiAitg.(i=1,2,3),μ1=gt2,μ2=t,

где t — некоторый параметр, а зависимость g(t) определяется уравнением
λ12(A1tg)2+λ22(A2tg)2+λ32(A3tg)2=1.

Гамильтониан H в точках, задаваемых уравнением (2.8), принимает значения
h=t22(λ12A1(A1tg)2+λ22A2(A2tg)2+λ32A3(A3tg)2).

Таким образом, для любых значений t и g, удовлетворяющих условию (2.9), мы получаем ровно одну точку (2.8), в которой градиенты функций f2 и H зависимы. Образ этой точки при отображении f2×H — это точка с координатами (g,h) из плоскости R2, где число h определяется формулой (2.10). Итак, нам надо построить кривую (g(t),h(t)) на плоскости R2(g,h), заданную неявно формулами (2.9), (2.10).

Из леммы 5.2 следует, что dhdt=tdgdt. Поэтому достаточно построить кривую (2.9) на Рис. 5.15 плоскости R2(t,g). Вводя полярные координаты t=rcosφ,g=rsinφ, из соотношения (2.9) получаем зависимость r от φ :
r=(λ12(A12+1)sin2(φφ1)+λ22(A22+1)sin2(φφ2)+λ32(A32+1)sin2(φφ3))12,

где tgφi=Ai(i=1,2,3). Кривая, определяемая этим соотношением, изображена на рис. 5.15. Она центрально симметрична относительно начала координат и имеет асимптоты g=Ait±λi, отмеченные на рис. 5.15 пунктиром.

Формула (2.10) задает отображение этой кривой в плоскость R2(g,h). Учитывая соотношение dhdt=tdgdt, легко понять, что локальным минимумам и максимумам функции g(t) (см. рис. 5.15) соответствуют точки возврата бифуркационной кривой, а локальным минимумам и максимумам обратной функции t(g) отвечают точки перегиба бифуркационной кривой. Если точка, двигаясь по кривой, изображенной на рис. 5.15, асимптотически стремится при t к прямой g=Ait±λi, то образ этой точки на бифуркационной кривой асимптотически стремится к параболе
h=(gλi)22Ai+12(λj2Aj(AjAi)2+λk2Ak(AkAi)2),

где {i,j,k}={1,2,3}.

Бифуркационная диаграмма отображения f2×H для гамильтониана случая Жуковского представляет собой объединение построенной кривой с отрезком
{h=0,|g|λ,λ12}.

Она изображена на рис.5.16. Диаграмма симметрична относительно прямой g=0. Точки касания отрезка с кривой имеют координаты
(±λ12+λ22+λ32,0)R2(g,h),

а точка пересечения двух ветвей кривой имеет координаты:
(0,λ122A1+λ222A2+λ322A3)R2(g,h).

Рис. 5.16
Топологический тип многообразия Qh3 в каждой из областей можно определить с помощью предложения 5.2 , рассматривая проекцию Qh3 на сферу Пуассона. То же самое можно получить и другим способом. Заменим λi на λi=αλi и устремим α к нулю. В пределе бифуркационная диаграмма, изображенная на рис. 5.16 , переходит в бифуркационную диаграмму для гамильтониана случая Эйлера, показанную на рис. 5.3. При этом неограниченные области на рис. 5.16, для которых 3-многообразие Qh3 гомеоморфно 2S3,S1×S2,RP3, переходят в соответствующие области на рис. 5.3. В треугольнике около начала координат на рис.5.16 многообразия Qh3 гомеоморфны S1×S2, а в граничащих в ним областях — это S3. Эти утверждения следуют из теории Морса-Ботта, поскольку в прообразе каждой точки отрезка лежит ровно одна критическая окружность, на которой функция H~=H|{f1=1,f2=g} достигает минимума, а в прообразе других точек границы образа TS2 — ровно одна критическая точка, на которой эта же функция тоже достигает минимума.
5.2.6. Случай Сретенского
Рассмотрим теперь гамильтониан
H=12(S12+S22+4(S3+λ)2)+R1.

Дополнительный интеграл для него существует лишь на одной 4-поверхности {f1=1,f2=0}R6(S,R). Поэтому для описания инвариантов этого интегрируемого случая нет необходимости изучать отображение f2×H. Топологический тип Q3 и молекула W, т.е. инвариант Фоменко, для Q3 будут в данном случае зависеть от значения параметра λ и значения h, определяющего 3 -поверхность Qh3={f1=1,f2=0,H=h}. Поэтому для гамильтониана (2.11) мы построим разделяющие кривые на плоскости R2(λ,h) и определим топологический тип многообразия Q3 в каждой из областей.

Топологический тип Q3 определим, проектируя Q3 на сферу Пуассона. Для гамильтониана (2.11) функция (2.2) имеет вид:
φ0(R)=2λ2R3243R32+R1.

Для того, чтобы определить вид областей
{φ0(R)h},

найдем критические точки функции (2.12), заданной на сфере {R12+R22+R32= =1}R3(R). Они определяются системой уравнений
φ0Ri=2μRi(i=1,2,3)R12+R22+R32=1.

Ее решениями являются две критические точки, существующие при любом значении λ :
R1=±1,R2=R3=0,μ=±12,

а также две критические точки, зависящие от значения λ :
R1=12ξ,R2=0,R3=±114ξ2,μ=ξ,

где λ2=(4ξ2+3)2128ξ3,ξ12.
В точках (2.14) значение функции (2.12) равно φ0(R)=±1. Так как число λ произвольно, то мы получаем на плоскости R2(λ,h) две прямые h=±1, разделяющие области с различным топологическим типом Q3. В точках (2.15) значение функции (2.12) равно
φ0(R)=16ξ4+40ξ2364ξ3.

Таким образом, эти точки определяют на плоскости R2(λ,h) разделяющую кривую, заданную параметрически:
h=16ξ4+40ξ2364ξ3,λ2=(4ξ2+3)2128ξ3,ξ12.

Объединяя прямые h=±1 и кривую (2.16), получаем бифуркационную диаграмму на плоскости R2(λ,h), изображенную на рис.5.17. Она симметрична относительно прямой λ=0. Координаты точек возврата кривой (2.16) таковы:
(±13,79)R2(λ,h).

Координаты точек касания прямой h=1 и кривой таковы:
(±1,1)R2(λ,h).

Кривая (2.16) асимптотически стремится к параболе h=λ2.

Каждой точке прямых h=±1 соответствует ровно одна критическая точка функции (2.12), а каждой точке кривой — ровно две критические точки. Вычислив их индексы, получим: для прямой h=1 индекс равен 0 . Для пря-
Рис. 5.17 мой h=1 индекс равен 2 на отрезке между точками касания с кривой, и индекс равен 1 на остальной части прямой. Для кривой индекс равен 1 на участках между точкой возврата и точкой касания с прямой, и индекс равен 2 на остальной части кривой. Зная индексы критических точек, несложно установить вид областей (2.13) на сфере Пуассона. Это, соответственно:
пустое множество ,
диск D2,
кольцо S1×R1, т.е. диск D2 с одной дыркой,
диск D2 с двумя дырками,
затем — вся сфера S2.
Соответствующие им 3 -поверхности Q3 указаны на рис. 5.17. Через K3 здесь обозначена та же поверхность, что и для гамильтониана (2.1).
5.2.7. Случай Клебша
Построим бифуркационные диаграммы отображения f1×H для
H=S122A1+S222A2+S322A3+ε2(A1R12+A2R22+A3R32),

где 0<A1<A2<A3,ε=±1. Условие
gradH=μ1gradf1+μ2gradf2,

определяющее критические точки отображения f2×H:S2×R3R2, можно записать в виде:
(A1100μ2000A2100μ2000A3100μ2μ200εA12μ1000μ200εA22μ1000μ200εA32μ1)(S1S2S3R1R2R3)=0.

Так как R12+R22+R32=1, то равенство (2.17) верно лишь в том случае, если определитель матрицы, входящей в это равенство, равен нулю. Матрица распадается на три блока размера 2×2, определители которых равны
Di=ε2μ1Ai1μ22(i=1,2,3).

Если Dieq0, то Si=Ri=0. Поэтому, если выполнено лишь одно из условий Di=0, то, учитывая равенство R12+R22+R32=1, мы получаем следующие критические точки:
(t1,0,0,±1,0,0),(0,t2,0,0,±1,0),(0,0,t3,0,0,±1) из R6(S,R),

где t1,t2,t3 — некоторые параметры.
Вычисляя значения функций f2 и H в этих точках, мы получаем на плоскости R2(g,h) три параболы
h=g22Ai+εAi2(i=1,2,3).

В прообразе каждой точки парабол (2.19) лежат ровно две критические точки отображения f2×H. Так как AieqAj при ieqj, то в случае ε=1 никакие два из определителей (2.18) не равны нулю одновременно. Поэтому для ε=1 мы получаем бифуркационную диаграмму, состонщую из трех непересекающихся парабол (2.19). Она изображена на рис.5.18.

Рассмотрим случай ε=1. Тогда из равенства нулю любых двух определителей (2.18) следует, что μ1=0,μ22=1. Поэтому равен нулю и третий определитель. Из (2.17) получаем следующую систему уравнений:
Si=±AiRi(i=1,2,3),R12+R22+R32=1.

Уравнения (2.20) определяют в R6(S,R) две двумерные сферы, целиком заполненные критическими точками отображений f2×H. Значения функций f2 и H на множестве (2.20) равны
g=±(A1R12+A2R22+A3R32) и h=A1R12+A2R22+A3R32.

Таким образом, при отображении f2×H двумерные сферы (2.20) отображаются в два отрезка
{h=|g|,A1|g|A3}R2(g,h).

Прообраз каждой точки этих отрезков — это две окружности, заполненные критическими точками.
Объединяя параболы (2.19) при ε=1 и отрезки (2.21), получаем бифуркационную диаграмму, изображенную на рис. 5.19. Отрезки (2.21) касаются всех трех парабол в точках (±Ai,Ai)R2(g,h). Параболы (2.19) пересекаются в точках
(±AiAj,Ai+Aj2)R2(g,h).

Топологический тип многообразия Qh3= ={f1=1,f2=0,H=h} в каждой из областей на рис. 5.18 и рис. 5.19 можно определить, рассматривая проекцию Qh3 на сферу Пуассона. Образ Qh3 при проекции определяется условием
g2z1+εz2h,A1zA3,

где через z обозначено выражение A1R12+A2R22+A3R32. График функции φ(z)=g2z1+εz изображен на рис. 5.20a при ε=1 и на рис. 5.20 b при ε=1. В случае ε=1 функция φ(z) имеет минимум при z=g. На рис. 5.20 b приведен один из возможных вариантов расположения g на оси z относительно A1,A2,A3. Несложно исследовать все возможные варианты и установить, какой вид имеет область (2.22) на сфере Пуассона в каждом случае. После этого топологический тип Q3 определяется стандартным способом, при помощи предложения 5.2. Список всех возникающих здесь 3 -поверхностей Q см. на рис.5.19. Через N3 на рис. 5.19 обозначено многообразие (S1×S2)#(S1×S2)#(S1×S2). Его проекция на сферу Пуассона — это диск с тремя дырками.
5.2.8. Случай Стеклова

Построение бифуркационных диаграмм отображения f2×H для случая Стеклова можно провести совершенно аналогично случаю Клебша. Как и в случае Клебша, бифуркационная диаграмма здесь является объединением трех парабол

Рис. 5.20
и отрезков общих касательных к ним. Но оси этих парабол, в отличие от случая Клебша, уже не совпадают. В результате получается большое количество различных вариантов расположения этих парабол в зависимости от соотношений между параметрами A1,A2,A3 гамильтониана. Описание этого случая, аналогичное предыдущим, становится уже достаточно громоздким. По этой причине мы его здесь не приводим. Все возможные молекулы W для случая Стеклова будут вычислены ниже.

1
Оглавление
email@scask.ru