Вестник СПбГУ. Прикладная математика. Информатика... 2024. Т. 20. Вып. 2
УДК 519.71, 521.1 MSC 70Q05
Управление орбитальным движением космического аппарата вблизи точки либрации с помощью изменения отражательной способности солнечного паруса*
А. С. Шмыров, В. А. Шмыров, Д. В. Шиманчук
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Шмыров А. С., Шмыров В. А., Шиманчук Д. В. Управление орбитальным движением космического аппарата вблизи точки либрации с помощью изменения отражательной способности солнечного паруса // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2024. Т. 20. Вып. 2. С. 298-306. https://doi.org/10.21638/spbu10.2024.214
Рассматривается стратегия удержания космического аппарата в окрестности неустойчивой точки Лагранжа (точки либрации) системы Солнце — Земля. Эта стратегия подразумевает использование специального управления, действующего коллинеарно линии Солнце — Земля и обеспечивающего устойчивость по Ляпунову орбитального движения. Реализация такого управления предлагается с помощью солнечного паруса. Обсуждаются преимущества и трудности данной стратегии.
Ключевые слова: задача трех тел, точка либрации, солнечный парус, управление, удержание, машинное обучение.
1. Введение. Несколько последних десятилетий исследователи в области небесной механики и управления активно разрабатывают проекты, связанные с колли-неарными точками либрации (точками Лагранжа) системы Солнце — Земля и системы Земля — Луна. Список таких проектов, как реализованных, так и планируемых, весьма внушителен. Первый проект ISEE-3 был запущен NASA в 1978 г. В 1994 г. запущена станция наблюдения Солнца SOHO (NASA, ESA). В рамках проекта Genesis (2005 г.) реализовано интересное маневрирование между окрестностями точек либрации L1 и L2. Космический аппарат Чанъэ-2 CNSA в 2011 г. совершил маневры в окрестности точки L2. Проект Спектр-РГ Роскосмоса и DLR был реализован в 2019 г. и посвящен построению схемы карты Вселенной в рентгеновском диапазоне.
В последние годы имеются практические достижения космической навигации, применяющие солнечный парус. Исследования космических полетов с солнечным парусом ведутся давно. Можно сказать, что разработана фундаментальная теоретическая база для полетов. Настоящие трудности обусловлены технической реализацией размещения солнечного паруса в космическом пространстве и управления им для маневрирования. Также к актуальным проблемам для современных полетов такого рода относится разработка подходящих систем управления, которые, с одной стороны, будут эффективно работать, а с другой, — просты и естественны. Потому следует отметить монографии [1, 2], содержащие фундаментальные исследования в этой области. Полетам с солнечным парусом, связанным с точками либрации, посвящено
* Работа выполнена при финансовой поддержке Санкт-Петербургского государственного университета (проект № 95438429).
© Санкт-Петербургский государственный университет, 2024
немало работ. В [3] проведены исследования в рамках линеаризованных уравнений круговой ограниченной задачи трех тел. В [4] предлагается с помощью сил светового давления управлять взаимным расстоянием группировки космических аппаратов в окрестности точки либрации ¿2. В [5] изучена задача изменения площади паруса и его ориентации для удержания в окрестности условно-периодической орбиты в окрестности точки либрации. В 2020-2023 гг. авторами настоящей работы были опубликованы статьи [6-8], посвященные решению задачи управления солнечным парусом в окрестности точки либрации, спецификой которой является использование специальной функции фазовых переменных, где рост модуля свидетельствует об уходе из окрестности точки либрации.
Коллинеарные точки неустойчивые [9], и, таким образом, задача управления нужна не только для маневрирования, но и для удержания космического аппарата (КА) в этих областях пространства.
В представленной работе изучается орбитальное движение в окрестности кол-линеарной точки либрации Ь\ системы Солнце — Земля. Здесь оказывается удобно применить специальную нелинейную аппроксимацию ограниченной круговой задачи трех тел, так называемую модель Хилла. Точки Лагранжа, или точки либрации, являются модельными понятиями круговой ограниченной задачи трех тел и ее частными решениями (рис. 1). При этом подразумевается, что система трех тел, т. е. КА и два массивных небесных тела, сохраняет свою геометрическую конфигурацию в течение времени. Будем рассматривать коллинеарную точку Лагранжа Ь\ системы Солнце — Земля. Эта точка является неустойчивой, и для длительного пребывания в ее окрестности требуется решать задачу управления. Окрестность Ь\ — притягательное место для реализации многих космических проектов.
ф Земля
О Солнце
¿3
Рис. 1. Положения точек либрации
Опубликовано много научных работ, посвященных движению в окрестности точек либрации, в частности обзор [10]. Обычно стратегия полета в окрестность точки либрации заключается в полете с околоземной орбиты на заранее просчитанную орби-
ту около точки либрации, так называемую гало-орбиту. На таких орбитах КА может находиться в неуправляемом режиме несколько месяцев с последующей коррекцией орбиты в виде малых импульсов. В этой статье рассматривается другая задача, когда подразумевается реализация непрерывного управления для удержания КА.
Перспективно использовать солнечный парус для решения такой задачи [1, 2]. Ранее нами была опубликована серия работ, посвященная применению солнечного паруса в задаче орбитального маневрирования в окрестности Li за счет изменения ориентации паруса [6-8].
В настоящей работе предполагается, что управляющим параметром является отражательная способность солнечного паруса. Решение дается в виде синтеза, причем управление зависит только от одной фазовой переменной.
2. Математическая модель движения КА. Уравнения неуправляемого движения КА в геоцентрических координатах в рамках круговой задачи трех тел можно записать в виде [9]
- Y'm1x ( Rl(t) — X 1 Л
* = " W + 7m2 (j|i«(i)-x||3 " Д2 Kt)) ' (1)
где X = (Xi, Х2, X3) — геоцентрические координаты КА; y — гравитационная постоянная; mi, m2 — массы Земли и Солнца; R — расстояние между Землей и Солнцем (1 а.е.); l(t) — единичный вектор, направленный по линии Земля — Солнце.
Если положить у = X, то уравнение (1) приводится к гамильтоновой системе с гамильтонианом
m*,y,t)-2yy рн д ^ц/ад-щ R )■
Далее, как показано в работе [11], раскладывая, следуя Хиллу [12], солнечный потенциал по величинам (которые малы в окрестности первой точки либрации), вводя специальную систему единиц и переходя во вращающуюся систему координат с помощью канонической замены переменных с производящей функции поворота
S(X, у, t) = yi(Xi cos wt + X2 sin wt) + y2(—Xi sin wt + X2 cos wt),
получаем аппроксимированные уравнения движения ограниченной задачи трех тел в виде, который принято называть хилловским приближением или уравнениями модели Хилла [10, 11]
x1=y1+x2, У1 = +2xi +У2,
Х2=У2~Х1, У2 = —[j^JT - ~ У1, (2)
здесь X = (Xi,X2,X3) — координаты КА во вращающейся системе координат, y = (У1,У2,Уз) — импульсы. Расстояние между центром масс Земли и точкой либрации (примерно 0.01 а.е.) берется в качестве единицы расстояния. Единица времени берется из соображения, что время обращения Земли вокруг Солнца (1 год) равно 2п единиц времени (w = 1). В этом случае единица скорости составляет 303.14 м/c, а единица ускорения — 5.93 • 10-5 м/с2.
Точка либрации Li в системе (2) является стационарной с координатами
X* = (1, 0,0), у* = (0,1,0).
Система (2) имеет гамильтонов вид с гамильтонианом
тт / \ 1 и 112 3 3 о ||x||
Hix, У) = 2IMI - - 2Ж! + + Х2Ш ~ XíV2'
Уравнения управляемого движения записываются в виде
х1=у1+х2, = +
Х2=У2-Х\, У2 = - Х1 - Ш +м2, (3)
3жз .
хз=Уз, Уз = -рург - Х3 + из-
В работе [10] был предложен закон управления
ui(xi) = a(xi — 1), а < —9,
ио = 0, (4)
из = 0.
Из (4) видно, что управляющее воздействие реализуется по линии Солнце — Земля. Таким образом, управление в системе (3) будет входить только при yi. В работах [11, 13] показано, что управлением ui(x) = a(xi — 1) при а < —9 обеспечивается устойчивость по Ляпунову стационарного решения (x* , y*) системы (2). Действительно, система (3) с управлением (4) имеет гамильтонову форму. Гамильтониан записывается следующим образом:
1 11 ..о 3 3 о |Ы| (X . .о
В линеаризованном случае гамильтониан H*(x,y) будет иметь вид
-1' \2 , -1/ i \2 , 3 2 /9 + aY ,n2 . 1 2 , о 2
Н\п(х,У) = ^(У2-х1У + -(у1+Х2Г + -х4- j(xi 1) +-yi + 2xi
откуда видно, что при а < —9 квадратичная форма H*n(x, y) положительно определенная. Полная производная от гамильтониана H*(x,y) равна нулю, и из положительной определенности квадратичной формы H*n(x, y) следует, что гамильтониан H*(x,y) будет функцией Ляпунова системы (3), (4) в некоторой окрестности точки либрации в фазовом пространстве. Таким образом, управление ui(xi) обеспечивает устойчивость по Ляпунову стационарного решения (x*,y*) управляемой системы (3), (4).
Приведем уравнения движения в случае, когда управление реализуется с помощью солнечного паруса [6-8]:
¿1 = ух +х2, yi = -ррг +2ж1 +У2- '¿кЩ- cos3 <р — (1 - Щf cosy, Х2 =У2 - XI, У2 = -руз -Х2-У1- 2к1— cos3 ysiny, ¿3=2/3, Уз = -рут - хз,
где m — масса паруса; S — его площадь; p — солнечное давление; k — коэффициент отражательной способности, k G [0,1]; у — угол ориентации паруса (рис. 2). При k =1 парус имеет идеальную отражательную поверхность, при k = 0 обладает абсолютно черной поверхностью.
Рис. 2. Модель солнечного паруса
Управляемое движение с солнечным парусом в случае, когда поверхность паруса перпендикулярна солнечным лучам (угол ^ = 0), описывается уравнениями
3Ж2
Х2=У2~Х1, У2 = —щр ~ Х2 ~ У1,
3хз
(5)
жз = Уз,
2/з =
- х3
в соответствии с законом управления (4). Коэффициент отражательной способности к будем рассматривать как элемент управления. Далее управление, с помощью которого обеспечивается устойчивость по Ляпунову, реализуется через давление солнечных лучей. Далее, будем считать, что отражательная способность паруса — это управление, реализующее закон ^(ж) = а(ж — 1) при а < -9. В связи с тем, что давление солнечных лучей направлено от Солнца, стационарная точка смещается, и, согласно уравнениям (3)-(5), закон управления запишется так:
— ( \ ( л ^
и\{х) = а{х\ — 1 — о) = —к---,
т т
(6)
где константа 5 > 0 связана с фотогравитационной точкой либрации при коэффициенте отражения к = 0.5.
3. Численное моделирование. Рассмотрим движение КА с солнечным парусом с массой т = 5 кг и площадью паруса $ = 32 м2. Приведем начальные данные для траектории движения КА, представленной на рис. 3:
ж(0) = (1.093,0.001,0), у (0) = (0,1.089,0).
Движение происходит в плоскости эклиптики. Рисунок 3 иллюстрирует управляемое движение КА в окрестности фотогравитационной точки либрации, на рисунке показана явная картина устойчивости по Ляпунову, временной интервал представленного графика — порядка 20 месяцев, а коэффициент в управлении а = —12. Управление (6) обеспечивает устойчивость по Ляпунову движения КА в окрестности фотогравитационной точки либрации. При определении фотогравитационной точки либрации коэффициент к = 0.5, тогда для рассматриваемого примера параметр
3
6 = 0.089116540. Координаты фотогравитационной точки либрации для этого примера имеют вид
ж* = (1.089116540,0,0), у* = (0,1.089116540,0).
а
I ! -1 1 ■! ! 1 |
1 \2 1 \4 \ 6 / \ 8 1 \ 1 г
б
Рис. 4- Графики управления (а) и управляющего параметра (б)
На рис. 4, а приведен график функции управления на временном промежутке порядка 20 месяцев.
На рис. 4, б показана зависимость коэффициента отражения k от времени. Значения коэффициента отражения в представленном моделировании не выходят за границы отрезка [0,1]. Однако при других начальных данных, удаленных от точки либрации более чем на несколько тысяч километров, уже наблюдается выход величины коэффициента отражения из отрезка [0,1]. Таким образом, возникает задача управляемости или построения окрестности Li в фазовом пространстве, в которой задача удержания КА решается с помощью солнечного паруса с изменяющейся отражательной способностью и заданной парусностью. Также, одним из перспективных подходов для построения таких управлений являются методы машинного обучения.
4. Заключение. В настоящей работе рассмотрена возможность стабилизации орбитального движения в окрестности коллинеарной точки либрации Li системы Солнце — Земля с помощью сил светового давления. При этом отказались от управления пространственной ориентацией солнечного паруса, поскольку на практике реализация такого управления технически весьма сложна. В описываемой постановке поверхность солнечного паруса расположена перпендикулярно направлению солнечных лучей, а управлением является отражающая способность солнечного паруса.
Основной вывод из численного моделирования следующий: КА с парусностью 64 см2/г можно удержать, изменяя отражательную способность, в окрестности фотогравитационной точки либрации на расстоянии несколько тысяч километров от этой точки. Положение самой фотогравитационной точки определяется параметром задающим ее отклонение от точки либрации Li.
Литература
1. Поляхова Е. Н. Космический полет с солнечным парусом. М.: Книжный дом «ЛИБРОКОМ», 2011. 320 с.
2. McInnes C. R. Solar sailing // Technology, dynamics and mission applications. Berlin; London: Springer, 1999. 296 p. (Springer-Praxis Series in Space Science and Technology)
3. Лукьянов С. С. Управление движением космического аппарата в окрестности коллинеарной точки либрации круговой задачи трех тел с помощью сил светового давления // Космические исследования. 1981. Т. 19. № 4. С. 518-527.
4. Назиров Р. Р., Эйсмонт Н. А., Чумаченко Е. Н., Данхэм Д. У., Логашина И. В., Фе-доренко А. Н. Управление группировкой космических аппаратов в окрестности солнечно-земных коллинеарных точек либрации с помощью солнечного паруса // Вестник машиностроения. 2013. № 2. С. 43-46.
5. Bookless J., McInnes C. Control of Lagrange point orbits using solar sail propulsion // Acta Astronautica. 2008. Vol. 62. P. 159-176.
6. Shymanchuk D. V., Shmyrov A. S., Shmyrov V. A. Controlled motion of a solar sail in the vicinity of a collinear libration point // Astronomy Letters. 2020. Vol. 46. N 3. P. 185-192.
7. Shymanchuk D. V., Shmyrov A. S., Shmyrov V. A. Orbital maneuvering in the vicinity of collinear libration points using light pressure forces // Astronomy Letters. 2021. Vol. 47. N 10. P. 710-719.
8. Shymanchuk D. V., Shmyrov A. S., Shmyrov V. A. Stabilization of the orbital motion of a solar sail in the vicinity of a collinear libration point through a change in reflectivity // Astronomy Letters. 2022. Vol. 48. N 8. P. 469-477.
9. Маркеев А. П. Точки либрации в небесной механике и космодинамике. М.: Наука, 1978. 312 с.
10. Gomez G., Llibre J., Martinez R., Simo C. Dynamics and mission design near libration points. Vol. 1. Fundamentals: The case of collinear libration points. Singapore: World Scientific, 2001. 462 p.
11. Шмыров В. А. Стабилизация управляемого орбитального движения космического аппарата в окрестности коллинеарной точки либрации Li // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. Сер. 10. 2005. Вып. 1-2. С. 192-198.
12. Hill G. W. Researches in the lunar theory // American Journal of Mathematics. 1878. Vol. 1. P. 5-26.
13. Shmyrov A. S., Shmyrov V. A. On controllability region of orbital motion near Li // Applied Mathematical Sciences. 2015. Vol. 9. P. 7229—7236.
Статья поступила в редакцию 27 января 2024 г. Статья принята к печати 12 марта 2024 г.
Контактная информация:
Шмыров Александр Сергеевич — д-р физ.-мат. наук, проф.; a.shmyrov@sbpu.ru Шмыров Василий Александрович — канд. физ.-мат. наук, доц.; v.shmyrov@spbu.ru Шиманчук Дмитрий Викторович — канд. физ.-мат. наук, доц.; d.shimanchuk@spbu.ru
Controlling the orbital motion of the spacecraft near the libration point by changing the solar sail reflectivity*
A. S. Shmyrov, V. A. Shmyrov, D. V. Shymanchuk
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Shmyrov A. S., Shmyrov V. A., Shymanchuk D. V. Controlling the orbital motion of the spacecraft near the libration point by changing the solar sail reflectivity. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2024, vol. 20, iss. 2, pp. 298-306. https://doi.org/10.21638/spbu10.2024.214 (In Russian)
The article considers the strategy of keeping the spacecraft in the vicinity of the unstable Lagrange point (libration point) of the Sun — Earth system. This strategy involves the use of special management. This control acts collinearly to the Sun — Earth line and ensures Lyapunov stability of the orbital motion. The implementation of such control is proposed using a solar sail. The advantages and difficulties of such a strategy are discussed.
Keywords: three bodies problem, libration point, solar sail, control, keeping, machine learning.
References
1. Polyakhova E. N. Kosmicheskiy polet s solnechnym parusom [Space flight with solar sail]. Moscow, Book House "LIBROKOM" Publ., 2011, 320 p. (In Russian)
2. McInnes C. R. Solar sailing Technology, dynamics and mission applications. Berlin, London, Springer, 1999, 296 p. (Springer-Praxis Series in Space Science and Technology)
3. Lukyanov S. S. Upravleniye dvizheniyem kosmicheskogo apparata v okrestnosti kollinearnoy tochki libratsii krugovoy zadachi trekh tel s pomoshch'yu sil svetovogo davleniya [Control of the motion of a spacecraft in the vicinity of the collinear libration point of a circular three-body problem using light pressure forces]. Space Research, 1981, vol. 19, no. 4, pp. 518—527. (In Russian)
4. Nazirov R. R., Eismont N. A., Chumachenko E. N., Danhem D. U., Logashina I. V., Fedo-renko A. N. Upravleniye gruppirovkoy kosmicheskikh apparatov v okrestnosti solnechno-zemnykh kollinearnykh tochek libratsii s pomoshch'yu solnechnogo parusa [Control of a grouping of spacecraft in the vicinity of the solar-terrestrial collinear libration points using a solar sail]. Vestnik of Mechanical Engineering, 2013, no. 2, pp. 43—46. (In Russian)
5. Bookless J., McInnes C. Control of Lagrange point orbits using solar sail propulsion. Acta, Astronautica, 2008, vol. 62, pp. 159-176.
6. Shymanchuk D. V., Shmyrov A. S., Shmyrov V. A. Controlled motion of a solar sail in the vicinity of a collinear libration point. Astronomy Letters, 2020, vol. 46, no. 3, pp. 185-192.
* This research was supported by Saint Petersburg State University, project N 95438429.
7. Shymanchuk D. V., Shmyrov A. S., Shmyrov V. A. Orbital maneuvering in the vicinity of collinear libration points using light pressure forces. Astronomy Letters, 2021, vol. 47, no. 10, pp. 710—719.
8. Shymanchuk D. V., Shmyrov A. S., Shmyrov V. A. Stabilization of the orbital motion of a solar sail in the vicinity of a collinear libration point through a change in reflectivity. Astronomy Letters, 2022, vol. 48, no. 8, pp. 469—477.
9. Markeev A. P. Tochki libratsii v nebesnoy mekhanike i kosmodinamike [The libration points in celestial mechanics and space dynamics]. Moscow, Nauka Publ., 1978, 312 p. (In Russian)
10. Gomez G., Llibre J., Martinez R., Simo C. Dynamics and mission design near libration points. Vol. 1. Fundamentals: The case of collinear libration points. Singapore, World Scientific, 2001, 462 p.
11. Shmyrov V. A. Stabilizatsiya upravlyayemogo orbital'nogo dvizheniya kosmicheskogo apparata v okrestnosti kollinearnoy tochki libratsii Li [Stabilization of the controlled orbital motion of a spacecraft in the vicinity of the collinear libration point Li]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes. Series 10, 2005, iss. 1—2, pp. 192—198. (In Russian)
12. Hill G. W. Researches in the lunar theory. American Journal of Mathematics, 1878, vol. 1, pp. 5-26.
13. Shmyrov A. S., Shmyrov V. A. On controllability region of orbital motion near Li. Applied Mathematical Sciences, 2015, vol. 9, pp. 7229—7236.
Received: January 25, 2024.
Accepted: March 12, 2024.
Authors' information:
Alexander S. Shmyrov — Dr. Sci. in Physics and Mathematics, Professor; a.shmyrov@spbu.ru
Vasiliy A. Shmyrov — PhD in Physics and Mathematics, Associate Professor; v.shmyrov@spbu.ru
Dmitry V. Shymanchuk — PhD in Physics and Mathematics, Associate Professor; d.shimanchuk@spbu.ru