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

Изобретение относится к определению массово-инерционных характеристик космических аппаратов (КА). Способ включает ориентацию КА и стабилизацию в инерциальной системе координат (ИСК) его строительной оси, ближайшей к оси максимального момента инерции. Далее выполняют закрутку КА вокруг этой оси с угловой скоростью не менее 2°/с. Измеряют в системе строительных осей КА направления на регистрируемые звезды и угловую скорость КА до определённого момента времени. Последний зависит от времени закрутки КА и интервала движения КА, слабо возмущенного действием гравитационного градиента и вычисляемого с некоторым коэффициентом надежности. Опознают указанные звезды и определяют в ИСК направления на них. Тензор инерции КА определяют по указанным направлениям на звезды и значениям угловой скорости КА. Технический результат изобретения заключается в повышении достоверности определении тензора инерции КА, в т.ч. при отсутствии на его борту инерционных исполнительных органов.

 

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

Тензор инерции любого твердого тела является важной характеристикой для управления его движением. Поэтому был разработан ряд способов для определения тензора инерции тела, описанных, например, в [1] (Способ определения тензора инерции и координат центра масс тела и устройство для его осуществления, патент RU 2348020 C1). В способе-аналоге [1] телу сообщается заданное движение и по измерениям параметров движения определяется тензор инерции тела. Главный недостаток способа [1] и других аналогичных способов заключается в отсутствии возможности их применения для определения тензора инерции КА в полете.

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

Для определения тензора инерции КА в полете был предложен способ [2] (Севастьянов Н.Н., Бранец В.Н., Банит Ю.Р., Беляев М.Ю., Сазонов В.В. «Определение тензора инерции геостационарных спутников «Ямал» по телеметрической информации. Препринт ИПМ им. М.В. Келдыша №17, 2006 г.). Предложенный способ [2], взятый авторами за прототип, включает инерциальную ориентацию и развороты КА и измерение суммарного кинетического момента маховиков. При изменении ориентации КА путем его разворотов, по измерениям суммарного кинетического момента маховиков (инерционных исполнительных органов) определяется тензор инерции КА в полете.

Недостаток способа-прототипа связан с низкой точностью определения тензора инерции КА и необходимостью использования измерений от инерционных исполнительных органов (ИИО) [2]. В то же время многие КА не имеют в своем составе ИИО. Например, транспортный грузовой корабль (ТГК) «Прогресс», являющийся основным транспортным грузовым кораблем в программе МКС, не имеет в своем составе ИИО. Вместе с тем, за счет перемещения грузов космонавтами внутри ТГК «Прогресс» и расхода большого количества топлива на ТГК, его тензор инерции меняется в процессе полета. Особенно важно знание углового рассогласования главных осей инерции ТГК и его строительных осей, т.к. двигатели ориентации и коррекции ТГК установлены относительно строительных осей корабля.

Задачей, на решение которой направлено настоящее изобретение, является определение тензора инерции КА в полете.

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

Технический результат достигается тем, что в способе определения тензора инерции космического аппарата в полете, включающем инерциальную ориентацию и развороты космического аппарата, ориентируют космический аппарат, стабилизируя в инерциальной системе координат его строительную ось, ближайшую оси, соответствующей максимальному моменту инерции, выполняют закрутку космического аппарата вокруг этой строительной оси с угловой скоростью Ω₂ не менее 2°/с, измеряют в строительной системе координат космического аппарата направления на регистрируемые звезды и угловую скорость космического аппарата до момента времени

T=T₀+Δt, где $$\Delta t=\frac{4{\Omega }_2\cdot R^3}{3\cdot {\mu }_{гр}\cdot K}$$ ,

где T₀ - момент времени выполнения закрутки космического аппарата;

Δt - интервал времени слабо возмущенного движения космического аппарата;

Ω₂ - угловая скорость закрутки вокруг строительной оси, ближайшей оси максимального момента инерции;

R - радиус орбиты;

µ_гр - гравитационный параметр Земли;

K - коэффициент надежности,

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

За счет выполнения предлагаемых действий определение тензора инерции КА осуществляется надежно и даже при отсутствии на борту КА ИИО. Действия способа обеспечивают слабо возмущенное движение КА на интервале времени Δt. Это позволяет надежно определять тензор инерции КА даже при отсутствии на его борту ИИО. На угловое движение КА оказывают влияние, в основном, гравитационный и аэродинамический возмущающие моменты, причем основное влияние на большинство КА оказывает гравитационный момент.

Соотношение для Δt получено для КА, имеющего вытянутую форму, с учетом действия на него гравитационного возмущающего момента. При выводе соотношения для Δt учитывается максимальное значение гравитационного момента, действующего вокруг поперечной оси КА. Для повышения надежности обеспечения слабо возмущенного движения на интервале времени Δt вводится специальный коэффициент надежности К. Коэффициент надежности может быть взят равным, например, 10. Для ТГК «Прогресс», например, Δt оказывается равным нескольким десяткам минут. На этом интервале времени угловое движение КА считается невозмущенным. Тензор инерции космического аппарата в этом случае по измеренным и определенным параметрам определяют минимизацией функционала

