О равновесиях и устойчивости спутника с системой двухстепенных силовых гироскопов в центральном гравитационном поле Текст научной статьи по специальности «Физика»

УДК 531.36

Н.И. Амелъкин1, А. В. Зыков2

1 Московский физико-технический институт (государственный университет)

2ОАО РКК «Энергия»

Ои и

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

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

Ключевые слова: спутник, силовой гироскоп, равновесие, устойчивость, асимптотическая устойчивость.

1. Введение

Вопросам применения силовых гироскопов в системах пассивной ориентации спутников посвящены работы [1-5]. Подробное исследование свойств положений относительного равновесия спутника с одним двухстепенным силовым гироскопом в центральном гравитационном поле проведено в работах [3-5]. Получено полное аналитическое решение задачи о положениях равновесия спутника с одним двухстепенным гироскопом для случая, когда ось прецессии гироскопа параллельна одной их главных осей инерции спутника [3], и для динамически симметричного спутника при произвольном расположении оси прецессии гироскопа в несущем теле [5].

Метод исследования устойчивости по Ляпунову положений равновесия спутника при наличии диссипации в осях рамок гироскопов изложен в работе [4]. Согласно этому методу вопрос об асимптотической устойчивости или неустойчивости положений равновесия систем с частичной диссипацией решается на основе результатов анализа вековой устойчивости и усеченных уравнений линейного приближения.

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

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

2. О положениях равновесия спутника и их устойчивости

Рассмотрим спутник, состоящий из несущего твердого тела и N двухстепенных силовых гироскопов. Обозначим через вк единичные векторы, указывающие фиксированные направления осей прецессии гироскопов в корпусе спутника. Текущие положения осей роторов будем задавать единичными векторами Нк (хк), где Хк — углы прецессии. Предполагается, что оси роторов ортогональны осям рамок, т.е. «^Г Н = 0, а скорости собственного вращения роторов постоянны. Тогда кинетические моменты собственного вращения роторов будут зависеть только от углов прецессии гироскопов:

Hк = Hк hk (хк), iifc = const > 0.

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

Исследуется вращательное движение спутника в центральном гравитационном поле в рамках ограниченной круговой задачи, т.е. в предположении, что центр масс спутника движется по кеплеровой круговой орбите. Орбитальный базис определяется взаимно ортогональными единичными векторами г, пит = п X г, направленными по радиусу орбиты, по нормали к плоскости орбиты и по касательной к орбите соответственно. Вектор угловой скорости орбитального базиса ^о = Шо п постоянен в инерциальном базисе.

Положение спутника относительно орбитального базиса определяется 3 + N независимыми переменными. Три из них задают ориентацию корпуса спутника относительно орбитального базиса, которая однозначно определяется через компоненты векторов г и п в связанном с корпусом базисе, а остальные переменные — углы хь поворота рамок гироскопов.

Уравнения движения спутника, несущего систему двухстепенных гироскопов, выводятся по аналогии с уравнениями движения спутника с одним двухстепенным гироскопом [4] и имеют следующий вид:

N

J (й + ^о X ш) + 1ь Хь + вь х Йь Xь)+

ь=1

N

+ ( ^0 + ш) X ^(ш + Эь 1ь Xь + Йь) =3 ^0 г X Jr,

ь=1

(1)

1ь (хь + (й + ^о X ш)) — (ш + ^о)Т (эь X Й ь ) = Мь, к = 1, ...,Ы. (2)

Здесь ш — угловая скорость корпуса спутника относительно орбитального базиса, ^о = ^о п — угловая скорость орбитального базиса, правая часть уравнений (1) представляет собой момент действующих на спутник гравитационных сил, Мь — моменты внутренних дисси-пативных и потенциальных сил в осях рамок гироскопов, 1ь - моменты инерции гироскопов относительно осей вь-

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

г = —и X г, п = —и X п. (3)

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

[3]

W = —пТ J п/2 — птЙ + 3 гт J г/2 + П(х), (4)

представляющей собой измененную потенциальную энергию, поделенную на величину ш2. Здесь

N

Й = ^ Й ь; Йь(хь) = Нь Кь(хь); Нь = Йь/^о (5)

ь=1

ции, И(ж) = П(ж)/^о — «приведенная» энергия действующих в осях рамок внутренних

«приведенный» суммарный кинетическии момент роторов с размерностью момента инер-

потенциальных сил.

Поскольку векторы п и г связаны равенствами

,т,__1 _ 1 „т

п п = 1, г г = 1, п г = 0, (6)

то, используя функцию Лагранжа с множителями

Ь = № + А1 пт п/2 + \2 гтг/2 + Аз птг, (7)

получим для положений равновесия систему уравнений

дЬ/дп = —3 п — Н + А1 п + Аз г = 0, (8)

дЬ/д г = 33 г + А2 г + Аз п = 0, (9)

дЬ/дхк = -Нк пт (8к х Нк )+дП/дхк = 0 , к = 1, 2 ,...,Ы. (10)

