close

Вход

Забыли?

вход по аккаунту

?

Определение пространственных стационарных точек типа polar-sitter в круговой задаче трех тел с использованием солнечного паруса.

код для вставкиСкачать
Известия ТулГУ. Технические науки. 2016. Вып. 5
Akinshin Oleg Nikolayevich, candidate of technical sciences, head of department,
rts@cdbae.ru, Russia, Tula, JSC Central Design Bureau of Apparatus Engineering,
Yesikov Dmitry Olegovich, postgraduate, rts@cdbae.ru, Russia, Tula, Tula State
University,
Akinshina Natalia Yuryevna, engineer, bauman@bmstu.ru, Russia, Moscow, Bauman Moscow State Technical University
УДК 629.154.4
ОПРЕДЕЛЕНИЕ ПРОСТРАНСТВЕННЫХ СТАЦИОНАРНЫХ
ТОЧЕК ТИПА POLAR-SITTER В КРУГОВОЙ ЗАДАЧЕ ТРЕХ ТЕЛ
С ИСПОЛЬЗОВАНИЕМ СОЛНЕЧНОГО ПАРУСА
С.М. Кабанов, Г. В. Фридлендер
Рассмотрена задача определения положений равновесия солнечного паруса в
круговой ограниченной задаче трех тел в системе Солнце-Земля-парус. Получены уравнения движения паруса, определены положения равновесия солнечного паруса. Найдено
параметрическое семейство оригинальных решений в виде асимптотического ряда в
области между Солнцем и Землей.
Ключевые слова: солнечный парус, точки либрации, устойчивость, polar – sitter, асимптотический ряд.
Идея полетов в космосе с использованием солнечного паруса возникла в 1920-е годы в России и принадлежит одному из пионеров ракетостроения Фридриху Цандеру, исходившему из того, что частицы солнечного света — фотоны — имеют импульс и передают его любой освещаемой
поверхности, создавая давление. Величину давления солнечного света
впервые измерил русский физик Пётр Лебедев в 1900 году. В наше время в
связи с высоким техническим прогрессом существуют необходимые материалы для технической реализации солнечногопаруса. Первое развёртывание солнечного паруса в космосе было произведено на российском корабле
«Прогресс М-15» 24 февраля 1993 года в рамках проекта «Знамя-2». Японское космическое агентство (JAXA) 21 мая 2010 года запустило ракетуноситель H-IIA, на борту которой находились космический аппарат
IKAROS с солнечным парусом и метеорологический аппарат для изучения
атмосферы Венеры. Также в конце 2014 года специалистами NASA готовится запуск самого большого в истории паруса Sunjammerплощадью
111,5 квадратных метров при весе в 31,7 килограмм. Все это подтверждает
актуальность использования солнечного паруса в виде движетеля космического аппарата.
116
Информатика, вычислительная техника, обработка и защита информации
Одна из первых книг [1] о солнечном парусе была издана в 1992 году и содержала обзоры основных работ по исследованию некеплеровых
орбит в задаче Солнце-парус. Первые работы, посвященные солнечному
парусу в задаче трех тел Солнце-Земля-парус, начали появляться с выходом книги [2], в которой, в частности, были приведены численные расчеты
точек либрации солнечного паруса для малых значений параметра паруса и
их возможное использование в практических целях. В настоящее время
упор делается на исследование гало-орбит относительно классических точек либрации и , как представляющих наибольший интерес для практического применения. В работах [3, 4] рассматривается стабилизация малой периодической орбиты паруса около точек , . В работе [5] исследована возможность импульсного управления парусом в плоскости эклиптики. Также есть статья [6], посвященная управлению парусом в задаче Земля-Луна-парус.
Исследование гало-орбит солнечного паруса в задаче трех тел –
большая и развивающаяся сфера в научном мире. Из свежих работ можно
по изучению практических миссий солнечного паруса в задаче трех тел
можно отметить [7, 8, 9].
В нашей работе не рассматривается динамика паруса, напротив, мы
хотим получить и исследовать статические положения солнечного паруса.
Подобные исследования проведены в работах [2,8,10]. Однако, в нашей работе мы больше делаем упор на поиск аналитического решения. Нам удается найти параметрическое семейство оригинальных решений в виде
асимптотического ряда в области между Солнцем и Землей, что является
большой редкостью для нелинейных уравнений такого типа и, таким образом, отличает от подобных работ [2, 8, 10] по исследованию точек либрации солнечного паруса.
1. Постановка задачи.
В качестве модели солнечного паруса, следуя пионерской работе
[2], возьмем плоскую поверхность, полностью отражающую солнечный
свет без потерь. Считаем, что масса такого паруса сосредоточена в одной
точке. Тогда, согласно сделанным предположениям, сила давления солнечной радиации направлена по внешней нормали к отражающей поверхности и описывается формулой
= ⊕ (⊙ ⋅ ) ,
⊙
!"
где ⊕ ≃ 5 × 10 Н/м – максимально возможное давление со стороны
солнечной радиации на парус, находящейся на расстоянии в одну астрономическую единицу; – площадь активной отражающей поверхности.
Будем рассматривать такой парус в круговой ограниченной задаче
трех тел Солнце-Земля-парус. Для этого введем две ортогональные системы координат, следующим образом:
117
Известия ТулГУ. Технические науки. 2016. Вып. 5
$%&' – инерциальная система координат, центр которой совпадает с
центром масс Солнца. Ось $% лежит в плоскости эклиптики, и в начальный
момент времени направлена на центр масс Земли. Ось $& дополняет систему до правой; (( , ( , () ) – орты соответствующих осей системы *%&'.
$+,- – неинерциальная система координат, совпадающая в начальный момент с системой $+,- и вращающаяся относительно оси $- с угловой скоростью . = /
0⊙
1
, равной угловой скорости обращения Земли во-
круг Солнца. Таким образом, Земля и Солнце всегда находятся на оси
$+;(2 , 2 , 2) ) – орты соответствующих осей системы $+,- (рис. 1).
Рис. 1. Солнечный парус во вращающейся системе координат 3456
2. Вывод уравнений движения
Получим уравнения движения солнечного паруса под действием
фотогравитациооных сил Солнца и Земли. Векторное уравнение движения
паруса массой 7 в инерциальной системе координат $%&' запишется в
виде
:⊙
:⊕
89 = ) 8 + ) < − 8 > + ?
,
(1)
8
⊕
где 8 = %( + &( + '() – радиус-вектор паруса в системе координат
$%&'; ?
= /7 – ускорение, сообщаемое парусу силой солнечного
давления.
Пусть во вращающейся системе координат $+,- = +2 + ,2 + -2)
– радиус" – вектор паруса относительно центра масс Солнца. Причем, очевидно, что для первой производной вектора 8 по времени во вращающейся
системе координат выполнено равенство
118
Информатика, вычислительная техника, обработка и защита информации
8A = A + .2) × ,
(2)
тогда для второй производной по времени справедлива цепочка равенств
B
89 = <A + .2) × > + .2) × <A + .2) × > =
(3)
BC
= 9 + 2.2) × A + .2) × (.2) × ).
Переписывая уравнение (1) с учетом выражения (3), получим векторное уравнение движения во вращающейся системе координат $+,:⊙
9 = −2.2) × A − .2) × (.2) × ) − ) −
(4)
:⊕
− ) ( − ) + ?
.
⊕
Перепишем уравнение (4) в скалярном виде. Для этого воспользуемся следующими векторными равенствами:
.2) × A = .2) × (+A 2 + ,A 2 + -A 2) ) = .+A 2) × 2 + .,A 2) × 2 =
= .+A 2 − .,A 2 ,
.2) × (.2) × ) = . 2) × (+2 − ,2 ) = . +2) × 2 − . ,2) × 2 =
= −. +2 − . ,2 .
Также запишем выражения для векторов и ⊕ :
= +2 + ,2 + -2) ,
⊕ = − = (+ − )2 + ,2 + -2) .
В итоге получим следующую систему скалярных уравнений движения паруса во вращающейся системе координат:
:⊙
:⊕
+9 = 2.,A + . + − ) + + ) ( − +) + ?F
,
⊕
:⊙
:⊕
,9 = −2.+A + . , − ) , − ) , + ?G
,
(5)
⊕
:⊙
:⊕
-9 = − ) - − ) - + ?H
,
⊕
где
⊕ =
< ⋅ > I ,
7 ⊙
= J+ + , + - , ⊕ = J(+ − ) + , + - .
Для удобства работы с уравнениями (5) перейдем к безразмерным
параметрам и переменным. Введем следующую каноническую систему
единиц: за единицу массы примем массу Солнца K⊕ = 1; за единицу расстояния примем расстояние между центрами масс Земли и Солнца = 1;
за единицу времени – время, необходимое Земле для описания дуги в 1 радиан вокруг Солнца. В такой системе единиц: :⊕ = 1 и. = 1. И введем
?
119
Известия ТулГУ. Технические науки. 2016. Вып. 5
безразмерный гравитационный параметр Земли как : = :⊕ /:⊙ . Что касается солнечного давления, то введем безразмерный коэффициент L следующим образом:
⊕ L=
.
(6)
7:⊙ /
Из равенства (6) видно, что безразмерный параметр паруса равен
отношению силы, действующей на парус со стороны давления солнечного
давления, к силе гравитационного притяжения Солнца на расстоянии в одну астрономическую единицу.
Теперь можно переписать уравнения (5) в безразмерном виде.
Здесь, в отличии от уравнений (5), координаты +, ,, - измеряются в терминах астрономической единицы. Обозначения +, ,, - сохраненыдля удобства.
+
:
< ⋅ >
+9 = 2,A + + − ) − ) (+ − 1) + L
IF ,
⊕
,
:
< ⋅ >
,9 = −2+A + , − ) − ) , + L
IG ,
⊕
где
(7)
:
< ⋅ >
-9 = − ) − ) - + L
IH ,
⊕
(8)
= J+ + , + - , ⊕ = J(+ − 1) + , + - .
Уравнения (7) можно переписать в более удобном виде, если ввести
безразмерную потенциальную функцию
1
1 :
M = (+ + , ) + + .
2
⊕
Таким образом система (7) примет вид
NM
< ⋅ >
+9 − 2,A −
=L
IF ,
N+
NM
< ⋅ >
,9 + 2+A −
=L
IG ,
N,
NM
< ⋅ >
-9 −
=L
IH .
N
Или еще короче – в векторной записи
⋅ >
<
9 + 22) × A − O M = L
.
120
(9)
(10)
Информатика, вычислительная техника, обработка и защита информации
3. Определение положений равновесия в плоскости 346
Для определения точек равновесия паруса положим нулю все производные по времени в системе уравнений (7) и получим систему уравнений
+
:
< ⋅ >
+ − ) − ) (+ − 1) + L
IF = 0,
⊕
,
:
< ⋅ >
,− )− ) ,+L
IG = 0,
⊕
(11)
:
< ⋅ >
− ) -+L
IH = 0.
) ⊕
С первого взгляда на систему (11) становится понятно, что определить аналитическое решение таких нелинейных уравнений затруднительно. Поступим следующим образом: будем искать решение в плоскости $+(рис. 2, а), полагая, что G = 0. Тогда второе уравнение системы (11) допускает решение , = 0.Таким образом, система уравнений перепишется в
проекции на плоскость $+
+
:
< ⋅ >
+ − ) − ) (+ − 1) + L
IF = 0,
⊕
(12)
:
< ⋅ >
−
+
L
IH = 0.
)
) ⊕
3.1. Задача Солнце-парус
Даже после такого существенного упрощения аналитическое решение не очевидно. Будем считать, что отношение массы Земли к массе
Солнца мало, то есть : ≪ 1. Попробуем искать решение в виде асимптотического ряда по целым степеням малого параметра :. Вкачестве порождающего решения возьмем решение при : = 0, что соответствует решению
фотогравитациоонной задачи в поле одного притягивающего центра,
Солнца (рис. 2, б). Будем называть такую задачу Солнце-парус. Итак, при
: = 0, с учетомтого, что ⋅ = cos U и
+
F = cos( U + V) = cos U − sin U ,
+
F = sin( U + V) = cos U + sin U ,
перепишем систему уравнений (12) в виде
L cos U
+
(+
++
cos
U
−
sin
U)
=
,
)
)
(13)
Lcos α
(- cos U + + sin U) = ) .
)
121
Известия ТулГУ. Технические науки. 2016. Вып. 5
a
б
Рис. 2. Солнечный парус в плоскости 346: а - с учетом гравитационного
поля Земли; б - без учета поля Земли, то есть при [ = \
Система уравнений (13) допускает точное решение
+ = ]1 − L ^_` U +
)
-=
L ^_` U `b U
+.
1 − L ^_` ) U
L ^_` U `b U
a
)
1 − L ^_` U
122
)
!
L ^_` a U `b U
+ 1
c (1 − L ^_` ) U)
(14)
Известия ТулГУ. Технические науки. 2016. Вып. 5
L cos U
+ :(+ − 1)
(+
++
cos
U
−
sin
U)
=
+
,
)
)
)
⊕
L cos U
- :(cos
U
+
+
sin
U)
=
+ ).
)
) ⊕
(15)
Рис. 4. Схематическое изображение кривых равновесия задачи
Солнце-Земля-парус
Решения системы (15) можно записать в виде асимптотического ряда по целым степеням малого параметра ::
+f (U) = +j (U; L) + :l (+j , -j ) + ⋯,
-f (U) = -j (U; L) + :n (+j , -j ) + ⋯,
(16)
где +j и-j – решение задачи Солнце-парус, определяемое формулами (14).
Асимптотическое решение задачи Солнце-Земля-парус мало отлично от
порождающего решения задачи Солнце-парус ввиду малости параметра
:,: = 3 × 10!" . Также решение (16) в виде ряда неверно описывает поведение кривых равновесия задачи Солнце-Земля-парус в окрестности Земли. Точнее сказать, что ряд расходится внутри сферы действия Земли [17],
так как отношение :/⊕ будет определяющим и им пренебрегать нельзя.
Поэтому будем считать, что ряд (16) сходится к точному решению внутри
области (+ − 1) + - > : /p .
124
Информатика, вычислительная техника, обработка и защита информации
Перед использованием численных методов для определения точного решения вблизи Земли системы Солнце-Земля-парус, решим систему
(15) для частных предельных значений углов U = 0 иU = d/2. При
U = d/2 имеем классическую задачу трех тел, в которой решения лежат на
оси $+; имеем три классические точки Лагранжа [11], представимых в виде
разложения по малому параметру :
+qr
: /)
: /)
= 1 − s t + ⋯ , +qu = 1 + s t + ⋯ ,
3
3
:
+q1 = −1 −
+ ⋯.
12
(17)
При U = 0 иL ≤ 1 решения также лежат на оси $+ и задаются уравнением
(+ − 1) (+ ) + L − 1) = :+ ,
(18)
для которого, после применения методов [11, 12], найдется три корня
:
1
+vqr = J1 − L −
+ ⋯,
1
3<J1 − L − 1>
2 : /
+vqu = 1 + w s t + ⋯ ,
L 2
:
1
+vq1 = −J1 − L −
+ ⋯.
1
3<J1 − L + 1>
(19)
Полученные аналитические выражения (19) позволяют оценить зависимость от безразмерного параметра L сдвига классических точек Лагранжа (рис. 4). Также эти точки наряду с точками Лагранжа дают хорошее
начальное приближение для численного поиска решения в области, где
аналитическое решение в виде ряда (16) несправедливо, то есть в окрестности Земли.
Используя численный метод продолжения по параметру U для решения системы (15), изложенный в [13, 14, 15], построим параметрическое
семейство кривых равновесия паруса в трех областях для различных параметров паруса: за Солнцем (кривые в окрестности точки ) ); между Солнцем и Землей (кривые в окрестности точки ); и за Землей (кривые в окрестности ) (рис. 5). Схематическое изображение кривых равновесия в
трех областях представлено на рис. 6, где синим обозначены решения, к
которым сходится асимптотический ряд (16) (то есть справедливо аналитическое решение); красным – решения, полученные в процессе численного решения, для которых ряд расходится.
Перейдем к исследованию устойчивости в линейном приближении
полученных точек равновесия.
125
Информатика, вычислительная техника, обработка и защита информации
4. Исследование устойчивости точек равновесия в линейном
приближении
Запишем уравнения движения солнечного паруса в плоскости $+- в
виде
+9 = ΩF − yH ,
(20)
-9 = ΩH + yF ,
где для удобства введены следующие обозначения
1 1 + L ^_` ) U :
Ω(+, -) = + +
+
, y(+, -)
2
⊕
(21)
L cos U `b U
=
.
Пусть (z, {) – вариации координат (+, -) относительно точек равновесия (+j , -j ), тогда уравнения в вариациях для системы (1) запишутся в
виде
z9 = z + |{,
(22)
{9 = ^z + B{.
Здесь
∗
∗
∗
= ΩFF
− yFH
, | = Ω∗FH − yHH∗ , ^ = Ω∗HH + yFH
,B
(23)
∗
∗
= ΩFH + yHH ,
где Ω∗~ =
€u 
‚
€€~ (F ,H )
ƒ ƒ
иy~∗ =
€u „
‚
€€~ (F ,H )
ƒ ƒ
.
Перейдем к фазовым координатам. Вводя вектор … = [z{zA {A]ˆ , перепишем систему (22) в матричном виде
(24)
…A = ‰….
Матрица ‰ задается в блочном виде
\ ‹
|
(25)
‰=Š
 , гдеŒ = Š
,
Π\
^ B
а матрицы \ и‹ – нулевая и единичная матрицы соответственно. Характеристический многочлен матрицы ‰ равен
(26)
‘(’) = ’a − tr Œ ⋅ ’ + det Œ,
∗
∗
где tr Œ = + B = ΩFF + ΩHH иdet Œ = B − |^.
Численное исследование для следа и определителя матрицы Œ в
точках равновесия дает цепочку неравенств:
(tr Œ) − 4 det Œ > 0,
(27)
tr Π> 0,
det Œ ≥ 0.
Исходя из неравенств (27), можно заключить, что характеристический многочлен (26) имеет пару действительных корней ’, = ±|.| и пару
чисто мнимых корней ’),a = ±b|š|. Однако, это справедливо не для всех
точек равновесия. В одной предельной точке (‘), находящейся в окрестности точки , определитель матрицы Œ обращается в нуль. Тогда харак127
Известия ТулГУ. Технические науки. 2016. Вып. 5
теристический многочлен имеет только действительные корни: ’, = 0 –
корень второй кратности и ’),a = ±√tr Œ. В таблице приведены значения
координат предельных точек для некоторых значений параметра паруса.
Численный анализ показывает, что первая предельная точка появляется
для параметра паруса L ≈ 0.16.
Таким образом, заключаем, что плоское движение относительно точек равновесия неустойчиво в линейном приближении для любых параметров паруса L ≤ 1 в отличии от точек в плоскости эклиптики [10, 16].
Рис. 6. Схематическое изображение кривых равновесия задачи
Солнце–Земля–парус
Значения предельных точек для некоторых параметров ž
Заключение
В выполненном этапе работы стало ясно, что нужные для практического применения положения равновесия типа polar-sitter существуют. Более того, в области между Солнцем и Землей существуют решения в виде
замкнутых кривых, параметризованных углом установки паруса, мало от128
Информатика, вычислительная техника, обработка и защита информации
личные от равновесных решений задачи Солнце-парус. Последние существуют в аналитическом виде. Таким образом, для общей задачи справедливы решения в виде асимптотического ряда по малому безразмерному гравитационному параметру Земли.
Однако вблизи Земли решение в виде ряда становится несправедливым ввиду сильного влияния планеты, тем самым сила притяжения Земли
становится определяющей. В этой области в процессе численного исследования плоских решений методом продолжения были найдены решения в
окрестности точки также в виде семейства замкнутых параметрических
кривых. Эти решения не характерны задаче Солнце-парус и никак не могут
быть определены аналитически даже в виде асимптотического ряда. Такие
решения обусловлены добавлением силы притяжения Земли, которая в
этой области становится определяющей.
Списоклитературы
1. Wright J., «Space Sailing». Routledge, 1992. 258 p.
2. McInnes C.R. Solar sailing: technology, dynamics and mission applications. Springer Praxis, 1999. P. 214 - 222.
3. Bookless J., McInnes C. Control of lagrange point orbits using solar
sail propulsion. In: Proceedings of 57th International Astronautical Congress,
IAC-05-C1.6.03, 2005.
4. Baoyin H., McInnes C. Solar sail halo orbits at the sun-earth
cal point. Celestial Mechanics and Dynamical Astronomy 94, 2006.
P. 155 - 171.
5. Ming Xu, Shijie Xu Nonlinear dynamical analysis for displaced orbits
above a planet. Celest. Mech. Dyn. Astron. 102, 2008. P. 327 - 352.
6. Keita Tanaka and Jun’ichiro Kawaguchi. A Design of Small Circular
Halo Orbits around the of the Earth-Moon System. AAS 12-539, 2008.
7. Heiligers J., Ceriotti M., McInnes C.R. and Biggs J.D. Design of Optimal Transfers Between North and South Pole-Sitter Orbits. 22nd AAS/AIAA
Space Flight Mechanics Meeting, Charleston, SC, 2012.
8. Jeannette Heiligers, Colin R. McInnes. Agile Solar Sailing in threebody problem: motion between artificial equilibrium points. In: Proceedings of
64th International Astronautical Congress, IAC-13-C1.8.3, 2013.
9. Waters T.J. and McInnes C.R. Periodic Orbits Above the Ecliptic in
the Solar-Sail Restricted Three-Body Problem. Journal of Guidance, Control,
and Dynamics 30 (3), 2007. P. 687-693.
10. Bombardelli C., Pelaez J. On the stability of artificial equilibrium
points in the circular restricted three-body problem. Celest. Mech. Dyn. Astron.
109(1), 2011. P. 13–26.
11. Маркеев А.П. Точки либрации в небесной механике и космодинамике. М.: Наука, 1978. С. 17-24.
129
Известия ТулГУ. Технические науки. 2016. Вып. 5
12. Вигдорович И.И., Алексин В.А. Асимптотические разложения
корней алгебраических уравнений: учебно-методическое пособие. М.:
МГИУ, 2007. 135 с.
13. Kuznetsov Yu.A. Elements of Applied Bifurcation Theory. SpringerVerlag, 1998.
14. D. Roose, B. De Dier and A. Spence. Aspects of continuation software. In: Continuation and Bifurcations: Numerical Techniques and Applications, NATO ASI series, Series C, Vol. 313, Kluwer, 1990. P. 261-268.
15. Yu.A. Kuznetsov, V.V. Levitin, A.R. Skovoroda. Continuation of
stationary solutions to evolution problems in CONTENT. Mathematisch Centrum, 1996.
16. Thomas Waters, Patrick Browne. Using AUTO on the Solar sail
CR3BP. School of Mathematics, NUI Galway, Ireland, 2010.
17. Балк М.Б. Элементы динамики космического полета. М.: Наука,
1965. С. 201-209.
Кабанов Сергей Михайлович, главный инженер-математик отдела разработки специального программного обеспечения, smkabanov@inbox.ru, Россия, Москва, АО
«Системная динамика»,
Фридлендер Григорий Владимирович, инженер-математик отдела разработки
специального программного обеспечения, fridlender.g.v@gmail.com, Россия, Москва, АО
«Системная динамика»
DETERMINATION OF THE SPATIAL STATIONARY POINTS OF TYPE POLAR-SITTER
IN A CIRCULAR THREE-BODY PROBLEM WITH THE USE OF SOLAR SAIL
S.M. Kabanov, G.V. Fridlender
This paper considers the problem of determination of equilibrium position of solar
sail in the circular restricted three – body problem in the system Sun – Earth – sail. Obtained
the equitation of the sails motion, determined the equilibrium position of solar sail. Founded
parameter family of original solutions in form of an asymptotic series in the region between
the Sun and Earth.
Key words: solar sail, libration point, steadiness, polar - sitter, an asymptotic series.
Kabanov Sergej Mihajlovich, chief engineer-mathematician of department of development of the special software, smkabanov@inbox.ru, Russia, Moscow, JSC "Sistemnaya dinamika",
Fridlender Grigorij Vladimirovich, engineer-mathematician of department of development of the special software, fridlender.g.v@gmail.com, Russia, Moscow, JSC "Sistemnaya
dinamika"
130
1/--страниц
Пожаловаться на содержимое документа