$${Ф}_{\Omega }={\displaystyle \sum _{n=1}^N{\displaystyle \sum _{i=1}^3{\left[{\Omega }_i^{\left(n\right)}-{\Omega }_i\left(t_n\right)\right]}^2}}$$ ,

на решениях системы уравнений (уравнений Эйлера, записанных в безразмерном виде)

$${\dot{\omega }}_1=\mu {\omega }_2{\omega }_3$$ , $${\dot{\omega }}_2=\frac{{\mu }^{\text{'}}-\mu }{1-\mu {\mu }^{\text{'}}}\cdot {\omega }_3{\omega }_1$$ , $${\dot{\omega }}_3=-{\mu }^{\text{'}}{\omega }_1{\omega }_2$$ ,

где: $$\mu =\frac{I_2-I_3}{I_1}$$ , $${\mu }^{\text{'}}=\frac{I_2-I_1}{I_3}$$ , $${\Omega }_i={\displaystyle \sum _{k=1}^3{в}_{iк}{\omega }_k\left(i=\mathrm{1,}\mathrm{2,}3\right)}$$ ,

ω₁, ω₂, ω₃ - компоненты угловой скорости на главные центральные оси инерции;

I₁, I₂, I₃ - моменты инерции космического аппарата;

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

$${\Omega }_i^{\left(n\right)}$$ - приближенные измеренные значения компонент угловой скорости в строительной системе координат.

Минимизация Ф_Ω является первым этапом определения искомых величин и осуществляется методом Гаусса-Ньютона.

Ф_Ω рассматривается как функция набора из восьми параметров ω_iо=ω_i(t_о) (i=1, 2, 3), µ, µ′, γ, α, β. Углы γ, α, β задают положение строительной системы координат оу₁у₂у₃ относительно системы координат ох₁х₂х₃, образованной главными центральными осями инерции КА.

Система оу₁у₂у₃ может быть переведена в систему ох₁х₂х₃ тремя последовательными поворотами: 1) на угол α вокруг оси оу₂, 2) на угол β вокруг новой оси оу₃, 3) на угол γ вокруг новой оси оу₁, совпадающей с осью ох₁.

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

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

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

Наиболее ценным для управления движением КА является точное знание элементов матрицы в_ik (т.е. углов γ, α, β). Это обеспечивается выполнением всей совокупности действий и приемов способа.

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

В настоящее время технически все готово для реализации предложенного способа, например, на ТГК «Прогресс» или других КА. На ТГК «Прогресс» отсутствуют ИИО. Однако система управления ТГК «Прогресс» позволяет выполнять инерциальную ориентацию, развороты и закрутку КА. Для измерения направлений на звезды может использоваться звездный датчик типа БОКЗ или ОЗД. Звезды, попавшие в поле зрения датчика, регистрируются в зависимости от заложенной в прибор яркости (могут регистрироваться звезды, например, до 6-ой звездной величины). Опознавание попавших в его поле зрения звезд осуществляется автоматически (по яркости звезд и угловому расстоянию между зарегистрированными звездами). На ТГК измеряются угловые скорости в строительной системе координат корабля, направление на Солнце (которое, строго говоря, является звездой). Для определения необходимых направлений и вычислений, ТГК снабжен бортовой вычислительной системой БВС.

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

ЛИТЕРАТУРА

1. Способ определения тензора инерции и координат центра масс тела и устройство для его осуществления, патент RU 2348020 С 1.

2. Севастьянов Н.Н., Бранец В.Н., Банит Ю.Р., Беляев М.Ю., Сазонов В.В. «Определение тензора инерции геостационарных спутников «Ямал» по телеметрической информации. Препринт ИПМ им. М.В. Келдыша №17, 2006 г.

Способ определения тензора инерции космического аппарата в полете, включающий инерциальную ориентацию и развороты космического аппарата, отличающийся тем, что ориентируют космический аппарат, стабилизируя в инерциальной системе координат его строительную ось, ближайшую к оси, соответствующей максимальному моменту инерции, выполняют закрутку космического аппарата вокруг этой строительной оси с угловой скоростью Ω₂ не менее 2 град/с, измеряют в строительной системе координат космического аппарата направления на регистрируемые звезды и угловую скорость космического аппарата до момента времени
T = T₀ + Δt, $$\Delta t=\frac{4{\Omega }_2\cdot R^3}{3\cdot {\mu }_{гр}\cdot K}$$ ,
где T₀ - момент времени выполнения закрутки космического аппарата,
Δt - интервал времени слабо возмущенного движения космического аппарата,
Ω₂ - угловая скорость закрутки вокруг строительной оси, ближайшей оси максимального момента инерции,
R - радиус орбиты,
µ_гр - гравитационный параметр Земли,
K - коэффициент надежности,
опознают зарегистрированные звезды, определяют в инерциальной системе координат направления на опознанные звезды и определяют тензор инерции космического аппарата по измеренным и определенным на интервале времени Δt направлениям на опознанные звезды и измерениям угловой скорости космического аппарата.