Способ восстановления ориентации орбитального космического аппарата

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

Известен способ управления ориентацией, приведенный в статье авторов О.В. Ермошиной, А.П. Крищенко «Синтез программных управлений ориентацией КА методом обратных задач динамики», Москва, МГТУ им. Н.Э. Баумана, 2010 г., 12 с. Недостатком способа является затянутость переходного процесса, связанного с выбором «гладкой» переходной функции, близкой к асимптотическому переходному процессу.

Известен способ ВО, приведенный в книге авторов В.Н. Бранец, И.П. Шмыглевский «Применение кватернионов в задачах ориентации твердого тела». Москва, Наука 1973 г., 320 с. (см. стр. 205-226), где также рассматриваются асимптотические процессы приведения ССК к ОСК.

Известен способ ВО, изложенный в статье «Система ориентации и стабилизации космического аппарата по информации с астродатчиков», Электронный журнал «Труды МАИ». Выпуск №38, который можно принять за прототип. Этот способ включает начальное демпфирование угловых скоростей КА относительно инерциальной - геоцентрической абсолютной системы координат (ИСК), расчет по данным аппаратуры спутниковой навигации (АСН) матрицы А положения орбитальной системы координат (ОСК) относительно ИСК, измерение с помощью астродатчика (АД) матрицы М_ro ориентации КА относительно ИСК с последующим расчетом матрицы ориентации S связанных осей КА относительно ОСК по формуле:

К недостатку прототипа также следует отнести повышенное время ВО. С целью сокращения времени на ВО известный способ дополняют новыми операциями. Для этого в задемпфированном положении КА относительно ИСК, фиксируют положение ОСК на момент времени t_Φ и рассчитывают проекции и модуль собственного вектора Е матрицы S на оси ОСК-Ф по формулам:

где Е_X, E_Y, E_Z - проекции собственного вектора ориентации КА на соответствующие оси ОСК;

S_ij - элементы матрицы S;

E_mod - модуль собственного вектора матрицы S,

рассчитывают угол Эйлера (угол конечного поворота) для приведения КА к ОСК-Ф, зафиксированной относительно ИСК на время t_Φ, по формуле:

где angV_O - угол Эйлера.

Задают программную скорость ω_ПР вращения КА вокруг собственного вектора до совпадения связанных осей КА с осями ОСК-Ф, при этом программные скорости по осям КА рассчитывают по формулам:

а время конечного поворота по формуле:

где К<1 - коэффициент упреждения (определяется экспериментально). В течение расчетного времени t_ПР поворачивают КА одновременно по каждой оси ориентации с угловыми скоростями , , и по истечении времени программного поворота t_ПР демпфируют угловые скорости КА до выполнения условий:

где ω_o - допустимое значение угловой скорости КА,

вновь определяют матрицу ориентации S, по элементам которой находят углы ориентации КА относительно ОСК, например, для последовательности поворотов ψ, ϑ, γ по формулам:

проверяют условия приведения углов ориентации КА к наперед заданным значениям:

где ψ_ЗАД, ϑ_ЗАД, γ_ЗАД - заданные значения углов ориентации КА,

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

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

На фиг. 1 показана кинематическая схема ВО, где обозначено:

1 - КА;

НП - направление полета;

X₀Y₀Z₀ - ОСК;

X_SY_SZ_S - ССК;

angV_O - угол Эйлера (угол конечного поворота);

Е_X, Е_Y, E_Z - проекции собственного вектора ориентации КА на оси ОСК, равные проекциям на оси ССК;

S_ij - элементы матрицы S;

E_mod - модуль собственного вектора матрицы S.

На фиг. 2 обозначено:

1 - участок демпфирования на маховиках;

2 - участок программного поворота;

3 - участок включения астрокоррекции;

4 - участок текущей астроориентации;

TSS - текущее состояние системы.

На фиг. 3 обозначено:

ψ, ϑ, γ - углы ориентации КА относительно ОСК по курсу, тангажу и крену соответственно.

При использовании АД для ВО необходимо снизить (задемпфировать) угловые скорости КА по связанным осям до скоростей ω_o, при которых АД становится работоспособным и выдает в систему ориентации положение КА относительно ИСК в виде кватерниона или матрицы - M_ro.

После того как КА приобретает скорости, меньшие ω_o, по показаниям АД и навигационной системы рассчитывают матрицу S ориентации КА относительно ОСК:

где A^T транспонированная матрица ориентации ОСК относительно ИСК, получаемая из навигационной системы КА. По элементам матрицы S_ij для фиксированного момента времени - t_Φ рассчитывают проекции собственного вектора (вектора конечного поворота) Е на связанные оси КА, которые буду равны также проекциям собственного вектора Е на оси ОСК, и угол конечного поворота для приведения КА к ОСК-Ф (Фиг. 1), т.е. к фиксированному положению ОСК относительно ИСК на момент времени t_Φ. Расчеты проводят по формулам (3, 5). После этого задаются величиной программной скорости , с которой будут вращать КА относительно вектора конечного поворота до совмещения ССК с ОСК-Ф.

Так как вращение осуществляется с помощью ИО, расположенных по связанным осям, то необходимо найти проекции на связанные оси КА, которые будут являться программными скоростями КА относительно каждой связанной оси КА. Это можно сделать, вычислив косинусы углов между и его проекциями на связанные оси КА (4). Время конечного поворота t_ПР рассчитывается по формуле (5), которое уменьшают на некоторую величину времени путем введения К (коэффициент упреждения), достаточную для естественного демпфирования (торможения) КА до заданных скоростей ω_o по условиям (6).

Поворачивают КА со скоростями (4), например, путем введения в закон стабилизации КА по угловой скорости смещений на величины программных скоростей (8):

По окончании времени t_ПР проверяют выполнение условий (8). Переходный процесс показан на фиг. 3. Если углы ориентации не будут находиться в заданных пределах (8), то операции ВО повторяют, пока не будут выполнены одновременно условия (6) и (8), после чего переводят систему в режим астроориентации КА по показаниям астродатчика.