Динамика быстро вращающихся малых спутников в геомагнитном поле тема диссертации и автореферата по ВАК РФ 01.02.01, кандидат физико-математических наук Ильин, Андрей Александрович

• Ильин, Андрей Александрович
• кандидат физико-математических наук
• 2006

Оглавление диссертации кандидат физико-математических наук Ильин, Андрей Александрович

Введение

Глава 1. Активная магнитная система ориентации спутника, стабилизируемого собственным вращением

1.1. Уравнения движения спутника

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

Управление ориентацией оси вращения спутника

Управление скоростью вращения спутника

Алгоритм гашения нутационного движения спутника 22 1.3 Алгоритм управления магнитными катушками, не требующий измерения угловой скорости спутника (второй алгоритм)

Управление ориентацией оси вращения спутника

Управление скоростью вращения спутника

Алгоритм гашения нутационного движения спутника 28 Определение скорости вращения спутника на основании измерений магнитометров

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

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

2.2. Осредненная модель геомагнитного поля и используемые системы координат

2.3. Уравнения возмущенного движения

2.4. Исследование алгоритмов управления ориентацией 46 Уравнения движения при управлении ориентацией спутника 47 Осредненные уравнения по переменной v|/ 50 Случай вращения спутника вокруг оси симметрии

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

Глава 3. Моделирование работы активной магнитной системы ориеитации на спутнике ТНС

3.1. О спутнике ТНС

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

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

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

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

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

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

Глава 4. Динамика спутника ТНС-0 относительно центра масс

4.1.0 цели запуска спутника ТНС

4.2. Пассивная магнитная система ориентации спутника ТНС

4.3. Задача определение параметров движения наноспутника ТНС-0 по измерениям солнечных датчиков

4.4. Обработка кратковременных измерений при быстром вращении спутника

4.5. Движение быстро вращающегося спутника с постоянным магнитом

4.6. Асимптотическая устойчивость вращающегося спутника с постоянным магнитом и гистерезисными стержнями

4.7. Определение параметров движения спутника по измерениям, сделанным в течение одного витка

Введение диссертации (часть автореферата) на тему «Динамика быстро вращающихся малых спутников в геомагнитном поле»

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

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

Первой рассматриваемой в настоящей работе задачей являлась задача обеспечения ориентации осесимметричного спутника, стабилизируемого собственным вращением, при помощи активной магнитной системы. Одним из наиболее простых и часто используемых способов обеспечения заданного углового движения спутника заключается в том, что аппарат, закрученный вокруг оси с наибольшим моментом инерции, сохраняет достаточно долго неизменную ориентацию оси закрутки в инерциальном пространстве [33]. Системы ориентации, построенные по такому принципу, рассматриваются например в работах [1,42]. Для ориентации оси вращения и регулирования скорости собственного вращения можно использовать магнитную систему, использующую токовые катушки для создания магнитного дипольного момента спутника. Магнитная система обладает тем преимуществом, что она не требует расхода рабочего тела, что особенно важно для малогабаритных недорогих спутников, связанных множеством ограничений. Кроме того, отсутствие подвижных частей значительно повышает надежность, а, следовательно, и продолжительность активного функционирования такого аппарата.

Однако, система управления, построенная по такому принципу, обладает как минимум двумя недостатками: магнитное торможение и потеря осью вращения первоначальной ориентации из-за ряда возмущений. Спутник, вращаясь в магнитном поле Земли, теряет часть своей кинетической энергии за счет индуцирования вихревых токов в токопроводящих материалах, гистерезисных потерь и потерь на вихревые токи в магнитном материале. Например, скорость вращения спутника "Vanguard" (1958 02) уменьшилась в течение года с 2.7 до 0.55 рад/сек в основном за счет индуцирования вихревых токов в алюминиевом корпусе [1]. За счет потерь на гистерезис скорость вращения спутника (1960, 01) упала от 0.785 до менее чем 0.002 рад/сек за 24 дня.

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

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

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

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

Ранее в системах управления вращением спутников использовалось несколько способов стабилизации скорости вращения. На двух аппаратах метеорологического АЕ-В (The Atmosphere Explorer satellite) и спутника прямых измерений ДМЕ-А (The Direct Measurements Explorer satellite) использовалась магнитная система функционирующая, следующим образом 3

