13 / 23
Методика формирования устойчивых околокруговых солнечно-синхронных орбит…
Инженерный журнал: наука и инновации
# 6·2017 13
сут
c З др
2
(
).
L
N Т
∆ = π − ω − Ω
Расчет фазового вектора КА в абсолютной системе координат
(АСК) на момент
0
t
в восходящем узле орбиты осуществляется по
формулам:
1
1
p x
e
=
+
;
0
y
=
;
0;
z
=
2
х
V
e
p
µ = −
;
1
(1 ) cos
у
V
e
i
p
µ = +
;
1
(1 ) sin
z
V
e
i
p
µ = +
.
(21)
Задача формирования параметров околокруговых программных
геометрически устойчивых кратных ССО формулируется следующим
образом. Задана кратность программной орбиты. Число учитываемых
в модели ГПЗ зональных гармоник
8.
n
≥
Требуется определить
в
восходящем узле начальные значения параметров программной гео-
метрически устойчивой ССО:
0 0 0 0
, , ,
a e i
ω
. В расчетах принимается
0
0.
Ω =
В дальнейшем требуемое значение долготы
0
Ω
восходящего
узла орбиты в момент времени
0
t
обеспечивается расчетным време-
нем запуска КА.
С учетом вышеизложенного
алгоритм расчета параметров гео-
метрически устойчивой ССО с заданной кратностью замыкания
трассы полета
принимает следующий вид.
1. По формуле (15) определяется номинальное значение
н
др
T
ССО
k
-суточной кратности с замыканием трассы через
l
витков.
2. По формуле (18) рассчитывается оскулирующее значение
большой полуоси
( )
j
a
в восходящем узле орбиты. Полагается, что
1 ( )
2 ( )
0.
= =
j
j
e
e
Номер итерации
1
j
=
.
3.
По формуле (17) вычисляется фокальный параметр
j
p
.
4.
По формуле (19) определяется наклонение
j
i
плоскости ССО.
5.
По формулам (13) находятся значения элементов
1 ( )
j
e
и
2 ( )
j
e
.
6. По формуле (16) рассчитывается
др( )
j
T
в
j
-м приближении.
7. Проверяется условие
н
др др( )
,
j
Т
Т Т Т
∆ = −
≤ ε
0, 0001
Т
ε =
с. (22)
Если условие выполняется, то
0 ( )
j
a a
=
,
0 ( )
j
i
i
=
и осуществляется
переход к п. 8, в противном случае рассчитывается новое значение
большой полуоси орбиты