CC BY
41
Механика и машиностроение
УДК 531.383 : 629.78 : 681.51
ОЦЕНКА РЕАЛИЗУЕМОСТИ ПОВОРОТНОГО МАНЕВРА КОСМИЧЕСКОГО АППАРАТА ПРИ НЕОПРЕДЕЛЕННОСТИ НАКОПЛЕННОГО КИНЕТИЧЕСКОГО МОМЕНТА СИЛОВОГО ГИРОКОМПЛЕКСА
© 2008 Е.И. Сомов, С.С. Мещеряков
Самарский научный центр РАН
Представляются методы априорной гарантирующей оценки реализуемости поворотного маневра космического аппарата при неопределенности накопленного кинетического момента силовых гироскопических комплексов кратных схем на основе четырех и шести гиродинов.
Введение
При среднесрочном наземном планировании работы КА информационного назначения (спутников землеобзора, метеорологии, связи, геодезии, навигации и т.д.) очень актуальной является проблема априорной оценки реализуемости очередного поворотного маневра (ПМ) космического аппарата (КА) с помощью силовых гироскопических комплексов (СГК) с учетом неопределенности вектора накопленного кинетического момента (КМ) [1]. В статье представляются методы решения этой проблемы применительно к СГК двух кратных схем на основе четырех и шести гиродинов.
Силовые гирокомплексы кратных схем
В классе гиродинных систем масса СГК будет минимальной, если при наименьшем числе m > 4 используемых гиродинов (ГД) он обеспечивает потребную область вариации суммарного КМ без непроходимых сингулярных состояний. Как известно [2], примене-
ние групп (в простейшем случае пар) ГД с коллинеарными осями подвеса дает очень важное преимущество при формировании схемы СГК - начиная с числа ГД m = 6 сингулярные состояния такой схемы во всех внутренних точках области S являются строго проходимыми. Коллинеарная пара безу-порных ГД в оригинальной работе J.W. Crenshaw (1973) названа Scissored Pair Ensemble (SPE), а избыточные кратные схемы на основе трех и двух коллинеарных пар ГД - как 3-SPE и 2-SPE соответственно. Некоторые внутренние сингулярные состояния схемы 2-SPE на основе четырех ГД являются лишь условно проходимыми, поэтому при применении этой схемы приходится существенно ограничивать используемую область вариации суммарного КМ либо привлекать сложные адаптивные законы настройки СГК. Достоинство кратной схемы СГК на основе 6 ГД, представленной двумя коллинеарными группами - по три ГД в каждой группе, заключается в ее отказоустойчивости, но она
799
Известия Самарского научного центра Российской академии наук, т. 10, №3, 2008
Рис. 3. Оболочка области вариации КМ схемы 3-SPE
“проигрывает” схемам 2-SPE и 3-SPE по сложности алгоритмов настройки. Поэтому в классе гиродинных систем наиболее рациональными схемами СГК являются схема 3-SPE на базе 6 ГД (рис. 1) и минимально избыточная схема 2-SPE (m = 4, рис. 2). На этих рис. представлены канонические структуры указанных схем, когда оси подвесов ГД параллельны осям ортогонального канонического гироскопического базиса (КГБ) Gc (Ox g yg zg), который будем считать совпадающим с базисом B(Oxyz) связанной с корпусом КА системе координат (ССК). Оболочки областей S вариации нормированного КМ данных схем представлены на рис. 3 и рис. 4.
Математическая постановка задач
Схема 2-SPE является частным случаем схемы 3-SPE (при отсутствии в ней третьей
Рис. 4. Оболочка области вариации КМ схемы 2-SPE
пары ГД), поэтому далее более подробно выполняется анализ схемы 3-SPE и указываются дополнительные аналитические результаты, которые удается получить для схемы 2-SPE. Свяжем с каждым p = 1 ^ m из m ГД, имеющих вектор КМ Hp = hg hp c одинаковым модулем hg , правый триэд его осей: орт h p (Р p) вектора КМ ГД, положение которого определяется углом Рp , фиксированный в ССК орт gp оси подвеса ГД и орт
p p (Р p) = h p (Р p) x g p.
Введем вектор-столбец P = {Р p}, составленный из углов поворота ГД относительно осей их подвеса, см. рис. 1 и рис. 2. Тогда вектор суммарного нормированного КМ СГК h(P) = х h (Р ), а парковые состояния СГК, когда вектор нормированного КМ СГК h(P) = 0 , представлены на рис. 5 и рис. 6. Орт mp(рp) = ^hp(рp)/дрp = gp x hp(pp) гирос-
Рис. 5. Парковое состояние схемы 3-SPE
800
Механика и машиностроение
копического момента p -го ГД всегда противоположен орту pp(Pp).
Введем кватернион Л = (Л,0,X),
X = {Хг., i = 1 - 3} ориентации связанного с корпусом КА базиса в относительно некоторого инерциального базиса I и вектор ю = {шг-} абсолютной угловой скорости КА. Все векторы и тензор инерции корпуса КА J представляются в ССК и при m > 4 определяется прямоугольная матрица Якоби
Ah (Р) = ^h / ЭР = [m g (Р p ), Р = 1 - m].
Пусть Vt e Tr = [t;, tf] заданы функции времени Л^), ю(t), s(t) = со(t) и
Ё(t) = е*(t) + o(t)XE(t) , представляющие программное угловое движение базиса B относительно базиса I . В рамках прецессионной теории силовых гироскопов при отсутствии внешнего возмущающего момента угловое движение КА удовлетворяет векторному интегралу нормированного КМ механической системы “КА+СГК” в виде
начальный момент времени t = t- считается известной, но направление этого вектора неопределенно, т.е. оно может быть произвольным. Модуль g вектора gi принимается постоянным и максимально допустимым, что справедливо при существенно меньшем темпе накопления КМ по сравнению с темпом поворотного маневра КА.
Решаемые в статье задачи состоят: в создании явных функций настройки СГК (распределения суммарного КМ СГК между гиродинами) для исключения избыточности и возможных сингулярных состояний;
в разработке метода оперативной бортовой гарантирующей оценки реализуемости очередного ПМ КА при известном (измеренном и вычисленном) векторе накопленного КМ;
в разработке метода априорной наземной гарантирующей оценки реализуемости последовательности ПМ КА на заданных смежных интервалах времени при неопределенности направления накопленного КМ.
g(0 = k(<)+И(Р(0) = Л() о gX Л(t), С1)
где нормированный КМ корпуса КА k(t) = Jo(t)/hg , а g) = Л(ti) о g(t;) о Л(^) -вектор накопленного КМ, постоянный в инерциальном базисе I. Угловое движение КА описывается дифференциальными уравнениями
Л = 2Л о ю;
k (t) + o(t) x g(t) = -m g (Р,Р). (2)
Здесь вектор k(t) = JO (t)/ hg и нормированный вектор управляющего момента СГК mg (Р,р) вычисляется в КГБ по формуле
mg (Р, р) = ^ mgp (Рp ) рp = Ah (Р) ug;
Р = ug; Ug = Р = vg,
Явные законы настройки СГК
Введем обозначения проекций ортов КМ каждого ГД на оси ортогонального КГБ G c(Ox g y gzg), см. рис. 1 и рис. 2:
x1 = C1 = cosP1; x2 = C2 = cosP2;
У1 = S = sin Pi; У2 = ^2 = sin P2;
X3 = S3 = sin P3; x4 = S4 = sin P4;
Z3 = C3 = cos P3; z4 = C4 = cos P4;
y5 = C5 = cos [З5; y6 = C6 = cos P6;
Z5 = S5 =sin P5; z6 = S6 =sin P6.
Тогда вектор-столбец нормированного суммарного КМ СГК h в ортогональном КГБ Gc и градиентная матрица Якоби Ah (Р) = 5h / Эр представляются для схемы 3-SPE в виде
где “управлением” считается вектор-столбец Р = ug = {ир}с компонентами up скоростей прецессии ГД, ограниченных по модулю заданным значением um = const. При этом вектор-столбец производных “управления” Ug = Р = vg = {vp} имеет компоненты vgp угловых ускорений ГД относительно осей подвеса, также ограниченные по модулю заданным значением vmg = const.
Верхняя оценка модуля gi =| gi |< g0 вектора gi = g(ti) накопленного КМ в ССК в
x C+c2+s3+s4 "
h = У = S1 + S2 + C5 + C6
z _C3 + C4 + S5 + S6 _
- S1 - S2 C3 C4 0 0
Ah(Р) = Q C2 0 0 -S5 -S6
0 0 -S3 -S4 C5 S6
а для схемы 2-SPE - в виде
801
Известия Самарского научного центра Российской академии наук, т. 10, №3, 2008
x 'C1 + C2 + S3 + S4"
h = У = S1+S2
z _ C3 + C4 _
"- S1 - S2 C3 Q
Ah (Р) = Q C2 0 0
0 0 - S3 - S4
Сингулярные состояния СГК возникает при таких положениях ГД, когда матрица Грамма G(h (Р)) = G(P) = Ah (Р)АТ(Р) теряет полный ранг, т.е. при G = det(G(P)) = 0 . Определитель G(P) матрицы Грамма представляет объем области вариации вектора нормированного управляющего гироскопического момента mg (Р,р) , соответствующего значению h (Р).
Угловые скорости и ускорения каждого ГД относительно оси подвеса ограничены по
m m
модулю заданными постоянными u и v , а именно
|р, (t)|< um ;|р, (t)|< vm t e Tr, (4) что соответствует ограниченным ресурсам его привода. Такие ограничения существенно нелинейным образом (3) “трансформируются” в ограничения на компоненты вектора нормированного управляющего момента СГК mg (Р,р) и производной от него по времени, т.е. моментные характеристики СГК зависят как от предыстории, так и текущего расположения векторов КМ гиродинов в его составе.
Вводятся обозначения
x12 = x1 + x2; x34 = x3 + x4;
'34
= y + y2; y56 = y5 +y6;
= z3 + z4; z56 = z5 +z6;
= x12 ~34 9 x34
1 1 1 J4 - z^4 ;
1 V; -- ~56 9 y 56
1 1 X J4 - z56 '
= z34 ~56 9 z56
V4 - x34 V4 - y56
Компоненты явного векторного закона настройки fp (p) = fp2, /р3} = 0 СГК схе-
мы 3-SPE принимаются в виде
34
12
fpx (Р) = /р1(Р) = ~12 - ~34 + р(~12~34 - j) ; fpy (Р) = fp2 (Р) = ~56 - ~12 + Р(^5б.~12 - 1 ; (5)
fpz (Р) = /р3 (Р) = ~34 - ~56 + Р(~34~56 - 1) ,
где постоянный параметр р удовлетворяет условию 0 < р < 1.
Для значении р = 2л/2/3 « 0.9428 парковое состояние схемы 3-SPE определяется аналитически: углы ГД в парах составляют + 45° относительно плоскости, которая имеет нормаль, совпадающую с главной биссектрисой координатных плоскостей КГБ, см. рис. 5. При обозначении a = 1^/2 такому парковому состоянию СГК соответствуют значения
x1 = x2 = a; y1 = y2 = -a ;
z3 = z4 = a ; x3 = x4 = -a;
y 5 = y 6 = a; z5 = z6 = -a.
При значении р « 0.87 в парковом состоянии СГК по схеме 3-SPE достигается глобальный максимум определителя G(h) = G(x, y, z), именно gl max G = Gm = G(0,0,0) = 8, однако даже в сечении z = 0 области S вариации нормированного КМ СГК проявляется существенно несимметричный и многоэкстремальный характер зависимости G(x, y,0), (рис. 7), в частности, имеется глубокий провал с дном G(2,-2,0)« 3.3.
При р = 0.65 в парковом состоянии такого СГК значение G(0,0,0) = 6, а в сечении z = 0 области S вариации КМ зависимость G(x, y,0) имеет вполне приемлемую несимметричность и многоэкстремальную “волнистость”, (рис. 8).
Важное достоинство векторного закона (5) настройки СГК заключается в возможности аналитического симметричного “размещения” внутри области S двух сфер, а именно сферы s g и сферы Sg, удовлетворяющих следующим условиям:
1. радиус rg сферы sg соответствует потребной области вариации нормированного КМ СГК с учетом возможного накопленного КМ;
2. радиус r0g = rg - g0 сферы Sg соответствует потребной области вариации нормированного КМ СГК без учета накопленного КМ;
3. для h (Р) edSg (на поверхности сферы sg ) определитель Грамма G(P)
802
Механика и машиностроение
Рис. 7. Значения определителя G матрицы Грамма в сечении z = 0 при р = 0.87
Рис. 8. Значения определителя G матрицы Грамма в сечении z = 0 при р = 0.65
принимает положительные значения G(h (P))=G (h (P))> 0, причем d Gm < Gg (h (P)) << Gm, где Gm = gl minG(P) 3h (P) e Sg и постоянный положительный коэффициент d << 1;
4. при h (P) e Sg / Sg (в шаровом слое Sg / Sg ) определитель Грамма G(h (P)) монотонно уменьшается с возрастанием модуля нормированного вектора h (P) КМ СГК.
Для компактного представления условий однозначной разрешимости вытекающего из
(1) векторного соотношения
h(P(t)) = h(t) = g(t) - k (t), (6)
где h(t) = {x(t), y(t), z(t)} является известной векторной функцией, относительно синусов и косинусов всех углов Р p (t) гиродинов, вводятся обозначения
Р12 = V4 (Х12)2 ' ^ 412 = V4 (У12 )2
Р34 = V4 - (Z34)2 : > 434 = V4 - (Х34)
803
Известия Самарского научного центра Российской академии наук, т. 10, №3, 2008
V4 (>56) ; q56 -V4 (z56) ;
x-А.
Р56 — V^56
x + А
xi —
a, - db b + da
; >i —
x — ________x x — x
Л12 ~ ; Л34
2
a + db
2
y +2А
Z+2A
2
> y-А>
>; >12 — > ■
x — -4 ' “l"! > — b1 diai , 2 2 ; >2 2 ;
вторая пара (ГД 3 и ГД 4): a - db
Z34 — ■
Z56 —
dx — 4l2 + Р34; d> — 456 + Р12;
dZ — 434 + Р56 .
Условия однозначной разрешимости уравнения (6) имеют вид
Аx — dx{1 - [1 - 4Р((412 - Р34)(x/2)
+ Р4Р34 - (X/2)2))/d2]1/2}/Р ;
А> — d>{1 - [1 - 4Р((456 - Р12 )(У /2)
+ р(456Р12 - (У/2)2))/d>2]1/2}/р ;
Аz — dz {1 - [1 - 4Р((434 - Р56 )(z / 2)
+ Р(?,4Р„ -(z/2)2))/d;]1'2}/р (7)
и при введении вектора-столбца A(t) — {Ах (t), А> (t), Аг (t)} очевидным образом преобразуются к нелинейному векторному уравнению A(t) — 0(h(t), A(t)) . Получить аналитическое решение этого векторного уравнения весьма затруднительно, но его численное решение достигается практически мгновенно по методу простой итерации - при рациональном выборе начального точки итерационного цикла достаточно лишь 1-2 итерации для получения результата с приемлемой точностью. Далее вычисление синусов и косинусов углов Р Р (t) всех шести гиродинов выполняется по явным аналитическим
2 ’ x — b2 + d2a2 z x — 2 ; Z3
z-А z
2 ; x — b2 - d2a2 z x — 2 ; Z4
2
a + db
соотношениям x + А
a1 —
x; b —
У-А > 2
третья пара (ГД 5 и ГД 6):
_a3 - d3b3 _b3 + d3a3
>5 — 2 ; Z5 — 2 ;
_a3 + d3b3 _b3 - d3a3
>6 — 2 ; Z6 — 2 .
При явном законе настройки (5) схемы 3-SPE принципиально отсутствуют сингулярные состояния СГК для всех внутренних точек области вариации S его суммарного вектора КМ h(t). Значения углов РР (t) е [0,2л] гиродинов представляют интерес только для графической интерпретации результатов, явное однозначное аналитическое представление этих углов получается очевидным образом.
Для схемы 2-SPE используются единственная скалярная функция распределения нормированного КМ СГК между первой и второй парами гиродинов
fpx (Р) — /р1 (Р) = ~12 - ~34 + Р(~12~34 - 1) — 0 (8)
и соотношения
x + А„
x12 — x1 —
x ■ V — x — .
34 2
x -А,
2 34 2 2 Р12 — V4 - (x12) ; q34 — V4 - (x34 ) ;
>12— у ; 412— 4у — V4 - (У)2;
Z34 — z ; Р34 — pz - л14 - (z)2 ; dx — 4у + Рz .
Это позволяет получить явную однозначную аналитическую зависимость
Х56 _
2
c1 —tJa2 + bj2 ; d1 — ^4 - c2 / c1;
z + А,
a2 —
2
; b—
2
x - А x 2
C2 '\la2 + b2 ; d2 V4 ~C2 / C2 ;
a3 —
У + A > 2
; b—
3
z-А
2
; d^ Cg” / C3 ;
первая пара (ГД 1 и ГД 2):
Аx — dx{1 - [1 - 4Р((4у - Рz)(x/2)
+ Р(4уpz - (x/2)2))/d2]1/2}/р
и затем аналитическое представление как синусов и косинусов углов Р Р (t), так и собственно значения этих углов всех четырех гиродинов. При законе настройки (8) СГК данной схемы внутри области вариации его суммарного нормированного КМ остается лишь множество сингулярных состояний Q >z (Р) — qp u Q Р - два не пересекающихся
804
Механика и машиностроение
между собой подмножества в виде двух полуэллипсов в сочетании с сингулярными направлениями вариации КМ пар гиродинов
QP = Q* ^ S*;
(9)
S, = {s = 0;| st |=| s21= 1, s = y, z;
Q* = {(х24/2р))2 + (z/2)2 = 1;Х34 < 0};
Q* = {(Х12 /2р))2 + (у/2)2 = 1;Х12 > 0}.
В парковом состоянии СГК схемы 2-SPE глобально оптимальной как по объему области вариации нормированного управляющего гироскопического момента mg, так и по гарантированной величине его модуля в произвольном направлении, является конфигурация
Ри = ±РP; р3,4 =-(я/2)±РP, где при значении р = 2V6 / 5 угол РP = arcco^/273 = 35°15', (рис. 6). В этом состоянии достигается глобальный максимум gl max G = Gm = 64/27^ 2.37 определителя G матрицы Грамма. Значение этого определителя монотонно уменьшается с возрастанием модуля нормированного КМ СГК и обращается в нуль на множестве сингулярных состояний (9).
Расчет характеристик движения ГД
Наличие векторного интеграла КМ (1) и применение известных правил дифференцирования вектора с учетом подвижности связанного базиса в относительно инерциального базиса I позволяет получить явное аналитическое представление вектора h(t) е S и его производных по времени t е Tr для произвольного вектора g(t) = A(t) ° g) ° Л(t):
h(t) = g(t) - k (t); h (t) = g1(t) -k(t); h(t) = g 2(t) - k(t)
(10)
где
g1(t) = -©(t) х g(t);
g 2(t) = -c6(t) х g(t) - ю(t) х g1(t), кватернион Л(t), векторы ю (t), CO(t) = s(t), со (t) = e (t) и k (t) = Jo(t)/hg, k (t) = Js(t)/ hg, k(t) = Je (t)/ hg
являются известными функциями времени. В результате получается явное аналитическое представлению вектора h(t) = {x(t), y(t), z(t)} и его производных в КГБ.
В свою очередь аналитические зависимости
h(t) = ah / ap p = Ah (P) p;
h(t) = Ah(P) p+[a(Ah(P) p)/ap ] p
с учетом производных компонентов (5) явного векторного закона настройки СГК по схеме 3-SPE для индексов s = 1 ^ 3 позволяют получить явные формулы для однозначного аналитического расчета векторов-столбцов p (t) и p(t). Здесь основная задача состоит в аналитическом обращении матрицы
(a/р, (p)/ ap , p) = 0
afps (p)/ap, p+ a([a/p, (p)/ ap] p)/ ap, p) = o,
A(p) =
Ah (p)
a/P1(p)/ ap a/P2(p)/ ap a/P3(p)/ ap
- У1 - У 2 Z3 Z 4 0 0
x1 X2 0 0 - Z5 - Z6
0 0 - x3 - X4 У5 У6
~41 ~~42 ~~43 ~44 0 0
~51 ~~52 0 0 ~~55 ~56
0 0 ~~63 ~64 ~65 ~66
с учетом 6 тождеств
X12 III + 1; X22 +У22 = 1;
x32 + z32 = 1; X42 +z 4; = 1;
У52 + Z52 = 1; у6 +z62 e 1.
При обозначениях
S12 = 412 - -px12; s34 = 434 pZ34 ;
S56 = 456- -РУ56; V12 = P12 + РУ12 :
V34 = P 34 +РХ34; V56 = P56 + PZ56
Г12 = 1 + У1У2 + X 1 X2;
805
Известия Самарского научного центра Российской академии наук, т. 10, №3, 2008
r34 = 1 + z3 z4 + x3 x4 ;
r56 =1 + z5 z6 + У 5 ye;
f (p, p)={(am„ (p)/ap] p)/ap, p), s = 1«■3}
получаются аналитические соотношения
4 y1 - r12 y12
«41 = -V34
«42 V34
p34(q12)
4 У 2 ~ r12 У12 _
p34(q12) 5
«43 S12
4 Z3 Г34 Z34
q12 (p34)
«44 S12
4 Z4 Г34 Z34
q12( p34) 5
«51 S56
4 X1 -r12 X12 q56 (p12)
«52 S56
4 X2 r12 X12
q56(p12) ’
«55 = -V12
4Z5 -r56 Z56 p12 (q56)
«56 = -V12
4 Z6 -r56 Z56 _ p12(q56) ’
«63 = - V56
4 x3 -r34 x34
p56(q34) ’
«64 = -V56
«65 S34
4 x4 -r34 x34
p56 (q34 ) ’
4У5 -r56 У56 _ q34 (p56 ) ’
«
66
= -s„
4У6 -r56 У56
^34( P56)3
матрица D(p) = (A(p)) 1 и искомое аналитическое представление векторов-столбцов
p(t) = D(p(t)){h(t),0} ; (11)
p(t) = D(p(t)){h(t),-f(p, p)}. (12)
Оценка реализуемости ПМ
Бортовая гарантирующая оценка реализуемости поворотного маневра КА при известных кинематических параметрах A(t), ю (t), s(t), ё (t) движения его корпуса и известном (измеренном и вычисленном) в ССК в начальный момент времени t = t. фиксированном векторе g. = g(t.) накопленного КМ для СГК
по схеме 3-SPE получается проверкой для t е Tr условия h(t) е Sg и далее по соотношениям (7) и (10) - (12) с явной проверкой выполнения условий (4) для каждого p -го гиро-дина, p = 1 ^ 6. Здесь все вычисления осуществляются по аналитическим соотношениям, за исключением численного решения векторного уравнения A(t) = 0(h(t), A (t)) (7), на сетке дискретных моментов времени из заданного интервала, причем такая сетка сгущается в окрестности моментов времени, когда компоненты вектора производной ускорения Ё (t) обращаются в нуль либо происходит конечный разрыв их непрерывности.
В случае применения СГК по схеме 2-SPE все вычисления для оперативной бортовой оценки реализуемости ПМ КА выполняются по явным аналитическим соотношениям.
Наземная априорная оценка реализуемости ПМ КА с учетом неопределенности направления вектора накопленного КМ осуществляется в два этапа. Согласно (10) наличие накопленного КМ приводит к смещению траектории вектора h(t) КМ СГК на величину вектора g(t) = A(t) ° g) ° Л(t) . Введем обозначение орта h(t) вектора h(t).
На этапе 1 определяется набор моментов времени TR = {t;. е Tr}, в которые достигаются максимальные значения модулей потребных угловых скоростей (3 p (t) и ускорений (3 (t) каждого ГД, когда вектор нормированного КМ СГК h(t) изменяется по закону h(t) = -k (t) , т.е. когда вектор накопленного КМ g. = g(t.) = 0. Здесь на сетке дискретных моментов времени t е Tr по указанной выше методике проверяется справедливость условия h(t) е Sg и определяются как множество моментов времени TR, так и орты h j = h(t;.) вектора h(t;.) е Sg , имеющих наихудшее влияние на реализуемость ПМ КА в отношении угловых скоростей |3 p (t) и ускорений 3 (t) гиродинов.
На этапе 2 выполняется повторный расчет максимальных значений модулей угловых скоростей и ускорений прецессии каждого ГД при t е Tr для набора нормированных векторов накопленного КМ gO, определенных в инерциальном базисе I соотношением
glj = go A(tj ) ° h(t; ) ° A(t; ) , tj е TR .
Если получаемые в результате такого
806
Механика и машиностроение
расчета максимальные значения модулей угловых скоростей P (t) и ускорений (3 (t) гиродинов при каждом значении вектора накопленного КМ g);. удовлетворяют ограничениям (4), то поворотный маневр КА считается гарантированно реализуемым, в противном случае нереализуемым - указываются конкретные “критические” направления максимального по модулю накопленного КМ, при которых нарушаются условия (4) для конкретных ГД в составе СГК.
Заключение
Разработаны и проанализированы явные функции настройки СГК кратных схем на основе четырех и шести гиродинов при их каноническом расположении на корпусе КА. Созданы конструктивные вычислительные алгоритмы оперативной бортовой и априорной наземной гарантирующих оценок реализуемости ПМ КА, в том числе при неопределенности накопленного КМ.
Работа поддержана РФФИ (07-08-97611, 08-08-99101, 08-08-00512), Президиумом РАН (программа фундаментальных исследований 22) и Отделением энергетики, механики, машиностроения и процессов управления РАН (программы 15 и 18)
СПИСОК ЛИТЕРАТУРЫ
1. Сомов Е.И., Герасин И.А. Оценка реализуемости поворотного маневра космического аппарата, управляемого избыточной системой гиродинов // Управление движением и навигация летательных аппаратов. Том 2. Самара: СГАУ, 1997.
2. Сомов Е.И., Бутырин С.А., Сорокин А.В., Платонов В.Н. Управление силовыми гирокомплексами космических аппаратов / / Труды X Санкт-Петербургской международной конференции по интегрированным навигационным системам. С.-Петербург: ЦНИИ “Электроприбор”, 2003.
ESTIMATION OF REALIZABILITY FOR A SPACECRAFT ROTATIONAL MANEUVER AT UNCERTAINTY BY AN ACCUMULATED ANGULAR MOMENTUM OF GYROMOMENT CLUSTER
© 2008 Ye.I. Somov, S.S. Meshcheryakov
Samara Science Centre of Russian Academy of Sciences
A rotational maneuver of a spacecraft with the gyromoment clusters by multiply schemes based on six and forth gyrodines, is considered. Methods for a priori guarantee estimation of such maneuver’s realizability at uncertainty by the clusters’ accumulated angular momentum, are presented.
807
