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

РАКЕТНО-КОСМИЧЕСКАЯ ТЕХНИКА

УДК 531.381:629.78

DOI 10.26732/2618-7957-2019-2-73-86

СОВРЕМЕННЫЕ АЛГОРИТМЫ АКТИВНОЙ МАГНИТНОЙ ОРИЕНТАЦИИ СПУТНИКОВ

М. Ю. Овчинников Д. С. Ролдугин

Федеральный исследовательский центр Институт прикладной математики им. М. В. Келдыша РАН, г. Москва, Российская Федерация

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

Ключевые слова: ориентация спутника, угловое движение, магнитная ориентация, магнитное управление, магнитная система ориентации, магнитная катушка.

Введение

Магнитные системы управления используются на спутниках с начала эры освоения космического пространства. Первым спутником с пассивной магнитной системой стал Transit 1B [1], запущенный 13 апреля 1960 года. Первая активная магнитная система использовалась на аппарате Tiros II [2], запущенном 23 ноября 1960 года. Реконструкция движения по данным измерений магнитного поля была впервые проведена для Спутника-3 [3], запущенного 15 мая 1958 года.

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

Н ovchinni@keldysh.ru © Овчинников М. Ю., Ролдугин Д. С., 2019

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

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

ТЕХНОЛОГИИЕ

гг

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

Из векторного произведения видно, что управляющий момент не создается вдоль вектора геомагнитной индукции. Система может показаться неуправляемой. Однако направление, вдоль которого недоступен управляющий момент, меняется в пространстве при движении спутника по орбите. В результате все направления с течением времени оказываются доступны. Строго управляемость спутника с магнитной системой ориента-74 ции показана в [4; 5]. В этих работах используется предположение о периодическом изменении вектора геомагнитной индукции. Это предположение вполне оправданно, т. к. геомагнитное поле обычно моделируется полем диполя на этапах предварительных исследований и синтеза управления. Более популярна модель наклонного диполя, указывающего примерно в магнитные полюса Земли. Изменение вектора индукции определяется движением аппарата по орбите и вращением Земли. Поскольку второй фактор приводит к заметно более медленному изменению, чем первый, им зачастую пренебрегают, рассматривая диполь, антипараллельный оси вращения Земли. В этом случае вектор индукции изменяется периодически с течением времени (по мере движения аппарата по круговой орбите). Применимость других моделей в задачах ориентации спутников рассмотрена в [6; 7], квадрупольные и модели большего порядка рассмотрены в [8]. В [9] показана управляемость линейной системы кроме случаев двух соотношений между моментами инерции, в том числе, если спутник не осесимметричный. В [10] представлено доказательство управляемости и предложена стратегия управления для малых отклонений от требуемой ориентации.

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

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

Том з

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

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

Вторая часть посвящена применению магнитной системы совместно с другими исполнительными элементами или с использованием пассивных подходов и свойств внешней среды. В первую очередь рассматривается спутник, стабилизированный вращением. Быстро вращающийся аппарат приобретает свойства гироскопа и сохраняет ориентацию оси вращения в инерциальном пространстве. Внешние возмущения приводят к медленному изменению направления вращения в инерциальном пространстве и возрастанию амплитуды нутационных колебаний. Магнитное управление используется для компенсации нутационных колебаний и ориентации оси вращения в требуемом направлении. Спутник с тангажным маховиком реализует похожий принцип. Ось вращения маховика, имеющего высокую скорость закрутки и/или массу и, как следствие, кинетический момент, ориентируется по нормали к плоскости орбиты. Магнитная система обеспечивает асимптотическую устойчивость такого движения и поворот в плоскости орбиты. Наконец, пассивная гравитационная ориентация может обеспечить важные условия устойчивости орбитальной ориентации аппарата. Описанные подходы, хорошо известные и проверенные на множестве аппаратов, рассматриваются кратко. Более полный обзор соответствующих алгоритмов приведен в [11].

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

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

Основное внимание в обзоре уделено работам последних двух десятилетий. Краткий обзор алгоритмов магнитного управления [12], подробные обзоры основных достижений в области ориентации спутников в первые тридцать лет освоения космического пространства [13; 14], обзор методов пассивной магнитной ориентации [15] и обобщающие результаты в области магнитного управления [16; 17] позволяют включить лишь наиболее важные и актуальные сегодня ранние результаты.

1. Гашение угловой скорости

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

Наиболее популярен алгоритм демпфирования Bdot, впервые опубликованный в [18]. Управление строится из требования отрицательности производной кинетической энергии и имеет вид:

m = kw х b,

где w - абсолютная угловая скорость спутника. Здесь и далее к - коэффициент усиления. Это управление преобразуется в алгоритм Bdot на основании связи производной вектора геомагнитной индукции db/dt в связанной и dB/dt в инерциаль-ной системах координат:

db dB

— =--ш х b.

dt dt

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

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

1 б/Ь IU

m = -к— = -кЪ dt

Управление (2) простое и удобное на предварительном этапе успокоения аппарата. Не применяя сложных методов обработки датчиковой информации, управление обычно реализуется в виде:

т = _к Ь *+1 ~ Ь * .

4+1 4

Для этого достаточно иметь два последовательных измерения магнитометра.

Скорость, до которой управление (2) способно успокоить аппарат, обусловлена отброшенным членом, и чуть меньше удвоенной орбитальной [19]. Управление (1) лишено этого недостатка и обеспечивает экспоненциальное стремление скорости к нулю [20-23] с быстродействием, зависящим от наклонения орбиты [24].

В управлении (1) можно использовать относительную угловую скорость движения спутника в орбитальной системе координат. Управление строится на основе функции Ляпунова, в качестве которой используется интеграл Якоби при движении под действием гравитационного момента [25]. В [26] вместо скалярного коэффициента усиления используется положительно-определенная матрица. Также можно использовать переменный коэффициент усиления для компенсации уменьшения величины скорости и поддержания величины управляющего момента. Во всяком случае, можно использовать две величины, изменяя коэффициент усиления при уменьшении скорости ниже определенного значения [27]. Для дальнейшего выигрыша в быстродействии можно использовать оптимальные законы, например, приведенный в [28]. Однако простота, надежность и быстродействие законов управления (1) и (2), а также обычно мягкие требования на время переходных процессов, делают такие подходы полезными лишь в редких случаях.

2. Использование магнитного управления совместно с другими исполнительными элементами и пассивными подходами

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

75

ТЕХНОЛОГИИ t

76

т

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

Одна из наиболее популярных схем управления быстро вращающимся аппаратом была предложена в [29]. Управление разделено на три основных этапа: гашение нутационных колебаний, раскрутка спутника, переориентация оси вращения. Было предложено релейное управление и соответствующие функции переключения. В [30] рассматривается управление:

m. =i

j

m

-mn

