CC BY
21
УДК 517.925+531.01
НЕКОТОРЫЕ УСЛОВИЯ ИНТЕГРИРУЕМОСТИ ДИНАМИЧЕСКИХ СИСТЕМ В ТРАНСЦЕНДЕНТНЫХ
ФУНКЦИЯХ
© 2013 Н.В. Походня,1 М.В. Шамолин2
Изучаются некоторые общие условия интегрируемости в элементарных функциях для систем на касательных расслоениях двумерной сферы. При этом приводится интересный пример трехмерного фазового портрета системы маятникового типа, которая описывает движение сферического маятника, помещенного в поток набегающей среды. Приводятся достаточные условия существования первых интегралов, выражающихся через конечную комбинацию элементарных функций, для многопараметрических систем третьего порядка.
Ключевые слова: динамическая система с переменной диссипацией, интегрируемость, трансцендентный первый интеграл.
1. Предварительные сведения
Для начала в качестве важного примера рассмотрим следующую динамическую систему:
в + Ъв cos в + sin в cos в - ф2^ =0,
1 1 г cos в '
1+cos2 в sin в cos в
0 (1Л)
ф + Ьф cos в + вф
на касательном расслоении T*S2 двумерной сферы S2{в, ф}. Данная система описывает сферический маятник, помещенный в поток набегающей среды (см. также [1-3]. При этом в системе присутствует консервативный момент
sin в cos в, (1.2)
а также момент силы, линейным образом зависящий от скорости с переменным коэффициентом:
Ь(ф )cosв (О)
Оставшиеся коэффициенты в уравнениях являются коэффициентами связности, а именно:
гв = _sintf = 1 +cos2 в (1 4)
М cos в1 ^ sin в cos в. ( . )
хПоходня Наталья Витальевна ([email protected]), кафедра математики и физики Московского государственного гуманитарного университета им. М.А. Шолохова, 109240, Российская Федерация, г. Москва, ул. Верхняя Радищевская, 16—18.
2Шамолин Максим Владимирович ([email protected]), Институт механики Московского государственного университета им. М.В. Ломоносова, 119192, Российская Федерация, г. Москва, Мичуринский пр., 1.
Рис. Относительно грубый фазовый портрет в трехмерной области
Система (1.1) фактически имеет порядок 3, поскольку переменная ф является циклической, при этом в систему входит лишь производная ф. Предложение 1. Уравнение
ф = 0 (1.5)
задает семейство интегральных плоскостей для системы (1.1).
Более того, уравнение (1.5) редуцирует систему (1.1) к уравнению, описывающему цилиндрический маятник, находящийся в потоке набегающей среды (см. также [4]).
Предложение 2. Система (1.1) эквивалентна следующей системе:
в = —¿2 + Ь эт в, ¿2 = эшвсоэв — ,
2 1 бШ 0 '
~ _ ~ ~ сса 0
¿1 = ¿^ЦПв, ф = г1Щ
Т 1 Я1Л 0
(1.6)
на касательном расслоении Т^Б2{¿1, ¿2,в,ф} двумерной сферы Б2{в,ф].
Отделение четвертого уравнения системы (1.6) также произошло по причине цикличности переменной ф.
О построении фазового портрета системы (1.1), изображенного на рисунке, см.
[5; 6].
Пример 1. Исследуем систему, возникающую в пространственной (3Д—) динамике твердого тела, взаимодействующего со средой (см. также [7; 8]):
а = —¿2 + Ь(^2 + ¿|) эт а + Ь эт а соэ2 а, ¿2 = эт а соэ а + Ьх2 (х^ + г"2) соэ а — Ьх2 эт2 а соэ а — х2 ¿1 = Ьх1(х2 + ¿1) соэ а — Ьх1 эт2 а соэ а + ¿1х2с°-а,
(1.7)
которая соответствует следующей системе с алгебраической правой частью:
йг2 = т+Ъ7-2(7-1 + г1)-Ъг2Г 2-г^/т йт = -22 + Ъг(г2 + г2)+Ът(1-т2) , йг1 = Ъг 1 (г 1+г'^)-Ъг 1 т 2+2 1X2/т йт = г2 + Ът(г2 + г2)+Ът(1-т2) '
(1.8)
(1.11)
Итак, мы по-прежнему рассматриваем пару систем: первоначальную систему (1.7) и соответствующую ей алгебраическую систему (1.8).
Аналогичным образом производится переход к однородным координатам uk, k = 1, 2 по формулам
Zk = Uk Т. (1.9)
Система (1.8) приводится с помощью замены (1.9) к виду
Т du2 I u = Т +bu2T3(uj+ul)-bu2T3-и\т
Т dT + U2 = -u2T + bT3(u2 +u2) + bT(1-T2) ' (1 10)
Т dui , u = buiT3(u2+u2)-buiT3 +uiu2T V ' '
Т dT + U1 = -u2T + bT3(u2 +u2) + bT(1-T2) ,
который, в свою очередь, соответствует уравнению
du2 1 — bu2 + u2 — u1 dui 2uiu2 — bui
Данное уравнение интегрируется в элементарных функциях, поскольку интегрируется тождество
, (1 — bu2 + u2\ ,
di-u ) + du1 =0, (1.12)
и имеет в координатах (т, z1,z2) первый интеграл следующего вида (ср. с [9; 10]):
z2 + z| — ez2T + Т2
-= const.
z1 Т
2. Некоторые обобщения
Зададим вопрос: каковы возможности интегрирования в элементарных функциях следующей системы более общего вида, включающей в себя рассмотренные выше системы, в трехмерных фазовых областях:
dz _ ax+by+cz+ciz2/x+C2Zy/x+C3y2/х
dx = dix+ey+f z , (2 1)
dy = gx+hy+iz+iiz2/x+i2zy/x+i3y2/x V ' '
dx dix+ey+fz '
имеющей особенность типа 1/x?
Ранее уже был получен ряд результатов по данному вопросу (см. также [9]). Приведем краткое резюме по данным результатам, а также дополним материал оригинальными рассуждениями.
Вводя, как и ранее, подстановки
y = ux, z = vx, (2.2)
получим, что система (2.1) приводится к следующей системе:
dv ax + bux + cvx + C1v2x + C2vux + c^u? x
xj~ + v =-1—T-"Г7-, (2.3)
dx d1x + eux + fvx
du gx + hux + ivx + ¿1v2x + i2vux + ¿3 u2x
x— + u = -----—-, (2.4)
dx d1x + eux + fvx
эквивалентной
(2.5)
dv ax + bux + (c — dAvx + (c1 — f)v x + (c2 — e)vux + C3u x
x— = -,
dx d1x + eux + fvx
du gx + (h — d1)ux + ivx + ¿1v2x + (¿2 — f )vux + (¿3 — e)u2x
x^- =-;-;-, (2.6)
dx d1x + eux + fvx
которой сопоставим следующее неавтономное уравнение с алгебраической правой частью:
dv a + bu + cv + civ2 + c2vu + C3U2 — v[di + eu + fv] du g + ftu + iv + ¿1 v2 + «2vu + ¿3U2 — u[di + eu + fv]
Интегрирование последнего уравнения сводится к интегрированию уравнения в полных дифференциалах:
[g + hu + iv + ¿iv2 + i2vu + ¿3u2 — diu — eu2 — fuv]dv =
= [a + bu + cv + civ2 + C2vu + C3u2 — div — euv — fv2 ]du. (2-8)
Имеем, вообще говоря, 15-параметрическое семейство уравнений вида (2.8). Для интегрирования в элементарных функциях последнего тождества как однородного уравнения достаточно наложить 6 соотношений:
g = 0, i = 0, ii =0, e = c2, h = c, i2 = 2ci — f. (2.9)
Введем 9 параметров ßi,...,ßg и рассмотрим их в качестве независимых:
ßi = a, ß2 = b, ß3 = c, ß4 = ci, ß5 = c2, (2 io)
ße = c3, ßr = di, ßg = f, ßg = ¿3. (-)
Таким образом, уравнение (2.8) при выполнении групп условий (2.9), (2-10) сводится к виду
dv ßi + ß2u + (ß3 — ßr)v + (ß4 — ßg)v2 + ßeu2
du (ß3 - ßr)u + 2(ß4 - ßs)«u +(ßg - ß5)u2 а система (2.5), (2.6), соответственно, к виду
(2.11)
= в1 + в2и + (вз — вг)« + (в4 — вв>2 + веи2 (2 12)
¿Ж в7 + вб^ + вв^ '
= (вз — вг)ц + 2(в4 — вв)уц + (во — вб )и2 (2 13) ¿X в7 + вби + вв^ '
после чего уравнение (2.11) интегрируется через конечную комбинацию элементарных функций.
Действительно, интегрируя тождество (2.8), получаем следующее соотношение:
d
(ß3 - ßr)«'
+ d
(ß4 - ßs)«
+ d[(ßg - ß5)«] + d
ßi u
u
-dtf2 ln |u|] - d[e6u]=0, (2.14)
которое для начала позволяет получить инвариантное соотношение:
(вз - вг)у , в - вв)У2 + (в в ) + в1
--1---+ (в9 - в5 )v +---
u u u
-в2 ln |u| - в6и = C1 = const, (2.15)
а затем в координатах (x,y,z) — первый интеграл в виде:
(в4 - es)z2 - вбУ2 + (вз - e?)zx + (в9 - в5)zy + в1Х2
-ß2 ln
yx
У
x
const. (2.16)
Таким образом, можно сделать вывод об интегрируемости в элементарных функциях следующей, вообще говоря, неконсервативной системы третьего порядка, зависящей от 9 параметров:
dz = ¡3ix+@2y+@3Z+@4Z2/x+@5 zy/x+¡36y2/x dx p7x+p5y+psz
dy = e3y+(2e4-es)zy/x+e9y2/x dx fÍ7x + fÍ5y+fÍ8Z '
(2.17)
Следствие. Система третьего порядка на множестве
81{атоа2п}\{а = 0, а = п} х В2{г1 ,г2}, (2.18)
зависящая от 9 параметров
а = вг sin а + + e»Z2, z2 = в1 sin a cos а + в2 z1 cos а + в3г2 cos а+ +в4 + esz^S^ + в^2^,
2 sin а 1 2 sin а 1 sin а '
Z1 = e3Z1 cos а + (2в4 - es)z1 z?^ + вог?^
(2.19)
обладает, вообще говоря, трансцендентным первым интегралом, выражающимся через элементарные функции:
(в4 - вв)г? - вбг? + (вз - вг)г2 sin а + (во - вб^! + в1 sin2 а2
-в2 ln
z1 sin а Z1
const. (2.20)
sin а
В частности, система (2.19) при в1 = 1, в2 = вз = в4 = в5 = вЭ = 0, вб = = вв = — 1, в7 = b имеет вид системы, составленной из первых трех уравнений формулы (1.6).
Для нахождения дополнительного первого интеграла неавтономной системы (2.1) используется найденный первый интеграл (2.16), выражающийся через конечную комбинацию элементарных функций.
Преобразуем для начала соотношение (2.15) следующим образом:
(в4 — вв>2 + [(вэ — в5)и + (вз — вг)] v + fi(u) = 0, (2.21)
где
fi(u) = в1 — вбu2 — в2и ln |u| — Ciu. При этом формально величину v можно найти из равенства
1
2(в4~ вв)
где 2
f2(u) = A1 + A2u + A3u2 + A4uln luí, Ai = (вз — вг)2 — 4в1(в4 — вв), A2 = 2(вэ — вб)(вз — вг) + 4С1(в4 — вв), Аз = (вэ — в5 )2 + 4вб(в4 — вв), A4 = 4в2(в4 — вв). Тогда искомая квадратура для поиска дополнительного, вообще говоря, трансцендентного первого интеграла (например, системы (2.12), (2.13) или (2.5), (2.6), при этом используется уравнение (2.13)) примет следующий вид:
f dx _ f [вг + в5u + ввV1,2(u)]du _
V1,2(u) = 2(в _ в ){(вб - во)и + (в7 - вз) ± у/ш} , (2.22)
(вз - вг)и + (во - вб)и2 + 2(в4 - вз)м^112(м) = Г [B1 + B2u + B3^f2(n)]du J B4UV/ f2 (u) '
(2.23)
Bk = const, к = 1, . . . , 4.
Искомая же квадратура для поиска дополнительного, вообще говоря, трансцендентного первого интеграла (например, системы (2.12), (2.13) или (2.5), (2.6), когда при этом используется уравнение (2.12)) примет следующий вид:
i dx I'_[в7 + вбЦ(у) + e&v]dv__(2 24)
J x J + ^2«(v) + (вз - £r)v + (в4 - es)v2 + ¡e6u2(v)' (. )
при этом функция u(v) должна быть получена в результате разрешения неявного уравнения (2.15) относительно u (что в общем случае не всегда очевидно).
Достаточные условия выражения интегралов в (2.24) через конечную комбинацию элементарных функций дает следующая Теорема 1. При A4 = 0, т. е. при
в2 =0 (2.25)
или при
в4 = в8 (2.26)
неопределенный интеграл в (2.24) выражается через конечную комбинацию элементарных функций.
Следствие. Система (2.19) при выполнении достаточных условий теоремы 1 (в данном случае выполнено свойство (2.25)), обладает полным набором первых интегралов, выражающихся через конечную комбинацию элементарных функций.
Рассматриваемые в данной работе динамические системы относятся к системам с переменной диссипацией с нулевым средним по имеющейся периодической координате [9]. Более того, такие системы часто обладают полным списком первых интегралов, выражающихся через элементарные функции.
Метод приведения исходных систем уравнений с правыми частями, содержащими полиномы от тригонометрических функций, к системам с полиномиальными правыми частями позволяет искать (или же доказывать их отсутствие) первые интегралы для систем более общего вида, а не только тех, которые обладают указанными симметриями (см. также [10; 11]).
Литература
[1] Шамолин М.В. К задаче о движении тела в среде с сопротивлением // Вестн. МГУ. Сер. 1. Математика. Механика. 1992. №1. С. 52-58.
[2] Шамолин М.В. Многообразие типов фазовых портретов в динамике твердого тела, взаимодействующего с сопротивляющейся средой // Доклады РАН. 1996. Т. 349. №2. С. 193-197.
[3] Шамолин М.В. Введение в задачу о торможении тела в сопротивляющейся среде и новое двухпараметрическое семейство фазовых портретов // Вестн. Моск. ун-та. Сер. 1. Математика. Механика. 1996. №4. С. 57-69.
[4] Шамолин М.В. Новые интегрируемые по Якоби случаи в динамике твердого тела, взаимодействующего со средой // Доклады РАН. 1999. Т. 364. №5. С. 627-629.
[5] Шамолин М.В. Пространственные топографические системы Пуанкаре и системы сравнения // Успехи матем. наук. 1997. Т. 52. Вып. 3. С. 177-178.
[6] Шамолин М.В. Классификация фазовых портретов в задаче о движении тела в сопротивляющейся среде при наличии линейного демпфирующего момента // Прикл. матем. и механ. 1993. Т. 57. Вып. 4. С. 40-49.
[7] 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. V. 114. №1. P. 919-975.
[8] Шамолин М.В. Интегрируемость по Якоби в задаче о движении четырехмерного твердого тела в сопротивляющейся среде // Доклады РАН. 2000. Т. 375. №3. С. 343-346.
[9] Shamolin M.V. Classes of variable dissipation systems with nonzero mean in the dynamics of a rigid body // Journal of Mathematical Sciences. 2004. V. 122. №1. P. 2841-2915.
[10] Шамолин М.В. Об одном интегрируемом случае уравнений динамики на so(4) х R" // Успехи матем. наук. 2005. Т. 60. Вып. 6. С. 233-234.
[11] Походня Н.В., Шамолин М.В. Новый случай интегрируемости в динамике многомерного тела // Вестник СамГУ. Естественнонаучная серия. 2012. № 9(100). С. 136-150.
Поступила в редакцию 18/XI//2013;
в окончательном варианте — 19/XII/2013.
CERTAIN CONDITIONS OF INTEGRABILITY OF DYNAMICAL SYSTEMS IN TRANSCENDENTAL
FUNCTIONS
© 2013 N.V. Pokhodnya? M.V. Shamolin4
Certain general conditions of integrability in elementary functions for the systems on the tangent bundle of two-dimensional sphere are studied. At that an interesting example of three-dimensional phase pattern of pendulum-like system which describes the motion of spherical pendulum, placed in an over-run medium flow. Sufficient conditions of existence of the first integrals expressed through the finite combination of elementary functions, for multi-parametric third order systems are presented.
Key words: variable dissipation dynamic system, integrability, transcendental first
integral.
Paper received 18/XI/2013. Paper accepted 19/ХД/2013.
3Pokhodnya Natalya Vitalievna ([email protected]), the Dept. of Mathematics and Physics, Sholokhov Moscow State University for the Humanities, Moscow, 109240, Russian Federation.
4Shamolin Maxim Vladimirovich ([email protected]), Institute of Mechanics, Moscow State University, Moscow, 119899, Russian Federation.