33]. Введем связанную со спутником систему координат Ох\х^х3. Ось Ох3 является осью вращения спутника. Векторные магнитометры измеряют напряженность магнитного поля по осям Ох\ и Ох2 спутника. Сигнал с магнитометра, установленного на оси Ох\, усиливается и подается на обмотку электромагнита оси 0x2. Подобным образом сигнал с магнитометра оси Ох2 усиливается, инвертируется и поступает на обмотку электромагнита оси Ох\.

Если проекция индукции магнитного поля Земли в плоскости 0х\х2 спутника - В12 то система, включающая в себя магнитометры и электромагниты, создаст магнитный диполь 1113, перпендикулярный В^. В зависимости от того, в какую сторону вращается спутник, момент 1113 будет увеличивать или уменьшать скорость вращения. При работе системы возникает гироскопический момент 1П3ХВ, вызывающий прецессию спутника. Если скорость вращения спутника велика, а произведение 1П3ХВ мало, то прецессионное движение спутника можно сделать медленным.

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

При полете спутников «АЕ-В» и «ДМЕ-А» требовалось, чтобы ось вращения была перпендикулярна плоскости орбиты. Так как эти аппараты двигались по полярным орбитам, напряженность магнитного поля Земли в требуемой плоскости достигает значительных величин, поэтому была возможность создавать заданную прецессию спутника.

Для функционирования большинства спутников требуется определенным образом ориентировать ось их вращения. Например, спутник "Tiros" ориентировался так, чтобы получить наилучшие условия наблюдения за метеорологических обстановкой [18]. Метеорологические спутники этой серии представляли из себя 18-угольный многоугольник около метра в диаметре и около полуметра высотой. Весил этот спутник 122 кг.

Решение о стабилизации аппарата посредством вращения было принято в процессе разработки первоначальной конструкции спутника "Tiros". Несмотря на ограничения, был выбран этот довольно простой метод. Телевизионные камеры были смонтированы на аппарате так, чтобы их оптические оси были направлены параллельно оси вращения, за которою была принята ось минимального момента инерции космического аппарата. Спутник был выведен на орбиту, вращающимся относительно оси, лежащей в плоскости орбиты, направленной по касательной в точке вывода с номинальной скоростью около 10 об/мин. При стабилизации аппарата выбранным способом его поведение подобно простому идеальному гироскопу, который при отсутствии возмущающих моментов сохраняет фиксированное направление в инерциальном пространстве.

Медленное снижение скорости вращения спутника, вызываемое влиянием магнитного поля Земли, возмещалось включением небольших твердотопливных двигателей, смонтированных на спутнике. Поэтому, в первом приближении, можно считать, что ось вращения спутника сможет поддерживать в инерциальном пространстве примерно постоянное направление, по крайней мере, в течение нескольких месяцев. Кинематика вывода и отделения спутника от носителя была выбрана такой, чтобы в течении первых нескольких недель ось вращения спутника при пролете над Северным полушарием была бы направлена всегда вниз. Такая ориентация позволяет получить информацию о состоянии облачного покрова и инфракрасном излучении Земли. Однако прецессия орбитальной плоскости со скоростью приблизительно 5 градусов в день уже через две недели так изменила положение плоскости орбиты, что ось вращения спутника перестала лежать в ее плоскости, а зона полезного оборота исчезла. Явление непрерывной прецессии плоскости орбиты приводит к периодическому восстановлению полезных условий, при которых снова можно получать полезные снимки. Таким образом, создавалась ситуация, при которой последовательно чередуются промежутки времени, когда, например, в течении 2-х недель имеется хорошая зона обзора затем зона плохого обзора и т.д.

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

Основные источники внешних моментов имеют гравитационную и магнитную природу [9]. Несимметричное твердое тело будет испытывать в гравитационном поле ощутимый момент, если главные моменты инерции тела не равны между собой. У спутника "Tiros I" главный момент инерции относительно оси собственного вращения больше, чем остальные главные моменты инерции. Оценка величины результирующего гравитационного момента показала, что, несмотря на очевидное влияние момента на спутник, этим моментом можно пренебречь. Выяснилось, что основное возмущающее воздействие является по своей природе магнитным.

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

