CC BY
59
УДК 531.1+521.1
И. В. Карпасюк
МОДИФИКАЦИЯ АЛГОРИТМА УПРАВЛЕНИЯ КОСМИЧЕСКИМ АППАРАТОМ С СОЛНЕЧНЫМ ПАРУСОМ НА НИЗКОШИРОТНОЙ ОКОЛОКРУГОВОЙ ОРБИТЕ
Введение
В настоящее время околоземное космическое пространство все более интенсивно и широко используется как в целях проведения многочисленных научных исследований, так и для обеспечения самых разнообразных нужд современной цивилизации. В связи с быстрым увеличением числа всевозможных тел искусственного происхождения, движущихся по геоцентрическим орбитам, космическое пространство в окрестности планеты Земля загрязняется отходами (называемыми «космический мусор»), появляющимися в процессе функционирования космических аппаратов.
Постоянно увеличивающаяся скорость прироста объема космического мусора (сосредоточенного большей частью в экваториальной плоскости Земли ввиду ее активной эксплуатации) несет в себе реальную угрозу космическим аппаратам, как уже функционирующим в околоземном пространстве, так и впервые выводимым на геоцентрическую орбиту. Космический мусор уже становился причиной повреждения спутников и орбитальных станций [1]. Неоднократно корректировалась орбита Международной космической станции, чтобы избежать ее столкновения с отслеженными крупными элементами космического мусора. Таким образом, задача уменьшения негативного влияния космического мусора на штатное функционирование космических аппаратов является весьма актуальной.
Решение задачи
Одним из возможных способов решения данной задачи является использование цилиндрических (называемых также широтными) орбит, параллельных экваториальной плоскости Земли и находящихся на определенном расстоянии от нее, над некоторой географической широтой. Для удержания космического аппарата на низкоширотной цилиндрической орбите в течение длительного времени возможно использование солнечного паруса. Солнечный парус представляет собой разновидность двигателей малой тяги. Он может применяться для осуществления долговременного управляющего малого воздействия на элементы орбиты, которого, тем не менее, вполне достаточно, чтобы поддерживать движение космического аппарата на заданной высоте над плоскостью экватора при определенном алгоритме управления ориентацией паруса.
В [2] такой алгоритм управления строился для периода быстрого колебания (суток), в связи с чем направление вектора солнечных лучей I предполагалось постоянным. Однако, ввиду необходимости долговременного характера действия солнечного паруса как управляющего устройства, являющегося источником возмущающей силы, имеющей ненулевую проекцию на нормаль к экваториальной плоскости, изменение направления вектора I необходимо учитывать [3].
Рассмотрим движение космического аппарата с солнечным парусом по околокруговой низкоширотной орбите на высоте Н над плоскостью экватора (предполагая тем самым, что значение Н будет намного меньше величины радиуса орбиты). Выберем декартову систему координат ОХУХ таким образом, чтобы ее центр О совпадал с центром Земли, ось ОХ была направлена в точку весеннего равноденствия, ось ОХ - по нормали к экваториальной плоскости Земли, а ось ОУ составляла с предыдущими осями правую тройку.
Возмущающее ускорение, сообщаемое космическому аппарату солнечным парусом, можно представить в виде [4]:
где и - вектор нормали к теневой стороне плоского паруса, непосредственно реализующий управляющее воздействие на движение космического аппарата; а - коэффициент, зависящий от физических характеристик паруса и являющийся малым параметром. Предполагается, что ||и||2=||/||2=1, откуда
м> = а(и • I)2 и,
(1)
Таким образом, модуль вектора определяется значением параметра а.
В течение периода быстрых колебаний будем считать направление вектора I неизменным, но в течение следующего периода это направление будет уже другим, выбранным с учетом движения Солнца по эклиптике.
Из формулы (1) следует, что меридиональная компонента возмущающего ускорения, соответствующая оси 02 введенной системы координат и обеспечивающая движение спутника по широтной орбите, имеет вид
w3 = а и3 (и • Г)2 = а и3 (u1l1 + u2l2 + u3l3)2, где I = (li, l2, l3) = const в течение периода быстрых колебаний. Введем параметр
и выразим и1, и2 через полярные координаты р, v по формулам:
и1 = р cos v, u2 = р sin v, р e (0, 1], v e [0, 2л). Подставляя (4), (5) и (2) в (3), находим
k = д/1 -р2 рZjcosv + pl2sinv + l3^1-p2
(3)
(4)
(5)
откуда получаем
4k
■ = l1 cos v +12 sin v +13
1 -P2
и
Введем функции:
- + -
4k
, + , • l3y¡1 -P2
l,cosv + l2sin v = ——-----— . ,----
P P-4^7
-4 -p2)+ Jk • V1 -p2
P-V1-P2
(6)
/ ч , , • ,4 -4(1 -р2)+Vk• 4/1 -р2
g(v) = l1 cos v + l2 sin v, д(р) =-^- i - .
рч 1 -р2
График функции д(р) представлен на рисунке. Он построен для предельных значений параметров k, l3.
2
Значение параметра k =
зТз
соответствует максимально возможному значению мери-
диональной составляющей вектора управляющего ускорения от солнечного паруса, обеспечивающему наибольшую высоту цилиндрической орбиты [2]. Реальные значения параметра к, приемлемые с точки зрения эффективности гашения вековых возмущений вектора Лапласа, будут меньше.
Значение параметра /3 »0,3987 соответствует максимально возможному углу наклона плоскости широтной орбиты по отношению к направлению на Солнце (в точках солнцестояний), что составляет примерно 23,5°.
2
2
P
-0,999 р 0,999
График функции д(р)
Очевидно, что в заданных границах изменения значений параметров (у) £ 1, |#(р)| є (0, ¥.
Найдем максимум функции g(v) для известных 11, 12, обозначим его gтах. Используя численные
методы, можно найти отрезки [Рі, Р2], для которых выполняется условие: |д(р)| £ gmax . Значения V,
соответствующие каждому значению р из найденных промежутков, определяются из уравнения (6). Дальнейшие же действия полностью соответствуют алгоритму, описанному в [2]. Этот алгоритм позволяет вычислять компоненты вектора и(ґ) управления ориентацией солнечного паруса для реализации широтного движения космического аппарата, минимизирующего вековые возмущения вектора Лапласа.
Выводы
Таким образом, для каждого периода движения космического аппарата с солнечным парусом по цилиндрической орбите определяется свой закон управления парусом с учетом движения Солнца по эклиптике.
СПИСОК ЛИТЕРАТУРЫ
1. Гаврилов В. Космический мусор: Обломки недавнего прошлого. [Электронный ресурс]. - Заглавие с экрана. Режим доступа: http://www.popmech.ru/article/479-kosmicheskiy-musor.
2. Карпасюк И. В., Шмыров А. С. Управление космическим аппаратом с солнечным парусом на низкоширотной околокруговой орбите // Вестн. С.-Петерб. ун-та. Сер. 1. - 1999. - Вып. 4, № 22. - С. 89-93.
3. Карпасюк И. В. Моделирование возмущенного движения в гравитационном поле // Вестн. С.-Петерб. ун-та. Сер. 1. - 1998. - Вып. 4, № 22. - С. 100-102.
4. Поляхова Е. Н. Космический полет с солнечным парусом: проблемы и перспективы. - М.: Наука, 1986. - 304 с.
Статья поступила в редакцию 15.12.2009
MODIFICATION OF CONTROL ALGORITHM OF A SPACESHIP WITH SOLAR SAIL ON LOW-LATITUDE NEAR-CIRCULAR ORBIT
I. V. Karpasyuk
The problem of algorithm development of the flat solar sail orientation control is considered in the paper. It is set on the spaceship in order to see fulfilled its long-term movements on low-latitude near-circular orbit. The change of sun rays direction, caused by the Sun ecliptic movements, is taken into consideration in the algorithm. Algorithm application can decrease the probability of spaceship with solar sail damage by space debris.
Key words: space debris, spaceship with solar sail, control algorithm of solar sail, latitude orbit, cylindrical orbit, the Sun ecliptic movements taken into consideration.