Уравнение (10) описывает условия равновесия гироскопов относительно корпуса спутника, а уравнения (8), (9) сводятся к одному векторному уравнению, не содержащему множителей Лагранжа:

п х (3 п + Н)=3 г х 3г. (11)

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

птйп = 0, гтйг = 0, гт йп + птйг = 0. (12)

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

йп = и г + ут , йг = —и п + т, (13)

а множители Лагранжа определяются из уравнений (8), (9) формулами

А1 = пт3 п + птН, А2 = —3гт3 г, Аз = —3пт3 г. (14)

При учете этих соотношений и вытекающего из (9) равенства тт3 г = 0 квадратичная форма й2Ь на множестве (12) сводится к следующей квадратичной форме от N+3 независимых переменных и, V, и> и хк = йхк:

Ф = й2Ь | (12) = (4 (пт3п — гт3г) + птН) и2 + (пт3п — тт3 т + птН) у2+ +3(тт3т — гт3 г) ад 2 — 6 пт3 т и,ш — 6пт3 г V -ш+

N

+

к=1

N д2П

птН к — 2(г и + т у)т(вк х Нк) гк + ^ дсс дх ^

3 = 1 3

(15)

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

3. Положения равновесия спутника с ко л линеарной системой двухстепенных гироскопов

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

вк = к = 1, 2,...,И. (16)

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

Положим для определенности, что модули приведенных кинетических моментов роторов удовлетворяют неравенствам

Н1 ^Н2 ^ ... ^НN. (17)

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

N N

К = ^Нк,р = Н\ — ^Нк. (18)

к=1 к=2

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

р < |Н| < П. (19)

Отметим, что для рассматриваемой системы р > 0 только в том случае, когда в системе есть гироскоп, для которого модуль кинетического момента ротора превышает сумму модулей кинетических моментов роторов всех остальных гироскопов. В противном случае множество всевозможных значений вектора Н представляет собой круг радиуса К, а неравенство (19) принимает вид |Н| ^ К.

В рассматриваемом случае уравнения (10) принимают вид

пт (в х Нк) = 0; к = 1, 2 ,...,М, (20)

а система (11), (20) имеет две группы решений. Для первой группы оси всех роторов взаимно параллельны, т.е.

Нх = Н, Нк = ак Н; к = 2,... ак = ±1, (21)

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

п х (I п + Н Н) = 3 г х Зг, пт(з х Н) = 0, зтН = 0, (22)

где

N

Н = Нх + ^ ак Нк. (23)

к=2

Здесь каждый из символов ак может принимать два значения +1 и —1, вследствие чего величина Н может принимать 2м-1 фиксированных значения.

Отметим, что система (22) описывает всевозможные положения равновесия спутника с одним гироскопом, имеющим кинетический момент ротора Н = Н Н(х) [4]. Поэтому

оМ— 1

задача определения первой группы решении сводится к решению 2" 1 задач о положениях равновесия спутника с одним гироскопом. При этом каждая из этих задач характеризуется одинаковым направлением оси прецессии гироскопа в в корпусе спутника, а различаются только значения величины кинетического момента ротора Н.

Вторая группа решений системы (11), (20) характеризуется условием п х в = 0 (ось ^ направлена по нормали к плоскости орбиты) и находится из системы уравнений

п = а в, в х З 8 + а Н) = 3 г х Зг, зтН = 0; а = ±1. (24)

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

зт(г х Зг) = 0, зтг = 0, гтг = 1. (25)

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

Второе уравнение задает плоскость, ортогональную оси s, а третье уравнение — единичную сферу.

Подробный анализ показал, что система (25) всегда имеет только четыре решения:

г = ±ri(s), г = ±r2(s), (26)

где решения ri(s) и г2 (s) определяются в явном виде из квадратного уравнения. Наиболее простой вид эти решения имеют в случае, когда ось s параллельна главной плоскости инерции спутника. Например, если

s = —е2 sin ,0 + е3 cos 0, (27)

то

r1(s) = е1, r2(s) = е3 sin, + е2 cos,. (28)

Решения H( s) находятся в явном виде из второго уравнения системы (24). Умножив это уравнение векторно на s, получим

