УДК 629.78 : 681.51
КОМПЬЮТЕРНАЯ ИМИТАЦИЯ И АНИМАЦИЯ УПРАВЛЯЕМОГО ДВИЖЕНИЯ ИНФОРМАЦИОННОГО СПУТНИКА
© 2018 Т.Е. Сомова
Самарский государственный технический университет
Статья поступила в редакцию 22.02.2018
Кратко представляются результаты применения разработанной программной среды для имитации и анимации пространственного движения информационных спутников, которая входит в состав системы автоматизированного проектирования систем управления космических аппаратов. Ключевые слова: космический аппарат, управляемое движение, имитация, анимация.
Работа поддержана РФФИ (грант 17-48-630637) и отделением ЭММПУРАН (программа фундаментальных исследований № 13).
ВВЕДЕНИЕ
При проектировании космических информационных систем применяются современные компьютерные технологии. Программная система SIRIUS-S разработана в отделе «Навигации, наведения и управления движением» НИИ ПНМС СамГТУ и предназначена для эскизного проектирования систем наведения, навигации и управления движением информационных спутников [1]. Данная программная система может применяться для анализа различных вариантов служебных систем космического аппарата (КА) и формирования требований к измерительным и исполнительным подсистемам бортового комплекса управления.
Программная система SIRIUS-S содержит диалоговый монитор, подсистемы моделирования, синтеза и анализа, а также технологические подсистемы анимации пространственного движения КА и документирования результатов. Управление работой SIRIUS-S выполняется с помощью системы иерархических меню. В результате инженер-разработчик получает функциональный облик системы управления проектируемого спутника радиолокационного или оптико-электронного наблюдения в отношении периодичности, производительности и оперативности, разрешения на местности, точности наведения и стабилизации КА с учетом возмущений, ограничений и др. факторов.
Отметим, что одним из наиболее наукоемких направлений развития космических систем в радиолокационном спектральном диапазоне является применение технологии синтеза апертуры антенны, основанной на приеме сигналов
Татьяна Евгеньевна Сомова, аспирантка, младший научный сотрудник отдела навигации, наведения и управления движением, НИИ Проблем надежности механических систем СамП'У. E-mail te_somova@mail.ni
при длительном временном интервале полета носителя радиолокатора. При этом происходит искусственное увеличение (синтез) апертуры антенны за счет движения радиолокатора и последующей цифровой обработки измерений. Радиолокаторы с синтезированной апертурой (РСА) космического базирования в настоящее время достигли высокого технического совершенства и являются важным средством мониторинга поверхности Земли. За счет особенностей электромагнитного излучения в радиодиапазоне такие системы в меньшей степени зависят от метеорологических условий и орбитального расположения, чем их аналоги в видимом и инфракрасном диапазонах. Применение РСА с активной фазированной антенной решеткой (АФАР) существенно упрощает реализацию различных режимов радиоэлектронной съемки.
В системе 8ШШ5-Б реализованы модели Земли, конструкции КА, его поступательного и углового движения как при различных видах сканирующей съемки, так и при поворотных маневрах. Эти модели позволяют выполнить точный расчет всех кинематических параметров пространственного движения КА и сформировать набор хш-файлов для подсистемы визуализации. Документируемые результаты представляются сценами маршрутов съемки на картографической основе, таблицами и графиками изменения координат движений КА и исполнительных органов в функции времени, значениями достигнутых характеристик качества наблюдения и критериев оптимальности.
Подсистема моделирования содержит следующие компоненты: модель Земли и электронные карты с базами данных для объектов на земной поверхности; структурная модель КА - геометрические и инерционные характеристики конструкции, параметры телескопов и антенн, исполнительных органов, измерительных
а
Рис. 1. Геостационарные спутники связи: а) Луч-5А/В; Ъ) Экспресс-АМ5/6
Рис. 2. Спутники оптико-электронного землеобзора: а) Ресурс-Л; Ъ) Сириус-1
подсистем и др.; баллистическая модель движения центра масс КА; модели углового движения КА - схемы обзора земной поверхности, методы расчета маршрутных движений и поворотных маневров (ПМ) с краевыми условиями общего вида и др.; модели внешней обстановки - методы расчета освещенности наблюдаемой поверхности, облачности и др.; модели приемного тракта - методы расчета скоростей движения изображения на матрицах оптико-электрон-ных преобразователей (ОЭП) при оптической сканирующей съемке, разрешающей способности и др.; модели целевой обстановки - базы данных задач и объектов наблюдения.
Подсистема анализа и синтеза предназначена для выполнения таких функций: отображение карт земной поверхности с объектами наблюдения; отображение орбиты и трассы полета КА, полосы обзора; расчет орбитального поступательного и углового движения КА, проверка их реализуемости при ограниченных ресурсах исполнительных органов; синтез алгоритмов определения ориентации, законов наведения, цифрового и широтно-импульсного управления КА; анализ устойчивости и качества переходных процессов в системе управления ориентацией (СУО) спутника при выполнении требуе-
мого движения; расчет энергобаланса системы электропитания, освещенности наблюдаемой местности, радио-видимости пунктов приема информации и др.; расчет пространственного разрешения на местности; анализ вариантов построения СУО КА по различным критериям; расчет эффективности функционирования информационного спутника по критериям периодичности, производительности и оперативности наблюдений на основе компьютерной имитации. В этой подсистеме алгоритмы имитации реализованы в среде MatLab.
Подсистема анимации пространственного движения КА является технологическим программным средством, созданным в среде Delphi 7 с применением графической библиотеки OpenGL. 3D-модель конструкции КА реализована в среде Blende, отображение ее элементов с учетом освещённости Солнцем выполняется средствами OpenGL. Здесь с использованием стандартной процедуры «наклейки» текстуры в виде карты Земли на её поверхность отмечаются объекты наблюдения, их геодезические координаты рассчитываются в системе координат (CK) текстуры и определяется положение КА в CK сцены съемочной камеры. Отображаются конструкция КА с текущей точкой трассы по-
лета, точкой пересечения линии визирования с поверхностью Земли и проекцией линейки ОЭП на эту поверхности, если в этот момент выполняется съемка. Программа визуализации позволяет изменять масштаб изображения и ракурс наблюдения сцены при сканирующей съемке заданного участка на поверхности вращающейся Земли.
Программная система SIRIUS-S успешно применялась при автоматизированном проектировании многих информационных спутников, например геостационарных спутников связи Sesat, Луч-5А/В, Экспресс-АМ5/6 разработки АО «Информационные спутниковые системы им. акад. М.Ф. Решетнева» (рис. 1), спутника оптико-электронного землеобзора Ресурс-П разработки АО РКЦ «Прогресс», а также перспективного спутника Сириус-1 для оптико-электронного наблюдения, рис. 2.
В последнее время проявился существенный интерес предприятий космической промышленности к применению программной системы SIRIUS-S для полетного сопровождения (поддержки) систем управления движением информационных спутников. Современное представление фактической ориентации КА с помощью кватернионов создает проблему при пространственном восприятии поступающей телеметрической информации для верификации бортового расчета соответствующих законов наведения КА. Вторая проблема состоит в том, что при возникновении отказов бортовой аппаратуры управления возникает существенное отклонение фактического углового движения КА от программного движения. При этом, несмотря на доступность телеметрической информации от работоспособных измерительных приборов о фактической ориентации спутника, возникает задача квалифицированной интерпретации аномальной ориентации КА относительно направлений на Землю, на Солнце и на другие внешние ориентиры. Поэтому наличие в центре управления полетами компьютерной среды для анимации пространственного движения КА по телеметрическим данным весьма желательно.
В статье кратко представляется компьютерная среда для анимации движения информационных спутников, работа которой демонстрируется для последовательности сканирующих маршрутов, включающей маршруты стереоскопической съемки заданных участков земной поверхности.
ПРОГРАММНАЯ СРЕДА АНИМАЦИИ
На рис. 3 показан пример визуализации конструкции спутника землеобзора Sirius-1 в среде Blender. Модель конструкции этого КА представлена системой твердых тел и содержит корпус
спутника, телескоп с подвижной крышкой объектива, кластер четырех реактивных двигате-лей-маховиков (ДМ), четыре звездных датчика системы определения углового положения (СОУП) в составе бесплатформенной инерци-альной навигационной системы (БИНС), двигательную установку на основе 8 реактивных двигателей, бортовые антенны, принимающие сигналы навигационных систем FJIOHACC/GPS, и четыре подвижные панели солнечных батарей (СБ), каждая панель в двухстепенном кардано-вом подвесе по схеме «вилка». Предусмотрена возможность отображения перемещений элементов конструкции спутника, что достигается за счет изменения их угловых координат по степеням подвижности. Отображение пространственного положения 31)-модели спутника с учетом освещённости Солнцем выполняется средствами OpenGL. При изменении положения центра масс и ориентации КА изображение его конструкции в текущей сцене «перерисовывается», что создает эффект анимации пространственного движения спутника относительно вращающейся Земли с отображением текущего положения важнейших объектов внешней космической среды - Солнца, Луны и др.
Анимация движения КА реализуется программным средством, созданным в среде Delphi 7 с применением графической библиотеки OpenGL, т.е. на основе свободно распространяемых программ, в отличие от известной системы STK (Analytical Graphics, Inc.. http://www.agi.com). Здесь
Рис. 3. Облик конструкции КА Sirius-1
программно выполняется процедура «наклейки» текстуры в виде карты Земли на её поверхность, отмечаются объекты наблюдения, рассчитываются их географические координаты и далее программно реализуется собственно бВ-анимация движения спутника.
Обработчик события OnCreate для генерации формы содержит следующие строки:
SetDCPixelFormat(Canvas. Handle); // формат пикселя
hrc: =wglCreateContext(Canvas. Handle); / / контекст воспроизведения
Здесь первая строка - обращение к пользовательской процедуре
procedure SetDCPixelFormat (hdc: HDC) ; var pfd: TPixelFormatDesçriptor; nPixelFormat: Integer; begin
FillChar (pfd, SizeOf (pfd), 0); nPixelFormat: = ChoosePixelFormat (hdc, @pfd);
SetPixelFormat (hdc, nPixelFormat, @ pfd) ; end;
Обработка события onpaint выглядит так:
wglMakeCurrent. (Canvas. Handle, hrc) ; // установить контекст
glClearColor (0.5, 0.5, 0.75, 1.0); // цвет фона
glClear (GL_COLOR_BUFFER_BIT); //очистка буфера цвета
wglMakeCurrent (0, 0); // освободить контекст
Необходимо отметить объекты информационного обслуживания на поверхности Земли. В качестве текстуры используется карта Земли (файл earth. bmp). В разделе var задаются переменные tdir:string; tx:pointer;
procedure texture(kl,k2: integer;p:poin ter) ;
begin
glTexParameteri(GL_TEXTURE_2 D, GL_
TEXTURE_MAG_FILTER, GL_linear);
glTexParameteri(GL_TEXTURE_2 D, GL_
TEXTURE_MIN_FILTER, GL_linear);
glTexImage2D ( G L TEXTURE_2D,0,3,kl,k2,0,GL_RGB,GL_ UNSIGNED_BYTE,p);
glTexEnvi(GL_TEXTURE_ENV, GL_TEXTURE_ ENV_MODE, GL_MODULATE);
glEnable(GL_TEXTURE_2D); end;
В процедуре glTexImage2D используются параметры kl, k2 - размер текстуры (т.е. разрешение загружаемой картинки с расширением bmp) ир- переменная, которая подгружает файл earth. bmp. Далее вызывается функция texture, которая зависит от трех параметров: texture(2048,1024,tx); tx:=readbitmap(tdir+'textures \EarthMap_2500xl250_l.jpg' , 2048, 1024); Этой командой описывается двумерная текстура и указывается ее разрешение.
Теперь мы имеем карту земной поверхности и необходимо нанести на нее объекты информационного обслуживания, которые хранятся в массиве. Здесь следует преобразовать географические координаты объектов в координаты среды OpenGL. В рассчитанном массиве мы имеем широту и долготу объектов обслуживания, причем широта ф изменяется в интервале (- к/2,+к/2 ), а долгота А,- в диапазоне (0, 71). Условия преобразования: если X < 71, то без изменений, иначе X := X - 2л. Формулы преобразования х = X/ к иу=ф /(2л)реализованы в упрощенном коде программы
glpointsize (2) ; // размер точек(2) for i:=l to num do begin //цикл рисования glcolor3f (1,0,0); // цвет точек (красный) glbegin(gl_points); if lm[i]<=180 then
glvertex3f(lm[i] /180, 0.5* fi[i] /90, 1. 001)
else givertex3f(-1+(lm[i]/-180)/180,
0.5* fi [i]/90, 1.001) ;
glend;
end;
Представленный код отличается от оригинального кода программы в среде OpenGL. Это сделано для того, чтобы не акцентировать внимание на условия и способ извлечения координат из массивов записей. Между командами glbegin ( gl_points ) и glend располагается код рисования с использованием описанной процедуры пересчета координат, где lm [ i ] и fi [ i ] - массивы значений долготы и широты соответственно, a num - длина массива. Для изменения масштаба используется команда glScalef с 3 аргументами в виде масштабных множителей для каждой оси.
Имитация движения КА землеобзора с цифровой системой управления выполнена в программной системе SIRIUS-S с применением среды Matlab. С использованием среды Blender отображаются конструкция КА, текущая точка трассы КА, точка пересечения линии визирования бортового телескопа с поверхностью Земли и проекция центральной линии матрицы ОЭП на эту поверхность, если в этот момент выполняется сканирующая съемка.
МАТЕМАТИЧЕСКИЕ МОДЕЛИ
Стандартным образом вводятся инерциаль-ная (ИСК) и геодезическая Гринвичская (ГСК) системы координат с началом центре Земли О , орбитальная (ОСК) и связанная с корпусом спутника (ССК) системы координат с началом в его центре масс О. Ориентация ССК в ИСК определяется кватернионом А = (А0Д), где X = \ Л: J, а в ОСК - углами Крылова ф, в последовательности 312.
ьЬ<7) \г .г
4D Л-1 Е'
Ж Y „ v Iг
Г^Л .у dJ2\ У
Ф :
hV а Ъ
Рис. 4. Схема GE (а) и оболочка ее КМ (Ь)
Схема General Electric (GE) кластера ДМ, рис. 4, обладает возможностью управлять ориентацией КА при отказе любого одного маховика. Здесь вССК Оxyz оси вращения четырёх ДМ располагаются на поверхности конуса с углом полу-рас-твора у. Далее используются стандартные обозначения col(-) = И, Нпе(-) = [•], (•)', <,) , [ах] и о5 • для векторов, матриц и кватернионов, Cy=cosy, S sin у, /' = 1,2,3 = 1-^3, применяется вектор модифицированных параметров Родрига (МПР) а = {а.} = е tg(<E>/4) с традиционными обозначениями орта Эйлера е и угла ф собственного поворота, который взаимнооднозначно связан с кватернионом А ориентации КА в ИСК прямыми О = А/(1 + А,0) и обратными =(1-о2)/(1 + а2), X = 2а/(1 + а2) соотношениями. Модель углового движения КА учитывает упругость его конструкции и имеет вид [2] -[4].
А = Л°оэ/2; А°{юДП} = {Fffl,Fq,Fn}, (1) где F® = -[© x]G + Mm+Md;
Fq = -AiJ(Vyq + Wtjq); Fn = M - Mf;
J D, ЛАТ
А0 = D; Aq 0
J, А t У 0 /Л _
A
Ау = S y -S y 0 0
0 0 S y
Здесь С = С +Б а является вектором кинетического момента (КМ) системы, где С = К + Н и К = Л© , столбцы Н = (И ;} и И = 'И ,, = }, р = \~4, представляют КМ кластера и отдельных ДМ, которые связаны соотношением Н = А И , где матрица А, составлена из ортов осей вращения ДМ в ССК; А9 = (Иа^цД; V, =
тчт&аёЩ)2}; Мш-{<}; М = {тр};
Mf={»7^}; вектор механического момента МП Мш = {да™ } = -L х В [6], где вектор электромагнитного момента (ЭММ) L = {/,} с ограниченными компонентами \I:j | < 1m и вектор индукции магнитного поля Земли В = b В с ортом b определены в ССК; векторы-столбцы
M = {тр} и М1 = \ту } представляют управляющие моменты и моменты сил сухого трения
по осям вращения ДМ, а вектор Md - внешние возмущающие моменты. Ресурсы каждого ДМ по управляющему и кинетическому моментам ограничены согласно соотношениям
| «7^(0 |< mm, | |< hm, = 1 ^4.
Вектор Мг управляющего момента кластера ДМ формируется в виде М! = — H , где (•) - символ локальной производной по времени. Если корпус КА считать твердым телом, то
G = G° и модель динамики его углового движения принимает вид
J<D+[a> x]G° =МГ +Mm +Md. (2)
Если задано программное угловое движение КА Л"(/),Ш'(/),0)"(0 = Е"(0 (закон наведения) в ИСК, то погрешность ориентации определяется кватернионом Е = (е0,е) = Ар о Л , вектором Эйлера Е = {е0,е}, где вектор е = {е.}, матрицей Се(Е) = I; 2[ex]Q^, где Qe = I3e0 + [ex] j и вектором 5ф = |0ф J = {2с ej , а погрешность стабилизации угловой скорости
ÔC0 = {8®,} = со - €ею ^(7).
Пусть дискретное измерение кватерниона ориентации A',"=A'"(t,) КА с помощью СОУП выполняется с периодом Т , где
A'" (t, ) = Л(/>) о А;, кватернион At" = Л"(/, ) представляет центрированный гауссовский шум, tM =t, +Тр, 1 g N0 = [0,1,2,3...), а дискретное измерение угловых скоростей вращения ДМ Q, =Q (tj) - в моменты времени | с периодом Tq, ts+1 =ts+Tq, s e N0. Будем считать, что в моменты времени tk с периодом Ти, tk+l = tk + Ти, k g N0, формируется цифро-
вое управление ДМ, а в моменты времени I, с
периодом У;;1 > Ти, , = % + 'Г , Г й -
цифровое управление МП, когда значения компонентов вектора ЭММ Ь = {/,} фиксируются
При балансе СУО по вектору КМ с условием С3 = 0 и отсутствии внешних моментов (Мт =0, Ма=0) модель (2) динамики углового движения КА принимает вид СО = £ с вектором углового ускорения 8 = .Р Мг, а модель углового движения КА - простое кинематическое представление
А = Л о ш/2; ш = 8 ; ё = в* = V. (3)
Модули векторов (0 (/), 8 (7) и 8 (/) ограничены, и (/) <о7, £ (/ ) и б (/) и , что обусловлено ограниченностью областей вариации векторов КМ Н и управляющего момента Мг = — Н кластера ДМ, а также допустимым темпом его изменения. Разработанный [4] аналитический метод синтеза закона наведения КА при пространственном ПМ спутника на интервале времени / е [// . /; | с краевыми ушовиями общего вида
Жр) = Л,;со(/') = ю ;еЩ) = е,;
(4)
Щ) = Лг; га(^) = со(; Цф = ег; е(^) = вг
основывается на необходимом и достаточном условии разрешимости классической задачи Дарбу. Здесь решение представляется как результат сложения трех одновременно происходящих элементарных поворотов «вложенных» базисов вокруг ортов к = 1-^3 осей Эйлера, положение которых определяется условиями (4).
Дтя кластера ДМ принципиальная проблема заключается в распределении векторов его кинетического Н и управляющего М =—Н* моментов между четырьмя ДМ. Дтя решения этой проблемы используется скалярная функция настройки, которая обеспечивает однозначное распределение векторов Н и Мг по явным аналитическим соотношениям.
В СУО применяются дискретная рекуррентная фильтрация сигналов СОУП с периодом Тр , цифровой закон управления кластером ДМ с периодом Ти, компенсационная схема разгрузки от накопленного КМ с помощью цифрового управления МП с периодом Т™ и компенсация влияния моментов сил сухого трения Ш по осям вращения ДМ с использованием дискретного наблюдателя Луенбергера. Здесь выполняется фильтрация значений вектора углового рассогласования €,= — 8фг = — 2е0/е,, /еД-,, с периодом Тр и формируются значения вектора
рассогласования е; . к е М0, которые используются в законе управления кластером ДМ
8м-1 = в8А + ; й, = К«; + Р е I; М1 = а, хс; + .1(С-Е;:' + [С 1(йрКх]щ +щ)и где С:=Се(Е,)
ЗАКОНЫ НАВЕДЕНИЕ ПРИ СТЕРЕОСКОПИЧЕСКОЙ СЪЕМКЕ
На рис. 5 приведена упрощенная схема стереоскопической съемки, которая содержит два маршрута сканирования заданного участка земной поверхности, причем для иллюстрации идеи эта поверхность считается неподвижной плоскостью. Здесь красной прямой линией отмечена орбита спутника с начальными (initial) i и конечными (final) f точками первого ( /' 1, голубой цвет) и второго ( / 2, бордовый цвет) маршрутов, красная пунктирная прямая линия представляет след орбиты. Если для наглядности описания стереосъемки с углом конвергенции 8 использовать стандартные понятия для направлений векторов в ОСК, то каждый маршрут сканирования начинается в фиксированной точке А в ГСК, заданной с превышением относительно поверхности общеземного эллипсоида. При этом в момент времени начала первого маршрута с дальностью D наблюдения точки А орт целеуказания vs получается отклонением от надира сначала по крену и затем по тангажу на угол 5, = 5/2. Первый стерео маршрут съемки должен выполняться с условиями D = const, 5j = const V/ е за
счет углового наведения КА в ОСК по углу тангажа 0, см. рис. 5, голубой цвет. Далее следует поворотный маневр КА на интервале времени
t G [/ta), t[2' ], В момент времени начала второго стерео маршрута съемки с той же дальностью D наблюдения указанной точки А в ГСК орт целеуказания vs формируется отклонени-
i - /(2>
ем от надира, достигнутого при I I , по крену и тангажу на угол ô-, = -5/2. Этот маршрут должен выполняться с условиями D = const, ô2 = const V/ G |X(2),/f2>] также за счет углового наведения КА в ОСК по углу тангажа 9, см. рис. 5, бордовый цвет.
При синтезе закона наведения КА для получения стереопары изображений заданного участка в полной модели земной поверхности на маршрутах сканирования с двумя ракурсами следует учитывать, что невозможно обеспечить строгое постоянство продольной скорости движения изображения (СДИ) в центре матрицы ОЭП при соблюдении основного условия D = const стереосъемки. Здесь приходится изменять текущие значения азимута А маршрута и продольной СДИ в центре матрицы ОЭП. Начало А первого стерео маршрута съемки (MCI) определяется из условия достижения значения угла 5, =5/2 между траверсом точки А на орбитальное положение КА и ортом целеуказания Vs на эту точку при заданном угле конвергенции ô . Это условие определяет также момент времени /¡(1) начала MCI и дальность 1У-]' наблюдения точки А. Выполняется расчет закона углового наведения КА, движущегося по известной орбите, при итерационной вариации значений азимута А маршрута и значения продольной СДИ так, чтобы в результате выполнения такого шага дальность D наблюдения достигнутой точки на земном эллипсоиде была равна /) = D. Далее описанная процедура повторяется с периодом Ts в очередной достигнутой точке MCI и процесс численного построения закона углового наведения КА на первом маршруте съемки завершается по превышению заданной длины MCI на земном эллипсоиде. Момент времени начала второго стерео маршрута съемки (МС2) определяется из условия достижения дальности D наблюдения точки А при угле 5 , = —5/2 между траверсом этой точки на КА и ортом целеуказания v на точку А. Для совпадения MCI и МС2 на земном эллипсоиде шаги построения закона углового наведения КА на втором маршруте выполняются также с периодом Ts , но в классе ортодро-мической съемки с обеспечением постоянства дальности наблюдения D при азимуте и совпадения следа оси визирования телескопа
на земном эллипсоиде с центральной линией MCI [5].
Аналитическое представление законов углового наведения спутника в ИСК при известных численных данных для двух маршрутов сканирующей стереосъемки основывается на интерполяции этих данных векторной функцией МПР а как векторного сплайна седьмого порядка с 8 трехмерными столбцами коэффициентов, которые вычисляются на основе с помощью аналитических соотношений. Краевые условия для поворотного маневра КА между этими маршрутами стереосъемки рассчитываются также на основе явных аналитических соотношений.
Законы наведения КА на солнечно-син-хронной орбите (ССО) высотой 720 км для
маршрутов съемки, в том числе стереосъемки с углом конвергенции ô = 30 град, представлены на рис. 6. Здесь при отчете времени от момента прохождения КА восходящего узла орбиты отмечены
• трассовая сканирующая съемка, маршрут Триполи, /е[539 25,549 25] с;
• ПМ маршрут Триполи => стерео MCI Тунис, / е [549.25,580.75] с;
• стерео MCI Тунис с азимутом А = -13.42 град,точки ¡пи Г , /е[580.75,591.25] с;
• ПМ стерео MCI Тунис => стерео МС2 Тунис,/£[591.25,647.25] с;
• стерео МС2 Тунис с азимутом А = -13.42 град, точки ir, и f12 / £ [647.25,658] с;
• ПМ стерео МС2 Тунис => стерео MCI Флоренция, /£[658,690.75] с;
• стерео MCI Флоренция с азимутом А = -14.18 град, точки i и f2J, /£[690.75,701.5] с;
• ПМ стерео MCI Флоренция => стерео МС2 Флоренция, /£ [701.5,756] с;
• стерео МС2 Флоренция с азимутом А = —14.18 град, точки i22 и f22, /£[756,766.75] с;
• ПМ стерео МС2 Флоренция => маршрут Франкфурт, / £ [766.75,810.75] с;
• съемка с выравниванием, марш-
рут Франкфурт с азимутом А = -45 град, /£[810.75,820.75] с.
На рис. 7 приведен синтезированный закон наведения спутника для выполнения маршрутов сканирующей съемки указанных видов и поворотных маневров между ними. Здесь учитывались ограничения 0) (/) 1.5 град/с и |е(/) |< 0.2 град/с2 на модули векторов угловой скорости и углового ускорения корпуса КА, обусловленные ограниченными ресурсами кластера ДМ по допустимой области вариации вектора его КМ и допустимым ускорениям ротора каждого ДМ.
Первые два набора графиков на рис. 7 представляют ориентацию ССК относительно ИСК (вектор МПР (5 = ¡а }) и относительно ОСК (углы крена ф1, синий цвет; рыскания ф2, зеленый цвет и тангажа ф3 = 0, красный цвет).
ИМИТАЦИЯ И АНИМАЦИЯ
Имитация выполнена для спутника землеоб-зорастензроминерции Л = Ла§{812; 587; 910} кгм2 на указанной ССО, в его модели (1) учтены упругие колебания каждой панели СБ с парциальными частотами первого тона О® » 0.6 р/с и декрементом 5 = 0.01. При шуме БИНС с СКО о" =1 угл. сек в процессе измерении угло-
и
7. Закон углового наведения КА для выполнения заданных маршрутов съемки
-0.5
о
-1
О) 20
ф
0
-е-
-20
ю 1 О)
-8 о
3 -1 0.2
N
W
О) п
a) U ~о
ы
-0.2
•1................Л..................1+................Г.................Г .............1'
550 600 650 700 750 800
Ъв
Рис. 8. Погрешности реализации закона наведения для заданных маршрутов съемки
Рис. 9. Фрагменты четырех кадров
анимации движения КА при выполнении задания
вого положения КА в ИСК с периодом Tq =1/8 с дискретная фильтрация измерений выполнялась с этим же периодом, а цифровое управление кластером ДМ формировалось по закону (5) с периодом Ти =1/4 с. Погрешности по углам 5фг и угловым скоростям ôû)„ полученные при компьютерном анализе реализации синтезированного закона углового наведения КА для выполнения заданных маршрутов съемки (см. рис. 7), представлены на рис. 8.
На рис. 9 приведены фрагменты четырех кадров анимации движения КА при выполнении указанного задания:
a) трассовая сканирующая съемка в надир, маршрут Триполи, ПМ, первый стерео MCI Тунис;
b) продолжение, ПМ, второй стерео МС2 Тунис, ПМ, первый стерео MCI Флоренция;
c) продолжение, ПМ, второй стерео МС2 Флоренция, поворотный маневр;
d) продолжение, съемка с выравниванием, маршрут Франкфурт, поворотный маневр.
Нетрудно убедиться, что только результаты анимации на рис. 9 позволяют интерпретировать (зрительно воспринять и объяснить) изменение ориентации спугника в процессе его пространственного движения, определенного законом углового наведения при выполнении заданных маршрутов съемки (рис. 7) с учетом орбитального движения центра масс космического аппарата.
ЗАКЛЮЧЕНИЕ
Кратко представлена разработанная компьютерная среда для анимации движения информационных спутников и продемонстрирована ее
работа для последовательности сканирующих маршрутов оптико-электронной съемки, включающей маршруты стереоскопической съемки заданных участков земной поверхности. Приведены результаты компьютерной имитации динамических процессов при угловом движении спутника землеобзора, управляемого минималь-но-избыточным кластером реактивных маховиков по сигналам астроинерциальной системы определения углового положения, в процессе реализации закона углового наведения КА, который синтезирован для выполнения заданных маршрутов сканирующей съемки.
СПИСОК ЛИТЕРАТУРЫ
1. SIRIUS-S software environment for computer-aided designing of attitude control systems for small information satellites / Ye.I. Somov, S.A. Butyrin, S.Ye. Somov, T.Ye. Somova // Proceedings of 20th Saint Petersburg International Conference on Integrated Navigation Systems. 2013. P. 325-328.
2. Somova T. Attitude guidance and control, simulation and animation of a land-survey mini-satellite motion I I Journal of Aeronautics and Space Technologies.
2016. Vol. 9, no. 2. P. 35-45.
3. Сомова Т.Е. Алгоритмы наведения и цифрового управления ориентацией маневрирующего спутника землеобзора // Известия Самарского научного центра Российской академии наук. 2016. Том 18. №4(6). С. 1173-1180.
4. Somova T.Ye. Guidance and digital attitude control of a maneuvering land-survey satellite // Proceedings of 24th Saint Petersburg International Conference on Integrated Navigation Systems. 2017. P. 564-568.
5. Сомов ЕЖ, Бутырин C.A. Наведение и гиросило-вое управление ориентацией спутника землеобзора при сканирующей стереоскопической съемке // Известия Самарского научного центра РАН.
2017. Том 19. № 4. С. 70-80.
COMPUTER SIMULATION AND ANIMATION CONTROLLED MOTION OF AN INFORMATION SATELLITE
© 2018 T.Ye. Somova
Samara State Technical University
We briefly present the results on application developed software environment for simulation and animation of a spatial motion of the information satellites, which is a part of the system for computer-aided designing of the spacecraft control systems. Keywords: spacecraft, controlled motion, simulation, animation
Tatyana Somova, Postgraduate Student, Associate Research Fellow of Department "Navigation, Guidance and Motion Control", Research Institute for Problems of Mechanical Systems Reliability, Samara State Technical University. E-mail te_somova@mail.ni