О стационарных движениях механической системы с частным интегралом Стеклова
Автор: Новиков М.А.
Журнал: Вестник Бурятского государственного университета. Математика, информатика @vestnik-bsu-maths
Рубрика: Теоретическая механика
Статья в выпуске: 3, 2024 года.
Бесплатный доступ
В статье методом Рауса - Ляпунова найдены стационарные движения механической автономной консервативной системы, для которой возможно существование дополнительного частного интеграла Стеклова. В зависимости от количества интегралов, участвующих по методу Рауса - Ляпунова в связке из первых интегралов, установлен ряд свойств. В частности показано, что при последовательном уменьшении числа интегралов в функции Лагранжа анализ решений стационарности проводится с меньшими вычислительными операциями. При всех способах составления связки интегралов получаются одинаковыми стационарные движения. Последним свойством установлена теорема, согласно которой функция Лагранжа с наименьшим количеством участвующих необходимых интегралов получается при не включении в нее интеграла Стеклова, но при этом выполняются возникающие в случае условий Стеклова частные интегралы: q = const, r = 0. При таком построении связки интегралов решения стационарности почти полностью совпадают со стационарными движениями.
Стационарное движение, частный интеграл, связка интегралов, решение стационарности
Короткий адрес: https://sciup.org/148329911
IDR: 148329911 | DOI: 10.18101/2304-5728-2024-3-19-30
Текст научной статьи О стационарных движениях механической системы с частным интегралом Стеклова
Отдельные динамические свойства механических автономных консервативных систем часто выявляются при исследовании стационарных движений. Хорошим объектом для этого может быть задача о вращении твердого тела вокруг неподвижной точки [1–3]. Движение твердого тела описывается известными дифференциальными уравнениями:
'AP = ( B - C ) qr + Mg ( y o Y 3 - z 0 Y 2 X Y 1 = r Y 2 - q / 3 ,
“ Bq = ( C - A ) rp + Mg ( z o / i - x o Y 3 ) y 2 = P Y 3 - r Y l , (1)
.Cr =(A - B)pq + Mg(xo/2 - Уo/i), Y3 = q/i - PY2, где A, B, C — моменты инерции твердого тела относительно главных осей Ox, Oy, Oz; p,q, r — проекции мгновенной угловой скорости на подвижные, связанные с телом оси; xo, yo, zo — координаты центра масс в подвижных осях; /1, /2, /3 — проекции ортов подвижных осей на неподвижную вертикальную ось OZ, направленную к притягивающему центру (углы Пуассона), M — масса тела; g — ускорение свободного падения.
Для системы (1) известны три общих интеграла [3]:
V , = Ap 2 + Bq 2 + Cr 2 - 2( x o / 1 + y o / 2 + z o / 3) = c o = const (интеграл энергии), V = Ap / 1 + Bq / 2 + Cr / 1 = c 1 = const (интеграл кинетического момента), V 2 = Y 1 2 + Y 2 + Y 3 = 1 (интеграл Пуассона).
Проведем изучение стационарных движений для системы (1), содержащей дополнительный частный интеграл Стеклова. В работах [3, 4] показано существование дополнительных частных интегралов системы (1) при определенных значениях xo = o, zo = o, (2)
а для моментов инерции выполняются дополнительные условия. При условиях (2) предварительно записываются движения тела в форме шести обыкновенных дифференциальных уравнений. Вводя обозначения:
x = ( p , q , r , / 1, / 2, / 3) ' — фазовые переменные; m = Mgy o/ A ; k1 = B / A ;
k 2 = C / A ( m ^ o; o < k 1 ^ 1; k 2 > o), уравнения движения задаются в виде:
p =(k1 -k2) qr + mY3, Y1 = rY2 -q/3, q = (k2 - 1) rp, Y2 = p/3 - Y (3)
. r = - [(1 - k 1 ) pq + m / 1 )]/ k 2 , Y 3 = q / 1 - p y 2 .
Для системы (3) два первых интеграла запишутся в более короткой форме:
V ю = p 2 + k 1 q 2 + k 2 r 2 - 2 m / 2 = const , V o = p / 1 + k 1 q / 2 + k 2 r Y 3 = const .
В системе с частным интегралом В. А. Стеклова дополнительно накладываются условия на моменты инерции.
-
1 Об известных решениях исследуемой задачи
Простейшая классификация частных решений системы (3) может быть проведена при обращении в нуль дифференциального равенства для q = o . В таком случае существует общий интеграл
V3 ( x ) = q = q o = const .
Это может достигаться в трех ситуациях: 1) r = o; 2) p = o; 3) A = C .
Как показано в [4], при условиях r = 0 и B = 2 A ( k 1 = 2 k 2) в системе (3) существует другой дополнительный интеграл В. А. Стеклова
V4(x ) = m Y - 2 q 2 у 2 = const . (1.1)
В этом случае решения для фазовых переменных [3, 4] записываются при помощи эллиптического интеграла. Последний интеграл может быть в системе (3) и при A = C , когда имеются другие решения. Поэтому отделение решений с частным интегралом (1.1) от других возможных решений, существующих без интеграла Стеклова, можно осуществить при условиях
C * A = B /2. (1.2)
Такое отделение необходимо для выделения стационарных движений, содержащих только дополнительный интеграл (1.1).
-
2 Некоторые свойства систем с частным интегралом Стеклова
Существование частного интеграла (1.1) в системе (3) обеспечивается условиями r = 0, k j = 2, k 2 * 1. Любое алгебраическое равенство в фазовых переменных можно рассматривать дополнительным множеством (условием), наравне первому интегралу. Поэтому введем в качестве частного интеграла
-
V 41 ( x ) = r = 0.
Нахождение стационарных движений традиционно проводится методом Рауса — Ляпунова [5-7]. Как в [2] этот метод опирается на известные первые интегралы или дополнительные условия. Его основу составляет отыскание фазовых переменных при поиске экстремума одного из первых интегралов движения системы (3), когда выполняются остальные интегралы. По общей теории экстремумов [8] условный экстремум можно свести к безусловному, вводя функцию Лагранжа, составленную связкой из всех известных первых интегралов и дополнительных условий. Тогда экстремум ищется по всем фазовым переменным, включая в их число вещественные множители Лагранжа, входящие в связку интегралов. Конечно, здесь возникает частично условная задача экстремума, так как общие интегралы, принимающие постоянные, но неопределенные значения, дополнительных числовых условий не доставляют.
Вначале проведем нахождение стационарных движений при полном наборе всех известных первых интегралов и дополнительных соотношений:
-
V )0 ( x ) = const , V 0 ( x ) = const , V , ( x ) = 1, V 3 ( x ) = const ,
V4 ( x ) = const , V 41 ( x ) = r = 0.
Искомая функция Лагранжа имеет вид
K 0 ( x , 2 ) = V ,0( x ) - Л V 0 ( x ) - Л 2 V 2 ( x ) - 2 , V 3 ( x ) - 2 4 V 4 ( x ) - 2 5 V 4X x ) с вещественными множителями 2 i ( i = 1,2, l ,5).
Уравнения стационарности составляются обращением в нуль частных производных от функции K0(x, А) по всем фазовым переменным и множителям Лагранжа. Притом учитываются только соотношения, соответствующие интегралам с фиксированными константами, в данном случае V,(x) = 1, V41( x) = 0. Особо следует выделить особенность, что в выраже- д V (x) д V (x)
ниях —4-- и —4-- уравнений стационарности величина q0 исполь- д/1 ду 2
зуется одинаково с переменной q. Поэтому уместно изучить вопрос о дV(x).
выражении —4 : считать его обычной производной при замене q0 на д q q , или полагать равным нулю. С этой целью изучим оба варианта.
Пусть вначале q 0 рассматривается переменной, равной q . С учетом этого уравнения стационарности запишутся восьмью алгебраическими уравнениями от одиннадцати переменных:
дK0( x, А)/дp = 2 p - AY1 = 0, дK 0( x, А)/д q = 4 q - 2А Y 2 - А + 4А4 qY 2 = 0, дK0 (x, А) / дr = 2k2 r - А k2 Y3 - А5 = 0, дK0(x,А) /дY1 = -(А P + 2А2 y1 + 2тА4 Y1) = 0, дK0(x, А) / дY2 = -2(т + А q + A Y2 - А4 q2) = 0,
(2.1)
дK 0( x, A)/ дY 3 =-( k2 Ar + 2A2 y 3) = 0, дK0(x,А)/дА =-(Yi2 + Y2 + Y32 -1) = 0, дK 0( x, А)/ дА5 = - r = 0.
После подстановки r = 0 в (2.1) получается система упрощенных уравнений
F = 2 P - А Y 1 = 0, F 2 = 4 q - 2 А Y 2 - А + 4 А 4 q Y 2 = 0, F 3 = А k 2 Y 3 + А 5 = 0, - F 4 = А р + 2 Y 1 ( A 2 + т А 4 ) = 0, Fs = m + A q + А Y 2 - А 4 q 2 = 0, F e = А 2 Y 3 = 0, F 7 = Y 12 + Y 22 + Y 32 - 1 = 0.
(2.2)
пре-пер-
Из F 3 = 0 можно определить y 3 = - А 5 / ( k 2 Y 2). Равенство F 6 = 0 доставляет две ветви решений стационарности: 1) А = 0, 2) y 3 = 0 .В вом случае А = 0 уравнения (2.2) становятся проще:
F 0 = 2 P - A Y i = 0, F w = A P + 2 m & 4 Y i = 0,
F 20 = 4 q - 2 А у 2 - А , + 4 А 4 q Y 2 = 0, F 50 = m + A q - А 4 q 2 = 0.
Выразим из F 0 = 0 величину p = Л 1 Y 1/ 2 и подставим в F 40 = 0 :
Y ( Л 2 + 4 m Л 4) = 0.
(2.3)
Одним из возможных решений стационарности отсюда является Y 1 = 0, из него можно заключить p = 0. Предполагая в общем случае Л 4 ^ 0, равенство F 50 = 0 имеет решение q * = ( Л ± 7 Л + 4 m Л 4 ) / (2 Л 4). При Л 4 = 0 решением будет q * = — m / Л 1 .
Равенство F 20 = 0 используем для определения величины Л 3. Легко видеть, что ввиду отличных от нуля произвольных параметров Л и Л 4 величина q * может принимать любые вещественные значения. Тогда одним из решений стационарности (2.2) можно записать:
( P.C . 1) p = 0, q = q * = const , r = 0, y = 0,0 < Y 2| < 7 1 - Y 3 2 , Y 3 = — Л 5 / ( k 2 Л ).
Конечно, не каждое решение стационарности должно быть стационарным движением. Согласно [2, 7] к стационарному движению выдвигается требование постоянства решений (или не зависимость от времени), что означает равные нулю правые части уравнений движения (3). Приведенные к простому виду они здесь будут такими:
fi = (2 - k2) qr + mY 3 = 0, f, = rp = 0, f, = pq + m^ = 0, ft = qY 3 = 0,
f 5 = p ^ = 0, f , = q Y i - p y 2 = 0.
(2.4)
Для выполнения f1 = 0 необходимо y 3 = 0, откуда следует Л 5 = 0. Тогда одним из видов стационарных движений будет
( 1 .1) p 0 = 0, q 0 = const , Г, = 0, Y 10 = 0, Y 20 = ± 1, Y 30 = 0.
Частным движением к нему будет состояние покоя
( 1 .0) p 0 = 0, q 0 = 0, r 0 = 0, Y 10 = 0, Y 20 = ± 1, Y 30 = 0, которое легко проверяется равенствами (2.4).
В другом случае выполнения равенства (2.3), когда y 1 ^ 0, Л 4 =- Л 1 2 / (4 m ), выражение F 5 приводится к к виду ( Л q + 2 m ) 2 = 0. Для нетривиальных значений y 2, q отсюда можно заключить Л 1 =— 2 m / q . Подставляя найденные в F 0 = 0 значения, получим pq + m Y 1 = 0 . Решением уравнений стационарности здесь будет:
( P.C .2) p = — m ^ -, q ^ 0, r = 0, Y 1 = ±7 1 — Y 22 - Y 3 ,0 < | y 21 ^ 1, Y 3 = - Л / ( k 2 Л ).
q
Здесь параметрами следует рассматривать q , y 2, Л 5. Значение Л 2 = 0 по предыдущей ветви решений, и из F 20 = 0 вычисляется Л 3 = 0 . При подстановке найденных решений в (2.4) получается более общий вид стационарных движений:
( I -2) p 0 = q 0 — , q 0 = ±4- m Y 20 , r o = 0, / 10 = ± л/ 1 - Y 220 ,0 < [ Y 20 [ ^ !> Y 30 = 0 • Y 20
При этом должно выполняться sign Y 20 = - signm • Движения ( 1 .1) и ( 1 .0) различны, и совпадают между собой только при q 0 = 0 .
Для другой возможности при y 3 = 0 из F 30 = 0 сразу следует 2 5 = 0 .
Исключение p из F 0 = 0 приводит F 40 = 0 к равенству:
Ш2 + 4а + m^4)] = 0.
Здесь так же одним из решений является y , = 0, которое сводится к стационарному движению вида ( 1 . 1).
При подстановке другого решения Л4 = -(22 + 422) / (4m) в F0 = 0 в предположении y ^ 0 получает равенство
( 2 1 q + 2 m ) 2 + 4 2 2 q 2 = 0.
(2.5)
Квадратное относительно q уравнение имеет вещественные решения
q ** = 2 m
( - Л ± 2 y- 2 J
2 12 + 4 2 2
Решением стационарности здесь будет:
( P.C .3) p = ^р q = q „, r = 0, y , = ±У 1 - Y 22 ,0 < | y 21 ^ 1, Y 3 = 0 .
При подстановке найденных решений в равенства (2.4) найдем Р = q Y , / Y 2. Тогда из равенства f , = 0 следует Y i ( q 2 + m Y 2) = 0.
При y = 0 уже получены стационарные движения вида ( 1 . 1), а при
Y ^ 0, в частности, можно получить 2 2 = 0, и затем стационарные движения вида ( 1 .2).
Стационарные движения вида (1.2) реализуются при следующих значениях множителей Лагранжа: 24 — вещественный параметр, У - mY, т
2 = ± 2 y ^^, 2 2 = — - m 2 4 , 2 3 = 4 q 0 2 4 y 2 0 , 2 = 0 .
Y 20 Y 20
Для другого вида стационарных движений ( 1 . 1) множители Лагранжа такие: 2 1 , 2 4 — свободные параметры, 2 2 = ( 2 4 q 2 - 2 1 q 0 - m ) / y 20;
2 3 = 2 y 20 (2 2 4 q 0 - 2 1 ); 2 5 = 0.
Движения видов ( 1 .1) и ( 1 .2), вообще говоря, разные, и совпадают между собой в случае, когда для ( 1 .2) значение y 10 = 0, а для ( 1 .1) определено значение q 0 = ±У- m Y 20 .
Стационарное движение вида (1.0) может быть осуществлено при множителях Лагранжа: 21,24 — свободные параметры, 22 =-m / Y20; 23 =-2y2021; 25 = 0.
Рассмотрим теперь нахождение решений стационарности и стационар- ных движений при допущении
d V 4( x ) ∂ q
= 0. Единственным отличием от
предыдущего будет замена выражения F 2 = 0 на F 21 = 4 q - 2 A 1 Y 2 - A 3 = 0.
Последнее равенство так же используется для определения A3. Проводя аналогичные вычисления , получим те же самые решения стационарно сти, соответственно — такие же стационарные движения.
Сопоставляя два способа задания
d V 4 ( x ) ∂ q
, можно сформулировать
следующее свойство.
С в о й с т в о 1
При отыскании стационарных движений методом Рауса — Ляпунова механической системы (3) с частным интегралом Стеклова при полной связке интегралов, составленные уравнения стационарности как для независимой переменной q 0 , так и для q 0 , используемой переменной q, одинаково приводят к стационарным движениям видов ( 1 .0) ( 1 .1) ( 1 .2) .
При различных представлениях значений
д V ( x ) .
—4— имеется общее за-∂q ключение для всех видов стационарных движений о равном нулю множителе Лагранжа А5. Фактически это означает отсутствие в связке интегралов V41 (x). Поэтому предварительно исключим множество V41 (x) из функции Лагранжа, и подставим r = 0 в общие интегралы
V )0 ( x ) = const , V 0 ( x ) = const .
При этом в уравнениях (2.2) не будут участвовать слагаемые с А5. Ко- личество переменных при таком подходе уменьшается на два, и на единицу уменьшается число равенств в системе (2.2). В этом случае порядок нахождения решений стационарности остается прежним. Конечно, упрощается отыскание решений стационарности от меньшего числа переменных . Так как решения стационарности составляются из тех же равенств F1 = 0 , F4 = 0 , F5 = 0 , то и стационарные движения остаются теми же (1.0) (1. 1) (1.2). При таком способе задания связки интегралов проводится значительно меньшее количество вычислений. Отсюда следует следующее свойство.
С в о й с т в о 2
При нахождении стационарных движений методом Рауса — Ляпунова механической системы (3) с частным интегралом Стеклова целесообразно предварительно выполнить подстановку r = 0 в исходные первые интегралы.
Можно показать, что и в получаемых решениях стационарности для (1.0) (1. 1) (1.2) получается A4 = 0. Следовательно, в связке интегралов не участвует интеграл Стеклова V4(x) = const, но условия (1.2) и r = 0 предполагаются выполненными.
Выпишем отдельно получаемые здесь уравнения стационарности, и для ясного преставления результатов вместо Л 1 , Л 2 л , введем новые множители Лагранжа д , д 2 д 3:
F 1 = 2 p - д Y 1 = 0, F 21 = 4 q - 2 д 1 Y 2 - д 3 = 0,
‘ F 41 = д 1 p + 2 m д 2 = 0, F 51 = m + д 1 q + д 2 у 2 = 0, (2.6)
F 61 = П 2 Y з = 0, F 7 = Y 2 + Y 22 + Y з2 = 1
Представленная система нелинейных уравнений состоит из шести уравнений, и зависит от восьми переменных.
Из F 1 = 0 выразим p = д Y 1 /2 и подставим в F 41 = 0:
Y ( m д 1 + 4 д 2) = 0.
При y = 0 в предположении д 1 ^ 0 из F51 = 0 найдем q = - ( m + д2 y 2)/ д. Равенство F 21 = 0 предназначено для вычисления Д з = 2( q — Ц 1 Y 2 ).
Решением стационарности в этом случае будет
( P .C .4) p = 0, q = - m / д, r = 0, Y 1 = 0, 0 < | y 2| < 1, Y 3 = ± 1 - — Y 22 •
Здесь д, y 2 рассматриваются параметрами. При подстановке найденных решений стационарности в (2.4) получается: Y 3 = 0, Y 2 =± 1. А ввиду неопределенных значений д е R полагаем q произвольным. Следовательно, в этом случае стационарными движениями получаются движения вида ( 1 . 1) .
При Y 3 = 0 в предположении д 1 ^ 0 решением стационарности будет ( P .C .5) p = Uy , q = - m - д Y 2 /4 , r = 0, Y 1 =± 1 -^ ,0 < Y 2I < 1, Y 3 = 0
2 П 1
Несложной проверкой равенств (2.4) для последнего решения стационарности получаются стационарные движения вида ( 1 .2) .
Таким образом, справедливо следующее свойство.
С в о й с т в о 3
Нахождение стационарных движений механической системы (3) с частным интегралом Стеклова методом Рауса — Ляпунова можно проводить при неполной связке интегралов, исключая интеграл Стеклова, V 4( x ) = const, сохраняя лишь условия его существования (1.2) и r = 0 .
Для такой функции Лагранжа выполняется еще меньше вычислений. В такой форме составленной связки интегралов стационарными движениями получаются только виды (1.0) (1. 1) (1.2). Наиболее общие стационарные движения вида (1.2) соответствуют следующим значениям множителей Лагранжа: д1 = ±^-m / y20; д2 = m; д3 = 0. В остальных стацио- нарных движениях так же р3 = 0. Последнее соответствует не включению общего интеграла V3 = q = const в состав связки интегралов. Тогда в уравнениях (2.6) следует изначально полагать р3 = 0 .Такая же ситуация может возникать и в других системах с аналогичными частными интегралами выше приведенного вида V3(x) = q = const. Поэтому проведем полные выкладки вычислений стационарных движений в этом конкретном случае. Уравнения стационарности здесь будут представлены системой (2.6), в которой будет изменен только полином F21 = 0, имеющий вид F22 = 2q - Pi Y2 = 0.
При р.2 = 0 для системы стационарности имеется такое решение
( P.с .6) p = r = y з = 0, q -'^ m Y , Y i = ±F i - y 22 ,0 <| y 21 ^ 1.
В результате подстановки найденных решений в (2.4) получаются стационарные движения вида p0 = 0, q0 =±V—mY20 / 2, Г0 = 0, Y10 = 0, Y20 =±1, Y30 = 0, являющиеся частным случаем (I.1) .
При y 3 = 0 исключение n i из равенств { F 1 = 0, F 22 = 0 } приводит к
P Y 2 = q Y 1 . Далее из { F 1 = 0, F 41 = 0 } исключение p составляет равенство
Y i ( P i + 4 P 2 ) = 0.
Полагая y 2 ^ 0 исключение п 1 из { F 22 = 0, F 51 = 0 } получает равенство 2 q 2 + m Y 2 + р 2 Y 2 2 = 0. Тогда уравнения стационарности приводятся к более простой системе:
p y 2 = q Y i ,
Y i ( P i2 + 4 P ) = 0,
< 2 q 2 + m Y 2 + P 2 Y 22 = 0, (2.7)
2 q - P i 2 2 = 0,
Y i2 + 2 + Y 32 = 1.
Второе равенство (2.7) при Y i = 0 выделяет часть решений стационарности. Этому соответствуют значения:
p = y 2 = 0, q = ±V- Y 2 ( m + P 2 Y 2 )/2 .
Из множества всех вещественных значений р 2 е ( -^ ; + « ) всегда можно подобрать такие, что для q получаются произвольные вещественные значения. Отсюда возникает решение стационарности ( P.C .7) p = r = Y i = Y 3 = 0, q = const , Y i =± i.
Легко видеть, оно включает в себя решения стационарности (P.C.6), и в точности соответствует стационарным движениям вида (I.1) .
Другое решение второго соотношения (2.7) ^ 2 --4 214 приводит третье равенство (2.7) к виду: 2 q 2 + m y 2 - ^ 2 у 2 /4 - 0. Последнее при подстановке р. 1 - 2 q / у 2 сводится к q 2 - m y 2 - 0. Окончательно получается решение стационарности
Р - q —, q -Ч 'm Y ^ , r - Y з = 0, Y i =± 7 1 - Y 22 ,0 < Y 2I ^ 1, Y 2
в точности совпадающее со стационарным движением вида ( 1 .2).
Как видно из приведенных выкладок, имеет место следующее свойство.
С в о й с т в о 4
Нахождение стационарных движений механической системы (3) с частным интегралом Стеклова методом Рауса — Ляпунова при неполной или не совсем полной связке интегралов необходима постоянства решений стационарности равенствами (2.4).
Вместе с тем справедлива следующая теорема.
Т е о р е м а
В механической системе (3) с частным интегралом Стеклова при нахождении стационарных движений методом Рауса — Ляпунова с неполной связкой интегралов, в которой участвуют только интеграл энергии, кинетического момента и Пуассона, получаемые решения стационарности почти полностью совпадают с стационарными движениями.
Кроме системы (3) имеются другие механические системы, содержащие общие интегралы вида V3 ( x ) - q - q 0 - const , где так же может оказаться p - const , или r — const . Для таких систем могут быть полезны изученные здесь свойства. Конечно, важный интерес представляют системы с четвертым дополнительным интегралом, у которых решения стационарности в точности совпадают с стационарными движениями.
-
3 Обсуждение полученных результатов
С учетом знаков y 20 стационарные движения видов ( 1 .0) и ( 1 . 1) каждое содержит два простых движения. Для стационарных движений вида ( 1 .2) ввиду двух знаков подкоренных выражений для q 0 различаются четыре простейших движения. В совокупности для систем с частным интегралом Стеклова допускается восемь простых стационарных движений.
Выявленные свойства для систем с частным интегралом Стеклова связаны с вычислительными вопросами применения метода Рауса — Ляпунова. В статье последовательно устанавливается тенденция уменьшения вычислительных операций: первоначально — при исключении множества r - 0, затем — при уменьшении числа первых интегралов в функции Лагранжа (сначала исключается частный интеграл Стеклова, далее исключается общий интеграл V3 (x) - q - const). При этом уточняется обязательное участие основных общих интегралов: энергии, кинетического момента, Пуассона. Показано предельно наименьшее число первых инте- гралов в связке, равное трем. Вместе с тем отмечается обязательное выполнение необходимых условий существования частного интеграла Стеклова, выраженных соотношениями (1.2) и r = 0.
Существует предположение, что полный набор стационарных движений выявляется при участии связки из наибольшего числа первых интегралов. Демонстрируемые в статье способы составления функций Лагранжа для системы с частным интегралом Стеклова показательны тем, их участие в методе Рауса — Ляпунова не добавляет стационарных движений, а только усложняет вычислительный процесс. Хотя это не противоречит исходному высказыванию о наибольшем количестве стационарных движений, но выражает другое предположение, что наибольшее число стационарных движений по методу Рауса — Ляпунова находится не обязательно при полной связке всех первых интегралов изучаемой системы.
В исследуемой задаче оказалось необязательным формальное выражение дополнительного интеграла Стеклова, а требовались лишь условия его существования. Такое обстоятельство не встречается в системах с другими дополнительными общими или частными интегралами, ранее изученными.
Заключение
В статье разными способами получены три вида стационарных движений для механической системы, в которой существует частный интеграл Стеклова. Установлено несколько свойств для таких систем. Прежде всего они составляют вычислительные аспекты применения метода Рауса — Ляпунова. Они показывают как можно упростить процесс нахождения стационарных движений в подобных и других системах. Основное свойство выражается последней теоремой, согласно которой все стационарные движения можно получить при неполной связке интегралов.
Список литературы О стационарных движениях механической системы с частным интегралом Стеклова
- Аппель П. Теоретическая механика. Москва: ГИФМЛ, 1960. Т. 2. 487 с.
- Уиттекер Э. Т. Аналитическая динамика. Ижевск: Изд-во Удмурт. ун-та, 1999. 584 с.
- Голубев В. В. Лекции по интегрированию уравнений движения тяжелого твердого тела около неподвижной точки. Москва: Регулярная и хаотическая динамика, 2002. 287 с.
- Стеклов В. А. Один случай движения тяжелого твердого тела, имеющего неподвижную точку. (сообщение в заседании Харьковского математического общества 5 марта 1893 г.) // Сочинения. Москва: Тип. М. Г. Волчанинова, 1896. 9 с.
- Routh E. J. A treatise on the stability of a given state of motion, particulary steadly motion. London: McMillan, 1877. 108 p.
- Routh E. J. The advanced part of a treatiseon the dynamics of a system of rigid bodies. London: McMillan, 1884. 343 p.
- Ляпунов А. М. О постоянных винтовых движениях твердого тела в жидкости // Собрание сочинений. Москва: Изд-во АН СССР, 1954. Т. 1. С. 276-319.
- Иоффе А. Д., Тихомиров В. М. Теория экстремальных задач. Москва: Наука, 1974. 480 с.