Уравнения получены на базе теоремы об изменении кинетического момента. При выводе уравнений был принят ряд упрощающих допущений. Используемые уравнения имели следующую форму [9]: = ~£ (s • n)(s х n) + //(b х s),

3<Ул I-J М V0 т > М = ---j,

COs I I(DS где s - единичный вектор скорости вращения спутника вокруг центра масс, coq - угловая скорость радиус-вектора, характеризующего положение спутника на орбите, (os - скорость вращения спутника вокруг его центра масс, / - момент инерции спутника относительно оси собственного вращения, J - поперечный момент инерции спутника, п - единичный вектор нормали к орбите, М -составляющая момента магнитного диполя спутника вдоль его оси вращения, Vq - магнитная постоянная для диполя в центре Земли, R - радиус орбиты, b -нормализованный вектор геомагнитного диполя.

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

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

Такая управляющая система использовалась на спутнике "Tiros-II" и всех последующих спутников, вплоть до "Tiros-VIII". Позднее такая же система стабилизации вращением была применена на связном спутнике-ретрансляторе "Relay" и спутнике "Telstar".

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

Система организации магнитного управляющего момента, посредством которого ось разворачивается на 90 градусов и аппарат переводится в режим «колеса», использует тот же самый тип катушки, что и в предыдущих конструкциях спутника, имеющих дипольный момент, направленный параллельно оси вращения. Если обеспечить реверсирование тока протекающего по катушке через каждую четверть орбиты (для круговой орбиты), то можно организовать прецессию в желаемом направлении. Такой прием позволил устранить ограничения, присущие схеме, используемой на ряде первых спутников. На "Tiros-IX" имеется дополнительная магнитная катушка, плоскость которой совпадает с осью вращения. Ее назначение состоит в том, чтобы создавать моменты которые можно было использовать для управления скоростью вращения спутника.

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

Поскольку система управления пространственным положением оси вращения содержит катушку, в которой производится изменение направления тока через каждую четверть орбиты, то ее принято называть сокращенно 8

QOMAC (the Quarter-Orbit Magnetic Attitude Control system). Аналогично система управления скоростью вращения спутника сокращенно называется MASC (the Magnetic Attitude Spin Control).

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

Магнитный момент в грубом режиме, который способен обеспечить прецессию со скоростью около 10 градусов за орбиту, оценивается величиной диполя равной 30 А м2. Величина диполя в точном режиме оценивается 6 Ам 2, а скорость, развиваемая моментом прецессии, около 2 градусов за орбиту. Более поздние спутники этой серии [15] использовали для ориентации помимо токовых катушек также маховики и жидкостные демпферы.

Приведем краткий обзор по системам ориентации, используемых на некоторых современных малых спутниках [56].

Oersted - датский микроспутник изготовлен кооперацией фирм и университетов во главе с CRI, имеет размеры 720x340x450 мм и массу 62 кг [63]. Корпус спутника имеет в разрезе форму буквы Н, где на одной стороне перекладины смонтированы блоки электроники, а на другой - механизм развертывания штанги гравитационной стабилизации. Система ориентации также включает: магнитные катушки, звездный датчик, солнечный датчик и магнитометр. Задачи, решаемые спутником: изучение и съемка магнитного поля Земли.

Наноспутник Берлинского университета Tubsat-N размерами 320x320x104 мм и массой 8.5 кг использует систему ориентации, состоящую из магнитных катушек, магнитометра, силового маховика и звездного датчика [57]. Основное назначение - слежение за перемещениями крупно- и среднеразмерными млекопитающими, угнанными машинами, сбор данных с радиобуев.

Sapphire - микроспутник Стентфортского Университета, имеет форму шестигранного цилиндра шириной 60 см и весит 18 кг [58]. Для его ориентации используется пассивная магнитная система, состоящая из постоянного магнита и гистерезисных стержней. Цель запуска спутника - проверить работоспособность в космосе нового инфракрасного датчика.

Sunsat - микроспутник ЮАР, создан в Стелленбосхском Университете, имеет размеры 45x54x62 см и массу 59 кг [64]. Система ориентации состоит из штанги гравитационной стабилизации, маховиков и магнитометров.

Starshine - микроспутник создан военно-морской исследовательской лабораторией США, имеет форму сферы диаметром 48 см и массой 39.3 кг [59]. Поверхность спутника покрыта 877 полированными алюминиевыми зеркалами диаметром 25 мм. Спутник ориентируется осью вращения вдоль линий магнитного поля с помощью постоянного магнита и закручивается посредством «фотонного флюгера» - пластин, части которых окрашены отражающей и поглощающей красками. Скорость вращения спутника около 1 об/мин.

Спутник TERRIERS разработан Бостонским университетом, служебный борт изготовлен Aero Astro Inc. [65]. Масса спутника 124 кг, служебный борт -75 кг. Спутник имеет прямоугольную форму. Его длина составляет 97 см, поперечный размер 51 см. Используется стабилизация вращением 1-10 об/мин. Система ориентации состоит из двух магнитометров, 6 грубых солнечных датчиков и магнитной катушки. Цель запуска - исследование ионосферы.

DLR-Tubsat - немецкий микроспутник [66]. Полезная нагрузка изготовлена институтом технологии космических датчиков, служебный борт -институтом аэронавтики и астронавтики технического университета Берлина. Имеет форму параллелепипеда 32x32x27 см. Его масса составляет 45 кг. Система ориентации включает звездный датчик, три волоконно-оптических лазерных гироскопа и три маховика. Когда спутник не ведет съемку, он стабилизируется собственным вращением.

ASUSat-1 - наноспутник Аризонского университета изготовлен в виде призмы с 14 боковыми гранями высотой 25.4 см и максимальным поперечным размером основания 34.3 см [67]. Масса спутника составляет 5.9 кг. Система определения ориентации состоит из фотодиодов (простейших солнечных датчиков), размещенных группами по 4 на каждой боковой грани корпуса (один датчик Солнца и три датчика Земли), и фотодиодов, размещенных на основаниях призмы. Система ориентации состоит из гравитационной штанги и жидкостного демпфера.

TiungSat-1 - малазийский микроспутник, построенный SSTL. Его масса составляет 54 кг [68]. Размеры спутника 69x36.6x36.6 см. Система ориентации имеет штангу гравитационной ориентации развертываемая в полете. Спутник оснащен четырьмя камерами для съемки поверхности Земли в разных диапазонах и с разным разрешением, цифровым процессором для обработки данных, оборудованием для мониторинга радиационной обстановки.

Munin - наноспутник, разработанный Шведским институтов космической физики и студентами университетов Умеа и Лулеа [2,60]. Масса спутника: 7.5 кг, размеры: 20x20x25 см. Этот спутник ориентируется по магнитному полю с помощью постоянного магнита. Спутник предназначен для изучения полярных сияний и «космической погоды» над северным и южным полушариями.

QuikToms - малый спутник, изготовленный компанией OSC для изучения и прогнозирования озонового слоя [61]. Масса спутника -162 кг. Спутник представляет собой два кольцевых модуля, расположенных рядом. Система ориентации включает датчик Земли, солнечный датчик, GPS приемник, гироскопический блок, три маховика, магнитометры и токовые катушки и обеспечивает стабилизацию по трем осям с точностью 0.5 град по направлению и 0.01 град/с по скорости, а также определение фактической ориентации с точностью 0.25 град.

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

11 магнитометра [62]. Спутник создан РЖИ РАН при участии школьников Обнинска и Сиднея. Масса спутника составляет 20.5 кг. Аппарат имеет форму шестиугольной призмы (высота 51.5 см, диаметр описанной окружности 37 см) с четырьмя раскрывающимися двусторонними панелями солнечных батарей. Система ориентации - магнитно-гравитационная.

Как видно из приведенного обзора, на малых спутниках довольно часто используются магнитные катушки или постоянные магниты, хотя другие подходы к обеспечению ориентации также имеют место. Помимо спутников, стабилизируемых собственным вращением, магнитное управление может быть использовано для спутников с маховиками [40]. На таких спутниках магнитные катушки часто используются для разгрузки маховичной системы [37]. В работе [41] рассматриваются возможности одновременного использования маховиков и магнитных катушек для минимизации расходуемой энергии при поддержании требуемой ориентации. С помощью активной магнитной системы можно осуществлять трехосную ориентацию спутника [44]. Для синтеза управления, осуществляющего трехосную ориентацию спутника, удобно рассматривать линеаризованную систему, как это сделано в работах [45,53]. В работе [54] вращательное движение спутника описывается линейной дискретной периодической системой, и осуществляется поиск оптимального управления.

Ранее на практике использовались различные способы управления токовыми катушками, для обеспечения заданного углового движения спутника, стабилизируемого собственным вращением [48]. Так, в [46] рассматривается способ управления, основанный на методе осреднения, в котором создаваемый механический момент осредняется на периоде управления, а закон управления не зависит от направления оси вращения. В этом случае закон управления изменяется так, чтобы управляющий магнитный момент обеспечивал заданную ориентацию. В [46] исследовались пять законов управления, причем для околополярных орбит выбиралось релейное управления с переключением через каждую четверть орбиты. В [38,52] разработан закон управления для круговой

12 орбиты. В [38] предложен критерий максимума величины магнитного момента и ошибки углового момента. В работе [50] использован фильтр Калмана для минимизации угловой ошибки и создан метод управления с минимальными энергетическими затратами. В работе [16] рассматривается возможность управления как скоростью вращения, так и ориентацией с помощью одной катушки.

Способ управления угловым движением спутника, стабилизируемого собственным вращением, который предложен в работах [35,47], основан на анализе величины рассогласования между номинальным (расчетным) и текущим кинетическим моментом спутника. Этот подход в настоящее время предполагается использовать для ориентации американского малого спутника LionSat [43] и бразильского спутника BDCS (Brazilian Data Collecting Satellite) [39]. В диссертационной работе также используется такой подход для управления ориентацией и скоростью вращения спутника [23]. Алгоритмы управления ориентацией и скоростью вращения являются частями автономного алгоритма управления угловым движением спутника, использующего информацию об угловой скорости спутника.

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

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

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

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

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

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

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

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

15 за определенные интервалы времени. Сеансы связи происходили тогда, когда спутник ТНС-0 находился в подходящей точке пространства и ориентировался нужным образом по отношению к спутникам системы GLOBALSTAR, через которую осуществлялась связь спутника с Землей.

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

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

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

Заключение диссертации по теме «Теоретическая механика», Ильин, Андрей Александрович

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

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

2. Исследовано влияние гравитационного момента на быстро закрученный спутник. Аналитически оценивается величина минимального дипольного момента необходимого для компенсации возмущения, производимого гравитационным моментом. Проведено аналитическое исследование динамики спутника, стабилизируемого собственным вращением, для предложенных алгоритмов управления. Объяснен возникающий эффект перерегулирования. Предложены способы преодоления этого эффекта. Проведено численное моделирование вращательного движения на примере наноспутника ТНС-1 с магнитной системой ориентации для обоих алгоритмов управления.

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

Список литературы диссертационного исследования кандидат физико-математических наук Ильин, Андрей Александрович, 2006 год

1. Артюхин Ю.П., Каргу Л.И., Симаев В.Л., Системы управления космических аппаратов, стабилизированных вращением, М.: Наука, 1979, - 295 с.

2. Барабаш С.В., Кирюшкин И.Ю., Норберг У., Овчинников М.Ю., Пеньков В.И., Методы определения углового движения наноспутника Munin, М.: Препринт ИПМ им.М.В.Келдыша РАН, 2003, № 24, 24 с.

3. Белецкий В.В., Движение искусственного спутника относительно центра масс, М.: Наука, 1965, 416 с.

4. Белецкий В.В., Движение спутника относительно центра масс в гравитационном поле, М.: Изд-во Моск. ун-та, 1975, 308 с.

5. Белецкий В.В., Новогребелъский А.Б., Существование устойчивых относительных равновесий искусственного спутника в модельном магнитном поле //Астрономический журнал, 1973, т. 50, с. 327-335.

6. Белецкий В.В., Хентов А.А., Вращательное движение намагниченного спутника, М.: Наука, 1985, 288 с.

7. Биркхольд Е., Лейб Ф., Стопфкюхен К, Магнитная система активной стабилизации ориентированных на Солнце вращающихся спутников, Управление в космосе, Труды ИФАК, вып.2, М.: Наука, 1972, с. 223-230.

8. Гехт Е., Мендлсер У.П., Магнитная система управления угловым положением спутников серии «Тайрос», В сб.: Проблемы ориентации искусственных спутников Земли, М.: Наука, 1966, с. 336-350.

9. Гончарский А.А., Хентов А.А., О некоторых режимах вращения намагниченного спутника гиростата в геомагнитном поле // Радиофизика, 1972. т. 15,№ 11, сЛ 631-1636.

10. Гребенников Е.А., Метод усреднения в прикладных задачах, М.: Наука, 1986,-256 с.

11. Ильин А.А., Куприянова Н.В., Овчинников М.Ю., Пеньков В.К, Селиванов А.С., Анализ вращательного движения первого российского наноспутника ТНС-0 по результатам летных испытаний, М.: Препринт ИПМ им.М.В.Келдыша РАН, 2006, № 18, 28 с.

12. Ильин А.А., Овчинников М.Ю., Пеньков В.И., Обеспечение ориентации малого спутника, стабилизируемого собственным вращением, М.: Препринт ИПМ им.М.В.Келдыша РАН, 2004, № 83, 28 с.

13. Ильин А.А., Овчинников М.Ю., Пеньков В.И., Алгоритмы магнитной системы ориентации малого спутника, стабилизируемого собственным вращением, М.: Препринт ИПМ им.М.В.Келдыша РАН, 2005, № 19, -32 с.

14. Ковалик X., Система управления вращением и ориентацией спутников ISIS-I и ISIS-B, Управление в космосе, Труды ИФАК, вып.2, М.: Наука, 1972, с. 183-193.

15. Куприянова Н.В., Овчинников М.Ю., Пеньков В.И, Селиванов А.С., Пассивная магнитная система ориентации первого российского наноспутника ТНС-0, Препринт ИПМ им.М.В.Келдыша, 2005, №46, -23 с.

16. Моисеев Н.Н., Асимптотические методы нелинейной механики. М.: Наука, 1971,-400 с.

17. Непобедимый С.П., Белов А.А., Ильин А.А., Калинин А.П., Овчинников М.Ю., Орлов А.Г., Родионов И.Д., Родионов А.И., Трехмерное техническое зрение на основе монофотонных технологий // Доклады Академии наук, 2006, т.406, № 3, с. 333-336.

18. Новожилов И.В., Фракционный анализ, М.: Изд-во механико-математического факультета МГУ, 1995, 224 с.

19. Нохрина Е.Е., Исследование возможности демпфирования колебаний спутника, В сб.: Обработка информации и моделирование, М.: МФТИ, 2002, с. 119-122.

20. Овчинников М.Ю., Пеньков В.И., Ильин А.А., Селиванов А.С., Алгоритмы управления ориентацией российского наноспутника ТНС-1 // Актуальные проблемы авиационных и аэрокосмических систем: процессы, модели, эксперимент, 2005, т. 10, № 1(19), с. 78-101.

21. Пивоваров М.Л., Жидкостное демпфирование колебаний спутника с большим магнитным моментом // Космич. исслед. 1990, т 28, № 6, с. 865-873.

22. Плотников В. А., Зверкова Т.С., Метод усреднения для систем стандартного вида с разрывными правыми частями. // Диффиренц. уравнения, 1982, т. 18, № 6, с. 1091-1093.

23. Раушенбах Б.В., Овчинников М.Ю., Лекции по динамике космического полета: учебное пособие, М.: МФТИ, 1997, 188 с.

24. Сарычев В.А., Овчинников М.Ю., Магнитные системы ориентации искусственных спутников Земли // Итоги науки и техники, Сер.: Исследование космического пространства. М.: ВИНИТИ, т. 23, 1985, -104 с.

25. Сидоренко В.В., Разрушение адиабатических инвариантов на резонансах: пример из динамики твердого тела, М.: Препринт ИПМ им.М.В.Келдыша РАН, 1995,№76,-24 с.

26. Сидоренко В.В., Об одном классе движений спутника, несущего сильный магнит // Космич. исслед. 2002, т 40, № 2, с. 147-155.

27. Тихонов А.А., Уточнение модели «наклонный диполь» в задаче об эволюции ротационного движения заряженного тела в геомагнитном поле // Космич. исслед. 2002, т 40, № 2, с. 171-177.

28. Фишелл Р.Э., Стабилизация вращения спутников, Автоматическое управления космическими летательными аппаратами, Труды I Международного симпозиума ИФАК по автоматическому управлению в мирном использовании космического пространства, М.: Наука, 1968.

29. Фрелих X, Ганзенхаузер Р., Калман X., Лейб Ф., Активная магнитная система управления для спутника «Аэрос», Управление в пространстве, Труды IV Международного симпозиума ИФАК по автоматическому управлению в пространстве, М.: Наука, 1973, с. 103-112.

30. Шегехара М, Выбор законов управления угловым положением вращающихся спутников, Управление в пространстве, Труды V Международного симпозиума ИФАК по автоматическому управлению впространстве, т.1, М.: Наука, 1975, с.238-249.144

31. Camillo P., MarkleyF., Orbit-Averaged Behavior of Magnetic Control Laws for Momentum Unloading // Journal of Guidance and Control, v.3, № 6, 1980, pp. 563-568.

32. Ergin E.I., Wheeler P.C., Magnetic Attitude Control of a Spinning Satellite // Journal of Spacecraft and Rockets, 1965, v.2, № 6, pp.846-850.

33. Ferreira D.D., Crus J.J., Attitude and Spin Rate Control of a Spinning Satellite Using Geomagnetic Field // Journal of Guidance, Control, and Dynamics, 1991, v. 14, № 1, pp.216-218.

34. Goel P.S., Rajaram S., Magnetic Attitude Control of a Momentum-Biased Satellite in Near-Equatorial Orbit // Journal of Guidance, Control, and Dynamics, 1979, v.2, № 4, pp.334-338.

35. Grassi M., Pastena M., Minimum Power Optimum Control of Microsatellite Attitude Dynamics // Journal of Guidance, Control, and Dynamics, 2000, v. 23, № s, pp.798-804.

36. Hughes P.C., Spacecraft Attitude Dynamics, New York, Dover Publications, 2004, 570 p.

37. Hur Ph., Melton R., Spencer D., Meeting Science Requirements For Attitude Determination and Control in a Low-Power, Spinning Nanosatellite, Paper IAC-04-IAF-A.4.05 at the 55th Congress IAF, Vancouver, Canada, 4-8 Oct., 2004, 9 p.

38. Lovera M., Astolfi A., Global Magnetic Attitude Control of Inertially Pointing Spacecraft // Journal of Guidance, Control, and Dynamics, 2005, v.28, № 5, pp. 1065-1072.

39. Psiaki M.L., Magnetic Torquer Attitude Control via Asymptotic Periodic Linear Quadratic Regulation // Journal of Guidance, Control, and Dynamics, 2001, v. 24, №2, pp.3 86-3 94.

40. Renard M.L., Command Laws for Magnetic Attitude Control of Spin-Stabilized Earth Satellites I I Journal of Spacecraft and Rockets, 1967, v.4, №2, pp. 156-163.

41. Shigehara M., Geomagnetic Attitude Control of an Axisymmetric Spinning Satellite, Journal of Spacecraft and Rocket, 1972, v.9, № 6, pp 391-398.

42. 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, v. 6, № 6, pp. 461-471.

43. Sidi M.J., Spacecraft Dynamics and Control. A practical Engineering Approach, UK, Cambridge University Press, 1997, 409 p.

44. Sorensen J.A., A Magnetic Three-Degree of Freedom Attitude Control System for an Axisymmetric Spinning Spacecraft // Journal of Spacecraft and Rockets, 1971, v. 8, №5, pp. 441-448.

45. Vallado D.A., McClain W.D., Fundamentals of Astrodynamics and Application., El Segundo California: Microcosm Press, 2001, 958 p.

46. Wheeler P.C., Spinning Spacecraft Attitude Control via the Environmental Magnetic Field // Journal of Spacecraft and Rockets, 1967, v. 4, № 12, pp 1631-1637.

47. Wisniewski R., Linear Time-Varying Approach to Satellite Attitude Control Using Only Electromagnetic Actuation // Journal of Guidance and Control, 2000, v. 23, № 4, pp. 640-647.

48. Wisniewski R., Stoustrup J., Periodic H2 Synthesis for Spacecraft Attitude Control with Magnetorquers // Journal of Guidance and Control, 2004, v. 27, №5, pp. 874-881.

49. Zajac E.E., Some simple solutions to magnetic attitude control of satellites. Proc. 4th U.S. Nat. Congr. Appl. Mech., Berkeley, Calif., 18-21 June, 1962,