CC BY
195
УДК 629.78 В01: 10.18287/2412-7329-2016-15-2-50-56
БОРТОВОЙ АЛГОРИТМ ДЛЯ СИСТЕМЫ ОРИЕНТАЦИИ И СТАБИЛИЗАЦИИ НАНОСПУТНИКА SAMSAT-218Д
© 2016 А. В. Крамлих1, М. Е. Мельник2
1Самарский национальный исследовательский университет имени академика С.П. Королёва 2Акционерное общество «РКЦ «Прогресс», г. Самара
В статье описываются алгоритмы определения ориентации и демпфирования угловых скоростей наноспутника SamSat-218Д, разработанного в Самарском государственном аэрокосмическом университете. В качестве источников информации в системе ориентации используются трёхосный магнитометр и датчики освещённости. Исполнительными органами, обеспечивающими демпфирование угловых скоростей на борту наноспутника, являются три магнитные катушки. В качестве алгоритма определения ориентации применяется известный метод векторного согласования с оригинальным критерием предварительной отбраковки измерений, поскольку метод векторного согласования имеет высокую чувствительность к ошибкам измерений. В качестве алгоритма демпфирования используется традиционный алгоритм B-dot. Изучено влияние начальных угловых скоростей на коэффициент магнитного момента в алгоритме управления B-dot. Поскольку непрерывный режим работы магнитных катушек на борту наноспут-ника в орбитальном полёте нецелесообразен из-за необходимости измерения вектора напряжённости магнитного поля Земли, то проведён анализ процесса демпфирования в дискретном режиме работы при различных величинах начальных угловых скоростей отделения наноспутника. Приведено сравнение быстродействия алгоритма в непрерывном и дискретном режимах работы при различных начальных угловых скоростях отделения.
Наноспутник, ориентация, демпфирование углового движения, трёхосный магнитометр, датчик Солнца, магнитная катушка, алгоритм.
Введение
В настоящее время в Самарском государственном аэрокосмическом
университете завершена работа по созданию наноспутника стандарта CubeSat3USamSat-218Д [1]. Основным назначением наноспутника является изучение динамики его движения, отработка алгоритмов определения ориентации и демпфирования углового движения с использованием встроенных измерительных и исполнительных средств, входящих в состав бортовых модулей [1].
Ограничения, накладываемые принятым составом измерительных и исполнительных средств, а также характеристиками бортового компьютера формируют жёсткие требования на вычислительную сложность алгоритмов определения ори-
ентации и демпфирования углового движения.
Системы координат, используемые для решения задачи определения ориентации и стабилизации
При решении задач определения ориентации и стабилизации используются следующие системы координат [1,2]:
1. Орбитальная система координат (ОСК) 0ХУ2 с началом в центре масс наноспутника; ось 02 направлена от притягивающего центра; ось 01 совпадает с направлением вектора производной истинной аномалии (направление по бинормали к траектории центра масс); ось 0Х дополняет СК до правой.
Цитирование: Крамлих А.В., Мельник М.Е. Бортовой алгоритм для системы ориентации и стабилизации наноспутника SamSat-218Д // Вестник Самарского государственного аэрокосмического университета имени академика С.П. Королёва (национального исследовательского университета). 2016. Т. 15, № 2. С. 50-56. DOI: 10.18287/2412-7329-2016-15-2-50-56
2. Связанная система координат (ССК) Oxyz с началом в центре масс наноспутника, расположенная таким образом, что центробежный момент инерции спутника Iyz = 0, ось Ox является продольной осью наноспутника.
Матрицу перехода от ОСК к ССК представим в виде:
Bi = bj,(i, j = 13), bii = eos a,
b12 = sinan sin Ya , b13 =- Sin an C0S Ya , Ь21 = Sinan Sin Pn ,
Ь22 = C0S Pn C0S Ya - C0S an Sin Pn Sin Ya , Ь23 = C0S Pn Sin Ya + C0S an Sin Pn C0S Ya , b31 = Sin an C0S Pn ,
b32 = - Sin Pn C0S Ya - C0S an C0S Pn Sin Ya , Ь33 = - Sin Pn Sin Ya + C0S an C0S Pn C0S Ya ■
Критерии отбраковки измерительной информации
В алгоритме определения ориентации наноспутника используются следующие критерии отбраковки измерительной информации.
Для отбраковки магнитометрических измерений необходим расчёт вектора напряжённости магнитного поля Земли (МПЗ) по выбранной модели. После расчёта модельных значений вектора напряжённости МПЗ отбраковываются измеренные векторы напряжённости с использованием соотношений:
8Иприб - инструментальная погрешность; 5Имод - погрешность модели МПЗ.
Для отбраковки информации о токосъёме воспользуемся критерием, аналогичным критерию отбраковки магнитометрических измерений:
Ризм SMOÓ - K2^S ,
5S =(5Sприб + 5Sмод ) =
где 8изм - измеренное значение тока в панелях СБ; 8мод - модельное значение тока в панелях СБ; 88 - допустимое отклонение; К2 - коэффициент запаса (К2 > 1); 88приб - инструментальная погрешность; 88мод - погрешность модели токосъёма.
Алгоритм определения ориентации
В качестве алгоритма определения ориентации выбран двухвекторный алгоритм определения углового положения наноспутника [2]. На каждый момент времени расчёта ориентации имеем составляющие векторов 8 и И в ССК и ОСК. Составляющие единичных векторов 8 и И в ССК и ОСК связаны соотношениями:
^x Sx Hx H X
Sy = B1 SY , Hy =B1 HY
^z SZ Hz Hz
Низм - HмоЛ - ^SH ,
5H = (5H
+ SH
приб мод
) ■
где Иизм - измеренное значение модуля вектора напряжённости МПЗ; Имод - модельное значение модуля вектора напряжённости МПЗ; 8И - допустимое отклонение; К1 - коэффициент запаса (К1 > 1);
Р
H
Введём в рассмотрение орты H х S х H H х S
q
H х S
H х S
Матрицы перехода М1, М2 от вспомогательной системы координат Opqr соответственно к осям ССК и ОСК имеют вид
r
M T
1
|H X S 4 a11 = Hx H x S|, a12 = Hy H x S|, a13 = Hz H x Si,
au j = 1,3)
a
a
22
a
23
= Sx - Hx (H, S),
= Sy - Hy (H, S),
= Sz - Hz (H, S),
a3i = HySz - HzSy, a32 = HzSx - HxSz , a33 = HxSy - HySx ,
m T
1
|H x S j c11 = HX (H x S|: c12 = HY |H x S|, c13 = HZ |H x SI,
cti ,(i, j = 1,3),
22
23
Запишем систему уравнений движения наноспутника в общем виде в осях ССК [3]:
dK0 dt
+ с x K0 = M0,
(1)
= SX - HX (Н S )-
= SY - HY (h, s),
= SZ - HZ (н, S),
С31 = HYSZ - HZSY , С32 = HZSX - HXSZ , С33 = HXSY - HY SX ■
Используя матрицы М1 и М2, найдём матрицу перехода от ОСК к ССК.
в 2 = МХМ Т.
Углы ап, уа находятся с помощью матриц в1 и в 2 по известным тригонометрическим соотношениям.
Модель движения относительно центра масс наноспутника
Рассмотрим силы и моменты, действующие на наноспутник.
При орбитальном полёте наноспут-ника по заданной орбите основное влияние на него оказывают гравитационный и аэродинамический возмущающие моменты [1].
где сс - абсолютная угловая скорость; K0 = IC - вектор кинетического момента;
I - тензор инерции; M0 - главный момент внешних сил.
Поскольку оси ССК являются главными центральными осями инерции нано-спутника для точки O, то векторное уравнение (1) в проекциях на данные оси запишется в виде:
IxO x + (Iz - Iy )CyCz = Mxe + Mxy , Iy<Cy + (Ix - Iz )<z<x = Mye + My ,
IzG)z + (Iy - Ix << = Mze + My ,
где Mxe, Mye, Mze - проекции главного
момента внешних сил на оси ССК; M ,
xy
Myy, My - проекции управляющего момента на оси ССК; «cx, а>у , «cz - проекции
угловых скоростей сс на оси ССК.
Учтём, что орбита круговая, и дополним динамические уравнения тремя кинематическими уравнениями, связывающими первые производные по времени Эйлеровых углов an, (n, уа (угол нутации - пространственный угол атаки ап ; угол собственного вращения - аэродинамический угол крена (п ; угол прецессии
- скоростной угол крена уа ) с проекциями вектора угловых скоростей на связанные оси «)x, са, G)z :
<x = Ya C0s an + (n + <0b21,
<y =/a SÏn (n SÏn an + an COS (n + <0b22 , < = Y COS (n SÍn an - an SÍn (n + <0b23 ,
С
21
где о0 - орбитальная скорость наноспут-ника.
В алгоритме демпфирования угловых скоростей наноспутника дипольный магнитный момент М формируется по закону [4; 5]
М = -кИ,
где к - постоянный положительный коэффициент, И - вектор производной напряжённости МПЗ.
Проекции дипольного магнитного момента на оси ССК имеют вид:
Му = - к (И у и г - И ж И у ), Муу = - к (И ж и х - И х и ж), Му = - к (И х иу - Иу и х),
где Их, И , И2 - проекции вектора
напряжённости МПЗ на оси ССК; Их, И ,
И 2 - проекции вектора производной напряжённости МПЗ на оси ССК.
Логика работы алгоритма заключается в следующем. По информации о направлении вращения вектора напряжённости МПЗ в ССК формируется управляющий момент, вращающий нано-спутник в обратном направлении.
Коэффициент к соответствует максимальному дипольному моменту магнитной катушки. В зависимости от начальных угловых скоростей отделения наноспутника коэффициент к может принимать различные значения (рис. 1), причём предполагается, что закрутки вокруг продольной оси наноспутника нет, т.е. (( = 0,1 °/с. При запуске наноспутника в
бортовой компьютер по умолчанию закладывается коэффициент для следующих начальных угловых скоростей: ах = 0,1 °/с, ( = 10°/с, = 10°/с. В дальнейшем он
У 2
может быть изменён в зависимости от фактических скоростей отделения нано-спутника.
Рис. 1. Величина коэффициента к при различных начальных угловых скоростях отделения наноспутника
При решении задачи демпфирования начальных угловых скоростей наноспут-ника необходимо учитывать следующее ограничение: на наноспутнике невозможна одновременная работа магнитометра и магнитных катушек. Поэтому существуют два варианта: использовать математические модели на наноспутнике в орбитальном полёте (модель МПЗ и модель движения) или попеременно включать магнитометр и катушки. Первый вариант соответствует непрерывному режиму работы алгоритма демпфирования начальных угловых скоростей наноспутника, второй -дискретному.
Под непрерывным режимом работы алгоритма демпфирования начальных угловых скоростей наноспутника понимается режим, при котором ток на магнитные катушки подаётся непрерывно.
Под дискретным режимом работы алгоритма демпфирования начальных угловых скоростей наноспутника понимается режим, при котором ток на магнитные катушки подаётся в течение 16 секунд с интервалом в одну секунду для проведения измерения напряжённости МПЗ [ 1].
На рис. 2 представлено сравнение быстродействия алгоритма в непрерывном и дискретном режимах работы при различных начальных угловых скоростях отделения. Как видно из графиков, выигрыш по времени для непрерывного режима работы несущественный - около 30 %, а для больших начальных угловых скоростей -10 %.
60
О -1-1-1-1-1-1-1-1—
5 10 15 20 25 30 35 40 45 Начальные угловые '/с
Рис. 2. Длительность процесса демпфирования в непрерывном и дискретном режимах работы при различных начальных угловых скоростях отделения наноспутника
Поэтому непрерывный режим работы алгоритма применим для проверки функционирования разработанного про-
граммного обеспечения, а в качестве бортового режима для наноспутника SamSat-218Д выбран дискретный режим работы алгоритма демпфирования угловой скорости.
Исследование выполнено при финансовой поддержке Министерства образования и науки Российской Федерации в рамках договора на проведение НИОКТР от 11.03.2013 г. № 122/12 «Создание высокотехнологичного производства мало-массогабаритных космических платформ, информационных технологий и программно-аппаратного комплекса приёма и обработки гиперспектральных данных».
Библиографический список
1. Мельник М.Е. Моделирование процесса демпфирования угловых скоростей наноспутника SamSat-218Д после выхода из транспортно-пускового контейнера // Вестник Самарского государственного аэрокосмического университета имени академика С. П. Королёва (национального исследовательского университета). 2016. Т. 15, № 1. С. 89-96. DOI: 10.18287/2412-7329-2016-15-1-89-96
2. Беляев М.Ю. Научные эксперименты на космических кораблях и орбитальных станциях. М.: Машиностроение, 1984. 264 с.
3. Belokonov I.V., Kramlikh A.V., Timbai I.A. Low-orbital transformable nanosatellite: research of the dynamics and possibilities of navigational and communication problems solving for passive aerodynamic stabilization // Proceedings of 2th IAA Conference on Dynamics and Control of Space System (DYCOSS). 2014. V. 153. 15 p.
4. Whitford Ch., Forrest D. The CATSAT Attitude Control System // 12th AIAA/USU Conference on Small Satellites.
http://digitalcommons.usu.edu/cgi/viewcontent.cgi?article=2252&context=smallsat.
5. Сёмкин Н.Д., Любимов В.В., Малышев В.И. Моделирование законов функционирования магнитных исполнительных органов при ориентации микроспутника по местной вертикали // Физика волновых процессов и радиотехнические системы. 2012. Т. 15, № 1. С. 103-108.
Информация об авторах
Крамлих Андрей Васильевич, кандидат технических наук, доцент, доцент кафедры космического машиностроения, Самарский национальный исследовательский университет имени академика С.П. Королёва. E-mail: [email protected]. Область научных интересов: навигация, ориентация и управление космическими аппаратами, ком-плексирование и обработка информации.
Мельник Мария Евгеньевна, инженер-конструктор второй категории АО «РКЦ «Прогресс». E-mail: [email protected]. Область научных интересов: алгоритмы управления угловым движением космического аппарата, комплексирование и обработка информации.
ON-BOARD ALGORITHM FOR SAMSAT-218D NANOSATELLITE ORIENTATION AND STABILIZATION SYSTEM
© 2016 A. V. Kramlikh1, M. E. Melnik2
1Samara National Research University, Samara, Russian Federation 2JSC Space Rocket Center «Progress», Samara, Russia Federation
Algorithms for determining the attitude and angular rate damping of the SamSat-218D nanosatellite developed at Samara State Aerospace University are described in the paper. The date of the prospective launching of the satellite into the Earth orbit from Vostochniy spaceport is April, 2016. A triaxial magnetometer and luminance sensors are used as sources of information in the orientation system. Three magnetic coils serve as actuators providing the angular rate damping onboard the nanosatellite. A known method of vector coordination with the original criterion of preliminary measurement rejection is used as the orientation determination algorithm since the vector coordination method is very sensitive to measurement errors. The traditional B-dot algorithm is used as the damping algorithm. The influence of the initial angular rates on the magnetic moment coefficient in the B-dot control algorithm is studied. Since the continuous mode of operation of magnetic coils on board the nanosatellite is unpractical due to the necessity of measuring variations in the Earth magnetic field intensity vector the damping process is analyzed in the discrete mode of operation at different values of angular rates of nanosatelite separation. The comparison of the algorithm processing speed in the continuous and discrete operating modes at various initial angular rates of separation is given.
Nanosatellite, attitude, angular motion stabilization, triaxial magnetometer, solar sensor, magnetic coil, algorithm.
References
1. Melnik M.E. Simulation of the process of angular rate damping of the SamSat-218D nanosatellite after separation from the launching-transporting container. Vestnik of the Samara State Aerospace University. 2016. V. 15, no 1. P. 89-96. (In Russ.) DOI: 10.18287/24127329-2016-15-1-89-96
2. Belyaev M.Yu. Nauchnye eksperimenty na kosmicheskikh korablyakh i orbital'nykh stantsiyakh [Scientific experiments on spacecraft and orbital stations]. Moscow: Mashi-nostroenie Publ., 1984. 264 p.
3. Belokonov I.V., Kramlikh A.V., Timbai I.A. Low-orbital transformable nanosatellite: research of the dynamics and possibilities of navigational and communication problems solving for passive aerodynamic stabilization. Proceedings of 2th IAA Conference on Dynamics and Control of Space System (DYCOSS). 2014. V. 153. 15 p.
4. Whitford C., Forrest D. The CATSAT Attitude Control System. 12th AIAA/USU Conference on Small Satellites. Available at:
http://digitalcommons.usu.edu/cgi/viewcontent.cgi?article=2252&context=smallsat.
5. Semkin N.D., Lyubimov V.V., Malyshev V.I. Modeling the operation laws of magnetic control devices with orientation a microsatellite along the local vertical. Physics of Wave Processes and Radio Systems. 2012. V. 15, no. 1. P. 103-108. (In Russ.)
Citation: Kramlikh A.V., Melnik M.E. On-board algorithm for SamSat-218D nanosatellite orientation and stabilization system. Vestnik of the Samara State Aerospace University. 2016. V. 15, no. 2. P. 50-56. DOI: 10.18287/2412-7329-2016-15-2-50-56
About the authors
Kramlikh Andrey Vasilyevich, Candidate of Science (Engineering), Associate Professor, Department of Space Mechanical Engineering, Samara National Research University, Samara, Russian Federation. E-mail: [email protected]. Area of Research: spacecraft navigation, orientation and control, information integrating and processing.
Melnik Maria Evgenievna, design engineer, JSC Space Rocket Center «Progress», Samara, Russian Federation. E-mail: [email protected]. Area of Research: spacecraft angular motion control algorithms, information integration and processing.