H = a [s х (s х J s — 3 r x Jr)j = a [s (sTJ s) — J s + 3 r(sTJr)].

Подставляя сюда решения (26), получим

H(1) = a [s (sTJ s) — J s + 3 ri(sTJn)], H(2) =a [s (sTJ s) — J s + 3 r2(sTJr2)]. (29)

Условия существования решений (26), (29) системы (24) определяются вытекающими из (19) неравенствами

р2 ^ H1 ^ R , р ^ H2 В случае (13) решения (29) и условия (30) имеют следующий вид:

H(1) = a Г2 (В — С) sin 0 cos 0; р2 < (В — С)2 sin2 0 cos2 0 ^R

(2)

р2 < H? < R2, р2 < H2 < R2. (30)

H(2) = 4a Г2 (В — С) sin 0 cos 0; р2 < 16(В — С)2 sin2 0 cos2 0 < R2

В зависимости от значений параметров спутника (главных центральных моментов инерции, расположения оси в, модулей кинетических моментов роторов) число равновесных ориентаций спутника {п, г}, определяемых системой (24), равно либо 0, либо 4, либо 8.

Отметим, что решения (29) дают значения для суммарного кинетического момента роторов Н(в). При этом значения кинетических моментов роторов отдельных гироскопов Н к(в), т.е. углы прецессии хк( в), определяются неоднозначно. В случае N = 2 каждому значению Н( 8) из интервала р < |Н| < К соответствует две комбинации углов Х1( 8), х2( 8), а при N ^ 3 каждому значению Н( 8) соответствует (^ — 2)-мерное множество значений вектора х = (х1, х2,..., х^).

Исследуем характер вековой устойчивости положений равновесия спутника. Начнем с первой группы решений, определяемых уравнениями (21), (22). Для них квадратичная форма (15) записывается в виде

Ф = (4 (nTJn — rTJr) + nTH) u2 + (nTJn — ttJ t + nTH) v2+ +3 (tt J T — rT J r)w2 — 6 nT J T uw — 6 nT J r V w+

N N

+nTakHkz2k — 2(ru + tv)t(s x h) ^ akHkzk. k=1 k=1

(31)

Здесь а1 = 1, ак = ±1; к = 2,..., N.

Определим сначала необходимые условия вековой устойчивости рассматриваемых решений. Покажем, во-первых, что для вековой устойчивости необходимо выполнение неравенства

пт к = 0. (32)

Действительно, из второго уравнения системы (22) следует, что вектор п ортогонален вектору s х h. Поэтому выполняется одно из неравенств гт(s х h) = 0 либо тт(s х h) = 0. Пусть для определенности выполняется первое из указанных неравенств. Тогда получаем, что при птh = 0 на множестве v = w = zs = 0 s = k, квадратичная форма (31) принимает вид

Ф(и, zk) = (4 (птJn - гт Jr) + птH) и 2 - 2г т(s х h) ak Hk uzk.

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

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

ak = 1; к = 1,...,N, hk = h; к = 1,...,N, птh > 0, (33)

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

N

Т1к,

H = R h, R = Е Нк, (34)

к=1

а квадратичная форма (31) принимает вид

Ф=(4 (птJn - гт Jr) + птH) и 2 + (пт Jn - тт J т + птH) v2 + 3 (ттJ т - гтJ г) w 2-

N N

-6 пт J т uw - 6 пт J г vw + птh ^ Нк - 2(r и + т v)T(s х h)^^ Нк zk.

к=1 к=1

Введем переменную z соотношением

N

(35)

Rz = Hk zk. (36)

к=1

Из неравенства (32) и неравенства Коши-Буняковского

2

22

N / N \ / / N \

Е н* z2 > Е Нк А Е

к=1 \к=1 / / \к=1 /

_ к г2к ^^ Ик хк\ / ( > Як )= Rz¿ (37)

к=1

получим

Ф ^ (4 (птЗп - гтЗг) + птН) ик + (птЗп - тт3 т + птН) Vк+

гр ггт с\ ггт ггт

+3 (т 3 т — г 3 г) ад — 6 п 3 т и,ш — 6 п 3 г (38)

+ пт к Rzк — 2(г и + т у)т (в х h)Rz = Фь

При этом на множестве переменных гк = г ^^^^^^^^таая форма Ф (35) совпадает с квадратичной формой Ф1 (38). Отсюда следует, что вековой устойчивостью обладают только те из рассматриваемых решений, для которых положительно определена квадратичная форма (37). Эта форма зависит от четырех переменных и, V, ад, г и определяет характер вековой устойчивости положений равновесия спутника с одним гироскопом, имеющим кинетический момент ротора Н = Rh.

Исследуем характер вековой устойчивости второй группы решений, определяемых уравнениями (24). Поскольку для этих решений п х в = 0, то птН^ = 0 Vк, и, кроме того, вектор п ортогонален каждому вектору 8 х Н&. Выше было показано, что при этих условиях

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

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

N

зовании одного гироскопа с кинетическим моментом ротора К = ^ Нк и осью прецессии е.

к=1

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

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

4. Заключение

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

Работа выполнена при финансовой поддержке Правительства Российской Федерации в

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

за № 2014/120. НИР № 2583.

Литература

1. Сарычев В.А., Лукапип К.В., Мирер С.А. Оптимальные параметры гравитационно-гироскопических систем ориентации спутников // Космич. исслед. — 1975. — Т. 13, № 3. - С. 311-321.

2. Сарычев В.А. Вопросы ориентации искусственных спутников // Итоги науки и техники. Сер. Исследование космического пространства. Т. 11. — М.: ВИНИТИ, 1978. — 223 с.

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

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

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

Поступим в редакцию 29.11.2012.