УДК 517+531.01
М.В. Шамолин1
СЛУЧАИ ИНТЕГРИРУЕМОСТИ, СООТВЕТСТВУЮЩИЕ ДВИЖЕНИЮ МАЯТНИКА В ЧЕТЫРЕХМЕРНОМ ПРОСТРАНСТВЕ2
В работе систематизируются результаты по исследованию уравнений движения динамически симметричного закрепленного четырехмерного твердого тела-маятника, находящегося в некотором неконсервативном поле сил. Его вид заимствован из динамики реальных закрепленных твердых тел, помещенных в однородный поток набегающей среды. Параллельно рассматривается задача о движении свободного четырехмерного твердого тела, также находящегося в подобном поле сил. При этом на данное свободное тело действует также неконсервативная следящая сила либо заставляющая во все время движения величину скорости некоторой характерной точки твердого тела оставаться постоянной во времени (что означает наличие в системе неинтегрируемой сервосвязи). Указаны нетривиальные механические и топологические аналогии.
Ключевые слова: четырехмерное твердое тело, неконсервативное поле сил, динамическая система, случаи интегрируемости.
1. Модельные предположения
Рассмотрим однородный трехмерный круговой диск V3 с центром в точке D, гиперплоскость которого в четырехмерном евклидовом пространстве E4 перпендикулярна державке OD. Диск жестко закреплен к державке, находящейся на (обобщенном) сферическом шарнире O, и обтекается однородным потоком среды. В этом случае тело представляет собой физический (обобщенный сферический) маятник. Поток среды движется из бесконечности с постоянной скоростью v = = 0, а державка сопротивления не создает (ср. с [1; 2]).
Предположим, что суммарная сила S воздействия потока среды на диск перпендикулярна диску 3, а точка N приложения этой силы определяется, по крайней мере, углом атаки а, измеряемым между вектором скорости vd точки D относительно потока и державкой OD, углами въ/®2, измеряемыми в гиперплоскости диска D3 (таким образом, (v,a, /31, (32) — (обобщенные) сферические координаты конца вектора vD), а также тензором приведенной угловой скорости w = lQ/vd, vd = |vd | (l — длина державки, Q — тензор угловой скорости маятника). Подобные условия обобщают модель струйного обтекания пространственных тел [3; 4].
Вектор e = OD/l определяет ориентацию державки. Тогда S = s(a)v'2De, где s(a) = si(a)sign cos а, где коэффициент сопротивления si ^ 0 зависит лишь от угла атаки а. В силу свойств осевой симметрии тела-маятника относительно точки D функция s(a) является четной.
Пусть DX1X2X3X4 — система координат, жестко связанная с телом, при этом ось Dxi имеет направляющий вектор e, а оси Dx2, DX3 и DX4 лежат в гиперплоскости диска D3.
Углами (С, П1, П2) мы определим положение державки OD в четырехмерном пространстве E4. При этом угол С будем измерять между державкой и направлением набегающего потока. Другими словами, вводимые углы являются (обобщенными) сферическими координатами точки D центра диска на трехмерной сфере постоянного радиуса OD.
Пространством положений такого (обобщенного) сферического (физического) маятника является трехмерная сфера
S3{(C,ni,т) е R3 : 0 < £,щ < п, П2 mod 2п}, (1.1)
а фазовым пространством — касательное расслоение трехмерной сферы
TS3{(£,rj 1, П2; С,т,т) е R6 : о < С,т < п, П2 mod 2п}. (1.2)
!© Шамолин М.В., 2017
Шамолин Максим Владимирович (shamolin@rambler.ru, shamolin@imec.msu.ru), Институт механики Московского государственного университета им. М. В. Ломоносова, 119192, Российская Федерация, г. Москва, Мичуринский пр., 1.
2Работа выполнена при финансовой поддержке РФФИ (грант 15-01-00848-а).
Тензор (второго ранга) П угловой скорости в системе координат 0x1x2x3x4 будем определять через косо-симметрическую матрицу
П =
/ 0 —ш6 ^5 — ^3 \
ш6 0 —Ш4
—^5 Ш4 0 — Ш1
\ и3 —Ш2 Ш1 0 )
П е 8о(4).
(1.3)
Расстояние от центра О диска V3 до центра давления (точки N) будет иметь вид | = г^ = = DN (а, ¡3\, ^2,1П/ь о), где т^ = {0,X2N ,x3м ,x4м} в системе 0x1x2x3x4 (волну над П опустим).
Сразу же заметим, что, также как и в двумерном и трехмерном случаях, используемая модель воздействия потока среды на закрепленный маятник аналогична построенной модели для свободного тела и в дальнейшем учитывает влияние вращательной производной момента силы воздействия среды по тензору угловой скорости маятника (см. также [5; 6]). Анализ задачи об обобщенном сферическом (физическом) маятнике в потоке позволит обнаружить качественные аналогии в динамике частично закрепленных и свободных четырехмерных тел.
2. Группа динамических уравнений на алгебре Ли 8с(4)
Пусть четырехмерное твердое тело © массы т с гладкой трехмерной границей д© находится под воздействием некоторого неконсервативного поля сил (а именно, это можно интерпретировать как движение тела в сопротивляющейся среде, заполняющей четырехмерную область евклидового пространства Е4). Предположим, что оно является динамически симметричным, при этом имеются две логические возможности представления его тензора инерции в случае наличия двух независимых равенств главных моментов инерции: либо в некоторой связанной с телом системе координат 0x1x2xзx4 оператор инерции имеет вид
&а8{/ь/2,/2,/2}, (2.1)
либо вид diag{/l,/1,/3,/3}. В первом случае в гиперплоскости 0x2x3x4 тело динамически симметрично (другими словами, ось 0x1 — ось динамической симметрии тела), а во втором случае двумерные плоскости 0x1 X2 и 0x3X4 являются плоскостями динамической симметрии тела.
Конфигурационным пространством свободного п-мерного твердого тела является прямое произведение пространства И" (определяющего координаты центра масс тела) на группу его вращений БО(п) (определяющую вращение тела вокруг центра масс) И" х БО(п) и имеет размерность п + п(п — 1)/2 = п(п + 1)/2. Соответственно, размерность фазового пространства равна п(п + 1).
В частности, если П — тензор угловой скорости четырехмерного твердого тела (а он является терзором второго ранга [7; 8; 9]), П е бо(4), то та часть динамических уравнений движения, которая отвечает алгебре Ли эо(4), имеет следующий вид [8; 9]:
ПЛ + ЛП + [П, ПЛ + ЛП]= М, (2.2)
Л = diag{Лl,Л2,Лз,Л4 },
—/1 + /2 + /3 + /4 л /1 — /2 + /3 + /4 Л2
22
/1 + /2 — /3 + /4 л /1 + /2 + /3 — /4
Л3 - - , Л4 - - ,
3 2 ' 4 2 '
М = Мр — момент внешних сил Е, действующих на тело в И4, спроектированный на естественные координаты в алгебре Ли во(4), [.,.] — коммутатор в во(4). Кососимметрическую матрицу (соответствующую данному тензору второго ранга) П е эо(4) будем представлять в виде (1.3), где ш1, и2, и3, и4, и5, и6 — компоненты тензора угловой скорости в проекциях на координаты в алгебре Ли эо(4). При этом, очевидно, выполнены следующие равенства: Л^ — Л^ = /^ — /^ для любых г,] = 1,..., 4.
При вычислении момента внешней силы, действующей на тело, необходимо построить отображение И4 х х И4 —> бо(4), переводящее пару векторов (В^ Е) е И4 х И4 из И4 х И4 в некоторый элемент из алгебры Ли эо(4), где DN = {0^^^4м}, Е = Е — внешняя сила, действующая на тело. При
этом строится соответствующая вспомогательная матрица
0 X2N X3N X4N
.
Тогда правая часть системы (2.2) примет вид
М = {М1, М2, М3, М4, М5, М6} =
= ^4 — X4N ^3 ,x4м ^2 — X2N F4, — X4N Fl,X2N Fз — X3N F2,X3N Fl, —X2N Fl},
( 0 -Мб М5 -М3
Мб 0 - М4 М2
-M5 М4 0 -М1
\ М3 -М2 М1 0
где М2, М3, М4, М5, Мб — компоненты тензора момента внешней силы в проекциях на координаты в алгебре Ли во(4),
М
В нашем случае закрепленного маятника реализуется случай (2.1). Тогда динамическая часть уравнений его движения примет следующий вид:
/ + 12)^1 =0, (II + 12)^2 =0, 212^3 + (/1 - /2+ = Х4М (а, @1, @2, ¡(а^2,
(II + /2)^4 = 0, (2.3)
2/2^5 + (/1 - 12)(шАш6 - Ш1Ш3) = -Х3И (а,в1,в2, fi/v) 8(аУи2, 2/2^6 + (/2 - /1)(ш4^5 + Ш2Ш3) = Х2м (а, @1, @2, ПМ 8(а^2,
поскольку момент силы воздействия среды определяется через следующую вспомогательную матрицу:
(
0 X2N X3N X4N
-s(a)v 2 0 0 0
где {—s(a)v2D,0, 0,0} — разложение силы S воздействия среды в системе координат DX1X2X3X4.
Поскольку размерность алгебры Ли so(4) равна 6, система уравнений (2.3) и составляет группу динамических уравнений на so(4).
Видно, что в правую часть системы уравнений (2.3) входят, прежде всего, углы a,@i,@2, поэтому данная система уравнений не является замкнутой. Для того, чтобы получить полную систему уравнений движения маятника, необходимо к динамическим уравнениям на алгебре Ли so(4) присоединить несколько групп кинематических уравнений.
2.1. Циклические первые интегралы
Сразу же заметим, что система (2.3), в силу имеющейся динамической симметрии
I2 = I3 = I4, (2.4)
обладает тремя циклическими первыми интегралами
= ш0 = const, ш2 = ш0 = const, = и0 = const. (2.5)
При этом в дальнейшем будем рассматривать динамику системы на нулевых уровнях:
,0 = = ш1
При условиях (2.4)—(2.6) система (2.3) примет вид незамкнутой системы трех уравнений:
2^3 = Х4м (а, @2, Q/v) 8(а^2, 2^5 = -хзм (а,в1,в2, tt/v) 8(а^2,
(2.6)
(2.7)
2I2Ú6 = X2N (а,@1,в2, &/v) s(a)v .
3. Первая группа кинематических уравнений
Для получения полной системы уравнений движения нам потребуется группа кинематических уравнений, связывающих скорости точки D (центра диска D3) и набегающего потока:
v D = vD • iv(a,pup2) = ñ 1 + (-vTO)iv(-tm,m), 1 = {1,0,0,0}, (3.1)
I cos a \
sin a cos Pi sin a sin в1 cos в2 У sin a sin в1 sin в2 J
Равенство (3.1) выражает теорему сложения скоростей в проекциях на связанную систему координат
DX1X2X3X4.
Действительно, в левой части равенства (3.1) стоит скорость точки D маятника относительно потока в проекциях на связанную с маятником систему координат DX1X2X3X4. При этом вектор iv(a,@i,@2) — единичный вектор вдоль оси вектора v d . Вектор iv (a, ,#1,^2) имеет (обобщенные) сферические координаты (1,a,el,в2), определяющие разложение (3.2).
iv (a,Pi,fi2) =
(3.2)
)
В правой части равенства (3.1) стоит сумма скоростей точки D при повороте маятника (первое слагаемое) и движения потока (второе слагаемое). При этом в первом слагаемом имеются координаты вектора l = OD = = {l, 0,0, 0} в системе координат Dxix2x3x4.
На втором слагаемом правой части равенства (3.1) остановимся подробнее. В нем имеются координаты вектора (—vTO) = { v^q, 0,0,0} в неподвижном пространстве. Чтобы его записать в проекциях на связанную систему координат Dxix2x3x4 необходимо произвести (обратный) поворот маятника на угол (—£), что алгебраически эквивалентно умножению величины (— vTO) на вектор iv ( — £, ni,n2).
Таким образом, первая группа кинематических уравнений (3.1) в нашем случае примет следующий вид:
vD cos а = —vcos £, vD sin a cos в = l<6 + vsin £ cos щ,
vD sin a sin в1 cos (32 = — l<5 + vsin £ sin щ cos n2, (3.3)
vD sin a sin в1 sin в2 = l<3 + vsin £ sin ni sin n2 •
4. Вторая группа кинематических уравнений
Нам также потребуется группа кинематических уравнений, связывающих тензор угловой скорости П и координаты £,П1,П2,£,П1,П2 фазового пространства (1.2) исследуемого маятника — касательного расслоения ,П1,П2; £ т,П2}.
Проведем рассуждения в стиле, допускающем любую размерность. Искомые уравнения получаются из следующих двух групп соотношений. Поскольку движение тела формально происходит в евклидовом пространстве Е",п = 4, сначала выражается набор, состоящий из фазовых переменных ^3,^5,^6, через новые переменные ^1,^2,^3 (из набора г). Для этого производится следующая композиция поворотов на углы П1,П2:
Ti,2(П2) 0 T2,3(ni)
(4.1)
T2,3(ni)
0
cos ni sin ni
0
— sin ni cos ni
Ti ,2 (П2)
Другими словами, справедливы соотношения
Zi Z2 Z3
cos n2 sin n2 0
<3
— sin n2 cos n2 0
T2,3 — ni) 0 Ti,2 ( П2) I <5
<e
т.е.
Zi = <3 cos ni + <5 sin n2, Z2 = —<¿3 cos ni sin n2 + <5 cos ni cos n2 + <6 sin ni, Z3 = <3 sin ni sin n2 — <5 sin ni cos n2 + <6 cos ni • Затем вместо группы переменных z подставляется следующая зависимость:
sin £ sin £
Z3 = £ , Z2 = —ni-т, Zi = n2-
cos £ cos £
sin щ-
(4.2)
Таким образом, две группы уравнений (4.1) и (4.2) дают вторую группу кинематических уравнений:
<3 = £ sin ni sin n2 + ni fon| cos ni sin n2 + n'2 CoS| sin ni cos n2,
<5 = —£ sin ni cos n2 — ni cos! cos ni cos n2 + n'2 C0S| sin ni sin n2, (4.3)
<e = £ cos ni — n* i CoS| sin ni-
Видно, что три группы соотношений (2.7), (3.3), (4.3) образуют замкнутую систему уравнений. В эти три группы уравнений входят следующие функции:
x2N (a,l3i,l32, — ) , x3N (a,Pi,P2, — ) , xiN (. V vD) \ vDJ V
a, fii,@2, — ) , s(a)-vD
При этом функция в считается зависимой лишь от а, а функции X2N,xзN,X4N могут зависеть, наряду с углами а, вь в2, вообще говоря, и от приведенного тензора угловой скорости 1П/ур.
5. Задача о движении свободного тела при наличии следящей силы
Параллельно рассматриваемой задаче о движении закрепленного тела, рассмотрим пространственное движение свободного динамически симметричного (случай (2.1)) четырехмерного твердого тела с передним торцом (круговым трехмерным диском Р3) в поле силы сопротивления в условиях квазистационарности [10; 11] с той же моделью воздействия среды.
Если (v,a,pl,p2) — сферические координаты вектора скорости центра D диска 3, лежащего на оси симметрии тела, — — тензор угловой скорости тела (см. (1.3)) в системе координат DX1X2X3X4, связанной с телом, при этом ось симметрии CD совпадает с осью Dxi (C — центр масс), а оси Dx2,Dx3,Dx4 лежат в гиперплоскости диска, Ii, I2, I3 = I2, I4 = I2, m — инерционно-массовые характеристики, то динамическая часть уравнений движения тела, при котором касательные силы воздействия среды на диск отсутствуют, примет вид:
V cos a — av sin a — w6v sin a cos в1 + w5v sin a sin в1 cos (32-
p
—w3v sin a sin в1 sin в2 + &(w6 + ш5 + wf) = —,
m
v sin a cos в1 + av cos a cos в1 — в iv sin a sin в1 + +шбv cos a — w4v sin a sin в1 cos в2 + w2v sin a sin в1 sin в2 — a(w4w5 + ш2ш3) — аШб = 0, v sin a sin в1 cos в2 + av cos a sin в1 cos в2 + вlv sin a cos в1 cos в2 — —e2v sin a sin в1 sin в2 — w5v cos a + w4v sin a cos в1 — wlv sin a sin в1 sin в2 — —a(—wiU2 + Ш4 Шб) + = 0, v sin a sin в1 sin в2 + av cos a sin в1 sin в2 + вlv sin a cos в1 sin в2 + (5.1)
+e2v sin a sin в1 cos в2 + w3v cos a — w2v sin a cos в1 + wiv sin a sin в1 cos в2+ +a(w2we + Ш1Ш5) — aÚ3 = 0,
(Ii + I2)Wi =0, (Ii + I2)¿5 =0,
2I2W3 + (Ii — I2)(W2W6 + wiW5) = X4N ^a, ei, в2, s(a)v2, (Ii + h)Ü4 = 0, 2I2W5 + (Ii — I2)(w4w6 — wiw3) = —X3N (a,ei, в2, — ) s(a)v2,
2I2W6 + (I2 — Ii)(w4W5 + Ш2Ш3) = X2N [a,ei,e2, — ^ s(a)v2,
где Fi = —S, S = s(a)v2, a = CD, при этом
(о, X2N (а, в1,в2, -V^j , X3N ^a, , ^ , x4N ^a, Pi, fc,
— координаты точки N приложения силы S в системе координат DxiX2X3X4, связанной с телом.
Первые четыре уравнения системы (5.1) описывают движение центра масс в четырехмерном евклидовом пространстве E4 в проекциях на систему координат DXiX2X3X4. Вторые же шесть уравнений системы (5.1) получены из (2.2).
Таким образом, фазовым пространством системы динамических уравнений (5.1) десятого порядка является прямое произведение Rl х S3 х so(4) четырехмерного многообразия на алгебру Ли so(4). При этом, поскольку сила воздействия среды не зависит от положения тела в пространстве, система динамических уравнений (5.1) отделяется от системы кинематических уравнений и может быть рассмотрена самостоятельно (см. также [12, 13]).
5.1. Циклические первые интегралы
Сразу же заметим, что система (5.1), в силу имеющейся динамической симметрии
I2 = I3 = I4, (5.2)
обладает циклическими первыми интегралами
ш1 = ш0 = const, ш2 = ш0 = const, ш4 = и0 = const. (5.3) При этом в дальнейшем будем рассматривать динамику системы на нулевых уровнях:
ш0 = ш0 = ш0 = 0. (5.4)
5.2. Неинтегрируемая связь
Если рассматривается более общая .задача о движении тела при наличии некоторой следящей силы T, проходящей через центр масс и обеспечивающей во все время движения выполнение равенства
v = const, (5.5)
то в системе (5.1) вместо Fi будет стоять величина T — s(a)v2.
В результате соответствующего выбора величины T следящей силы можно формально добиться во все время движения выполнения равенства (5.5). Действительно, формально выражая величину T в силу системы (5.1), получим при cos а = 0:
Т = Tv (а, в1, в2, —) = ma(wf + Ш + ) + s{a)v2
ma sin а
1 —--1 ^
2I2 cos а
(а,в1,в2, П)
Г„
^а,вьв2, — ^ = X4N ^
X4N ( а, Pi, Р2 ) sin в1 sin в2 +
П
+X3N а, Pi, Р2, — sin Pi cos P2 + X2N а
)si
(а,в1 ,в2, — )
— I в
— cos p1.
v
(5.6)
На данную процедуру можно посмотреть с двух позиций. Во-первых, произошло преобразование системы при помощи наличия в системе следящей силы (управления), обеспечивающей рассмотрение интересующего нас класса движений (5.5). Во-вторых, на это все можно посмотреть как на процедуру, позволяющую понизить порядок системы. Действительно, система (5.1) в результате действий порождает независимую систему шестого порядка следующего вида:
áv cos а cos в1 — в 1 v sin а sin в1 + w6v cos а — аш6 = 0, áv cos а sin в1 cos в2 + в 1v sin а cos в1 cos в2 — p2v sin а sin в1 sin в2 — w5v cos а + аШ5 = 0, аи cos а sin в1 sin в2 + в 1v sin а cos в1 sin в2 + в2 v sin а sin в1 cos в2+
+^3v cos а — ашз = 0,
(5.7)
212Ш3 = X4N а
(а,в1 ,в2, —)
— ) s(а)v2 v
212Ш5 = —X3N а
[а,в1 ,в2, — )
-v^ s(а)v2, 2I2Ш6 = X2N (уа, в1,в2, — ^ s(а)v~
в которой к постоянным параметрам, указанным выше, добавляется параметр v. Система (5.7) эквивалентна системе
аv cos а + v cos а [ш6 cos в1 — ш5 sin в1 cos в2 + ш3 sin в1 sin в2 ] +
+а [—Ш6 cos в1 + ш5 sin в1 cos в2 — Ш3 sin в1 sin в2] = 0, e1v sin а — v cos а [ш5 cos в1 cos в2 + ш6 sin в1 — ш3 cos в1 sin в2] + +а [Ш5 cos в1 cos в2 + ш6 sin в1 — Ш3 cos в1 sin в2] = 0, e2v sin а sin в1 + v cos а [ш3 cos в2 + ш5 sin в2] + а [—Ш3 cos в2 — ш5 sin в2] = 0,
(5.8)
v2 —
Ш3 = 212 X4N (а, в1,в2, v ) в(а),
Ш5 = — 2I" X3N 1 а Введем новые квазискорости в системе:
^а,в1,в2, -— ^ «(а), Шб = 2I- X2N (а
v
v2
2I2
(а,в1,в2, — )
— ) «(а).
Z1
Т1,2(в2) ◦ Т2,3(в1) | Z2 Z3
1 0 0
Т2,3(в1)= | 0 cos в1 — sin в1 0 sin в1 cos в1 Другими словами, справедливы соотношения
Т1,2(в2)
Z1 Z2 Z3
cos в2 — sin в2 0 sin в2 cos в2 0 0 0 1
Ш3
Т2,3( — в1) ◦ Т1,2(—в2) | Ш5
(5.9)
(5.10)
(5.11)
Ш6
2
v
т.е.
zi = Ш3 cos ei + w5 sin в2 , z2 = —w3 cos ei sin в2 + w5 cos ei cos в2 + шб sin ei, (5.12)
z3 = w3 sin ei sin в2 — w5 sin ei cos в2 + Шб cos ei.
Как видно из (5.8), на многообразии
Oi = {(a,ei,e2,w3,ш5,шб) G R6 : a = Пk, ¡i = nl, k,l G zj (5.13)
нельзя однозначно разрешить систему относительно а, / 1, /?2. Формально, таким образом, на многообразии (5.13) происходит нарушение теоремы единственности. Более того, при к четном и любом I неопределенность возникает по причине вырождения сферических координат а, /1, /2), а при к нечетном происходит явное нарушение теоремы единственности, поскольку при этом первое уравнение (5.8) вырождается.
Из этого следует, что система (5.7) вне и только вне многообразия (5.13) эквивалентна системе
av s(a) ( —Л
—Z3 + ^—— М a, ei, в2,-
2i2 cos a V v
v2 — 2 2 cos a
Z3 = ТГТs(a)i v a,¡íi,¡í2,— — z + z^---
2I2 \ v J sin a
av s(a) . ( „ „ av s(a) _ ( „ „
-z2Aj a,ei,e2,-) + zi©j a,ei,e2,-),
2I2 sin a \ v J 2I2 sin a \ v J
v2 — cos a 2 cos a cos в1
z2 = —wr s(a)Av a,¡íi,¡í2,— + z2 z3~--+ z^--:—— +
2I2 \ v J sin a sin a sin вi
av s(a) Л ( n n — \ av s(a) cos ¡i „ ( „ — A „
z3Aj a,ei,e2,- ~—TTzi©v [a,eue2,-), (5.14)
2I2 sin a \ v J 2I2 sin a sin ¡i \ v J
v2 — cos a cos a cos в1
zi = t~¡-s(a)©v a,ei,¡2,— + ziz3 ;--ziz2~.--——
2I2 \ v J sin a sin a sin ei
av s(a) — av s(a) cos в1 —
2I2 sin a \ v J 2I2 sin a sin ¡i \ v J
■ cos a av s(a) л ( „ „
в i = z^-+ — Av [a,ei,e2,-),
sin a 2I2 sin a \ v J
cos a av s(a) „ / — \
¡2 = — z^——- + — . .\©v (a,ei,e2,-)
sin a sin ei 2I2 sin a sin pi \ v J
^a, ¡i, ¡2, —v^ = X4N ^a, ¡i, ¡2,—— ^
Av ( a, ¡i, ¡2, -v j = X4N (a, ¡i, ¡2, — ) COS вi sin ¡2 +
+X3N (a,ei,e2, —^ cos вi cos e2 — X2N [a, ¡i, ¡2, —^ sin вi, (5.15)
vv
©v ( a, ¡i, ¡2, v I = X4N \a,ei,e2, v I cosв2 —
—X3N [a, ¡i, ¡2, —^ sin в2, (5.16)
а функция rv (a,ei,e2, — /v) представляется в виде (5.6).
Здесь и далее зависимость от групп переменных (a, ¡i, ¡2, — /v) понимается как сложная зависимость от (a, ¡i, ¡2, zi/v, z2/v, z3/v) в силу (5.12).
Нарушение теоремы единственности для системы (5.8) на многообразии (5.13) при нечетном k происходит в следующем смысле: почти через любую точку из многообразия (5.13) при нечетном k проходит неособая фазовая траектория системы (5.8), пересекая многообразие (5.13) под прямым углом, а также существует фазовая траектория, полностью совпадающая во все моменты времени с указанной точкой. Но физически это различные траектории, так как им отвечают разные значения следящей силы.
5.3. Постоянная скорость центра масс
Если рассматривается более общая .задача о движении тела при наличии некоторой следящей силы T, проходящей через центр масс и обеспечивающей во все время движения выполнение равенства
Ус = const (5.17)
(Ус — скорость центра масс), то в системе (5.1) вместо должна стоять величина, тождественно равная нулю, поскольку на тело будет действовать неконсервативная пара сил: Т — в(а)у2 = 0. Очевидно, что для этого нужно выбрать величину следящей силы Т в виде
Т = Ту(а,вьв2, = s(a)v2, Т =-8. (5.18)
Случай (5.18) выбора величины Т следящей силы является частным случаем возможности отделения независимой подсистемы четвертого порядка после некоторого преобразования системы (5.1). Действительно, пусть выполнено следующее условие на величину Т:
Т = Т;(а,в1,в2, —) = n,j (а,в1,в2, — jwjWj = Т1 (а,вьв2, — ) v2, шо = v.
i,j=0, i^j ^ / \ /
Введем для начала новые квазискорости (5.9)—(5.11). Систему (5.1) в случаях (5.2)—(5.4) можно переписать в виде
2 2 2 2 —
v + a(z-, + Z2 + Z3) cos а — а——в(а) sin а ■ rv а, вь в2, —
2I2 v
Т1 (а,в1,в2, V) v2 — s(а)v2
m
2 2 2 v2 — —v + Z3v — a(z2 + z| + z3) sin а — а-«(а) cos а ■ rv I а, в1, в2, ~ ) =
2I2 V v)
s(а)v2 — Т1 (а,в1,в2, V)
а,в1,в2, V) v2 ,
v sin а,
в 1 sin а — z2 cos а--«(а) ■ Av (а, в1, в2, — ) =0. (5.19)
v
2I2
с
2I2
в2 sin а sin в1 + z1 cos а — -"в(а) ■ ©v ( а, в1, в2, — ) =0.
v
v2 —
Ш3 = тгр-X4N а, в1,в2,— «(а), 2I2 v
v2 — v2 —
Ш5 = —— X3N а,в1,в2,— «(а), Шб = — X2N а, вьв2, — «(а). 2I2 v 2I2 v
Вводя далее новые безразмерные фазовые переменные и дифференцирование по формулам zk = n1vZk, к = 1, 2, 3, < ■ >= n1v <'>, П1 > 0, n1 = const, система (5.19) приведется к следующему виду:
v' = vЪ(а,в1,в2,Z), (5.20)
а' = —Z3 + an1(Z1 + Z| + Z|) sin а +----«(а) cos а ■ Tv (а, вь в2, n1Z) —
2I2n1
Т1 (а,вl,в2,nlZ) — «(а) .
mn1
sin а, (5.21)
Z3 = «т^ ■ Tv (а,в1,в2^) — (Z2 + Z2) — — -^Z2■ Av (а^ьв^^) +
2I2n1 sin а 2I2n1 sin а
а я ( а)
-М- ■ ©v (а, в1,в2, nxZ) — Z3 ■ Ф (а, в1, в2, Z), (5.22)
2I2 n1 sin а
/ «(а) cos а
Z2 = — ^—2 ■ Av (а,в1,в2,n1 Z) + Z2Z3---+
2í2n\ sin а
2 cos а cos в1 а «(а)
+Zt —--—— + —-Z3--■ Av {а,вl,в2,nlZ) —
sin а sin в1 2I2n1 sin а
а «(а)
— ■ ©v (а, в1,в2, nZ) — Z2 ■ Ф (а, в1, в2, Z), (5.23)
2I2 n1 sin а
Z1 = ^ ■ ©v (а, в1, в2, n1Z) + Z1Z3— — 2I2n1 sin а
cos а cos в1 а «(а)
—Z1Z2---—— — —--:-—— ■ ©v (а,вl,в2,nlZ) X
sin а sin в1 2I2n1 sin а sin в1
X [Z3 sinв1 — Z2 cosв1] — Z1 ■ Ф (а, вь в2, Z), (5.24)
cos а
m
= Z2^ + ^^ . Av (a,e1,e2,n1Z), (5.25)
sin a 212n1 sin a
в2 = -Zi . COSa в + -f- . <a) в ■ ©v (а,в1 ,e2,niZ), (5.26)
sin a sin pi -I2n1 sin a sin p1
Ф(а,вьв2,Z) = —an1(Z? + Z2 + Z3)cos a + f s(a)sin a ■ rv (a,e1,e2,n1Z) +
2i2«1
T\ (a, n1Z) - s(a)
+--cos a,, (5.27)
mn1
а функции rv {а,в1,в2, , Av (а,в1,в2, , ©v (а,в1,в2, представляются в виде (5.6), (5.15), (5.16), соответственно.
Видно, что в системе седьмого порядка (5.20)—(5.26) может быть выделена независимая подсистема шестого порядка (5.21)—(5.26), которая может быть самостоятельно рассмотрена на своем шестимерном фазовом пространстве.
В частности, при выполнении условия (5.18) только что рассмотренный прием выделения независимой подсистемы шестого порядка также возможен.
6.
6.1.
Случай зависимости момента неконсервативных сил от угловой скорости
Введение зависимости от угловой скорости
Данная работа посвящена динамике четырехмерного твердого тела в четырехмерном пространстве. Но, поскольку данный раздел посвящен исследованию случая движения при наличии зависимости момента действующих сил от тензора угловой скорости, введем такую зависимость с более общих позиций. К тому же, данная точка зрения поможет нам вводить эту зависимость и для многомерных тел.
Пусть х = (хш,х2^) — координаты точки N приложения неконсервативной силы (воздействия среды) на трехмерный диск Р3, Q = (^1,^2,^3,^4) — компоненты, не зависящие от угловой скорости. Будем вводить зависимость функций (х1ш,х2ш,хзN,х4ш) от тензора угловой скорости О лишь линейным образом, поскольку само данное введение априори не очевидно [3, 10].
Итак, примем следующую зависимость: х = Q+Д, где К = (Д1, Д2, Д3, Л4) — вектор-функция, содержащая тензор угловой скорости О. При этом зависимость функции Д от тензора угловой скорости — гироскопическая:
R
R1 \ / 0 —U6 U5
R2 1 ш6 0 —u4
Rs VD — U5 U4 0
R4 J \ Us — U2 U1
—us\ ( h-1 \ h2 hs
\h4 J
Ш2
0
/
Здесь (Л-1, Л-2, Нз, Ъ,4) — некоторые положительные параметры.
Теперь, применительно к нашей задаче, поскольку х1ш = хш = 0, то x2N = Q2 — , х3ш = Q3 +
+ Н1Ш5/ув, х4ш = Q4 — Ъ,1Ш3/ю.
Таким образом, функция гш выбирается в следующем виде (диск Р3 задается уравнением х1ш = 0):
rN
iN
iv(-,в1,в2) , h
/ 0
X2N
XSN
X4N
h1
h2
hs
h4
R(a)iN--l h,
VD
íí
0 —U6 U5 —Us \
U6 0 — U4 U2
—Us U4 0 —U1
us —U2 U1 0 )
(6.1)
(6.2)
iN
(6.3)
(см. (1.3), (3.2)). В нашем случае
0
cos j31 sin в1 cos в2 У sin j31 sin в2 J
Таким образом, выполнены равенства X2N = R(a) cos — h1U6/vD, xsn = R(a) sin cos в2 + + h^s/vD, X4N = R(a) sinв1 sinв2 — h1Us/vD, убеждающие нас о том, что в рассматриваемой системе присутствует также еще и дополнительный демпфирующий (а в некоторых областях фазового пространства и
разгоняющий) момент неконсервативной силы (т.е. присутствует зависимость момента от тензора угловой скорости).
Итак, для построения силового поля также используется пара функций R(a), s(a), информация о которых носит качественный характер. Подобно выбору аналитических функций типа Чаплыгина динамические функции s и R примем в следующем виде:
R(a) = A sin a, s(a) = B cos a, A, B > 0. (6.4)
6.2. Приведенные системы
Теорема 6..1 Совместные уравнения (2.3), (3.3), (4-3) при выполнении условий (2.4)-(2.6), (6.1), (6.4) редуцируются к динамической системе на касательном расслоении (1.2) трехмерной сферы (1.1). Действительно, если ввести безразмерные параметры и дифференцирование:
AB h\B
b* = lno, n0 = —p-, Hi* = —-, < • >= novж <'>, (6.5)
¿12 ¿hno
то полученные уравнения будут иметь следующий вид:
+ (b* - Hi*)£' cos e + sin e cos e - [n? + n22 sin2 ni] COé, =0, < + (b* - Hi*)n' cos e + e'n'l С+ОЁг - n22 sin ni cos ni = 0, (6.6)
n2 + (b* - Hi *)n2 cos e + e 'ric+oiá + ¿n = 0, b* > 0, Hi * > 0.
После же перехода от переменных z (о переменных z см. (4.2)) к промежуточным безразмерным переменным zk = n0vTO(1 + b*Hi*)Zk, k = 1, 2, z3 = n0vTO(1 + b*Hi*)Z3 -n0v^b* sine, система (6.6) будет эквивалентна системе
e ' = (1 + b*Hi*)Zs - b* sine, (6.7)
Z'3 = - sinecose + (1 + b*Hi*)(Z2 + Z22)cos| + Hi*Z3 cos e, (6.8)
sin e
Z2 = -(1 + b*Hi*)Z2Z3cos| - (1 + b*Hi*)Z2cos|^^ + Hi*Z2 cose, (6.9) 2 sin e sin e sin n
Z[ = -(1 + b*Hi*)ZiZ3cos| + (1 + b*Hi*)ZiZ2cos|^^ + Hi*Zi cose, (6.10) sin e sin e sin n
nl = -(1 + b*Hi *)Z2, (6.11) sin e
n2 = (1 + b*Hi*)Zi . cose , (6.12) 2 sin e sin n
на касательном расслоении
T*S3{(Z3,Z2,Zi;e,ni,m) e R6 : 0 < e,ni < n, П2 mod 2n} (6.13)
трехмерной сферы S3{(e,ni,n2) e R3 : 0 < e,ni ^ n, n2 mod 2n}.
Видно, что в системе шестого порядка (6.7)—(6.12) по причине цикличности переменной П2 выделяется независимая подсистема пятого порядка (6.7)—(6.11), которая может быть самостоятельно рассмотрена на своем пятимерном многообразии.
6.3. Полный список первых интегралов
Перейдем теперь к интегрированию искомой системы шестого порядка (6.7)—(6.12) (без всяких упрощений — при наличии всех коэффициентов).
Аналогичным образом, для полного интегрирования системы (6.7)—(6.12) шестого порядка необходимо знать, вообще говоря, пять независимых первых интегралов. Однако после замены переменных
+ z2, w i = Z2
'п
система (6.7)-(6.12) распадается следующим образом:
e' = -(1 + ь* Hi *)w3 - ь* sin e, 1
w3 = sinecos e - (1 + b*Hi+ Hi*W3 cos e,(6.15) w2 = (1 + b*Hi*)W2W3 0 + Hi*W2 cose,
w1 = di(w3,w2,wi; e, ni,m)1 + Wi cosПi, n'i = di(w3,w2,wi; e,ni,m), (6.16)
w sin n
I--Z
W3 = -Z3, W2 = ^Z2 + Z2, wi = Z2, (6.14)
П2 = ¿2(ws,w2,w1; £,П1,П2), (6.17)
d1 (ws,w2,w1; е,П1,П2) = —(1 + b,Hu)Z2(ws,W2,W1)fOf| = T-^ Ш,
¿2(ws ,w2,w\; £,П1,П2) = (6.18)
= (1 + hHu )Z1(ws,w2,w1) . 1os.e = . |os.e ,
^ 1 * ^ s' 2' ^ sin | sin ni \/1+wf sln Isln ni
при этом Zk = Zk (ws ,w2,w1), k = 1, 2, 3, — функции в силу замены (6.14).
Видно, что в системе шестого порядка (6.15)-(6.17) выделяется независимая подсистема третьего порядка (6.15), которая может быть самостоятельно рассмотрена на своем трехмерном многообразии, независимая система второго порядка (6.16) (после замены независимой переменной), а также (по причине цикличности переменной п2) уравнение (6.17) на п2 отделяется.
Таким образом, для полной интегрируемости системы (6.15)-(6.17) достаточно указать два независимых первых интеграла системы (6.15), один — для системы (6.16) и дополнительный первый интеграл, "привязывающий" уравнение (6.17) (т.е. всего четыре).
Для начала сопоставим системе третьего порядка (6.15) неавтономную систему второго порядка
dw3 _ sln I cos cos I/ sln I+Hi„ W3 cos I
dl = -(1+btH1t)w3-bt sin I ,
dw2 _ (1 + 6,Hi,)w2W3 cos I/ sln I + HitW2 cos I
dI —(1+bt Hit )w'3—bt sln I '
Используя замену т = sin C, перепишем ее в алгебраическом виде
dws т — (1 + b*H1*)w2/T + Huws dw2 (1 + b*H1*)w2ws/T + H1*w2
(6.19)
¿т — (1 + — ЬФ т ' ¿т — (1 + — Ь*т
Далее, вводя однородные переменные по формулам ^3 = и2т, ^2 = «1т, приводим систему (6.19) к следующему виду:
т и = (1+Ь»Нц,)(«2—"2) + (Ь* + Нц,)«2 + 1 т и = 2(1+Ь»Нц)цщ2 + (Ь»+Нц)щ ^ —(1+Ь»Н1»)«2-Ь» , ^ —(1+Ь»Н1»)«2-Ь» '
Сопоставим системе второго порядка (6.20) неавтономное уравнение первого порядка
¿«2 = 1 — (1 + ЬфН1ф)(и2 — «2) + Ъ + Ни)«2 ¿«1 2(1 + ЬфН1Ф)«1«2 + (ЬФ + Ни)«1
которое несложно приводится к полному дифференциалу:
¿ ((1 + ьФН1Ф)(«2 + «1) + Ъ + Н1Ф)«2 +1 \ = 0
(6.20)
(6.21)
Итак, уравнение (6.21) имеет следующий первый интеграл:
(1 + b*H1*)(w| + ul) + (b* + H1*)u2 + 1
= C1 = const, (6.22)
u1
который в прежних переменных выглядит как
©1(ws, w2; C)
„2 + w2) + (b + H )w„ C + s;n2 C
C = const' (6.23)
(1 + b*H1*)(w2 + w2) + (b* + H1*)ws sinC + sin2 C
w2 sin C
Замечание 6..1 Рассмотрим систему (6.15) с переменной диссипацией с нулевым средним, становящейся консервативной при b* = H1* :
C' = —(1 + b*)ws — b* sinC, w's = sinCcosC — (1 + b*)w2coI + b*ws cosC, (6 24)
w2 = (1 + b2*)w2wsconI + b*w2 cosC'
Она обладает двумя аналитическими первыми интегралами вида
(1 + b2*)(w| + w2)+2b*ws sinC + sin2 C = C* = const, (6.25)
w2 sin C = C2* = const. (6.26)
Очевидно, что отношение двух первых интегралов (6.25), (6.26) также является первым интегралом системы (6.24). Но при b* = H1* каждая из функций
(1 + b*H1*)(w| + w2) + (b* + H1*)ws sinC + sin2 C (6.27)
и (6.26) по отдельности не является первым интегралом системы (6.15). Однако отношение функций (6.27), (6.26) является первым интегралом системы (6.15) при любых b*,H1*.
Далее, найдем явный вид дополнительного первого интеграла системы третьего порядка (6.15). Для этого преобразуем для начала инвариантное соотношение (6.22) при и1 = 0 следующим образом:
2 / ^ \ 2
U2 +
b* + Hi *
+
+ u i -
Ci
(b* - Hi*)2 + C2 - 4
2(1 + Ь* Н *)) 'У' 2(1 + ЬН1 *)) 4(1 + Ъ*Н\ *)2
Видно, что параметры данного инвариантного соотношения должны удовлетворять условию
(Ь* - И*)2 + С2 - 4 > 0,
(6.28)
(6.29)
и фазовое пространство системы (6.15) расслаивается на семейство поверхностей, задаваемых равенством (6.28).
Таким образом, в силу соотношения (6.22) первое уравнение системы (6.20) примет вид
<1и2 = 2(1 + Ь*Нг*)и22 + 2(Ь* + И*)и2 +2 - СгПг(Сги) ¿т -Ь* - (1 + Ь*Ни)и2
Ui(Ci,U2)
1
-{Ci ± Ü2(Ci,U2)},
2(1 + Ь*Нг*)
и2(С1,и2) = ^С2 - 4(1 + Ь*Н1*)(1 + (Ь* + Н1*)и2 + (1 + Ь*Н1*)и2),
при этом постоянная интегрирования С1 выбирается из условия (6.29).
Поэтому квадратура для поиска дополнительного первого интеграла системы (6.15) примет вид
[ ¿т Г (-Ь* - (1 + Ь*Н1*)и2)<и2
*
т J Qi
2
(6.30)
Qi = 2(1 + (b* + Hu)u2 + (1 + b*Hi*)u2) - Ci{Ci ± U2(Ci,u2)}/(2(1 + b*Hi*)). Левая часть (с точностью до аддитивной постоянной), очевидно, равна ln | sinЕсли
b* + Hi*
U2 +
= ri, b2 = (b* - Hi*)2 + C2 - 4,
2(1 + b*Hi*) то правая часть равенства (6.30) примет вид
1 í d(b2 - 4(1 + b*Hi*)r2)
4 J (b2 - 4(1 + b*Hi*)r2) ± Ci^Jb2 - 4(1 + b*Hi*)r2
dri
+
+№. -H„)(i + №,)/
-
(b2 - 4(1 + b*Hi*)r2) ± Ciy/bi - 4(1 + b*Hi*)r2 Vхb2i - 4(1 + b*Hi*)r2
Ci
1
-b* + Hi* ±--ii,
Ii
dr3
vw-4 (гз ± Ci)
, гз = ^bi - 4(1 + b*Hi*)r2.
(6.31)
(6.32)
При вычислении интеграла (6.32) возможны три случая. I. Ib* - Hi* I > 2.
/1 = - 1
V(b* - Hi*)2 - 4 '
x In
V(b* - Hi*)2 - 4 +
гз ± ci
+
±
Ci
V(b* - Hi*)2 - 4
+
V(b* - Hi*)2 - 4'
x In
V(b* - Hi*)2 - 4 -^/Ь-Л
i - r3
II. \b* - Hi* \ < 2.
III. \b* - Hi*\ = 2.
гз ± ci 1
T
Ci
V(b* - Hi*)2 - 4 . ±Cir3 + b2i
+ const.
V4 - (b* - Hi*)2 bi(r3 ± Ci)
+ const.
Ii = T
Vb—2
Ci(r3 ± Ci)
+ const.
(6.33)
(6.34)
(6.35)
1
1
3
Возвращаясь к переменной
= ^ , ь* + Н1* (6 36)
Г1 = в1пе + 2(1 + Ь*Ни), (6.36) имеем окончательный вид для величины /1:
I. \Ь* — Ни \ > 2.
/■ = — 1
2V(b* - Hu)2 - 4 '
х ln
y/(b« - Hb)2 - 4 ± 2(1 + b«Hi«)ri C\
Vхb2 - 4(1 + 6H1OM ± Ci y/(b« - Hu)2 - 4 1
+
+
V(b* - Hu)2 - 4 '
x ln
II. \b« - Hu \ < 2.
■sj(b, - Hu)2 - 4 t 2(1 + b«Hu)ri ^ Ci
vb - 4(1 + bHu)2r2 ± Ci V(b« - Hu)2 - 4
. iC^Vbi-iö^b^^T2 + b1
+ const. (6.37)
I1 = —. arcsin-. „ —---+ const. (6.38)
V4 - (b« - H1«)2 bi^b2 - 4(1 + b«Hi«)2r2 ± Ci) v ;
III. \b« - Hu\ = 2.
r ^ 2(1 + b«Hi« )ri t
Ii = T-, 0 =--+ const. (6.39)
^Ci(Vb2 - 4(1 + b«Hu)2r'2 ± Ci) v ;
Итак, только что был найден дополнительный первый интеграл для системы третьего порядка (6.15) — предъявлен полный набор первых интегралов, являющихся трансцендентными функциями своих фазовых переменных.
Замечание 6..2 В выражение найденного первого интеграла формально необходимо вместо Ci подставить левую часть первого интеграла (6.22).
Тогда полученный дополнительный первый интеграл имеет следующий структурный вид:
©2(^3,^2; Z) = g(sinZ, W:, = C2 = const. (6.40)
V SinZ sinZJ
Итак, найдены два первых интеграла (6.23), (6.40) независимой системы третьего порядка (6.15). Осталось указать один первый интеграл — для системы (6.16) и дополнительный первый интеграл, "привязывающий" уравнение (6.17).
Действительно, искомые первые интегралы имеют следующий вид:
©3(wi; ni) = Л//1 + W2 = C3 = const, (6.41)
sin ni
cos ni
e4(wi; ni, n2) = П2 ± arctg^^ -= C4 = const, (6.42)
C32 sin2 ni - 1
при этом в левую часть равенства (6.42) вместо C3 необходимо подставить интеграл (6.41).
Теорема 6..2 Система (6.15)-(6.17) шестого порядка обладает достаточным количеством (четырьмя)
независимых первых интегралов (6.23), (6.40), (6.4I), (6.42).
Итак, в рассматриваемом случае система динамических уравнений (6.15)—(6.17) имеет четыре первых интеграла, выражающихся соотношениями (6.23), (6.40), (6.41), (6.42) (при этом используются выражения (6.36)—(6.39)), являющихся трансцендентными функциями фазовых переменных (в смысле комплексного анализа) и выражающихся через конечную комбинацию элементарных функций.
Теорема 6..3 Три группы соотношений (2.3), (3.3), (4.3) при условиях (2.4)-(2.6), (6.1), (6.4) обладают четырьмя первыми интегралами (полным набором), являющимися трансцендентными функциями с точки зрения комплексного анализа, выражающимися через конечную комбинацию элементарных функций.
6.4. Топологические аналогии
Предъявим далее еще две группы аналогий, связанных с системой (5.1), описывающей движение свободного твердого тела при наличии следящей силы.
Первая группа аналогий снова касается случая наличия в системе неинтегрируемой связи (5.5). В данном случае динамическая часть уравнений движения при некоторых условиях приводится к системе (5.14).
При выполнении условий (6.1), (6.4) система (5.14) примет вид
а = -(1 + bHi)Z3 + bsin а, (6.43)
cos a
Z'3 = sin a cos a - (1 + bHi)(Z21 + Z22)---H2Z3 cos a, (6.44)
sin a
Z2 = (1 + bHi)Z2Z3 — + (1 + bHi)Z¡ — ^ - HiZ2 cos a, (6.45) sin a sin a sin pi
Zi =(1 + bHi)ZiZ3 — - (1 + bHi)ZiZ2 — ^^ - HZ cos a, (6.46)
sin a sin a sin в2
cos a
el = (1 + bHi)Z2 --, (6.47)
1 sin a
cos a
в2 = -(1 + bHi)Zi-.-—, (6.48)
sin a sin pi
если ввести безразмерные параметры, переменные и дифференцирование по аналогии с (6.5):
2 AB h1 B
b = ano, = ——, Hi = —-, zk = novZk, k = 1, 2, 3, < • >= nov <'> . (6.49)
2I2 2l2'no
Теорема 6..4 Система (6.43)-(6.48) (для свободного тела) эквивалентна системе (6.7)-(6.12) (для закрепленного маятника).
Действительно, достаточно положить
e = a, ni = Pi, П2 = в2, b* = -b, Hi* = -Hi, (6.50)
а также сопоставить переменные Zk ^ —Zk, k =1,2, 3.
Для полного интегрирования системы (6.43)-(6.48) необходимо, вообще говоря, знать пять независимых первых интегралов. Однако после следующей замены переменных
W3 = Z3, W2 = Jz%+zi, wi = ^Г, (6.51)
Z1
система (6.43)-(6.48) распадается следующим образом:
a = -(1 + bH1)w3 + b sin a, w'3 = sin a cos a — (1 + bH^w^^n^ — Hiw3 cos a,(6.52) w'2 = (1 + bHi)w2w3C°na — Hiw2 cos a,
1 i cos в i
w[ = d i(w',w2,w i; a, в i, в2)-—tt , в'\ = d i(w',W2 ,w i; a,ei,e2), (6.53)
w i sin р2
в'2 = d2(w',w2,w 1; a, в 1,в2), (6.54)
d 1 (W3 ,w,w 1; a,ei,e2) = (1 + bHi )Z2 (w',w,w 1) fff = sna,
d2(w3,w2,wi; a, вi, в2) = (6.55)
= -(1 + bHi )Zi (w',w2,w 1) . ^ я = T-r^i ■ c°sa я ,
^ 1 > 1 ^ 3 2' l^sm a sin pi Wl+wf sln a sln Pi
при этом Zk = Zk (W3 ,W2,wi), k = 1, 2, 3, — функции в силу замены (6.51). Система (6.52)-(6.54) рассматривается на касательном расслоении
T*S3{(w3,w2,wl; a, в2, в2) e R6 : 0 < a, в2 < п, в2 mod 2п} (6.56)
трехмерной сферы S3{(a,вl,в2) e R3 : 0 ^ a, вт ^ п, в2 mod 2п}.
Видно, что в системе шестого порядка (6.52)-(6.54) выделяется независимая подсистема третьего порядка (6.52), которая может быть самостоятельно рассмотрена на своем трехмерном многообразии, независимая система второго порядка (6.53) (после замены независимой переменной), а также (по причине цикличности переменной в2) уравнение (6.54) на в2 отделяется.
Таким образом, для полной интегрируемости системы (6.52)-(6.54) достаточно указать два независимых первых интеграла системы (6.52), один — для системы (6.53) и дополнительный первый интеграл, "привязывающий" уравнение (6.54) (т.е. всего четыре).
Следствие 6..1 1. Угол атаки a и углы вт, в2 для свободного тела эквивалентны соответственно углам отклонения e и n1, n2 закрепленного маятника.
2. Расстояние a = CD для свободного тела соответствует длине державки l = OD закрепленного ма-
ятника.
3. Первые интегралы системы (6.52)-(6.54) могут быть автоматически получены через равенства (6.23), (6.40), (6.41), (6.42) после подстановок (6.50):
©1(ws,w2; а) = (1 + bH1)(wg + w2) — (b + H1)ws sin а + sin2 а
w2 sin а
C1 = const. (6.57)
©2(ws, w2; а) = G (sin а, ——, —— ) = C2 = const. (6.58)
2 sin а sin а
©s(wb £1) = Л//1 + w = Cs = const, (6.59)
sin e1
cos
©4(w1; в1,в2) = в2 ± arctg—. = C4 = const, (6.60)
yC sin2 в1 — 1
при этом в левую часть равенства (6.60) вместо Cs необходимо подставить интеграл (6.59).
Вторая группа аналогий касается случая движения с постоянной скоростью центра масс тела, т.е. когда выполнено свойство (5.17). В данном случае динамическая часть уравнений движения при некоторых условиях приводится к системе (5.20)-(5.26). Тогда, в силу условий (5.17), (6.1), (6.4), (6.49) преобразованная динамическая часть уравнений движения (система (5.21)-(5.26)) примет вид аналитической системы
Zs = sin а cos а — (1 + bH1)(Z? + Z|)C0Sа + bZs(Z? + Z2 + Z|)cos а — bZs sin2 а cos а+ s 1 2 sin а 1 2 s
= —Zs + b(Z2 + Z22 + Z|) sin а + b sin а cos2 а — bH1Zs cos2 а, (6.61)
2 2)
+bH1Z| sin а cos а — H1 Zs cos а, (6.62)
V (Л , hU ^ 7 7 cos а , n , ъи w2 cos а cos e1 ,
Z2 = (1 + bH1)Z2Zs---+ (1 + bH1 )ZX---——+
sin а sin а sin p1
+bZ2 (Z2 + Z22 + Z|) cos а — bZ2 sin2 а cos а + bH1Z2Zs sin а cos а — H1Z2 cos а, (6.63)
Z1 = (1 + bH1)Z1Zs ^ — (1 + bH1 )Z1Z2 ^ ^+ sin а sin а sin p1
+bZ1(Z1 + Z22 + Z|) cos а — bZ1 sin2 а cos а + bH1Z1Zs sin а cos а — H1Z1 cos а, (6.64)
cos а
в1 = (1 + bH1)Z2 c-, (6.65)
1 sin а
cos а
в2 = —(1 + bHx)Zx~.-—, (6.66)
sin а sin p1
при этом выбирая постоянную «4 следующим образом: «1 = no.
Для полного интегрирования системы (6.61)-(6.66) необходимо, вообще говоря, знать пять независимых первых интегралов. Однако после замены переменных (6.51) система (6.61)-(6.66) распадается следующим образом:
а' = —ws + b(w2 + w|) sin а + b sin а cos2 а — bH1ws cos2 а, w's = sin а cos а - (1 + bHi )w2cosa+
s V 1 / 2 sln a
+bws(w2 + w|)cos а — bws sin а cos а + bH1w| sin а cos а — H1ws cos а, (6.67)
w2 = (1 + bH1)w2wscosa + bw2(w2 + w|) cos а — bw2 sin2 а cos а+ +bH1w2ws sin а cos а — Hw2 cos а,
w1 = d1(ws,w2,w1; а,^1,в2)—1 c°sel, в = d1(ws,w2,w1; а, ^1,^2), (6.68)
w1 sin p1
в2 = d2(ws,w2,w1; а,в1 (6.69)
где выполнены условия (6.55).
Система (6.67)-(6.69) рассматривается на касательном расслоении (6.56) трехмерной сферы S3{(а,, в1, в2) е Rs : 0 < а, в1 < п, в2 mod 2п}.
Видно, что в системе шестого порядка (6.67)-(6.69) выделяется независимая подсистема третьего порядка (6.67), которая может быть самостоятельно рассмотрена на своем трехмерном многообразии, независимая система второго порядка (6.68) (после замены независимой переменной), а также (по причине цикличности переменной ^2) уравнение (6.69) на в отделяется.
Таким образом, для полной интегрируемости системы (6.67)—(6.69) достаточно указать два независимых первых интеграла системы (6.67), один — для системы (6.68) и дополнительный первый интеграл, "привязывающий" уравнение (6.69) (т.е. всего четыре).
Если вопрос о первых интегралах системы (6.43)—(6.48) (или (6.52)—(6.54)) решается с помощью следствия 6..1, то аналогичный вопрос для системы (6.61)-(6.66) (или (6.67)-(6.69)) решает следующая теорема 6..5.
Сначала отметим, что один из первых интегралов системы (6.67) имеет следующий вид:
e'/(w3,w2; а) _
_ (1 + bH1)(w'2 + w%) - (b + Hi)w3 sin а + sin2 а _ ^ = cQngt (6 7g)
W2 sin а
Далее, изучим вопрос дополнительного первого интеграла системы третьего порядка (6.67), используя при этом первый интеграл (6.70). Для этого введем следующие обозначения и новые переменные:
т _ sin а, W3 _ U2т, W2 _ uiT, p _ -1. (6.71)
т 2
Тогда вопрос о явном виде искомого первого интеграла сводится к решению линейного неоднородного уравнения:
dp _ 2((1 + bHi)u2 - b)p + 2b(1 - H1U2 - u2 - U¡(Ci, U2)) (6 72)
du2 1 - (b + Hi)u2 + (1 + bHi)u2 - (1 + bHi)U2(Ci,U2) ' (. )
Ui(Ci,u2) _ 1 {Ci ± ^C2 - 4(1 + bHi)(1 - (b + Hi)u2 + (1 + bHi)u2)},
при этом постоянная интегрирования C1 выбирается из условия
(b - Hi)2 + C¡ - 4 > 0. (6.73)
Последний факт означает, что может быть найден еще один трансцендентный первый интеграл в явном виде. При этом общее решение уравнения (6.72) зависит от произвольной постоянной C2. Полные выкладки в данном месте приводить не будем, отметив лишь для примера, что общее решение линейного однородного уравнения, полученного из (6.72), даже в частном случае
1 A4 i
Ib - Hi| _ 2, Ci _TTa4, A i _2(b + Hi)
имеет следующее решение:
VC - 4A1(1 - A1 U2)2 ± C\
p = P0(U2)= C[1 - A i Ui]2/( 1+Al)
VC1 - 4A2(1 - A 1 U2)2 T Ci
±Al/( 1+A\)
2(Ai - b) ,
x exp —-д . A-—, C = const.
F (1 + A4)Ai(Aiu2 - 1)'
Тогда искомый дополнительный первый интеграл имеет следующий структурный вид:
/ W3 W2 \
©2'(w3, w2; а) = G (sin a,-,-) = C2 = const, (6.74)
V sin a sin a J
используя при этом обозначения и замены (6.71).
Итак, найдены два первых интеграла (6.70), (6.74) независимой системы третьего порядка (6.67). Осталось
указать один первый интеграл — для системы (6.68), и дополнительный первый интеграл, "привязывающий"
уравнение (6.69).
Действительно, искомые первые интегралы совпадают с первыми интегралами (6.59), (6.60), а именно:
es(wi; ¡3i) = ^ + W2 = C3 = const, (6.75)
Sin pi
cos в i
©4(wi; вьвг) = в2 ± arctg— = C4 = const, (6.76)
yC sin2 ei - 1
при этом в левую часть равенства (6.76) вместо C3 необходимо подставить интеграл (6.75). Теорема 6..5 Четыре первых интеграла (6.70), (6.74), (6-75), (6.76) системы (6.67)-(6.69) являются трансцендентными функциями своих фазовых переменных и выражаются через конечную комбинацию элементарных функций.
Теорема 6..6 Четыре первых интеграла (6.70), (6.74), (6.75), (6.76) системы (6.67)-(6.69) эквивалентны четырем первым интегралам (6.57), (6.58), (6.59), (6.60) системы (6.52)-(6.54).
Действительно, пары первых интегралов (6.70), (6.57), (6.75), (6.59) и (6.76), (6.60) совпадают. Осталось формально отождествить фазовые переменные w^, к = 1,2,3, для системы (6.67)—(6.69) с фазовыми переменными Wk, к = 1,2,3, для системы (6.52)—(6.54). Аналогичные рассуждения, касающиеся пары первых интегралов (6.74), (6.58), не приводим ввиду громоздкости изложения.
х
Заключение
Итак, мы имеем следующие топологичекие и механические аналогии. (1) Движение закрепленного на (обобщенном) сферическом шарнире физического маятника в потоке набегающей среды (неконсервативное поле сил при учете дополнительной зависимости момента сил от угловой скорости). (2) Движение четырехмерного свободного твердого тела в неконсервативном поле сил со следящей силой (при наличии неинте-грируемой связи и при учете дополнительной зависимости момента сил от тензора угловой скорости). (3) Сложное движение четырехмерного твердого тела, вращающегося вокруг центра масс, движущегося прямолинейно и равномерно, и находящегося в неконсервативном поле сил при учете дополнительной зависимости момента сил от тензора угловой скорости.
О более общих топологических аналогиях см. также [9; 10].
Литература
[1] Шамолин М.В. Случаи интегрируемости, соответствующие движению маятника на плоскости // Вестник СамГУ. Естественнонаучная серия. 2015. № 10(132). С. 91-113.
[2] Шамолин М.В. Случаи интегрируемости, соответствующие движению маятника в трехмерном пространстве // Вестник СамГУ. Естественнонаучная серия. 2016. № 3-4. С. 75-97.
[3] Shamolin M.V. New integrable cases and families of portraits in the plane and spatial dynamics of a rigid body interacting with a medium // Journal of Mathematical Sciences. 2003. Vol. 114. №1. P. 919-975.
[4] Шамолин М.В. Многообразие случаев интегрируемости в динамике маломерного и многомерного твердого тела в неконсервативном поле // Итоги науки и техники. Сер.: Современная математика и ее приложения. Тематические обзоры. T. 125. Динамические системы. 2013. C. 5-254.
[5] Походня Н.В., Шамолин М.В. Некоторые условия интегрируемости динамических систем в трансцендентных функциях // Вестник СамГУ. Естественнонаучная серия. 2013. № 9/1(110). С. 35-41.
[6] Шамолин М.В. Многообразие типов фазовых портретов в динамике твердого тела, взаимодействующего с сопротивляющейся средой // Доклады РАН. 1996. Т. 349. № 2. С. 193-197.
[7] Шамолин М.В. Динамические системы с переменной диссипацией: подходы, методы, приложения // Фунд. и прикл. мат. 2008. Т. 14. Вып. 3. С. 3-237.
[8] Арнольд В.И., Козлов В.В., Нейштадт А.И. Математические аспекты классической и небесной механики. М.: ВИНИТИ, 1985. 304 с.
[9] Трофимов В.В. Симплектические структуры на группах автоморфизмов симметрических пространств // Вестн. Моск. ун-та. Сер. 1. Математика. Механика. 1984. №6. C. 31-33.
[10] Трофимов В.В., Шамолин М.В. Геометрические и динамические инварианты интегрируемых гамильтоновых и диссипативных систем // Фунд. и прикл. мат. 2010. Т. 16. Вып. 4. С. 3-229.
[11] Шамолин М.В. Методы анализа динамических систем с переменной диссипацией в динамике твердого тела. М.: Изд-во "Экзамен", 2007. 352 с.
[12] Shamolin M.V. Classes of variable dissipation systems with nonzero mean in the dynamics of a rigid body // Journal of Mathematical Sciences, 2004, 122, no. 1, P. 2841-2915.
[13] Шамолин М.В. Некоторые модельные задачи динамики твердого тела при взаимодействии его со средой // Прикл. механика. 2007. Т. 43. №10. С. 49-67.
References
[1] Shamolin M.V. Sluchai integriruemosti, sootvetstvuiushchie dvizheniiu maiatnika na ploskosti [Cases of integrability corresponding to the pendulum motion on the plane]. Vestnik SamGU. Estestvennonauchnaia seriia [Vestnik of Samara State University. Natural Sciences Series], 2015, no. 10(132), pp. 91-113 [in Russian].
[2] Shamolin M.V. Sluchai integriruemosti, sootvetstvuiushchie dvizheniiu maiatnika v trekhmernom prostranstve [Cases of integrability corresponding to the pendulum motion on the three-dimensional space]. Vestnik SamGU. Estestvennonauchnaia seriia [Vestnik of Samara State University. Natural Sciences Series], 2016, no. 3-4, pp. 75-97 [in Russian].
[3] Shamolin M.V. New integrable cases and families of portraits in the plane and spatial dynamics of a rigid body interacting with a medium. Journal of Mathematical Sciences, 2003, Vol. 114, no. 1, pp. 919-975 [in English].
[4] Shamolin M.V. Mnogoobrazie sluchaev integriruemosti v dinamike malomernogo i mnogomernogo tverdogo tela v nekonservativnom pole [Variety of cases of integrability in dynamics of lower-, and multi-dimensional body in nonconservative field]. Itogi nauki i tekhniki. Ser. "Sovremennaia matematika i ee prilozheniia. Tematicheskie obzory" [Results of science and technics. Series: "Contemporary Mathematics and its Applications. Subject reviews"], Vol. 125, "Dynamical Systems", 2013, pp. 5-254 [in Russian].
[5] Pokhodnya N.V., Shamolin M.V. Nekotorye usloviia integriruemosti dinamicheskikh sistem v transtsendentnykh funktsiiakh [Some conditions of integrability of dynamics system in transcedent functions]. Vestnik SamGU. Estestvennonauchnaia seriia [Vestnik of Samara State University. Natural Sciences Series], 2012, no. 9/1(110), pp. 35-41 [in Russian].
[6] Shamolin M.V. Mnogoobrazie tipov fazovykh portretov v dinamike tverdogo tela, vzaimodeistvuiushchego s soprotivliaiushcheisia sredoi [Variety of types of phase portraits in the dynamics of a rigid body interacting with a resisting medium]. Doklady RAN [Doklady Mathematics], 1996, Vol. 349, no. 2, pp. 193-197 [in Russian].
[7] Shamolin M.V. Dinamicheskie sistemy s peremennoi dissipatsiei: podkhody, metody, prilozheniia [Dynamical Systems With Variable Dissipation: Approaches, Methods, and Applications] Fund. i prikl. mat. [Fundamental and Applied Mathematics], 2008, Vol. 14, issue 3, pp. 3-237 [in Russian].
[8] Arnold V.I., Kozlov V.V., Neyshtadt A.I. Matematicheskie aspekty klassicheskoi i nebesnoi mekhaniki [Mathematical aspects in classical and celestial mechanics]. M.: VINITI, 1985, 304 p. [in Russian].
[9] Trofimov V.V. Simplekticheskie struktury na gruppakh avtomorfizmov simmetricheskikh prostranstv [Symplectic structures on symmetruc spaces of automorphysm groups]. Vestnik Moskovskogo Universiteta. Ser. 1. Matematika. Mekhanika [Moscow University Mathematics Bulletin], 1984, no. 6, pp. 31-33 [in Russian].
[10] Trofimov V.V., Shamolin M.V. Geometricheskie i dinamicheskie invarianty integriruemykh gamil'tonovykh i dissipativnykh sistem [Geometrical and dynamical invariants of integrable Hamiltonian and dissipative systems]. Fund. i prikl. mat. [Fundamental and Applied Mathematics], 2010, Vol. 16, no. 4, pp. 3-229 [in Russian].
[11] Shamolin M.V. Metody analiza dinamicheskikh sistem s peremennoi dissipatsiei v dinamike tverdogo tela [Methods of analysis of various dissipation dynamical systems in dynamics of a rigid body]. M.: Izd-vo "Ekzamen", 2007, 352 p. [in Russian].
[12] Shamolin M.V. Classes of variable dissipation systems with nonzero mean in the dynamics of a rigid body. Journal of Mathematical Sciences, 2004, Vol. 122, no. 1, pp. 2841-2915 [in English].
[13] Shamolin M.V. Nekotorye model'nye zadachi dinamiki tverdogo tela pri vzaimodeistvii ego so sredoi [Some model problems of dynamics for a rigid body interacting with a medium]. Prikl. mekhanika [International Applied Mechanics], 2007, Vol. 43, no. 10, pp. 49-67 [in Russian].
M.V. Shamolin3
CASES OF INTEGRABILITY CORRESPONDING TO THE PENDULUM MOTION IN FOUR-DIMENSIONAL SPACE
In this article, we systemize some results on the study of the equations of motion of dynamically symmetric fixed four-dimensional rigid bodies-pendulums located in a nonconservative force fields. The form of these equations is taken from the dynamics of real fixed rigid bodies placed in a homogeneous flow of a medium. In parallel, we study the problem of the motion of a free four-dimensional rigid body also located in a similar force fields. Herewith, this free rigid body is influenced by a nonconservative tracing force; under action of this force, either the magnitude of the velocity of some characteristic point of the body remains constant, which means that the system possesses a nonintegrable servo constraint. We also show the nontrivial topological and mechanical analogies.
Key words: four-dimensional rigid body, non-conservative force field, dynamical system, case of integra-bility.
Статья поступила в редакцию 11/77/2017. The article received 11/77/2017.
3Shamolin Maxim Vladimirovich (shamolin@rambler.ru, shamolin@imec.msu.ru), Institute of Mechanics, Lomonosov Moscow State University, 1, Leninskie Gory, Moscow, 119192, Russian Federation.