h.(, h.(

e j x

e j x

b )> 0 b )< 0

где mj - дипольныи момент вдоль j-й оси спутника с ортом ej, h - вектор кинетического момента аппарата. Если ось вращения третья, то выбор j = 1 или j = 2 отвечает управлению скоростью закрутки (3). При j = 3 реализуется гашение нутационных колебании.

Три простых закона управления для описанных стадиИ предложены в [31; 32]. Для гашения нутационных колебании применяется управление:

, db m = —k—e,

3 dt 3

соответствующее алгоритму Bdot, реализованному одноИ катушкой. В результате не гасится скорость вдоль оси симметрии. Детально этот алгоритм исследован в [33]. Закрутка обеспечивается управлением:

m = k(b2, -Ьь 0).

Популярный алгоритм закрутки «Y-Thom(p) son» был предложен в [34]. Закрутка вокруг второй оси y задается управлением:

mi = k(2 - <»2re/)sign(b3):

где ю2ге/ - требуемая скорость вращения. Алгоритм переориентации оси симметрии:

m3 = kAh-(e3 х b),

где Ah = hre/ - h - разница между требуемым hre/ и текущим h вектором кинетического момента.

В [35] идея использования ошибки ориентации обобщена как:

Aj = Hre/ - h - Jw, A2 = hre/ - h - Jw.

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

Том з

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

Три простых закона управления предложены в [36]:

