Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2 УДК 517.925.51:517.93:531.36 МБС 93С10, 34Н15
Динамика твердого тела от уравнений Эйлера до управления угловым движением ИСЗ в трудах ученых СПбГУ. Ч. 2*
А. А. Тихонов
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Тихонов А. А. Динамика твердого тела от уравнений Эйлера до управления угловым движением ИСЗ в трудах ученых СПбГУ. Ч. 2 // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2. С. 259302. https://doi.org/10.21638/spbu01.2024.203
Работа является продолжением обзора, посвященного 300-летию Санкт-Петербургского государственного университета (СПбГУ) и представляющего собой попытку анализа научных достижений Санкт-Петербургской школы механики в области динамики твердого тела. Вторая часть обзора посвящена 50-летнему периоду, завершающемуся в 2023 г. Основное внимание в ней уделяется общетеоретическим исследованиям, выполненным учеными СПбГУ и посвященным как неуправляемому, так и управляемому движениям твердого тела.
Ключевые слова: твердое тело, космический аппарат, искусственный спутник Земли, динамика, вращательное движение, управление, стабилизация, колебания.
1. Введение. В первой части обзора [1] были представлены некоторые результаты научных исследований, выполненных учеными СПбГУ (далее — Университет) в области динамики твердого тела в первые 250 лет истории университета. Рассмотрены не только качественные и количественные результаты анализа движения тел, но и методы исследования, позволившие получить эти результаты.
Данная вторая часть обзора посвящена сравнительно короткому, но чрезвычайно богатому на научные результаты 50-летнему периоду истории Университета, завершающемуся в 2023 г. В этот период времени динамика твердого дела в Ленинградском — Санкт-Петербургском университете развивалась как базовый раздел механики, необходимый для моделирования актуальных прикладных задач. Поэтому в первую очередь в данной части обзора рассматриваются общетеоретические исследования, посвященные анализу движения твердого тела. Основное внимание концентрируется на вопросах динамики. При этом рассматриваются как неуправляемое, так и управляемые движения твердого тела. Вместе с тем систематическое изучение вопросов динамики твердого тела не могло идти в полном отрыве от изучения вопросов кинематики твердого тела. Поэтому попутно рассматриваются также и актуальные вопросы кинематики твердого тела.
* Первую часть статьи см.: Тихонов А. А. Динамика твердого тела от уравнений Эйлера до управления угловым движением ИСЗ в трудах ученых СПбГУ. Ч. 1 // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2023. Т. 10 (68). Вып. 3. С. 457—487. https://doi.org/10.21638/spbu01.2023.303
© Санкт-Петербургский государственный университет, 2024
2. Лагранжев и гамильтонов формализм в динамике твердого тела.
Канонические уравнения Гамильтона обладают определенными преимуществами при использовании методов теории возмущений [2]. Естественно предположить, что эти преимущества могут сыграть роль и в динамике твердого тела. В связи с этим в работе [2] В. С. Новоселов ставит вопрос о приведении дифференциальных уравнений механики к каноническому виду. Известно, что в случае действия потенциальных сил движение механической системы может быть описано каноническими уравнениями. С целью указать другие случаи приведения уравнений движения к каноническим автор предполагает, что имеются уравнения движения, разрешенные относительно обобщенных ускорений:
дг = Яг(Чз (г,з = 1, 2,...,в). (1)
Поскольку с математической точки зрения представление уравнений в гамиль-тоновом виде равносильно записи этих уравнений в форме уравнений Лагранжа
с1 дЬ дЬ Л дщ дцг '
то для существования решения поставленной задачи необходимо и достаточно выполнения равенств
Я Л ди +ди (3)
И дцг дцг
в
на траекториях системы (1). Функция Лагранжа будет равна Ь = ^ ^ д? + £/. Для
г=1
тождественного выполнения равенств (3) во всем фазовом пространстве имеются необходимые и достаточные условия Гельмгольца:
дЯг 1 дЯэ п дЯг дЯэ 1 (I (дЯг дЯ
дщ дщ ' 2 Л \ дщ дщ ) ' ^^ ' ' ' ^ ^ ^
Функция и, удовлетворяющая равенствам (4), может быть найдена с помощью известной процедуры, изложенной в [3], и в усовершенствованном виде — в [4]. Если же потребовать выполнения равенств (3) не во всем фазовом пространстве, а лишь на траекториях системы (1), то допустимо преобразование правых частей уравнений (1) в виде
в
Яг = Яг + - Я»), (5)
м=1
где Ргм — подлежащие определению функции времени, обобщенных координат и обобщенных скоростей. Далее автор выписывает необходимые и достаточные условия Гельмгольца для представления функции Яг в виде (3) во всем фазовом пространстве:
дЯг дЯз „ дЯг , дЯз _ 0 1 дЯг дЯг дЯз _ 1 1 ( дЯг дЯз \
и' + ---"Е---Уь>
дс[з дцг ' дц^ дцг Л дц^ дцз дцг 2 И \ дцз дд.
и находит следующие условия, которым должны удовлетворять неопределенные множители Ргц для приведения уравнений движения к гамильтонову виду:
Ргз = Рзг, (7)
dPjj -A dQM _ / dPjj -A dQ
di . . - I df + ад ^
dQi , 1 A dQi dQM 1 d dQ
Р Р -Р Р Р - -и ± _
¿1 - ^ - дд. +4^ 2<Й • 1 ;
Несложно проверить, что если выполнены условия (4), то уравнения (7)—(9) допускают нулевое решение для множителей Р^. Далее рассмотрены следующие частные случаи приведения уравнений к каноническому виду: Qi = Qi(t), Qi = Qi(qj), Qi = Qг(Чj). Для каждого из них приведены примеры, включая такие, в которых условия (4) не выполняются.
В работе [5] В.С.Новоселов рассматривает метод минимизации функционала действия для нахождения характеристик нелинейных колебаний. Ранее такой подход предлагался в [6] и применялся в [7] для решения уравнения Дуффинга:
х + х + /х3 = 0 (10)
с малым параметром Для дифференциального уравнения более общего вида
х + ¡3^)х + а(х^)=0 (11)
в [6] предложено заменить его решение исследованием минимума функционала:
42
ti \о /
dQ dR ( Г п. . 1
x(ti) = xi, хШ) = Х2,---= a(x,t)exp / р(т)ат
dt dt \ J
dQ dR dt dt
для которого (11) является уравнением Эйлера — Лагранжа. В.С.Новоселов в работе [5] устраняет ошибку, а также дает уточнение и развитие результатов работы [7]. Идея работы [7] заключается в использовании интеграла механической энергии на траекториях уравнения (10) и в обращении в ноль множителя Q при t в вековом члене выражения I на аппроксимирующем колебании x = xi cos wt. Однако в действительности минимизация функции действия требует не выполнения равенства Q = 0, а минимизации Q по w. Поэтому искомые значения w получаются иными. Показано, что если удерживать члены порядка л в выражении для частоты w, то и в аппроксимирующем выражении для x также следует удерживать члены порядка ¡ . Сделан вывод о том, что метод минимизации функционала действия может практически применяться для определения характеристик нелинейных колебаний. Но нужно учитывать дополнительные трудности, которые возникают вследствие того, что членам порядка л в фазовых координатах могут отвечать члены порядка ¡л2 в минимизируемом функционале.
Таким образом, проблема сведения задачи исследования решения дифференциального уравнения (11) к изучению функционала (12) равносильна замене уравнения (11) каким-либо эквивалентным уравнением в форме Лагранжа или к представлению уравнения (11) в каноническом виде. Отмечается, что приведение уравнения
(11) к каноническому виду может быть выполнено непосредственно и рассматривается пример. Далее в [5] обсуждается задача приведения к каноническому виду для диссипативной или возбуждаемой механической системы с п степенями свободы. Если с помощью неособенного преобразования привести уравнения колебательного движения такой системы к виду
Щ + 2ЪгХг + • = г = 1, П,
дхг
то, вводя канонические переменные по формулам цг = хгвь^, рг = Хгвь^, можно представить эту систему в виде
ди 2Ь.г
Рг = -ОгРг - — е г , Чг=Рг + Мг-ддг
Далее несложно проверить, что Н = {^¡Рг + ЬчгРг + е26^ о представ-
ляет собой функцию Гамильтона.
Одним из важных и фундаментальных вопросов кинематики твердого тела, а вслед за ней и динамики твердого тела, является вопрос о нахождении наилучшего описания ориентации твердого тела. Удачный выбор координат для задания ориентации твердого тела позволяет существенно упростить математическую модель и избежать таких неприятностей, как сингулярности и сложные нелинейности кинематических дифференциальных уравнений. И наоборот, плохой выбор координат может искусственным образом ограничить возможности системы управления движением твердого тела путем требования функционирования в области, не содержащей сингулярностей. Поэтому неудивительно, что вопрос о выборе наилучших переменных (иногда называемых параметрами ориентации) для описания ориентации и углового движения твердого тела изучался такими известными учеными, как Эйлер, Якоби, Гамильтон, Кэйли, Клейн, Родриг, Гиббс, Хопф и привел к появлению многих выдающихся результатов. Тем не менее вопрос о возможности выбора таких обобщенных координат, которые были бы минимальны по числу (соответствовали бы трем степеням свободы твердого тела) и не имели бы точек вырождения, оставался открытым на протяжении 200 лет после появления первых работ Эйлера по кинематике твердого тела. Лишь в 1964 г. был дан отрицательный ответ на поставленный вопрос: невозможно «улучшить» известные методы параметризации вращений твердого тела без введения избыточных переменных [8].
Отчасти подытоживая 200-летние изыскания предыдущих поколений ученых, отчасти приближая изложение методов теоретической механики к изложению математических дисциплин в Университете, С. Н. Кирпичников и В. С. Новоселов опубликовали книгу [9], которая оказалась весьма своевременной из-за отсутствия в литературе достаточно полной математической формализации методов кинематики твердого тела, а также в связи с тем, что математическое описание кинематических методов стало явно уступать изложению методов динамики. В этой книге теорию перемещений твердого тела, теорию его движения и кинематику сложного движения предлагается изучать с помощью единообразных групповых, операторно-матричных и топологических методов. Большое внимание уделяется определению и изучению структуры пространств конфигураций твердого тела, а также вопросам задания ориентации тела, как локальной, так и глобальной. Дано достаточно полное математически формализованное изложение кинематики твердого тела
и кинематики сложного движения, в котором естественным образом большое внимание уделено вопросам топологии. Примененная в [9] операторная запись формул позволила эффективно использовать качественные методы исследования (в отличие от координатной записи, которая часто ввиду своей громоздкости маскирует математическую сущность записываемых соотношений). При этом последовательное применение операторных и топологических методов в кинематике твердого тела оказалось весьма плодотворным и позволило помимо единообразия изложения, математической ясности и большей строгости доказательств получить и некоторые новые результаты. Например, в [9] впервые показано, что введение угловой скорости твердого тела соответствует установлению структуры тривиального расслоения в фазовом пространстве твердого тела.
В работах Н. Н. Поляхова, С. А. Зегжды, М. П. Юшкова (см., например, [10]) показано, что математический аппарат теоретической механики, разработанный применительно к системе материальных точек, может быть применен и к механическим системам, состоящим из твердых тел. Это означает, что при исследовании динамики системы твердых тел могут быть использованы уравнения Лагранжа как I, так и II рода. При этом задача отыскания простейших по форме уравнений Лагранжа I рода для системы твердых тел сводится к отысканию новой формы записи уравнений движения одного твердого тела. Одна из таких форм уравнений динамики твердого тела предложена в [11]. Положение тела с массой т в пространстве можно задать с помощью избыточных координат, в качестве которых выбираются радиус-вектор р центра масс тела и тройка ортов г, р, к, направленных по главным центральным осям инерции тела Схуг. В силу абсолютной твердости тела выполняются условия ортонормированности /1 = г — 1 = 0, /2 = ] — 1 = 0, /3 = к — 1 = 0, /4 = г • р = 0, 15 = Р •к = 0, /6 = к • г = 0, рассматриваемые как уравнения голономных связей1. С введением обозначений 1Х = х21т, 1у = у21т, 1г = г2вт удвоенная кинетическая энергия тела примет вид
• 2 • 2 • 2 г 2 2Т = тр + 1x1 + 1у• + . (13)
Пусть к телу в точках Ми (хи ,у1/, ги) приложены силы 1Ри. Виртуальная работа этих сил записывается в виде
5А = ^2• (5 р + хр5г + у„5р + г„5к) = Ср • 5 р + СРр • 5г + Я] • 5] + ^ • 5к,
V
где Ср IV, Ср =]Т XV, С] = yvIV, С% = ^ гv¥„. (14)
V V V V
Тогда движение тела можно описать четырьмя векторными уравнениями Лагранжа I рода с неопределенными множителями Лагранжа А^, г = 1,6:
-—" —(15)
й дТ дТ Л др др
<1 дТ дТ Р А А
-— - —; = д?+ V = д?+ 2А1г + Л4;+ Авк, (16)
11 дг дг К=1 дг
1 Скалярное произведение векторов здесь обозначается знаком •. В дальнейшем этот знак опус-
кается. Векторное произведение векторов обозначается X.
Л дТ дТ А , д]к ^ ^ ^
= + + (17)
Л д] д3 3 ^ д]
Л дТ дТ А ,
-— - Т7 = ^ + Е = У* + + Лб* + (18)
^ дк дк дк
где приняты обозначения
д _ д ^ д -> д _ ^ ^ 77"? = я /1 г + а л 3 + ял'5' ^ = + -^-уЗ + Агк.
дА дЛх дАу дАг С учетом этих обозначений и равенств (14) уравнения (15)—(18) запишутся в виде
ш/ = ^ , (19)
V
У = Е XV + 2Л° +А4о +Лбк, (20)
V
1уЗ = Е yv ^ +2Л2/ +Л5к + Л4°, (21)
V
/у/ = Е ZV ^ + 2Лзк + Л° +Л5о. (22)
V
Двукратное дифференцирование по времени уравнений связей и исключение множителей Лагранжа из уравнений (20)—(22) дает:
Г = (23)
1х + 1у 1 + 1х 1х + 1у 1 + 1х
Ч- ГЛ2^ 2-Гг (о ■ к) ^ 21х(г ■ р) Ьх ^ Ьх
3 = 3~ т I г г I г г+Г1ГА:~Г1Г г (24)
Iу + /г Iх + Iу Iу + ^ Iх + 1хУ
к = -{к к-^—г1~^ТГг + ТТТг~ТТТ]- (25)
Здесь Ъх, Ьу, Ъг — проекции главного момента Ъ = ^V(х] + yVj + zVk) х FV внешних сил относительно центра масс. В системе уравнений (19), (23)—(25) уравнение (19) соответствует движению центра масс, а уравнения (23)—(25) — вращению тела относительно центра масс.
Во многих задачах небесной механики и космодинамики большой интерес представляет анализ углового движения твердого тела (естественного или искусственного происхождения) в окрестности устойчивого положения равновесия, обусловленного влиянием гравитационного момента. Пусть центр масс тела движется по круговой кеплеровой орбите в центральном гравитационном поле, тензор инерции тела Л = diag(A, В, С) в главных центральных осях инерции Cxyz с ортами г,р,к удовлетворяет неравенствам В > А > С. Тогда устойчивым будет такое положение тела, в котором орты г,р, к совпадают соответственно с ортами Оо,Оо, Оо орбитальной системы координат, вращающейся с угловой скоростью в инерциальном пространстве, причем ось £ направлена по касательной к орбите в сторону движения, ось ц —
по нормали к плоскости орбиты, а ось Z — вдоль радиус-вектора центра масс тела относительно притягивающего центра [12]. Если задать ориентацию тела в орбитальной системе координат с помощью «самолетных» углов ф, р, в, то указанному устойчивому положению будут соответствовать значения ф = р = в = 0. Такое положение тела часто называют положением гравитационной ориентации, а тело, совершающее либрационное движение в окрестности этого положения, — гравитационно-ориентированным телом. Либрационное движение тела происходит, вообще говоря, не только под действием гравитационного момента [12] Mg = 3w2Mo х (JM0), но и под действием других моментов различной природы (обозначим их главный момент Mb ), вызывающих возмущения в колебаниях гравитационно-ориентированного тела. Исследование возмущенных колебаний является нелинейной задачей, описываемой динамическими уравнениями Эйлера
— Ul5)+L5XUL5)=Mg + Mb, (26)
dt
и для приближенного аналитического решения этой задачи, с учетом условия \Mb | ^ |mg|, естественно использовать методы возмущений. Однако методы возмущений обычно связаны с переходом к новым переменным, и часто бывает непросто построить преобразованную систему дифференциальных уравнений движения. В некоторых случаях эти трудности можно в значительной степени преодолеть, используя метод вариации канонических произвольных постоянных. Этот метод удобен для построения дифференциальных уравнений возмущенного движения механических систем, состояние которых можно описать каноническими уравнениями. Например, в работах [12, 13] метод вариации канонических произвольных постоянных в сочетании с методом усреднения успешно применялся при Mb = 0 для исследования нелинейных колебаний гравитационно-ориентированного спутника под действием гравитационного момента, взятого в квадратичном приближении.
В Университете метод канонической вариации произвольных постоянных получил развитие и был обобщен на механические системы с неканоническими возмущениями. В работе [14] автором построены дифференциальные уравнения возмущенного либрационного движения гравитационно-ориентированного тела, обобщающие известные уравнения, использованные в [12] и [13], на случай, когда наряду с потенциальными возмущающими моментами действуют непотенциальные возмущающие моменты, взятые в квадратичном приближении.
Для получения этих уравнений вводятся в рассмотрение канонические переменные Ф,ф,в, Р1,Р2,РЗ, где pi = Р2 = Щ, Рз = §§ - Blvq, L = + Вш2у + Сш\) — — СУ/1 + (В — Cyfá] — функция Лагранжа, в которой шх, шу, wz, 71,72 выражены через ф, р, в, ф, р, в на основании кинематических формул wx = р — ф sin в + w0 sin ф cos в, uy = в cos р + ф sin р cos в + w0 (cos р cos ф + sin р sin ф sin в), = ф cos р cos в — в sin р+Wi^cos р sin ф sin в— sin р cos ф), 71 = — sin в, 72 = sin р cos в. После введения в рассмотрение функции Гамильтона H = Р1ф + Р2р + Рзв — L дифференциальные уравнения (26) могут быть записаны в виде
дН
Pl = — -7— + Мву sin Ср COS в + Mbz COS Ср COS в — Мвх sin дф
дН дН
Р2 = --Q^ + МВх, РЗ = --QQ + МВу COS <р - MBz sin <p, (27)
ф-дн ■ - дн в - дн
dpi' ^ др2' дрз'
Далее вводится новая безразмерная независимая переменная u = шоЬ и новые безразмерные обобщенные импульсы:
Pi = (Лшо)-1Р2 - ф, P2 = (Buo)-1ps, P3 = (Cuo)-1pi + (1 - EC-1 )у. (28)
Идея такой замены заключается в том, что формулы (28) совпадают с линеаризованными и деленными на шо выражениями для проекций относительной угловой скорости тела и могут рассматриваться как формулы перехода от обобщенных импульсов Р1,Р2,Рз к новым обобщенным импульсам Pi,P2,Рз того же порядка малости. В новых переменных функция Гамильтона принимает вид H(ф, у, 6, P1,P2,P3) = Но + Hl + H2, где Но содержит члены 2-го порядка малости, а Hi — члены 3-го порядка малости:
1 13
Н0 =-(АР2 + BP2 + С Pi) + 2 (В - СV2 + - (В- А)-ф2 + С)в\
B
Н1=(В-С)^фв+(В-С)Р2^2 + (А-С)Рз0ф + -Р2ф2 + (С-В)Р2Рз^+АР1Рз0.
В системе (27), записанной для Н = Но + Hi, осуществляется переход к новым импульсам по формулам (28). Путем разрешения системы (27) относительно производных Pi, P2, P3, ф', у', 6, взятых по безразмерной переменной u, образуется система
P1 = Fi(Pi,P2,P3, ф, 6, у), у' = Gi(Pi,P2,P3, ф, 6, у),
P2 = F2P1 ,P2 ^3,ф,6,у), 6' = G2 (Pi, P2, P3, ф, 6, у), (29)
P3 = F3 (Pi ,P2 ,P3,ф,6,у), ф' = G3 (Pi ,P2 ,P3 ,ф,6,у),
учитывающая члены 2-го порядка малости включительно. Дальнейшее ее преобразование опирается на переход к новым переменным, связанным с рассмотрением уравнений (29) в линейном приближении, соответствующем гамильтониану Ho. При отсутствии возмущающего момента Miß общее решение уравнений линейного приближения имеет вид [12]
Р\ = ~ ^к^гл/зйвтп, Р2 = л/2x3В-1 соетз, Р3 = ^ Ьацл/зйсоей, (30)
г=1 г=1
2 2 (р = '^2(T2i^/x¡cosтi, в = -—у/2хзВ-1вттз, ф = ^ацу^втп. (31)
- -I к3 -
1=1 1=1
Здесь т = (г = 1, 2, 3); х^, VI — канонические произвольные постоянные; к —
невозмущенные собственные частоты колебаний гравитационно ориентированного тела:
2 = Ja/2TVа2/4—ï>, к3 = ^3(1 - e)/ö, ^и = \Hr, a2i = ~ Т\ГТ '
ki V Ai ki V Ai
e = C/A, S = В/A, a = 1 + 3(6 - e) + b/4, b = 4(6 - l)(S - e)/e,
_ (-1 yBe(ki-kï) _ h(l-ô + e) 6(6-l-ek2) ' ô-l-ek2'
С учетом членов 2-го порядка малости решение системы (29) отыскивается в виде (31), где хг(и), v¿(u) — некоторые неизвестные функции и, подлежащие определению. Формулы (30) и (31) рассматриваются как формулы перехода к новым переменным хг(и), Vi(u). После подстановки этих формул, а также результата их дифференцирования по и, в систему (29) и разрешения этой системы относительно производных хг' и VI (г = 1, 2, 3), получена система дифференциальных уравнений вида
Хг' = хн'(хн^г), VI = Vi(xi,Vi), (г = 1,2,3), (32)
построенных с учетом членов 2-го порядка малости. Аналогичным образом преобразуются слагаемые в функции Н. В частности, для Н получено выражение
#1 = жхл/жз(/г.]1'1 совтз + ссе(2т 1 - т3) + соз(2т1 + т3)) +
+ Х2\/хз(1г^ соэтз + сов(2т2 - т3) + Ь^ соз(2т2 + т3)) +
+ л/ж!Ж2Жз (Н^ соэ(т1 - т2 - т3)+ + к^ ео8(т1 - Т2 + тз) + 41} ео8(т1 + Т2 - тз) + кЦ еов(т1 + Т2 + тз)),
коэффициенты которого, зависящие только от инерционных параметров тела, приведены в [14]. В результате система (32) имеет вид
V ai ао /
dvi у ' у \ аi Bx а2
/ . (-1)* а1,3-гТГ(3) . a2,3-i7J(3) . \ /оо\
** = ^ + M-COSTi + —M-sm4 (33)
*з' = -Ц + V^ Mg cosr3, ,3' = Ц - ^Mg sinr3,
где Mb* = Mßx(A^02)-1; MBy = Мву(Вш1)-1 ^Мв = Mbz(Cu20)-1.
При отсутствии возмущающего момента Мв дифференциальные уравнения (33) являются каноническими и совпадают с известными уравнениями в канонических вариациях [12]. Уравнения (33) имеют вид, позволяющий непосредственно, без громоздких выкладок, применять для их исследования асимптотические методы нелинейной механики. Поэтому они оказалась удобными для исследования и позволили получить новые результаты. Далее этим результатам будет посвящен специальный раздел обзора. Здесь лишь упомянем первые по времени появления две статьи [15] и [16], в которых на базе уравнений (33) исследованы нелинейные колебания гравитационно-ориентированного заряженного спутника в суперпозиции гравитационного и магнитного (в дипольном приближении) полей Земли. В области гравитационной ориентации были обнаружены и исследованы 22 параметрических и четыре внутренних нелинейных резонанса, обусловленных квадратичными нелинейными членами.
В дальнейшем на основе уравнений (33), полученных в [14], был развит общий метод исследования либрационного движения твердого тела под действием возмущающего момента общего вида в классе нелинейных квадратичных функций относительно малых углов ориентации тела и их производных по времени. В [17] была выяснена спектральная структура этого возмущающего момента, что позволило обнаружить возможные параметрические и внутренние резонансы. В [18] найдены 63 резонансные комбинации невозмущенных собственных частот тела, реализу-
емых в области гравитационной ориентации, и построены соответствующие резонансные кривые в плоскости инерционных параметров тела. Также было выявлено 520 кратных резонансов в области гравитационной ориентации тела и найдены соответствующие значения инерционных параметров. В работах [19] и [20] система дифференциальных уравнений возмущенного движения тела усреднялась по явно входящему безразмерному времени в нерезонансном случае и в условиях выявленных резонансов. Выяснена структура полученных усредненных дифференциальных систем, получены коэффициенты усредненных уравнений в зависимости от параметров задачи и выведена унифицированная форма записи усредненных дифференциальных уравнений возмущенного движения тела для всех резонансных случаев. Исследованы колебания тела в нерезонансном случае и в условиях всех возможных резонансов. В рамках принятой квадратичной аппроксимации получены условия затухающих, ограниченных или нарастающих колебаний тела. Эти результаты позволяют исследователям конкретных прикладных задач динамики гравитационно-ориентированных спутников выяснить, какие резонансы могут возникать при различных возмущениях в классе квадратичных функций и какие динамические эффекты могут вызвать эти резонансы.
В то же время многие примеры известных задач нелинейной механики демонстрируют важность учета нелинейных членов 3-го порядка малости. В полной мере это справедливо применительно к задачам нелинейных колебаний в условиях внутренних резонансов. В работе [21] выведены такие дифференциальные уравнения возмущенного либрационного движения гравитационно-ориентированного тела, которые позволяют учесть влияние не только квадратичных, но и кубических нелинейных членов:
х3
дуг V а1 Вх Ст2
, « А(72 З—г ,- ТГ(2)
(-1) _——^х1хз М^^совпвттз -
2{-1)гВ(7^ ' л/ж7(а21л/д:1со8г1 + 0"22л/^2 соя совп,
д(#1 + #2) , (-1)г (а1,3-г ~(3) а2,3-г^{3) . \ ,
-я--I" -—— —-м^ сое п н--:-МУВ> 81П п +
дхн к1у/х1\ СГ1 Вх <72 Вг )
пг, ^ А(72 з—г^/хз —(2) . . , .
+ у2( —1) _ ; Мув;8ШП8ттз+ (34)
+ ( -1)' {а 21 л/жТ СОЭ п + 022\[Х2 СОвТ2)М^ ЯПТ;,
-—--Ь \J2BX3 МКВу соэтз -
/— 3 [?21 со8 Т1 + °"22 ф>2 СОЭ Г2 ) МВ1 СОИ Т3,
V В
Э(Я1+Я2) /В 77(3) . ^
+ ^ ^2Вх (СГ21л/Д71со8г1 +СГ22л/^2со8Т2)М^;;8тТз.
268 Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2
х
3
V
Здесь г = 1, 2; А, В, С — главные центральные моменты инерции тела; гамильтониан Н1 содержит члены 3-го порядка малости, Н2 содержит члены 4-го порядка малости; Ы{2, мВУ, мВ1 — выражения для безразмерных проекций возмущающего момента Мвх, Мву, Мвг соответственно, содержащие члены до второго порядка малости по л/жг", у/хз, а М^, М^у, М^ — аналогичные выражения, содержащие члены до 3-го порядка малости включительно. Легко видеть, что дифференциальные уравнения (33) вытекают из уравнений (34), если в последних нелинейные члены взять в квадратичной аппроксимации.
Уравнения (34) удобны как для выявления резонансных комбинаций 2, 3 и 4-го порядков (именно эти три набора являются наиболее важными в подавляющем большинстве случаев [22]), так и для дальнейшего аналитического исследования резонансных колебаний тела с помощью асимптотических методов. Так, для случая, когда возмущающий момент Мв отсутствует, получены следующие резонансные комбинации 4-го порядка для собственных частот гравитационно-ориентированного тела:
2к1 - 2к3 = 0, 3к1 - к2 = 0, 2к2 - 2к3 = 0, к1 + к2 - 2к3 = 0, к1 - к2 + 2к3 = 0. (35)
Все они могут реализоваться в области, где выполняются достаточные условия устойчивости В > А > С [12] гравитационной ориентации тела. В то же время выполнение достаточных условий устойчивости не гарантирует пребывания амплитуд колебаний тела (в частности, ИСЗ) в допустимом с практической точки зрения диапазоне. Поэтому особенно актуальным в задачах такого рода является анализ резонансных эффектов и вопрос об амплитудах колебаний тел в условиях возможных резонансов.
В работе [23] с использованием метода усреднения проанализированы дифференциальные уравнения возмущенных колебаний тела в нерезонансном случае и для всех резонансов гравитационного происхождения, включая резонансы (35). Построены полные наборы первых интегралов для усредненных дифференциальных уравнений в каждом из резонансных случаев. С помощью этих интегралов дана интерпретация решений соответствующих дифференциальных систем в амплитудно-фазовых плоскостях. Выявлено, что амплитуды резонансных колебаний ограничены. Этот вывод проверен численно путем компьютерных экспериментов, в которых на основе исходной нелинейной дифференциальной системы найдены верхние границы амплитуд колебаний тела для всех резонансных случаев.
3. От параметризации вращательного движения твердого тела к новым динамическим уравнениям. Уравнения динамики твердого тела являются математической базой для решения большого разнообразия прикладных задач из разных областей науки и техники. В зависимости от специфики этих задач, порой ярко проявляющейся в характере движения тел и целях выполняемых исследований, вначале применяются те или иные приближенные методы аналитических исследований, а затем, как правило, — методы компьютерного моделирования. При этом успех предпринимаемых аналитических исследований во многом зависит от того, насколько удачно выбраны переменные, используемые при построении математической модели задачи.
Методы компьютерного моделирования также чувствительны к выбору переменных, используемых при построении математической модели задачи. Например,
известны трудности использования традиционных углов Эйлера или различных вариантов «самолетных» углов для компьютерного моделирования вследствие того, что правые части дифференциальных систем в нормальной форме оказываются сложными трансцендентными функциями углов и сингулярными в области определения переменных, задающих ориентацию тела. Схема Эйлера — Пуассона не связана с сингулярностями в кинематических уравнениях, однако также является неудобной для компьютерного моделирования вследствие слишком высокого порядка соответствующих дифференциальных систем.
Одним из широко распространенных вариантов преодоления указанных трудностей является использование четырех параметров Родрига — Гамильтона Ао,А1,А2,Аз, представляющих собой компоненты кватерниона и подчиненных одному условию нормировки А2 + А'2 + А2 + Х^ = 1. Эти параметры не имеют точек вырождения и характеризуют наиболее естественный способ задания положения твердого тела с помощью плоского поворота. В силу своей наглядности и удобства использования как в аналитических исследованиях, так и в численных, параметры Родрига — Гамильтона широко применяются в прикладной теории гироскопов [24], в динамике искусственного спутника Земли (ИСЗ) [25] и в других областях механики [26].
При этом использование параметров Родрига — Гамильтона в динамике ИСЗ может стать еще более эффективным, если, учитывая специфику вращательного движения ИСЗ, опираться не на обычные уравнения Эйлера, замыкаемые кинематическими уравнениями в параметрах Родрига — Гамильтона, а на специально выведенные дифференциальные уравнения, учитывающие ротационный характер движения ИСЗ, т. е. наличие сравнительно большой угловой скорости ИСЗ, при которой кинетическая энергия вращательного движения ИСЗ намного больше работы главного момента M внешних возмущающих сил, действующих на ИСЗ. Такие дифференциальные уравнения получены в работе [27]. При выводе этих уравнений используются две системы координат с началом в центре масс тела: кенигова система координат CXYZ и система координат CL1L2L3, ось L3 которой направлена вдоль кинетического момента L тела. Ориентация тела или, что то же, системы его главных центральных осей инерции Cxyz, относительно системы координат CL1L2L3 задается с помощью параметров Родрига — Гамильтона А0 ,А1, А2 ,А3. Тогда матрица А = (Aij) (г = 1,3, j = 1,3) направляющих косинусов углов между этими осями имеет вид [26]
(А2 + А2 - А2 - АЗ3 2(А1 А2 - А0А3) 2АА3 + А0А2) \ A = I 2(А1 А2 + А0А3) А2 + А2 - А22 - А3 2(А2Аз - А0А1) | . V 2(А1 АЗ - А0А2) 2(А0А1 + А2А3) А0 + А2 - А2 - А2)
Вводится абсолютная угловая скорость тела ш = ше + шг, где ше — угловая скорость системы координат CL1L2L3 относительно кениговой системы CXYZ, а шг — угловая скорость тела относительно системы координат CL1L2L3. Тогда теорема об изменении кинетического момента приводит к дифференциальным уравнениям:
L = Мз, Ше2 = Mi/L, Ше1 = -M2/L, (36)
где ше1, ше2, шез — проекции угловой скорости ше системы координат CL1L2L3 относительно кениговой системы координат CXYZ на оси Li,L2,L3 соответственно, а Mi, М2, М3 — проекции момента М на оси Li, L2, L3 соответственно.
Кинематические уравнения в параметрах Родрига — Гамильтона
1\ Т (А31Л1 А32Л2 А33Л3 \ , , , , ,
2Л0 = -Ь I —----1--—--1--I + {ШехМ + ШеуМ +
1\ Т (А31Л0 А32Л3 А33Л2 \ , , , л \
¿М — ь I ——---—--1--—— I — (шехло — шеулз + шехЛ2),
1\ Г ( А31Л3 А32Л0 А33ЛЛ , , , л \
¿Л2 — Ь I ——--1--—---—— I — (ШежЛз + и)еуЛ§ — ШегЛ1),
о \ т ( А31Л2 А32Л1 А33ЛоА , , , л \
2Л3 = Ь I--—--1--—--1--I - (^е2Л0 + - шех\2)
совместно с динамическими уравнениями (36) образуют замкнутую дифференциальную систему. Эта система может принимать тот или иной конкретный вид в зависимости от того, какие переменные использовать для задания ориентации системы координат СЬ^2Ь3 относительно кениговой системы координат CXYZ. Если, следуя В. В. Белецкому [12], использовать для этой цели угловые переменные р и а, то из (36) получаются известные динамические уравнения:
Ь = М3, р = М1/Ь, а = Ы2/(Ь этр), (38)
относительно «медленных» переменных Ь, р, а, зависящие, однако, в общем случае, также и от «быстрых» переменных, определяющих тем или иным образом ориентацию главных центральных осей инерции тела Схуг относительно системы координат СЬ1Ь2Ь3. Если в качестве «быстрых» переменных использовать углы Эйлера, то из (37) получаются уравнения
. 1 эш2 ф соэ2 ф \ М1 соэ ф + М2 эт ф ^ = ---—)+-—--,
г .'эт2 ф соэ2 ф \ М1 соэф + М2 этф М2 ■Ф = Ь\ —— + —д— )--1- ^ ~ — Р,
(39)
11
= Ь I —■ — — ) 8И1 'в эш ф соэ ф + А В
М2 соэ ф — М1 эт ф Ь '
образующие вместе с (38) известные уравнения Белецкого — Черноусько [28], удобные тем, что в условиях ротационного движения, т. е. при выполнении неравенства \Ь|/А0 ^ \М\/\Ь|, где А0 — характерный момент инерции тела, в дифференциальной системе (38), (39) переменные Ь,р,а изменяются медленно, а ф,ф,д — быстро, что позволяет применять к этой системе приближенные асимптотические методы исследования движения.
Очевидные недостатки уравнений (39), заключающиеся в наличии трансцендентных функций угловых переменных ф,ф,д и сингулярности при д = 0, не возникают, если пользоваться уравнениями (37) непосредственно в переменных Родри-
га — Гамильтона:
'А31А1 , А32А2 , А33АЛ М1л М2
2А0 = -Ь + + + -^А2 - -^(А! - Аз <АёР),
т(А31А0 А32А3 А33А2 \ М, М2 , . , .
I ~4---в— ) ~ Т Т^
( ' (40)
т(Аз1Аз А32А0 А33АЛ М1 М2,л , . , .
= ^ + -в---—)-—х°- —(Лз + Л1
г ( А31А2 А32А1 А33А0N МЧ М2.л 2Аз = Ь (--— + —— + I + —А1 - — (А2 -\octgp).
В уравнениях (40) остается еще сингулярность при р = 0, но и эта неприятность может быть преодолена. Для этого достаточно вместо модуля кинетического момента и двух углов, определяющих направление вектора Ъ, ввести другие три медленно изменяющиеся переменные — проекции Ъх, Ъу, вектора Ъ на оси кениговой системы координат.
Действительно, если определить ориентацию тела, т. е. системы координат Схуг, относительно системы координат CXYZ с помощью матрицы А = А(Л0, А1; А2, Ад) направляющих косинусов, выразив их, в свою очередь, через параметры Родрига — Гамильтона А® (г = 0,3), то кинематические уравнения в параметрах А4
2А0 = -Ах (АпЬх + А21Ьу + АзгЬг) /А - А2 {А12ЬХ + А22Ьу + А32Ьг) /В— Аз (АгзЬх + А2зЬу + АззЬг) /С, 2АХ = А0 {АпЬх + А21ЬУ + А31Ьг) ¡А- Ад (А12ЬХ + А22Ьу + А32Ьг) /В+
+ Л2 (А13ЬХ + А23Ьу + АззЬг) /С, 2А2 = Аз (АпЬх + А21Ьу + Аз1Ьг) /А + А0 {А12ЬХ + А22Ьу + А32Ьг) /В- А1 (А13ЬХ + А23Ьу + АззЬг) /С, 2Аз = -А2 {АцЬХ + А21Ьу + А31Ьг) ¡А + Ах (А12ЬХ + А22Ьу + А32Ьг) /В+ + А0 (А13ЬХ + А23Ьу + АззЬг) /С
оказываются свободными от каких-либо сингулярностей или вырождений. Остается записать динамические уравнения в виде
(.Ьх Ьу Ьг)т = А (Мх Му Мг)т, (42)
чтобы получить замкнутую дифференциальную систему, удобную для описания ротационного движения тела, не имеющую особенностей и не содержащую тригонометрических функций [27]. Система (41), (42) удобна и для численного моделирования на ЭВМ. Недостаток, связанный с избыточностью переменных для описания ротационного движения тела, легко преодолевается путем введения в систему дифференциальных уравнений дополнительных слагаемых, обеспечивающих асимптотическую устойчивость решения и автоматическую коррекцию нормы кватерниона [26].
Тем не менее избыточность переменных Родрига — Гамильтона естественным образом мотивирует к использованию тех или иных вариантов модифицированных переменных Родрига — Гамильтона, число которых минимально и совпадает с числом степеней свободы тела, находящегося в сферическом движении.
(41)
Первая широко известная публикация на эту тему [30] предлагает следующий вариант введения модифицированных переменных Родрига — Гамильтона s?,s2,s3, называемых далее для краткости s-параметрами:
Si = Л</(1 - Ло) (i = 1, 2, 3). (43)
Формулы (43) позволяют рассматривать s-параметры как результат стереографического проектирования четырехмерной сферы, представляющей множество кватернионов Л = Лого + Лх«1 + Л2«2 + Лзгз, подчиненных условию нормировки ||Л|| = ЛО + Л? + Л2 + Л3 = 1, на трехмерную гиперплоскость s?,s2,s3, ортогональную оси ОЛо. В течение последующего десятилетия s-параметры оставались еще малоизвестными, а область их использования ограничивалась преимущественно демонстрационными примерами описания либрационного движения твердого тела [31, 32]. В работе [31] s-параметры упоминаются как компоненты модифицированного вектора Родрига.
Отчасти отсутствие интереса к практическому использованию s-параметров могло быть спровоцировано содержащимся в обширном обзоре [31] утверждением о том, что, если угол поворота тела не ограничен, например, как в случае вращающегося космического аппарата или очень маневренного самолета, сингулярность или разрыв непрерывности модифицированного вектора Родрига является неизбежным («.. .in the case of a spinning spacecraft or a very agile aircraft, the singularity or the discontinuity of the modified Rodrigues vector is unavoidable»).
В действительности s-параметры имеют лишь одну (!) точку вырождения и в этом плане они не только не хуже, а намного лучше, чем любые широко используемые углы конечных поворотов, имеющие плоскости вырождения. Как известно, наличие хотя бы одной точки вырождения при использовании любых трех переменных для задания ориентации твердого тела гарантируется общей теоремой, доказанной в [8]. Однако из наличия такой точки вырождения еще не следует неизбежность попадания точки четырехмерной единичной сферы, изображающей ориентацию твердого тела, в полюс отображения, находящийся в точке пересечения гиперсферы ||Л|| = 1 с осью ОЛо. Наоборот, попадание изображающей точки в полюс отображения является исключительной ситуацией, реализуемой лишь при специальном подборе параметров тела и начальных условий его движения.
Если же все-таки в процессе вычислений будет наблюдаться критический рост модуля ^ вектора s-параметров, то всегда можно будет переключиться с вектора s на вектор s * = —s/(s2 + s2 + s3), описывающий, как показано в [30], ту же самую ориентацию твердого тела.
С учетом вышеизложенных соображений в СПбГУ был выполнен анализ эффективности использования s-параметров для описания динамики ротационного движения твердого тела. Практическая важность этих исследований для кос-модинамики обусловлена тем, что во многих прикладных задачах приходится иметь дело с космическими аппаратами (КА), получившими значительную угловую скорость, например с целью создания искусственной силы тяжести на борту, или для стабилизации КА вращением, или в результате толчка, имевшего место вследствие срабатывания механизма отделения КА от последней ступени ракеты-носителя.
По аналогии с системой уравнений (36), (37) в работе [27] получены дифференциальные уравнения в s-параметрах. При этом динамические уравнения (36)
замыкаются кинематическими уравнениями:
¿\ (А31/А\ /и>ел
¿2 | = ЬВ I А32/В | — Бт I ^ | , (44)
\А33/С) \ше3)
в которых элементы матрицы А теперь выражены через ¿-параметры:
1 / и — 8(¿2 + ¿2) 4(2¿l¿2 — ¿3М1) 4(2¿l¿з + ¿2«1)\ А=- 4(2в1в2+ в3М1) + 4(2в2«з-«1М1) , (45)
\4(2sis3 - S2U1) 4(2s2s3 + S1U1) u — 8(s2 + s\) где u = (|s|2 + 1)2, |S|2 = s2 + s2 + s2, u1 = |s|2 — 1, а матрица B такова:
-1
S
1 — s2 — s2 + 1 2(S1S2 + S3) 2(S1S3 — S2)
В = (Щ) = —\ 2 (S1S2-S3) -si + sl-s^ + l 2(S2S3 + si) |. (46)
2(S1S3 + S2) 2(S2S3 — S1) —s2 — s2 + s2 + 1y
Далее, если, следуя [12], использовать «медленные» переменные L, р и а для изучения эволюции ротационного движения тела, то из (36), как и ранее, получаем динамические уравнения (38), а из (44) — кинематические дифференциальные уравнения:
М31/М ( —M2 \ S2 I = LB ( A32/В I — BT ( M1 I L-1. (47)
^3/ \A33/Cj \M2 ctg р)
Дифференциальная система (38), (47), в отличие от системы (38), (40), обладает минимальной размерностью. Содержащаяся в ней сингулярность при р = 0 может быть преодолена, как показано в [27], путем введения других «медленных» переменных, а именно проекций Lx ,Ly, Lz вектора L на оси кениговой системы координат вместо модуля L кинетического момента и двух углов р и а, определяющих направление вектора L. В качестве «быстрых» переменных при этом выступают s-параметры, определяющие ориентацию тела относительно системы координат CXYZ с помощью матрицы А = А(s 1, s2, S3). В результате выведены следующие кинематические уравнения в s-параметрах:
(1/А 0 0 \ /Ьх\ = В О 1 /В О Ат \Ly . (48)
V 0 0 1 /С) \Lzj
Система (42), (48) представляет собой замкнутую дифференциальную систему, удобную для описания ротационного движения тела. Эта система не имеет особенностей и не содержит тригонометрических функций и поэтому более предпочтительна, чем система (38), (47), для численного моделирования на ЭВМ.
В продолжение анализа эффективности использования s-параметров для описания динамики ротационного движения твердого тела в работе [27] проведен анализ структуры конфигурационного s-пространства уравновешенного твердого тела (M = 0). Для этого построен первый интеграл
2Т _ 1 64
Iß ~ с +
(s2 + s2)2 г( s2 s2\ /1 1 ... _
+ sl 4 + 4 + Ö-T (l*T-l)siS2S3
C Л A В \В A
+
16/, -и2 -,s2{ S2 s2 Л 16 S? + s2
имеющий место в условиях постоянства кинетической энергии Т и вектора кинетического момента Ъ и задающий неявным образом поверхность в пространстве в-параметров. Анализ поверхности (49), представляющей собой конфигурационное в-пространство уравновешенного твердого тела, позволил установить, что при произвольном фиксированном значении параметра 2Т/Ь2 эта поверхность представляет собой пару замкнутых ограниченных поверхностей (соответствующих двум эквивалентным наборам в-параметров: в и вв*), топологически эквивалентных тору. Следовательно, при вращении уравновешенного твердого тела абсолютные величины в-параметров ограничены. Неограниченный рост абсолютных величин в-параметров возможен только в случае, когда 2Т/Ъ2 = 1/С и начальные значения в-параметров лежат на оси Св3. Это неудобство, однако, легко преодолимо: достаточно сменить направления главных центральных осей инерции х и г или у и г на противоположные. При этом кинетическая энергия тела сохранит свое значение, а изображающая точка будет двигаться по окружности единичного радиуса с центром в точке С, лежащей в плоскости Св1в2. Такова практическая рекомендация, позволяющая избежать возможности неограниченного нарастания модуля в-вектора.
Для случая динамически симметричного твердого тела (А = В = С) в работе [27] произведено полное аналитическое исследование конфигурационного пространства. Установлено, что это пространство, определяемое соотношением
2 Т
Т2
1 ю. 2 2
16 i2
|s|2 - 1)2 + 4s2, 4(sl + s22) — u
A
+
C
представляет собой два правильных тора, вложенных один в другой. Каждый из торов можно запараметризовать углами а и ß в виде
cos ß sm ß
si =-(1 + cos во cos a), S2 =-(1 + cos Oq cos a), S3 = ctg Oq sina, (50)
sin в0 sin во
где а € [—п; п], ß € [0; 2п]. Параметр во € (0; п/2), определяемый соотношением
2T/L2 — 1/C
sin2 2во
1/A — 1/C
принимает два различных значения, лежащих в пределах от 0 до п/2 и определяющих форму торов (50). Найдены зависимости углов а и в от времени:
а = 2 arctg ■
'1 + cos $0 1 — cos во
tg
Lki
(t + to)
ß = ßo + (51)
Здесь
= 2 sin2 в0
+
C
cos2 во sin2 во = 2 ( ----h
A
C
1
C'
Тем самым для динамически симметричного твердого тела выполнено аналитическое интегрирование уравнений ротационного движения тела в в-параметрах. Из формул (51) следует, что угол а изменяется во времени периодически с периодом Та = 4п/(Ък\). Угол в тоже изменяется периодически с периодом Тр = 4п/(Ък2), поскольку в € [0; 2п] и при достижении значения 2п угол в скачком переходит к значению 0.
4
1
В работе [29] доказано, что поскольку периоды Та и Т^, вообще говоря, несоизмеримы, то траектория изображающей точки вращающегося твердого тела покрывает конфигурационное пространство всюду плотно. Дальнейшие исследования показали, что это свойство является важным в динамике ротационного движения тела. Оно позволяет упростить процедуру усреднения непрерывных на конфигурационном в-пространстве функций в-параметров, зависящих от ориентации тела (например, возмущающего момента М), по невозмущенному движению Эйлера — Пуансо путем замены усреднения по времени на усреднение (с некоторым весом) по конфигурационному в-пространству. В качестве весовой функции усреднения выбирается относительное время пребывания изображающей точки в пределах элемента ¿Б конфигурационного пространства, равное пределу Ит^-юо где Т — время движения, а ¿Т — суммарное время пребывания изображающей точки в пределах элемента ¿Б за время Т. В результате среднее значение некоторой функции I(в!,в2,в3) по невозмущенному движению Эйлера — Пуансо определяется выражением
где в 1, в2, в3 — известные на основании (50) функции углов а и в. Таким образом, на базе установленных эргодических свойств конфигурационного в-пространства обосновано представление об в-параметрах как идеальных «быстрых» переменных и предложен практически удобный вариант модификации метода усреднения, предполагающий переход от усреднения по невозмущенному движению к усреднению по конфигурационному пространству. В [29] данный подход применен в задаче об эволюции ротационного движения заряженного твердого тела в геомагнитном поле при движении по околоземной орбите, регрессирующей вследствие полярного сжатия Земли.
Завершая рассмотрение вопроса об эффективности использования в-параметров для описания динамики ротационного движения твердого тела, отметим, что в [27] проведен анализ особенности в-параметризации. Как отмечалось выше, эта особенность представляет собой единственную точку. Однако попадание изображающей точки в полюс отображения в процессе движения является исключительной ситуацией, требующей особых начальных условий, причем множество этих начальных условий в пространстве в-параметров имеет нулевую меру. Действительно, существенные неприятности, связанные с большими значениями модуля в-вектора при данном выборе систем координат возможны лишь в том случае, когда кинетическая энергия близка к максимальному значению (в невозмущенном движении это соответствует 2Т/Ъ2 = 1/С), а начальные значения в-параметров близки к точкам, лежащим на оси Св3. Однако указанное обстоятельство легко преодолимо. Достаточно систему в-параметров ввести так, чтобы ориентация, соответствующая полюсу отображения в-параметров, находилась по возможности дальше от начальной ориентации твердого тела. На практике это означает, что достаточно поменять обозначения главных центральных осей инерции тела. Следовательно, ротационное движение твердого тела, как близкое к невозмущенному, можно также описывать при помощи в-параметров, не опасаясь при этом, что изображающая точка может попасть в полюс.
Другим важным следствием установленной в [27] ограниченности конфигурационного пространства вращающегося твердого тела является снятие необходимости
п 2п
(Лз = ! J 1 (в1(а,в),в2 (а,в),в3(а,в))
эш в0 ¿айв
4п2(1 + соэ во соэ а)'
использования предлагаемого в [30] переключения с одного набора в-параметров на другой, соответствующий той же самой ориентации тела, но лежащий внутри единичной сферы в пространстве в-параметров (соответственно в и в *). Предложенное в [30] переключение было введено с целью предотвратить возможный неограниченный рост модуля в-вектора. Вместе с тем такое переключение может оказаться полезным, поскольку избавляет от необходимости заботиться о выборе полюса в-параметров.
В целом после публикации [27] были введены в широкий практический обиход предложенные в ней новые варианты математических моделей динамики твердого тела, построенные на базе трехмерных модифицированных параметров Родрига — Гамильтона ( в-параметров).
Продолжение общетеоретических исследований, начатых в [27], содержится в работах [33] и [34]. В первой из них установлена возможность использования в-пара-метров как обобщенных переменных для построения лагранжева формализма в динамике твердого тела. Соответствующие дифференциальные уравнения движения построены в форме канонических уравнений Гамильтона. Для этого вначале получено следующее выражение кинетической энергии тела с тензором инерции Л = ¿1^(А,В,С):
128 т
Т= — (¿1, ¿2, ¿з ) В Л Вт ( ¿1, ¿2, ¿з ) • (52)
Здесь матрица В, как ранее, определяется формулой (46). Затем введены канонические переменные в1, в 2, вз и соответствующие им обобщенные импульсы Р1,Р2,Рз. В предположении, что тело находится под воздействием только потенциальных сил и моментов и потенциальная энергия П не зависит от обобщенных скоростей в1, в 2, вз,
дТ
найдены обобщенные импульсы р^ = ——, (г = 1, 2, 3) в виде
д в^
, ,т 256 т , ,т
[Р1, Р2, Рз ) = —ВЛВ ( вь в2, «з ) •
Последнее равенство может быть разрешено относительно обобщенных скоростей
(в1, в2, вз )т = ВЛ-1 Вт( Р1, Р2, рз )т (53)
и подставлено в (52). Тогда функция Гамильтона примет вид
Н=(р1, Р2, Рз ) (¿1, ¿2, ¿з )Т —Т+И= ^ р2, рз ) В Л ХВТ (ръ р2, рз )Т+П, а уравнения (53) могут быть переписаны следующим образом:
( Р1\-( — — —V (54)
¿3 У у II ) V дР1' дР2' дрз )
Совместно со второй группой канонических дифференциальных уравнений
( дН дН дН \Т
= 4^? дГз) (55)
они образуют замкнутую систему канонических дифференциальных уравнений, описывающих динамику вращательного движения твердого тела в -параметрах.
В статье [34] выведена новая форма дифференциальных уравнений для описания ротационного движения твердого тела. В качестве вспомогательной системы координат вместо вышеупомянутой Ь1,Ь2, Ь3 используется похожая по сути, но более наглядная система координат Ш1,Ш2,Ш3, ось Ш3 которой направлена вдоль угловой скорости тела. Модуль угловой скорости ш и два угла а, в, задающие направление вектора угловой скорости в инерциальной системе координат, образуют множество «медленных» переменных, определяющих угловую ориентацию вспомогательной системы координат относительно инерциальной. В качестве «быстрых» переменных используются три в-параметра, определяющих угловую ориентацию тела (или, что то же, системы главных центральных осей инерции Схуг) относительно вспомогательной системы координат Ш1,Ш2,Ш3. Соответствующая матрица направляющих косинусов А = (Л{у) формально определяется тем же равенством (45), но смысл в-параметров в ней, естественно, другой. Выведенные в [34] динамические уравнения определяются равенствами
и? \ ( Mx/A ша sin в ) = A I My/В | + ш / \ Mz/С
и2 ABC
A13
(A-B)(B-C)(A-C)A31A32 I A23 | +
A
33
AC (A - С )A2i A31 + BC(B - С A22A32 + | -AC(A - C)AiiA3! - BC(B - CA2A32
0
(56)
Совместно с кинематическими уравнениями
¿Л ( а sin в
¿2 ) = Вт I -в | (57)
s3 J \ ш - а cos в
они образуют замкнутую дифференциальную систему (56), (57).
Проанализировано конфигурационное пространство твердого тела, находящегося под уравновешенным внешним воздействием (M = 0). Для этого построен первый интеграл
16 .о . 16
u
^(2sis3 - s2Ul)2(A2 - AD) + —(2s2s3 + slUl)2(B2 - BD)+ 22
+ {1~-u(sl + sl?j (C2-CD)= 0, (58)
где D = L2 / (2T) = const, имеющий место в условиях постоянства кинетической энергии T и вектора кинетического момента L и задающий неявным образом поверхность в пространстве s-параметров. Эта поверхность представляет собой конфигурационное s-пространство уравновешенного твердого тела. Поскольку качественная зависимость интеграла (58) от переменных S1,S2,S3 совпадает с аналогичной зависимостью в интеграле (49), то результаты анализа конфигурационного s-пространства уравновешенного твердого тела совпадают с изложенными выше. В случае динамически симметричного твердого тела (A = B = C) третье из уравнений (56) ожидаемо дает результат ш = const, а первые два уравнения системы (56) принимают вид
р \ = ш{А - С) f А21А31 + А22А32 \ asin/3 J А у —АцА31 — Ai2A32 )' [ >
Показано, что в этом случае конфигурационное пространство представляет собой правильный тор, а его параметризация дана в явном виде в [34].
Далее, на базе дифференциальных уравнений (56), (57) исследована эволюция ротационного движения твердого тела в среде с сопротивлением, создающим момент М = — Нш, где к > 0. Построены траектории конца вектора ш и траектории конца вектора в, наглядно подтверждающие удобство использования переменных ш, а, в, в1, в2, вз для качественного и количественного анализа ротационного движения твердого тела.
4. Стационарные движения двух гравитирующих тел. В работах С. Н. Кирпичникова [35-37] рассмотрены общие вопросы устойчивости стационарных движений двух гравитирующих тел с различными инерционными характеристиками. Найдены множества всех стационарных относительных движений и исследована их орбитальная устойчивость в общем случае движения двух гравитирую-щих тел: осесимметричного и сферически симметричного. Задача изучается в общем виде: используется точное выражение для потенциала гравитационного взаимодействия и никаких предположений, кроме предположения о его осесимметричности, не делается. Существенную упрощающую роль при исследовании сыграл специальный выбор лагранжевых позиционных координат, отвечающих некоторой «полуподвижной» декартовой системе координат.
Относительное движение рассматриваемых тел, при котором центр сферически симметричного тела во все время движения неподвижен в каком-либо пространстве, жестко скрепленном с осью симметрии осесимметричного тела, названо стационарным. Доказана теорема о двух типах стационарных движений, которая позволила впервые описать множества всех возможных стационарных движений и дать классификацию этих движений, связанную с их характерными свойствами. При стационарных движениях осесимметричное тело совершает регулярную прецессию, а его центр масс описывает круговую орбиту, причем центр сферически симметричного тела у стационарных движений I типа, вообще говоря, не лежит в плоскости этой орбиты, а у стационарных движений II типа всегда является ее центром.
Найденные стационарные движения двух типов охватывают в качестве частных случаев все известные движения («спица», «стрела», «квазистрела», «поплавок» и др.), определенные ранее [12] при тех или иных предположениях для некоторых конкретных тел, как правило, простого вида. Отказ от общепринятых весьма частных моделей гравитационного поля осесимметричного тела привел к качественно новым выводам. Например, в отличие от известных результатов доказано, что движения «стрела» возможны и у несимметричных, а «квазистрела» — у симметричных относительно плоскости экватора осесимметричных тел.
Разработан новый метод исследования орбитальной устойчивости компланарных (ось симметрии осесимметричного тела и центр сферически симметричного движутся в одной неподвижной плоскости) стационарных движений, заключающийся в установлении связи между их устойчивостью в рамках компланарной и некомпланарной постановок задач. Определено, когда метод позволяет свести исследование устойчивости к анализу существенно более простых условий устойчивости компланарной задачи.
При исследовании орбитальной устойчивости стационарных движений большое внимание уделено качественным аспектам: определению числа и типов соответствующих точек либрации, их возможным бифуркациям, нахождению параметров, от-
вечающих переходам от устойчивости к неустойчивости, и т. д. Доказан ряд теорем об орбитальных устойчивости и неустойчивости стационарных движений в некомпланарной задаче. В качестве примера приведем лишь один результат: если масса сферически симметричного тела меньше определенного критического значения, то невозможна орбитальная устойчивость стационарных движений с минимумом эффективной потенциальной энергии. Изучена устойчивость стационарных движений, близких к «поплавку», как для симметричных относительно экваториальной плоскости осесимметричных тел, так и для несимметричных.
Исследованы компланарные стационарные движения. Доказаны теоремы, устанавливающие связь орбитальной устойчивости в рамках компланарной и некомпланарной постановок задач. Так, области орбитальных устойчивости и устойчивости в линейном приближении для неколлинеарных (центр сферически симметричного тела не лежит на оси симметрии осесимметричного) стационарных движений всегда сохраняются, а для коллинеарных — могут уменьшиться при переходе к некомпланарной постановке задачи.
Движение простых модельных тел определяется более простыми соотношениями, допускающими более глубокий анализ. Поэтому полученные результаты были проанализированы в случае аппроксимации поля тяготения осесимметричного тела полем тяготения нескольких материальных точек, среди которых могут быть точки с положительными и отрицательными массами, а также пары точек с комплексно сопряженными массами. Для случая аппроксимирующей гантели в рамках компланарной задачи изучены множества стационарных движений, их орбитальная устойчивость, влияние на нее сжатия и асимметрии гантели, возможные бифуркации, определены критические расстояния между телами, которым отвечает переход от устойчивости к неустойчивости и наоборот.
Опишем лишь один из результатов С. Н. Кирпичникова. Широко известна и используется на практике система пассивной гравитационной стабилизации вращательного движения круговых спутников планет в рамках обычной спутниковой модели. Учет конечности массы спутника в рамках поступательно-вращательного его движения приводит к качественно новому явлению: соответствующие стационарные движения становятся орбитально неустойчивыми, если расстояние между телами меньше определенного критического значения. При отношении массы осе-симметричного тела к массе сферически симметричного тела порядка единицы это значение имеет порядок характерного размера осесимметричного тела. С ростом указанного отношения масс соответствующее критическое расстояние неограниченно увеличивается.
Поэтому выводы о возможных стационарных относительных движениях взаимно гравитирующих тел и их устойчивости имеют важные конкретные приложения в астрономии и космодинамике. Так, они обосновывают возможность длительного существования разного рода двойных космических объектов, например двойных астероидов, комет, звездных систем, систем «спутник-планета», и определяют устойчивые относительные конфигурации при стационарных движениях образующих их тел. Другая область практических приложений относится к космодинамике. При достаточно общих предположениях дан ответ на важный уже в практике современных космических исследований вопрос: на каких относительных орбитах космический аппарат (КА) может образовывать устойчивую двойную систему с астероидом или спутником планеты, в частности в каких точках КА может в течение длительного времени пассивно «зависать» над поверхностью астероида или спутника? Важно,
что расстояние между телами такого двойного объекта может иметь порядок линейных размеров самих тел, а их массы могут быть сравнимы, и хорошо развитые теории для классических спутниковых моделей здесь не применимы.
5. Управление ориентацией твердого тела. В конце рассмотренного в первой части обзора 250-летнего периода развития научных исследований в СПбГУ [1] в динамике твердого тела наметилась явная тенденция к переходу от изучения классических задач к изучению задач динамики управляемого движения твердого тела, а проблема управления ориентацией твердого тела стала одной из важнейших. Вызванное запросами практики и, в частности, необходимостью активного управления угловым движением летательных и космических аппаратов, это новое направление в динамике твердого тела потребовало для своего развития привлечения последних достижений из теории управления, теории устойчивости движения, теории нелинейных колебаний и других областей фундаментальной и прикладной математики. В то же время передовые достижения математики нашли широчайшее поле для своего применения в динамике управляемого движения твердого тела, что отражается в огромном количестве публикаций, число которых продолжает расти. Некоторые результаты этих двух взаимообогащающих тенденций рассматриваются в данном разделе.
В. И. Зубов рассмотрел задачу активного управления ориентацией твердого тела в общей нелинейной постановке [38-47]. Он предложил новые подходы к синтезу управляющих моментов, обеспечивающих стабилизацию заданных ориентаций тела. Кроме того, им проведено полное исследование качественного поведения решений соответствующих замкнутых систем управления.
Рассмотрим подход В. И. Зубова к задаче об одноосной стабилизации тела [46]. Пусть твердое тело вращается вокруг своего центра масс О с заданной угловой скоростью ш. Обозначим через Л тензор инерции тела в главных центральных осях инерции тела Охуг. Вращательное движение тела под действием управляющего момента М описывается уравнениями Эйлера:
Пусть заданы единичные векторы в и г, причем вектор в неподвижен в инерциаль-ном пространстве, а вектор г неподвижен в системе координат, жестко связанной с телом. Тогда вектор в поворачивается относительно системы координат Охуг с угловой скоростью —ш. Следовательно,
Для дифференциальной системы, состоящей из динамических уравнений Эйлера (60) и кинематических уравнений Пуассона (61), требуется построить управляющий момент М, обеспечивающий одноосную стабилизацию тела: система (60), (61) должна допускать асимптотически устойчивое положение равновесия
Лш + ш х Лш = ММ.
(60)
в
ш х в.
(61)
ш = 0,
в = г.
(62)
В. И. Зубов предложил выбирать управляющий момент в виде
М = —ш + кг х в,
(63)
где к = const > 0. Следует отметить, что при таком управлении система (60), (61) наряду с (62) имеет также положение равновесия
ш = 0, s = —r. (64)
С помощью функции Ляпунова
V = Joj + k\\s - г)\2) (65)
доказана следующая теорема [46].
Теорема 1. Пусть управляющий момент определяется формулой (63). Тогда положение равновесия (62) асимптотически устойчиво, а положение равновесия (64) неустойчиво. Кроме того, любое движение .замкнутой системы, отличное от положения равновесия (64), при соответствующем выборе коэффициента к обладает свойством ш ^ 0, r ^ s при t ^ж.
В. И. Зубов рассматривал также задачу сканирования оси тела по заданной программе [44, 46]. Предполагалось, что единичный вектор so(t) вращается в инерци-альном пространстве с заданной угловой скоростью wo(t). Управляющий момент должен обеспечивать выполнение условий uj(t) ^ ¿b(t), r(t) ^ so(t) при t ^ ж>. Доказано, что искомое управление можно выбрать в виде [44, 46]
M = w0(t) — ш + Лш0(t) + w0(t) х Лш + кг х s0(t), к = const > 0.
Аналогичные подходы были развиты для задачи трехосной стабилизации твердого тела [40, 43, 44].
5.1. Строгие функции Ляпунова в динамике твердого тела. Результаты В. И. Зубова по управлению ориентацией твердого тела основаны на построении слабых функций Ляпунова. Производные этих функций в силу изучаемых систем только неотрицательны. Известно, что такие функции Ляпунова плохо подходят для анализа робастности, так как их отрицательные полуопределенные производные вдоль траекторий могут становиться положительными при сколь угодно малых возмущениях динамических уравнений. Это обстоятельство мотивировало развитие методов построения строгих функций Ляпунова, т. е. функций с отрицательно определенными производными. Е. Я. Смирнов предложил способ преобразования слабых функций Ляпунова, построенных В.И.Зубовым, в строгие [48].
Например, в задаче об одноосной стабилизации твердого тела он предложил модификацию функции Ляпунова (65)
Vi = ^(ojTJoj + k\\s - г]|2) -7сЗтЛР(г х s), (66)
где y — положительный параметр; P — постоянная матрица. Доказано, что если матрица (PT + P) — положительно определенная, то для сравнительно малых значений Y производная функции V в силу системы (60), (61), (63) является отрицательно определенной относительно переменных ш и r — s.
На базе подхода, разработанного Е. Я. Смирновым, в работах его научной группы [48-50] был получен ряд новых результатов. В частности, разработаны методы
робастного управления угловым движением твердого тела для случаев, когда тензор инерции тела и моменты, действующие на тело, заданы с некоторой погрешностью [50].
Строгие функции Ляпунова строились также в работах А. Ю. Александрова и А. А. Тихонова [51-55] для решения задач об одноосной и трехосной стабилизации твердого тела. В работе [51] эти задачи рассмотрены для тела с изменяющимися во времени моментами инерции и найдены достаточные условия асимптотической устойчивости заданных ориентаций. Наряду с применением специальных конструкций строгих функций Ляпунова использовался метод дифференциальных неравенств. В работах [52, 54] с помощью метода сравнения для функций Ляпунова проанализирована стабилизация твердого тела по части переменных. Получены достаточные условия устойчивости и асимптотической устойчивости программных движений твердого тела по части переменных. Подход Е. Я. Смирнова к построению строгих функций Ляпунова [48] можно использовать только для случая линейного восстанавливающего момента. Поэтому в работе [55], посвященной использованию линейных и нелинейных управлений в задаче о стабилизации тела в условиях воздействия нестационарных возмущений с нулевыми средними значениями, предложена другая конструкция строгой функции Ляпунова.
5.2. Эволюционирующее нелинейное управление. В работах [51, 52, 54, 56-58] выполнен ряд исследований, развивающих задачу реактивного управления угловым положением твердого тела в направлении использования восстанавливающего и диссипативного моментов, зависящих явным образом не только от фазовых переменных, но и от времени. Постановка вопроса о возможности реализации таких систем управления, в которых восстанавливающая или диссипативная составляющая управляющего момента стремится к нулю с течением времени, является важной и закономерной в связи с ограниченностью ресурсов реактивных систем управления. С математической точки зрения варианты исчезающего с течением времени управления известны и общепризнаны в мировой науке [59, 60] как наиболее сложные в проблеме анализа устойчивости механических систем с нестационарным параметром. Рассмотрены варианты линейного и существенно нелинейного управления. С помощью прямого метода Ляпунова и метода сравнения получены условия, при которых можно гарантировать устойчивость заданных положений равновесия тела, несмотря на эволюцию управляющих моментов.
В работах [56, 57] рассматривается задача о трехосной стабилизации тела, с которым неизменно связаны орты ri,M2,M3, в базисной системе координат с ортами si, S2, S3, принимаемой за инерциальную. Для дифференциальной системы, состоящей из динамических уравнений Эйлера (60) и кинематических уравнений Пуассона
si = —ш х si, i = 1, 2,3, (67)
решается задача обеспечения асимптотической устойчивости положения равновесия
ш = 0, si = ri, i = 1,2,3, (68)
с помощью суперпозиции восстанавливающего момента Mr = —aisi х ri — 0,2 S2 х r2 и эволюционирующего со временем диссипативного момента Add = —h(t)giadW(ш), где W(ш) — непрерывно дифференцируемая при ш G R3 положительно определенная однородная порядка р +1 функция, р ^ 1. Таким образом, уравнения (60)
принимают вид
Лш + ш х Лш = — h(t)gradW(ш) — a1s1 х r1 — a2s2 х r2,
где h(t) — положительная и непрерывно дифференцируемая при t ^ 0 функция. В работе [56] предполагается, что h(t) ^ 0 при t ^ Доказана лемма об оценке снизу нормы восстанавливающего момента в окрестности стабилизируемого движения твердого тела, а также две теоремы об асимптотической устойчивости стабилизируемого движения тела. Первая из них относится к случаю линейного диссипа-тивного момента. В этом случае функция W(ш) представляет собой квадратичную форму W(ш) = ojtDoj с постоянной симметричной положительно определенной матрицей D, а уравнения (60) имеют вид
Лш + ш х Лш = —2h(t)Dш — aisi х ri — a2S2 х r2.
Теорема 2. Пусть h(t) ^ 0 при t ^ +ж. Если существует число v ^ 1 такое, что \h(t)\ < Lh^(t) при t > 0, L = const > 0 и f0 ЪУ{т)(1т —> +oo при t —> +oo; mo положение равновесия (68) системы (60), (67) асимптотически устойчиво.
Следствием из теоремы 2 является асимптотическая устойчивость положения равновесия (68) системы (60), (67), если h(t) = (t + 1)-^ и в ^ 1. Доказано также, что неравенство в ^ 1 представляет собой не только достаточное, но и необходимое условие асимптотической устойчивости положения равновесия.
Следующая теорема, доказанная в [56], относится к модели управления с существенно нелинейным диссипативным моментом, в котором W (ш) — непрерывно дифференцируемая при ш G R3 положительно определенная однородная порядка ¡i +1 функция, причем ¡л > 1.
Теорема 3. Пусть h(t) ^ 0 при t ^ +ж. Если существует число v ^ 1 такое, что \h(t)\ < Lh1+~(t) при t > 0, L = const > 0 и ЬУ{т)<1т —> +оо при t —> +оо; то положение равновесия (68) системы (60), (67) асимптотически устойчиво.
Сформулированное выше следствие из теоремы 2 остается в силе при управлении с существенно нелинейным диссипативным моментом. Этот факт указывает на то, что найденные в теоремах 2 и 3 достаточные условия асимптотической устойчивости близки к необходимым. Таким образом, получены критерии, позволяющие гарантированно решать задачу угловой стабилизации тела в условиях экономии ресурсов реактивных систем управления.
В работе [57] рассматривается в некотором смысле противоположная по отношению к [56] задача о трехосной стабилизации тела. Предполагается, что h(t) ^ при t ^ +ж. Доказаны две теоремы об асимптотической устойчивости.
В работах [52, 61] решена задача об одноосной стабилизации твердого тела в положении равновесия (62) с помощью суперпозиции эволюционирующего со временем восстанавливающего момента и не зависящего явно от t линейного [61] или нелинейного [52] относительно угловой скорости диссипативного момента. В случае линейного диссипативного момента динамические уравнения (60) управляемого движения тела имеют вид
Лш + ш х Лш = —Dш — ah(t)\\s — х r, (69)
где О — постоянная симметричная положительно определенная матрица; Н(Ь) — положительная и непрерывная при Ь ^ 0 функция, а> 0, р ^ 1. С помощью метода декомпозиции доказана теорема.
Теорема 4. Пусть Н(Ь) ^ 0 при Ь ^ Если /0 Н(т)йт ^ при Ь ^ то положение равновесия (62) системы (69), (61) асимптотически устойчиво.
Следствием из теоремы 4 является устойчивость положения равновесия (62) системы (69), (61), если Н(Ь) ^ 0 при Ь ^
В работе [52] та же задача об одноосной стабилизации решена при нелинейном относительно угловой скорости диссипативном моменте в уравнениях движения:
Лш + ш х Лш = (ш) — аН(г)\в — х г.
Рассмотрены два подхода к построению функции Ляпунова. На основании первого из них с учетом обозначения ^(Ь) = шах |0; к(Ь)/к(Ь)^ сформулирована и доказана следующая теорема.
Теорема 5. Если существует постоянная Ь > 0, такая что /0 т)йт ^ Ь для Ь ^ 0, то положение равновесия (62) системы (69), (61) устойчиво по отношению к ш.
Доказана также теорема об асимптотической устойчивости положения равновесия по отношению к ш.
Другой подход к построению функции Ляпунова позволил найти достаточные условия асимптотической устойчивости не только по отношению к ш, но и по отношению к в. С использованием метода дифференциальных неравенств для функций Ляпунова доказана следующая теорема.
Теорема 6. Если Н(€) ^ 0 для Ь ^ 0 и существуют положительные числа в и Ь такие, что < ЬЬ1-^^) для £ > 0, /liЯгr(í) ¡0 1гв+^+1 (т) (1т ->• +аэ
при Ь ^ то положение равновесия (62) системы (69), (61) асимптотически
устойчиво.
Проблема трехосной стабилизации углового положения твердого тела с использованием эволюционирующих во времени нелинейных восстанавливающих моментов изучалась в работах [54, 58]. Анализировалась возможность построения системы управления, в которой восстанавливающий момент стремится к нулю с течением времени, и остается единственный управляющий момент — диссипативный (линейный или нелинейный), не зависящий явно от времени. Рассматривались случаи линейного и существенно нелинейного восстанавливающих моментов. С помощью прямого метода Ляпунова и метода сравнения получены условия, при которых можно гарантировать устойчивость или асимптотическую устойчивость равновесного положения тела, несмотря на исчезновение восстанавливающего момента.
В работе [54] применяется метод сравнения и метод декомпозиции в варианте, разработанном в [62, 63]. Сначала рассмотрен вариант линейного восстанавливающего момента. Управляемое движение тела в этом случае описывается уравнениями:
Лш + ш х Лш = О ш — И,(Ь)(а181 х т\ + а2в2 х г2), (70)
где О — постоянная симметричная отрицательно определенная матрица. Вводится в рассмотрение функция ф(Ь) = шах |2Л*, к(Ь)/к(Ь)^, где Л* — наибольшее собственное значение матрицы Л ^ОЛ. Доказаны следующие две теоремы.
Теорема 7. Пусть функция Н(€) положительна и непрерывна для Ь € [0, и к(Ь) ^ 0 при Ь ^ Если существует постоянная Ь > 0 такая, что
/0 ф(т)йт ^ Ь при Ь ^ 0, то положение равновесия (68) системы (70), (67) устойчиво по отношению к ш.
Теорема 8. Пусть функция Н(€) положительна и непрерывна для Ь € [0, и Н(€) ^ 0 при Ь ^ Если /0 ф(т)йт ^ при Ь ^ то положение равновесия (68) системы (70), (67) асимптотически устойчиво по всем переменным.
Затем рассмотрено управляемое движение тела с нелинейным восстанавливающим моментом:
Лш + ш х Лш = Ош — (в1, в2)(а,1в1 х т\ + а2в2 х Т2), (71)
где п(в\,в2) = (а1\\в\ — т1\2 + а2\\в2 — г2\\2)/2 и V > 0. Доказаны две теоремы — одна об устойчивости по отношению к ш, а другая — об асимптотической устойчивости по всем переменным.
В работе [58], в отличие от [54], рассмотрено управление с нелинейным демпфированием:
Лш + ш х Лш = (ш) — к(Ь)пХ(в1, вв2)(а1в1 х г1 + а2в2 х г2), (72)
где функция к(Ь) положительна и непрерывно дифференцируема для Ь € [0, и к(Ь) ^ 0 при Ь ^ Применены развитые в [52] подходы к построению функций Ляпунова. Доказано, что если р = 2Л +1, то при выполнении условий теоремы 6 положение равновесия (62) системы (72), (61) асимптотически устойчиво по всем переменным.
5.3. Развитие приема усреднения для задач динамики твердого тела. При построении математической модели, описывающей динамику твердого тела, мы неизбежно ограничиваемся учетом основных действующих сил и моментов, относя все остальные к категории возмущающих. В ряде случаев структура возмущающих моментов известна и их влияние можно учесть с помощью методов теории возмущений, которые эффективны для широкого круга динамических задач. В частности, они широко и успешно используются в задачах динамики искусственных и естественных небесных тел [21], о чем упоминалось выше в разделе 1. В этом случае наибольшая трудность может возникнуть из-за изменяющихся во времени возмущений. Мощным инструментом анализа устойчивости нестационарных систем является прием усреднения, позволяющий свести исследование устойчивости нестационарных систем к исследованию соответствующих стационарных усредненных систем с вытекающими из этого существенными упрощениями. Однако применение приема усреднения хорошо развито лишь для случая, когда рассматриваемые системы являются быстро изменяющимися во времени. В этом случае техника усреднения широко применяется в различных задачах динамики и управления [21, 29].
В работах [64, 65] предложен оригинальный подход к построению функций Ляпунова и на его основе найдены новые условия устойчивости нелинейных нестационарных систем. По сравнению с известными условиями устойчивости, полученными
с помощью различных методов, использующих усреднение, принципиальная новизна результатов [64, 65] заключается в том, что для обеспечения асимптотической устойчивости рассматриваемых систем на правые части дифференциальных уравнений не накладывается требование быстрого изменения во времени.
В работе [53] исследуются динамика и управление угловым положением твердого тела, находящегося под воздействием нестационарного возмущающего момента. Никаких ограничений на скорость изменения во времени возмущающего момента не налагается. Исследуется задача одноосной стабилизации твердого тела. Управляющий момент состоит из линейной диссипативной составляющей и нелинейной
восстанавливающей составляющей. Возмущающий момент М представлен в виде линейных комбинаций однородных функций с изменяющимися во времени коэффициентами: М = Ф(г)<?(V — г). Здесь Ф(£) € К3хк — непрерывная и ограниченная для £ ^ 0 матрица, а компоненты вектора Я (и) € — непрерывно дифференцируемые для и € К3 однородные функции порядка р. Предполагается также, что
1 Гг+т
- / Ф(т) <1т ->• 0 при Т ->• +оо (73)
равномерно по £ ^ 0.
Последнее условие выполняется для периодических, квазипериодических и почти периодических функций с нулевыми средними значениями. Важно отметить, что на величины амплитуд этих возмущений никаких ограничений не накладывается. Возмущающие моменты такого типа часто встречаются в математических моделях углового движения ИСЗ, взаимодействующих с геофизическими полями [66].
Тогда дифференциальные уравнения вращательного движения тела имеют вид
Лш + ш х Лш = Шш — а||в — х г + Ф(£)<?(V — г). (74)
Здесь к > 0, р ^ 1, а > 0, О — постоянная симметричная отрицательно определенная матрица. Таким образом, предполагается, что порядки однородности возмущающего и восстанавливающего моментов совпадают, а коэффициенты восстанавливающего момента имеют нулевые средние значения.
В случае р =1, соответствующем линейным восстанавливающему и возмущающему моментам, Ф(£) € К3х3 и система (74) упрощается:
Лш + ш х Лш = КО ш — ая хг + Ф(£)(в — г). (75)
С использованием метода функций Ляпунова на базе подходов, развитых в работах [62, 64, 67], доказана следующая теорема, утверждающая, что дестабилизирующее влияние нестационарного возмущения можно компенсировать выбором достаточно большого множителя при диссипативном моменте.
Теорема 9. Существует число к0 > 0 такое, что положение равновесия (62) системы (75), (61) асимптотически устойчиво для всех К ^ К0.
Кроме того, доказано, что при существенно нелинейном восстанавливающем моменте (р > 1) не нужно подбирать параметр К > 0 достаточно большим для стабилизации положения равновесия (62): его асимптотическая устойчивость гарантируется при любом К > 0.
В работе [68] рассматривается более широкий по сравнению с [53] класс нестационарных возмущений, для которых неприменимы подходы, использовавшиеся в [53]. При этом получены условия, гарантирующие решение задачи одноосной стабилизации для более широкого класса нестационарных возмущений с нулевыми средними значениями, т. е.
1 Гт
- Ф(г) (¿г —^ 0 при Т —>• +оо. (76)
Во-первых, не требуется, чтобы предельное соотношение (73) выполнялось равномерно при Ь ^ 0 и, во-вторых, порядок однородности а возмущений может не только совпадать с порядком однородности р восстанавливающего момента, но и быть меньше его. Заметим, что если а > р, то асимптотическую устойчивость положения равновесия (62) системы (74), (61) гарантируют известные результаты об устойчивости по линейному и нелинейному приближениям [48, 69]. Доказана следующая теорема.
Теорема 10. Пусть р > 1 и существует число в € (0,1] такое, что 1 Т
/0 Ф(т) ¿т —> 0 при Т —> +оо. Тогда, если 2а ^ р + 1 + ¡3(р — 1); то положение равновесия (62) системы (74), (61) асимптотически устойчиво.
В качестве следствия из теоремы 10 вытекает, что если р > 1 и выполнено условие (76), то для асимптотической устойчивости положение равновесия (62) системы (74), (61) достаточно выполнения неравенства а ^ р. Это следствие обобщает результат работы [53] на более широкий класс нестационарных возмущений. Однако, в отличие от [53], где установлена равномерная асимптотическая устойчивость положения равновесия, теорема 10 не позволяет сделать вывод о том, что асимптотическая устойчивость будет равномерной.
Практическая значимость теорем, доказанных в [53, 68], будет рассмотрена в одной из следующих частей обзора. Здесь же продолжим анализ общетеоретических результатов, являющихся развитием исследований, выполненных в [53, 68], но относящихся к другой задаче — задаче о трехосной стабилизации твердого тела в условиях воздействия нестационарных возмущений, рассмотренной в [55, 70].
В работе [70] динамические уравнения имеют вид
Лш + ш х Лш = Б ш — от)"(81, вв2)(а181 х Т1 + а2в2 х Г2) + Ф(Ь)(((в1 — гь в2 — Т2), (77)
где сохраняется смысл обозначений, использованных ранее в уравнении (71); о — положительная постоянная; матрица Ф(Ь) — непрерывная, ограниченная и удовлетворяет условию (76); компоненты вектора (( — непрерывно дифференцируемые однородные функции порядка 2v +1. Известно, что если восстанавливающий и возмущающий моменты линейны (V = 0), то возмущения рассматриваемого типа могут разрушить устойчивость положения равновесия. В работе [53] показано, что если порядок неоднородности рассматриваемой системы больше, чем 1, то нестационарные возмущения с нулевыми средними значениями не разрушают асимптотическую устойчивость. Однако в работе [53] при рассмотрении похожей задачи на функцию Ф(Ь) накладывалось более жесткое ограничение (73). Поэтому результаты работы [53] в данной задаче неприменимы и используются другие подходы, основанные на методе декомпозиции, принципе сравнения и технике усреднения. Доказана теорема.
Теорема 11. Если V > 0, то положение равновесия (68) системы (77), (67) асимптотически устойчиво.
В работе [55] рассматривается та же задача, что и в [70], но с другими предположениями относительно возмущающих и управляющих моментов. Вначале рассматривается случай линейного восстанавливающего момента (V = 0) и линейного возмущающего момента Мр = Фх(£)(в1 — г;[) + Ф2(£)(«2 —М2), где Ф].(£), Ф2(£) € К3х3 — непрерывные и ограниченные при £ € [0, матрицы. Динамические уравнения имеют вид
Лш + ш х Лш = КБМ — аМ х г1 — а2^2 х М2 + Фх(Ь)М — п) + Ф2(£)М — М2), (78)
где К — положительный параметр. Доказана теорема.
Теорема 12. Пусть матрицы /^ФДт) ¿т (г = 1, 2) ограничены для £ ^ 0. Тогда существует число к0 такое, что положение равновесия (68) системы (78), (67) равномерно асимптотически устойчиво при любом К ^ К0.
г
Следует заметить, что условие ограниченности матриц §ФДт) ¿т (г = 1, 2) вы-
0
полняется, например, если элементы этих матриц являются периодическими с нулевыми средними значениями. Для более широкого класса возмущений доказана следующая теорема.
Теорема 13. Пусть матрицы Ф].(£) и Ф2(£) удовлетворяют условиям (73). Тогда существует число К0 такое, что положение равновесия (68) системы (78), (67) равномерно асимптотически устойчиво при любом К ^ К0.
Заметим, что условия (73) выполняются для матриц Фх(£) и Ф2(£), если их элементы являются почти периодическими с нулевыми средними значениями.
Далее рассматриваются уравнения (77) с существенно нелинейным восстанавливающим моментом (V > 0) и нелинейным возмущающим моментом того же порядка однородности. Установлено, что в этой ситуации асимптотическая устойчивость положения равновесия (68) может быть гарантирована при менее обременительных условиях, чем в рассмотренном выше линейном случае. Доказана теорема.
Теорема 14. Пусть матрица Ф(£) удовлетворяет условию (73). Тогда положение равновесия (68) системы (78), (67) равномерно асимптотически устойчиво при любом К ^ 0.
5-4- Управление с распределенным запаздыванием. В серии статей А. Ю. Александрова и А. А. Тихонова [71-74] изучаются задачи одноосной и трехосной стабилизации твердого тела при наличии распределенного запаздывания в системе управления.
В [71] рассматривается линейная механическая система с большим параметром при векторе скоростных сил и распределенным запаздыванием в позиционных силах. С помощью метода декомпозиции получены условия, при выполнении которых задача анализа устойчивости исходной системы дифференциальных уравнений второго порядка может быть сведена к исследованию устойчивости двух вспомогательных подсистем первого порядка. При этом одна из вспомогательных подсистем не содержит запаздывания, а для второй подсистемы, содержащей распределенное запаздывание, условия устойчивости формулируются в терминах разрешимости
систем линейных матричных неравенств. Для обоснования указанной декомпозиции используется прямой метод Ляпунова. Предложены специальные конструкции функционалов Ляпунова — Красовского. Разработанный подход применяется в задаче (60), (61) об одноосной стабилизации твердого тела в положении (62). Уравнения Эйлера имеют вид
JU + ш x Ju = -hD ш - a s(t) xr + bi s(£) xrd£, (79)
Jt-T
где D — постоянная положительно определенная матрица; h > 0, a > 0, b = const = 0; т — постоянное положительное запаздывание. Доказана теорема.
Теорема 15. Если выполнено неравенство a > т|b|, то существует число h0 > 0 такое, что при всех h ^ h0 положение равновесия (62) системы (60), (61) асимптотически устойчиво.
В работе [72] также решается задача об одноосной стабилизации твердого тела. При этом, в отличие от [71], используется управление по типу ПИД-регулятора, но с существенно нелинейной восстанавливающей компонентой и с распределенным запаздыванием (интегральный член) в отличие от классического ПИД-регулятора. Уравнения (60) управляемого движения тела имеют вид
JU + ш х JU = —D ш - a||s(t) - r\\»-1s(t) х r + bi ||s(£) - г|Г-1«Ш x rd£, (80)
t-T
где p ^ 1. Будем считать, что начальные функции х(0 для системы (61), (80) принадлежат пространству С([-т, 0], R6) непрерывных функций с равномерной нормой ||x||T = max£e[-T,0] Относительно условий асимптотической устойчивости по-
ложения равновесия (62) системы (61), (80) доказана следующая теорема.
Теорема 16. Пусть p > 0. Тогда при выполнении неравенства a > тЬ положение равновесия (62) системы (60), (61) асимптотически устойчиво.
Сравнение теорем 15 и 16 показывает, что при использовании управления по типу ПИД-регулятора, но с существенно нелинейной восстанавливающей компонентой и с распределенным запаздыванием одноосная стабилизация тела может быть обеспечена при менее жестких ограничениях на параметры системы. В частности, ослабление требований к величине диссипативного момента (не требуется наличия большого параметра в качестве множителя при диссипативном моменте) играет важную роль в задачах угловой стабилизации ИСЗ, поскольку создание диссипа-тивных моментов в космических условиях представляет непростую задачу. Кроме того, в отличие от линейного случая [71], из теоремы 16 следует, что при использовании существенно нелинейного управляющего момента одноосную стабилизацию тела можно гарантировать при a = 0, b < 0, т. е. когда на тело не действует восстанавливающий момент.
В работах [73, 74] решается задача (60), (67) о трехосной стабилизации твердого тела в положении равновесия (68) с использованием управления, построенного по типу ПИД-регулятора, но с распределенным запаздыванием (интегральный член). Уравнения (60) управляемого движения тела имеют вид [73]
Juu + ш х JU = -hD ш - aisi xri - a2S2 xr2 + b / (aisi(£) xri + a2&2(£) ХГ2) d£, (81)
t-T
где b = const. Пусть начальные функции х(£) для системы (67), (81) принадлежат пространству C([—г, 0], R12) непрерывных функций х(£) : [—r, 0] ^ R12 с равномерной нормой ||х||т = max5e[-r,0] ||х(£)||. Относительно условий асимптотической устойчивости положения равновесия (68) системы (67), (81) доказана следующая теорема.
Теорема 17. Пусть выполняется неравенство \b\r < 1. Тогда существует число h0 > 0 такое, что при всех h ^ h0 положение равновесия (68) системы (67), (81) асимптотически устойчиво.
В работе [74] тоже решается задача о трехосной стабилизации твердого тела с использованием управления по типу ПИД-регулятора с распределенным запаздыванием. При этом, в отличие от [73], впервые применяется нелинейный однородный восстанавливающий момент с распределенным запаздыванием. Кроме того, в отличие от подходов, применявшихся в [71-73], в работе [74] для получения условий асимптотической устойчивости вместо метода декомпозиции используется специальная конструкция функционала Ляпунова — Красовского. Уравнения (60) управляемого движения тела под действием диссипативного момента Add = — D Ш, где D — постоянная симметричная положительно определенная матрица, и восстанавливающего момента Air = —nv (Si, S2)(aisi xri + a2S2 xr2), где n(si, ¿2) = (ai ||si — ri||2 + a2||s2 — r2||2)/2, ai, a2 — положительные постоянные, v — неотрицательный параметр, имеют вид
ЛШ + ш x Лш = —D ш + Mr(t) + bi Mr(£) d£. (82)
Jt-r
Здесь b — постоянный коэффициент; r — постоянное положительное запаздывание.
Будем считать, что начальные функции х(£) для системы (67), (82) принадлежат пространству C([—r, 0], Ri2) непрерывных функций х(£) : [—r, 0] ^ Ri2 с равномерной нормой ЦхЦт = тах^-^о] ||х(01.
Задача нахождения условий асимптотической устойчивости положения равновесия (68) системы (67), (82) решена в работе [73] в случае линейного восстанавливающего момента (v = 0). Согласно теореме 17, стабилизацию положения равновесия можно гарантировать, если выполнено неравенство \b\r < 1, а перед вектором диссипативного момента имеется достаточно большой положительный множитель.
В работе [74] доказано, что при использовании существенно нелинейного восстанавливающего момента (v > 0) стабилизация имеет место при значительно менее жестких ограничениях на параметры системы, а именно справедлива теорема.
Теорема 18. Пусть v > 0. Если выполнено неравенство 1 + br > 0, то положение равновесия (68) системы (67), (82) асимптотически устойчиво.
В целом обращает на себя внимание простота условий асимптотической устойчивости в теоремах 15-18.
Практическая значимость доказанных теорем будет обсуждаться в следующей части обзора, а здесь приведем иллюстративный пример. Пусть твердое тело имеет главные центральные моменты инерции A = 5, B = 6, C = 4 (все численные значения параметров имеют размерности в единицах системы СИ).
Ставится задача стабилизации тела в положении (68), которому соответствуют нулевые значения «самолетных» углов ориентации тела в кениговой системе координат: ф = 0, ф = 0, в = 0. Управляющий момент берется в виде правой части
уравнений (82), где Б = diag(0.5, 0.5, 0.5), а\ = а2 = 2, Ь = 1.3. Рассматривается процесс стабилизации углового положения тела при разных значениях запаздывания т и показателя V степени нелинейности управления для одних и тех же начальных значений углов и угловых скоростей у(0) = 0.5, ф(0) = —0.2, 0(0) = 0.6, шж(0) = —0.055, шу (0) = —0.045, (0) = 0.05. В случае линейного управления (V = 0) при отсутствии распределенного запаздывания (т = 0) процесс стабилизации тела показан на рис. 1
Рис. 1. Процесс стабилизации при V = 0, т = 0.
Видно, что процесс стабилизации практически завершается при £ = 100. Если ввести в управление интегральную составляющую, выбрав т = 0.8, и снова рассмотреть поведение тела при тех же начальных условиях, то, как видно из рис. 2, стабилизация тела отсутствует.
„■-к : 11 и II 11 N. .. ----ч 1 \ \ \ \ И к/ |/ ,4 1 V .1 /I
\ \ 1 ' я ' I
20 ----1 ^ \ я \ ч \ 1 1 N А 'Л / /1 ' А | / Л1-'' н I / м I / ч А/ \
|— <г>(0 — у(<) — е(01
Рис. 2. Отсутствие стабилизации при V = 0, т = 0.8 (|с|т > 1).
Этот результат полностью соответствует теореме 17, гарантирующей асимптотическую устойчивость нулевого решения лишь при \Ь\т < 1. В данном случае \Ь\т = 1.3 • 0.8 = 1.04 > 1. Однако теорема 18 позволяет надеяться на то, что при нелинейном восстанавливающем моменте процесс стабилизации может состояться даже при Ь = 1.3 и т = 0.8, поскольку неравенство 1 + Ьт > 0 выполняется. Действительно, при V = 4 тело успешно стабилизируется, как показано на рис. 3.
При этом важно заметить, что, во-первых, тело в процессе стабилизации совершенно не испытывает колебаний, вызванных работой системы управления, и, во-вторых, процесс стабилизации завершается намного быстрее, чем в случае линейного управления без распределенного запаздывания (см. рис. 1). Наглядная демонстрация этих двух важных качеств предложенного нелинейного управления с распределенным запаздыванием свидетельствует о его перспективности для решения многочисленных практически важных задач, требующих не только быстродействия, но и соблюдения условий гладкости переходных режимов.
Рис. 3. Процесс стабилизации при v = 4, т = 0.8.
Из доказательства теоремы 18 следует также, что если восстанавливающий управляющий момент вовсе отсутствует, то диссипативный момент и интегральный член с запаздыванием позволяют обеспечить асимптотическую устойчивость стабилизируемого режима углового движения тела, причем наличие большого параметра при диссипативном моменте не обязательно. Таким образом, из сравнения процессов стабилизации при линейном и нелинейном восстанавливающем моментах (теоремы 17 и 18) следует, что если в линейном случае стабилизация тела происходит за счет восстанавливающего момента, а интегральная составляющая эффективно сглаживает переходные процессы, выполняя роль своего рода гасителя колебаний, вызванных работой системы управления, то при нелинейном восстанавливающем моменте интегральная составляющая не только стабилизирует тело, но и одновременно гасит его колебания.
5.5. Управление с неполной обратной связью. Пусть твердое тело вращается относительно центра масс С с угловой скоростью С. Кроме того, пусть в инерциальном пространстве имеется ортонормированный базис ¿1,М2,Мз, вращающийся относительно инерциального пространства с угловой скоростью с?1, а в теле имеется ортонормированный базис т\,г2,гз, вращающийся относительно тела с угловой скоростью с?2, где С и с2 — непрерывные и ограниченные вектор-функции. Тогда в системе Схуг главных центральных осей инерции тела движение тела и изменение ортов описывается дифференциальными уравнениями Эйлера — Пуассона:
ЛсС + С х Лсс = М + Мр, (83)
= (сс1 - С) х п = сс2 х г., г = 1,2,3, (84)
где М — управляющий момент; Мр — возмущающий момент.
Практический интерес представляют две основные задачи.
Задача 1. Выбрать управляющий момент М так, чтобы либо в. = г., либо в. ^ г.\ при £ ^то (г = 1, 2, 3).
Задача 2. Выбрать управляющий момент М так, чтобы обеспечить асимптотическую устойчивость вращательного движения тела вокруг одной из его главных центральных осей инерции с заданной угловой скоростью сс = мо.
При наличии полной обратной связи, когда имеется три датчика угловой скорости, измеряющих ее проекции на выбранные оси, обе эти задачи решены В.И. Зубовым [40]. В то же время А. М. Летов поставил задачу построения такого управления с использованием меньшего числа датчиков [75]. Задача Летова, продолжая
и развивая результаты В. И. Зубова, рассматривалась в ряде работ ученых СПбГУ. Некоторые из них упомянуты ниже.
1. В работе [76] рассматриваются задачи 1, 2 в предположении о том, что измеряются не все фазовые переменные тела, а именно компоненты ортов si, s2, ¿3 в системе координат Cxyz и проекции вектора ш на два неколлинеарные направления с ортами 61 и 62, которые, не ограничивая общности, можно считать ортогональными. С использованием обозначений ш* = — Ш2, 63 = 61 х 62, Bi = (61,62,63), B2 = Bi diag(1,1, 0), w = BJoj, w* = Bjш* управляющий момент строится в виде
3
M = лш* + ш* х Лш* + ш* х JB2 (w — w*) — B1D1(w — w*) + ^ aM х s, — Mp, (85)
i=i
где D = (dkl), du > 0, dud22 — №1 + d^)2/4 > 0, dk3 = 0, d3m = w*(6m х J63), a1 > 0, a,2 > 0, аз ^ 0 (k,l = 1, 2, 3; m =1, 2). Тогда справедлива теорема.
Теорема 19. Пусть v > 0. Если на программном движении задачи 1 выполнено неравенство M*(Л63 х 63) > 0, то это движение при управляющем моменте (85) асимптотически устойчиво по Ляпунову. Если ш* = 0 и соотношения Ь3 х Л63 = 0, 63 х r3 =0 одновременно не выполняются, то программное движение задачи 1 при управляющем моменте (85) асимптотически устойчиво по Ляпунову; если же ш* = 63 х Л63 = 63 х r3 = 0, то это движение лишь устойчиво.
Доказательство теоремы 19 проводится путем исследования поведения функции V = (ш — ш*)Л(ш — ш*) + 3=1 — ri)(S1 — r,) на решениях системы (83)—(85).
Аналогичным образом в [76] получены и доказаны достаточные условия, при выполнении которых построенные управляющие моменты осуществляют программную ориентацию тела в задаче 2.
2. При ш0 = 0 задача 1 решена в [77], а задача 2 — в [78].
3. При ш0 = 0 задача 2 решалась в локальной постановке в работах [79, 80] для случая, когда вектор шо направлен вдоль главной оси инерции тела.
4. В случае двух датчиков задача 2 решена в [81] для почти всех тензоров инерции.
5. В случае одного датчика задача 2 решена в [82] в линейном приближении с использованием динамического регулятора.
6. Задача 2 в случае одного датчика нелокально исследовалась в [83] с учетом нелинейных членов. Исследование основано на методе функций Ляпунова.
7. В [84] исследована задача 2. Искомые управления предлагается выбирать как линейные функции относительно отклонений проекций векторов угловых скоростей от заданных значений. Получены условия, при которых такое управление обеспечивает асимптотическую устойчивость вращения вокруг большой, малой или средней оси инерции тела.
6. Заключение. Последние 50 лет истории Университета отличаются от первых 250 лет взрывным ростом количества публикаций по динамике твердого тела. Причина этого явления — быстрое расширение круга актуальных задач, для которых методы динамики твердого тела, разработанные ранее и развитые в рассматриваемый период, являются надежной теоретической базой. Рассмотрение некоторых из этих методов, развитых учеными Университета, дает представление о достигнутых в рассматриваемый период достижениях, относящихся как к неуправляемому,
так и к управляемому движениям твердого тела. О достижениях ученых Университета в области прикладных исследований, опирающихся на динамику твердого тела, в том числе и о достижениях в области динамики космических аппаратов, будет рассказано в последующих частях данной серии статей.
Литература
1. Тихонов А. А. Динамика твердого тела от уравнений Эйлера до управления угловым движением ИСЗ в трудах ученых СПбГУ. Ч. 1. Вестник Санкт-Петербургского университета.. Математика.. Механика.. Астрономия 10 (68), вып. 3, 457—486 (2023). https://doi.org/10.21638/spbu01.2023.303
2. Новоселов В. С. Приведение дифференциальных уравнений механики к каноническому виду. Вестник Ленинградского университета. Сер. 1. Математика. Механика.. Астрономия 3, 100— 105 (1972).
3. Суслов Г. К. О кинетическом потенциале Гельмгольца. Математический сб. 19, 1—14 (1896).
4. Новоселов В. С. Аналитическая механика систем с переменными массами. Ленинград, Изд-во Ленингр. ун-та (1969).
5. Новоселов В. С. Об исследовании нелинейных колебаний методами гамильтоновой механики. Вестник Ленинградского университета. Сер. 1. Математика. Механика. Астрономия 4, 123-125 (1976).
6. Leitmann G. Some remarks on Hamilton's principle. Trans. ASME E 85, 623 (1963).
7. Jones S. A class of minimal variational principles for some nonconservative systems. Trans. ASME E 40, 1144-1146 (1973).
8. Stuelpnagel J. On the parametrization of the three-dimensional rotation group. SIAM Review 6 (4), 422-430 (1964).
9. Кирпичников С. Н., Новоселов В. С. Математические аспекты кинематики твердого тела. Ленинград, Изд-во Ленингр. ун-та (1986).
10. Поляхов Н. Н., Зегжда С. А., Юшков М.П. Теоретическая механика. Москва, Юрайт (2015).
11. Поляхов Н. Н., Зегжда С. А., Юшков М. П. Специальная форма уравнений динамики системы твердых тел. Доклады Академии наук СССР 309 (4), 805-807 (1989).
12. Белецкий В. В. Движение спутника относительно центра масс в гравитационном поле. Москва, Изд-во Моск. ун-та (1975).
13. Breakwell J.V., Pringle R.Jr. Non-linear resonances affecting gravity-gradient stability. Astrodynamics, 305-325. Paris, Gauthier-Villars (1966).
14. Тихонов А. А. Об одной форме дифференциальных уравнений возмущенного движения гравитационно-ориентированного твердого тела. Вестник Ленинградского университета. Сер. 1. Математика. Механика. Астрономия 1 (1), 71-75 (1990).
15. Тихонов А. А. О нелинейных резонансных колебаниях заряженного твердого тела в условиях нестационарного магнитного поля. Ч. 1. Вестник Санкт-Петербургского университета. Сер. 1. Математика. Механика. Астрономия 3 (15), 60-65 (1992).
16. Тихонов А. А. О нелинейных резонансных колебаниях заряженного твердого тела в условиях нестационарного магнитного поля. Ч. 2. Вестник Санкт-Петербургского университета. Сер. 1. Математика. Механика. Астрономия 4 (22), 86-91 (1992).
17. Тихонов А. А. Резонансные явления в колебаниях гравитационно-ориентированного твердого тела. Ч. 1. Спектральная структура возмущающего момента. Вестник Санкт-Петербургского университета. Сер. 1. Математика. Механика. Астрономия 2 (8), 82-91 (1997).
18. Тихонов А. А. Резонансные явления в колебаниях гравитационно-ориентированного твердого тела. Ч. 2. Реализуемые резонансы. Вестник Санкт-Петербургского университета. Сер. 1. Математика. Механика. Астрономия 3 (15), 85-94 (1997).
19. Тихонов А. А. Резонансные явления в колебаниях гравитационно-ориентированного твердого тела. Ч. 3. Усредненные уравнения. Одночастотные резонансы. Вестник Санкт-Петербургского университета. Сер. 1. Математика. Механика. Астрономия 3 (15), 116-126 (1998).
20. Тихонов А. А. Резонансные явления в колебаниях гравитационно-ориентированного твердого тела. Ч. 4. Многочастотные резонансы. Вестник Санкт-Петербургского университета. Сер. 1. Математика. Механика. Астрономия 1 (1), 131-137 (2000).
21. Kosjakov E. A., Tikhonov A.A. Differential equations for librational motion of gravity-oriented rigid body. International Journal of Non-Linear Mechanics 73, 51-57 (2015).
22. Куницын А. Л. О резонансной стабилизации одного класса неустойчивых систем. Прикладная математика и механика 75 (5), 727—730 (2011).
23. Kosjakov E. A., Partalin T. A., Tikhonov A. A. On Nonlinear Resonances in Satellite Oscillations. Scientific publications of the State University of Novi Pazar. Ser. A: Appl. math., inform. and mech. 7 (1), 19-24 (2015).
24. Кошляков В. Н. О применении параметров Родрига — Гамильтона и Кейли — Клейна в прикладной теории гироскопов. Прикладаня математика и механика 29 (4), 729-733 (1965).
25. Белецкий В. В., Голубков В. В., Лавровский Э.К., Трушин С. И., Хацкевич И. Г. Определение ориентации и вращения искусственных спутников по данным измерений. Космические исследования 5 (5), 686-702 (1967).
26. Бранец В. Н., Шмыглевский И. П. Применение кватернионов в задачах ориентации твердого тела. Москва, Наука (1973).
27. Петров К. Г., Тихонов А. А. Уравнения ротационного движения твердого тела, основанные на использовании кватернионных параметров. Известия РАН. Механика твердого тела 3, 3-16 (2002).
28. Тихонов А. А. О вековой эволюции ротационного движения заряженного ИСЗ на регрессирующей орбите. Космические исследования 43 (2), 111-125 (2005).
29. Белецкий В. В. Движение искусственного спутника относительно центра масс. Москва, Наука (1965).
30. Marandi S.R., Modi V.J. A preffered coordinate system and the associated orientation representation in attitude dynamics. Acta Astronaut. 15 (11), 833-843 (1987).
31. Shuster M. D. A survey of attitude representations. Journal of the astronautical sciences 41 (4), 439-517 (1993).
32. Schaub H., Junkins J. L. Stereographic orientation parameters for attitude dynamics: a generalization of the Rodrigues parameters. The Journal of the Astronautical Sciences 44 (1), 1-19 (1996).
33. Dolicanin C. B., Tikhonov A. A. On dynamical equations in s-parameters for rigid body attitude motion. International Conference on Mechanics-Seventh Polyakhov's Reading. February 2-6, 2015, St. Petersburg. https://doi.org/10.1109/P0LYAKH0V.2015.7106723
34. Tikhonov A. A., Dolicanin C., Partalin T. A., Arandjelovic I. A new form of equations for rigid body rotational dynamics. Technical Gazette 21 (6), 1221-1227 (2014).
35. Кирпичников С. Н. Об устойчивости стационарных движений двух гравитирущих тел: осе-симметричного и сферически-симметричного. Некомпланарный случай. Астрономический журнал 66 (3), 612-621 (1989).
36. Кирпичников С. Н., Кокорев A. A. О первом типе стационарных движений двух гравити-рующих тел: осесимметричного и сферически-симметричного. Вестник Ленинградского университета. Сер. 1. Математика. Механика. Астрономия 2 (8), 115-123 (1991).
37. Кирпичников С. Н. О втором типе стационарных движений двух притягивающихся тел: осесимметричного и сферически-симметричного. Вестник Ленинградского университета. Сер. 1. Математика. Механика. Астрономия 3 (15), 92-97 (1991).
38. Зубов В. И., Ермолин В. С., Сергеев С. Л., Смирнов Е. Я. Управление вращательным движением твердого тела. Ленинград, Изд-во Ленингр. ун-та (1978).
39. Зубов В. И. Консервативные численные методы интегрирования дифференциальных уравнений в нелинейной механике. Доклады РАН 354 (4), 446-448 (1997).
40. Зубов В. И. Об активном управлении вращательным движением твердого тела. Диффе-ренц. уравнения 6 (11), 2086-2087 (1970).
41. Зубов В. И. Об управлении вращательным движением. Дифференц. уравнения 7 (7), 13201322 (1971).
42. Зубов В. И. Управляемое вращение и динамика относительного движения. Доклады Академии наук СССР 219 (3), 565-566 (1974).
43. Зубов В. И. Аналитическая конструкция закона регулирования, управляющего вращательным движением твердого тела при непрямом регулировании. Дифференц. уравнения 10 (10), 18981899 (1974).
44. Зубов В. И. Лекции по теории управления. Москва, Наука (1975).
45. Зубов В. И. Оптимальное демпфирование вращательного движения. Прикладная механика 16 (5), 113-121 (1980).
46. Зубов В. И. Динамика управляемых систем. Москва, Высшая школа (1982).
47. Зубов В. И., Ермолин В. С., Иголкин В. Н. Чернецкий В. И. Динамика свободного твердого тела и определение его ориентации в пространстве. Чернецкий В. И. (ред.). Ленинград, Изд-во Ленингр. ун-та (1968).
48. Смирнов Е. Я. Некоторые задачи математической теории управления. Ленинград, Ле-низдат (1981).
49. Смирнов Е. Я. Управление вращательным движением твердого тела в условиях неопределенности. Дифференциальные уравнения 11 (5), 928—929 (1975).
50. Смирнов Е. Я., Павликов В.Ю., Щербаков П. П., Юрков А. В. Управление движением механических систем. Ленинград, Изд-во Ленингр. ун-та (1985).
51. Александров А.Ю. Об устойчивости положений равновесия нелинейных неавтономных механических систем. Прикладная математика и механика 71 (3), 361—376 (2007).
52. Aleksandrov A. Yu., Aleksandrova E. B., Tikhonov A. A. Monoaxial attitude stabilization of a rigid body under vanishing restoring torque. Nonlinear Dynamics and Systems Theory 18 (1), 12—21 (2018).
53. Aleksandrov A. Yu., Tikhonov A. A. Rigid body stabilization under time-varying perturbations with zero mean values. Cybernetics and Physics 7 (1), 5—10 (2018).
54. Aleksandrov A. Yu., Tikhonov A. A. Attitude stabilization of a rigid body under the action of a vanishing control torque. Nonlinear Dynamics 93 (2), 285—293 (2018).
55. Aleksandrov A. Yu., Tikhonov A. A. Attitude stabilization of a rigid body under disturbances with zero mean values. Acta Mechanica 233, 1231-1242 (2022).
56. Александров А.Ю., Тихонов А. А. Стабилизация вращательного движения твердого тела в условиях убывающей диссипации. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 4 (62), вып. 3, 631-640 (2017). https://doi.org/10.21638/11701 /spbu01.2017.411
57. Aleksandrov A.Yu., Tikhonov A. A. On the Triaxial Stabilization of a Rigid Body under an Evolution of Dissipative Torque. Proceedings of 2018 14th International Conference Stability and Oscillations of Nonlinear Control Systems (Pyatnitskiy Conference), STAB 2018 (July, 6), 1-4 (2018).
58. Aleksandrov A.Yu., Aleksandrova E. B., Tikhonov A. A. Attitude Stabilization of a Rigid Body under Vanishing Restoring Torque. Advances in Stability and Control Theory for Uncertain Dynamical Systems. Cruz-Hernandez C., Martynyuk A. A., Mazko A.G. ^ds). Stability, Oscillations and Optimization of Systems 11, 213-228 (2021).
59. Rouche N., Habets P., Laloy M. Stability theory by Liapunov's direct method. New York, Springer (1977).
60. Хатвани Л. О действии демпфирования на свойства устойчивости равновесий неавтономных систем. Прикладная математика и механика 65 (4), 725-732 (2001).
61. Aleksandrov A.Yu., Aleksandrova E. B., Tikhonov A. A. On the Monoaxial Stabilization of a Rigid Body under Vanishing Restoring Torque. AIP Conference Proceedings 1959, 080001 (2018). https://doi.org/10.1063/1.5034718
62. Александров А. Ю., Косов А. А., Чэнь Я. Об устойчивости и стабилизации механических систем с переключениями. Автоматика и телемеханика 6, 5-17 (2011).
63. Aleksandrov A.Yu., Aleksandrova E. B. Asymptotic stability conditions for a class of hybrid mechanical systems with switched nonlinear positional forces. Nonlinear Dynamics 83 (4), 2427-2434 (2016).
64. Александров А. Ю. Об устойчивости равновесия нестационарных систем. Прикладная математика и механика 60 (2), 205-209 (1996).
65. Александров А.Ю. Об устойчивости сложных систем в критических случаях. Автоматика и телемеханика 9, 3-13 (2001).
66.Антипов К. А., Тихонов А. А. Электродинамическое управление в задаче о стабилизации космического аппарата в геомагнитном поле. Космические исследования 52 (6), 512-520 (2014).
67. Kosov A. A. The exponential stability and stabilization of non-autonomous mechanical systems with non-conservative forces. J. Appl. Math. Mech. 71 (3), 371-384 (2007).
68. Александров А. Ю., Тихонов А. А. Одноосная стабилизация вращательного движения твердого тела при наличии возмущений с нулевыми средними значениям. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 6 (64), вып. 2, 270-280 (2019). https://doi.org/10.21638/11701/spbu01.2019.209
69. Зубов В. И. Устойчивость движения. Москва, Высшая школа (1973).
70. Aleksandrov A. Yu., Tikhonov A. A. Nonlinear Control for Attitude Stabilization of a Rigid Body Forced by Nonstationary Disturbances with Zero Mean Values. Journal of Applied and Computational Mechanics 7 (2), 790-797 (2021).
71. Александров А. Ю., Тихонов А. А. Анализ устойчивости механических систем с распределенным запаздыванием на основе декомпозиции. Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления 17 (1), 13-26 (2021).
72. Александров А. Ю., Тихонов А. А. Анализ устойчивости механических систем с существенно нелинейными позиционными силами при наличии распределенного запаздывания. Автоматика и телемеханика 1, 3-22 (2023).
73. Aleksandrov A. Yu., Tikhonov A. A. On the attitude stabilization of a rigid body under control with distributed delay. Mechanics Based Design of Structures and Machines 51 (4), 2241-2250 (2023).
74. Александров А. Ю., Тихонов А. А. Нелинейное управление с распределенным запаздыванием для угловой стабилизации твердого тела. Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия 9 (67), вып. 4, 653-664 (2022). https://doi.org/10.21638/spbu01.2022.408
75. Летов А.М. Динамика полета и управление. Москва, Наука (1969).
76. Смирнов Е. Я. Управление вращательным движением твердого тела. Дифференц. уравнения 8 (12), 2264-2266 (1972).
77. Мейлахс А. М., Хитров Г. М. О гашении вращения твердого тела. Дифференц. уравнения 12 (8), 1520-1522 (1976).
78. Мейлахс А. М. Об ориентации инертного твердого тела. Дифференц. уравнения 12 (11), 2101-2120 (1976).
79. Лебедев Д. В. Об управлении движением твердого тела при неполной информации об угловой скорости. Автоматика и телемеханика 12, 5-11 (1978).
80. Ермолин В. С. Стабилизация вращательного движения динамически симметричного твердого тела. Управление, надежность и навигация 3, 132-137 (1976).
81. Мейлахс А.М. О стабилизации вращательного движения твердого тела. Дифференц. уравнения 16 (10), 1888-1890 (1980).
82. Лебедев Д. В. Об управлении вращательным движением твердого тела при неполной информации о векторе угловой скорости. Прикладная математика и механика 41 (2), 219-224 (1977).
83. Мейлахс А. М. Об управлении вращательным движением твердого тела. В: Математические методы исследования управляемых механических систем. Ленинград, Изд-во Ленингр. ун-та (1982).
84. Джумагалиева М. Б. Об асимптотической устойчивости вращательного движения твердого тела в случае Эйлера. В: Проблемы механики управляемого движения. Нелинейные динамические системы.. Пермь, Пермский университет (1984).
Статья поступила в редакцию 2 октября 2023 г.;
доработана 4 ноября 2023 г.; рекомендована к печати 9 ноября 2023 г.
Контактная информация:
Тихонов Алексей Александрович — д-р физ.-мат. наук, проф.; a.tikhonov@spbu.ru
Rigid body dynamics from the Euler equations to the spacecraft attitude control in the works of scientists from Saint Petersburg State University. Part 2*
A. A. Tikhonov
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Tikhonov A. A. Rigid body dynamics from the Euler equations to the spacecraft attitude control in the works of scientists from Saint Petersburg State University. Part 2. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11 (69), issue 2, pp. 259-302. https://doi.org/10.21638/spbu01.2024.203 (In Russian)
*See first part: Tikhonov A. A. Rigid body dynamics from the Euler equations to the spacecraft attitude control in the works of scientists from Saint Petersburg State University. Part 1. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2023, vol. 10(68), issue 3, pp. 457— 487. https://doi.org/10.21638/spbu01.2023.303 (In Russian)
This article represents the second part of a survey dedicated to the 300th anniversary of Saint Petersburg State University (SPbSU). The article is intended on analysis the scientific achievements of the Saint Petersburg School of Mathematics and Mechanics in the field of rigid body dynamics in the last 50 years. The main attention in the article is paid to theoretical research in the field of uncontrolled and controlled motions of a rigid body. Keywords: rigid body, spacecraft, artificial Earth satellite, dynamics, attitude motion, attitude control, stabilization, oscillations.
References
1. Tikhonov A.A. Rigid-Body Dynamics from the Euler Equations to the Attitude Control of Spacecraft in the Works of Scientists from Saint Petersburg State University. Part 1. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 10 (68), iss. 3, 457—486 (2023). https://doi.org/10.21638/spbu01.2023.303 (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mathematics 56, iss.3, 322-340 (2023). https://doi.org/10.1134/S1063454123030081].
2. Novoselov V. S. Reducing differential equations of mechanics to canonical form. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 3, 100-105 (1972). (In Russian)
3. Suslov G.K. O kineticheskom potentsiale Helmholtsa [On the Helmholtz kinetic potential]. Matematicheskii sbornik — Mathematical collection 19 (1), 197-210 (1896). (In Russian)
4. Novoselov V. S. Analytical mechanics of systems with variable masses. Leningrad, Leningrad University Press (1969). (In Russian)
5. Novoselov V. S. On the study of nonlinear oscillations by methods of Hamiltonian mechanics. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 4, 123-125 (1976). (In Russian)
6. Leitmann G. Some remarks on Hamilton's principle. Trans. ASME E 85, 623 (1963).
7. Jones S. A class of minimal variational principles for some nonconservative systems. Trans. ASME E 40, 1144-1146 (1973).
8. Stuelpnagel J. On the parametrization of the three-dimensional rotation group. SIAM Review 6 (4), 422-430 (1964).
9. Kirpichnikov S.N., Novoselov V. S. Mathematical aspects of rigid body kinematics. Leningrad, Leningrad University Press (1986). (In Russian)
10. Polyakhov N.N., Zegzhda S.A., Yushkov M. P. Theoretical mechanics. Moscow, Iurait Publ. (2015). (In Russian)
11. Polyakhov N.N., Zegzhda S.A., Yushkov M.P. A Special Form of Dynamical Equations of a Rigid Bodies System. Doklady Akademii nauk SSSR 309 (4), 805-807 (1989).
12. Beletsky V. V. Motion of a satellite about its center of mass in field of gravity. Moscow, Moscow University Press (1975). (In Russian)
13. Breakwell J.V., Pringle R.Jr. Non-linear resonances affecting gravity-gradient stability. Astrodynamics, 305-325. Paris, Gauthier-Villars (1966).
14. Tikhonov A. A. A form of differential equations for disturbed motion of gravity-oriented rigid body. Vestnik of Leningrad University. Ser. 1. Mathematics. Mechanics. Astronomy 1 (1), 71-75 (1990). (In Russian) [Eng. transl.: Vestnik Leningrad University. Mechanics. Bulletin 1, 15-18 (1990)].
15. Tikhonov A. A. Nonlinear Resonance Vibrations of a Charged Rigid Body in a Nonstationary Magnetic Field. Part 1. Vestnik of Saint Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy 3 (15), 60-65 (1992). (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mechanics. Bulletin 3, 12-19 (1992)].
16. Tikhonov A. A. Nonlinear Resonance Vibrations of a Charged Rigid Body in a Nonstationary Magnetic Field. Part 2. Vestnik of Saint Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy 4 (22), 86-91 (1992). (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mechanics. Bulletin 4, 28-34 (1992)].
17. Tikhonov A. A. Resonance Phenomena in Oscillations of a Gravity-Oriented Rigid Body. Part I. The Spectral Structure of a Perturbation Moment. Vestnik of Saint Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy 2 (8), 82-91 (1997). (1992). (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mechanics. Bulletin 2, 15-23 (1997)].
18. Tikhonov A. A. Resonance Phenomena in Oscillations of a Gravity-Oriented Rigid Body. Part II. The Realizable Resonances. Vestnik of Saint Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy 3 (15), 85-94 (1997). (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mechanics. Bulletin 3, 14-22 (1997)].
19. Tikhonov A.A. Resonance Phenomena in Oscillations of a Gravity-Oriented Rigid Body. Part III. The Averaged Equations. Single-Frequency Resonances. Vestnik of Saint Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy 3 (15), 116—126 (1998). (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mechanics. Bulletin 3, 42—51 (1998)].
20. Tikhonov A.A. Resonance phenomena in oscillations of a gravity-oriented rigid body. Part 4: multifrequency resonances. Vestnik of Saint Petersburg University. Ser. 1. Mathematics. Mechanics. Astronomy 1 (1), 131—137 (2000). (In Russian) [Eng. transl.: Vestnik St. Petersburg University. Mechanics. Bulletin 1, 131-137 (2000)].
21. Kosjakov E. A., Tikhonov A.A. Differential equations for librational motion of gravity-oriented rigid body. International Journal of Non-Linear Mechanics 73, 51-57 (2015).
22. Kunitsyn A. L. The resonance stabilization of a class of unstable systems. Journal of Applied Mathematics and Mechanics 75 (5), 508-510 (2011). (In Russian)
23. Kosjakov E. A., Partalin T. A., Tikhonov A.A. On Nonlinear Resonances in Satellite Oscillations. Scientific publications of the State University of Novi Pazar. Ser. A: Appl. math., inform. and mech. 7 (1), 19-24 (2015).
24. Koshlyakov V.N. On the application of the Rodriguez — Hamilton and Cayley — Klein parameters in the applied gyroscope theory. Applied mathematics and mechanics 29 (4), 729-733 (1965). (In Russian)
25. Beletsky V. V., Golubkov V. V., Lavrovskii E. K., Trushin S. I., Khatskevich I. G. Determination of the orientation and rotation of artificial satellites from measurement data. Kosmicheskie issledovaniia
5 (5), 686-702 (1967). (In Russian) [Eng. transl.: Cosmic Research 5 (5), 585-598 (1967)].
26. Branets V. N., Shmyglevsky I. P. Application of Quaternions to the Problems of Rigid Body Orientation. Moscow, Nauka Publ. (1973). (In Russian)
27. Petrov K. G., Tikhonov A. A. Equations of Rotational Motion of a Solid Body Based on the Use of Quaternion Parameters. Izvestiia RAN. Mekhanika tverdogo tela 3, 3-16 (2002). (In Russian)
28. Tikhonov A. A. Secular Evolution of Rotary Motion of a Charged Satellite in a Decaying Orbit. Kosmicheskie issledovaniia 43 (2), 111-125 (2005). (In Russian) [Eng. transl.: Cosmic Research 43 (2), 107-121 (2005)].
29. Beletskii V. V. Motion of an Artificial Satellite with Respect to the Center of Mass. Moscow, Nauka Publ. (1965). (In Russian)
30. Marandi S.R., Modi V.J. A preffered coordinate system and the associated orientation representation in attitude dynamics. Acta Astronaut 15 (11), 833-843 (1987).
31. Shuster M. D. A survey of attitude representations. Journal of the astronautical sciences 41 (4), 439-517 (1993).
32. Schaub H., Junkins J. L. Stereographic orientation parameters for attitude dynamics: a generalization of the Rodrigues parameters. The Journal of the Astronautical Sciences 44 (1), 119 (1996).
33. Dolicanin C. B., Tikhonov A. A. On dynamical equations in s-parameters for rigid body attitude motion. International Conference on Mechanics-Seventh Polyakhov's Reading. February 2-6, 2015, St. Petersburg. https://doi.org/10.1109/POLYAKHOV.2015.7106723
34. Tikhonov A. A., Dolicanin C., Partalin T. A., Arandjelovic I. A new form of equations for rigid body rotational dynamics. Technical Gazette 21 (6), 1221-1227 (2014).
35. Kirpichnikov S. N. The stability of the steady-state motions of two gravitating bodies: Axisymmetric and spherically symetric. The noncomplanar case. Astronomicheskii zhurnal 66, 612-621 (1989). (In Russian)
36. Kirpichnikov S. N., Kokorev A. A. The first type of steady motions of two gravitating bodies: Axisymmetric and spherically symmetric. Vestnik Leningrad University Mathematics. Mechanics. Astronomy 2 (8), 115-123 (1991). (In Russian)
37. Kirpichnikov S. N. The second kind of steady motions of two gravitating bodies — Axisymmetric and spherically symmetric. Vestnik Leningrad University Mathematics. Mechanics. Astronomy 3 (15), 92-97 (1991). (In Russian)
38. Zubov V. I., Ermolin V. S., Sergeev S. L., Smirnov E. Ya. Control of rotational motion of a rigid, body. Aleshkov Yu. Z. (ed.). Leningrad University Press (1978). (In Russian)
39. Zubov V. I. Conservative numerical methods for integrating differential equations in nonlinear mechanics. Doklady RAN 354 (4), 446-448 (1997). (In Russian)
40. Zubov V. I. On the active control of the rotational motion of a rigid body. Differential Equations
6 (11), 2086-2087 (1970). (In Russian)
41. Zubov V.I. About rotary motion control. Differential Equations 7 (7), 1320-1322 (1971). (In Russian)
42. Zubov V.I. Controlled rotation and relative motion dynamics. Doklady AN USSR 219 (3), 565-566 (1974). (In Russian)
43. Zubov V. I. Analytical design of the control law governing the rotational motion of a rigid body under indirect control. Differential equations 10 (10), 1898-1899 (1974). (In Russian)
44. Zubov V.I. Lectures on control theory. Moscow, Nauka Publ. (1975). (In Russian)
45. Zubov V.I. Optimal damping of rotational motion. Applied mechanics 16 (5), 113-121 (1980). (In Russian)
46. Zubov V.I. Dynamics of Controlled Systems. Moscow, Vysshaya shkola Publ. (1982). (In Russian)
47. Zubov V. I., Ermolin V. S., Igolkin V.N. Chernetsky V.I. Dynamics of a free rigid body and determination of its orientation in space. Chernetsky V.I. (ed.). Leningrad, Leningrad University Press (1968). (In Russian)
48. Smirnov E. Ya. Some Problems of Mathematical Control Theory Leningrad, Lenizdat Publ. (1981). (In Russian)
49. Smirnov E. Ya. Control of rotational motion of a rigid body under conditions of uncertainty. Differential equations 11 (5), 928-929 (1975). (In Russian)
50. Smirnov E. Ya., Pavlikov V.Yu., Shcherbakov P.P., Yurkov A.V. Control of Motion of Mechanical Systems. Leningrad, Leningrad University Press (1985). (In Russian)
51. Aleksandrov A.Yu. The stability of the equilibrium positions of non-linear non-autonomous mechanical systems. Journal of Applied Mathematics and Mechanics 71 (3), 324-338 (2007).
52. Aleksandrov A.Yu., Aleksandrova E. B., Tikhonov A.A. Monoaxial attitude stabilization of a rigid body under vanishing restoring torque. Nonlinear Dynamics and Systems Theory 18 (1), 12-21 (2018).
53. Aleksandrov A.Yu., Tikhonov A.A. Rigid body stabilization under time-varying perturbations with zero mean values. Cybernetics and Physics 7 (1), 5-10 (2018).
54. Aleksandrov A.Yu., Tikhonov A.A. Attitude stabilization of a rigid body under the action of a vanishing control torque. Nonlinear Dynamics 93 (2), 285-293 (2018).
55. Aleksandrov A. Yu., Tikhonov A. A. Attitude stabilization of a rigid body under disturbances with zero mean values. Acta Mechanica 233, 1231-1242 (2022).
56. Aleksandrov A. Y., Tikhonov A. A. Attitude stabilization of a rigid body in conditions of decreasing dissipation. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 4 (62), iss.3, 631-640 (2017). https://doi.org/10.21638/11701/spbu01.2017.411 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 50, 384-391, iss.3 (2017). https://doi.org/10.3103/S1063454117040021].
57. Aleksandrov A.Yu., Tikhonov A.A. On the Triaxial Stabilization of a Rigid Body under an Evolution of Dissipative Torque. Proceedings of 2018 14th International Conference Stability and, Oscillations of Nonlinear Control Systems (Pyatnitskiy Conference), STAB 2018 (July, 6), 1-4 (2018).
58. Aleksandrov A.Yu., Aleksandrova E. B., Tikhonov A.A. Attitude Stabilization of a Rigid Body under Vanishing Restoring Torque. Advances in Stability and Control Theory for Uncertain Dynamical Systems. Cruz-Hernandez C., Martynyuk A.A., Mazko A.G. (eds). Stability, Oscillations and Optimization of Systems 11, 213-228 (2021).
59. Rouche N., Habets P., Laloy M. Stability theory by Liapunov's direct method. New York, Springer (1977).
60. Hatvani L. The effect of damping on the stability properties of equilibria of non-autonomous systems. J. Appl. Math. Mech. 65 (4), 707-713 (2001).
61. Aleksandrov A.Yu., Aleksandrova E. B., Tikhonov A.A. On the Monoaxial Stabilization of a Rigid Body under Vanishing Restoring Torque. AIP Conference Proceedings 1959, 080001 (2018). https://doi.org/10.1063/1.5034718
62. Aleksandrov A.Yu., Kosov A.A., Chen Y. Stability and stabilization of mechanical systems with switching. Automation and Remote Control 72 (6), 1143-1154 (2011).
63. Aleksandrov A.Yu., Aleksandrova E. B. Asymptotic stability conditions for a class of hybrid mechanical systems with switched nonlinear positional forces. Nonlinear Dynamics 83 (4), 2427-2434 (2016).
64. Aleksandrov A. Yu. The stability of equilibrium of non-stationary systems. Prikladnaia matematika i mekhanika 60 (2), 205-209 (1996). (In Russian) [Eng. transl.: Journal of Applied, Mathematics and Mechanics 60 (2), 199-203. https://doi.org/10.1016/0021-8928(96)00026-3].
65. Aleksandrov A. Y. Stability of Complex Systems in Critical Cases. Avtomatika i telemekhanika 9, 3-13 (2001). (In Russian) [Eng. transl.: Automation and Remote Control 62, 1397-1406 (2001). https://doi.org/10.1023/A:1011676123661].
66. Antipov K. A., Tikhonov A. A. Electrodynamic control for spacecraft attitude stability in the geomagnetic field. Kosmicheskie issledovaniia 52 (6), 512-520 (2014). (In Russian) [Eng. transl.: Cosmic Res. 52 (6), 472-480 (2014). https://doi.org/10.1134/S001095251406001X].
67. Kosov A. A. The exponential stability and stabilization of non-autonomous mechanical systems with non-conservative forces. J. Appl. Math. Mech. 71 (3), 371-384 (2007).
68. Aleksandrov A. Y., Tikhonov A. A. Uniaxial Attitude Stabilization of a Rigid Body under Conditions of Nonstationary Perturbations with Zero Mean Values. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 6 (64), iss. 2, 270-280 (2019). https://doi.org/10.21638/11701/spbu01.2019.209 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 52, iss. 2, 187-193 (2019). https://doi.org/10.1134/S106345411902002X].
69. Zubov V.I. Stability of Motion. Moscow, Vysshaia shkola Publ. (1973). (In Russian)
70. Aleksandrov A. Yu., Tikhonov A. A. Nonlinear Control for Attitude Stabilization of a Rigid Body Forced by Nonstationary Disturbances with Zero Mean Values. Journal of Applied and Computational Mechanics 7 (2), 790-797 (2021).
71. Aleksandrov A.Yu., Tikhonov A.A. Stability analysis of mechanical systems with distributed delay via decomposition. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes 17, iss. 1, 13-26 (2021). (In Russian)
72. Aleksandrov A.Yu., Tikhonov A.A. Stability Analysis of Mechanical Systems with Highly Nonlinear Positional Forces under Distributed Delay. Avtomatika i telemekhanika 1, 3-22 (2023). (In Russian) [Eng. transl.: Automation and Remote Control 84 (1), 1-15 (2023)].
73. Aleksandrov A. Yu., Tikhonov A. A. On the attitude stabilization of a rigid body under control with distributed delay. Mechanics Based Design of Structures and Machines 51 (4), 2241-2250 (2023).
74. Aleksandrov A.Yu., Tikhonov A.A. Nonlinear Control with Distributed Delay for the Angular Stabilization of a Rigid Body. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy 9 (67), iss.4, 653-664 (2022). https://doi.org/10.21638/spbu01.2022.408 (In Russian) [Eng. transl.: Vestnik St Petersburg University. Mathematics 55, iss. 4, 426-433 (2022). https://doi.org/10/1134/S1063454122040021].
75. Letov A.M. Flight Dynamics and Control. Moscow, Nauka Publ. (1969). (In Russian)
76. Smirnov E. Ya. Control of rotational motion of a rigid body. Differential Equations 8 (12), 2264-2266 (1972). (In Russian)
77. Meilakhs A. M., Khitrov G. M. On damping the rotation of a solid. Differential Equations 12 (8), 1520-1522 (1976). (In Russian)
78. Meilakhs A. M. On the orientation of an inert solid. Differential Equations 12 (11), 2101 (1976). (In Russian)
79. Lebedev D. V. On the control of the motion of a rigid body with incomplete information about the angular velocity. Automation and Remote Control 12, 5-11 (1978). (In Russian)
80. Ermolin V. S. Stabilization of the rotational motion of a dynamically symmetric rigid body. Control, Reliability and Navigation 3, 132-137 (1976). (In Russian)
81. Meilakhs A. M. On the stabilization of the rotational motion of a solid. Differential Equations 16 (10), 1888-1890 (1980). (In Russian)
82. Lebedev D. V. To control the rotational motion of a solid with incomplete information about the angular velocity vector. Prikladnaia matematika i mekhanika 41 (2), 219-224 (1977). (In Russian)
83. Meilakhs A. M. On the control of the rotational motion of a solid. In: Mathematical Methods for the Study of Controlled Mechanical Systems, 121-126 (1982). Leningrad, Leningrad University Press (In Russian)
84. Dzhumagalieva M. B. On the asymptotic stability of a rigid body rotation motion in the Euler case. In: Problems of Mechanics of Controlled Motion. Nonlinear Dynamic Systems, 66-73 (1984). Perm, Perm State University. (In Russian)
Received: October 2, 2023 Revised: November 4, 2023 Accepted: November 9, 2023
Author's information:
Alexey А. Tikhonov — a.tikhonov@spbu.ru