Стационарное движение твердого тела, несущего силовые гироскопы, и их устойчивость тема автореферата и диссертации по механике, 01.02.01 ВАК РФ

На правах рукописи

Амелькин Николай Иванович

СТАЦИОНАРНЫЕ ДВИЖЕНИЯ ТВЕРДОГО ТЕЛА, НЕСУЩЕГО СИЛОВЫЕ ГИРОСКОПЫ, И ИХ УСТОЙЧИВОСТЬ

Специальность 01.02.01 - теоретическая механика

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математических наук

Москва-2011

1 СЕН 2011

4852588

Работа выполнена на кафедре теоретической механики Государственного образовательного учреждения высшего профессионального образования «Московский физико-технический институт (государственный университет)»

Официальные оппоненты:

академик РАН, доктор физико-математических наук профессор Журавлев Виктор Филиппович

доктор физико-математических наук

профессор Карапетян Александр Владиленович

доктор физико-математических наук

профессор Степанов Сергей Яковлевич

Ведущая организация: Учреждение Российской академии наук Институт Прикладной математики им. М.В. Келдыша РАН

Защита состоится 6 октября 2011 года в 1500 часов на заседании диссертационного совета Д002.240.01 при Институте Проблем Механики им. А.Ю. Ишлинского РАН по адресу: 119526, Москва, проспект Вернадского, д. 101, корп. 1.

С диссертацией можно ознакомиться в библиотеке ИПМех РАН.

Автореферат разослан 11 августа 2011 года.

Ученый секретарь диссертационного совета Д002.240.01 при ИПМех РАН

кандидат физико-математических наук ^ысоева Е.Я.

Общая характеристика работы

Актуальность темы. Системы управления и стабилизации вращательного движения твердого тела, использующие в качестве исполнительных элементов силовые гироскопы или роторы (маховики), разделяются на активные и пассивные. К активным относятся системы, в которых движение гироскопов или маховиков относительно несущего тела является управляемым. В таких системах управление, обеспечивающее требуемое вращательное движение несущего тела, строится на принципах обратной связи и осуществляется с помощью активных моментных устройств, устанавливаемых на осях рамок гироскопов или осях маховиков. В качестве примеров активных систем можно указать системы ориентации космических станций, в которых используются гиродины, и системы ориентации геостационарных спутников, использующие управляемые маховики.

В пассивных гиросиловых системах взаимодействие между несущим телом и гироскопами обеспечивается только за счет реакций связей и устанавливаемых в осях рамок гироскопов пассивных моментных устройств, например, пружин и демпферов. В принципах их работы помимо гироскопических свойств вращающихся тел используются также свойства внешней среды (моменты гравитационных и аэродинамических сил, магнитное поле и др.). Гироскопы (или роторы) применяются в таких системах с целью получить дополнительные восстанавливающие моменты, а также новые стационарные движения, отличные от стационарных движений твердого тела. Кроме того, использование гироскопов с демпфированием в осях рамок обеспечивает во многих случаях асимптотические свойства системы, что особенно важно для практики.

В диссертации рассматриваются задачи пассивной гироскопической и гравитационно-гироскопической стабилизации вращательного движения твердого тела. В связи с этим исследуются свойства стационарных движений твердого тела, несущего силовые гироскопы, в однородном внешнем поле (при отсутствии внешних моментов сил) и в центральном гравитационном поле.

К настоящему времени достаточно подробно изучена задача о стационарных движениях твердого тела, несущего осесимметрич-ные силовые роторы с постоянной скоростью собственного вращения. Для случая однородного внешнего поля полное аналитическое

решение этой задачи для всевозможных вариантов установки оси ротора в несущем теле приведено в книге Й. Виттенбурга «Динамика систем твердых тел».

Большое число работ (В.В. Румянцев, В.Н. Рубановский, Р.В. Лонгман, P.E. Роберсон, В.А. Сарычев, С.Я. Степанов и др.) посвящено исследованию стационарных движений спутника, несущего силовой ротор, в центральном гравитационном поле в рамках ограниченной круговой задачи. Подробное аналитическое решение «прямой» задачи, в которой определяется зависимость стационарных движений от величины кинетического момента ротора и главных центральных моментов инерции спутника, получено для случаев, когда ось ротора параллельна главной оси, либо главной плоскости инерции спутника (В.А. Сарычев, С.А. Мирер, A.A. Дегтярев).

Проведен сравнительный анализ стационарных движений спутника с силовым ротором в ограниченной и неограниченной задаче (В.В. Румянцев, С.Я. Степанов).

Стационарные движения твердого тела, несущего силовые гироскопы, до последнего времени были изучены в меньшей степени. В опубликованных работах на эту тему (В.А. Сарычев, С.А. Мирер, A.B. Исаков, P.C. Суликашвили и др.) в основном исследовались отдельные стационарные движения систем с конкретной схемой установки гироскопов в несущем теле, либо анализ был ограничен узким диапазоном значений параметров системы.

В диссертации основное внимание уделено решениям «прямой» задачи, в которой определяется зависимость стационарных движений от таких параметров, как величина суммарного кинетического момента системы (для случая однородного внешнего поля), значения кинетических моментов роторов гироскопов, углы установки осей рамок гироскопов в несущем теле, значения главных центральных моментов инерции системы.

Целью работы является изучение свойств стационарных движений твердого тела, несущего двухстепенные и трехстепенные силовые гироскопы.

Задачами исследования являются:

- разработка эффективных методов определения стационарных движений и методов анализа их устойчивости.

- поиск решений «прямой» задачи о стационарных движениях твердого тела, несущего силовые гироскопы, в однородном внешнем поле и в центральном гравитационном поле.

- анализ предельных движений системы при наличии диссипации в осях рамок гироскопов.

- анализ стационарных конфигураций систем, состоящих из произвольной совокупности связанных тел, в однородном и центральном гравитационном поле.

Методы исследования. В диссертации применяются методы теории дифференциальных уравнений и теории устойчивости, в том числе первая теорема Рауса и ее модификации, теоремы Ляпунова и Четаева, теоремы Барбашина и Красовского.

Научная новизна. В диссертации получены следующие основные новые результаты, которые выносятся на защиту:

1. Уравнения движения системы «несущее тело - силовые гироскопы» и алгоритм исследования устойчивости стационарных движений в однородном внешнем поле.

2. Теорема о предельных движениях гиростатов в однородном внешнем поле при наличии диссипации в осях рамок гироскопов.

3. Аналитическое решение «прямой» задачи о стационарных движениях твердого тела, несущего один двухстепенной силовой гироскоп, в однородном внешнем поле для различных вариантов установки оси прецессии гироскопа в несущем теле.

4. Аналитическое решение «прямой» задачи о стационарных движениях в однородном внешнем поле твердого тела, несущего систему двухстепенных силовых гироскопов, установленных по коллинеарной схеме.

5. Аналитическое решение «прямой» задачи о стационарных движениях в однородном внешнем поле твердого тела, несущего произвольное число трехстепенных силовых гироскопов в кардано-вых подвесах.

6. Алгоритм анализа устойчивости положений относительного равновесия спутника, несущего двухстепенные силовые гироскопы с диссипацией в осях рамок, в центральном гравитационном поле.

7. Аналитическое решение «прямой» задачи о положениях равновесия на круговой орбите спутника, несущего двухстепенной силовой гироскоп с осью рамки, параллельной одной из главных центральных осей инерции спутника.

8. Аналитическое решение «прямой» задачи о положениях равновесия на круговой орбите динамически симметричного спутника с двухстепенным силовым гироскопом при произвольном расположении оси прецессии гироскопа в корпусе спутника.

9. Алгоритм определения стационарных конфигураций и характера их устойчивости для произвольной системы связанных тел в однородном внешнем поле и в центральном гравитационном поле на круговой орбите.

Достоверность и обоснованность результатов работы обеспечивается корректностью постановок задач, наличием полных и строгих доказательств утверждений. Полученные в работе аналитические выводы подтверждаются результатами численного моделирования.

Практическая ценность работы. Полученные в работе результаты могут быть использованы при проектировании систем гироскопической и гравитационно-гироскопической стабилизации вращательного движения спутников.

Апробация результатов работы. Результаты диссертационной работы докладывались на заседаниях Всероссийского научного семинара имени академика Румянцева В.В. (Москва, МГУ им. Ломоносова, 2008 и 2011 гг.), на заседаниях Всероссийского научного семинара «Механика систем» имени академика А.Ю. Ишлинского (Москва, Институт Проблем механики, 2010 г.), на международном конгрессе «Нелинейный динамический анализ» (Санкт-Петербург 2007), на IX Всероссийском съезде по теоретической и прикладной механике (Нижний Новгород, 2006), на X международном семинаре им. Е.С. Пятницкого «Устойчивость и колебания нелинейных систем управления» (Москва, ИПУ РАН, 2008), на X Международной конференции «Устойчивость, управление и динамика твёрдого тела» (Донецк, 2008), на XI международном семинаре им. Е.С. Пятницкого «Устойчивость и колебания нелинейных систем управления» (Москва, ИПУ РАН, 2010), на X Крымской международной математической школе «Метод функций Ляпунова и его приложения» МРЬ-2010 (Крым, Алушта, 2010), на XXX академических чтениях по космонавтике (Москва, 2006), на ХЫХ научной конференции МФТИ (Долгопрудный, 2006), на заседаниях семинара кафедры теоретической механики МФТИ (Долгопрудный, 2005-2011).

Публикации. Основные результаты диссертации опубликованы в 22 научных работах, в том числе 14 статей [1-14] опубликовано в журналах из Перечня ведущих научных журналов и изданий ВАК и 9 тезисов докладов [15-22].

Личный вклад соискателя. Все исследования, результаты которых изложены в диссертационной работе, проведены лично соискателем в процессе научной деятельности.

Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения и списка литературы. Работа изложена на 249 страницах, содержит 53 иллюстрации. Библиография включает 115 наименований.

Содержание работы

Во введении обосновывается актуальность темы, сформулированы цели и задачи исследования, дается краткий исторический обзор и краткое содержание диссертации.

Глава 1 посвящена методам исследования стационарных движений твердого тела, несущего N двухстепенных силовых гироскопов с постоянными скоростями вращения роторов (рис. 1).

Рассматриваются системы, не обладающие в общем случае свойством гиростата, и используются следующие обозначения:

СЕ ( е, е ) - базис с началом в центре масс системы «несущее тело - гироскопы», параллельный жестко связанному с несущим телом базису, С1 (С1112 ¡3) - базис Кенига, параллельный

инерциальному базису 0\ (0\{\2 ¡,), $ {к=\,Ы) - единичные векторы, определяющие фиксированные направления осей рамок гироскопов в базисе СЕ, хк - углы прецессии гироскопов, х - И-

мерный вектор, составленный из углов хк, Ь к {хк) - единичные

векторы, указывающие текущие направления осей роторов, Н^х^,) = Нк\\ к - кинетические моменты собственного вращения

роторов, Н(х) - вектор суммарного кинетического момента роторов, (*,,.) - тензор инерции к-го гироскопа (ротора вместе с рамкой) в базисе С^Е , получаемом параллельным переносом базиса СЕ в центр масс гироскопа Ск, «Ц*) - тензор инерции системы в базисе СЕ, <й - угловая скорость несущего тела (базиса СЕ) относительно базиса Кенига, заданная проекциями на оси базиса СЕ.

В разделе 1.1 дается постановка задачи и сформулированы основные допущения о геометрии системы и внешнем силовом поле.

