ТЕОРЕТИЧЕСКАЯ МЕХАНИКА
УДК 531.36
DOI: 10.18101/2304-5728-2018-3-22-39
ОБ УСТОЙЧИВОСТИ СТАЦИОНАРНОГО ДВИЖЕНИЯ МЕХАНИЧЕСКОЙ КОНСЕРВАТИВНОЙ СИСТЕМЫ1
© Новиков Михаил Алексеевич
доктор физико-математических наук, старший научный сотрудник,
Институт динамики систем и теории управления
им. В. М. Матросова СО РАН
Россия, 664033, г. Иркутск, ул. Лермонтова, 134,
E-mail: nma@icc.ru
В статье исследуется устойчивость стационарного движения нелинейной механической консервативной автономной системы, описывающей вращение твердого тела вокруг неподвижной точки. Для исследуемой системы известны три первых общих интеграла: энергии, момент количества движения, Пуассона. При равенстве Аппельрота, связывающем моменты инерции тела с координатами центра масс, допускается частный интеграл Гесса. Исследуемое стационарное движение имеет место и при существовании интеграла Гесса. Исследование устойчивости проведено по уравнениям линейного приближения возмущенного движения. Оно опирается на существование только нулевых и чисто мнимых корней характеристического уравнения с соответствующими им простыми элементарными делителями. Это выражается системой трех неравенств от коэффициентов характеристического уравнения, притом двукратный нулевой корень имеет простые элементарные делители. Анализ трех неравенств, выражающих чисто мнимые простые корни характеристического уравнения, позволил выделить семь областей решений. Отдельно рассмотрены случаи вырождения характеристического уравнения: появление дополнительных нулевых и кратных чисто мнимых корней.
В частности, установлена неустойчивость в линейном приближении при условии существования частного интеграла Гесса. Показана необходимость применения системы аналитических вычислений.
Ключевые слова: устойчивость движения; интеграл уравнений движения; характеристическое уравнение: элементарный делитель.
Введение
При проектировании механических объектов часто проводится исследование динамических свойств. В качестве таких свойств могут рассматриваться, в частности, выделение стационарных или установившихся движений и установление их устойчивости [1, 2]. Исследование устойчивости опирается на известные методы Ляпунова [1-4]. В настоящей статье рассматривается механическая система,
1
Работа выполнена при частичной поддержке Совета по грантам Президента Российской Федерации для государственной поддержки ведущих научных школ (проект № НШ-8081.2016.9).
описываемая дифференциальными уравнениями: Ap = (B-C)qr + z0у2, уj = ry2 - qy3,
< Bq = (C-A)rp + x0y3-z0yu у2 = ру3 -Ни (1)
Cr = (A- B)pq - х0у2, Уз = ЯУ\ ~ РУ2 ■
Здесь х0 ^ 0 ^ z0, А,В,С — моменты инерции твердого тела относительно главных осей Ox, Оу, Oz\ — проекции мгновенной
угловой скорости на подвижные оси; х0, z0 — координаты центра масс в тех же осях; У\,У2>Уз — проекции ортов подвижных осей на неподвижную вертикальную ось OZ .
Для системы (1) известны первые общие интегралы [5]
V0 = Ар2 +Bq2 +Сг2 +2х0уг +2 z0y3 = const, Fj = Арух +Bqy2 +Cry3 = const, V2=y2l+y22+y23=\.
При условии Аппельрота
AC(x2 +z2) = B(Ax2 + Cz2) (2)
система (1) допускает частный линейный интеграл Гесса [5] V3 = Ах0р + Cz0r = 0. В связи с этим представляет интерес нахождение
условий устойчивости некоторых стационарных движений и в случаях, когда не выполняется равенство Аппельрота.
1. Постановка задачи
Рассмотрим одно из стационарных движений для системы (1)
tjA2xl+C2zl V
С (Ах2 +Cz2) -Cz0
= 0, Гю
Га =■
A(C-A)x0z0> ' /10 Ja2x2+C2z2'
^ (1.1)
А(Ах0 +Cz0) Ах0
720 - Узо
Здесь рассматриваются движения, для которых С > А, х0г0 > 0. Ставится цель исследования устойчивости стационарного движения (1.1). Оно имеет место и при существовании частного интеграла Гесса.
Связки из первых интегралов возмущенного движения не смогли составить знакоопределенную квадратичную форму. Вторым методом Ляпунова в этом случае не удалось установить устойчивость стационарного движения (1.1). Ввиду сложности вычислительной задачи предложим получение достаточных условий устойчивости (1.1) простейшей линеаризованной системы уравнений движения.
2. Необходимые условия устойчивости по линейному приближению
Проведем исследование устойчивости (2) по уравнениям первого приближения. Составим отклонения х1 = р — р0~, х2 = q —
= Р~Ро, Х2 Х3 = Г~Г0, ХА=Ух ~Ую> =У2 -У20;
х6 = у3 — у30. Необходимые условия устойчивости будем получать по
характеристическому уравнению, составленному из дифференциальных уравнений линеаризованной системы возмущенного движения,
х = Дх, (2.1)
где матрица линейной части имеет вид
( 0 (В -С)г0/А 0 0 z0/A 0 ^
(С-А)г0/В 0 (<С-А)р0/В -zJB 0 xJB
0 04- -B)pJC 0 0 -xJC 0
0 ~ Узо 0 0 ro 0
У 30 0 -Ую ~го 0 Po
1 0 Ую 0 0 -Po 0 J
Характеристическое уравнение для матрицы 1)} будет таким:
/(Я) = det(D, -Щ = X2fx (Я) = 0, (2.2)
где (Я) = Я4 + k2?i + к0 при обозначениях
к2 =-1 ={А2[С2 -С(А + В) + 2АВ] х4 +
ABC(C-A)x20z20yjA2x20 +C2z20
+ АС[2(С - А)2 +B(A + C)]x2z2 +С2[А2 -A(B + C) + 2BC]z*j, k _ (Ах2 + Cz2 )s0 (xQ, Zq )^(x0 , Zq ) 0 A2BC2(C ~A)(A2xg +C2zg) ' ^(x0, z0) = Л3х04 + 4^4C(C - A)x2z2 - C3z04,
(x0, z0) = AC(x2 +z2)- B(Ax2 + Cz2).
Здесь и далее требуется огромное число различных вычислений, связанных с алгебраическими операциями над матрицами, факторизацией полиномиальных выражений. Это выполняется системой аналитических вычислений «Mathematica» на персональных компьютерах.
При непростых элементарных делителях нулевого корня (2.2) независимо от остальных корней уравнения fx (Я) = 0 стационарное движение (2.1) будет неустойчивым [3, 4, 6]. Известно [7], что
24
элементарные делители матрицы 1)} и любой другой матрицы, полученной из элементарными преобразованиями над строками и столбцами, будут одинаковыми. Для этого умножим 1)} слева и справа на соответствующие невырожденные матрицы 1\ и 0] . Полагаем, например, матрицы преобразования такими:
Г1 0 х0/г0 14 0
0 1 0 0 0
0 0 1 0 0
0 0 »43 1 0
0 и 152 0 0 1
1° 0 0 - Сг0/(Ах0) 0
Г 1 0 0 0 0
0 1 0 0 0 0
а = 0 0 0 1 0 1 0 0 1 0
0 0 0 1 0
а(1) \Чб\ 0 а(1) Убз а(1) ^ 64 0 ъ
где обозначено:
С = ^
х0 + С г0
(С-А)(Ах20+С220)
(^Х Q 2 Q
г(1) =
»43
^А
х0 +С
1
АС(Ах2 +Сг2)
(С А)х0г0
Г(1) =
'52
В
^А
х0 + С г0
=
Ах0 (2 Ах^ +
С12
2 _3
Чзл г0
АХд (АХд И" )
С5(С-А)г
В-А
№
х0 + С г0
1
С*0(Ах2 +С22)
А(С - А)хд
Ч{61 =
Р2Х20+С2220 \ 25
А3 (С - А)Хд (Ах2 +Сг20)
СЧ5
О) _ (Ах0 + 2Сг0) С г0
4А
Х0 + С
С (С - А)г0 (Ах1 +Сг1)
Ахп
В результате получается
ада = д =
Го 0 0 0 0 0
0 0 0 0 0 х0/В
0 0 0 0 - х0/С 0
0 ё(2) 42 0 0 0 0
0 0 ё(2) 53 0 0 0
0 0 0 0 1
где
ЛГ) _ _(Х0 ; 20 )
42 I-
(С - А)х0 у]А2х2 + С
V
о
— С2 2 ^(2) = _^ 20
Ах^ А х^ ~С ^
2 _2 0
Такое представление матриц Д, Тх, (7, обусловлено выбором базовых строк и столбцов матрицы Д. Вместе с тем операцию «сжатие-растяжение» можно выполнить как левосторонним, так и правосторонним преобразованием матриц. Матрица Д имеет равный двум дефект. В предположении лп (х0, гп) ^ 0 тогда двукратный нулевой корень уравнения (2.2) имеет простые элементарные делители [7].
3. О достаточных условиях устойчивости
Исследование достаточных условий устойчивости линейных консервативных систем сводится к необходимости существования чисто мнимых или нулевых корней характеристического уравнения с всеми простыми элементарными делителями [8]. При ранее установленных простых нулевых корнях уравнения (2.2) для дальнейшего исследования необходим анализ корней уравнения (Я) = 0 , который сводится здесь к неравенствам
к2>0,к0>0, т2=к2-4к0>0. (3.1)
Свободный член уравнения (Я) = 0 обращается в нуль, поэтому удобнее осуществить замену переменных
В =
АС(х2 + г20) Ах0 + С^о
+ а
(3.2)
и вещественный параметр а может принимать значения разных знаков.
Выполняя подстановку (3.2) и умножая т2 на
АВ2С(С - А)х2г2 (А2х2 + С2г2) > О, получим Д = Ъ2а2 + Ъха + Ь0, где
Ъ2 = (.Ах2 + Сг2)2[АС(х2 + г20)2 +12 (С-А)2х^2),
\ = 2 (Ах2 + Сг2){АС(Ах2 + Сг2)2 +
+ (С -А)2 [(6 А2 + 4АС)х2 + (4 АС + 6С2)г2]},
Ь0 = АС[(Ах2 + Сг2)2 + 4(С - А)2 х2г2 ]2. Тогда система неравенств (3.1) запишется в виде
Д>0, Д>0, Д>0. (3.3)
Условия устойчивости оказались сведенными к системе неравенств, и поэтому прежде изучим неравенство > 0.
4. Дополнительные требования к условиям устойчивости
Рассмотрим возможность существования системы неравенств
А2-А(В + С) + 2ВС<0, С2-С(А + В) + 2АВ>0. (4.1) Вместе с тем должны учитываться основные свойства твердого тела А + В > С, А + С > В . Последнее свойство В + С > А выполняется тождественно ввиду С > А. Предварительно разрешим первое неравенство (4.1) относительно А :
(В + С-^)/2<А<(В + С + ^)/2, (4.2)
где = [(5 — С)2 — 4ВС] — дискриминат первого квадратичного неравенства (4.1). Для существования решений квадратного неравенства здесь необходимо (1Х > 0. Учитывая строгость неравенства в (4.2), должно быть (1Х Ф 0. Последнее условие с/1 > 0 приводит к необходимости анализа двух решений
\.В>(3 + 2л[2)С, (4.3)
2.0 < В <(3- 2у[2)С. (4.4)
В предположении решения (4.3) проверим выполнение свойства А + С > В. Составим выражение А + С -В, которое должно быть положительным. Подставляя здесь ограниченные значения В,
вычисления получают А + С-В<А + С-(Ъ + 2л/2)С =
= А — 2(42. + 1)С. Ввиду начально заданного С>А последнее выражение только отрицательно, чего не должно быть в твердом теле. Для решений В из неравенства (4.4) проверим выполнение свойства
А + В > С. Аналогично составим А+В — С < (В+ С + +
+ В — С — [3В — С + yfd^]/2 — Sí. Для оценки знака Sí сравним их рациональные и иррациональные части. Здесь выполняется ЗВ -С <3(3- 2у[2)С -С = 2(4 - 3^2)С < 0 ввиду 3^2>4. Далее имеет место (ЗВ - С)2 —dl= 8В2 > 0 . Отсюда заключаем S] < 0. что
соответствует А + В < С. Так как это недопустимо в твердом теле, то система неравенств (4.1) несовместна. Отметим, что также несостоятельны системы неравенств А2 — А(В + С) + 2ВС < 0, С2 - С (А + В) + 2АВ > 0 и А2-А(В + С) + 2ВС<0, С2-С(А + + В) + 2АВ > 0. Система уравнений А2 - А(В + С) + 2ВС = 0, С2 — С (А + В) + 2АВ = 0 несовместна, так как она может существовать только в двух случаях: 1.В = 0, 2.В = С . Здесь второй случай приводит
к А = 0, и оба эти случая соответствуют плоскому телу. Аналогично рассмотрим вторую систему неравенств
А2 - А(В +С) + 2ВС > 0, С2 - С(А + В) + 2АВ < 0. (4.5) Вычитая из второго неравенства первое, получим необходимое для
выполнения (4.5) условие (С - А)(А + С - ЗВ) < 0) . При С > А отсюда следует
А + С-ЗВ<0. (4.6) Второе неравенство (4.5) разрешим относительно С :
(A + B-JT2)I2<C<(A + B + JT2)I2, (4.7)
где d2 = (В — А)2 — ЛАВ, которое должно быть положительным ввиду строгого второго неравенства в (4.5). Здесь могут существовать два возможных решения:
l.B>(3 + 2yÍ2)A, (4.8)
2.0<В<(3-2л[2)А. (4.9)
Для первого решения (4.8) проверим выполнение свойства А + С > В . Составим выражение A + C-B<A + (A + B + 4d~2)l 2-
-B = [3A-B + y[d^]/2 = S2. Так как здесь 3A-B<3A-(3 + < 0, то знак рациональной части S2 будет отрицательным. Далее (ЗА -В)2 -d2= 8 А2 > 0, откуда следует S2 <0 . Следовательно, решение (4.8) для системы (4.5) не состоятельно.
Другое решение В из (4.9) должно, в частности, удовлетворять условию (4.6). Здесь следует оценка А + С-ЗВ > А + (А +В-
- -Jd^)/2 - ЗВ = [ЗА - 5В - sfd^y2 = S3. Легко показать, что здесь ЗА-5В > ЗА-5(3-2-Í2)A = 2(5*j2 -6)А > 0. Для проверки знака S3
составим разность (ЗА - 5В)2 — d2 = 8(А2 - ЗАВ + 3В2) = S4 . Так как квадратичная форма относительно А, В в последнем выражении положительно определена, то S4 > О, и отсюда следует S3> 0, что
приводит к А + С — 3В > О . Но последнее не согласуется с (4.6) и поэтому система (4.5) несовместна. Аналогичные выкладки приводят к невозможности существования системы неравенств А2— А(В + С) + + 2ВС > О, С2 -С(А+В) + 2АВ<0.
Рассмотрим следующую систему: А2 - А(В + С) + 2ВС < О, С - С (А + В) + 2АВ < 0. Складывая оба неравенства, получим (С — А)2 + В(С + А ) <0. что невозможно для вещественных положительных А,В,С . Таким образом, для выполнения 1'\ > 0 следует рассматривать систему
А2-А(В + С) + 2ВС>0, С2-С(А + В) + 2АВ>0. (4.10) При этом может допускаться обращение в нуль только одного из неравенств (4.10).
5. Полный набор требований к получению условий устойчивости
Составим наиболее полный набор для выполнения системы неравенств (3.3). Введем в рассмотрение уравнение
/1(л/^) = Л(^) = 0. (5.1)
С учетом ранее установленных в предыдущем разделе свойств полный набор возможных требований будет таким:
1.С2 -С(А + В) + 2АВ > 0 (коэффициент при х4 в выражении 1'\ должен быть положительным),
2. А2 -А(В + С) + 2ВС > 0 (коэффициент при z4 в выражении 1'\ должен быть положительным),
3.А + В — С > 0 (свойство моментов инерции твердого тела),
4.А + С — В > 0 (свойство моментов инерции твердого тела),
5.-Fj >0 (сумма корней уравнения (5.1) должна быть отрицательной),
6.F0>0 (произведение корней уравнения (5.1) должна быть положительным),
7.F2 >0 (корни уравнения (5.1) только вещественны).
Применим здесь подстановку (2.4), тогда в том же порядке требования запишутся:
1.а(Ах2 + Cz¡)(2A - С) + С[А2х2 + (2А2 - 2АС + C2)z¡ ] > 0,
2.а(Ах2 + Cz¡)(2С -А) + А[(А2 - 2АС + 2C2)x¡ + C2z¡] > 0,
3.а(Ах2 + Cz2) + A2x¡ + С(2А-C)z2 > 0,
4.а(Ах2 +Cz2)-(A2x2 +C2z¡)> 0,
5 ,a(Ax2 +Cz2)[2(A2X2 +C2Z2)-AC(X2+Z2)] + AC[(AX2 +CZ2)2 + + 4(C - A)2 x2z2 ] > 0, 6.a(p(x0,z0)< 0, l.b2a2 +bla + b0 > 0.
Дискриминант биквадратного уравнения ф(х0, z0) = 0 обозначим
о > о >
б/, = А2С2 (4А2 - 7АС + 4 С2) > 0. Дискриминант квадратного
относительно параметра а уравнения Ъ2а +Ь1а + Ь0= 0 равен
Й?4 = -36(С - А)2 х0\>Оо, ).
Шестое требование оказалось факторизованным. Поэтому основное разделение условий при анализе семи перечисленных требований проведем для положительных и отрицательных значений параметра а. Для краткости изложения введем обозначения
_ С[А2х2 + (2А2 - 2АС + С2)г2] (Ах1+Сг1)(С-2А) '
-А[(А2 - 2АС + 2С2)х2 + C2z¡] < Q
(Ах2 + С22)(2С-А)
-А2х2 + С(С - 2А)г1 А2х2 + С%2 ^ п
аъ=--—-----, а4=-\-^->0,
д Сг0 -1- С^д
-ЛС[(Лх2 + Сг2)2 + 4(С - Л)2х0%2] °5 (Лх2 + Сг2)[ Д2Л - С)х2 + С(2С - Л)г2]'
Вместе с тем введем точечные значения относительного параметра / = > 0. Наиболее часто используемой величиной при анализе является больший положительный корень уравнения ф(х0, г0) = 0 :
Д2(С - А) + ^4А2 - 7АС + 4С2 ] п t%--—->0. В пятом и третьем
_ А(С - 2 А)
требованиях будут применяться величины t5--,
С (2С А)
t3 =
С (С - 2 А)
Перечисленные величины ai (/ = 1,...,5), / , ( / = 3,5,8) введены для
выбора из них наибольших или наименьших значений при решении
30
систем неравенств. В зависимости от знака (С —2А) некоторые из них участвуют в оценках нижних и верхних точных границ.
6. Получение решений системы неравенств (3.3)
Проведем параметрический анализ решений системы (3.3) в терминах А,, и вначале положим а > О. Из шестого требования тогда необходимо (р{х0, г0) < 0, что приводит к выбору значений t > ts . Далее
рассмотрим первое распределение моментов инерции тела А < С < 2А . При любых положительных а тождественно выполняются первое, второе, третье и пятое требования. При этом седьмое требование заведомо выполнено при всех положительных слагаемых многочлена
Ъ2а2 + Л,а + Ь0. Четвертое требование ограничивает сверху значения
параметра а. При отсутствии других ограничений одну часть решений системы неравенств (3.3) можно записать:
А<С<2А, t>ts, 0<а<а4. (6.1)
Для другой возможности распределения моментов инерции С > 2А тождественно выполняются второе и седьмое требования. В пятом
требовании знак выражения [/1(2/1 - С)х2 + С(2С — А)г2 ] зависит от отклонения (t — t5). Сравним значения t5,t8 :
_ А[(2А2 - ЛАС + ЗС2) + (2С - А)^4А2 -1АС + 4С2]
8~5 С2(2С-А) '
В последнем выражении рациональная часть представляет положительно определенную относительно А, С квадратичную форму, и
положительную иррациональную часть. Поэтому < ^, и при
положительных а величину t5 можно не учитывать. Следовательно, в
этом случае пятое требование выполняется тождественно. Тогда можно получить нетривиальную область решений системы (3.3) при выполнении условий
t>ts, С>2А, шах[0,а3] < а < т[а„а4].
Для сравнения величин а1,а4 составим разность ах—а4 = 2А2С — 2 А, которая в данном случае положительна. Следовательно, верхней гранью является а4 . При оценке знака а3 используем :
_ - А[4А2 - 7АС + 2С2 + (С - 2А)л]4А2 -1 АС + 4С2 ] 3 8 (С-2А)
Здесь рациональная часть выражения в квадратных скобках является знакопеременной квадратичной формой, а иррациональная часть положительна. Отрицательный знак этого выражения может оказаться в
31
единственном случае, когда
С>2А, 4А2 -7АС + 2С2 <0,
14А2 - 7АС + 1С21> (С - 2А)-\14А2 - 7АС + 4С2.
Легко подсчитать
(4А2 -7АС + 2С2)2 -(С-2А)2(4А2 -7АС + 4С2) = = -АС(С - А)(5С -12А).
Решением последних неравенств является 2А < С <1215А.
Следовательно, при 2А < С < 12/5 А выполняется t3> ts, а при С > 12/5 А будет t3 <ts. Тогда при значениях 1<13 получим а3 < 0, а для 1> 13 следует а3> 0. Окончательно при 2А < С < 12/5 А величина а3< 0 на отрезке I е (/8;73), и а3 >0 при С > 12/5^4 .
Следующие решения системы неравенств (3.3) можно записать: 1Л3 >Ц : 2А<С < 12/5 А, [¿8 <t<t3; 0 < а < > 13\а3 <а<а4],
(6.2)
2.t3<ti■.C>12/5A,a3<a<a4. (6.3)
Далее рассмотрим другую часть значений параметра а : а < 0. Шестое требование здесь приводит к (р{х0, гп) > 0. откуда областью определения следует I е (0; ). Ввиду отрицательности дискриминанта ¿/4 и Ъ2 > 0 седьмое требование выполняется тождественно. Четвертое требование можно не принимать в расчет вследствие а < 0 < а4 (что заведомо выполняется). Здесь также выделим два распределения моментов инерции тела:
\.А<С<2А, 2.С>2А.
Для первого распределения сразу понадобится установить требования на параметр: а> ах, а> а2, а> а3 (аг < 0, а2 < 0) . В пятом требовании имеется две возможности:
1)^ < (где должно выбираться а <а5 (а5 > 0)),
2)t > t5 (где выбирается а> а5 (а5 < 0)).
Первая возможность может не учитываться для отрицательных значений а, где ^ получается суммой двух положительных слагаемых. Решение системы (3.3) тогда может существовать при условиях:
А<С<2А, тах[а1,а2,а3]<а<0, и
I е (/5;78), А<С < 2 А, тах[а1,а2,а3,а5]< а< 0.
Сразу установим соотношения
_ 2(С-А)(А2 - АС + С2) _ _ 2А(С-А) (С - 2А)(2С - А) 1 3 С-2А
а2-а3 =
_-2 (С-А)
2 С-А
При первом распределении А < С < 2А их легко упорядочить: ах< а2< а3< 0 . Для значений I е (I-; ) нужно учитывать а5:
_ - (С - А)[2А3хд + АС(ЗС - А)х2г2 +2 С2 (С- ] < Л5 а3 (^х02+Сг02)[Д2^-С)х02+С(2С-^)г02]
(здесь все слагаемые в скобках положительны). Следовательно, на отрезке (0;^8) для рассматриваемых А < С < 2А точной нижней гранью
значений параметра а является а3. Тогда можно записать решение системы неравенств (3.3):
А<С < 2А, 7 е (0, ), а3 < а <0. (6.4)
Рассмотрим другой случай распределения моментов инерции С > 2А . Здесь прежде всего следует учитывать только те значения параметра t, при которых а3< 0. Если при этом > . то гарантируется значение
а3< 0. Это выполняется, как ранее установлено, при 2А < С < 12/5А.
_-А(С-А)(С-5А)
Далее опорной точкой является t — --.
5 3 5 С (С - 2А)(2С - А) Отсюда можно заключить, что > I- при 2А < С < 5 А, и t3<t5 при С > 5А. В промежутке (0; I-) пятое требование выполняется тождественно при всех отрицательных значениях а как сумма положительных слагаемых. Тогда возможные решения системы (3.3) могут быть при выполнении условий тах[а2,а3] < а < тт[0,а1]. Точно
так же в промежутке (^О условиями будут тах[а2,а3,а5] <а< тт[0,а1].
При С > 2А легко установить: ал > 0, си > а2 > а;. Тогда решения системы неравенств (4.5) в этом случае можно записать:
\Л3 >Г8 : 2А<С < 12/5А, ^е (0;78), а3<а< 0, (6.5)
2Л3 е(?5;?8): 12/5А<С < 5А, te(0;t3), а3<а< 0, (6.6)
ЪЛ3 е(0;Г5):С >5А, Ге(0;Г3), а3<а< 0. (6.7)
Интересно отметить, что первое и второе требования хотя и участвовали в анализе, но не отразились в условиях устойчивости.
7. Исследование вырожденных случаев
Рассмотрим возможные особые случаи при анализе системы неравенств (3.3). Они возникают при обращении в нуль свободного члена уравнения (Я) = 0, когда получается кратный нулевой корень
(Я) = 0, и при существовании кратных чисто мнимых корней. При этом кратный нулевой корень уравнения (Я) = 0 существует в двух случаях: \.а = 0 (50(х0,г0) = 0), 2.<р(х0,г0) = 0.
В первом случае при а = 0 выполняется равенство Аппельрота (1.2) и возникает четырехкратный нулевой корень уравнения (2.2). Легко видеть, здесь дефект матрицы /А равен трем. Следовательно, кратному нулевому
корню при а = 0 соответствуют не все простые элементарные делители, что приводит к растущему решению системы (2.1) [3, 4, 6] и неустойчивости стационарного движения (1.1).
Для анализа других случаев вырождения будем предполагать $0(х0,г0)^0, и составим Я - матрицу (Д - XI С). В результате элементарных операций, что сводится к линейным операциям над строками и столбцами и попросту выражается умножением (Д — XI С) слева на матрицу Г(Х) и справа на матрицу О(Х), получим
С
Д(Я) =
0 с!(3) 12 0 0 0 0
0 0 0 0 0 х0/В
0 0 0 0 -х0/С 0
< 0 0 0 0 0
0 0 0 < 0 0
0 0 ё{ъ) 63 0 0 1
л
где
с!(3> =
41
^(3) _ Оо 7 2о )
12 I-
Р
(2 2 0
б/(3) = 54
+сч2 V -1
Ах0 +
А3С(С-А)х3020
/(3) = 16ъ
^А2х20+С2г20 V
ЗА
(Ах20+Сг20У
АС(С-А)х30г0
^2х02 + С%2 V А
С(С-А)(Ах20+С220У _
А хп
= Ъ6ъ\
Ахп
3^(х0,^^А14Тс14]1С(С-А)230(АХ20+С220)
г-ф4(Л) = Ь41Л-ф4(Л),
Левостороняя матрица преобразования здесь имеет вид
Т(Х) =
Г 1 0 Сг0/(Ах0) 0 0
0 1 0 0 0 0
0 0 1 0 0 0
¿41 ¿42 ¿43 ¿44 ¿45 ¿46
0 ¿52 -С/х0 0 1 0
- ВХ/х0 ¿63 — г0/х0 0 к
где
¿41
А (-V
{-А2ВСх0г0 х
(х0,)^2х02+С\2 ЦС- А)2о(Ах2 + Сг2)
х(Ах2 +Сг2) + Ах0Л/А2х2 +С2г2 [А3Вх4 + АС(ЗА2 - АВ- ЛАС + + ЗВС + С2)х20г20 +С2(2А2 -2АВ-2АС+ ЗВС)г4]Х2-2г0(Ах2 +
+ Сг2)*0(х0,г0)(А2х2+С2г2)},
¿42
-ВХ
3 (С-А)20(Ах20+С220У
~[АС(С-А)^А2х1+С2г1 х0г0 X2 + А3 х4 +
+ АС(А + С)х2г2 +С3г4],
143
-1
У! Схл
{АБС2 (С- А) х
3*0 (х0, )р2х2,+С221 ЦС- А)г0 (Ах2 + Сг2 )5 х(Ах2 +Сг2)х0г214 +Сгоу/А2х2 +С2г2 [А2 (ЗАВ -2АС -2ВС + 2С2)х04 +
+ АС(А2 + ЗАВ - ЛАС - ВС + ЗС2 )х2г2 + ВС3г4]Х2 +2х0(Ах2 +Сг2 )х
+СЧ2)},
1
¿44
3 (С-А)(Ах2 + Сг2 у
-[АСл/А2х2 +С2г2 х0г0Х2 + А2 (ЗС -А)х4 +
+ 2АС(2С-А)х2г2 + С3г4],
¿45 = Ь
АС (С А)х0г0 \ 22
/ л 2 , 2 \ 3 V О
¿52
¿46
(Ах2+С22) 1
р2х20+С2г20 \
хп
2 2
>
сг0 (Ах2+сг2)
А (С ~А)хд
3(С - А)г0(Ах1 + Сг2})'
~[АСр2х1 +С2г02 х0г0 X2 + А3 х4 +
+ 2АС(2А - С)х2г2 +С2(ЗА- 2], 1
^61 =
АС(С-А)г0
$ о (хо >2 о х0 + С г0
+ 2оЯо
х0 (Ах0 + Сг0)
[АВх0^¡А2х2 +С2г20 X2 +
163
1
АС(С-А)г0
(х0,)^А2х20+С2220 \ хо(Ах2 + сг2) — х050(х0,г0)].
Правосторонняя матрица О запишется:
[ВСх0у/А2х2 +С2г2 X2 -
№) =
'Чи 0 Сг0/(Ах0) 0 0
421 1 0 0 0 0
031 0 1 0 0 0
041 042 043 1 0 0
051 052 - С^/х0 0 1 0
4061 062 ^63 г0/х0 0 к
где 0п
Зя0(х0, г0)
-ЗВ + С) х2г2 + С2 (А - В) г40 ],
_ Х$]А2х20+С2г20
(Ах20 +Сг20)[-АВСх0г0 ^А2х2 +С2г2 X2 + АС(2А -
(хо >2о)
А х„
А3С(С-А)х30г0
^4x0 + С'2 0
031
(Ах0 + С20) (х0, )
[АВСх0г0 л)а2х20+С2г20 X2 +
+ А2 (С - В)х4 + АС(А -ЗВ + 2С)х2г2 +2 С2(А-В)г4], -1
041
(х0, г0 х0 + С г0
АС (С А)х0 {а2всх222{ах2 +
2о (Ахо + Сг0)
+ Сг20 )Х4+А х2 ^2х2 +С2г20 [А2 (С- В)х40 + АС(2С +В- 2 А)х20г20 ] X2 + + г0 (4АХд - Сг20) (А2х2 + С2г20 )г0э0 (х0,20)(Ах2 + Сг2)},
143
С3(С-А)2о
л]а2х1 + С2г1 \ Ах0(Ах2 + Сг2)
[-АхоЛ]А2х2 +С2г20 X2 +
+ 2 г0(Лх02+Сг02)],
АХ
0 42
С(А-2В)г0Х2
А х0 + С 2.0
32о$о(хо,2о)(АХ0 +Сг0)
[~АВСх0г0 у1А2х20+С2г2 X2 +
+ А2 (В- С)х40 - 2 АС (С - А)х20г20 + С2 (А- В)г4],
В-А
р2х20+С2220 \ -1
С20(Ах2+С22)
А(С - А)х30
Зя0 (х0, г0 х0 + С г0 ^
А3(С-А)х1
2 {АБС 2х2г2 (А2х2
Сг^+Сг^)
+ С2г20)Х4 +Сг^А2х1 +С2г20 [А2(С-2В)х40 -АС(2А + В - 2С)х2г2 + + С2(В - А)г4]Х2 -г0(Ах2-4 Сг20)(А2х20 + С\2 К(х0,г0)}, _ (АВх20 +С(А-В)г1)Х
062 ~д 2 2ч '
х0 (Лх0 + Сг0 ) 1
Г 63
р2х20+С2220 \
АС(С-А)г0
х0 (Ах0 + Сг0)
[-Сг0Л/^2х02 + С2г2 X2 +
+ 2 х0(Ах2 +Сг2)].
Перестановкой строк и столбцов с приведением коэффициентов ап,а54,Ь43,Ь63 к единице можно представить матрицу /А(л) в более простой диагональной форме, которую проще записать в терминах инвариантных множителей [7]. В случае ф(хп, г0) = 0 инвариантные
множители будут такими: 1,1,1,1, Л, Л1 (X2 +/?,). где
= (Ах2 + Сг2 )[(Ах2 + Сг2)2 + 4(С - А)2 х2г2 ] х
х АС (С - А)х0г0 (х2 + г2 )^А2х2 +С2г20 .
Легко видеть, что при полученных здесь инвариантных множителях четырехкратному нулевому корню соответствуют непростые элементарные делители. Следовательно, при ф(х0, гп) = 0 стационарное движение (1.1) неустойчиво.
В случае кратных чисто мнимых корней уравнения (Я) = 0 при
Х1/2 = +/д//?,/2 инвариантные множители будут: 1,1,1,1, Л, Х(Х2 + И]!2)1. Здесь элементарные делители для Я^ тоже не являются простыми.
Таким образом, в линейном приближении [3, 4] стационарное движение (1.1) является неустойчивым при всех вырождениях уравнения (Я) = 0 .
Учитывая, что условия устойчивости по начальному приближению нелинейной системы имеют значительно большие ограничения, чем условия устойчивости линеризованной системы, можно заключить о неустойчивости стационарного движения (1.1) при существовании интеграла Гесса. При этом неустойчивость движения (1.1) не зависит от знака z0 .
Заключение
В результате проведенных вычислений составлены области устойчивости (6.1)—(6.7) линеаризованной системы дифференциальных уравнений возмущенного движения для механической системы (1). Параметрический анализ устойчивости проведен в зависимости от трех параметров: а (отклонение от равенства Аппельрота), t (положительный
параметр z20lx20), относительный положительный параметр С!А . При
условии Аппельрота (2) стационарное движение (1.1) для линеаризованной системы (1) оказалось неустойчивым. Следует отметить, что даже в случаях, когда центр масс расположен ниже плоскости Оху, стационарное движение (1.1) системы (1) неустойчиво.
Всюду при анализе требовался большой объем вычислений, особенно при факторизации алгебраических выражений. С применением системы аналитических вычислений удалось значительно ускорить и упростить вычислительные операции.
Литература
1. Ляпунов А. М. Общая задача об устойчивости движения. Собрание сочинений. М.-Л.: Изд-во АН СССР, 1956. Т. 2. С. 7-263.
2. Четаев Н. Г. Устойчивость движения. Работы по аналитической механике. М.: Изд-во АН СССР, 1962. 535 с.
3. Каменков Г. В. Устойчивость движения, колебания, аэродинамика. М.: Наука, 1971. Т. 1.255 с.
4. Каменков Г. В. Устойчивость и колебания нелинейных систем. М.: Наука, 1972. Т. 2. 213 с.
5. Голубев В. В. Лекции по интегрированию уравнений движения тяжелого твердого тела около неподвижной точки. М.: Регулярная и хаотическая динамика, 2002. 287 с.
6. Малкин И. Г. Теория устойчивости движения. М.: Наука, 1966. 530 с.
7. Гантмахер Ф. Р. Теория матриц. М.: Наука, 1967. 576 с.
8. Бахтин А. Б., Брюно А. Д., Варин В. П. Множества устойчивости многопараметрических гамильтоновых систем // Прикладная математика и механика. 2012. Т. 76, вып. 1. С. 80-133.
ON STABILITY OF STATIONARY MOTION OF MECHANICAL CONSERVATIVE SYSTEM
Mikhail A. Novikov
Dr. Sci. (Phys. and Math.), Leading Researcher, Institute of Systems Dynamics and Theory of Management named after V. M. Matrosov SB RAS 134 Lermontova St., Irkutsk, 664033, Russia E-mail: nma@icc.ru
The article studies the stability of stationary motion of a nonlinear mechanical conservative autonomous system that rotates a rigid body around a fixed point. For the system under study, the first three common integrals are known: energy, moment of momentum, Poisson. At Appelrot equality that connects the moments of body inertia with the coordinates of the center of mass, a particular Hess integral is allowed. The steady-state motion under research is also valid for the existence of the Hess integral. The research of stability is carried out by the equations of linear approximation of perturbed motion. It is based on the existence of only zero and purely imaginary roots of the characteristic equation with the corresponding simple elementary divisors. This is expressed by a system of three inequalities from the coefficients of the characteristic equation, moreover, the double zero root has simple elementary divisors. The analysis of three inequalities expressing purely imaginary simple roots of the characteristic equation, enabled us to identify seven areas of solutions. The cases of degeneration of the characteristic equation are separately considered: the appearance of additional zero and multiple purely imaginary roots.
Thus, the instability in the linear approximation is revealed under the condition of the existence of a particular Hess integral. The necessity of using the system of analytical computations is shown.
Keywords: stability of motion; integral of motion equations; characteristic equation; elementary divisors.
References
1. Lyapunov A. M. Ob shay a zadacha ob ustoychivosti dvizheniya. Sobraniye sochineniy [General Problem of Stability of Motion. Collected Works]. Moscow: Leningrad: USSR Acad. Sci. Publ. 1956. V. 2. Pp. 7-263.
2. Chetayev N. G. Ustoychivost dvizheniya. Raboty po analyticheskoy mekhanike. [Stability of Motion. Works in Analytical Mechanics], Moscow: USSR Acad. Sci. Publ. 1962. 535 p.
3. Kamenko G. V. Ustoychivost dvizheniya, kolebaniya, aerodynamika [Stability of Motion, Oscillations, Aerodynamics]. Moscow: Nauka Publ. 1971. V. 1. 255 p.
4. Kamenkov G. V. Ustoychivost i kolebaniya nelimeinykh sistem [Stability and Oscillations of Nonlinear Systems]. Moscow: Nauka Publ. 1972. V. 2. 213 p.
5. Golubev V. V. Lektsiyi po inyegrirovaniyu uravneniy dvizheniya tyazhelogo tverdogo tela okolo nepodvizhnoy tochki [Lectures on Integration of Equations of Motion for a Rigid Body about a Fixed Point], Moscow, 2002. 287 p.
6. Malkin I. G. Teoriya ustoychivosi dvizhenia [Theory of Stability of Motion], Moscow: Nauka Publ. 1966. 530 p.
7. Gantmacher F. R. Teoriya matrits [Theory of Matrices]. Moscow: Nauka Publ., 1967. 576 p.
8. Bakhtin A. B., Bruno A. D., Varin V. P. Mnozhestva ustoichivosti mnogoparametricheskikh gamil'tonovykh sistem [Sets of Stability of Multiparametric Hamiltonean Systems]. Prikladnaya matematika i mekhanika - Applied mathematics andMtchanics. 2012. V. 76, Iss. 1. Pp. 80-133.