CC BY
132
УДК 531.36 + 629.78
ПОВЫШЕНИЕ ТОЧНОСТИ ОРИЕНТАЦИИ КОСМИЧЕСКИХ АППАРАТОВ МОНИТОРИНГА ЗЕМЛИ В ПОЛЁТЕ
© 2012 Р. Н. Ахметов, А. В. Соллогуб, В. П. Макаров, А. А. Головченко, А. В. Филатов
Г осударственный научно-производственный ракетно-космический центр «ЦСКБ-Прогресс», г. Самара
Рассмотрены вопросы реализации полётной калибровки бесплатформенной инерциальной системы ориентации в части оценки взаимного положения осей измерителя угловой скорости и астроизмери-тельного базиса системы управления ориентацией, а также вопросы определения и учёта погрешностей измерителя угловой скорости. Предложен алгоритм полётной геометрической и параметрической калибровки этого измерителя.
Измеритель угловой скорости, геометрическая и параметрическая полётная калибровка, взаимное положение осей, астроизмерительный базис, бортовая аппаратура, система управления ориентацией.
Введение
При управлении автоматическими космическими аппаратами мониторинга Земли (КА МЗ) одной из центральных является задача определения углового положения в пространстве (ориентации) осей КА МЗ в обеспечение требуемой высокой точности наведения целевой аппаратуры (ЦА) на наблюдаемые наземные объекты. На многих современных КА эта задача решается бесплатформенной инерциальной системой ориентации (БИСО). Принцип работы БИСО заключается в формировании данных об угловом положении КА по информации измерителя угловой скорости (ИУС) об абсолютной угловой скорости КА и периодической её коррекции по информации от астродатчиков (АД). Погрешности измерений, которые даёт бортовая аппаратура (БА), составляют единицы угловых секунд, а погрешности ориентации КА МЗ - на порядок хуже вследствие рассогласования осей БА в полёте и её систематических составляющих погрешностей. Это обусловлено как погрешностями изготовления конструкции КА и взаимной установки БА (угловые минуты), так и эксплуатационными условиями на участке выведения КА (удары и вибронагрузки от раке-
ты) и периодическими (витковыми) изменениями температурных градиентов в полёте. Поэтому повышение точности ориентации КА связывается не только с созданием более совершенных его конструкций, но и с разработкой способов (прямых и алгоритмических) определения в полёте текущего взаимного положения осей БА, систематических составляющих погрешностей и их учёта в процессе управления.
Наиболее известные подходы к решению задачи калибровки базируются на проведении серии плоских калибровочных разворотов КА по каждому каналу на углы от 180° до 360°, получении измерений угловой скорости по информации ИУС и АД, оценивании их невязок с последующим вычислением погрешностей
[1]. Их общим недостатком является большая длительность процесса, что не только снижает производительность КА, но и обуславливает зависимость результатов калибровки от температурной нестабильности конструкции КА.
В настоящей работе предлагается алгоритм полётной геометрической и параметрической калибровки центрального элемента современных БИСО - измерителя угловой скорости, отличающийся малой длительностью процесса калибровки, достигаемой вследствие уменьшения не-
обходимых углов поворотов (до ±20 гра- ляет совмещать калибровку с целевой ра-дусов), что, в свою очередь, также позво- ботой КА.
АД-№1г
АД-№21
АД-№ЗЯ
ИУС-ОЫг і,\азпІ лек
ИУС-Ж?*"
Алгоритм4" съёма-информации ■*! с-АД*~
Алгоритм*" съёма ■ информации-*" с-ИУС1"
::ск
сск
Алгоритм-
определения1
ориенгаиии-
асгроизме-*"
ригельного-
оазиса*"
лііск
Алгоритм*"
определения1
ориенгаиии-
КА1"
Рис. 1. Структурная схема БИСО
Общая схема бесплатформенной системы ориентации
Вариант структурной схемы БИСО КА типа «Ресурс-ДК» представлен на рис. 1.
АД обеспечивает прецизионное определение положения осей 0ХУ2 своей приборной системы координат (ПСКАд) в инерциальной системе координат (ИСК), которая физически реализуется бортовым каталогом звёзд (прямое восхождение и склонение звёзд) на эпоху 7"а. Внутренняя система координат АД (ВнСКАд) определяется положением плоскости фотоприёмника изображений звёзд. Угловые расстояния между направлениями на центры зафиксированных изображений звёзд сравниваются с углами между направлениями на звёзды по бортовому каталогу, что приводит к их опознаванию (с погрешностью около 1 угл. с). В результате становится известным кватернион ориентации осей ВнСКАд относительно ИСК. Приборная система координат АД фиксируется посадочной площадкой прибора.
ИУС предназначен для измерения абсолютной угловой скорости КА при его угловых движениях. Конструктивно ИУС представляет собой моноблок, жёстко закреплённый на корпусе КА и включающий в себя три и более чувствительных элемента (ЧЭ) - одноосных измерителей
угловой скорости. Наиболее простым является ИУС с тремя ЧЭ, оси чувствительности которых ортогональны и совпадают с осями приборной системы координат (ПСК) ИУС.
Алгоритм съёма информации с АД обеспечивает получение данных от АД в соответствии с протоколом информационного обмена.
Алгоритм съёма информации с ИУС обеспечивает получение данных (угловой скорости) от ИУС в соответствии с протоколом информационного обмена и проецирование её на оси связанной системы координат (ССК) КА.
Алгоритм определения ориентации астроизмерительного базиса (АИБ) - ортогональной системы координат, построенной по информации «виртуального» АД
[2], обеспечивает формирование информации об его угловом положении в ИСК (Л^ ) по информации об угловом положении ПСК АД в ИСК ( Л ПСКАД ).
Алгоритм определения ориентации КА обеспечивает формирование информации об угловом положении КА в ИСК с частотой работы системы определения ориентации. Текущее угловое положение КА ( Л 4) определяется выражением:
(1)
где Л t и Л t _At - кватернионы углового положения КА на момент окончания текущего и предыдущего периода работы системы;
АЛи - приращение кватерниона, определяющее изменение углового положения КА за текущий период работы системы At, по измерениям ИУС.
Приращение кватерниона (АЛи ) определяется по информации о проекциях средней за период работы системы абсолютной угловой скорости вращения:
М ‘At
АЯЛ = cos- *
(2)
2
нх . |н| М,
ДЯ. = Т^781П —-’
Н 2
„ ну Н ■ Д1
ДЯ = 7-^1^— ---’
Н 2
ДЯ н2 . |н| Д^
ДЯ, = г-^81П —-’
Н 2
где ДЯ0, ДЯ., ДЯ,, ДЯ, - элементы кватерниона ДЛ;
Нх, Ну, Нъ - проекции вектора средней за период работы системы угловой скорости на оси ПСК ИУС.
Для расчёта кватерниона (1), определяющего текущее угловое положение КА, необходимо знание его начального значения.
Начальному значению кватерниона Л {_Д присваивается значение кватерниона
ориентации (ЛИ™), полученное алгоритмом построения астроизмерительного базиса по информации с астродатчиков.
Следует отметить, что АД определяют ориентацию своих приборных осей относительно ИСК непосредственно, а ИУС - косвенно, путём интегрирования измеряемых значений угловых скоростей. При астрокоррекции в качестве эталонных значений углового положения КА используются данные с АИБ.
Наибольшие точности определения положения КА могут достигаться при использовании информации, получаемой с двух синхронно работающих АД, установленных на КА под разными углами. Но временной диапазон, при котором возможна синхронная работа двух АД,
ограничен из-за засветки их небесными телами - Солнцем, Луной, Землёй. Причём в процессе полёта КА число участков, на которых возможна синхронная работа с двумя АД, значительно меньше числа участков, на которых возможна работа с одним из АД. В этой связи возникла идея построения и использования «виртуального» АД (его математической модели), с помощью которого обеспечивается получение точности ориентации КА, близкой к точности, достигаемой при двух синхронно работающих АД, по аналогии с [2]. В этом случае, определив невязки между показаниями реальных и «виртуального» АД, можно осуществлять управление КА по измерениям одного АД, корректируя их с учётом невязок.
Повышение точности ориентации КА на основе полётной калибровки БА
Приращение ( ДЛи ) кватерниона, определяющего ориентацию КА по выражениям (2), вычисляется с использованием данных измерений угловой скорости КА прибором ИУС. Поскольку эти данные включают в себя инструментальные ошибки ИУС и погрешности его установки, то погрешность определения кватерниона Л{ возрастает со временем. Для
устранения ошибок определения углового положения КА проводится коррекция кватерниона Л{ с использованием информации, формируемой алгоритмом построения астроизмерительного базиса по измерениям астродатчиков.
Коррекция кватерниона Л{ проводится каждый раз при наличии информации с астродатчиков (режим «коррекция»). При этом информация с АД поступает в систему управления с максимальной частотой 4 Гц, если в их полях «зрения» отсутствуют Солнце, Луна или Земля. При нахождении в полях «зрения» АД Солнца, Луны или Земли информация об астроизмерениях не поступает и кватернион (ЛД определяющий ориентацию КА, не корректируется, а система управления работает в режиме «память», то
есть кватернион, определяющий ориентацию КА, формируется только по информации с измерителя угловой скорости.
Коррекция кватерниона Л; производится в соответствии с выражением:
Л; = Л ;°ДЛкор, (3)
где
ДЛкор (Я0кор, ягр ),
ДЯскор=41 _Х .(ДЯГ )2,
ДЯкор = т ■ДЯ1,1 = 1,3, т - коэффициент фильтрации (т < 1),
ДЛ =Л;, °ЛА,
ЛА = ЛАИБ -'Ч ИСК’
Л'; - кватернион откорректированного текущего углового положения КА.
В настоящее время разработаны алгоритмы калибровки, совмещающие в себе геометрическую и параметрическую калибровку. Под геометрической калибровкой ИУС понимается определение погрешности углового положения осей чувствительности ИУС относительно астро-измерительного базиса.
Под параметрической калибровкой ИУС понимается определение компонентов погрешности измерений ИУС, а именно:
- погрешностей масштабных коэффициентов;
- погрешностей постоянных составляющих ошибки измерения угловой скорости.
Схема взаимного положения систем координат (и соответствующие кватернионы), используемых в БИСО, отражающих связь КА (ССК), целевой аппаратуры (ВСК - визирная система координат), АИБ, ИУС и их компонентов (ВОГ -волоконно-оптический гироскоп) с ИСК, приведена на рис. 2.
Погрешность взаимного положения
Л АИБ
пЛ’-гт
) обусловлена погреш-
АИБ и ССК ностями:
п А- установки АД (° АД );
установки фотоприёмника
уст
АД
* аУСТ (оЛ*»);
измерений АД ( тЩ );
- установки ИУС (^Чіус);
гь луст
- установки ВОГ ИУС (°*^бог);
- измерений ИУС (^шиус).
Рис. 2. Взаимное положение систем координат, используемых в БИСО
Описание алгоритма калибровки ИУС
Задачей процесса калибровки ИУС является определение (рис. 3):
- погрешностей углового положения одноосных измерителей ИУС относительно астроизмерительного базиса (^н;,
1=1,2,3);
- погрешностей масштабных коэффициентов одноосных измерителей ИУС (Ж^, ¡=1,2,3);
- погрешностей постоянных составляющих ошибок измерения одноосных
измерителей ИУС (00іу*- і , ¡=1,2,3).
оч
Рис. 3. Положение осей чувствительности (04) ПУС в ССК На рис. 3 Д - углы между измерительными осями одноосных измерителей и их проекциями на плоскость 2сскXсск, которые отсчитываются от проекций в сторону + Усск;
_ углы между проекциями измерительных осей одноосных измерителей на
плоскость 2сскXсск и осью 2сск, которые отсчитываются от оси + 2сск в сторону
+Хсск •
Для выполнения калибровки ИУС КА выбираются четыре участка (/=1,2,3,4) длительностью Т каждый, где производятся вращения КА со скоростями: юх,72=0; ю#0, «72=0; юуфО, «х,2= 0; «2^0, юх,7=0.
На момент начала и окончания каждого участка запоминается информация с ИУС и АД:
- !н1 (1:к1) - время привязки информации на момент начала (окончания) участка;
- ЛАИб (ЛАИ ) - кватернион ориентации АИБ в ИСК (Л )на момент начала
(окончания) участка;
т ОЧ н1 / т ОЧ к1 \
- !^иус (!®иУс ) - интегралы проекций абсолютной угловой скорости вращения КА на оси чувствительности ИУС на момент начала (окончания) участка.
По завершении работы КА на всех четырёх участках (/ = 1,2,3,4) с учётом трёх осей чувствительности (ОЧ) ИУС (/' = 1,2,3) формируется система 12 алгебраических уравнений с 12 неизвестными для определения погрешностей масштабных
коэффициентов (Якг), уходов ; ) и
углового положения осей ИУС относительно базовых осей АД ,):
6к: ЮАД;/+ бСи^ е + ЮАД(г+1)/ 51П(^«Г. ) -
- ШАД(,+2)/51п(ад.)=ШАД,/-ШИ''С,/;
где г - «уход» ИУС по оси /; соЛД//, ЮАД(/+1)/, ЮАД(/+2>; - проекции вектора угловой скорости КА на /-м участке, определённого по измерениям АД (ЛАИщ , ЛА1ИБ );
юиусг/ - проекции того же вектора, определённого по измерениям ИУС
/т^ОЧ н1 т_,ОЧ к1 \
(1®ИУС , !®ИУС ).
Решение системы не представляет больших трудностей, так как матрица системы является разреженной. Например, я оч
параметр £и>'::г может быть вычислен по
окончании самого первого участка, для которого Юххг = 0, так как в данном случае разность показаний по АД и ИУС обусловлена только уходами ИУС по осям. В связи с этим приходим к системе из 9 уравнений с 9 неизвестными.
После завершения расчётов полученные результаты калибровки передаются в ИУС для использования при формировании выходной информации.
Алгоритм проведения калибровки ИУС представлен на рис. 4.
В первом блоке осуществляется съём с ИУС измеренной информации об угловой скорости вращения КА (®Пск) и
её проецирование на оси ССК («и) с использованием матрицы направляющих косинусов номинального положения осей
ПСК ИУС в ССК (мПСКиус). В случае использования в приборе результатов калибровки (аклб = 1), проецирования его выходной информации на оси ССК не проводится.
Во втором блоке осуществляется вычисление за период работы алгоритма Д! приращения кватерниона ориентации КА (ДЛи )и расчёт текущего неоткоррек-тированного значения кватерниона ориентации КА ( Л!).
В третьем блоке при наличии информации с АД (ЛА™ ) определяется погрешность (ДЛкор) кватерниона Л! на текущем интервале Дt и выполняется его коррекция.
Четвёртый блок является подготовительным для калибровки ИУС. Здесь запоминаются исходные данные для расчёта результатов калибровки на момент начала каждого участка (Л^ , 1®Иж_н1, 1н1,
1 = 14).
В пятом блоке запоминаются исходные данные для расчёта результатов калибровки на момент окончания каждого участка калибровки ( ЛАИБ , 1«ИУск1, !к1,
1 = 14).
С Вход 'У
— изм лАИБ “пск^ “иск
_ длПСКиуС нЗМ ш _ тсск опск
И — „НЗМ
(тг = т
гтск
ДЛ^ДЛ^Дл^ Дл“ =Л-.3т-^^г Л/..0 = 1 - ¡(Л^)3. А, = Л^^ЛЛ".
нет
Л$ = Л^МДЙ
ЛЛ = Л^Л?, ЛЛкор(/™!,рД?‘Л>-Г! = т-л^, д; ™р = 1 - ¿л>™?-)2. л; = Л,аДЛкор
К Л Б
Обход Н У
Задаййе
начальных
условий
А1ГА — дАИБ _ 1
“аив ~ Фискэ 1
1, _ I ОЧлзм Гн1 _ *.
*°Н\УС 16)ИУС ■ с жг
Накопление информации
о погрешностях Парасчет
ИУСна/<Т да Т погрешностей ИУС
ДК1 — АЛ
ал/№ — /1КЩ’
,ОЧ_*( _ т ОЧлэы *и1 _ >
1(*>ИУС. УС ' С — £
1=1 + 1, и^^Х + Т
нст
i !а задание результатов калибровки в ИУС
,Д4
,04 „
г
Л.УС
;
Рг/с. 4. Алгоритм калибровки ПУС, совмещённый с острокоррекцией
Собственно калибровка осуществляется в шестом блоке, где по накопленным данным на всех участках (/=1,...,4) для каждой оси чувствительности ( = 1,2,3) формируется и решается система 12 алгебраических уравнений с 12 неизвестными.
Решения полученных уравнений
(8а;, 8Д, 8к1, ¡=1,2,3) передаются в
ИУС для их использования при формировании его выходной информации.
Результаты моделирования процесса калибровки ИУС
Точность результатов калибровки зависит от длительности измерительных участков и от погрешностей выходной информации ИУС и АД.
Длительности участков измерения выбраны равными 240 с. Погрешности измерений АД и ИУС, а также значения оцениваемых в процессе калибровки параметров приведены в таблицах 1, 2, 3, 4. Влияние инструментальных погрешностей БА выражается в увеличении погрешности результатов калибровки, но точность калибровки остаётся достаточной для определения оцениваемых величин.
Таблица 1. Случайная (О) и предельная суммарная погрешности измерения углового положения осей приборной системы координат АД в ИСК
Параметр Угловая скорость, ° /с Предельная погрешность, угл. с Случайная составляющая (О ), угл. с
Направление оси 1 ПСК 0,0 7 0,5
0,067 9 1,0
1,2 25 5
Направление осей X, У ПСК 0,0 41 2
0,067 50 5
1,2 95 20
Таблица 2. Погрешности измерения интегралов проекций абсолютной угловой скорости прибором ИУС
Параметры Значение, угл. с
Погрешность приращения интегралов проекций угловой скорости за любые:
5 с 0,6
10 мин. 7,5
20 мин. 15
40 мин. 30
24 часа 300
Цена младшего разряда выдаваемой информации по приращению
интегралов 1-10"3
Таблица 3. Предельные значения оцениваемых параметров
8< , 0/с 8а , угл.мин. 80 , угл.мин. 8к, безразмерная
2,8 10-5 1,5 1,5 0,5 10-4
Таблица 4. Погрешность результатов калибровки
Оцениваемый параметр /с,0 о £ «1 8а , угл.мин. 80, угл.мин. 8к, безразмерная
Значение параметра 2,8 10-5 1,5 1,5 110-4
Ошибка оценки с учётом погрешностей БА 3,96 -10-6 0,184 0,09 1,85 -10"5
В процессе штатной эксплуатации КА МЗ типа «Ресурс-ДК1» успешно отработан первый участок калибровки, оценивающий погрешности постоянных составляющих ошибок измерения ИУС (доО4 )• Результаты оценки постоянных составляющих ошибок измерения ИУС (<5шОЧ) представлены на рис. 5.
8о^. 7с
-0.0001
-0.0002
-0.0003
<Ч\
Рис. 5. Результаты оценки £^°Ч изд. «Ресурс-
ух
ДК1». N - порядковый номер участка оценки ухода
Заключение
Предложенный алгоритм проведения геометрической и параметрической калибровки центрального звена бескар-данной инерциальной системы ориента-
ции - измерителя угловой скорости - позволяет не только достичь необходимые для работы КА МЗ точности ориентации, но и уменьшить время проведения калибровки за счёт совмещения процесса калибровки со штатной работой КА.
Библиографический список
1. Кирилин, А. Н. Методы обеспечения живучести низкоорбитальных автоматических КА зондирования Земли: математические модели, компьютерные технологии [Текст] / А. Н. Кирилин, Р. Н. Ахметов, А. В. Соллогуб, В. П. Макаров.
- М.: Машиностроение, 2010. - 384 с.
2. Аншаков, Г. П. Квазинепрерывная астрокоррекция инерциальных датчиков памяти спутников дистанционного зонди-р ования Земли [Текст] / Г. П. Аншаков, Р. Г. Залялова, Я. А. Мостовой, В. А. Ти-пухов // XV Международная конференция по интегрированным навигационным системам. - СПб.: Изд. ЦНИИ «Электроприбор», 2008. - С. 304-311.
IMPROVEMENT OF THE FLIGHT ATTITUDE CONTROL ACCURACY OF THE EARTH OBSERVATION SPACECRAFT
© 2012 R. N. Akhmetov, A. V. Sologub, V. P. Makarov, A. A. Golovchenko, A. V. Filatov State Research and Production Space Centre “TsSKB-Progress”, Samara
The paper presents the problems of assessing the relative position of angular accelerometer sensitivity axes and astromeasuring basis of the attitude control system as part of flight calibration of a platform-independent inertial system. The inaccuracy of the angular accelerometer is defined and taken into account. An algorithm of in-flight geometrical and parametric calibration of the meter is suggested.
Angular accelerometer, geometrical and parametric calibration, relative position of axes, astromeasuring basis ,onboard equipment, attitude control system.
Информация об авторах Ахметов Равиль Нургалиевич, доктор технических наук, генеральный конструктор, ФГУП ГНПРКЦ «ЦСКБ - Прогресс». E-mail: [email protected]. Область научных интересов: проектирование и испытания ракетно-космической техники.
Соллогуб Анатолий Владимирович, доктор технических наук, главный научный сотрудник, ФГУП ГНПРКЦ «ЦСКБ - Прогресс». E-mail: [email protected]. Область научных интересов: моделирование ракетно-космических систем.
Макаров Валентин Павлович, доктор технических наук, научный советник, ФГУП ГНПРКЦ «ЦСКБ - Прогресс». E-mail: [email protected]. Область научных интересов: системы управления движением космических аппаратов.
Головченко Анатолий Алексеевич, ведущий инженер-конструктор, ФГУП ГНПРКЦ «ЦСКБ - Прогресс». E-mail: [email protected]. Область научных интересов: измерительные средства и системы управления ориентацией космических аппаратов.
Филатов Артём Викторович, ведущий инженер-конструктор, ФГУП ГНПРКЦ «ЦСКБ - Прогресс». E-mail: [email protected]. Область научных интересов: измерительные средства и системы управления ориентацией космических аппаратов.
Akhmetov Ravil Nurgalievich, doctor of engineering science, General Designer, Space Rocket Center “TsSKB-Progress”. E-mail:[email protected]. Area of research: designing and testing of space rocket systems.
Makarov Valentin Pavlovich, doctor of engineering science, scientific adviser, Space Rocket Center “TsSKB-Progress”. E-mail: [email protected] of research: spacecraft motion control systems.
Sollogub Anatoliy Vladimirovich, doctor of engineering science, Chief scientific adviser, Space Rocket Center “TsSKB-Progress”. E-mail: [email protected] of research: modeling of space rocket systems.
Golovchenko Anatoly Alexeevich, leading design engineer, Space Rocket Center “TsSKB-Progress”. E-mail: [email protected]. Area of research: measuring instruments and spacecraft attitude control systems.
Filatov Artyom Victorovich, leading design engineer, Space Rocket Center “TsSKB-Progress”. E-mail: [email protected]. Area of research: measuring instruments and spacecraft attitude control systems.