В разделах 1.2-1.5 рассматриваются системы со статически уравновешенными гироскопами. В разд. 1.2 определяется зависимость тензоров инерции 3 и J от текущей конфигурации х. Получены формулы для производных от тензоров инерции:

Здесь - кососимметрическая матрица, представляющая собой

матричный оператор векторного умножения: 8^г=84хг. Из (1)

следует, что система «несущее тело - гироскопы» будет гиростатом, т.е. ее тензор инерции J будет неизменным, в том случае, когда помимо статической уравновешенности каждый гироскоп динамически симметричен относительно своей оси рамки . Для

гиростатов неизменны все тензоры 3 к и выполняются равенства

где I - моменты инерции гороскопов относительно осей ьк, г -

произвольный вектор из трехмерного пространства.

Для системы со статически уравновешенными гироскопами кинетический момент К относительно центра масс С имеет вид

дЗ/дхк -дЗк/дхк =ЪкЗк .

(1)

к к'

(2)

K^Jto+XJ^ + H; Н = xk=skxk. (3)

k=i *=i Кинетическая энергия системы относительно базиса Кенига выражается в виде суммы Т = Т2 + ^+Т0 квадратичной, линейной и нулевой формы скоростей ю,х , а обобщенная энергия V = T2~TQ + П с точностью до аддитивной постоянной Т0 (Т0= const вследствие постоянства угловых скоростей фк) определяется формулой

