CC BY
64
УДК629.78 : 681.51
РАЗГОН ГИРОРОТОРОВ И ВКЛЮЧЕНИЕ ИЗБЫТОЧНОГО КЛАСТЕРА ГИРОДИНОВ В КОНТУР УПРАВЛЕНИЯ ОРИЕНТАЦИЕЙ ИНФОРМАЦИОННОГО СПУТНИКА
© 2017 С.Е. Сомов
Самарский научный центр Российской академии наук
Статья поступила в редакцию 31.08.2017
Предлагается и исследуется метод разгона гиророторов и приведения кластера гиродинов в парковое состояние без создания возмущающих моментов при угловой стабилизации корпуса спутника с помощью двигательной установки. Приводятся результаты исследования процессов при разгоне гиророторов, приведения в парк и включения силового гироскопического кластера в контур управления ориентацией спутника на солнечно-синхронной орбите. Ключевые слова: космический аппарат, управление ориентацией, кластер гиродинов, разгон роторов.
Работа поддержана РФФИ (гранты 17-08-01708,17-48-630637) и отделением ЭММПУ РАН (программа фундаментальных исследований № 13).
ВВЕДЕНИЕ
После отделения информационного спутника от ракеты-носителя и механического раскрытия его солнечных батарей (СБ) космический аппарат (КА) начинает вращаться с вектором угловой скорости ш произвольного направления в связанной с корпусом КА системе координат (ССК). Для ориентации КА в орбитальной системе координат (ОСК) выполняются начальные режимы, где обычно выделяют этапы успокоения (остановки вращения), определения углового положения КА по внешним ориентирам, поворотных маневров для приведения ориентации КА к заданной в ОСК и поддержания такой ориентации при подготовке бортового оборудования системы управления ориентацией (СУО) для последующего выполнения основных режимов. Например, начальные режимы ориентации геостационарного спутника связи SESAT [1] выполнялись с помощью датчиков угловой скорости, приборов ориентации на Солнце и Землю, двигательной установки (ДУ) на основе реактивных двигателей (РД) малой тяги с широтно-импульсным управлением для совмещения ССК с ОСК посредством двух одноосных поворотных маневров [2] и последующей угловой стабилизации КА в ОСК при разгоне ротора силового гироскопического стабилизатора в виде реактивного двигателя-маховика (ДМ) в двухстепенном кардановом подвесе [3].
В СУО маневрирующих информационных КА (в частности, спутников землеобзора) применяются электромеханические исполнительные органы в виде избыточных кластеров крат-
Сомов Сергей Евгеньевич, научный сотрудник отдела «Динамики и управления движением».
E-mail [email protected]
ных схем 3-SPE (группа 3 ножничных пар - 3 Scissored Pair Ensemble) или 2-SPE на основе соответственно трех либо двух пар двухстепенных силовых гироскопов - гиродинов (ГД) с колли-неарными осями подвеса [4,5]. Энергетические возможности системы электропитания КА в начале его миссии не позволяют выполнить одновременный разгон роторов всех ГД составе силовых гироскопических кластеров (СГК) указанных схем, доступен разгон гиророторов до
номинального значения h их собственного ки-
g
нетического момента (КМ) лишь для пар ГД.
В работах [6 - 9] представлен метод разгона роторов ГД в составе таких СГК, когда оси гиро-роторов и угловые положения ГД относительно осей их подвеса зафиксированы в корпусе КА с помощью арретиров сразу в «парковом» состоянии (парке) СГК, которому соответствует значение вектора его КМ H = Zh = 0 при h p =| h p |= hg , где индексом p отмечен номер ГД. Здесь применяется пошаговый разгон гиророторов в циклах последовательности пар ГД, при завершении каждого такого цикла выполняется условие H = 0. При этом проявляются инерционные возмущающие моменты из-за разгона гиророторов пар ГД и гироскопические связи каналов углового движения спутника, стабилизируемого в ОСК с помощью ДУ, естественно с расходом топлива РД.
В работе [10] предложена и исследована схема выполнения начальных режимов ориентации информационного спутника с применением бесплатформенной инерциальной навигационной системы (БИНС), корректируемой сигналами кластера звездных датчиков и навигационных систем ГЛОНАСС/GPS, двигательной установки на основе 8 РД с широтно-им-пульсным управлением, кластера четырех ДМ
и магнитного привода с цифровым управлением. Здесь успокоение КА выполняется в инер-циальной системе координат (ИСК) с помощью ДУ, далее проводятся инициализация БИНС с автономным определением как ориентации, так и расположения КА в ИСК, инициализация и включение кластера ДМ в контур управления ориентацией КА; этот режим завершается угловой стабилизацией КА в ИСК и расчетом на борту КА параметров потребного пространственного поворотного маневра КА на заданном интервале времени с краевыми условиями общего вида, который реализуется СУО с последующей угловой стабилизации спутника в ОСК.
В данной статье впервые предлагается метод разгона роторов ГД и приведения СГК кратных схем в парковое состояние без создания инерционных и гироскопических возмущающих моментов при угловой стабилизации корпуса спутника в ИСК с помощью двигательной установки. Приводятся результаты исследования процессов при разгоне роторов ГД, приведения в парк и включения СГК по схеме 3-БРЕ в контур управления ориентацией спутника на солнечно-синхронной орбите.
МОДЕЛИ И ПОСТАНОВКА ЗАДАЧИ
Стандартным образом вводятся ИСК, ОСК и ССК Охуг, корпус КА с фиксированным положением панелей СБ считается твердым телом с тензором инерции J. Вводятся также ортогональная каноническая система координат СГК ОхСуСгС, для простоты совпадающая с ССК Охуг, и отсчет углов ГД р относительно осей их подвеса, рис. 1. Используются обозначения со1(-) = {•}, Нпе(-) = [•], (•)', [ах] и °,г для векторов, матриц и кватернионов, индексы . = 1,2,3 = 1 - 3 для осей ССК, применяется вектор модифицированных параметров Родрига (МПР) а = {Сг-} = е tg(Ф /4) с обозначениями орта Эйлера е и угла ф
собственного поворота, который однозначно связан с кватернионом А ориентации КА в ИСК прямыми О — Х/(1 + ^0) и обратными 10 = (1 -ст2)/(1 + а2), Х = 2ст /(1 + ст2) соотношениями.
Модель углового движения КА принимается в виде
А = А ° ш/2; (1)
Лш = -[© х^ + М8 + Ме + Мй.
Здесь С = Лю + Н является вектором КМ электромеханической СУО, где столбец Н = {Нг.}
представляют КМ СГК; Ме - вектор момента двигательной установки на основе восьми РД [9], которая обладает возможностью одновременно создавать векторы внешних сил и моментов произвольного направления в ССК за счет ши-ротно-импульсной модуляции (ШИМ) тяги РД;
вектор М8 = {М8} управляющего момента СГК
формируется в виде М8 = — Н , где (•) - символ локальной производной по времени, а вектор Мй представляет внешние возмущающие моменты.
Каждый РД с номинальным значением тяги Рт, одинаковым для всех РД, имеет
шим тяги рр (г) = ртр№м(г -г;
V,
)
V/ е , /г+1), р = 1 - 8, с периодом Т^ и временным запаздыванием Т^и. Здесь ург является дискретным входным сигналом и функции
PWM(Г, хг, гт ,у рг) -
г (?т ) =
р г->
О
^пурГ г е { 0 г е + I <т
Г р г' гг+\)
р^т' И(Т;,|ург I)
р г\
1 I
р г1
где Хг = гТ:, ^ = Хг + Т:; г е N0 . [0,1,2,3...) . Вектор момента ДУ формируется по соотношению Ме = 2 [рр х] рр (г), р = 1 - 8, где р (г) = -рр (г)ер с ортом оси сопла е и рр -
Рис. 1. Каноническая система координат СГК и отсчет углов ГД в схемах 3-БРЕ (а) и 2-БРЕ (Ь )
р г
вектор тяги p -го РД и радиус-вектор точки ее приложения в ССК, соответственно.
Для схемы 3-SPE (см. рис. 1 a) введем обозначения проекций ортов КМ каждого ГД на оси
системы координат O xg yg zg СГК
xl = Cl = cos ßj; x2 = C2 = cos ß2; yl = Sl = sin ßt;
y2 = S2 = sin ß2; x3 - S3 s sinß3; x4 = S4 = sin ß4;
z3 - C3 = cosß3; z4 = C4 = cosß4; y5 - C5 = cosß5;
y6 = C = cosß6; г5 - S5 = sinß5; z6 - S6 = sinß6.
Тогда вектор нормированного КМ такого СГК h(ß) = Н(Р)/\, где столбец ß = {ßp}, p = 1 -6, и градиентная матрица Ah(ß) = 3h / 3ß представляются в виде
x rCi + C2 + S3 +S4"
h = y = Si +S2 +C5 + C6 У
z C3 + C4 +S5 + S6 _
Г" Si - S2 C3 C4 0 0
Ci C2 0 0 - S5 - S6
0 0 - S3 - S4 C5 S6
x C + c2 + s3 + s4 "
h = y = Si + S2
z _ C3+C4 _
Г- Si - S2 C3 C4
Ah(ß)= Q C2 0 0
0 0 - S3 - S4
законах их настройки; (И) разработке алгоритмов наведения и управления ориентацией КА с помощью ДУ при переводе его углового положения из произвольного, достигнутого при завершении режима успокоения, в такое угловое положение в ИСК, при котором нормаль к плоскости панелей СБ направлена на Солнце; (ш) разработке метода разгона роторов ГД и приведения СГК кратных схем в парковое состояние без создания инерционных и гироскопических возмущающих моментов при угловой стабилизации корпуса спутника в ИСК с помощью ДУ; (¡у) анализе динамических процессов при включении СГК в контур цифрового управления ориентацией спутника.
ПАРКОВЫЕ СОСТОЯНИЯ СГК
Вводятся обозначения
Х12 — Х1 ^ Х2 ' Х34 Х3 Х4 ' У\2 ~ У1 ^ У 2 '
у56 = у5+у6; ^34 = ^з+^56 = +
Ah(ß)=
Для схемы 2-SPE (см. рис. 1 b ) столбец ß = {ßp}, p = 1 - 4 и
x
12
■y2r.
, x34
34
А
; y i2 =
yi2
- z.
34
- x
12
y56 =
y56
• 'z — Z34 . ~ _ z56
' >Z34 I-7" 'Z56
- Z
56
- x
34
y526
При цифровом управлении СГК и®(/) = {ирк (/)}
с периодом Ти, где ирк О) = ирк Vt е У к, tk+1) , ^+1 = + Ти и к е К0, когда управлением считается вектор (3 = иг скоростей ГД относительно осей их подвеса, кусочно-непрерывный управляющий момент СГК формируется по соотношениям
М|(0 = -\Ли (Р(0 <(t); (3(t) = и| (t). (2) Сингулярные состояния этой схемы возникают при таких угловых положениях ГД, когда матрица Грамма G(P) = Ль (()Л^ (() теряет полный ранг, т.е. при G = det(G(P) = 0. Важной проблемой управления СГК избыточной структуры является выбор функции распределения вектора КМ СГК между ГД - закона настройки СГК. Наиболее рациональными [4,5] являются явные законы настройки, когда все характеристики движения ГД получаются по аналитическим соотношениям.
Задача статьи состоит в (¡) определении парковых состояний СГК обеих схем при явных
и компоненты явного векторного закона на-
стр°йки = {/Р1,/Р2,/рэ} = 0 СГК схемы
3-БРБ принимаются в виде [5]
/р\ (Р) = Х12 _ Х34 + Р(Х12Х34 _ 1) '
/>2(Р) = йб - й2 + р(Зу5бЗУ12 -1);
/>з(Р) = ^34 - *56 + Р(234256 - !) » где постоянный параметр р удовлетворяет условию 0 < р < 1. Данный закон настройки обеспечивает отсутствие сингулярных состояний СГК для всех внутренних точек области вариации вектора его суммарного КМ. Для представления условий однозначной разрешимости уравнения ) = h , где h = {х, у) z} является известным вектором-столбцом, относительно синусов и косинусов углов ( всех шести гиро-динов, вводятся обозначения
Рп =74-(Х12)2 ; д12 4-(у12)2 ;
Рз4 = 74^0^; д34 =л!4-(Х34)2;
Р56 = л/4-(У56)2 ; Ъб = л/4-(г56)2 ;
x + А x 2
y-л
' x34
y ■ Z =
' 34
x ~A x 2
z + Д.
' ' x56
y + A y 2
z-A z
y'2 ' 2 ' 2 ' J" 2 dx = Qn+PÍA', dy = q¡6+Pi2; dz = q34+Ps6-
(3)
Условия разрешимости векторного уравнения ) = ^) имеют вид
Ах = Лх{\ - [1 - 4Р((?12 - Рза)(Х/2) _ + р(д12р34 - (х/2)2))/Л2]1/2}/ р ' А, - -[1 - 4р((д56 -рп)(у/2)
+ Р(?56Р12 - (У/2)2))/<]1/2}/р; А- = {1 - [1 - 4р((д34 - р56)(г/2)
+ Р(^34р56 - (г/2)2))/^Г}/Р
и при введении столбца Д = {Ах, А , А } очевидным образом преобразуются к нелинейному векторному уравнению Д = Ф(И, Д). Получить аналитическое решение этого уравнения
г\ *
для определения углов р ГД при парковом состоянии СГК весьма затруднительно, но его численное решение достигается практически мгновенно по методу простой итерации - при рациональном выборе начального точки достаточно лишь 2-3 итераций для получения результата с приемлемой точностью.
Далее при введении обозначений
X + Лх . У -А У а = х; ь=-
У .
Л =-\] 4 - С]2 / ех\ а2 =-
ь =
х -А.
У + Л У
2 ^ 2 с2 = д^а2 ^ ь" 5 Л2 = д^4 с2" / с2 5 а3 =
-——; с3 ; Л3 = д/ 4 - с2 / с
Ь3 =
х = а! + Л|Ь, . у = Ь1 -2 2 '
2-я пара (ГД3 и ГД4):
Ь2 + Л2а2
Ь ^ 2 а 2
4 2 4 3-я пара (ГД5 и ГД6):
а 3 Ьз
у5 = ' 3; - =
Уб =
2
а3 + Л3Ь3
; - =
2 '
а2 - Л2Ь2
2
а2 + Л 2Ь2
2
Ь3 1 Л
2 ;
Ь3 Л
вычисление синусов и косинусов углов р всех шести ГД выполняется по явным аналитическим соотношениям 1-я пара (ГД1 и ГД2):
х = а- ¿А ; у = Ь1 + ¿а ; х1 2 ; у1 2 '
2 й 2
Положение нечетных (р = 1,3,5) и четных (р = 2,4,6) ГД в трех парах I = 1 ^ 3 удобно представить углами = (Р21-1 + Р21)/2 и 5. = (Р2.-1 - Р2.) / 2, 1 = 1 ^ 3 , которые поясняют применяемое выше наименование «ножничной пары» - угол ai определяет центральную линию а1 «ножниц», а углы ± 5^ - положения векторов КМ И2г-1 и h2г относительно линии а1 в 1 - ой паре ГД.
Парковое состояние СГК по схеме 3-БРЕ приведено на рис. 2 а. Здесь векторы КМ всех 3 пар ГД = И21-1 + И21 с одинаковыми модулями и концами в точках А, В и С принадлежат плоскости Р, которая содержит начало системы координат СГК, причем векторы И направлены по линиям а^.
Для СГК по схеме 2-БРЕ при скалярном законе настройки /р1(Р) = х12 - х34 + р(х12х34 -1) = О аналитическое решение, представленное первым соотношением в (3), зависит параметра р . Парковое состояние этой схемы приведено на рис. 2 Ь , где не показаны центральные линии ах и а2 двух ножничных пар, которые определяются углами а1 = 0 и а2 = к, а также углы 51 = 52 в парах ГД.
Рис. 2. Парковые состояния СГК кратных схем 3-БРЕ (а) и 2-БРЕ (Ь )
2
2
г .
НАВЕДЕНИЯ И УПРАВЛЕНИЕ КА ПРИ ОРИЕНТАЦИИ СБ НА СОЛНЦЕ
После окончания режима успокоения КА (остановки его вращения с точностью ~ 0,05 град/с), включения БИНС с определением ориентации КА, его расположения в ИСК, а также положения Солнца, замыкания контура стабилизации углового положения КА с широтно-импульсным управлением ДУ и завершения переходных процессов далее на некотором интервале времени t е [tp, tp ] должен выполняться разворот спутника в такое положение, при котором нормаль к плоскости панелей СБ направлена на Солнце. Соответствующий закон пространственного углового наведения КА в ИСК синтезируется на основе простой кинематической модели
Лp =Лp о юр/2; сЪр =еp
с заданными краевыми условиями
Л p (tp ) = Лр; ю p (tp ) = 0; Лp (tp) = Л p; ю p (tp) = 0
(4)
и ограничениями на модули вектора угловой скорости шр (/) и кусочно-постоянного вектора углового ускорения 8р (/). В результате получается угловое наведение КА по кратчайшему пути - при кватернионе разворота
Л* = X*) = Лр о Лр программное угловое движение КА представляется его вращением
вокруг орта е* = А* / sin(ф*/2) оси Эйлера на
угол ф* = 2агссоБ(А0).
Если задано программное угловое движение КА Лр(/),Юр(У),8р(0 = Юр(0 - закон его углового наведения в ИСК, то кватерниону рассогласования Е = (е0, е) = Лр ° Л с вектором е = {е} соответствует вектор параметров Эйлера Е = {е0, е}, матрица Се (Е) = 13 - 2[ех^е, где Qe = 13е0 + [ех], и вектор погрешности ориентации 5ф = } = {2е0е; }.
После дискретной фильтрации измеренных БИНС с периодом Тч значений вектора углового рассогласования е1 = -2е01 е1, I е N, формируются значения вектора 8 Г, г е N, дискретного
закона широтно-импульсного управления ДУ с
7
и •
5Г + 1
- Bgr + Cel; m r - Kgr + Pel;
(5)
М; = юг X !юг + !(Сегер + [С>рХ]юг + тг).
Здесь матрица Се = Се(Е) и используются диагональные матрицы К, В, С и Р. Далее импульс вектора потребного момента ДУ «пересчитывается» в столбец длительностей хг = {трг} включения РД V/ е [/г, tr+1), кото-
рые используются при ШИМ тяги всех восьми РД с формированием вектора управляющего
момента Ме(/) = 2 [р х] р (/) двигательной установки.
РАЗГОН РОТОРОВ И ПРИВЕДЕНИЕ СГК В ПАРКОВОЕ СОСТОЯНИЕ
Для СГК по схеме 3-БРЕ в плоскости изменения КМ каждой / -ой паре ГД введем линии Ь, ортогональные линиям а1, см. рис. 3 а, и будем считать, что направления осей роторов нечетных (р = 1,3,5) и четных (р = 2,4,6) ГД в парах / = 1 ^ 3 зафиксированы в корпусе КА с помощью арретиров по линиям при углах ГД относительно осей их подвеса Р2-1 = Р2 — л, так как согласно применяемой стратегии отсчета углов гиродинов нечетные ГД в парах всегда перемещаются впереди четных, см. рис. 1 и рис. 3 а.
При последовательном разгоне роторов в парах ГД с одинаковыми ускорениями вплоть до номинального значения их собственных КМ инерционные возмущающие моменты таких гиророторов компенсируются. Поэтому в результате КМ роторов всех шести ГД принимают значения без влияния их инерционных возмущающих моментов на угловое движение корпуса КА и при отключении арретиров гиродинов СГК становится готовым для приведения в парк.
При развороте двух ГД с противоположными векторами КМ в / -ой паре в разные стороны с одинаковыми скоростями относительно осей подвеса создаваемые ими гироскопические моменты складываются и вектор КМ ^ (/) этой пары ГД изменяется вдоль центральной линии , см. рис. 3 Ь . Поэтому в процессе приведения СГК в парк с условиями Н ) ^ 0 и М® ) = — Н ) ^ 0 рационально выполнять указанные развороты векторов КМ гиродинов во всех трех парах одновременно в следующей последовательности: 1) все шесть ГД разворачиваются в соответствующих направлениях с одинаковыми постоянными угловыми скоростями до
значений их углов р2, —1 = Р2,-—1 — X , Р2; = Р2,- + X
в парах / = 1 ^ 3 с заданным постоянным углом X, далее принимается % = 1 град; 2) задается потребный управляющий гироскопический момент СГК Мg = 0 и включается закон его настройки, который автоматически приводит СГК точно в парковое состояние.
Для СГК по схеме 2-БРЕ разгон гиророторов и приведение к парковому состоянию на рис. 2 Ь выполняется аналогичным способом. Здесь линия Ь, ортогональная центральной линии аг первой пары, параллельна оси ОуС , а линия Ь2 параллельна оси Ог® системы координат СГК
О ^ усС .
ВКЛЮЧЕНИЕ КЛАСТЕРА ГД В КОНТУР УПРАВЛЕНИЯ
После дискретной фильтрации значений
углового рассогласования вектор 8\ применяется в цифровом законе управления СГК
»к+1
= в§ к + Сек; тк = К§к + Рек;
(6)
М« = X СI +1 (с: ерК + [с: < х]шк + йк),
где G0 = Jшk + Нк. Далее вектор потребного
цифрового управления СГК М| с помощью явного закона распределения команд между шестью ГД «пересчитывается» в столбец и ® = {ирк } сигналов управления ГД, которые фиксируются на полуинтервалах цифрового управления СГК с периодом Ти для формирования управляющего момента М | (t) (2), детали представлены в [11].
РЕЗУЛЬТАТЫ ИМИТАЦИИ
При компьютерной имитации рассматривался информационный спутник с мас-
сой m = 1000 кг и тензором инерции J = diag{812; 587; 910} кг м2 на солнечно-синхронной орбите (ССО) высотой 720 км. Пусть при широтно-импульсном управлении с периодом TU = 4 с каждый из восьми РД в составе ДУ
имеет максимальное значение тяги Pm = 0.5 Н, а каждый из шести ГД в составе СГК имеет модуль собственного КМ hg = 100 Нмс и период цифрового управления Tu = 1 0.25 с. Будем считать, что к моменту времени t = tl = 7000 с завершен режим успокоения КА, с помощью БИНС определяются ориентация КА и его расположение в ИСК, а также положение Солнца, и на интервале времени t е [7000,7600] с выполняется угловая стабилизация КА в ИСК при широтно-импульсном управлении двигательной установкой. Закон углового наведения КА на интервале времени t е [t;p, tfP ] = [7600,8548] с для совмещения нормали к панелям СБ с направлением на Солнце при итоговом развороте корпуса КА вокруг оси Эйлера на угол ф* =76.6 град представлен на рис. 4. Здесь приведены программные значения компонентов векторов углового
Рис. 3. Разгон роторов в 3 парах ГД (а) и приведение СГК схемы 3-БРЕ в парковое состояние (Ь )
Рис. 4. Начальное положение и программное угловое движение КА при наведении панелей СБ на Солнце
Рис. 5. Ошибки стабилизации начального положения и движения КА при наведении панелей СБ на Солнце
ускорения 5р (г) и угловой скорости ш р (г), а также компонентов вектора модифицированных параметров Родрига Ор (г) . Достигаемые при этом ошибки стабилизации движения КА приведены на рис. 5.
В момент времени г = г/ = 8548 с начинается длительный режим угловой стабилизации КА с панелями СБ, ориентированными на Солнце. Здесь на интервале времени г е [8548, 14428] с попарно разгоняются роторы шести ГД без возникновения реактивных возмущающих моментов от этих гиророторов при стабилизации углового положение корпуса спутника в ИСК с помощью широтно-импульсного управления ДУ. В момент времени г = 14428 с положение векторов КМ с модулями = 100 Нмс в парах ГД определяется углами _ = 45 град и = -135 град, / = 1 ^ 3 .
При значении параметра р = 0.65 и парковом состоянии СГК по схеме 3-БРЕ углы р ГД имеют значения р2. _1 = 15.661816459787 град, р*2г =_105.661816459787 град, I = 1 -3 , что приводит к симметричным значениям углов а. = -45 град и 5г- = 60.661816459787 град для ножничных пар.
На интервале времени г е [14428, 14488] с все шесть ГД разворачиваются в соответствующих направлениях с одинаковыми постоянными угловыми скоростями,
что дает соблюдение условий Н (г) ^ 0 и М8 (г) = 0 . В результате углы р ГД принимают значения р2._1 = 14.661816459787 град, р2, = -104.6618164-59787 град, I = 1 3, которые сохраняются до момента времени г = 14528 с, когда включается закон настройки СГК с условием М8 = 0. В течение 20 секунд на интервале времени г е [14528, 14548] с указанный закон настройки приводит СГК точно в парковое состояние.
Начиная с момента времени г = 14548 с широтно-импульсное управление двигательной установки отключается и в контур угловой стабилизации спутника в ИСК включается СГК с цифровым управлением. Следует отметить, что структуры законов широтно-импульсного управления ДУ (5) и цифрового управления СГК (6) близки, но аналитические соотношения для расчета диагональных элементов матриц К, В , С и Р существенно отличаются.
Изменения углов р нечетных ГД (р = 1,3,5) при завершении приведения СГК в парк и включении его в контур управления представлены на рис. 6, а углов р четных ГД (р = 2,4,6) - на рис. 7.
Ошибки стабилизации КА в ИСК по углам 5ф. и угловым скоростям 5юг- на интервале времени г е [14548, 14648] с при включении СГК в контур цифрового управления приведены на рис. 8.
Рис. 6. Углы ГД при завершении приведения СГК в парк и включении его в контур управления
-104 8 -105 -105.2
О)
-g -105 4 -105 6 -1058 -106
1 _____
—h - -
\ :
¡
i
14 520 14 540 14560 14 580 14 600 14 620 14 640 Рис. 7. Углы ГД при завершении приведения СГК в парк и включении его в контур управления
14 530
14630
о.оз
0 02
О)
О)
"О 001
-0 01
-—
-5(Pj
o<o3
i i i i t
14 540
14 560
14 590
14 600
14 620
14 640
ts
Рис. 8. Ошибки угловой стабилизации КА при включении СГК в контур управления ЗАКЛЮЧЕНИЕ
Разработан новый метод разгона (увеличения скоростей вращения) роторов гироди-нов и приведения силовых гироскопических кластеров (СГК) двух кратных схем в парковое состояние без создания инерционных и гироскопических возмущающих моментов при угловой стабилизации корпуса спутника в ИСК с помощью двигательной установки. В этом методе учитывается ограниченность энергетических ресурсов системы электропитания КА в начале его миссии, гиророторы разгоняются до номинального значения их собственного КМ последовательно по парам ГД. При этом возникает большая длительность разгона роторов всех гиродинов (несколько часов) при угловой
стабилизации КА с панелями СБ, ориентированными на Солнце. Отсутствие возмущающих моментов в процессе такой инициализации СГК приводит к экономичности работы ДУ в отношении расхода топлива при длительной угловой стабилизации КА в ИСК.
Представлены методика определения парковых состояний СГК обеих схем при явных законах их настройки, разработанные алгоритмы углового наведения и управления ориентацией КА с помощью ДУ при переводе ориентации спутника в такое состояние, при котором нормаль к плоскости панелей СБ направлена на Солнце. Приведены результаты компьютерной имитации динамических процессов при включении СГК в контур цифрового управления ориентацией спутника.
СПИСОК ЛИТЕРАТУРЫ
1. Nonlinear dynamics of gyromoment attitude control system at communication satellite SESAT / Ye.I. Somov, S.A. Butyrin, V.A. Rayevsky, G.P. Titov, A. Baiget, R.C. Rogers // Proceedings of 5th IFAC Symposium "Nonlinear Control Systems". 2001. Vol. 3. P. 13991404.
2. Сомов С.Е. Анализ колебаний конструкции спутника при наведении на Солнце и Землю с широт-но-импульсной модуляцией управления двигателями // Известия Самарского научного центра РАН. 2007. Т. 9. № 3. С. 847-858.
3. Сомов С.Е. Динамика начальных режимов упругого геостационарного спутника с широтно-им-пульсным управлением двигателями и разгоном ротора силового гироскопического стабилизатора // Гироскопия и навигация. 2006. № 4(55). С. 100.
4. Platonov Steering the spacecraft control moment gyroscope clusters / Ye.I. Somov, S.A. Butyrin, A.V. Sorokin, V.N. // Proceedings of 10th Saint-Petersburg International Conference on Integrated Navigation Systems. 2003. P. 403-419.
5. Сомов Е.И. Анализ сингулярных состояний и синтез явных законов настройки гирокомплексов кратных схем // Гироскопия и навигация. 2013. № 1(80). С. 134-148.
6. Somov S.Ye. Guidance and robust pulse-width control of an information satellite at initial modes // Proceedings of 18th IFAC Symposium on Automatic Control in Aerospace. 2010. P. 302-307.
7. Somov S., Butyrin S. Guidance and robust pulse-width control of an information satellite at initial, service and emergency modes // Proceedings of 18th IFAC World Congress. 2011. P. 2072-2077.
8. Somov S. Spacecraft attitude pulse-width control at initial, service and emergency modes // AIP Conference Proceedings. 2012. Vol. 1493. P. 933-940.
9. Somov S. Pulse-width attitude control of large-scale information satellite at service modes // Proceedings of 19th IFAC Symposium on Automatic Control in Aerospace. 2013. P. 7-12.
10. Сомов Е.И., Бутырин С.А. Наведение и управление ориентацией информационного спутника в начальных режимах // Аналитическая механика, устойчивость и управление. Труды XI международной Четаевской конференции. Том 3. Секция 3. Управление. Ч. III. Казань: КНИТУ-КАИ. 2017. С. 41-52.
11. Сомов Е.И., Бутырин С.А. Полетная идентификация и силовая гироскопическая стабилизация слабо демпированной конструкции крупногабаритного спутника // Проблемы управления. 2013. № 2. С. 51-57.
ACCELERATION OF GYROROTORS AND SWITCHING ON A REDUNDANT GYRODINE CLUSTER IN A LOOP FOR ATTITUDE CONTROL OF INFORMATION SATELLITE
© 2017 S.Ye. Somov
Samara Scientific Centre, Russian Academy of Sciences
We have suggested and studied the method for acceleration of gyrorotors and to bring the gyrodine cluster to its park state without the creation of disturbing torques during angular stabilization of the satellite body by a propulsion unit. We have presented the research results on processes of acceleration the gyrodine rotors, the bringing the gyro moment cluster into the park and its switching on in a loop for the attitude control of an information satellite in sun-synchronous orbit. Keywords: spacecraft, attitude control, а gyrodine cluster, acceleration of rotors.
Sergey Somov, Research Fellow at the Dynamics and Motion Control Department. E-mail s [email protected]
