CC BY
19Нелинейная динамика. 2015. Т. 11. № 2. С. 353-376. Полнотекстовая версия в свободном доступе http://nd.ics.org.ru
ПЕРЕВОДНЫЕ СТАТЬИ
УДК: 531.36
MSC 2010: 34D20, 70F25, 70Q05
Динамика систем с сервосвязями. I
■*
В. В. Козлов
В работе обсуждается динамика систем с сервосвязями Бегена, когда связи реализуются посредством управляемых сил. Классические неголономные системы представляют важный частный случай. Особое внимание уделено исследованию движения на группах Ли с левоинвариантной кинетической энергией и левоинвариантной связью. Наличие симмет-рий позволяет свести динамические уравнения к замкнутой системе дифференциальных уравнений с квадратичными правыми частями на алгебре Ли. В качестве примеров рассмотрено вращение твердого тела с левоинвариантной сервосвязью — проекция угловой скорости на некоторое фиксированное в теле направление равна нулю (обобщение неголо-номной задачи Суслова), а также движение саней Чаплыгина с сервосвязями определенного вида. Динамика систем с сервосвязями Бегена богаче и разнообразнее по сравнению с более привычной динамикой неголономных систем.
Ключевые слова: сервосвязи, симметрии, группы Ли, левоинвариантные связи, системы с квадратичными правыми частями
1. Введение
Работа посвящена динамике систем с сервосвязями. Основы теории систем со связями, реализуемыми с помощью управляемых сил, обсуждаются в диссертации Бегена [1]; черты общей теории приданы Аппелем [2]. Связь теории сервосвязей Бегена с основными принципами динамики изложена в работе [3]. В современной литературе по управляемым системам сервосвязи принято называть программой движения.
* Перевод статьи "The Dynamics of Systems with Servoconstraints. I", опубликованной в журнале Regular and Chaotic Dynamics, 2015, vol. 20, no. 3, pp. 205—224.
Получено 10 февраля 2015 года После доработки 05 марта 2015 года
Работа выполнена за счет средств гранта Российского научного фонда (проект 14-50-00005).
Козлов Валерий Васильевич kozlov@pran.ru
Математический институт им. В. А.Стеклова РАН 119991, Россия, г. Москва, ул. Губкина, д. 8
Уравнения движения с сервосвязью
. = п (
Ф(х,х,$ = 0 (||/0) (1.1)
имеют вид
(1 дТ дТ
(М дх дх ^ ^ ^'
Здесь х = (х1,хп) — обобщенные координаты, х = (х1,..., хп) — скорость системы, Т — кинетическая энергия, ^ — внешняя сила, N — заранее заданное ковекторное поле (поле сил), Л — множитель Лагранжа. Уравнения (1.1) и (1.2), конечно, следует рассматривать совместно.
Условие реализации сервосвязи имеет следующий вид:
(1.8) Здесь
дх 2
— положительно определенный оператор инерции. Условие (1.3) означает возможность найти параметр Л как функцию от х, х и Ь (не решая уравнения движения). Обсуждение условия (1.3) можно найти в [3].
Здесь следует дать некоторые пояснения. Общую теорию управляемых механических систем со связями не следует отождествлять с теорией сервосвязей Бегена (как это сделано, например, в работах [14, 15]).
Вообще, системы со связями, зависящими от управляющих параметров, рассматривались и до классической работы Бегена [1]. Упомянем работы Я. И. Грдины по динамике живых организмов (см. [17], где можно найти ссылки на его более ранние работы). В самом
простом варианте уравнения с управляющим параметром имеют вид
= ф(*>ж>*>Л) = °-
Функция Ь м Л(Ь) (управление) находится в результате совместного решения этой дифференциально-алгебраической системы. Такие же уравнения изучались в работах [14, 15]. Если
д*, о
9 А ^ '
то из уравнения связи (по крайней мере локально) параметр Л находится как функция от состояния х, х и времени Ь. Подставляя эту функцию в выражение для силы, получим замкнутую систему дифференциальных уравнений второго порядка, из которой находим движение Ь м х(Ь), а тем самым и управление Ь м Л(Ь). Если же уравнение связи не содержит управляющего параметра, то дифференцируя его по времени и используя уравнение Лагранжа, получим:
где Ф — известная функция от х, х и Ь. Если
то (снова по теореме о неявной функции) параметр Л находится в виде функции от состояния системы и времени. В частности, если в некоторый момент Ф = 0, то это равенство справедливо при всех значениях Ь. Наконец, когда внешняя сила также не зависит от параметра Л, то дифференциально-алгебраическая система, как правило, несовместна и для реализации движения со связью надо вводить дополнительную силу (реакцию связи). Зависящие от управляющих параметров связи были названы в работе [14] условными связями. Сила Е считается внешней силой, а реакция условной связи не вводится. Точнее, она считается равной нулю и поэтому не совершает работы на возможных перемещениях 5х, определяемых уравнением Четаева
Такие же возможные перемещения использовались и в работах [15-17].
Более того, в указанных статьях уравнения движения неоправданно усложняются введением дополнительного множителя Лагранжа:
При этом утверждается, что для замыкания системы необходимо задать (или из каких-либо дополнительных условий определить) закон управления
Л = Л(Х, х, Ь).
На наш взгляд, в этих уравнениях надо положить ц = 0 (как это и сделано в [14]).
В теории Бегена связи не зависят от управления, а действующая на систему сила равна сумме внешней силы Е (также не зависящей от параметра) и реакции сервосвязи ЛМ. Чтобы сервосвязь была идеальной, возможные перемещения удовлетворяют другому уравнению —
М, 6х) = 0.
Хотя это обстоятельство ясно выражено Бегеном (а также повторено в классическом учебнике Аппеля [2]), оно не всегда оказывается понятым с точки зрения основных принципов аналитической механики. Вообще, определение возможных перемещений — это ключевая аксиома динамики систем со связями. Изменение этой аксиомы (при сохранении других аксиом — принципа освобождаемости и условия идеальности связей) приводит к иной динамике системы со связью (Т, Е, Ф) (обсуждение этого круга вопросов см. в [3]). Сама собой возникает аналогия со знаменитым «пятым постулатом» евклидовой геометрии.
Ясно, что сервомоторы реализуют движение механической системы со связью с некоторой погрешностью. Анализ этих погрешностей, а также вопросы их минимизации обсуждаются, например, в работах [15, 18]. Для динамических систем общего вида некоторые подобные вопросы рассмотрены в книге [19].
Реализацию движения систем с сервосвязями можно назвать активной; она осуществляется с помощью регулируемых устройств. Если положить
ох
то условие реализации (1.3) заведомо выполнено и уравнения (1.1)-(1.2) будут описывать классическое неголономное движение. С другой стороны, как хорошо известно, неголоном-ное движение можно реализовать и пассивным способом, используя дополнительные силы
анизотропного вязкого трения. Эта идея восходит к Каратеодори. Обзор результатов по пассивной реализации связей и их обсуждение можно найти в книге [20].
Задача точного интегрирования уравнений (1.1)—(1.2) существенно сложнее соответствующей задачи для уравнений Гамильтона и неголономных уравнений (когда N = дФ/дх). В типичной ситуации здесь уже не справедлива известная теорема об изменении кинетической энергии
гЛ М
дх
= (Е,х),
(1.4)
из которой выводится интеграл энергии в стационарном случае и когда внешние силы потенциальны. Соотношение (1.4) справедливо лишь в том случае, когда для всех «действительных» скоростей х (удовлетворяющих (1.1))
х) = 0.
Другими словами, действительные скорости должны лежать среди возможных перемещений. Это условие выполняется лишь в частных случаях, например, когда
Ф =
где ц — ненулевая функция в фазовом пространстве. Если ковекторное поле N не зависит от х, то получаем однородную неголономную связь, для которой соотношение (1.4) заведомо выполняется.
С другой стороны, нетрудно сформулировать обобщенную теорему Нётер, дающую линейные по скорости первые интегралы. Пусть 'ш(х) — векторное поле на конфигурационном пространстве, поток которого (однопараметрическое семейство преобразований конфигурационного пространства) сохраняет кинетическую энергию Т(х, х) (после естественного продолжения этих преобразований на фазовое пространство). Если
(N,^=0, (1.5)
то
1(1'") = <*»■
Соотношение (1.5) означает, что поле симметрий в каждой точке конфигурационного пространства есть возможное перемещение нашей системы. В частности, если система движется по инерции Е = 0, то уравнения движения допускают интеграл
(§"«') =сопв^
С качественной точки зрения системы с сервосвязями ведут себя сложнее, чем более привычные нам неголономные системы. Обсуждение этого круга вопросов, а также рассмотрение динамики систем с сервосвязями на группах Ли составляют содержание настоящей работы.
2. Системы на группах Ли с левоинвариантной кинетической энергией и левоинвариантной сервосвязью
Пусть G — n-мерная группа Ли, g — ее алгебра Ли. Пусть T — левоинвариантная риманова метрика на G. С точки зрения динамики эту метрику можно интерпретировать как кинетическую энергию некоторой механической системы. Согласно Пуанкаре [4], динамические уравнения (уравнения Лагранжа) отделяются и представляют собой замкнутую систему дифференциальных уравнений на алгебре g (или на коалгебре — двойственном линейном пространстве g*):
m k = c)kmiWj, k = 1,...,n. (2.1)
Здесь wi,...,wn (декартовы координаты на алгебре g) — квазискорости рассматриваемой системы, а Ш\,...,mn (координаты на коалгебре g*) — импульсы. Они связаны линейными соотношениями
mk = ^ IkjWj, 1 ^ k ^ n, (2.2)
где \\IkjУ — постоянный тензор инерции. Постоянные j = -cikj — структурные постоянные алгебры Ли; они удовлетворяют известным тождествам Якоби. Кинетическая энергия
Т = \(1и,и) = \(Г1т,,т)
— первый интеграл уравнений (2.1). Для полного описания движения системы к уравнениям (2.1) следует добавить кинематические уравнения. Подробности см. в [5].
Замечание. Скорость системы ш = (wi,..., шп) — касательный вектор к группе G в единице. Его компоненты преобразуются по контрвариантному правилу. Поэтому (с последовательной точки зрения тензорного анализа) номера компонент скорости ш следовало бы записывать сверху. Однако это создает неудобства при рассмотрении конкретных примеров.
Пусть
Ф = (а, ш) = ^ aiwi = 0 (ai = const) (2.3)
— левоинвариантная связь, линейная по скоростям. Постоянный ненулевой ковектор а — элемент из коалгебры g*. Предположим (согласно Бегену), что эта связь реализуется с помощью управляемой силы
ХЬ = (ХЬ\,..., Xbn),
где X(t) — неизвестная функция времени (управление), а Ь Е g* \ {0}. Тогда динамика системы с левоинвариантной сервосвязью принимает следующий вид:
mk = ^ cjkmiWj + Xbk (1 ^ k ^ n), ^ aiwi = 0. (2.4)
Если а и Ь коллинеарны, то система уравнений (2.4) задает неголономное движение со связью (2.3). В этом случае уравнения (2.4) допускают интеграл энергии T = const. Если же ковекторы а и Ь не коллинеарны, то энергия, как правило, не сохраняется. Условие реализации связи (2.3) сводится к неравенству
(I-1а, Ь) = 0. (2.5) _НЕЛИНЕЙНАЯ ДИНАМИКА. 2015. T. 11. №2. С. 353-376_^
Оно означает, что ковекторы а и Ь не ортогональны во «внутренней» метрике на д*, индуцированной исходной левоинвариантной метрикой на группе С. Условие (2.5) — это общее условие реализации сервосвязи (1.3), представленное в групповых переменных.
После исключения множителя Л из (2.4), систему на алгебре д (или на коалгебре д*) можно свести к замкнутой системе из п — 1 дифференциальных уравнений с однородными квадратичными правыми частями. Этого можно добиться разными способами. Укажем один из них. Для этого представим уравнения (2.4) как систему на коалгебре д* (другими словами, обратим линейную систему (2.2) и заменим в (2.4) и\,...,ип линейными функциями от Ш\,..., тп). Тогда уравнение связи примет вид
т = о, (2.6)
где А = (Л\,..., Ап) — ненулевой элемент из алгебры д. Линейной заменой переменных этот вектор приводится к виду (0,..., 0,1), то есть уравнение связи (2.6) примет вид тп = 0. Следовательно, последнее уравнение в динамической системе (2.4) сведется к алгебраическому, из которого множитель Л находится как квадратичная форма от т1,... ,тп. Подставляя это выражение для Л в первые п — 1 уравнений, получаем требуемое.
Замкнутую систему из р = п — 1 дифференциальных уравнений с квадратичными правыми частями будем называть приведенной.
3. Общие свойства приведенных систем
Итак, рассмотрим замкнутую систему дифференциальных уравнений в Мр = {х1,..., хр }
ху = Уу (хь...,хр), 1 ^ з ^ р, (3.1)
причем ы,. ..,Ур — однородные квадратичные многочлены от р независимых переменных.
Случай р =1 тривиальный. При р = 2 система (3.1) интегрируется в квадратурах. Это следствие более общего результата об интегрируемости автономных систем на плоскости с однородными правыми частями (восходящего еще к Лейбницу). Поскольку случай р = 2 будет представлять для нас особый интерес, приведем формулы, позволяющие, с одной стороны, проинтегрировать систему (3.1), а с другой — дать качественный анализ поведения ее траекторий. Для этого положим г = х2/х\. Тогда
/у»2 х1
где
f (г) = У2(1,г) — гы(1,г) (3.3)
— многочлен от г третьей степени.
Предложение 1. Если Ь ^ х^) — решение системы (3.1), то Ь ^ —х(—Ь) также решение той же системы.
Это — простое следствие квадратичности правых частей и теоремы единственности решений. Для обычных линейных систем с постоянными коэффициентами это не так: здесь Ь ^ —х(Ь) также будет решением. Ясно, что —х(—Ь + с) — решение системы (3.1) при всех вещественных с.
Заметим, что нетривиальные состояния равновесия системы (3.1) (то есть точки х0 = 0, где v(x0) = 0) не могут быть изолированными. Если v(x0) = 0, то v(Лx0) = 0 для всех вещественных Л. Качественные свойства решений системы (3.1) существенно зависят от существования нетривиальных равновесных состояний.
Теорема 1. Если v(x) = 0 при х = 0, то система (3.1) допускает инвариантную прямую х = Л£, где £ — ненулевой вектор из Мр, а параметр Л изменяется со временем согласно дифференциальному уравнению
Л = сЛ2, с = 0. (3.4)
Это утверждение и его обобщения играют существенную роль в первом методе Ляпунова для сильно нелинейных систем [6]. Укажем кратко основные моменты его доказательства.
Снабдим Мр «стандартной» евклидовой структурой и рассмотрим гауссово отображение
г/ \ -
х | ► 1 (ж; — ,
1Кх)Н
х £ Бр-1 (единичная сфера в Мр с центром в нуле). Имеем гладкое отображение единичной сферы в себя. Индекс векторного поля V в точке х = 0 равен степени отображения Г. Степень Г будет четной, поскольку каждый элемент из 1т Г имеет четное число прообразов. Следовательно, степень Г отлична от ( — 1)р-1 и поэтому (по известной теореме из дифференциальной топологии) Г имеет неподвижную точку х = £. Но тогда система дифференциальных уравнений (3.1) допускает инвариантную прямую х = Л£, Л £ М. Подстановка х = Л£ в уравнения (3.1) дает
Л £ = v(Л£) = ЛМ0 = Л2|К£)||£.
Отсюда получаем уравнение (3.4), где с = |^(£)|| = 0. Что и требовалось.
Из теоремы 1 вытекает ряд важных следствий.
1°. Если х = 0 — единственная особая точка квадратичного векторного поля V, то она неустойчива. Более того, имеется решение со сколь угодно малым начальным условием, которое уходит на бесконечность за конечное время. Последнее вытекает из вида решения уравнения (3.4):
А«> ^ т^Ь А» ^ А<°>-
2°. Фазовый поток такой системы определен не на всей числовой оси времени Ь. Однако если равновесие х = 0 квадратичного векторного поля V не изолировано, то оно может быть устойчивым и при этом фазовый поток определен на всей оси М = {Ь}. Вот простой пример:
х 1 = х1х2, х2 = —х2.
Здесь вся ось х1 =0 состоит из положений равновесия.
3°. В условиях теоремы 1 система (3.1) не допускает положительно определенного первого интеграла (в частности, в виде положительно определенной квадратичной формы по х1,...,хр). В противном случае равновесие х = 0 было бы устойчивым. В частности, на решениях исходной динамической системы (2.4) на алгебре Ли энергия не сохраняется.
Забегая вперед, скажем, что имеются примеры систем на группах Ли с левоинвариант-ными связями Бегена, приведенные системы которых удовлетворяют условию теоремы 1. Наоборот, равновесные состояния приведенных неголономных систем на группах Ли никогда не бывают изолированными.
Теорема 2. Если ковекторы а и Ь коллинеарны, то тривиальное равновесие х = 0 приведенной системы устойчиво и имеется проходящая через начало координат прямая, целиком состоящая из положений равновесия.
Действительно, в этом случае левоинвариантная связь (2.3) будет неголономной и поэтому динамическая система допускает интеграл энергии — положительно определенную квадратичную форму на алгебре д. Ее ограничение на гиперплоскость (2.3) даст нам положительно определенный квадратичный интеграл приведенной системы. Следовательно, тривиальное равновесие х = 0 устойчиво по Ляпунову. Наличие прямой равновесий вытекает из теоремы 1.
Замечание. Не следует думать, что если ковекторы а и Ь не коллинеарны, то приведенная система имеет лишь тривиальное равновесие х = 0. Ответ на этот вопрос зависит (кроме всего прочего) от строения группы Ли С. Например, если группа С абелева (ск = 0), то при всех а и Ь приведенная система вырождается в тривиальную: х = 0. Все точки Мр являются состояниями равновесия, притом устойчивыми.
При р = 2 теорема 1 сразу же выводится из формул (3.2) и (3.3), поскольку многочлен третьей степени над полем вещественных чисел всегда имеет вещественный нуль. Здесь следует дать пояснения. Если ось х1 =0 не является инвариантной, то многочлен У1(1,г) имеет вторую степень. Тогда многочлен f из (3.3), действительно, третьей степени. Кроме того, в правой полуплоскости (где х1 > 0) отношение г = х2/х1 может принимать все вещественные значения.
Мы можем утверждать больше: система (3.1) на плоскости с единственным тривиальным равновесием может иметь либо одну, либо три инвариантных прямых (считая кратности). Это наблюдение можно распространить на многомерный случай. Справедлива
Теорема 3. Пусть снова у(х) = 0 при х = 0. Если количество инвариантных прямых (считая кратности) системы (3.1) конечно, то оно нечетно и не превосходит 2р — 1.
При р = 2 получаем, что система (3.1) имеет либо одну, либо три (с кратностями) инвариантных прямых. Теорема 1, очевидно, содержится в теореме 3 (хотя используется при ее доказательстве).
В отличие от топологического доказательства теоремы 1, доказательство теоремы 3 алгебраическое. Запишем в явном виде систему алгебраических уравнений для нахождения инвариантных прямых:
у1(х1,..., хр) — ¡х1 = 0,..., Ур(х1,..., хр) — ¡хр = 0. (3.5)
Это — система р однородных уравнений степени 2 относительно р +1 переменных х1,...,хр и ¡. Согласно теореме Безу, количество различных комплексных непропорциональных ненулевых решений системы (3.5) (считая кратности) в точности равно 2р. Одно из вещественных решений системы (3.5)
х1 = ... = хр = 0, л = 1 (3.6)
не дает никакой инвариантной прямой. Следовательно, число инвариантных прямых не превосходит 2р — 1. Поскольку у(х) = 0 при х = 0, то остальные ненулевые решения алгебраической системы (3.5) непропорциональны решению (3.6). Одно из вещественных решений алгебраической системы (3.5) существует по теореме 1. Остальные решения либо вещественные, либо пары комплексно-сопряженных решений. В любом случае общее количество вещественных непропорциональных решений (считая их с кратностями) нечетно. Что и доказывает теорему 3.
В заключение этого параграфа сделаем несколько замечаний о системах с квадратичными правыми частями на плоскости, допускающих нетривиальные положения равновесия. Такие системы, конечно, не типичны. Оказывается, после подходящей замены времени они сводятся к линейным системам.
Теорема 4. Пусть р = 2 и система (3.1) имеет целую прямую ах\ + Ьх2 = 0 (а2 + Ь2 =0), состоящую из положений равновесия. Тогда .заменой времени £ м- т,
система (3.1) приводится к линейной системе.
Прямой равновесий без ущерба для общности можно считать, например, ось Х2. Этого всегда можно добиться подходящей линейной заменой переменных, не меняющей структуры и вида системы (3.1). При таком выборе переменных уравнения (3.1) принимают следующий вид:
очевидно, состоят из траекторий исходной системы. Только направление движения в исходной системе зависит еще и от знака переменной Х\.
В качестве примера рассмотрим линейную систему (3.8) с положением равновесия в виде устойчивого узла (рис. 1а). На рисунке 1Ь показан фазовый портрет исходной нелинейной системы (3.7). В отличие от линейной системы, тривиальное равновесие нелинейной системы уже неустойчиво. Стоит подчеркнуть, что смена устойчивости при переходе от линейной системы к нелинейной происходит не всегда.
Теорема 5. Пусть р = 2 и система (3.1) допускает первый интеграл Е в виде невырожденной квадратичной формы. Тогда эта система имеет целую прямую равновесных состояний.
йт = (ах1 + Ьх2)йЬ,
Х1 = х1(Ах1 + Вх2), Х 2 = х1(Сх1 + Бх2).
Замена времени йт = х1(М приводит эту систему к линейной. Траектории линейной системы
х'1 = Ах1 + Вх2, х'2 = Сх1 + Бх2,
(3.7)
(3.8)
(а)
(Ь)
Рис. 1
Если интеграл Е положительно (отрицательно) определен, то это утверждение отмечено ранее в работе [7]. В этом случае теорема 5 — простое следствие теоремы 1. Приведем доказательство теоремы 5 для квадратичного интеграла с сигнатурой (1,1). Невырожденной линейной заменой переменных интеграл Е можно привести к виду
х1х2. (3.9)
Пусть в новых переменных система (3.1) имеет следующий вид:
х 1 = а1х1 + Ь1 х1х2 + с1х2, х 2 = а2х2 + Ь2х1х2 + с2х2. (3.10)
Условие Е = 0 эквивалентно следующим соотношениям на коэффициенты этой системы:
а2 = С1 = 0, а1 + Ь2 = 0, Ь1 + С2 = 0.
Таким образом, (3.10) принимает следующий явный вид:
х 1 = х1(а1х1 + Ь1х2), х 2 = —х2(а1х1 + Ь1х2). (3.11)
Значит, прямая а1х1 + Ь1х2 = 0 целиком состоит из положений равновесия. Если а1 = Ь1 = 0, то система вырождается в тривиальную х 1 = х2 = 0, все решения которой — равновесные
Отметим, что после замены времени
йт = (а1х1 + Ь1х2)йЬ (3.12)
система (3.11) переходит в линейную гамильтонову систему с гамильтонианом Е. Фазовый портрет системы (3.11) изображен на рисунке 2 (направление движения по траекториям, конечно, может быть противоположным). Портрет фазовых траекторий для случая положительно определенного квадратичного интеграла хорошо известен (см. [8], где рассмотрена известная неголономная задача Суслова).
4. Замечания о фазовых портретах на плоскости
Рассмотрим приведенную систему на плоскости (р = 2); она имеет вид (3.10). В типичной ситуации система (3.10) имеет единственное равновесное состояние х1 = х2 = 0. Более точно, это свойство имеет место при выполнении следующих неравенств:
(а1С2 — С1 а2 )2 + а1С1 Ь2 + Ь{а2С2 — а^Ь2С2 — Ь^а2Ь2 = 0, (а2 + а2)(с? + с2) = 0.
Как выглядят фазовые портреты таких систем (с единственным тривиальным равновесием)? Как мы уже видели в §2, система (3.10) имеет либо одну, либо три (считая кратности) инвариантных прямых, проходящих через начало координат — положение равновесия.
состояния. Что и требовалось.
Используя формулы (3.2)—(3.3) и их модификации, можно дать классификацию возможных видов фазовых портретов системы (3.10). Более того, используя проективные координаты, можно даже получить информацию об уходе фазовых траекторий в бесконечность (как это делалось, например, в [9] для систем на плоскости с полиномиальными правыми частями). Однако автору неизвестно, была ли проделана эта работа. Обычно исследуются более сложные системы с наличием линейных по х1 и х2 слагаемых (см. [9]).
Тем не менее, исчерпывающее представление о возможных видах фазовых портретов дают рисунки 3а-3е, на которых изображены фазовые портреты системы
х 1 = х2 + х2, х 2 = (к + 1)х1 + (2 — к)х1 х2 — х2, (4.1)
содержащей вещественный параметр к.
Рис. 3
Система (4.1), очевидно, удовлетворяет условию теоремы 1. При всех k она имеет инвариантную прямую xi = Х2. Если k > 0, то других инвариантных прямых нет, а при k < 0 добавляются две инвариантные прямые
х2 = (-1 + </^к)хг и х2 = (-1 - V^k)xi. (4.2)
При k = 0 («бифуркационное» значение параметра) эти прямые совпадают. Рисунки 3a-3e дают ясное представление о перестройке фазового портрета при прохождении параметра через точку бифуркации.
При положительных значениях k фазовые траектории будут двоякоасимптотическими к инвариантной прямой xi = x2. Они имеют вид «волны», горб которой неограниченно увеличивается при k ^ +0. При k < 0 появляется новая пара инвариантных прямых (4.2), между которыми все фазовые траектории в левой полуплоскости стремятся к положению равновесия xi = x2 = 0 при t ^
5. Вращение твердого тела с левоинвариантной сервосвязью
Это простой, но важный пример развиваемой теории; его рассмотрение намечено еще в работе [3]. Здесь группа Ли G — трехмерная группа SO(3) вращений трехмерного евклидова пространства. Наличие евклидовой метрики на алгебре SO(3) позволяет отождествлять ковекторы и векторы.
Уравнения (2.4) на алгебре SO(3) принимают следующий вид:
Itjj + и х Iu = Xb, (а,и) = 0. (5.1)
Символ «х» обозначает векторное произведение, I — тензор инерции твердого тела, постоянные ненулевые векторы a и b задают направления фиксированных в теле осей l и n (которые упоминаются в [3]).
Оси подвижного ортогонального трехгранника выберем так, чтобы третья ось совпала с осью l. В этих осях вектор a имеет компоненты 0, 0, 1. Обозначим Ij компоненты тензора инерции, а ui, и2, j — проекции вектора угловой скорости и на выбранные оси. В этих обозначениях из = 0 будет уравнением связи, а уравнения (5.1) примут следующий явный вид:
Ill jj 1 + Ii2 и 2 + (Ii3 и 1 + I23 jj 2)U2 = Xbi,
I12jj 1 + I22jj2 - (I13jj 1 + I23jji)wi = Xb2, (5.2)
Il3 jj 1 + I23 jj2 + (I12 j1 + I22U2)U1 - (I11 j1 + I12 W2)W2 = Xb:i.
Здесь мы уже учли уравнение связи j = 0. Чтобы получить приведенную систему на плоскости, надо из (5.2) исключить множитель X.
Если b1 = b2 =0 (а b3 = 0), то получим уравнения классической задачи Суслова из него-лономной механики. Ее фазовый портрет хорошо известен (см. [8]). Чтобы почувствовать разницу в качественном поведении системы в общем случае (когда b2 + b2 =0), положим I13 = I23 = 0. Из последнего уравнения найдем множитель Лагранжа
X = I12и2 + (I22 - I11 )и1и2 - I12w|. (5.3)
В общем случае эта квадратичная форма — не тождественный нуль. Из первых двух уравнений (5.2) легко выводится линейный интеграл
(b2111 - b1112)и1 + (b2I12 - b1I22)j2 = const.
Рис. 4
Этот интеграл не вырождается при естественном условии Ь\ + &2 = 0 (поскольку /ц/22 — /12 > 0). Его линии уровня прямые. Фазовый портрет изображен на рисунке 4. Две жирные прямые (состоящие из положений равновесия) — это прямые, на которые распадается уравнение Л = 0. Подчеркнем, что Л — нейтральная квадратичная форма от ш\ и С2. Указанное направление движения по прямолинейным участкам фазовых траекторий соответствует положительному значению произведения /ц (61/22 — &2/12). Если /ц = 0, то положения равновесия (где Л = 0) будут заполнять координатные оси ш\ = 0 и С2 = 0.
Видно, что фазовый портрет на рисунке 4 радикальным образом отличается от портрета неголономной задачи Суслова (который при /13 = /23 = 0 сплошь состоит из положений равновесия). В частности, если /ц = /22, то равновесие = 02 = 0 всегда неустойчиво.
Перейдем теперь к обсуждению общего случая. Начиная с этого момента мы будем предполагать, что
61 = 0.
Этого всегда можно добиться поворотом ортогональной системы координат на подходящий угол вокруг третьей оси (совпадающей с осью I).
Согласно [3], условие реализуемости связи
из = 0 (5.4)
сводится к условию, что неподвижные в теле оси I и п не ортогональны в метрике, задаваемой кинетической энергией тела (точнее, в метрике на двойственной алгебре (зо(3))*, которая задается тензором /_1). Это же условие можно получить по-другому, исключая из двух последних уравнений (5.2) множитель Лагранжа (первое уравнение не содержит множителя Л, так как 61 = 0 по предположению). После этого получаем систему двух уравнений относительно и 02, и условие реализуемости связи (5.4) сводится к условию разрешимости этой системы относительно производных Со 1 и Со2. В итоге приходим к неравенству
63(/11 /22 — /12) — 62(/11 /23 — /12/13) = 0.
(5.5)
Это и есть искомое условие реализуемости связи (5.4). В (5.5), конечно, &2 + 63 = 0.
Если 62 = 0, а 63 = 0, то условие (5.5) заведомо выполнено. Это наблюдение — частный случай более общего результата о реализуемости обычных неголономных связей. Наоборот, если 63 = 0, а 62 = 0 (оси l и n ортогональны относительно евклидовой метрики SO (3), индуцированной метрикой исходного трехмерного евклидова пространства), то условие реализуемости (5.5) принимает вид
/11/23 = /12/13.
Это условие, конечно, может нарушаться при некоторых специальных распределениях массы твердого тела.
Найдем все равновесные состояния приведенной системы. Для этого запишем в явном виде уравнения (5.2) после исключения множителя Лагранжа:
/11 0 1 + /12 С02 + (/13 СО 1 + /23 02)02 = 0, (5.6)
(63 /12 — 62 /13)01 + (63 /22 — 62 /23 )о2 — [(63 /13 + 62/12 )о1 + (63 /23 + 62/22)02 ]01 +
+ 62 (/1101 + /1202)02 = 0.
(5.7)
Из уравнения (5.6) вытекает, что в состоянии равновесия либо j = 0, либо
I13 и1 + I23 и2 = 0. (5.8)
Рассмотрим сначала первую возможность. Уравнение (5.7) дает в этом случае соотношение
(b3I13 + b2I12 )и2 = 0.
Таким образом, если
(b3I13 + 62I12) = 0, (5.9)
то приведенная система не допускает нетривиальных положений равновесия. Наоборот, если
(b3I13 + 62I12) = 0,
то вся ось и>2 =0 целиком состоит из положений равновесия. В этом случае приведенная система (5.6)-(5.7) сводится к линейной системе после замены времени
dr = j2dt.
Этот результат предсказывается теоремой 4.
Как выглядит фазовый поток соответствующей линейной системы? Мы ответим на него в частном случае, когда I12 = I13 = 0. Линейная система имеет следующий вид:
I11 j + I23 и2 =0,
1 , (5.10)
(b3I22 - b2123)j + (b2I11 - b2I22 - b3^3)j = 0.
Она допускает квадратичный интеграл
I11 [b2(I11 - I22) - b3I23\j2i - I23(b3122 - b2123)= 0. (5.11)
Фазовые траектории системы (5.10) (как и исходной нелинейной системы) состоят из линий уровня этой функции. Пусть b3I23 = 0. Тогда всегда найдутся интервалы изменения отношения b2/b3, на которых квадратичная форма (5.11) определенная (положительно или отрицательно), а также целые интервалы, когда форма (5.11) нейтральная. В первом случае фазовый поток схож с фазовым потоком неголономной задачи Суслова, а во втором — имеет вид, изображенный на рисунке 2.
Рассмотрим теперь вторую возможность, когда выполнено равенство (5.8). Тогда (в предположении b2 = 0) из уравнения (5.7) вытекает соотношение
I12и2 + (I22 - I11)W1W2 - I12W2 = 0. (5.12)
Общее решение уравнения (5.8) имеет вид
и1 = kI23, и2 = -kI13 (k e R). Подставляя в (5.12), получим, что при условии
I12(I23 - Ii23) - (I22 - I11 )I13123 = 0 (5.13)
имеется только тривиальное равновесие (поскольку k = 0).
Сказанное можно суммировать в следующий вывод:
• если 62 =0 и выполнены неравенства (5.9) и (5.13), то приведенная система не допускает нетривиальных положений равновесия.
В частности, при этих условиях тривиальное равновесие с = С = 0 системы уравнений (5.6)-(5.7) неустойчиво. К исследованию этой системы можно применить общие результаты §§2 и 3. Отметим, в частности, что если вместо неравенства (5.13) имеем равенство, то система нелинейных уравнений (5.6)-(5.7) допускает целую прямую равновесий и, согласно теореме 4, также сводится к линейной системе после подходящей замены времени.
В заключение этого параграфа скажем кратко об условиях существования интегрального инварианта (инвариантной меры) с гладкой положительной плотностью у приведенной системы (5.6)-(5.7). Согласно [10], такой интегральный инвариант существует тогда и только тогда, когда фазовый поток такой системы сохраняет стандартную меру Лебега на плоскости. Другими словами, если представить систему (5.6)-(5.7) в стандартной форме
СС1 = ^1(^1,^2), СС 2 = ^2(^1,^2), то это условие сведется к условию бездивергентности:
^ + = 0. (5.14)
дс1 дс2
Поскольку функции г>1 и г>2 квадратичны по с и С2, то левая часть (5.14) — линейная однородная функция от С 1 и С 2. Следовательно, в итоге (5.14) сведется к условиям равенства нулю коэффициентов при С 1 и С 2.
Опуская элементарные вычисления, укажем итоговый результат. Приведенная система (5.6)-(5.7) допускает интегральный инвариант с гладкой положительной плотностью только при выполнении следующих двух условий:
62(2/22 — /11/12 + /?1 + /?з) + 6з(/12 /13 — /11 /23) = 0, (5.15)
62 (/12 /22 + /13 /23) + 63 (/22/13 — /23 /12) = 0. (5.16)
Таким образом, системы с интегральным инвариантом образуют подмножество коразмерности 2 в пространстве всех приведенных систем.
В частном случае, когда /ц = /13 = 0 (плоскость, натянутая на векторы а и 6, будет инвариантной для тензора инерции твердого тела), условие (5.16) выполнено автоматически, а условие (5.15) упрощается:
62 (/22 — /11) + 63/23 = 0.
При условиях (5.15) и (5.16) приведенная система допускает первый интеграл в виде ненулевого однородного многочлена третьей степени. Действительно, согласно Эйлеру, дифференциальная 1-форма
г2йш1 — г1йш2 (5.17)
будет полным дифференциалом первого интеграла исходной системы. После известных последовательных интегрирований формы (5.17) получаем однородный интеграл третьей степени И. Поскольку г1 = —дИ/дс2, г2 = дИ/дс1, то приведенная система принимает гамильтонову форму:
... _ дН . _ дн
— — т;—, Сс>2 — т;—• дС 2 дС 1
Замечание. Динамические системы с наличием нетривиальных равновесий, фазовые портреты которых изображены на рисунках 1Ь, 2 и 4, не допускают интегрального инварианта с гладкой положительной плотностью из-за того, что все их траектории являются асимптотическими (при Ь ^ или Ь ^ -ж). Более того, они даже не допускают инвариантной меры, абсолютно непрерывной относительно обычной меры Лебега на плоскости. Тем же свойством обладает и приведенная система в неголономной задаче Суслова (обсуждение см. в [8]). Но не стоит думать, что все системы с квадратичными правыми частями (не обращающимися в тождественный нуль), допускающие нетривиальные положения равновесия, обладают этим свойством. Вот простой контрпример: х 1 ^ х2, X2 = 0. Укажем еще пример системы с единственным тривиальным равновесием, фазовый поток которой сохраняет обычную меру Лебега: X1 — Х2 , X2 — X . Эта система является гамильто-новой; гамильтонианом служит функция
Н
Множество {Н — 0}, очевидно, инвариантно. Оно сводится к единственной прямой XI — Х2, существование которой гарантирует теорема 1.
3
3
1
2
6. Сани Чаплыгина как система с сервосвязью
Сани Чаплыгина — это твердое тело, движущееся по горизонтальной плоскости с неин-тегрируемой связью: скорость некоторой его точки (обозначим ее О) всегда ортогональна фиксированной в теле горизонтальной оси I. Классический неголономный вариант этой задачи хорошо изучен (см., например, [7, 11, 12]). Мы рассмотрим динамику саней Чаплыгина в предположении, что упомянутая неинтегрируемая связь реализуется управляемыми силами (по Бегену). Для краткости такую систему будем называть сервосанями.
Положение твердого тела на плоскости задается тремя параметрами: декартовыми координатами х, у выделенной точки О и углом поворота ф. Введем подвижную систему отсчета с началом в точке О, причем одну из этих осей направим по выделенной оси I. Пусть и, V — проекции скорости точки О на эти оси, и — угловая скорость. Уравнение неинтегрируемой связи имеет вид
V — 0. (6.1)
Мы будем рассматривать самый простой вариант, когда центр масс твердого тела совпадает с точкой О. Неголономная динамика такой системы тривиальна: точка О движется по окружности или по прямой (окружности бесконечного радиуса) с постоянной скоростью. Динамика саней с сервосвязью более содержательная.
Конфигурационное пространство сервосаней — группа движений плоскости Е(2). Кинетическая энергия равна
Т = Ш^ + У2) + 1-(6.2)
где т — масса тела, а I — момент инерции относительно вертикальной прямой, проходящей через точку О. И кинетическая энергия (6.2), и связь (6.1), очевидно, левоинвариантны. Поэтому мы можем применить общий подход, указанный в § 3.
Замечание. Ввиду известной ретракции группы БО(3) к группе Е(2), задачу о сервосанях можно рассматривать как предельный случай задачи о вращении твердого тела со связью из § 5 (эти вопросы обсуждаются, например, в [7]). Но мы проведем независимый анализ этой задачи, дополнив его геометрическим рассмотрением движения саней по плоскости.
Уравнения (2.4) с учетом связи (6.1) имеют следующий вид:
mil = шХа, 0 = —шиш + Х6, /0о = Хе.
Здесь Ха, Х6 — «следящие» силы, направленные вдоль подвижных осей, Хе — дополнительный момент сил, относительно вертикальной оси, проходящей через точку O; а, 6, е = const. Если а = е = 0, 6 = 0, то имеем классическое неголономное движение. Условие реализации сервосвязи в общем случае также сводится к неравенству 6 = 0.
Приведенная система описывается совсем простыми дифференциальными уравнения-
ii = аий, й = @иш ( а =
в =
тс Ib
(6.3)
Координатные оси и = 0 и 0 = 0 сплошь состоят из равновесных состояний. Система (6.3) имеет линейный интеграл
ви — aw = y = const. (6.4)
Мы будем обсуждать общий случай, когда обе постоянные a и в отличны от нуля. Фазовый портрет системы (6.3) изображен на рисунке 5. Расположение инвариантных прямых (6.4) и направление движения по траекториям соответствуют случаю, когда a > 0 и в > 0. Фазовый портрет на рисунке 5 есть частный случай более общей ситуации, изображенной на рисунке 4.
Уравнения (6.3) легко интегрируются в явном виде с использованием линейного первого интеграла (6.4). Исключая переменную 0 рис. 5 из первого уравнения, получим уравнение для отыскания линейной скорости точки O:
ii = и(ви — Y).
Предполагая, что константы в и Y положительны, укажем решение с траекторией в виде интервала 0 < и < Y/в:
i(t) =
Y
в + е
if
(6.5)
Это решение двоякоасимптотическое: при t — —с и при t — оно стремится, соответственно, к Y/в и 0. Решения вида (6.5) отвечают траекториям системы (6.3), расположенным в четвертом квадранте на рисунке 5 (где и > 0, а й < 0). При этом угловая скорость й также монотонно убывает от 0 до —j/а.
Простой анализ приведенной системы позволяет получить явное представление о скольжении сервосаней по горизонтальной плоскости. Мы по-прежнему будем считать, что все константы а, в и Y положительны.
Декартовы координаты x, y выделенной точки O, а также угол поворота ф саней могут быть найдены из следующих кинематических соотношений:
Х = и cos ф, y = и sin ф, ф = й.
Мы предполагаем, что угол ф отсчитывается от неподвижной оси х в положительном направлении; это угол между осью х и осью подвижного репера, не совпадающей с выделенной осью I.
Сосчитаем кривизну траектории точки O:
I¿У ~ ху\ (.х2+у2)3/2
1
а
/Зи — y
и
eYt(e + eYt)
aY
Важно отметить, что кривизна не зависит от угла поворота саней. Кривизна всегда положительна, и при Ь ^ —ж она экспоненциально быстро стремится к нулю, а при Ь ^ +ж, наоборот, экспоненциально возрастает. Траектория точки О изображена на рисунке 6. Ее длина от любой точки влево бесконечна, а вправо — наоборот, конечна. Это — простое следствие формулы для длины дуги
¿2 ¿2
/v^W = /«М* ti ti
и формулы (6.5).
Таким образом, в «отдаленном прошлом» тело практически движется по прямой не вращаясь, затем тело начинает вращаться с увеличивающейся угловой скоростью и траектория точки О начинает закручиваться по спирали. При бесконечно больших положительных значениях времени точка О становит-рис 6 ся неподвижной, а тело в итоге переходит во вращение с посто-
янной скоростью вокруг неподвижной вертикальной оси, проходящей через точку О. Конечно, точка О может двигаться по «завитку», изображенному на рисунке 6, в противоположном направлении. При этом сани совершают асимптотический переход от вращения вокруг вертикальной оси к поступательному движению по прямой. Такие движения отвечают начальным данным из второго квадранта фазовой плоскости приведенной системы. Вообще, завитки получаются друг из друга преобразованиями подобия и зеркального отражения. Это — простое следствие свойства однородности системы (6.3). При зеркальном отражении сани закручиваются в противоположном направлении.
Движения, уходящие в бесконечность (или приходящие из бесконечности) за конечное время, представляют меньший интерес, и их анализ здесь не приводится.
7. Сервосани на наклонной плоскости
Усложним задачу и рассмотрим динамику сервосаней на наклонной плоскости. В отличие от § 6, теперь на твердое тело еще действует сила тяжести. Уравнения движения усложняются, и теперь динамические уравнения не отделяются от кинематических соотношений (как это было в §3).
В обозначениях § 6 уравнения движения сервосаней имеют следующий вид:
mil = Xa — p sin p, 0 = —шиш + Л — p cos p, ICC = Xc. (7.1)
Здесь мы уже учли уравнения связи v = 0 и, кроме того, положили b =1 (это всегда можно сделать в предположении реализуемости связи v = 0). Угол поворота p отсчитывается от горизонтальной прямой на наклонной плоскости; p > 0 — вес тела (точнее, произведение веса тела на синус угла между наклонной и горизонтальной плоскостями).
Второе уравнение (7.1) служит для нахождения множителя Л. Подставляя формулу для Л в первое и третье уравнения (7.1), получим:
mil = шаиш + ap cos p — p sin p, I Co = mcuo + cp cos p. (7.2)
Добавляя кинематическое соотношение
ф = и
(7.3)
получим замкнутую систему автономных дифференциальных уравнений в трехмерном пространстве, 2п-периодическую по угловой переменной ф.
Если a = c = 0, то получим интегрируемую (и хорошо изученную) неголономную систему; для почти всех начальных данных неголономные сани не будут спускаться по наклонной плоскости, двигаясь в среднем по горизонтали. При c = 0 система (7.2)—(7.3) также легко интегрируется в квадратурах. Однако при c = 0 система (7.2)—(7.3) вряд ли относится к числу интегрируемых. В связи с этим становится актуальной задача ее качественного исследования. Первый шаг в этом направлении — поиск частных решений и исследование их устойчивости.
Если c = 0, то имеется семейство частных решений, когда ф = фо = const: сани спускаются по наклонной прямой с постоянным ускорением (если a cos фо = sin фо, то сани движутся с постоянной скоростью). Однако все они неустойчивы.
В общем случае, когда c = 0, имеется только одно такое движение, при котором ф = ±п/2: сани движутся по прямой наискорейшего спуска. Исследуем устойчивость такого движения, полагая для определенности ф = п/2 (случай ф = —п/2 ничем не отличается от этого).
Теорема 6. Если c > 0, то движение сервосаней по вертикальной прямой (когда ф = п/2) устойчиво по отношению к £ = п/2 — ф и
При c ^ 0 такое движение неустойчиво. Обсудим физический смысл условия теоремы. Если íp = то й = —р/т. Следовательно, при больших значениях времени линейная скорость саней отрицательна. Пусть сани движутся по инерции; тогда справа во втором уравнении (7.2) надо отбросить второе слагаемое. Если u(t) < 0, то положительный управляющий момент (c > 0) будет стабилизировать прямолинейное движение саней.
Для доказательства теоремы 6 нам потребуется
Лемма. Пусть t — u(t) — решение уравнения
с начальным условием и(0) ^ 0. Тогда и(Ь) < 0 в интервале времени (0, т), где ^(Ь) < 0. Действительно, функция п(-) (как решение (7.4)) выражается формулой
Эти значения отрицательны, коль скоро отрицательны значения функции ц(-).
Теорема 6 доказывается непрерывной индукцией по времени. Сначала заметим, что первое уравнение (7.2) имеет вид (7.4), где
здесь g = p/m — ускорение свободного падения, Л и ц — сложные функции времени (они зависят от ф и ф, которые, в свою очередь, суть функции времени как решения рассматриваемой системы). С другой стороны, второе уравнение (7.2) можно представить в виде
U = X(t)u + n(t)
(7.4)
Л = aw, ^ = g(a cos ф — sin ф),
(7.5)
уравнения движения частицы по прямой R = {£} в потенциальном поле и сопротивляющейся среде:
/£ - (mcu)i + Щ = 0, V = cp( 1 - cos О-
Так как c > 0, то потенциал V имеет строгий локальный минимум в точке £ = 0. Если с некоторого момента времени u(t) < 0, то равновесное состояние £ = 0, £ = 0 устойчиво в силу известного неравенства
А dt
т£2
~Y + V(0 = (mcu)Z2 < 0.
(7.6)
Если сани движутся по вертикальной прямой (£ = ^ — р = 0), то, очевидно, в некоторый момент времени скорость и будет отрицательной (сани будут двигаться вниз). По непрерывности это свойство имеет место при малых начальных возмущениях £ и £.
Согласно (7.5), функция ц отрицательна при малых значениях £. Таким образом, если £(Ь) мало, то (по лемме) и(Ь) < 0, начиная с некоторого момента Ь. Но если и(Ь) отрицательно, то значения £ (и £) остаются малыми согласно неравенству (7.6). Применяя индукцию по времени, отсюда выводим заключение теоремы.
По-видимому, при с > 0 имеет место асимптотическая устойчивость по отношению к £ и £. По крайней мере это верно для линеаризованного уравнения, которое имеет вид
£ + at£ + а£ = 0,
Pc
а = у > 0.
(7.7)
Для доказательства асимптотической устойчивости равновесия £ = £ = 0 линейной системы можно воспользоваться разложениями в ряды решений уравнения (7.7). Представим (7.7) в виде системы
х(АЬ + В)х, х = (£,£ )т,
(7.8)
A = 0 0 , в = 0 1
0 -а -а 0
Поскольку собственные значения матрицы А (числа 0 и —а) различны, то (как показано в [13]) решения системы (7.8) можно искать в форме следующих рядов:
x(t) = t* ( XQ + ^ + |f + ... ) exp(P? + 7i).
(7.9)
Здесь Х0,Х1,... — векторы из М2, 2@ — собственное число матрицы А, а константы к и 7 подлежат нахождению.
Так как а > 0, то собственному значению 2в = —а отвечает решение (7.9), которое сверхэкспоненциально быстро стремится к нулю. Нулевому собственному значению матрицы А (случай в = 0) отвечает решение системы (7.8) в виде ряда по обратным степеням Ь (к = —1, 7 = 0):
1
/ Х1 t \ жо + — + •••
xj е R2.
Оно также стремится к нулю при Ь ^ Таким образом, равновесие линейной систе-
мы (7.7), действительно, асимптотически устойчиво.
8. Вращение твердого тела с левоинвариантной сервосвязью в поле силы тяжести
Вращение твердого тела вокруг неподвижной точки с левоинвариантной линейной сервосвязью в осесимметричном силовом поле с потенциалом V (7) описывается следующей системой дифференциальных уравнений, обобщающих систему (5.1):
1и + и х 1и = 7 х V'1 + \Ъ, А/ + и х 7 = 0, (а, и) = 0. (8.1)
Здесь 7 — единичный вектор вертикали в подвижной системе отсчета, связанной с твердым телом. Если тело вращается в поле силы тяжести, то
V (7 ) = (с,7), (8.2)
где постоянный вектор с равен произведению веса тела на радиус-вектор его центра масс в подвижной системе отсчета. Без ущерба для общности можно считать, что левоинвари-антная связь имеет вид из = 0.
Предложение 2. Если выполнены условия (5.15) и (5.16), то фазовый поток системы (8.1) сохраняет стандартную меру в М6 = М3 {и} х М3{7}.
Доказательство состоит в непосредственной проверке следующего факта: дивергенция правой части системы дифференциальных уравнений (8.1) равна нулю, если выполнены условия (5.15)—(5.16). В отсутствие силового поля это установлено в §5. Оказывается, добавление потенциального силового поля не влияет на величину дивергенции.
В неголономной модели (когда Ъ2 = 0, а Ъз = 0) условия (5.15) и (5.16) эквивалентны равенствам 1\з = 12з = 0. Это в точности означает, что вектор а есть собственный вектор оператора инерции. Результат о сохранении меры в этом случае отмечен в [8].
Рассмотрим теперь некоторые аспекты задачи о вращении твердого тела в однородном силовом поле с потенциалом (8.2).
Теорема 7. Пусть 1а = ра (р £ М), а вектор с ортогонален векторам а и Ъ. Тогда уравнения (8.1) допускают линейный по скорости интеграл
Е = (1и, с). (8.3)
В неголономной модели векторы а и Ъ коллинеарны. В этом случае имеем одно условие (а, с) = 0; интеграл (8.3) получен Е.И.Харламовой в работе [21]. Подчеркнем, что если а — собственный вектор оператора инерции, то в общем случае (когда Ъ2 = 0) уравнения движения не допускают инвариантной меры с гладкой положительной плотностью. В условиях теоремы 7 неголономные уравнения интегрируются в квадратурах, поскольку кроме инвариантной меры и интеграла (8.3) уравнения движения допускают еще три первых интеграла
+ 7), (а, со), (7,7). (8.4)
Свойство интегрируемости вытекает из классической теоремы Эйлера-Якоби об интегрирующем множителе.
Для сервосвязей общего вида интегралами являются (8.3) и две последние функции (8.4), что явно недостаточно для интегрируемости уравнений вращения твердого тела.
Линейный интеграл (8.3) имеет вид нётерова интеграла:
Покажем, что его действительно можно вывести из обобщенной (по сравнению с формулировкой из введения) теоремы Нётер. Для этого введем векторное поле
w = с (8.5)
и покажем, что оно будет обобщенным полем симметрий. Оно представляет собой лево-инвариантное векторное поле на группе БО(3). Напомним, что поток любого правоинва-риантного векторного поля сохраняет левоинвариантную кинетическую энергию. Потоки левоинвариантных полей в общем случае таким свойством не обладают.
Поскольку (Ъ, с) = 0, то поле (8.5) есть поле возможных перемещений; оно удовлетворяет уравнению (Ъ, 5х) = 0. Далее, поток векторного поля (8.5) сохраняет потенциальную энергию. Действительно, поле w = с порождает следующую систему дифференциальных уравнений в М3 = {^}:
7' = 7 х с.
Следовательно, полная производная потенциала (8.2) в силу этой системы равна
ЩЛ) = (с, 7 х с) = 0.
Осталось показать, что поток векторного поля w, продолженный в фазовое пространство, сохраняет кинетическую энергию (на гиперплоскости связи). Полная производная от кинетической энергии (функции на алгебре вв(3)) вдоль продолженного поля w равна
Эта квадратичная форма относительно и, конечно, не тождественный нуль (если оператор инерции не пропорционален единичному оператору). Однако она будет равна нулю на плоскости (а, и) = 0 при выполнении условий теоремы 7 (в этом и заключается незначительное обобщение теоремы Нётер). Действительно, вектор и ортогонален вектору а (согласно уравнению связи). Далее,
(1и, а) = (и, 1а) = р(и, а) = 0.
Следовательно, вектор и х 1и будет коллинеарен вектору а. Но тогда правая часть (8.6) равна нулю, так как (с, а) = 0 по условию теоремы. Что и требовалось.
9. Некоторые задачи
В качестве заключения мы приведем перечень задач, которые возникают в связи с теорией движения систем с сервосвязями Бегена.
1°. Найти условия сохранения фазового объема приведенной системы (2.4) в случае когда группа О унимодулярна (то есть когда ^ скк =0 для всех 1 ^ г ^ п; в частности, все компактные группы Ли унимодулярны). Этот вопрос упирается в вычисление дивергенции однородного векторного поля. Для классических уравнений Эйлера-Пуанкаре эта дивергенция обращается в нуль, если группа О унимодулярна [8].
2°. Было бы интересным найти многомерные обобщения теорем 4 и 5. Вот один из вариантов: если система (3.1) имеет целую плоскость равновесий ^ aiXi = 0, то заменой времени
dr = ^^ aix^ dt
она приводится к линейной системе. Больший интерес представляют многомерные аналоги теоремы 5.
3°. Дать классификацию фазовых портретов динамических систем на плоскости с однородными квадратичными правыми частями (с учетом анализа ухода траекторий на бесконечность).
4°. В рамках этой классификации указать все типы фазовых портретов в задаче о вращении волчка с левоинвариантной сервосвязью. Первый шаг в этом направлении — найти условия, при которых соответствующая приведенная система на плоскости имеет либо ровно одну, либо три инвариантных прямых.
5°. Изучить движение сервосаней по горизонтальной плоскости в общем неуравновешенном случае, когда центр масс твердого тела не лежит на вертикальной прямой, проходящей через точку контакта.
6°. Рассмотреть движение сервосаней по наклонной плоскости и выяснить общие условия, при которых сани не спускаются вниз, а движутся в среднем по горизонтали.
7°. Исследовать интегрируемость уравнений (8.1) с потенциалом (8.2), описывающих вращение тяжелого твердого тела с сервосвязью Суслова.
Автор дружески благодарит А. В. Борисова, И. С. Мамаева за полезные обсуждения.
Список литературы
[1] Béghin M. H. Etude théorique des compas gyrostatiques Anschütz et Sperry. Paris: Impr. nationale, 1921. 132 pp. См. также: Беген А. Теория гироскопических компасов Аншютца и Сперри и общая теория систем с сервосвязями. Москва: Наука, 1967. 172 с.
[2] Appel P. Traite de Mecanique rationnelle: Vol. 2. Dynamique des systemes. Mecanique analytique. 6th ed. Paris: Gauthier-Villars, 1953. 584 pp. См. также: Аппель П. Теоретическая механика: Т. 2: Динамика системы. Аналитическая механика. Москва: Физматгиз, 1960. 487 с.
[3] Козлов В. В. Принципы динамики и сервосвязи // Вестн. Моск. ун-та. Сер. 1. Матем. Механ., 1989, №5, с. 59-66.
[4] Poincare H. Sur une forme nouvelle des equations de la Mecanique // C. R. Acad. Sci., 1901, vol. 132, pp.369-371.
[5] Kozlov V. V. General theory of vortices. (Encyclopaedia Math. Sci., vol.67.) Berlin: Springer, 2003. 184 pp. См. также: Козлов В. В. Общая теория вихрей. Москва-Ижевск: НИЦ «Регулярная и хаотическая динамика», Ижевский институт компьютерных исследований, 2013. 324 с.
[6] Козлов В. В., Фурта С. Д. Асимптотики решений сильно нелинейных систем дифференциальных уравнений. Москва-Ижевск: НИЦ «Регулярная и хаотическая динамика», Ижевский институт компьютерных исследований, 2009. 312 с.
[7] Kozlov V. V. Exchange of stabilities in the Euler-Poincare-Suslov systems under the change of the direction of motion // Nonlinear Dynamics & Mobile Robotics, 2014, vol. 2, no. 2, pp. 199-211.
[8] Kozlov V. V. On the integration theory of equations of nonholonomic mechanics // Regul. Chaotic Dyn., 2002, vol.7, no. 2, pp. 161-176.
[9] Андронов А. А., Леонтович Е. А., Гордон И. М., Майер А. Г. Теория бифуркаций динамических систем на плоскости. Москва: Наука, 1967. 487 с.
[10] Козлов В. В. Об инвариантных мерах уравнений Эйлера - Пуанкаре на алгебрах Ли // Функц. анализ и его прил., 1988, т. 22, №1, с. 69-70.
[11] Неймарк Ю.И., Фуфаев Н. А. Динамика неголономных систем. Москва: Наука, 1967. 520 с.
[12] Борисов А. В., Мамаев И. С. Динамика саней Чаплыгина // ПММ, 2009, т. 73, №2, с. 219-225.
[13] Coddington E. A., Levinson N. Theory of ordinary differential equations. New York: McGraw-Hill, 1955. 429 pp. См. также: Коддингтон Э. А., Левинсон Н. Теория обыкновенных дифференциальных уравнений. Москва: ИИЛ, 1958. 475 с.
[14] Киргетов В. И. О движении управляемых механических систем с условными связями (серво-связями) // ПММ, 1967, т. 31, №3, с. 433-446.
[15] Голубев Ю. Ф. Механические системы с сервосвязями // ПММ, 2001, т. 65, №2, с. 211-224.
[16] Румянцев В. В. О движении управляемых механических систем // ПММ, 1976, т. 40, №5, с. 771-781.
[17] Грдина Я. И. Заметки по динамике живых организмов. Екатеринослав: Екатеринославск. горн. ин-т, 1916.
[18] Blajer W., Seifried R., Kolodziejczyk K. Servo-Constraint realization for underactuated mechanical systems // Arch. Appl. Mech., 2015 (DOI 10.1007/s 00419-014-0959-2).
[19] Уткин В. И. Скользящие режимы и их применения в системах с переменной структурой. Москва: Наука, 1974. 272 с.
[20] Арнольд В. И., Козлов В. В., Нейштадт А. И. Математические аспекты классической и небесной механики. Москва: Едиториал УРСС, 2009. 416 с.
[21] Харламова-Забелина Е.И. Быстрое вращение твердого тела вокруг неподвижной точки при наличии неголономной связи // Вестн. Моск. ун-та. Сер. 1. Матем. Механ., 1957, №6, с. 25-34.
The dynamics of systems with servoconstraints. I
Valery V. Kozlov
Steklov Mathematical Institute, Russian Academy of Sciences
Gubkina st. 8, Moscow, 119991, Russia
kozlov@pran.ru
The paper discusses the dynamics of systems with Beghin's servoconstraints where the constraints are realized by means of controlled forces. Classical nonholonomic systems are an important particular case. Special attention is given to the study of motion on Lie groups with left-invariant kinetic energy and left-invariant constraints. The presence of symmetries allows one to reduce the dynamic equations to a closed system of differential equations with quadratic right-hand sides on a Lie algebra. Examples are given which include the rotation of a rigid body with a left-invariant servoconstraint — the projection of the angular velocity onto some direction fixed in the body is equal to zero (a generalization of the nonholonomic Suslov problem) — and the motion of the Chaplygin sleigh with servoconstraints of a certain type. The dynamics of systems with Beghin's servoconstraints is richer and more varied than the more usual dynamics of nonholonomic systems.
MSC 2010: 34D20, 70F25, 70Q05
Keywords: servoconstraints, symmetries, Lie groups, left-invariant constraints, systems with quadratic right-hand sides
Received February 10, 2015, accepted March 05, 2015
Citation: Rus. J. Nonlin. Dyn., 2015, vol. 11, no. 2, pp. 353-376 (Russian)