У(ч>,х,х)=±(ягЗ(о+ат^5кхк +jfx Jkxk +77, (4)

Z к=\ ^ Ы

где П - потенциальная энергия всех действующих сил. Функция (4) представляет собой полную энергию движения системы относительно базиса Кенига, вычисленную при «замороженных» роторах гироскопов.

В разделе 1.3 выводятся уравнения движения системы относительно базиса Кенига. Они разделяются на уравнения Пуассона, описывающие движение базиса СЕ относительно базиса Кенига, и уравнения для фазовых переменных ю, х, х . Последние записываются в виде

Jtb+юхJto+X[Jк {хк +сохх^ )+(ш+х4 М Jt\ +Н к )+хк XJt<0]=M. (5)

t=i _

s^JЛ(o>+xt +ft>xxjt)+s[((ox(Jjt(o)+x(t)+H/l.))=Mi; k=\,N. (6)

Здесь M - момент действующих на систему внешних сил относительно ее центра масс, Мк - момент действующих на гироскоп сил

относительно оси sk. Для гиростатов уравнения (5), (6) принимают

следующий упрощенный вид:

Jci»+caxJcB+f;(/txi+fi>x(/tit+HJt)+iixHi)=M. (5*)

*=i

Iktf<b+xk)+*rk{<axHk) = Mk-, k=Uf. (6*)

В разделе 1.4 исследуются движения системы в однородном внешнем поле (при отсутствии внешних моментов сил). В рассматриваемой задаче М = 0, а моменты Мк обусловлены только внут-

ренними потенциальными и диссипативными силами, действующими в осях рамок гироскопов, и определяются выражениями

Мк = М'' + М", где М" = -дП/дхк, М°хк < 0 .

В этом случае уравнения движения (5), (6) образуют замкнутую систему порядка 2/У+З и имеют первый интеграл

~ /V

¿=1

а обобщенная энергия (4) удовлетворяет соотношению

(8)

¿ы

Определение. Стационарными движениями системы в однородном внешнем поле будем называть движения, которые характеризуются постоянными значениями фазовых переменных (о,х.

Уравнения для стационарных движений получаются из уравнений (4), (5) при соэ0,х=0 и имеют вид

сох^о+Н) = 0, 8[(юх(^м+НА))=-аЯ/&4; (9)

Теорема 1.1. Стационарные движения системы соответствуют условно-стационарным точкам функции (4) на многообразии (7) и в пространстве переменных со,х совпадают с условно-стационарными точками «усеченной» функции

Г(ю,х) = К(сй,х,0)=о)^ю/2+Я(х) (Ю)

на «усеченном» многообразии

/(ю,х)=/(ю,х,0)=(-1со+Н)2=^2. (11)

На основании теоремы 1.1 все стационарные движения могут быть найдены как стационарные точки функции Лагранжа с множителем:

1(ш,х) = Г+Я//2=о>^со/2+Я(х)+А(1а)+Н)2/2. (12)

В разделе 1.5 излагается метод исследования устойчивости стационарных движений в однородном внешнем поле.

Теорема 1.2. В окрестности стационарной точки {ю\х\0} функция V (4) на многообразии (7) представляется в виде

у\ = Р+\У{<а,х,К2)\ 1¥(ы,х,К2) = У\

(П)'

(13)

где У7 - строго положительно определенная квадратичная форма скоростей х, а \У{т,х,К2) -усеченная обобщенная энергия (10), вычисленная на усеченном многообразии (11).

Из теорем 1.1 и 1.2 следует, что усеченная обобщенная энергия V (10), вычисляемая на усеченном многообразии (11), служит аналогом эффективной потенциальной энергии. Стационарные точки, в которых функция (10) имеет строгий условный минимум, устойчивы по переменным со , х, х независимо от того, действуют в осях рамок гироскопов диссипативные моменты, или нет. В дальнейшем, следуя Кельвину, все такие точки будем называть устойчивыми в вековом смысле. В свою очередь, стационарные точки, в которых функция V не имеет условного минимума, будем называть неустойчивыми в вековом смысле. Характер вековой устойчивости определяется на основе анализа квадратичной формы с1гЬ на множестве вариаций ¿Лд, с!х, связанных уравнением

Если (¡гЬ является строго положительно определенной на множестве (14), то стационарное движение устойчиво, а если она принимает отрицательные значения на некотором подмножестве множества (14), то движение неустойчиво в вековом смысле.

При наличии диссипации в осях рамок гироскопов возможен более детальный анализ устойчивости стационарных движений с использованием теоремы Барбашина-Красовского:

Если стационарное движение {<а*,х*,0} удовлетворяет условию

1°. В достаточно малой окрестности точки г*={ш*,х*,0} на многообразии (7) нет других предельных движений (целых

траекторий Ф г*, таких что У(г°(!)) =0 ), то в случае вековой устойчивости это движение условно асимптотически устойчиво, а в случае вековой неустойчивости — неустойчиво по Ляпунову.

Здесь под условной асимптотической устойчивостью понимается асимптотическая устойчивость по отношению к возмущениям,

для которых постоянная К2 в интеграле (7) не меняется.

с!/=£/(1со+Н)2 =0.

(14)

Для случая, когда диссипативные моменты действуют в осях рамок всех гироскопов, предельные движения характеризуются постоянными значениями всех углов прецессии гироскопов и описываются уравнениями

Эти уравнения образуют переопределенную систему из N+3 дифференциальных уравнений относительно трех компонент вектора со , а углы хк (/с = 1, /V) выступают в них в роли параметров.

Система (15), (16) имеет тривиальные решения со=ам^, которые соответствует стационарным движениям и описывается уравнениями (9). Для выполнения условия 1° необходимо и достаточно,

чтобы рассматриваемая точка {о/,х*} была изолирована от других стационарных точек на многообразии (11), и чтобы переопределенная система (15), (16) не имела нетривиальных решений в некоторой окрестности точки {со*,х*} на многообразии (11).

Для анализа нетривиальных решений системы (15), (16) используются первые и частные интегралы этой системы.

Уравнение (15) имеет два первых интеграла

Кроме того, для совместности системы (15), (16) уравнение (15) должно иметь N частных интегралов следующего вида:

&ткЗкГ\(ох(3(о+Щ=ьгк(а)><(Зк(о+П1с))+дП/дхк-, к = 1,Ы. (19)

Первые интегралы (17) и (18) задают два эллипсоида в пространстве вектора со, а каждый из частных интегралов (19) задает поверхность второго порядка, которая в частных случаях может вырождаться в пару пересекающихся плоскостей.

Теорема 1.3. Для существования нетривиальных решений системы (15),(16) необходимо и достаточно, чтобы при некоторых

фиксированных значениях углов хк (к — \,Ы) и постоянных К2 и Т*

(15)

(16)

о7\1со/2 = !Г*=соп51, (1Ю + Н)2=К2=СОП81.

(17)

(18)

поверхности (17)-(19) имели в качестве общей линии пересечения «полноценную» кривую (континуум) в пространстве вектора со.

Теорема 1.3 применяется в разделе 2.1 при исследовании нестационарных предельных движений (нетривиальных решений системы (15), (16)) для гиростатов.

В разделе 1.6 дано обобщение результатов, полученных в разделах 1.2 —1.5, на случай статически неуравновешенных гироскопов.

В разделе 1.7 рассматриваются произвольные системы связанных тел (без гироскопов), не обладающие свойством гиростата. Показано, что стационарным движениям в однородном внешнем поле соответствуют конфигурации системы, доставляющие стационарные значения трем функциям

IVк (х) = К2/ик (х) + Щх); к = 1,2,3, (20)

где П(х) - потенциальная энергия внутренних консервативных сил, Лк{х) - главные центральные моменты инерции системы, х - вектор, задающий конфигурацию системы. Достаточным условиям устойчивости удовлетворяют точки строгого локального минимума функции IV3 = К2/Ыъ + П, где - тах.1 к.

В разделе 1.8 излагается метод исследования стационарных движений (относительных равновесий) спутника, несущего двухстепенные силовые гироскопы, в центральном гравитационном поле в рамках ограниченной круговой задачи. Орбитальный базис задается единичными взаимно ортогональными векторами

г=К/Я, п=ш0/(У0, т=пхг, где И - радиус-вектор, соединяющий центр притяжения с центром масс спутника, <х>() - вектор угловой скорости орбитального базиса. Обобщенная энергия спутника определяется выражением

У = Т2+Ж; (Зг^г-пг,1п—2пгН-й\1)/2+77, (21)

где Т2 - строго положительно определенная квадратичная форма скоростей £1,х, представляющая собой кинетическую энергию спутника относительно орбитального базиса, вычисленную при «замороженных» роторах гироскопов, Л - угловая скорость корпуса спутника относительно орбитального базиса, IV - измененная

потенциальная энергия, зависящая только от углов прецессии гироскопов и взаимной ориентации орбитального и связанного с корпусом базиса, Н(х) = Н/си0 - «приведенный» суммарный кинетический момент роторов, П(х) - потенциальная энергия внутренних консервативных сил. Все фигурирующие в формуле (21) векторы и тензор инерции спутника I заданы в базисе СЕ.

Обобщенная энергия спутника удовлетворяет условию (8), а стационарные движения соответствуют стационарным точкам

функции Ш и при учете равенств пгп = 1, ггг=1, п г=0 могут быть найдены как стационарные точки функции Лагранжа

1=}¥1(о1+1{птп12+А1гтг12+Луг (22)

из системы уравнений

дЦдхк =0 Ш\дхк =0; к = ГД

Э£/ап = 0 => -(1п + Й) + А|П + /Цг = 0, (23)

дЬ/дг = 0 =>31г + Д2г + Я3п = 0. Неизвестными в этой системе являются углы прецессии гироскопов хк(к=1,Ы) и три переменные, определяющие ориентацию орбитального базиса относительно корпуса спутника. Для гиростатов первая группа уравнений (23) принимает вид

-пг(84хЙ4)+аЛГ/а^=0; к = П=П/(о20 . (23*)

Анализ вековой устойчивости стационарных движений спутника сводится к исследованию второго дифференциала с12Ь функции (22) на множестве вариаций, связанных уравнениями

пг<Л1=0, гтск=0, птск+гтсЬ=0. (24)

Вариации векторов орбитального базиса выражаются через независимые вариации и,у,и> равенствами

й?п=г/г+ут, ¿/г^-ип+уит, (25)

а множители Лагранжа определяются из уравнений (23) формулами Л1 =пг(Лп + Н), Л2=-Згг1г, Я3=-Зпг1г. (26)

С учетом соотношений (25), (26) квадратичная форма с12Ь на множестве (24) сводится к квадратичной форме от Лг+3 независимых вариаций и,у, и' и гк =с1хк. Для гиростатов эта квадратичная форма

принимает следующий вид:

Ф=С/21|(24)=[4(Пг1П-Гг1Г)+ПгЙ]М2+(П^П-Т^Т+ПгН)У2 + +3(тг1т-гг.1г)и'2-6пг.1тг/и'-6пг.1гуи,+ (27)

к=1 ¡л ) к В разделе 1.9 проводится анализ равновесных конфигураций произвольной системы связанных тел, не обладающей свойством гиростата, в центральном гравитационном поле на круговой орбите. Показано, что равновесным конфигурациям соответствуют стационарные точки следующих шести функций:

W¡jk(x) = (2Jr2Jj-Jk)|2 + П■, ¿Ф]Фк= 1,2,3, (28)

а достаточным условиям устойчивости удовлетворяют точки строгого локального минимума функции №132, где П{х)- потенциальная энергия внутренних консервативных сил, Jl(x)<J2(x)<Jг(x) -главные центральные моменты инерции системы, х - вектор, задающий конфигурацию системы.

В главе 2 исследуются предельные движения системы при наличии диссипации в осях рамок гироскопов.

В разд. 2.1 проводится анализ предельных движений гиростатов в однородном внешнем поле.

Теорема 2.1. 1) Если система не обладает динамической симметрией (АФ В фС), то все множество ее предельных движений исчерпывается стационарными вращениями.

2) Если А = ВфС, то при специальном подборе значений параметров к,Нк{к-\,Щ система допускает нестационарные

предельные движения - регулярные прецессии, изолированные от стационарных движений на многообразии (7).

Следствие из теоремы 2.1. Для гиростатов условие Г теоремы Барбашина-Красовского сводится к изолированности стационар-

ного движения {ю*,х*} от других стационарных движений на многообразии (11) при К2 = [.Тсо* + Н(х*)]2 .

В силу теоремы 2.1 анализ устойчивости по Ляпунову стационарных движений гиростатов в однородном внешнем поле при наличии диссипации в осях рамок гироскопов сводится в основной части к анализу вековой устойчивости, поскольку проверка изолированности стационарных движений на многообразии (11) не представляет принципиальных трудностей.

В разделе 2.2 излагается основанная на результатах работ Пожа-рицкого Г.К. и Лилова Л.К. методика исследования предельных движений системы с частичной диссипацией в малой окрестности положения равновесия с использованием усеченных уравнений линейного приближения.

Рассматривается автономная система, положение которой задается (и+т)-мерным вектором обобщенных координат q , где

чТ=(хТ,уТ), х^(х,,...,х), У'=(У[,...,УП). Предполагается, что известна функция фазовых переменных ^(<7, я), для которой в силу уравнений движения системы выполняется условие У<0. Рассматривается случай, когда диссипация действует только по координатам х{,.. ,,х;, т.е.

У<0 при х^0, Г=0прих = 0. (29)

Запишем уравнения движения системы в виде, разрешенном относительно старших производных:

¿ = /(г,*,х), х = (р(г,х,х). (30)

Здесь гГ =(ут, ут), а / п <р - вектор-функции фазовых переменных размерности 2т и я, соответственно.

Обозначим через с{°(!) целые траектории предельных движений системы, т.е. движений, удовлетворяющих условию У(д°(()) = 0. Пусть система допускает положение равновесия д = 0 (г = 0, х = 0). Поскольку в силу (29) все предельные движения характеризуются тождественным равенством х(/) = 0, для отсутствия в

окрестности точки г = 0, х = 0 целых траекторий д°(() необходи-

мо и достаточно, чтобы в этой окрестности получаемая из (30) при х = 0 система уравнений

¿ = /(2,0,*), 0 = р(2,0,х) (31)

имела только тривиальное решение гзО,

Система (31) представляет собой переопределенную систему дифференциальных уравнений относительно переменных г , а переменные х могут принимать только фиксированные значения, т.е. выступают в роли параметров. Проведя линеаризацию уравнений (31) в окрестности положения равновесия, получим систему

¿ = АГг + Вх, РГ2+Сх = 0 (32)

где матрицы выражаются через производные от вектор-функций / и ср в точке 2 = 0, х = 0 формулами

Аг=а//д2г, Ъ = д//дхт, Р т=д(р/д2т, С = д<р/дхт.

Теорема 2.2 (Пожарицкий Г.К., Лилов Л.К.). Для отсутствия нетривиальных решений системы (31) в малой окрестности положения равновесия г = 0, х = 0 достаточно выполнения условий

АТ В

Рг С

*0, гапк Р АР А2Р ... А2"'-'Р =2т. (33)

Первое из условий (33) обеспечивает изолированность точки # = 0 от других положений равновесия линейной системы (32) и означает, что в рамках нелинейной системы (30) положение равновесия <7 = 0 не является точкой ветвления. Второе из условий (33) означает, что усеченная система

¿ = АГ2, Рг2 = 0, (34)

получаемая из уравнений (32) при х = 0, имеет только тривиальное решение ггО.

Для спутника, несущего силовые гироскопы с диссипацией в осях рамок, функцией У(д, д) служит обобщенная энергия (21), переменными х - углы прецессии гироскопов, а переменными у -углы ориентации спутника относительно орбитального базиса.

В разделе 2.3 выводятся линеаризованные в окрестности положений равновесия уравнения движения спутника на круговой орби-

те и приводится алгоритм анализа устойчивости положений равновесия при наличии диссипации в осях рамок гироскопов.

Нелинейные уравнения движения спутника-гиростата относительно орбитального базиса выводятся из уравнений (5*), (6*) и имеют следующий вид:

n

/ы (-зь)

+(со0 п+£1) X [ 1(ю0 п+П)+Н]=(г X 1г),

/Д8[(П+СУ0пх£2)+Ха.]-(С(;0П+П)7,(5(1.ХН/1.) = М/,; (36)

Здесь <у0 - угловая скорость орбитального базиса, £1 - угловая скорость корпуса спутника относительно орбитального базиса, векторы хк и хк определяются выражениями хк-§кхк, хк = ъкхк, а

Мк - моменты внутренних потенциальных и диссипативных сил в

осях рамок. В записанных уравнениях векторы орбитального базиса задаются в базисе, связанном с корпусом спутника, а их производные по времени выражаются равенствами

п = пх,П, г = гх£2. (37)

В предположении, что моменты Мк обусловлены силами пружин и демпферов, установленных на осях рамок, линеаризованные в окрестности положений равновесия уравнения движения спутника записываются в следующем виде:

•ги'+£^4+ЕЧх(н,-/,п)+пх|;(х,хн,) =г(и)+ад, (38)

4=1 к=1 к=1

/,(8[и'+х^=(54хЙ,)ги-/?А.4-(пгН,+г,)х,; к=Ц/, (39) \|/'=и+\(/хп. (40)

Здесь \|/ - вектор малого поворота корпуса спутника относительно орбитального базиса, хк - отклонения углов прецессии от их значений в положении равновесия, штрихом обозначены производные по безразмерному времени т=соЛ,

н*=н*/®0> Мк=мк/О}0, ск=ск/й)20- к = (41)

- «приведенные» кинетические моменты роторов, коэффициенты вязкого трения и коэффициенты жесткости пружин, соответственно, а векторы Г (у) и Е(11) определяются формулами

Г(Ч/)=3[(гх\|/)хЛг+гх1(гх\|/)], (42)

Р(Ц)=-их(1п+Н)-пхЛ1. (43)

В уравнениях (38)-(40) переменными являются компоненты векторов \(/, и и углы хк, а векторы п, фиксированы и равны их

значениям в положении равновесия.

При х = 0 система (38)-(40) принимает вид

313'и)+1ХчО, (44)

\|/=и+\|/хп, (45)

/,8,ги'=(8,хН,)ги; к=Т^ (46)

и представляет собой переопределенную систему из N+6 скалярных дифференциальных уравнений для шести переменных.

В силу теоремы 2.2 при наличии диссипации в осях рамок гироскопов в малой окрестности положения равновесия спутника заведомо нет целых траекторий других предельных движений, если это положение не является точкой ветвления, и если система (44)-(46) не имеет нетривиальных решений. Отсюда следует, что для анализа устойчивости положений равновесия спутника, несущего двухстепенные силовые гироскопы с диссипацией в осях рамок, применима теорема Барбашина-Красовского в следующей формулировке:

Теорема 2.4. Если для положения равновесия, не являющегося точкой ветвления, усеченная линеаризованная система (44)-(46) имеет только тривиальное решение Ц/ = 0, 11=0, то в случае вековой устойчивости это положение асимптотически устойчиво, а в случае вековой неустойчивости — неустойчиво по Ляпунову.

Систему (44)-(46) можно привести к виду (34), а отсутствие у нее нетривиальных решений описать вторым из условий (33), если разрешить уравнение (44) относительно переменной и' и подставить полученное выражение в уравнения (46).

В главе 3 приводятся подробные решения «прямой» задачи о стационарных движениях твердого тела с двухстепенными силовы-

ми гироскопами в однородном внешнем поле. Здесь рассматриваются системы, обладающие свойством гиростата, т.е. гироскопы статически уравновешены и обладают динамической симметрий относительно осей рамок. Предполагается также, что оси роторов ортогональны осям рамок, а в осях рамок могут действовать только диссипативные силы. Через с|,е2,е3 обозначается базис главных

центральных осей инерции, а через А,В,С - соответствующие главные моменты инерции системы. Уравнения для стационарных движений таких систем выводятся из условий стационарности функции Лагранжа (12) и имеют следующий вид:

51/5(0=0 => Л(1<»+Н) = 0=>с)+/1(.1га+Н)=0, (47)

дЬ/дхк= 0 => Л(1<о+Н)г(8,хНА.)=0; (48)

Для рассматриваемой системы величина суммарного кинетического момента ^=|К| является первым интегралом, а «прямая» задача состоит в том, чтобы из уравнений (47), (48) и (11) найти значения вектора со и углов хк (векторов И ) в зависимости от К .

Для стационарных движений векторы ш и К связаны формулой К=+Н, вследствие которой вектор си однозначно определяется по значениям векторов К и х. Поэтому в большинстве случаев определяемые ниже стационарные движения описываются вектором {К,х\, где вектор К задается своими компонентами

К1,К2,К3 в базисе в! е2е3. Такое описание удобно для нахождения

числа стационарных движений в зависимости от величины К , а также точек ветвления, в которых меняется характер устойчивости.

Заключение о характере вековой устойчивости стационарных движений делается на основе анализа квадратичной формы

с]2Ь = \т]у+Л\т52\+2Л\т^(8кхЯк)гк +

(49)

/,/Ы 4=1

на множестве вариаций V=<Ло, гк =скк, связанных уравнением

с// = 0 => (50)

к=\

Для стационарных движений с ненулевой угловой скоростью (при Лф0) уравнение (50) принимает вид угЛС=0, т.е. вариации гк в

квадратичной форме (49) независимы.

Если квадратичная форма (49) на множестве (50) является строго положительно определенной, то стационарное движение удовлетворяет достаточным условиям устойчивости (устойчиво в вековом смысле). Если же квадратичная форма (49) принимает отрицательные значения на некотором подмножестве множества (50), то стационарное движение неустойчиво в вековом смысле. При наличии диссипации в осях рамок гироскопов в силу теоремы 2.1 для изолированных на многообразии (11) стационарных движений из вековой устойчивости следует условная асимптотическая устойчивость, а из вековой неустойчивости - неустойчивость по Ляпунову.

В разд. 3.1-3.5 исследуются стационарные движения системы с одним гироскопом. Система допускает стационарные движения с нулевой угловой скоростью несущего тела (неподвижные состояния). Им соответствуют решения системы (47),(48) при Л-0, а угол х может принимать любое значение от 0 до 2л. Неподвижные состояния существуют при К~Н, а в пространстве вектора К они образуют окружность радиуса Я , ортогональную оси прецессии «. Как показал анализ, каждое неподвижное состояние устойчиво только по части переменных, а именно, по переменным ш и х . При наличии диссипации в оси рамки гироскопа неподвижные состояния обладают условной асимптотической устойчивостью по указанным переменным.

Стационарные движения с ненулевой угловой скоростью несущего тела определяются решениями системы (47),(48) при Лф0 и зависят от расположения оси рамки 5 в базисе е ,е2,е3.

В разд. 3.1 на основании уравнений (47),(48) получены формулы

о) = Я—^ (пх.Гп); п(х) = <ЭЬ/бх = 8х11(л:) , (51) (пх

К=^ш + Н = Я <-пх3пУ** (Пх1п), (52)

(п х «Г п) «Г п

дающие параметрическое представление стационарных движений системы через угол х поворота рамки гироскопа. В случае, когда

ось рамки гироскопа s не параллельна ни одной из главных плоскостей инерции системы, формулы (51) и (52) определяют все множество стационарных движений системы с ненулевой угловой скоростью. Если же ось s параллельна главной плоскости инерции, то в точках х*, для которых вектор п(х*) параллелен главной плоскости или главной оси инерции, в правых частях этих формул возникает неопределенность типа 0/0. Все соответствующие указанным точкам х* стационарные движения определяются непосредственно из системы (47)-(48).

В разд. 3.2 исследуется система, в которой ось рамки гироскопа параллельна одной из главных центральных осей инерции: s=e3,

т.е. вектор h определяется выражением h = et cos х + е2 sinx. Полученные для этой системы стационарные движения с ненулевой угловой скоростью описывается следующими решениями:

Ь=5,е1(*=0,я-), К2= 0, К3= 0, K{=S2K. (53)

Ь=3^2(Х=тг/2-Л/2)„ К}= 0, К3-0, K2=S2K. (54)

h = ^ep К2= 0, К{ =5Х СН/(С—А), K3=sjK2-K2 . (55)

hК= 0, К 2 -5{ СИ/(С-В), K^djK2-K; . (56)

Здесь и далее для записи многозначных решений используются символы <5j, i5"2,..., каждый из которых может принимать два значения +1 и -1. Все эти решения образуют прямые линии в пространстве вектора К, а число решений {К,х} в зависимости от К составляет либо 8, либо 12, либо 16.

В рассматриваемом случае без ограничения общности можно положить А>В. Тогда полученные на основе анализа квадратичной формы (49) на множестве (50) результаты исследования вековой устойчивости формулируются следующим образом (здесь и далее перечисляются только стационарные движения, удовлетворяющие достаточным условиям устойчивости, а остальные движения по умолчанию неустойчивы в вековом смысле):

1). Если ось рамки гироскопа установлена параллельно оси наименьшего или среднего момента инерции ( А > С ), то для каждого значения К> 0, кроме точки К=Н, вековой устойчивостью

обладают два движения, определяемые решениями (53) при

2). Если ось рамки установлена параллельно оси наибольшего момента инерции (А<С), то в диапазоне 0<K<K = CHj(C-A) система имеет два устойчивых движения (53) при S =S =±\, а в

диапазоне К>К - четыре устойчивых движения (55).

В разд. 3.3 исследуется система, в которой ось рамки гироскопа параллельна одной из главных центральных плоскостей инерции:

s=e3 cos/?-e2sin/?, т.е. h = cosx + (e2 cos /? + е3 sin/?)sin;t.

В этом случае множество стационарных движений с ненулевой угловой скоростью описывается следующими решениями:

h=Slel(x=0,ír), К2= О, К3= 0, Kl=S2K. (57)

h = Sl (е cos р+е sinр) (х=я/2 ,-л"/2), Кх = О, К3 (CHcos Р + 5у (В-С)Кг) = К2 ВН sin fi.

(58)

V, V ИР ( А-С , В~АХ\

К{Х)=(А-В)(А-С){С^Х+В^[С^ + С^)} (59)

Р=С(В-А)со$,2р + В(С-А)?,т2р.

Решения (57) образуют прямые линии, проходящие через начало координат, решения (58) - две гиперболы, расположенные в плоскости К2,Кг, а решения (59) - эллипс с полуосями

а2 -К2(0)=Н2Р2/(А~В)2(А-С)2,

Н2Р2

b =К (я/2)=-

Г(С-А)2 {А-В)2Л

(.А-ВУ(А-СУ{С-ВУ\ cosL р sin1 р плоскость которого проходит через начало координат и ортогональна вектору

^=e3(^-C)sin/?+e2(,4-S)cos/?.

В зависимости от величины К система имеет либо 8, либо 12, либо 16 стационарных движений. По результатам анализа квадратичной

формы (49) на множестве (50) установлено, что характер вековой устойчивости зависит от параметров системы следующим образом:

1) А>В>С. Для любого К>0, кроме точки К-Н, вековой

устойчивостью обладают только два решения (57) при 8 = 82 =±1.

2) В>А>С, Р<0. В диапазоне 0<К<а достаточным условиям устойчивости удовлетворяют два решения (57) при 8 = 81 =±1, в

диапазоне а<К<Ь - четыре решения (59), в диапазоне К>Ь - два решения (58) при 31К2 >-НР/(В-С)(В-А)оо&р.

3) В>А>С, Р>0. Для любого К> 0 вековой устойчивостью обладают только два решения (58) при 8^К2 > 0.

4) В>С>А, 1В(С-А)г. Для этой области параметров достаточным условиям устойчивости удовлетворяют те же движения, что и для области 3).

5) В>С>А, tg4/3>C(B—A)3/B(C-A)3 . Во всем диапазоне 0<К<со достаточным условиям устойчивости удовлетворяют два решения (58) при 8{К2 >0. Кроме того, в диапазоне

Я2((5мп^)2/3 +{Ссо$/3)т)ъ!{В-С)2 <К2 <Ьг вековой устойчивостью обладают два решения (58) из интервалов

В разд. 3.4 исследуется динамически симметричная система (А=В^С), в которой ось рамки гироскопа расположена под углом Р к оси симметрии е3. В рассматриваемом случае главные оси

е,,е2 можно выбрать так, чтобы вектор 5 был параллелен плоскости е2,е3,т.е. ъ^е^оър-е^тр, а стационарные движения определить на основе решений, полученных в разд. 3.3. Установлено, что для рассматриваемой системы стационарные движения исчерпываются решениями (57) и (58), где А = В . При этом, если А>С, то во всем диапазоне 0<К<оо достаточным условиям устойчивости удовлетворяют только два решения (57) при 8 =82=±\, а в случае А<С - только два решения (58) при 0<К28х <СНойърЦС-А).

В разд. 3.5 исследуется система, в которой все три проекции оси рамки гироскопа на главные центральные оси инерции отличны от нуля. В этом случае все стационарные движения с ненулевой угловой скоростью описываются формулой (52). С использованием этой формулы «прямая» задача определения зависимости стационарных движений {К,л:} от величины К сводится к нахождению корней многочлена восьмой степени

(и -6, )2 (и -62 )2 Р4 - (и - л, )2 (и-а2)2(и2 +\)К2 j Н2 =0, (60) в котором коэффициенты а{,а ,bvb2 зависят от главных центральных моментов инерции системы А,В,С и проекций s оси s на главные оси е ,е ,е3, а Р4 - многочлен четвертой степени по переменной и = tgx, не имеющий вещественных корней. При этом каждому вещественному корню многочлена (60) соответствует два стационарных движения, вследствие чего для каждого значения К система имеет не более 16 стационарных движений.

Подробный анализ показал, что при любом значении К> 0 уравнение (60) имеет, по крайней мере, 4 вещественных корня и, следовательно, система имеет, по крайней мере, 8 стационарных движений. При этом для сравнительно малых значений величины К/Н число вещественных корней равно 4, т.е. система имеет 8 стационарных движений, а при достаточно больших значениях К/Н число вещественных корней равно 6, т.е. число стационарных движений равно 12.

Более детальное исследование зависимости стационарных движений от величины К проводилось с использованием численного анализа. Исследовались корни х(К) полученного из (52) уравнения

Ply п\ п\ п\ 1 1

Здесь n(jc)=sxh(jt), nk =п-е^ k = 1,2,3, D=(A-B)(B-C){C-A). Функция F{х) (61) является периодической с периодом п, а число стационарных движений системы равно удвоенному числу решений х(К) уравнения (61) на полуинтервале [0,тг).

По результатам численного анализа установлено, что поведение функции (61) зависит от положения оси рамки б следующим образом. Октант, образованный осями е1,е2,е3, делится на четыре области 5 (рис. 2). Если вектор я лежит в центральной области 50, то функция (61) имеет на полуинтервале [0,/г) единственное ненулевое стационарное значение К1=К(х[), вследствие чего в диапазоне О <К<К{ система имеет 8 стационарных движений, а в диапазоне К>К{ - 12 стационарных движений.

Рис. 2

Если вектор в лежит в одной из областей 5'1,5'2,5'3, примыкающих к главным осям инерции, то функция (61) имеет на полуинтервале [0, л) три ненулевых стационарных значения К] = К(х{),

К2=К(х2), К3 = К(х}). При этом в диапазоне О<К<К{ система имеет 8 стационарных движений, в диапазоне К]<К<К2 - 12 стационарных движений, в диапазоне К2<К<К3 - 16 стационарных движений, в диапазоне К>КЪ - 12 стационарных движений. Установлено, что при А>В>С и для любого значения

А'>0 система имеет только два симметричных друг другу устойчивых в вековом смысле стационарных движения

{х*(Д0, К*(К)}, {х*(К)+л, -К*(К)}. (62)

Если же 8е53, то в диапазоне К2<К<К3 помимо (62) система

имеет еще два симметричных друг другу устойчивых в вековом смысле стационарных движения.

Из результатов проведенного в разд. 3.2-3.5 анализа следует, что независимо от расположения оси рамки гироскопа все стационарные движения с ненулевой угловой скоростью представляют собой изолированные точки на многообразии (11). Поэтому в силу теоремы 2.1 при наличии диссипации в оси рамки стационарные движения, обладающие вековой устойчивостью, будут условно асимптотически устойчивыми, а все неустойчивые в вековом смысле стационарные движения будут неустойчивы по Ляпунову.

В разд. 3.6 исследуется стационарные движения твердого тела, несущего систему из N>2 двухстепенных силовых гироскопов, в однородном внешнем поле. Здесь подробно рассматривается случай, когда гироскопы установлены в несущем теле по коллинеарной

схеме, т.е. оси рамок параллельны друг другу: я к = 1,]У.

В указанном случае стационарным движениям с нулевой угловой скоростью ю (неподвижным состояниям), соответствуют решения системы (47),(48) при Л=О, для которых углы прецессии гироскопов хк могут принимать любые значения. Полагая для

определенности //, >Н2>..>Н^ , получим, что неподвижные состояния определены в диапазоне

/?<|К|<7? , где Д = Р = Н-^Нк. (63)

Ш к=2

При р<О область определения неподвижных состояний описывается неравенством |К|<Л. В пространстве кинетического момента К множеству неподвижных состояний системы соответствует расположенное ортогонально к оси в кольцо с внешним радиусом Л и внутренним радиусом р .

Стационарные движения с ненулевой угловой скоростью описываются двумя группами решений. Первая группа имеет вид

N

К=±К$, н=]Гн,=-тг-к. (64)

в Л в

Эти решения представляют собой вращения вокруг оси в и определены в диапазоне

При этом для внутренних точек указанного диапазона каждому значению вектора Н, определяемого соотношением (64), соответствует для случая N=2 два значения вектора х (две комбинации векторов Н1;Н2), а в случае N>3 - (N-2) -мерное множество значений вектора х.

Для второй группы решений оси роторов параллельны, т.е.

Н,=Я,Ь, = ^=±1, (65)

а вектор <э и угол х (вектор И) определяются из системы

«+Я(1и+Н)=0, Я^ю+Н)г(8х11)=0, (66)

которая описывает исследованные в разделах 3.2-3.5 стационарные движения твердого тела, несущего один двухстепенной гироскоп с

N

кинетическим моментом ротора Н=//И , где Н-Н^+^8кНк .

к=2

Установлено, что в диапазоне (63) вековой устойчивостью (по переменным ш и х) обладают только неподвижные состояния, образующие в совокупности устойчивое инвариантное множество. При наличии диссипации в осях рамок всех гироскопов это множество асимптотически устойчиво, а система будет приходить в состояние покоя из любого начального состояния, для которого величина кинетического момента принадлежит отрезу /?<|К|<Й .

При |К|>7? вековой устойчивостью обладают только те из решений системы (65),(66), для которых кинетические моменты роторов направлены в одну сторону НА = Н Ь; к=\,Ы, и которые соответствуют устойчивым стационарным движениям системы с одним гироскопом, имеющим кинетический момент ротора .

При |К|</9 вековой устойчивостью обладают только решения,

для которых Н =//И, Н^ =-Нк И; к=2,Ы, соответствующие

устойчивым стационарным движениям системы с одним гироскопом, имеющим кинетический момент ротора Н=р\\.

Из полученных ранее результатов (разд. 3.2-3.5) следует, что в диапазонах |К|>Л и |К|</? система может иметь либо два, либо четыре устойчивых в вековом смысле стационарных движения.

Аналитическое исследование стационарных движений системы при произвольной неколлинеарной схеме установки гироскопов в несущем теле оказывается намного сложнее, чем в случае коллине-арной схемы. В разд. 3.6 рассмотрен пример системы с двумя гироскопами, оси рамок которых взаимно ортогональны и параллельны главным осям инерции системы. В зависимости от величины |К| такая система имеет от 16 до 48 стационарных движений с ненулевой угловой скоростью, а число устойчивых в вековом смысле движений равно либо два, либо четыре, либо восемь. Установлено, что система обладает тем же свойством, что и исследованные выше системы с коллинеарной схемой установки гироскопов, т.е. при наличии диссипации в осях рамок гироскопов система асимптотически стремится к состоянию покоя, если значение модуля суммарного кинетического момента системы принадлежит отрезку возможных значений модуля суммарного кинетического момента роторов.

В главе 4 исследуются стационарные движения (относительные равновесия) спутника с одним двухстепенным силовым гироскопом в центральном гравитационном поле в ограниченной круговой задаче. Предполагается, что спутник обладает свойством гиростата, а ось рамки ортогональна оси ротора. В случае, когда потенциальные силы в оси рамки отсутствуют, система (23), описывающая стационарные движения, принимает вид

Здесь через Н = //Ь обозначен «приведенный» кинетический момент ротора, имеющий размерность момента инерции.

«Прямая» задача состоит в том, чтобы при заданном положении оси рамки в определить из уравнений (67)-(70) положение оси ротора Ь (угол х) и взаимную ориентацию орбитального и связанного с корпусом спутника базиса в зависимости от величины Н .

-1п-Н + Я,п + Я3г = 0, 3.1Г + /12Г-1-/!.3П = 0 ,

8г(пхН) = 0, 8ГЬ = 0.

(67)

(68)

(69)

(70)

Для определения решений системы (67)-(70) использовался следующий алгоритм. Были введены скалярные функции вектора г

/ = гт5 г, ё = 4гт32г-(гтЛг)2 . (71)

Из уравнений (68) получены выражения

Л2=-3/, ^=аЪё, п = о-(/г-аг)/г; о- = ±1, (72)

с учетом которых вектор Н определяется из (67) формулой

П = а{{32г-/Зг)1ё + Ъёг-Ах{Зг-/г)1ё). (73)

После подстановки соотношений (72), (73) в уравнения (69), (70) получены следующие уравнения:

*т32 г- (/ + Я,) 3 г + (3 + А, />гг = 0. (74)

5г(ЛгхЛ2г)-/87'(гхЛ) + (/2-Зя2)8г(гхДг) = 0. (75)

Уравнение (75) содержит только переменную г. Поэтому задача определения множества всех стационарных движений системы сводится к нахождению множества всех решений г уравнения (75), принадлежащих единичной сфере гтг = 1. Соответствующие значения векторов п и Н можно определить как двухзначные функции вектора г с помощью формул (72), (73), выразив множитель Я, из

уравнения (74). При этом решение «прямой» задачи может быть найдено на основе анализа поведения величины Н на решениях уравнения (75).

В разд. 4.1 исследованы положения равновесия спутника в случае, когда потенциальные силы в оси рамки гироскопа отсутствуют, а ось рамки параллельна одной из главных центральных осей инерции: я = е3, т.е. кинетический момент ротора выражается формулой

Н=ЯЬ=Я(е1со8х+е281пх). (76)

Положения равновесия в этом случае определяются следующими группами решений:

г=^е3, п=32ек, Ь=<У3е,; ¿=1,2, Яе(0,+оо). (77)

г=б1е3_к, п=32ек, Ь=33ек; к=1,2, Яе(0,+оо). (78)

г = 5. е, ,, h = <£,е . , n -S-. е, cosy + S^ eJsin у|,

, i, (79)

cosy = -Hjak, Яб(0,|а{|]; i = l,2.

h = S3 ck, r=8l (S} e3 cos у - д2 ck |sin.y|), n=дъ ek cos у+ёгс3 |sinj>|,

со&у=-НЦАак), He{0,^ak\]; к=1,2. ^

Здесь ax-a- A-C, a2 = b = B-C.

Решения (77),(78) определяют в совокупности 32 «прямых» положения равновесия (главные оси инерции спутника параллельны осям орбитального базиса) для каждого значения Н> 0. В свою очередь, каждое из решений (79),(80) задает 16 «косых» положений равновесия, области существования которых ограничены полуинтервалами значений Н. В правых граничных точках указанных полуинтервалов решения (79),(80) ответвляются от соответствующих решений (77),(78).

При близких к нулю значениях Н число положений равновесия, описываемых формулами (77)-(80), равно 64, а по мере увеличения параметра Н при переходе через каждую точку ветвления уменьшается на 8 единиц, достигая минимального числа 32.

Если выполняется неравенство s=a/b>0, т.е. когда ось рамки

гироскопа параллельна оси наибольшего или наименьшего моментов инерции спутника, помимо (77)-(80) спутник имеет следующие положения равновесия:

cos2х = [£Я2 -ь4ab(£ +1)- 5\a\^H2+4ab]/(£~l)H2

■ n = h cos у + е3 sin у; cosy = (.H/^H2 + 4ab)sign(C-A) (81)

r = ±(Jn + H- n(nrJn + nrH))/|Jn + H- n(nrJn + nrH)|

Диапазон значений кинетического момента Н, для которых определены решения (81), задается неравенствами

Vi6b2 -4ab <H<Á6a2 -4ab при a/6<4,

л/a2 —4ab <H<^ll6a2 —4ab при a/A>4.

Каждой точке H из диапазона существования рассматриваемых решений соответствует 16 различных положений равновесия.

В рассматриваемой задаче без ограничения общности можно положить А>В. Тогда полученные на основе анализа квадратичной

формы (27) результаты исследования вековой устойчивости формулируются следующим образом:

1). Если ось рамки гироскопа установлена параллельно оси наименьшего или среднего момента инерции (А>В>С или А>С>В), то вековой устойчивостью обладают только четыре решения (77) при к = 1, 82=д3 для всех значений Н> 0.

2). Если ось рамки параллельна оси наибольшего момента инерции (С>А>В), то в диапазоне Н>С-А вековой устойчивостью обладают только четыре решения (77) при к=1, S2 , а в диапазоне 0<Н<С-А - восемь решений (79) при & = 1.

Исследован характер устойчивости по Ляпунову найденных положений равновесия спутника при наличии диссипации в оси рамки гироскопа. Анализ проводился с использованием теоремы 2.4 и других известных теорем о неустойчивости.

Установлено, что «прямые» положения равновесия (77),(78) не удовлетворяют теореме 2.4, а диссипация в оси рамки не приводит к асимптотической устойчивости тех из решений (77), которые устойчивы в вековом смысле. В то же время показано, что все неустойчивые в вековом смысле положения равновесия из множества (77) и (78), за исключением отдельных точек ветвления, неустойчивы по Ляпунову при наличии диссипации.

«Косые» положения равновесия (79)—(81) за исключением отдельных точек Н* удовлетворяют теореме 2.4. Следовательно, при наличии диссипации устойчивые в вековом смысле решения из множества (79) асимптотически устойчивы, а все остальные решения (79)—(81) для значений НфН" неустойчивы по Ляпунову.

Из результатов проведенного анализа следует, что для рассматриваемых вариантов установки гироскопа в корпусе спутника (ось рамки параллельна одной из главных осей инерции спутника) асимптотически устойчивые положения равновесия реализуются только в случае, когда ось рамки гироскопа параллельна оси наибольшего момента инерции спутника.

В разд. 4.2 исследуются положения относительного равновесия динамически симметричного спутника, несущего двухстепенной силовой гироскоп с осью рамки, расположенной под углом /? ^ 0 к оси динамической симметрии е3. В этом случае кинетический момент ротора выражается через угол прецессии х формулой

Н = Ни = Н[е1 созх + (е2со8/?+е3 эт /?)51п х],

а положения равновесия при отсутствии пружины в оси рамки определяются следующими тремя группами решений:

11=£,е,, п=£2е,, г=63ек, Л:=2,3. (82)

Р=(С-А)/Н,

(83)

п=е.,со5зн-е351пу, г=<? е .

(84)

Ь = ^(е2соз/?+е35ту9), /{у)=Ад^Р, п=е2со57+е35ту, г=£ (е3со8>>-е вт^)-Здесь через /{у) обозначена функция

7=(81п2/3/?+со52/3/?)3/2 (85)

- ее стационарные значения. Уравнения /{у) = 8^Р и /{у) = 45хР,

из которых определяются значения переменной в (83) и (84), сводятся к алгебраическим уравнениям четвертой степени.

Решения (82) задают 16 «прямых» положений равновесия, которые определены для всех значений Н>0. Решения (83) задают 8

положений равновесия в диапазоне |Р|</ и 16 положений равновесия в диапазоне |Р| > / , а решения (84) - 8 положений равновесия в диапазоне 4[Р|</ и 16 положений равновесия в диапазоне 4|Р| > / . Значения НХ-\С-А[// и #2=4|С-Л|// являются точками ветвления положений равновесия. При близких к нулю значениях Н число положений равновесия равно 48, а по мере увеличения параметра Н при переходе через каждую точку ветвления уменьшается на 8 единиц, достигая минимального числа 32.

Для симметричного спутника при близких к нулю значениях Н число положений равновесия меньше, чем для несимметричного, потому, что область определения части решений вырождается в точку Н = 0 .

Установлено, что при А>С, т.е. для «вытянутого» вдоль оси симметрии спутника, вековой устойчивостью обладают только четыре решения (82) при к=3, 3 =32, определенные для любого

значения Н>0. При С>А, т.е. для «сплюснутого» вдоль оси симметрии спутника, вековой устойчивостью обладают те из решений (83), которые удовлетворяют условиям

Неравенства (86) выделяют из решений (83) четыре устойчивых положения равновесия, определенных для всех значений Н >0.

Исследован характер устойчивости по Ляпунову найденных положений равновесия симметричного спутника при наличии диссипации в оси рамки гироскопа. Установлено, что для устойчивых в вековом смысле «прямых» положений равновесия из множества (82) диссипация в оси рамки не приводит к асимптотической устойчивости, а все неустойчивые в вековом смысле решения из множества (82), кроме, быть может, отдельных точек Н*, неустойчивы по Ляпунову при наличии диссипации.

Решения (83),(84) за исключением отдельных точек Н* удовлетворяют теореме 2.4. Отсюда следует, что при наличии диссипации устойчивые в вековом смысле решения из множества (83) обладают асимптотической устойчивостью, а все остальные решения (83),(84) для значений Н^Н* неустойчивы по Ляпунову.

Для асимптотически устойчивых положений равновесия из множества (83) определены оптимальные по быстродействию значения параметров спутника. Для этого исследовались корни характеристического уравнения линеаризованной системы (38)-(40) и определялась степень устойчивости взятое с обратным знаком значение вещественной части самого правого корня. Численным способом с применением метода случайного поиска найдено максимальное значение £тах=0.42, достигаемое при следующих значениях параметров:

Здесь через Н и ¡и обозначены «приведенный» кинетический момент ротора и «приведенный» коэффициент вязкого трения.

Исследованы положения равновесия динамически симметричного спутника с пружиной в оси рамки гироскопа, удовлетворяющей условиям

3=\, р<у<?г/2 и 5=-\, х+Р<у<Ъя12.

(86)

>0=42.7°, т4/С=0.58, #/С=1.53, /¿/С=1.2.

(87)

Здесь с - «приведенный» коэффициент жесткости пружины, а х°=7г/2 - угол прецессии, соответствующий ненапряженному состоянию пружины. Для этого случая положения равновесия описываются следующими решениями:

h =(е2 cos ft+е3 sin Д), n =е2 cos j>+е3 siny,

r=Siei,f(y)=P.

h =(e2 cos/?-f e3 sin/?), n =e2 cos^+e3 sin^,

r=¿3(e3cos.}'-e2sin>0, f(y) = 4P.

В совокупности решения (89), (90) определяют 8 положений равновесия в диапазоне 4|Р|</ , 12 положений равновесия в диапазоне f ¡4<\P\<f и 16 положений равновесия в диапазоне \Р\> f .

При С>А (для «сплюснутого» спутника) вековой устойчивостью обладают только решения (89), принадлежащие интервалам

P<y<n¡2, Зф<у<у2 = 2/Т-circtg((tg/?)1,3).

Эти решения определяют четыре устойчивых положения равновесия в диапазоне О <H<H={C-A)U и два устойчивых положения равновесия в диапазоне Н>Н{.

При А > С (для «вытянутого» спутника) вековой устойчивостью обладают только решения (90) из интервала 0 <у</3, которые определяют два устойчивых положения равновесия для каждого значения параметра #>0.

Установлено, что все решения (89),(90), за исключением отдельных точек Н*, удовлетворяют теореме 2.4. Поэтому, если II^ Н*, то при наличии диссипации в оси рамки из вековой устойчивости положения равновесия следует его асимптотическая устойчивость, а из вековой неустойчивости неустойчивость по Ляпунову.

Из результатов проведенного анализа следует, что при отсутствии пружины на оси рамки гироскопа реализация асимптотически устойчивых положений равновесия возможна только для «сплюснутого» вдоль оси симметрии спутника. Установка пружины на ось рамки дает возможность получить асимптотически устойчивые положения равновесия и для «вытянутого» спутника. Но максимально

35

(89)

(90)

достижимое значение степени устойчивости в этом случае почти на порядок меньше, чем для «сплюснутого» спутника.

В главе 5 исследуются свойства стационарных движений твердого тела, несущего N трехстепенных силовых гироскопов в кар-дановых подвесах, в однородном внешнем поле.

Обозначим через i , s и единичные векторы, указывающие направления осей внешних рамок, внутренних рамок и роторов, соответственно. Оси внешних рамок фиксированы в несущем

теле, положения осей внутренних рамок sk определяются углами прецессии у/к, а положения осей роторов h^. - углами прецессии у/к и нутации вк. Предполагается, что для каждого гироскопа ось внутренней рамки ортогональна к оси внешней рамки и к оси ротора, т.е. ik -sk =0, sk -hj, =0, а кинетические моменты вращения роторов Нд. = Нк Ьк постоянны по величине, т.е. Нк - |Hfc| = const > 0.

Предполагается также, что система «несущее тело - гироскопы» обладает свойством гиростата. Указанное свойство можно обеспечить соответствующим подбором моментов инерции роторов и рамок и имеет место при выполнении следующих условий:

1) ротор каждого гироскопа статически уравновешен и динамически симметричен относительно оси ht.

2) Система «внутренняя рамка - ротор» статически уравновешена и динамически симметрична относительно оси s^. Момент

инерции этой системы относительно оси %к обозначим через Вк.

3) весь гироскоп статически уравновешен и динамически симметричен относительно оси внешней рамки ¡А . Момент инерции

гироскопа относительно оси обозначим через Ак.

В предположении, что в осях рамок могут действовать только моменты внутренних диссипативных сил, движение системы в фазовых переменных о), i//(_, фк, 6к, 6к описывается уравнениями

+п 'a- Чл) * н *)=

к=1

/lt(i>+^)+^(itxsi)7'(Ht-fit«))+i^(<axHi) = M1/t; k = \N. (92) ЯД+s>+^(i,xsi)7'со]-(8лхН,)г(со+^) = М2/;; (93)

где Мц и М2к~ моменты дисеипативных сил в осях прецессии и нутации, соответственно. Эти уравнения имеют первый интеграл

7=[Jco+H+X(^ ук\к +Вк eksk)f =К2 = const, (94)

k=1

а обобщенная энергия V = Г? -TQ, которая, как и в задаче с двухстепенными гироскопами, с точностью до аддитивной постоянной Г0

совпадает с кинетической энергией движения относительно базиса Кенига, вычисленной при «замороженных» роторах гироскопов, удовлетворяет условию

F((o,V,0,4/,8) = f (М1кфк+М2квк)<0 . (95)

k=I

Здесь V)/ и 0 — iV-мерные векторы, составленные из углов прецессии и нутации, соответственно. Уравнения для стационарных движений получаются из уравнений (91)—(93) приравниванием к нулю

производных cb,\j/,O,V|),0 и имеют следующий вид:

cox(Jto+H)=0, (96)

¡Аг(<ахН,)=0;£=У^, (97)

S£(coxHa.)=0 -k=hN. (98)

Доказанные в главе 1 теоремы справедливы и для рассматриваемой системы с трехстепенными гироскопами. По теореме 1.1 стационарные движения соответствуют условно-стационарным точкам

«усеченной» кинетической энергии V-Gt

на многообразии /(w,i|/,G) = (Jtfl+H)2 =К2 =const, (99)

где К^со+Н - «усеченный» кинетический момент системы. В свою очередь, по теореме 1.2 точки строгого условного минимума функции V удовлетворяют достаточным условиям устойчивости, а точки, в которых функция V не имеет условного минимума, неустойчивы в вековом смысле.

Условно-стационарные точки функции V на многообразии (99) могут быть найдены как стационарные точки функции Лагранжа

L = V + À//2 = corJra/2+A(Jco + H)2/2 (100) из системы уравнений

dL/e<o = Jw+;iJ(Jco + H) = 0 => ю = (101)

дЬ/ду/к = ЛKr(ifcxHA) = 0; k=l,...,N, (102)

dL/eek = ÀKT(skxnk) = 0-, k=\,...,N, (103)

которая эквивалентна системе (96)-(98).

Решениям системы (101)—(103) при Л=0 соответствуют стационарные движения с нулевой угловой скоростью несущего тела (неподвижные состояния). Для них углы прецессии и нутации гироскопов могут принимать любые значения. В предположении, что модули кинетических моментов роторов удовлетворяют неравенствам Нх>Нг>...>Нм, область значений АГ = |К|, в которой определены неподвижные состояния, описывается неравенством

р<К<Я; я = Р = Н,-^Нк. (104)

к=I к=2

При р < 0 неравенство (104) принимает вид К < Я .

Стационарные движения с ненулевой угловой скоростью описываются двумя группами решений. Первая из них имеет вид

К = аКе,, Нк=ЗкНке.; а=± 1, к=ЦУ, / = 1,2,3, 5к=±1, (105)

где е), е2, е3 - главные центральные оси инерции системы. Эти стационарные движения соответствуют вращениям несущего тела вокруг главных осей инерции с угловой скоростью

ю^^-ХХЯ^е,.//,.; / = 1,2,3, (106)

к=1

где J. - момент инерции системы относительно главной оси е..

Каждому значению Кф0 соответствует 6-2м решений (105).

Вторая группа стационарных движений описывается решениями

Н Г8.Н] ¡.; }=\~т\ т> 1, Н,=£,ЯА.К7|Г|; к=т+\,Ы,

~„ , ~ ^ т (Ю7)

К=К*(1+Я/|К*|); К*=(Е+//1)_1Н*; Я=

к=т+1 7=1

Здесь каждый из символов <5.(г"=1,Я) может принимать два значения +1 и -1, Е - единичная матрица размера 3x3 , а параметр /л может принимать любое значение от -оо до -ьх>. Угловая скорость несущего тела для стационарных движений (107) определяется

формулой со=-//(Е+;иЛ)-'Н*.

Для решений (107) рамки первых т гироскопов «сложены», т.е. внутренняя рамка ] -го гироскопа параллельна внешней. Общее число стационарных движений вида (107) составляет не менее 2-2^(2^—1), а для достаточно больших значений К достигает 6-2м(2м -V).

Установлено, что решения (107) неустойчивы в вековом смысле, а из решений (105) вековой устойчивостью обладают только следующие стационарные движения:

К = стКе3,Нк=аНке3; сг = ± 1; К>Е, (108)

К=стКс3, Я{=стНхег Н^ =-сгЯ(Ье3; а=±\, к=2Л; К<р , (109) где е3 — ось наибольшего момента инерции системы. В диапазоне значений кинетического момента р<К<Я, где нет устойчивых в вековом смысле стационарных движений с ненулевой угловой скоростью, устойчивы по части переменных (по переменным та, у, 6) неподвижные состояния системы. Все они в совокупности образуют устойчивое инвариантное множество.

При наличии диссипации в осях рамок всех гироскопов решения (108), (109) и неподвижные состояния условно асимптотически устойчивы, а все остальные стационарные движения системы неустойчивы по Ляпунову. При наличии диссипации система асимптотически приходит в состояние покоя из любого начального состо-

39

яния, для которого величина К суммарного кинетического момента системы принадлежит отрезку [р, ß] возможных значений величины |Н| суммарного кинетического момента роторов.

В заключении перечислены основные результаты и выводы диссертационной работы.

Основные выводы диссертации:

1. Разработанные в диссертации теоретические положения дают в сочетании с теоремой Барбашина-Красовского сравнительно простые и эффективные алгоритмы исследования устойчивости стационарных движений гиросиловых систем.

2. Использование силовых гироскопов в системах пассивной гироскопической и гравитационно-гироскопической стабилизации дает возможность соответствующим подбором параметров получить более широкий диапазон устойчивых стационарных движений, чем для твердого тела, а диссипацией в осях рамок обеспечить их асимптотическую устойчивость.

3. Система из двух и более гиродинов с диссипаций в осях рамок обеспечивает в определенном диапазоне начальных условий гашение угловой скорости космического аппарата до нулевого значения без использования активной системы управления.

Основные результаты диссертации опубликованы в работах

1. Амелъкин Н.И. О движениях твердого тела, содержащего двухстепенные силовые гироскопы с диссипацией в осях рамок // Изв. РАН. МТТ. 2006. №4. с. 19 - 30.

2. Амелъкин Н.И. Анализ устойчивости стационарных вращений твердого тела, несущего двухстепенные силовые гироскопы с диссипацией в осях подвеса рамок // Изв. РАН. МТТ. 2007. №4. С. 26-40.

3. Амелъкин Н.И. О предельных движениях гиросиловой системы с внутренней диссипацией в однородном поле тяжести // Изв. РАН. МТТ. 2008. №3. С. 23 - 32.

4. Амелъкин Н.И. Об устойчивых стационарных движениях свободного твердого тела с двухстепенным силовым гироскопом, ось прецессии которого параллельна главной плоскости инерции // Изв. РАН. МТТ. 2009. №1. С. 3 - 16.

5. Амелькин Н.И. О стационарных движениях спутника с двухстепенным силовым гироскопом в центральном гравитационном поле и их устойчивости // ПММ. 2009. №2. С. 236 - 249.

6. Амелькин Н.И. Анализ устойчивости равновесий спутника, несущего двухстепенной силовой гироскоп с диссипацией в оси рамки // ПММ, 2010. Т. 74. № 4. С. 567-581.

7. Амелькин Н.И. О стационарных движениях твердого тела с двухстепенным силовым гироскопом при произвольном расположении оси прецессии гироскопа в несущем теле // Изв. РАН. МТТ, 2010. №5. С. 5-18.

8. Амелькин Н.И. О равновесиях и устойчивости динамически симметричного спутника с двухстепенным силовым гироскопом //ПММ, 2010. Т. 74. № 5. С. 718 -733.

9. Амелькин Н.И. Системы гиродинов с дополнительной степенью свободы. // Космич. исслед. 2003. том 41. №3. С.295 - 300.

10. Амелькин Н.И. Показатели локальной управляемости систем силовых гироскопов //Космич. исслед. 2004. Т.42. №4. С. 424 -430.

11 .Амелькин Н.И. Об оценках уходов уравновешенного гироскопа и погрешностях формулы Магнуса // Изв. РАН. МТТ. 2009. Т.6. С. 9 - 20.

12. Амелькин Н.И. Об асимптотических свойствах движений спутников в центральном поле, обусловленных внутренней диссипацией // ПММ, 2011. Т. 75. № 2. С. 204-223.

13. Амелькин Н.И. О свойствах стационарных движений твердого тела, несущего систему двухстепенных силовых гироскопов // ПММ, 2011. Т. 75. №3. С. 355-369.

14.Амелькин Н.И. О стационарных движениях твердого тела, несущего трехстепенные силовые гироскопы, и их устойчивости // Изв. РАН. МТТ, 2011. №3. С. 3-17.

15. Амелькин Н.И. Анализ устойчивости стационарных движений пассивных гиросиловых систем // IX Всероссийский съезд по теоретической и прикладной механике (Нижний Новгород, 2329 августа 2006 г). Аннотации докладов. Т. 1, стр. 13.

16.Амелькин Н.И. Об устойчивости стационарных движений системы связанных тел, содержащих вращающиеся роторы, в центральном поле // Международный конгресс Нелинейный динамический анализ -2007. Тезисы докладов, Санкт-Петербург, Россия, 4 - 8 июня 2007. С. 185.

17.Амелькин Н.И. О предельных движениях гиросиловой системы с внутренней диссипацией в однородном поле тяжести и их устойчивости. // X международный семинар «Устойчивость и колебания нелинейных систем управления». Москва. ИПУ РАН. 3-6 июня 2008г. Тезисы докладов. С. 20-21.

18.Амелькин Н.И. Об устойчивости стационарных вращений твердого тела с двухстепенным силовым гироскопом // «Устойчивость, управление и динамика твёрдого тела» / Тезисы докладов

X Международной конференции (5-10 июня 2008 года). - Донецк: Институт прикладной математики и механики НАНУ, 2008. С. 6-7.

19.Амелькин Н.И. Асимптотические свойства движений спутников в центральном гравитационном поле в неограниченной задаче //

XI международная конференция «Устойчивость и колебания нелинейных систем управления». Москва. ИПУ РАН. 1-4 июня 2010г. Тезисы докладов. С. 15-16.

20.Амелькин Н.И. О стационарных движениях твердого тела, несущего силовые гироскопы, и их устойчивости // X Крымская международная математическая школа «Метод функций Ляпунова и его приложения» (MFL-2010). Крым, Алушта, 13—18 сентября 2010 г. Тезисы доклада. С. 9.

21. Амелькин Н.И. Анализ свойств устойчивости пассивных гиро-силовых систем ориентации в рамках неограниченной постановки задачи // Материалы XXX академических чтений по космонавтике. Тезисы доклада. Москва, январь 2006 г. С. 94-95.

22.Амелькин Н.И. О сравнении двух подходов к анализу устойчивости стационарных движений гиросиловой системы // Труды XLIX научной конференции МФТИ, 2006 г. С. 24-25.

Амелькин Николай Иванович

СТАЦИОНАРНЫЕ ДВИЖЕНИЯ ТВЕРДОГО ТЕЛА, НЕСУЩЕГО СИЛОВЫЕ ГИРОСКОПЫ, И ИХ УСТОЙЧИВОСТЬ

Подписано в печать 20.06.2011. Формат 60 х 84 '/16. Усл. печ. л. 2,3. Тираж 100 экз. Заказ № 57

Государственное образовательное учреждение высшего профессионального образования «Московский физико-технический институт (государственный университет)» 141700, Московская обл., г. Долгопрудный, Институтский пер., 9

Издательский сектор оперативной полиграфии 141700, Моск. обл., г. Долгопрудный, Институтский пер., 9

Государственное образовательное учреждение высшего профессионального образования Московский физико-технический институт (государственный университет)

На правах рукописи

05201151951

Амелькин Николай Иванович

СТАЦИОНАРНЫЕ ДВИЖЕНИЯ ТВЕРДОГО ТЕЛА, НЕСУЩЕГО СИЛОВЫЕ ГИРОСКОПЫ, И ИХ УСТОЙЧИВОСТЬ

Специальность 01.02.01 - теоретическая механика

Диссертация на соискание ученой степени доктора физико-математических наук

Москва-2011

Оглавление

стр.

Введение 4

Глава 1. Стационарные движения. Метод исследования 12

1.1 Постановка задачи 12

1.2 Тензор инерции системы и его производные. Кинетический момент

и кинетическая энергия 14

1.3 Уравнения движения системы 18

1.4 Стационарные движения в однородном поле 20

1.5. Метод исследования устойчивости стационарных движений^

в однородном поле 27

1.6 Обобщения на случай статически неуравновешенных гироскопов 33'>

1.7 Стационарные конфигурации произвольной системы связанных

тел в однородном внешнем поле 44'

1.8 Стационарные движения в центральном гравитационном поле

в рамках ограниченной круговой задачи 47

1.9 Равновесные конфигурации системы связанных тел на круговой орбите 55

Глава 2. Предельные движения.системы с двухстепенными -

гироскопами »при наличии диссипации в осях рамок 60

2.1 Предельные движения гиростатов в однородном поле 60

2.2 О предельных движениях систем с частичной диссипацией

в окрестности положения равновесия 79

2.3 Анализ предельных движений спутника в центральном гравитационном поле на круговой орбите 85

Глава 3. Анализ стационарных движений гиростатов

в,однородном поле 89

3.1 Параметрическое представление множества стационарных движений системы с одним двухстепенным гироскопом 89

3.2 Стационарные движения при параллельном расположении

оси прецессии гироскопа к главной оси инерции системы 95

3.3 Стационарные движения при параллельном расположении оси прецессии гироскопа к главной плоскости инерции системы 102

3.4 Стационарные движения системы с динамической симметрией 117

3.5 Стационарные движения при произвольном расположении оси прецессии гироскопа в несущем теле 122

3.6 Стационарные движения твердого тела, несущего систему двухстепенных силовых гироскопов 139

Глава 4. Стационарные движения спутника, несущего

двухстепенной силовой гироскоп, на круговой орбите 161

4.1 Стационарные движения при параллельном расположении

оси прецессии гироскопа к главной оси инерции спутника 161

4.2 Стационарные движения динамически симметричного спутника

с двухстепенным силовым гироскопом на круговой орбите 189

Глава 5. Стационарные движения твердого тела, несущего

трехстепенные силовые гироскопы, в однородном поле 213

5.1 Постановка задачи 213

5.2 Уравнения движения 215

5.3 Стационарные движения 217

5.4 Анализ вековой устойчивости 224

5.5 Анализ устойчивости по Ляпунову 232

Заключение 235

Литература 240

Введение

Системы управления и стабилизации вращательного движения твердого тела, использующие в качестве исполнительных элементов силовые гироскопы или роторы (маховики), разделяются на активные и пассивные. К активным относятся системы, в которых движение гироскопов* или маховиков'относительно несущего тела является управляемым. В* таких системах управление, обеспечивающее требуемое вращательное движение несущего тела, строится на принципах обратной связи и осуществляется с помощью активных моментных устройств, устанавливаемых на осях рамок гироскопов или осях маховиков. В качестве примеров активных систем можно указать системы* ориентации' космических станций, в которых используются гиродины, и системы ориентации геостационарных спутников, использующие управляемые маховики.

В пассивных гиросиловых системах взаимодействие между несущим телом и гироскопами обеспечивается только за счет реакций связей и устанавливаемых в осях рамок гироскопов пассивных моментных устройств, например, пружин и демпферов. В > принципах их работы помимо гироскопических свойств вращающихся тел используются также свойства внешней среды (моменты гравитационных и аэродинамических сил, магнитное поле и др.). Гироскопы (или роторы) применяются в таких системах с целью получить дополнительные восстанавливающие моменты, а также новые стационарные движения, отличные от стационарных движений твердого тела. Кроме того, установка демпферов« на осях рамок гироскопов обеспечивает во многих случаях асимптотические свойства системы, что особенногважно для практики.

В диссертации рассматриваются задачи, имеющие отношение к пассивным способам гироскопической и гравитационно-гироскопической стабилизации вращательного движения твердого тела.

Диссертация посвящена исследованию стационарных движений твердого тела, несущего силовые гироскопы, в однородном внешнем поле и в центральном гравитационном поле. Под «однородностью» внешнего поля подразумевается допущение, что в каждой точке области пространства, занимаемой системой, величина и направление напряженности поля одинаковы; при этом величина и направление вектора напряженности может меняться со временем произвольным образом. Используемое в дальнейшем свойство такого поля состоит в том, что^оно не влияет на движение системы относительно базиса Кенига. Для> спутников) допущение об« однородности внешнего поля* в. указанном выше смысле эквивалентно игнорированию гравитационных, аэродинамических и других моментов'

В случае однородного' внешнего поля под стационарными' движениями будут пониматься движения; для которых вектор угловой скорости несущего- тела постоянен, а рамки гироскопов неподвижны относительно несущего' тела. В' случае центрального гравитационного поля под стационарными подразумеваются» движения, для которых центр масс системы движется равномерно по круговой орбите вокруг центра притяжения, рамки гироскопов неподвижны, относительно несущего тела, а угловая скорость несущего тела совпадает с угловой скоростью орбитального базиса (положения равновесия системы относительно орбитального базиса).

Для исследования стационарных движений применяются, методы теории, дифференциальных уравнений и теории устойчивости, в том числе первая теорема Рауса и ее модификации [22, 35, 36, 39], теоремы-Ляпунова и Четаева [28, 42], теоремы Барбашина и Красовского [3, 25, 41, 47].

К настоящему времени достаточно подробно изучена задача о стационарных движениях твердого тела, несущего осесимметричные силовые роторы с постоянной скоростью собственного вращения. Для случая однородного внешнего поля полное аналитическое решение этой задачи для всевозможных вариантов установки оси ротора в несущем теле приведено в [11].

Большое число работ посвящено исследованию стационарных движений спутника, несущего силовой ротор, в центральном гравитационном поле в рамках ограниченной круговой задачи ([39, 54, 68, 69; 82, 83, 84] и др.). Подробное аналитическое решение «прямой» задачи, в которой определяется зависимость стационарных движений от величины кинетического момента ротора и главных центральных моментов инерции спутника, получено только для случаев, когда ось ротора параллельна главной« оси, либо главной плоскости!инерции1 спутника [54, 82, 83].

Для; спутника с силовым ротором проведен также сравнительный анализ стационарных движений-в ограниченной и неограниченной задаче [39, 84]. Показано, что уравнения, определяющие стационарные движения, и достаточные условия устойчивости в ограниченной и1 неограниченной задаче совпадают с точностью до малых величин, равных квадрату отношения ^характерного размера спутника,к радиусу орбиты.

Стационарные движения твердого тела, несущего силовые гироскопы, до последнего времени были изучены в меньшей степени. В опубликованных работах на эту тему в основном исследовались отдельные стационарные движения систем с конкретной схемой установки* гироскопов в несущем теле, либо анализ был ограничен узким диапазоном значении параметров системы.

Например, в [38, глава 6] рассматривались полупассивные системы ориентации спутников, построенные с использованием силовых гироскопов. В них перемещения рамок гироскопов относительно несущего тела конструктивно ограничены узким диапазоном углов, а анализ стационарных движений проводился ;в ограниченном диапазоне значений суммарного кинетического момента системы.

В [63, 78] рассматривались системы, в которых установленный на спутнике двухстепенной силовой гироскоп используется для демпфирования нутационных колебаний вращающегося спутника при отсутствии внешних моментов сил. В указанных работах рассматривались конкретные

схемы установки гироскопов в несущем теле и для отдельных стационарных вращений спутника определялись условия устойчивости и оптимальные значения? параметров, при которых достигается максимальное значение степени устойчивости.

В' [74, 75, 77] исследовались стационарные движения?(положения; относительного < равновесия) в центральном гравитационном? поле спутников, несущих два одинаковых, двухстепенйых силовых гироскопа с пружинами в. осяхч рамок; установленные по-'определенной^схеме в.корпусе спутникаи При^этомш [74]; для>исследования^устойчивостшприменялась теорема?Бар-башина—Красовского, что» позволило?: получить, условия? асимптотической« устойчивости более простым способом, чем с: помощью - критерия Рауса — Г'урвица. •

Стационарные движения, спутника, несущего1 одшI трехстепенной силовой- гироскоп в. кардановом подвесе, в центральном гравитационном поле в рамках ограниченной круговойзадачи*исследовались в: [45; 46; 86]1 Подробное решение задачи приведено в [86]. Здесь рассматривался5случай; когда система «корпус спутника — гироскоп» не" обладает свойством гиростата, но предполагалось, что гироскоп статически уравновешен, а ось. внешней рамки: гироскопа параллельна одной« из главных центральных осей инерции корпуса, спутника.

В диссертации основное внимание* уделено исследованию стационарных движений'твердого тела; несущего; двухстепенные силовые гироскопы с постоянной скоростью вращения роторов и неограниченным диапазоном-углов; прецессии:, Рассмотрена и решена также задача о стационарных движениях твердого, тела, несущего; произвольное число трехстепенных силовых гироскопов в кардановых подвесах, в однородном внешнем поле (для случая, когда система «несущее тело - гироскопы» обладает свойством гиростата).

В рассматриваемых задачах стационарные движения системы, представляют собой множества; зависящие от таких параметров, как значения

кинетических моментов роторов, углы установки осей рамок гироскопов в несущем теле, значения главных центральные моменты инерции системы. В случае однородного внешнего поля стационарные движения зависят также от величины (модуля) суммарного кинетического момента системы, которая является первым интегралом уравнений движения.

При, исследовании конкретных систем основной целью являлся поиск решений «прямой» задачи, т.е. определение зависимости переменных, характеризующих стационарные движения, от указанных выше параметров и< анализ свойств устойчивости'найденных решений:

Стационарные движения системы воднородномвнешнем< поле определялись на- основании доказанной в главе 1 теоремы 1.1 как стационарные точки? усеченной' обобщенной энергии на многообразии, задаваемом усеченным интегралом кинетического момента. Стационарные движения в центральном поле в рамках ограниченной круговой »задачи -находились как стационарные точки измененной потенциальной энергии:

Анализ устойчивости стационарных движений проводился' в два этапа. На первом этапе определялся, характер вековой устойчивости. На втором этапе в предположении, что в осях рамок гироскопов действуют дис-сипативные силы, проводился детальный анализ устойчивости по Ляпунову с применением теоремы Барбашина-Красовского. Предварительно были получены конструктивные способы проверки условия теоремы Барба-шина - Красовского об отсутствии целых траекторий движения1 без диссипации в малых окрестностях стационарных движений.

В диссертации рассмотрены также системы, составленные из произвольной совокупности связанных тел без гироскопов, и проведен качественный анализ конфигураций, соответствующих стационарным движениям таких систем в однородном внешнем поле и в центральном гравитационном поле.

Диссертация состоит из введения, пяти глав, заключения и списка литературы. Во введении обосновывается актуальность темы, дается краткий исторический обзор и краткое содержание диссертации.

Глава 1 посвящена методам исследования стационарных движений твердого тела, несущего N двухстепенных силовых гироскопов с постоянными скоростями вращения роторов. Для системы, не обладающей свойст-вом< гиростата, выводятся формулы, определяющие зависимость, тензора инерции и его производных по1 углам! прецессии гироскопов^от текущей конфигурации системы (разд.Л .2, 1.6). Выводятся уравнения движенияI относительно базиса Кенига, при произвольных действующих силах (разд. 1.3) и уравнения, определяющие стационарные движения системы в однородном внешнем.поле (разд. 1.4, 1.6). Доказывается теорема 1.1 о взаимном соответствии между стационарными движениями системы в (однородном внешнем поле- и< условно-стационарными точками усеченной- обобщенной? энергии» на многообразии, определяемом усеченным1 интегралом кинетического момента (разд. 1.4', 1.6). Доказывается теорема 1.2, согласно которой характер вековой устойчивости определяется характером стационарности усеченной обобщенной энергии на многообразии, определяемом усеченным интегралом кинетического момента (разд. 1.5, 1.6).

В* разд. 1.7 рассматривается произвольная система связанных тел и проводится качественный анализ ее стационарных конфигураций в однородном внешнем поле.

В разд. 1.8 выводятся уравнения движения спутника, несущего двухстепенные силовые гироскопы, в центральном гравитационном поле в рамках ограниченной круговой задачи, и уравнения, определяющие стационарные движения спутника (относительные положенияфавновесия).

В разд. 1.9 проводится качественный анализ равновесных конфигураций произвольной системы связанных тел без гироскопов на круговой орбите.

В главе 2 излагается методика исследования устойчивости по Ляпунову стационарных движений системы, при наличии диссипации в осях рамок гироскопов. Методика основана на теореме Барбашина-Красовского. Приводятся конструктивные способы проверки условия 1° этой теоремы об изолированности: стационарных движений от других предельных движений системы (целых траекторий движения без диссипации).

В разд. 21Г проводится анализ предельных, движений гиростатов в' однородном- внешнем поле. Доказывается? теорема;2Л, согласно которой указанное:выше-условие 1 ° сводится к изолированности рассматриваемого стационарного движения от других стационарных движений на многообразии, определяемом? усеченным интегралом кинетического момента системы. ■ '

Дляг исследования; предельных движений гиростатов в центральном гравитационном поле на круговой орбите: используются результаты работ Пожарицкого Г.К. и Лилова Л.К, сформулированные теоремой. 2.2. Эта теорема определяет для« систем частичной диссипаций достаточные условия изолированности стационарных, движений от других предельных движений на! основе усеченных уравнений линейного приближения (разд. 2.2).

В разд. 2.3' выводятся линеаризованные в окрестности положений равновесия« уравнения движения спутника, несущего двухстепенные силовые гироскопы, на круговой орбите и усеченные уравнения: линейного приближения, фигурирующие в теореме 2.2. Приводится теорема 2.4, с помощью которой определяется характер устойчивости по Ляпунову по-ложений.равновесия спутника при наличии диссипации в осях рамок гироскопов. . '

В главе 3 приводятся решения «прямой» задачи о стационарных движениях твердого тела с двухстепенными силовыми гироскопами в однородном внешнемшоле. В разд. '3.1-3.5 исследуются стационарные движения системы с одним гироскопом для различных вариантов расположения оси рамки гироскопа в несущем теле. В'разд. 3.6 изучаются стацио-

нарные движения системы, содержащей произвольное число N гироскопов, установленных по коллинеарной схеме в несущем теле. Рассмотрен также пример системы с двумя гироскопами, оси рамок которых не параллельны друг другу (взаимно ортогональны и направлены по главным осям инерции системы).

В главе 4 исследуются стационарные движения спутника с двухстепенным силовым гироскопом в центральном гравитационном поле в рамках ограниченной круговой задачи. Приводятся решения «прямой» задачи для случая, когда ось рамки гироскопа параллельна одной из главных центральных осей инерции спутника (разд. 4.1), и случая динамичес