т = Ь х [^1(ш1, ©2, 0) + ¿2(0, 0, Из - А3/ + МЩ.

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

В [37] был предложен популярный алгоритм управления углами у, ф (углы у, ф, 9 задают последовательность поворотов 1-2-3):

т2=к1Ь1у-к2Ь2,

и показана асимптотическая устойчивость требуемой ориентации в зависимости от коэффициентов усиления. Таким образом, разработана конструктивная методика подбора коэффициентов, расширенная с помощью теории Флоке и численного моделирования с учетом многих возмущающих факторов. Более общее управление для вращения вокруг третьей оси предложено в [38]:

т3 = К! [А

Функция / должна удовлетворять условиям х/х) > 0 и /(0) = 0. Она задает обобщение управления с обратной связью и в большинстве случаев /х) = х дает удовлетворительный результат. В [39] строится кусочно-постоянное управление, основанное на простых ошибках по углам ориентации и скоростям.

2.2. Спутник с двойным вращением

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

В [18] предложено линеаризованное управление (4):

т2 = АгДф - к2

В [41] показания магнитометра в (3) заменены на данные инфракрасного датчика вертикали. Простое управление поворотом в плоскости орбиты предложено в [40]:

м

0, к sin

^3 у

при этом маховик установлен вдоль второй связанной оси.

В [42] управление (6) преобразовано к виду:

ш = -А1

до этого предложенному для гравитационно ориентируемого аппарата [43]. Это управление, как и линейно-квадратичный регулятор [44], было успешно проверено на аппарате Gurwin. В последней работе возмущающие моменты аппроксимировались простыми периодическими выражениями с орбитальной и удвоенной орбитальной частотами, для геомагнитного поля использовалась модель упрощенного прямого диполя. Результаты построения управления, заметно упрощенного в постановке линеаризованных уравнений с периодическими коэффициентами, оказались применимыми для более полной модели движения аппарата.

В [45] приведено сравнение предложенных ранее законов управления и рассмотрено их обобщение в виде:

т2 = ККыа (ЬхФ + к2Ь3V)-

-к3 (й3ф - кАЪ$) - к5 (63Ф - к6Ь^) ' где необходим подбор шести коэффициентов усиления. При этом к2 = к5 = 0, к4 = 1 дает управление (6). Параметры к5 = 0, к2 = к4 = 1 задают управление, аналогичное (5). В [46] рассмотрено гашение нутационных колебаний в случае эллиптической орбиты. В [47] добавлена интегральная часть в управление с обратной связью.

Спутники с двойным вращением и стабилизируемые вращением имеют высокий кинетический момент. Его изменение с помощью магнитной системы ориентации неизбежно медленно. Однако для аппарата с двойным вращением сохраняется возможность быстрого разворота вокруг оси установки маховика. Поэтому в этой области при построении управления применяются оптимизационные методы, в частности, #2-управление [48] и ^-управление [49; 50]. При этом уравнения движения линеаризуются и сводятся к уравнениям с периодическими коэффициентами, как и при построении линейно-квадратичного регулятора.

Скользящее управление было рассмотрено в [51] для аппарата на полярной орбите. В результате движение вне плоскости орбиты отделяется, и скользящая поверхность для соответствующих уравнений движения имеет вид:

5 = ^ф + к2\\> + £3ф + к4\у.

Коэффициенты усиления к меняются со временем из-за их зависимости от вектора геомагнитной индукции. Управление релейное:

m„

0,

-mn

s (t) > s -8 < S (t)< S S (t) < -8

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

2.3. Пассивная гравитационная стабилизация

Пассивная гравитационная стабилизация - один из надежных и простых подходов [13]. Однако малые аппараты в силу небольших размеров и массы позволяют создавать лишь слабый управляющий гравитационный момент. Установка штанги представляет проблему из-за обычно плотной компоновки аппарата и ее высокой стоимости. В результате магнитная система используется чтобы «усилить» гравитационный момент, особенно для вращения вокруг местной вертикали, и для обеспечения асимптотической устойчивости с помощью демпфирования угловой скорости (1) [52].

В [43] предложено управление (7), построенное на основе (4). Заменив вектор геомагнитной индукции и его производную на углы ориентации и их производные, можно добиться управления вращением вокруг местной вертикали. Управление (7) простое и может использовать необработанные показания магнитометра в переходном режиме для достижения требуемой ориентации с невысокой точностью [53]. После этого управление передается более точному алгоритму

Мге/ = -

В [54] проведен анализ устойчивости аппарата под действием управления (8) и гравитационного момента. Показана устойчивость при небольшой угловой скорости. В противном случае применяется лишь демпфирующая часть в (8).

2.4. Другие подходы

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

В случае использования одного реактивного маховика управляющий момент можно представить в виде [56]:

f о о -b2 Y m

m =

-b 1 -b

V Ь2

0 0

Л

wheel

m

77

з у

ТЕХНОЛОГИИЕ

гг

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

78 Возможна реализация одноосной ориента-

ции аппарата без обязательного требования высокой скорости закрутки. Так, в [59] предложено управление:

т = К (е3 х Ь),

реализуемое двумя катушками, так что К = diag(k, к, 0). В [60] требуемое управление, которое затем реализуется в виде проекции, имеет вид:

Мге/ = кх5е3 + ¿2е1 х 8тах. Здесь 5 - нормализованная разница токосъема двух противоположных солнечных панелей, задача первого слагаемого - уменьшение этой разницы. Вектор 8тах задает нормаль к панели с максимальным токосъемом, задача второго слагаемого - максимизация токосъема за счет разворота орта ех связанной системы на Солнце. В [61; 62] предложено управление:

т = к • cosа(w х 8).

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

В [63] для ориентации на Солнце используется управление:

т = -2^щ + £2 , т2 = -2£ (1 + Хц) ю3 / X,

т3 = -( 2£,ю3 + £2 ) / X, где X = А / С; ц = (В - С) / А; А, В, С - моменты инерции; £ - параметр алгоритма.

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

ления. Так, в [64] п т = к1

зедложена модификация (7):

ге/

-к2(ь

ге/

Том з

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

Среди других экзотических подходов, используемых совместно с магнитной системой ориентации, можно выделить силу Лоренца [65], солнечное давление [66], кольца с проводящей жидкостью [67].

3. Трехосная магнитная система ориентации

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

3.1. Простое локальное управление

Управление с обратной связью

Наиболее прямолинейный подход - использование обычного управления с обратной связью (8). Этот подход был реализован в [68; 69], в результате дипольный момент, задающий проекцию требуемого момента на плоскость, перпендикулярную вектору геомагнитной индукции, имеет вид:

т = -кт В х ю - kq В х q = -кт В х ю - к'аВ х S,

где 8 = (а23 - а32, а3Х - ав, ап - а2х) задает эквивалентную кватерниону формулировку закона управления через элементы матрицы направляющих косинусов аИсходное управление строится на основе функции Ляпунова. Реальное управление (9), реализуемое магнитной системой, не обладает теми же свойствами глобальной асимптотической устойчивости. Необходим аккуратный подбор коэффициентов усиления и общее ограничение величины управляющего момента [70]. Поскольку управление поворачивает спутник лишь примерно в нужном направлении, происходит медленное, итерационное приближение его к требуемой ориентации. Большая величина управляющего момента на такте управления может привести к ошибке ориентации, большей, чем была до его отработки. При этом демпфирующая часть в (9) - управление (1), которое обеспечивает асимптотическое стремление скорости к нулю. Поэтому ограничение на величину этой части управления мягче.

Выбор коэффициентов усиления в (9) - центральный момент в исследовании этого управления. Простая схема, основанная на аналитических

результатах и теории Флоке, предложена в [71] для аппарата, стабилизируемого в инерциальной системе координат, и в [72] для стабилизации в орбитальных осях. Управление особенно чувствительно к точности знания тензора инерции [73], и даже небольшая ошибка может привести к его неработоспособности, тогда как другие источники возмущений имеют в основном количественный характер. Однако при наличии ограничений на величины моментов инерции подбор коэффициентов усиления возможен [74]. В [75] подбор коэффициентов проводится на основе минимизации времени стабилизации. Принято во внимание ограничение на дипольный момент катушек, так как оптимизационный подход неизбежно приводит к завышенному значению дипольного момента. Более того, так как скорость стабилизации зависит от начальных условий, строится минимаксная задача с целью поиска коэффициентов, работоспособных в худших ситуациях. В [26] вместо скалярных коэффициентов используются матрицы.

В [42] приведены результаты безуспешной работы алгоритма (9) на борту спутника Gurwin. На аппарате TANGO миссии PRISMA управление (9) с подбором коэффициентов усиления с помощью линейно-квадратичного регулятора успешно решило поставленные задачи с точностью ориентации около 10-20° [76; 77].

Модификация (9), не требующая знания угловой скорости, предложена в [68]:

5 = а([<?0, q]-sX5),

u = -s2 (+ k2аШ(q0,q)([q0, q]-sXb)).

Здесь угловая скорость заменена новой переменной 5, вводятся два дополнительных параметра управления а и X, W - матрица кинематических уравнений.

В [78] рассмотрена реализация управления (9) на бортовом компьютере. Управление кусочно-постоянное. При этом вместо простой дискретизации исходного управления проводится дискретизация на заранее заданном интервале управления для нескольких шагов управления. Предложенная схема выбора параметров управления требует коротких промежутков измерения и расчета управления. Аналогичные результаты для интервалов измерений, сравнимых с интервалами управления, приведены в [79]. Задержка в выдаче управляющего воздействия рассмотрена в [80].

Основные результаты по подбору коэффициентов усиления верны для аппаратов с соизмеримыми моментами инерции. Управление (9), реализуемое в виде проекции ляпуновского управления, имеет низкое быстродействие для аппарата с моментами инерции, отличающимися в несколько раз [81]. Рассмотрим проекцию на плоскости, перпендикулярную вектору индукции, как мини-

мизацию евклидовои нормы разницы между требуемым и доступным управлением,

Чтобы учесть значительную разницу в моментах инерции, в [81] предложена норма с весовой матрицей:

J = (Mref - Mf Q (Mref - M)-

min

ге/ ) ^ V, г4

и при оптимизации принимается во внимание ограничение на направление управляющего момента М-Ь = 0. Управляющий момент имеет вид:

м=мге/ - О"1 ь (ьто-1 ь)-1 ьтмт.

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

В [82; 83] для расчета реализуемого момента используется типичный подход - небольшая добавка к обращаемой матрице векторного произведения, в результате ее ранг повышается до трех:

т = (в^ В + АЕ)"1 В^ Мге/. Здесь Вх - кососимметрическая матрица векторного произведения.

Управление (8) можно обобщить на любую функцию ошибки ориентации и ее производную вместо использования кватерниона и угловой скорости. В линейном приближении можно использовать углы ориентации и их производные. В [84] предложено управление:

M

ref

е = е — Ь

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

При построении управления зачастую предполагается, что спутник сферически симметричный. За счет этого из уравнений движения исключается слагаемое ю х ^ и динамические уравнения принимают вид х = и, где и - ограниченное по направлению управление. В результате появляется возможность подключить дополнительные инструменты теории управления. Такой подход был использован в [85; 86] для построения обратной связи без данных об угловой скорости, в [87; 88] для построения интуитивно понятного управления при стабилизации по местной вертикали и по вектору геомагнитной индукции и в [89] для построения управления на основе двух простых поворотов относительно вектора индукции.

Скользящее управление

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

79

ТЕХНОЛОГИИЕ

8о

гг

требуемой ориентации для сферически симметричного аппарата. При построении скользящей поверхности используется только угловая скорость, что применимо только к симметричному аппарату. В [91] используется классическая поверхность:

8 = Jw + Кц.

Управляющий момент должен быть параллелен вектору поверхности, и магнитная система реализует его проекцию. Коэффициент К - положительно определенная матрица. Однако в приведенной работе, как и в большинстве аналогичных работ, используется фактически скалярный коэффициент в виде К = кЕ. В [92] используется скользящая поверхность без «весовой» матрицы для компонент угловой скорости. В [93; 94] предложена поверхность с добавкой интегрального слагаемого в (10):

s = ю

К^ + К21 qdт + J1

ьь'

-и

dт.

Первый дополнительный член представляет собой среднюю ошибку ориентации, второй показывает близость требуемого управляющего момента к плоскости, перпендикулярной вектору геомагнитной индукции. Реализуется идея построения такой скользящей поверхности (и, соответственно, траектории углового движения аппарата), которая может быть реализована магнитной системой ориентации. Более того, предпринята попытка неявного учета возмущений, действующих на аппарат в требуемой ориентации. Если построение пути для переходного процесса возможно, средств прямого парирования возмущающего момента, направленного вдоль вектора индукции, нет. Однако наличие интегрального члена, зависящего от точности ориентации, позволяет опосредованно учесть возмущающие факторы. В [94] проведено моделирование работы такого управления и показано успешное парирование различных возмущений. В [95] строится скользящая поверхность, полностью доступная к реализации магнитной системой. Для этого используется зависящая от времени матрица К и итерационный подход к ее построению.

Нелинейная скользящая поверхность была предложена в [96] в виде:

s = Л« + К \д0q|а sign(д^)-

К 2

Том з

3.2. Оптимизационные методы

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

Один из наиболее популярных и простых методов оптимального управления - линейно-квадратичный регулятор (ЛКР). Ориентация спутника зачастую рассматривается в малой окрестности требуемого движения, так что линеаризация уравнений вполне естественна. Более того, так как вектор геомагнитной индукции меняется практически периодически, возникает линейная система с периодическими коэффициентами с хорошо разработанным аппаратом теории управления. ЛКР был использован в [98] для линеаризованных уравнений:

х = Ах + В(?)т + В/№

с периодической (в модели прямого диполя) матрицей В и постоянной матрицей возмущений В^ Функционал имеет вид:

J =

1 /(

+ mRm ) dt +

с диагональными матрицами и параметром ае(0, 1). Магнитное управление, однако, строится с помощью проекции. Тем не менее, это управление успешно справляется с возмущениями, в том числе с неточностью знания тензора инерции. В [97] для расчета наилучшего управляющего момента, реализуемого магнитной системой, используется метод наименьших квадратов, дающий лучший результат, чем простая проекция.

с весовыми матрицами Q, Я, Р, задающими влияние ошибки ориентации в ходе движения на одном периоде изменения магнитного поля, величины управления и конечной ошибки ориентации. Несмотря на то, что ЛКР разрабатывается для линейных уравнений движения, моделирование показало его работоспособность при ошибке ориентации около 30°. В [99] внешние возмущения аппроксимируются периодическими функциями. В [100] предложен алгоритм решения периодического уравнения Риккати в случае использования магнитной системы ориентации. В [101] применение ЛКР расширено на нелинейную систему с помощью уравнения Риккати, зависящего от вектора состояния системы (SDRE).

В [102] строятся ошибки ориентации для двух заданных направлений и используется ошибка в угловой скорости, предложенная в [103] и аналогичная алгоритму Bdot. Управление строится на основе усредненной за несколько витков ошибки. Таким образом, используется вращение вектора геомагнитной индукции, обеспечивающее управляемость «в среднем». При оптимизации ошибки используются генетические алгоритмы.

Н2-управление было предложено в [104], в [105] Н„-управление применяется для подбора коэффициентов в управлении (9) для линеаризованных уравнений движения. Мера эффективности управления задается выражением:

z = ( Е, 03х3, о К) х.

Здесь х - вектор состояния, включающий три компоненты угловой скорости и векторную часть кватерниона. Управление имеет вид обратной связи и = Кх, коэффициент о играет роль штрафа за величину управляющего момента. Указанная мера может использоваться для построения линейно-квадратичного функционала 3 = гтг. Однако для парирования возмущений используется норма и функционал имеет вид:

J ( K )= sup

T

J = J dt + k q2 + k2ю2

x ( k + N ) = AN x ( k ) +

A N-1 • A N -1B • AN -2 B2...

...[u(k),u(k +1),u(k + 2),...,u(k + N -1)].

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

N

J = Х x

j=i

с весами Q и R для ошибки ориентации и величины управления. Отличие от ЛКР заключается в горизонте предсказания. Движение сравнивается с требуемым на заданном интервале, вместо всего времени движения или периода коэффициентов линейных уравнений. Более того, на каждой итерации происходит пересчет управления. Такой подход предоставляет большую свободу в построении управления. Например, можно потребовать постоянного управления на определенном числе шагов, с последующим пересчетом на новое значение, снова постоянное на том же числе шагов, что отражает работу реальной системы управления. Можно упростить расчет правых частей уравнений движения, считая вектор геомагнитной индукции и возмущающие воздействия постоянными в течение нескольких итераций. Так, в [12; 49] используется метод наименьших квадратов для аппроксимации геомагнитного поля постоянной и периодическими частями с орбитальной и удвоенной орбитальной частотами. При этом, выбирая небольшой горизонт предсказания, можно ограничиться постоянными величинами в модели поля и возмущениях [109]. В [110] строится задача квадратичного программирования для обобщенного вектора параметров, включающего вектор состояния и управляющий момент.

Нелинейный MPC с непрерывным функционалом используется в [111] для оптимизация демпфирования угловой скорости. Функционал:

t +T

J ( t) = j f (x (s ), u (s )) ds + Ekm [x (t + T )]

t

рассматривается на горизонте предсказания T. Квадратичная функция f используется совместно с кинетической энергией после окончания процесса демпфирования. Показано, что последнее ограничение приводит к значительному возрастанию числа итераций для поиска управления. Введение условия остановки:

J ( x ( t ), u ( t ), t )<

< J (x (t - 5t), u (t - 5t), t - 5t)- yx2 (t) сокращает время расчета. Скорость изменения функционала на интервале St и величина скорости в начале итерации в момент t показывают, при каком изменении функционала можно считать, что процесс сошелся. Изменение параметра у позволяет задать требуемую «скорость» изменения функционала. Предложенный подход показал выигрыш в быстродействии по сравнению с неограниченным MPC для начальных скоростей бо-

За счет такого функционала появляется возможность оптимизации эффективности парирования возмущающих моментов Т^. Нш-управление также используется в [106] для дискретной системы, включающей модели измерений и возмущений.

В работе [107] используется классический функционал - время совершения маневра с дополнительным штрафом на конечную точность ориентации и скорость:

Время маневра T подлежит определению. Поэтому вводятся дополнительные переменные, отражающие масштабирование уравнений движения со временем, и уравнения рассматриваются на фиксированном интервале. Построенная угловая траектория отслеживается с помощью управления с предсказанием динамики (model predictive control, MPC).

3.3. Model predictive control

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

x ( k +1 ) = Ax ( k ) + Bu ( k ).

На k + N шаге:

8i

z

2

2

ТЕХНОЛОГИИ t

82

m

лее 5 °/с. В [112] предложен функционал, использующий свойства управления (1):

J =

J [kg2 + mTRmys + Ekin [x(t + T)],

где g = 2т0 / п • агСаи [К(Ь х га)] . Если обычный квадратичный функционал использует величину угловой скорости напрямую, здесь используется отклонение от траектории, задаваемой управлением (1). Эта траектория рассматривается как первое приближение.

3.4. Адаптивные методы

Адаптивные методы, как и классические методы оптимизации, позволяют учитывать ограничения. Основная проблема при их использовании - высокая вычислительная сложность. Вопрос реализации такого управления на борту реального аппарата, однако, обычно не рассматривается. При этом зачастую используется комбинация классического управления с обратной связью и адаптивных методов. Так, в [113] динамическая нейронная сеть используется для подбора параметров управления в (8). Результаты моделирования показывают эффективность подобранных коэффициентов, в особенности при парировании возмущений неизвестной природы и величины.

Заключение

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

Том з

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

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

Список литературы

[1] Fischell R. E. Magnetic damping of the angular motions of Earth satellites // Journal of the American Rocket Society, 1961, vol. 31, no. 9, pp. 1210-1217.

[2] Grasshoff L. H. A method for controlling the attitude of a spin-stabilized satellite // Journal of the American Rocket Society, 1961, vol. 31, no. 5, pp. 646-649.

[3] Белецкий В. В., Зонов Ю. В. Вращение и ориентация третьего советского спутника // Сборник "Искусственные спутники Земли". АН СССР, 1961. № 7. С. 32-55.

[4] Bhat S. P. Controllability of nonlinear time-varying systems: applications to spacecraft attitude control using magnetic actuation // IEEE Transactions on Automatic Control, 2005, vol. 50, no. 11, pp. 1725-1735.

[5] Bhat S. P., Dham A. S. Controllability of spacecraft attitude under magnetic actuation // 42nd IEEE International Conference on Decision and Control. Maui, HI, USA, 2003, vol. 3, pp. 2383-2388.

[6] Navabi M., Barati M. Mathematical modeling and simulation of the earth's magnetic field: A comparative study of the models on the spacecraft attitude control application // Applied Mathematical Modelling, 2017, vol. 46, pp. 365-381.

[7] Ovchinnikov M. Yu., Penkov V. I., Roldugin D. S., Pichuzhkina A. V. Geomagnetic field models for satellite angular motion studies // Acta Astronautica, 2018, vol. 144, pp. 171-180.

[8] Тихонов А. А., Петров К. Г. Мультипольные модели магнитного поля Земли // Космические исследования. 2002. Т. 40. № 3. С. 219-229.

[9] Yang Y. Controllability of spacecraft using only magnetic torques // IEEE Transactions on Aerospace and Electronic Systems, 2016, vol. 52, no. 2, pp. 954-961.

[10] Smirnov G., Ovchinnikov M., Miranda F. On the magnetic attitude control for spacecraft via the epsilon-strategies method // Acta Astronautica, 2008, vol. 63, no. 5-6, pp. 690-694.

[11] Ovchinnikov M. Yu., Roldugin D. S. A survey on active magnetic attitude control algorithms for small satellites // Progress in Aerospace Sciences, 2019. doi: 10.1016/j.paerosci.2019.05.006.

[12] Silani E., Lovera M. Magnetic spacecraft attitude control: a survey and some new results // Control Engineering Practice, 2005, vol. 13, no. 3, pp. 357-371.

[13] Сарычев В. А. Вопросы ориентации искусственных спутников / Итоги науки и техники. Серия "Исследование космического пространства". М. : ВИНИТИ, 1978. T. 11. 221 с.

[14] Shrivastava S. K., Modi V. J. Satellite attitude dynamics and control in the presence of environmental torques - a brief survey // Journal of Guidance, Control, and Dynamics, 1983, vol. 6, no. 6, pp. 461-471.

[15] Сарычев В. А., Овчинников М. Ю. Магнитные системы ориентации искусственных спутников Земли / Итоги науки и техники. Серия: "Исследование космического пространства". М. : ВИНИТИ, 1985. 106 с.

[16] Коваленко А. П. Магнитные системы управления космическими летательными аппаратами. М. : Машиностроение, 1975. 248 с.

[17] Овчинников М. Ю., Пеньков В. И., Ролдугин Д. С., Иванов Д. С. Магнитные системы ориентации малых спутников. М. : ИПМ им. М. В. Келдыша, 2016. 366 с.

[18] Stickler A. C., Alfriend K. T. Elementary magnetic attitude control system // Journal of Spacecraft and Rockets, 83 1976, vol. 13, no. 5, рр. 282-287.

[19] Ovchinnikov M. Yu., Roldugin D. S., Tkachev S. S., Penkov V. I. B-dot algorithm steady-state motion performance // Acta Astronautica, 2018, vol. 146, рр. 66-72.

[20] Сарычев В. А., Сазонов В. В. Оценка влияния диссипативного магнитного момента от вихревых токов на быстрое вращение спутника // Космические исследования. 1982. Т. 20, № 2. С. 297-300.

[21] Lovera M. Magnetic satellite detumbling: The B-dot algorithm revisited / Proceedings of the American Control Conference. Chicago, 2015. pp. 1867-1872.

[22] Драновский В. И., Яншин А. М. Влияние диссипативных моментов от вихревых токов на ориентацию спутника, стабилизированного вращением // Космические исследования. 1975. T. 13, № 4. С. 487-493.

[23] Сазонов В. В., Сарычев В. А. Влияние диссипативного магнитного момента на вращение спутника относительно центра масс // Изв. АН СССР, Мех. тв. тела. 1983. T. 2. С. 3-12.

[24] Овчинников М. Ю., Пеньков В. И., Ролдугин Д. С., Карпенко С. О. Исследование быстродействия алгоритма активного магнитного демпфирования // Космические исследования. 2012. Т. 50, № 2. С. 176-183.

[25] Avanzini G., Giulietti F. Magnetic detumbling of a rigid spacecraft // Journal of Guidance, Control, and Dynamics, 2012, vol. 35, no. 4, pp. 1326-1334.

[26] Wisniewski R., Blanke M. Fully magnetic attitude control for spacecraft subject to gravity gradient // Automatica, 1999, vol. 35, no. 7, pp. 1201-1214.

[27] Семкин Н. Д., Любимов В. В., Малышев В. И. Моделирование законов функционирования магнитных исполнительных органов при ориентации микроспутника по местной вертикали // Физика волновых процессов и радиотехнические системы. 2012. T. 15, № 1. С. 103-108.

[28] Любимов В. В., Подклетнова С. В. Оптимальные законы управления для уменьшения кинетического момента и демпфирования угловой скорости наноспутников и микроспутников с магнитными катушками на борту // Вестник самарского государственного аэрокосмического университета им. акад. С. П. Королёва. 2016. T. 15, № 2. С. 57-67.

[29] Shigehara M. Geomagnetic attitude control of an axisymmetric spinning satellite // Journal of Spacecraft and Rockets, 1972, vol. 9, no. 6, pp. 391-398.

[30] Vega K., Auslander D., Pankow D. Design and modeling of an active attitude control system for CubeSat class satellites / AIAA Modeling and Simulation Technologies Conference. Reston, Virigina, 2009. AIAA 2009-5812.

[31] Овчинников М. Ю., Ролдугин Д. С., Пеньков В. И. Исследование связки трех алгоритмов магнитного управления угловой скоростью и ориентацией спутника, стабилизируемого вращением // Космические исследования. 2012. T. 50, № 4. С. 326-334.

[32] Ovchinnikov M. Yu., Roldugin D. S., Penkov V. I. Asymptotic study of a complete magnetic attitude control cycle providing a single-axis orientation // Acta Astronautica, 2012, vol. 77, pp. 48-60.

[33] Zavoli A., Giulietti F., Avanzini G., de Matteis G. Spacecraft dynamics under the action of Y-dot magnetic control law // Acta Astronautica, 2016, vol. 122, pp. 146-158.

[34] Thomson W. T. Spin stabilization of attitude against gravity torque // The Journal of the Astronautical Sciences, 1962, vol. 9, no. 1. AAS 9-1-31-3.

[35] Avanzini G., de Angelis E. L., Giulietti F. Spin-axis pointing of a magnetically actuated spacecraft // Acta Astronautica, 2014, vol. 94, no. 1, pp. 493-501.

[36] de Ruiter A. A fault-tolerant magnetic spin stabilizing controller for the JC2Sat-FF mission // Acta Astronautica, 2011, vol. 68, no. 1-2, pp. 160-171.

[37] Alfriend K. T. Magnetic attitude control system for dual-spin satellites // AIAA Journal, 1975, vol. 13, no. 6, pp. 817-822.

[38] Wheeler P. C. Spinning spacecraft attitude control via the environmental magnetic field // Journal of Spacecraft and Rockets, 1967, vol. 4, no. 12, pp. 1631-1637.

[39] Ergin E. I., Wheeler P. C. Magnetic attitude control of a spinning satellite. // Journal of Spacecraft and Rockets, 1965, vol. 2, no. 6, pp. 846-850.

ТЕХНОЛОГИИ t

84

[40

[41 [42 [43 [44 [45 [46 [47 [48

[49 [50 [51 [52 [53 [54 [55 [56 [57 [58 [59 [60 [61 [62 [63 [64 [65 [66 [67

^ Том 3

Овчинников М. Ю., Пеньков В. И., Ролдугин Д. С., Варатарао Р., Рябиков В. С. Движение спутника, оснащенного тангажным маховиком и магнитными катушками, в гравитационном поле // Космические исследования. 2017. T. 55, № 3. С. 218-225.

Goel P. S., Rajaram S. Magnetic attitude control of a momentum-biased satellite in near-equatorial orbit // Journal of Guidance, Control, and Dynamics, 1979, vol. 2, no. 4, pp. 334-338.

Guelman M. M., Waller R., Shiryaev A., Psiaki M. Design and testing of magnetic controllers for Satellite stabilization // Acta Astronautica, 2005, vol. 56, no. 1-2, pp. 231-239.

Martel F., Pal P. K., Psiaki M. Active magnetic control system for gravity gradient stabilized spacecraft / Proceedings of the 2nd Annual AIAA/USU Conference on Small Satellites. Logan, USA, 1988. 19 p.

Pittelkau M. E. Optimal periodic control for spacecraft pointing and attitude determination // Journal of Guidance, Control, and Dynamics, 1993, vol. 16, no. 6, pp. 1078-1084.

Hablani H. B. Comparative stability analysis and performance of magnetic controllers for bias momentum satellites // Journal of Guidance, Control, and Dynamics, 1995, vol. 18, no. 6, pp. 1313-1320.

Hablani H. B. Magnetic precession and product-of-inertia nutation damping of bias momentum satellites // Journal of Guidance, Control, and Dynamics, 1995, vol. 18, no. 6, pp. 1321-1328.

Tsuchiya K., Inoue M. New control schemes for a magnetic attitude control system // IFAC Proceedings Volumes, 1983, vol. 16, no. 11, pp. 221-225.

Trégouët J.-F., Arzelier D., Peaucelle D., Ebihara Yo., Pittet C., Falcoz A. Periodic H2 synthesis for spacecraft attitude control with magnetorquers and reaction wheels / IEEE Conference on Decision and Control and European Control Conference. Orlando, USA, 2011. pp. 6876-6881.

Lovera M. Periodic H-inf attitude control for satellites with magnetic actuators // IFAC Proceedings Volumes, 2000, vol. 33, no. 14, pp. 631-636.

Calloni A., Corti A., Zanchettin A. M., Lovera M. Robust attitude control of spacecraft with magnetic actuators / 2012 American Control Conference, 2012. pp. 750-755.

Wang P., Shtessel Y. Satellite attitude control via magnetorquers using switching control laws // IFAC Proceedings Volumes, 1999, vol. 32, no. 2, pp. 8021-8026.

Игнатов А. И., Сазонов В. В. Стабилизация режима гравитационной ориентации искусственного спутника Земли электромагнитной системой управления / Препринты ИПМ им. М. В. Келдыша. 2016. № 28. 32 c. Grassi M. Attitude determination and control for a small remote sensing satellite // Acta Astronautica, 1997, vol. 40, no. 9, pp. 675-681.

Lovera M., Astolfi A. Global magnetic attitude control of spacecraft in the presence of gravity gradient // IEEE Transactions on Aerospace and Electronic Systems, 2006, vol. 42, no. 3. pp. 796-805.

Pulecchi T., Lovera M. Attitude control of spacecraft with partially magnetic actuation // IFAC Proceedings Volumes, 2007, vol. 40, no. 7, pp. 609-614.

Desouky M. A. A., Prabhu K., Abdelkhalik O. On spacecraft magnetic attitude control / Space Flight Mechanics Meeting. Reston, Virginia, 2018. AIAA 2018-0205.

Forbes J. R., Damaren C. J. Geometric approach to spacecraft attitude control using magnetic and mechanical actuation // Journal of Guidance, Control, and Dynamics, 2010, vol. 33, no. 2, pp. 590-595. de Ruiter A. Magnetic control of dual-spin and bias-momentum spacecraft // Journal of Guidance, Control, and Dynamics, 2012, vol. 35, no. 4, pp. 1158-1168.

Jan Y. W., Tsai J.-R. Active control for initial attitude acquisition using magnetic torquers // Acta Astronautica, 2005, vol. 57, no. 9, pp. 754-759.

Kim J., Worrall K. Sun tracking controller for UKube-1 using magnetic torquer only // IFAC Proceedings Volumes, 2013, vol. 46, no. 19, pp. 541-546.

Ovchinnikov M. Yu., Roldugin D. S., Tkachev S. S., Karpenko S. O. New one-axis one-sensor magnetic attitude control theoretical and in-flight performance // Acta Astronautica, 2014, vol. 105, no. 1, pp. 12-16. Karpenko S. O., Ovchinnikov M. Yu., Roldugin D. S., Tkachev S. S. One-axis attitude of arbitrary satellite using magnetorquers only // Cosmic Research, 2013, vol. 51, no. 6, pp. 478-484.

Игнатов А. И., Сазонов В. В. Стабилизация режима солнечной ориентации искусственного спутника Земли электромагнитной системой управления // Космические исследования. 2018. T. 56, № 5. С. 375-383. Psiaki M. L. Nanosatellite attitude stabilization using passive aerodynamics and active magnetic torquing // Journal of Guidance, Control, and Dynamics, 2004, vol. 27, no. 3, pp. 347-355.

Антипов К. А., Тихонов А. А. Параметрическое управление в задаче о стабилизации космического аппарата в магнитном поле земли // Автоматика и телемеханика. 2007. № 8. С. 44-56.

Modi V. J., Pande K. C. Magnetic-solar hybrid attitude control of satellites in near-equatorial orbits // Journal of Spacecraft and Rockets, 1974, vol. 11, no. 12, pp. 845-851.

Nobari N. A., Misra A. K. A hybrid attitude controller consisting of electromagnetic torque rods and an active fluid ring // Acta Astronautica, 2014, vol. 94, no. 1, pp. 470-479.

[68] Lovera M., Astolfi A. Spacecraft attitude control using magnetic actuators // Automatica, 2004, vol. 40, no. 8, pp. 1405-1414.

[69] Lovera M., Astolfi A. Global magnetic attitude control of inertially pointing spacecraft // Journal of Guidance, Control, and Dynamics, 2005, vol. 28, no. 5, pp. 1065-1072.

[70] Rossa F. D., Bergamasco M., Lovera M. Bifurcation analysis of the attitude dynamics for a magnetically controlled spacecraft / 51st IEEE Conference on Decision and Control. Maui, HI, USA, 2012. pp. 1154-1159.

[71] Ovchinnikov M. Yu., Roldugin D. S., Penkov V. I. Three-axis active magnetic attitude control asymptotical study // Acta Astronautica, 2015, vol. 110, pp. 279-286.

[72] Ovchinnikov M. Yu., Roldugin D. S., Ivanov D. S., Penkov V. I. Choosing control parameters for three axis magnetic stabilization in orbital frame // Acta Astronautica, 2015, vol. 116, pp. 74-77.

[73] Ivanov D. S., Ovchinnikov M. Yu., Penkov V. I., Roldugin D. S., Doronin D. M., Ovchinnikov A. V. Advanced numerical study of the three-axis magnetic attitude control and determination with uncertainties // Acta Astronautica,

2017, vol. 132, pp. 103-110. 85

[74] Celani F. Robust three-axis attitude stabilization for inertial pointing spacecraft using magnetorquers // Acta Astronautica, 2015, vol. 107, pp. 87-96.

[75] Bruni R., Celani F. A robust optimization approach for magnetic spacecraft attitude stabilization // Journal of Optimization Theory and Applications, 2017, vol. 173, no. 3, pp. 994-1012.

[76] Chasset C., Noteborn R., Bodin P., Larsson R., Jakobsson B. 3-Axis magnetic control: flight results of the TANGO satellite in the PRISMA mission // CEAS Space Journal, 2013, vol. 5, no. 1-2, pp. 1-17.

[77] Bodin P., Larsson R., Nilsson F., Chasset C., Noteborn R., Nylund M. PRISMA: An In-Orbit Test Bed for Guidance, Navigation, and Control Experiments // Journal of Spacecraft and Rockets, 2009, vol. 46, no. 3, pp. 615-623.

[78] Celani F. Spacecraft attitude stabilization with piecewise-constant magnetic dipole moment // Journal of Guidance, Control, and Dynamics, 2016, vol. 39, no. 5, pp. 1140-1146.

[79] Celani F. Spacecraft attitude stabilization using magnetorquers with separation between measurement and actuation // Journal of Guidance, Control, and Dynamics, 2016, vol. 39, no. 9, pp. 2184-2191.

[80] Xu C., Luo W., Yang X. Three-axis magnetic attitude control by delayed output feedback / Chinese Automation Congress. Jinan, China, 2017. pp. 3544-3549.

[81] Wood M., Chen W. Attitude control of magnetically actuated satellites with an uneven inertia distribution // Aerospace Science and Technology, 2013, vol. 25, no. 1, pp. 29-39.

[82] Sugimura N., Kuwahara T., Yoshida K. Attitude determination and control system for nadir pointing using magnetorquer and magnetometer / IEEE Aerospace Conference. Big Sky, MT, USA, 2016. 12 p.

[83] Sakai Sh, Fukushima Yo., Saito H., Kaneda R. Studies on magnetic attitude control system for the REIMEI microsatellite / AIAA Guidance, Navigation, and Control Conference and Exhibit. Reston, Virigina, 2006. AIAA 2006-6450.

[84] Gulmammadov F., Kahraman O., Yavuzyilmaz C., Tufekci C. S., Subasi Y. Magnetorquers only attitude maintaining using dynamic attitude simulator environment / AIAA Modeling and Simulation Technologies Conference. Toronto, 2010. AIAA 2010-8106.

[85] Reyhanoglu M., Ton C., Drakunov S. Attitude Stabilization of a nadir-pointing small satellite using only magnetic actuators // IFAC Proceedings Volumes, 2009, vol. 42, no. 19, pp. 292-297.

[86] Reyhanoglu M., Drakunov S. Attitude stabilization of small satellites using only magnetic actuation / 34th Annual Conference of IEEE Industrial Electronics, 2008. pp. 103-107.

[87] Bushenkov V. A., Ovchinnikov M. Yu., Smirnov G. V. Attitude stabilization of a satellite by magnetic coils // Acta Astronautica, 2002, vol. 50, no. 12, pp. 721-728.

[88] Smirnov G. V Attitude determination and stabilization of a spherically symmetric rigid body in a magnetic field // International Journal of Control, 2001, vol. 74, no. 4, pp. 341-347.

[89] Inamori T., Otsuki K., Sugawara Yo., Saisutjarit Ph., Nakasuka Sh. Three-axis attitude control by two-step rotations using only magnetic torquers in a low Earth orbit near the magnetic equator // Acta Astronautica, 2016, vol. 128, pp. 696-706.

[90] Wang P., Shtessel Y. Satellite attitude control using only magnetorquers / Proceedings of the Thirtieth Southeastern Symposium on System Theory. Morgantown, West Virginia, 1998. pp. 500-504.

[91] Wisniewski R. Sliding mode attitude control for magnetic actuated satellite // IFAC Proceedings Volumes, 1998, vol. 31, no. 21, pp. 179-184.

[92] Sofyali A., Jafarov E. M. Purely magnetic spacecraft attitude control by using classical and modified sliding mode algorithms / 12th International Workshop on Variable Structure Systems. Mumbai, 2012. pp. 117-123.

[93] Sofyali A., Jafarov E. M. Integral sliding mode control of small satellite attitude motion by purely magnetic actuation // IFAC Proceedings Volumes, 2014, vol. 47, no. 3, pp. 7947-7953.

[94] Sofyali A., Jafarov E. M., Wisniewski R. Robust and global attitude stabilization of magnetically actuated spacecraft through sliding mode // Aerospace Science and Technology, 2018, vol. 76, pp. 91-104.

[95] Ovchinnikov M. Yu., Roldugin D. S., Penkov V I., Tkachev S. S., Mashtakov Y. V. Fully magnetic sliding mode control for acquiring three-axis attitude // Acta Astronautica. 2016. vol. 121. pp. 59-62.

TEXHOJIOMIt

86

m

TOM 3

96] Janardhanan S., Nabi M., Tiwari P. M. Attitude control of magnetic actuated spacecraft using super-twisting algorithm with nonlinear sliding surface / 12th International Workshop on Variable Structure Systems, 2012. pp. 46-51.

97] Schlanbusch R., Kristiansen R., Nicklasson P. J. Spacecraft magnetic control using dichotomous coordinate descent algorithm with box constraints // Modeling, Identification and Control, 2010, vol. 31, no. 4, pp. 123-131.

98] Psiaki M. L. Magnetic torquer attitude control via asymptotic periodic linear quadratic regulation // Journal of Guidance, Control, and Dynamics, 2001, vol. 24, no. 2, pp. 386-394.

99] Lovera M., de Marchi E., Bittanti S. Periodic attitude control techniques for small satellites with magnetic actuators // IEEE Transactions on Control Systems Technology, 2002, vol. 10, no. 1, pp. 90-95.

100] Yang Y. An efficient algorithm for periodic Riccati equation with periodically time-varying input matrix // Automatica, 2017, vol. 78, pp. 103-109.

101] Rodriguez-Vazquez A. L., Martin-Prats M. A., Bernelli-Zazzera F. Full magnetic satellite attitude control using ASRE method / 1st IAA Conference on Dynamics and Control of Space Systems. Porto, Portugal, 2012. AAS 12-347.

102] Colagrossi A., Lavagna M. Fully magnetic attitude control subsystem for picosat platforms // Advances in Space Research, 2018, vol. 62, no. 12, pp. 3383-3397.

103] Avanzini G., de Angelis E. L., Giulietti F. Acquisition of a desired pure-spin condition for a magnetically actuated spacecraft // Journal of Guidance, Control, and Dynamics, 2013, vol. 36, no. 6, pp. 1816-1821.

104] Wisniewski R., Stoustrup J. Periodic H2 synthesis for spacecraft attitude determination and control with a vector magnetometer and magnetorquers // IFAC Proceedings Volumes, 2001, vol. 34, no. 12, pp. 119-124.

105] Zanchettin A. M., Lovera M. H-inf attitude control of magnetically actuated satellites // IFAC Proceedings Volumes, 2011, vol. 44, no. 1, pp. 8479-8484.

106] Kulkarni J., Campbell M. An approach to magnetic torque attitude control of satellites via H-inf control for LTV systems / 43rd IEEE Conference on Decision and Control, 2004. pp. 273-277.

107] Liang J., Fullmer R., Chen Y. Time-optimal magnetic attitude control for small spacecraft / 43rd IEEE Conference on Decision and Control, 2004. pp. 255-260.

108] Wood M., Chen W., Fertin D. Model predictive control of low earth orbiting spacecraft with magneto-torquers / 2006 IEEE Conference on Computer Aided Control System Design, 2006. pp. 2908-2913.

109] Wood M., Chen W. Regulation of magnetically actuated satellites using model predictive control with disturbance modelling / 2008 IEEE International Conference on Networking, Sensing and Control, 2008. pp. 692-697.

110] Krogstad T., Gravdahl J.T., Tondel P. Explicit model predictive control of a satellite with magnetic torquers / Proceedings of the 2005 IEEE Mediterrean Conference on Control and Automation Intelligent Control, 2005. pp. 491-496.

111] Ahmed S., Kerrigan E. C. Suboptimal predictive control for satellite detumbling // Journal of Guidance, Control, and Dynamics, 2014, vol. 37, no. 3, pp. 850-859.

112] Böhm Ch., Merk M., Fichter W., Allgöwer F. // Lecture Notes in Control and Information Sciences: vol. 384, Nonlinear Model Predictive Control / ed. Magni L., Raimondo D. M., Allgower F. 2009. pp. 511-520.

113] Das S., Sinha M., Misra A. K. Dynamic neural units for adaptive magnetic attitude control of spacecrafts // Journal of Guidance, Control, and Dynamics, 2012, vol. 35, no. 4, pp. 1280-1291.

RECENT ADVANCES IN THE ACTIVE MAGNETIC CONTROL OF SATELLITES

M. Yu. Ovchinnikov, D. S. Roldugin

Keldysh Institute of Applied Mathematics of RAS, Moscow, Russian Federation

The paper covers main recent results in the active magnetic attitude control of satellites. Three main implementation situations are outlined. Angular velocity damping opens the problem as the auxiliary control task. Next, implementation with other actuators and passive stabilization concepts is considered. Magnetic attitude control is restricted in the direction: control torque cannot be applied along the magnetic induction vector. Other actuators or environmental properties may enhance the control, providing control authority along the restricted axis. This comes at the cost of restricted attitude motion. Passive gravitational stabilization, spin stabilization and dual spin satellites present main cases. The satellite may acquire the local vertical and one axis inertial attitude that represent important cases. The most challenging and practically promising situation is the fully magnetic three axis attitude control. This reduces the hardware requirements for the attitude control system to the minimum. However, this also comes at the cost of a restriction on the control torque vector and low attitude accuracy and time-response. Feedback law with proper control gains tuning, sliding control and optimization techniques are considered for this problem.

Keywords: satellite attitude, angular motion, magnetic attitude control, magnetorquer.