Cодержание

1. ВВедение

2. Траектории "прямого" полета ка от земли к луне

3. БИЭЛЛИПТИЧЕСКИЕ траектории полета ка к луне

    в системе земля-луна

    3.1. переходы между орбитами в центральном гравитационном поле

    3.2. Биэллиптический полет к Луне в модельной системе Земля-Луна

4. АНАЛИЗ "ОБХОДНОГО" ПОЛЕТА КА  К ЛУНЕ В СИСТЕМЕ

    ЗЕМЛЯ-ЛУНА-СОЛНЦЕ

    4.1. Схема и основные особенности "обходного" полета

    4.2. Оценка влияния гравитации Солнца на перигейное

           расстояние орбиты КА

     4.3. Оценка влияния гравитации Земли на гашение энергии

             селеноцентрического движения КА

    4.4. Оценка влияния гравитации Земли при эволюции на

           конечную эллиптическую орбиту у Луны

5. ЧИСЛЕННЫЕ ХАРАКТЕРИСТИКИ ТРАЕКТОРИЙ "ОБХОДНОГО"

     ПОЛЕТА К ЛУНЕ

    5.1. Модель расчета траектории

    5.2. Характеристики "обходного" полета к Луне

6. ВЫВОДЫ

7. ЛИТЕРАТУРА

CONTENTS

1.     iNTRODUCTION

2.     TRAJECTORIES OF DIRECT SPACE FLIGHT FROM THE EARTH

     TO THE MOON

3.     TRAJECTORIES OF BIELLIPTICAL SPACE FLIGHT TO THE MOON

      IN THE EARTH-MOON SYSTEM

3.1.          Orbital transfers in central gravitational field

3.2.          Bielliptical space flight to the Moon in model Earth-Moon system

4.     ANALYSIS OF “BYPASS” SPACE FLIGHT TO THE MOON IN THE

     EARTH-MOON-SUN SYSTEM

4.1.          Scheme and principal peculiarities of “bypass” flight

4.2.          Evaluation of the Sun gravity influence on a spacecraft orbit perigee

4.3.          Evaluation of the Earth gravity influence on a spacecraft celenocentric energy extinguishing

      4.4.   Evaluation of the Earth gravity influence  for evolution to final

            celenocentric elliptical orbit

5.     NUMERICAL CHARACTERISTICS OF TRAJECTORY OF “BYPASS”  FLIGHT TO THE MOON

5.1.          model of a trajectory calculation

5.2.          Some characteristics of  “bypass” flight to the Moon

6.     Conclusions

7.     references .

1.     Введение

В работе делается анализ траекторий полета к Луне нового типа, для которых характерны большие (> 1 млн. км) расстояния отлета от Земли, большое (~ 100 сут.) время полета, но и уменьшенные энергетические затраты на торможение у Луны за счет эллиптичности орбиты подлета к Луне. Сначала, в п. 2, дан анализ обычных траекторий "прямого" полета к Луне. Затем, в п. 3, исследованы характеристики биэллиптических траекторий полета к Луне в модельной системе Земля-Луна. Нa основе этого в п. 4 описана схема и выделены особенности построения "обходных" траекторий в реальной системе Земля - Луна - Солнце, даны оценки гравитационных возмущений для этих траекторий, доказывающие реализуемость временного захвата КА Луной. В п. 5 приведена модель численного расчета и представлены основные численные характеристики некоторых полученных траекторий, на которых осущест-вляется захват КА Луной и уменьшается энергетика торможения при переходе на орбиту спутника или при посадке на ее поверхность.

2. Траектории “прямого” полета КА от Земли к Луне

Исследование характеристик космического полета от Земли к Луне важно как для космонавтики, так и для небесной механики. Это явилось стимулом появления множества работ, посвященных данной проблеме. Особенно важными, как с практической, так и с теоретической точек зрения, оказались работы В.А. Егорова [1 – 4 и др.], в которых выполнен очень полный анализ траекторий полета от Земли к Луне.

Важными являются также результаты В.А. Егорова по захвату точки Луной, ибо даже в случае временного захвата точки оказывается возможным пассивное, без затрат топлива создание спутника Луны, а в случае торможения затраты топлива будут меньше, чем при гиперболическом подлете к Луне. Анализ В.Г. Фесенкова [5] и В.А. Егорова [1 –4]  показал, что вопрос о полном захвате Луной точки является достаточно сложным и до конца нерешенным. Для случая движения в рамках круговой ограниченной задачи трех тел оказался невозможным захват Луной КА, запущенного с Земли на первом витке его траектории. Для произвольных траекторий была доказана теорема о невозможности захвата точки лишь для случая достаточно малого отношения масс главных тел (< 10^-4). Однако В.А. Егоров (1957г.) отметил, что возможен временный захват точки Луной и дал идею построения траекторий временного захвата, проходящих вблизи точки либрации L₁. Это было затем подтверждено результатами ряда исследований, см., например, M.C.Davidson [6], B.Себехей [7], E. Belbruno [8]. Правда, эти траектории захвата проходили далеко от Земли.

В связи с выполненным В.А. Егоровым и в большинстве других работ анализом траекторий КА от Земли к Луне следует отметить два обстоятельства. В.А. Егоров [1]  сделал оценку возмущения Солнцем DF_S движения точки в сфере действия Земли, при r < 930 тыс. км, и получил, что это возмущение мало по сравнению с притяжением Земли F_E, DF_S/ F_E  < 0,14. При движении точки в пределах орбиты Луны это возмущение будет еще меньше, DF_S/ F_E  < 0,01. Поэтому В.А. Егоров свой анализ выполнил для движения точки (КА) в рамках системы Земля – Луна. Во-вторых, построенные траектории полета от Земли к Луне лежали, в основном, в пределах орбиты Луны и, если удалялись за ее  пределы, то недалеко.

Выполненный В.А. Егоровым анализ траекторий полета от Земли к Луне в рамках указанной модели – системы Земля – Луна – имел фундаментальный характер и послужил основой определения практически всех траекторий полета от Земли к Луне. Эти траектории имеют ряд общих черт и иногда сейчас называют траекториями “прямого” полета. Для примера на рис. 1 приведена схематично траектория полета к Луне советского КА «Луна-9»,  запущенного 31.01.1966 г. и осуществившего 3.02.1966 г. первую мягкую посадку на Луну [9].

Для этой траектории время полета на Луну составляет ~ 3.5 сут. КА удаляется от Земли на расстояние до ~ 405 тыс. км.  Подлет к Луне  осуществляется по гиперболической селеноцентрической орбите, скорость «на бесконечности» V_¥ » 1 км/c, скорость торможения при посадке DV₂ » 2600 м/c.

На рис. 2 представлена схема траекторий полета советского КА «Луна-10», запущенного 31.03.1966 г. и выведенного впервые на орбиту искусственного спутника Луны (ИСЛ) 3.4.1966 г. [9]. Эта траектория имеет примерно те же характеристики, что и траектория КА «Луна-9», только после торможения получается не нулевая скорость, а эллиптическая для орбиты с периселением H_p = 350 км, апоселением H _a = 1017 км, поэтому скорость торможения меньше, DV₂ » 0,75 км/с.

В таблице 1 в первой колонке (прямой полет) приведены основные характеристики данного класса траекторий полета к Луне. Они  осуществляются в гравитационном поле системы Земля – Луна, расстояния до Земли меняется до ~ 400 тыс. км, время полета до Луны не превышает ~ 5,4 сут. Захвата Луной нет, КА подлетает к Луне по гиперболе со скоростью V_¥ » 0,8 –

Табл. 1.   ОCНОВНЫЕ ХАРАКТЕРИСТИКИ

НЕСКОЛЬКИХ СХЕМ ПОЛЕТА К ЛУНЕ

1,1 км/с. Не считая импульсов коррекции, на траектории сообщается два импульса (в импульсном приближении): при разгоне с орбиты ожидания у Земли и при переходе на орбиту спутника или посадке. Первый импульс меняется мало, DV₁ = 3,126 – 3,224 км/с, при увеличении начальной энергии от минимальной до нулевой параболической. Во второй снизу и в последней строках для данной схемы приведена величина тормозного импульса DV_f для перехода на круговую орбиту спутника высотой 100 км, DV_f » 0,81 – 0,93 км/с для V_¥ = 0,8 – 1,1 км/с. Оптимальным без учета коррекций будет плоский касательный перелет в апогей Лунной орбиты. По аналогии с обычным касательным двухимпульсным переходом в центральном поле такой перелет к Луне сейчас иногда называется также Гомановским перелетом.

Отметим теперь, что если бы был осуществлен захват КА Луной и подлет к Луне проиcходил по эллипсу с апоселением, например, r_a = 80 тыс. км, то скорость в периселении и импульс торможения уменьшатся, DV_f » 0,65 км/с (или ~2,4км/с для посадки), что меньше на ~ 0,2 км/с импульса торможения для обычного «прямого» перелета. Это может дать заметную экономию топлива и увеличение полезной массы. Например, пусть начальная масса перед посадкой m₀ = 1000 кг, скорость истечения c = 2,9 км/с, масса конструкции m₁ ~ 308 кг. Тогда при DV₂ = 2,6 км/с масса топлива m₂ ~ 592 кг, полезная масса m_f = 100 кг. При уменьшении импульса до DV₂ = 2,4 км/с масса топлива ~ 563 кг, и, при прежней конструкции, полезная масса m_f » 129 кг, т.е. больше на ~ 30%.

Поэтому важно для практики и интересно с теоретической точки зрения найти и исследовать траектории полета к Луне с захватом точки на эллиптическую орбиту. Предварительно выполним, аналогично [10 - 12], анализ еще одной схемы полета к Луне в рамках системы Земля –Луна.

3. БИЭЛЛИПТИЧЕСКИЕ ТРАЕКТОРИИ ПОЛЕТА К ЛУНЕ В СИСТЕМЕ ЗЕМЛЯ – ЛУНА

3.1. Межорбитальные переходы в центральном поле

Для прояснения ситуации сначала рассмотрим основные типы плоских переходов между круговыми орбитами в центральном ньютоновском гравитационном поле.

Задача оптимизации перехода между орбитами встала перед учеными на заре космонавтики в связи с анализом проблемы перелета c Земли на Марс и Венеру. Практически одновременно, в 20-х годах XX века, немецкий ученый В. Гоман [13] и советский ученый Ф.А. Цандер [14] показали, что в классе двухимпульсных плоских переходов между круговыми орбитами оптимальным по сумме величин импульсов является перелет по касательному эллипсу, см. рис. 3. Работы Ф.А. Цандера стали известны широкой научной общественности значительно позже, поэтому данный переход обычно называют Гомановским. В случае перелета между эллиптическими орбитами T₀, T_f со свободной ориентацией линий апсид переходная орбита T соединяет перицентр внутренней орбиты, имеющей менее удаленный апоцентр, и апоцентр внешней орбиты. Если r_a0 < r_af, то r_p = r_p0, r_a = r_af.

А.А. Штернфельд затем открыл, что при большом отношении радиусов исходных круговых орбит T₀, T_f указанный «классический» перелет по касательному эллипсу может уступать в энергетике т.н. биэллиптическому перелету, см. рис. 4. В этом случае используются две эллиптические орбиты: сначала КА отлетает с начальной орбиты T₀ по первой переходной ор-бите T₁ на расстояние r_a, большее радиусов обеих исходных орбит (r_a > r₀, r_a > r_f), затем, после сообщения там некото-рого импульса DV₂, летит по второй орбите T₂ к конечной орбите T_f. Такой перелет А.А. Штернфельд получил сначала (в 1934 г.) для задачи спуска с орбиты [15], затем - для перехода между орбитами [16]. В случае перехода между эллиптическими орбитами r_a > max(r_a0, r_af), переходные орбиты T₁ и T₂ соединяют перицентры исходных орбит и удаленный апоцентр: r_p1 = r_p0, r_a1 = r_a2 = r_a, r_p2 = r_pf. А.А. Штернфельд назвал такие перелеты «обходными». В приложении к задаче перелета от Земли к Луне оптимальный прямой перелет к Луне будет моделироваться Гомановским касательным перелетом с околоземной круговой орбиты на орбиту Луны, при учете ее эллиптичности – в ее апогей. При этом второй импульс этого перехода DV₂ будет играть роль скорости на бесконечности (V_¥= - DV₂) для гиперболической орбиты селеноцентрического движения КА.

Рассмотрим характеристики перелета от Земли к Луне при использовании биэллиптических «обходных» траекторий.

3.2. Биэллиптический полет к Луне в модельной системе
Земля - Луна

Для оценки характеристик биэллиптического полета к Луне полагаем, что сначала осуществляется перелет с круговой околоземной орбиты высотой H₀ = 200 км в перигей орбиты Луны в центральном поле притяжения Земли. Последний импульс скорости этого перехода DV₃ физически не сообщается, а станет, по величине, скоростью «на бесконечности», V_¥ = - DV₃, для гиперболической орбиты дальнейшего селеноцентрического подлета КА к Луне. Полагаем, что в периселении этой орбиты происходит переход на круговую орбиту спутника Луны высотой H_f = 100 км, импульс этого маневра торможения

DV_f = (V_¥² + 2 m_M / (R_M + H_f))^1/2 - (m_M / (R_M + H_f))^1/2,                               (3.1)

здесь R_M=1738 км–радиус Луны, m_M=4902,65км³/c² – ее гравитационная постоянная. Максимальное расстояние r_a, достигаемое в точке приложения промежуточного импульса DV₂, будет параметром задачи. В таблице 1, во втором столбце приведены основные характеристики данной схемы полеты для  r_a = (1,5 млн. км; 4 млн. км; ¥). Величина второго импульса DV₂ = (0,27; 0,11; 0) км/с, величина конечного импульса DV_f = (0,69; 0,70; 0,71) км/с, при этом скорость V_¥ уменьшается до (0.25; 0.34; 0.4) км/c. Учитывая, что увеличение  r_a достигается небольшим изменением начального импульса DV₁ и, более того, может быть обеспечено, в основном, пассивно, с помощью гравитационного маневра у Луны, будем сравнивать переходы этого семейства по сумме величин двух последних импульсов:

w_f = DV₂ + DV_f.                                                                                            (3.2)

Получаем, что w_f = (0,96; 0,81; 0,71) км/с, т.е. данный перелет экономичнее лучшего прямого перелета при  r_a > 4 млн. км. Время перелета для данной схемы очень велико, Dt= (90; 350; ¥) сут для данных вариантов, не считая времени селеноцентрического полета.

Таким образом, биэллиптическая схема полета к Луне, в принципе, при большом расстоянии отлета от Земли позволяет уменьшить энергетику полета. Однако, так как при таких больших расстояниях от Земли велико влияние возмущения от Солнца, в чистом виде эта схема полета нереальна, следует рассмотреть движение КА Земля - Луна - Солнце.

4. Анализ «обходного» полета Ка в системе
земля – луна - солнце

4.1. Схема и основные особенности «обходного» полета
от Земли к Луне

Рассматривая идею биэллиптического перелета к Луне в рамках реальной системы Земля – Луна – Солнце, получаем прежнюю схему полета, рис. 4. Сначала КА разгоняется с начальной околоземной орбиты и выходит на сильно вытянутую эллиптическую орбиту (при этом может быть гравитационный маневр у Луны для уменьшения импульса DV₁). После пролета апогея, где, может быть, сообщается промежуточный импульс DV₂, поднимающий, вместе с Солнечными возмущениями, перигей до Лунной орбиты, КА движется к Лунной орбите. Затем осуществляется сближение с Луной, и, если нужно, торможение для перехода на орбиту спутника Луны или для посадки.

Рассматривая в рамках реального силового поля эту схему полета с учетом полученных выше, в п. 3, характеристик модельной схемы полета, отметим некоторые особенности данного «обходного» полета.

Во-первых, интересно понять, можно ли увеличить перигейное расстояние орбиты КА до радиуса лунной орбиты с помощью, в основном, солнечной гравитации. При положительном ответе второй импульс исчезает (или существенно уменьшается), тогда значительно уменьшается энергетика полета (на ~ 0,2 – 0,3 км/с).

Далее, характерным для биэллиптического перелета является существенно меньшая (по сравнению с прямым полетом) скорость V_¥. Встает вопрос, нельзя ли погасить гравитацией Земли энергию селеноцентрического гиперболического движения КА и осуществить захват КА Луной, пусть даже временный. Это еще уменьшит энергетику полета и снимет ряд геометрических ограничений на ориентацию орбиты подлета к Луне. Отметим, что с этой точки зрения выгодно иметь орбиты с минимально возможным апогеем  r_a, это уменьшит скорость V_¥ и время полета.

Сделаем оценки указанных влияний гравитационных возмущений Солнца и Земли.

4.2. Оценка влияния гравитации Солнца на перигей орбиты КА

Для оценки влияния солнечных возмущения на перигей орбиты КА применим результаты М.Л. Лидова [17]. В этом методе определяются возмущения Dq_ij, причем i = 1 и i = 2 соответсвуют линейному и квадратичному членам в разложении возмущающего ускорения, а j = 1 и j = 2 соответствуют неподвижному и переменному положению возмущающего тела.

Аналогично [17, 18], оценим изменение перигея геоцентрической орбиты КА за промежуток времени от отлета с начальной орбиты до подлета к Лунной орбите в рамках первого приближения, как возмущение от Солнца за виток орбиты КА:

Dr_p = - a De₁₁ » (15 / 2) p (m_S / m_E) (a / r_S)³ a e (1 – e²)^1/2 b₃^'.                        (4.1)

Здесь: m_E, m_S – гравитационные параметры Земли и Солнца; r_S – расстояние до Солнца; a, e – элементы орбиты КА;

b₃^'= cos² g sin 2a,                                                                                           (4.2)

g - угол наклона радиус-вектора Солнца к плоскости орбиты КА, a - угол между проекцией на эту плоскость и направлением на перицентр орбиты КА, |b₃^'| £ 1. Получаем, аналогично [17, 18, 12], что для увеличения r_p должно быть sin 2a > 0, 0 < a < p/2, p < a < 3p/2. Преобразуем (4.1), учитывая, что 1- e = r_p / a, полагая для данной орбиты e » 1, 1- e² » 2 r_p / a и принимая в качестве r_p среднее значение

r_p = (2r₀ + Dr_p) / 2 » Dr_p / 2.                                                                           (4.3)

Тогда

Dr_p » ((15 / 2) p (m_S /m_E) b₃^')² a⁷ / r_S⁶,                                                             (4.4)

или

a » [Dr_p r_S⁶ / ((15/2) p (m_S / m_E)b₃^')²]^1/7.                                                          (4.5)

Для оценки значения большой полуоси орбиты КА, приводящей к необходимому поднятию ее перигея, полагаем Dr_p = 500 тыс. км. Тогда a » 0,71 млн. км, r_a = 1,2 млн. км при b₃^' = 1; a » 0,87 млн. км, r_a = 1,5 млн. км при b₃^' = 0,5.

Таким образом, оценка влияния гравитации Солнца на перигей орбиты КА показывает, что отлет КА на расстояние r_a » (1,2 – 1,5) млн. км при подходящей ориентации Солнца относительно линии апсид обиты КА позволяет осуществить пассивное увеличение перигейного расстояния орбиты КА до радиуса лунной орбиты. Численные расчеты движения КА соответствуют эти оценкам.

Заметим, что оценка дополнительного изменения перигея за счет члена D₁₂a дает существенно меньшую величину, чем (4.4).

4.3. Оценка влияния гравитации Земли на гашение энергии селеноцентрического гиперболического движения КА

Для оценки влияния гравитации Земли на гашение энергии селеноцентрического гиперболического движения КА при его подлете к Луне рассмотрим простейшую модель радиального движения КА. На рис. 5

приведена схема этого движения КА к Луне. Здесь E, M – центры масс Земли,

Луны, EM  = r_M – радиус-вектор Луны относительно Земли, P – положение КА, MP = r, V – селеноцентрические радиус-вектор и вектор скорости КА. Полагаем, что скорость V направлена к Луне и dr / dt = - V < 0.

Возмущающее ускорение Земли

a =a_P–a_M = -(m_E / (r_M + r)²)((r_M + r) / (r_M + r)) + (m_E / r_M²)(r_M / r_M).          (4.6)

Здесь a_p, a_M – ускорения притяжения Землей КА и Луны. Так как a_M > a_P, то возмущение a будет направлено противоположно скорости V и будет тормозить движение КА. Величина возмущения

a = a_M – a_P = m_E / r_M² - m_E / (r_M + r)².                                                           (4.7)

Энергия движения относительно Луны

E = V² / 2 - m_M /r.                                                                                          (4.8)

Ее изменение определяется уравнением:

dE /dt º (dE / dr) (dr / dt) º - V dE / dr = V dV/ dt = -V a,

т.е.

dE / dr = a = m_E / r_M² - m_E / (r_M + r)².                                                            (4.9)

Полагаем r_M = const, тогда

E = (m_E /r_M²) (r - r_E) + m_E / (r_M + r) -  m_E / (r_M + r_*).                                  (4.10)

Здесь r_* - расстояние до Луны при захвате:

E(r_*) = 0.                                                                                                      (4.11)

Разлагая (4.10) в ряд по степеням (r - r_*), получим в квадратичном приближении:

E » a (r - r_*) – a_P (r -r_*)²/ (r_M + r).                                                            (4.12)

Отсюда следует, что

r_* @ r - (r_M + r) [(a/2 a_p) – ((a/2 a_p)² – E / a_P / (r_M + r))^1/2].                        (4.13)

Рассмотрим численный пример. Пусть V_¥ = 0,41 км/с, r = 161 тыс. км, r_M = 390 тыс. км. Тогда (4.13) дает r - r_* = 74 тыс. км, а r_* = 87 тыс. км. При этом a » 1,3×10^-6 км/с². Время достижения нулевой энергии Dt_* » 3,6 сут. «Точный» численный расчет траектории движения КА в системе Земля –Луна – Солнце показал, что захват Луной происходит даже быстрее: r_* » 105 тыс. км, а Dt_* »2,7 сут.

Таким образом, гравитация Земли позволяет для данного класса орбит осуществить гашение гиперболической скорости движения КА относительно Луны и осуществить захват Луной КА.

4.4. Оценка влияния гравитации Земли при эволюции на конечную эллиптическую орбиту у Земли

Оценим влияние гравитации Земли на уменьшение энергии движения КА относительно Луны от нулевой до значения, соответствующего конечной эллиптической орбите. Пусть в конце перелета КА находится у Луны в периселении орбиты, для которой высота периселения H_pf = 100 км. Оскулирующее расстояние в апоселении r_af » 70 тыс. км, т.е. большая полуось a_f » 35 тыс. км, энергия E_f » - 0.07 км²/с².

Согласно [17], в линейном приближении изменение большой полуоси за виток орбиты КА

D₁₂ a = - (3/2) (m_E / m_M) (a / r_M)³ a T b₁₂.                                                     (4.14)

Здесь b₁₂ = e (e + 4) (db₂ / dt - db₆ / dt / 3) + (2 e² + 4 e -1) (db₁ / dt - db₂/dt), b₁ =  x₁² D_M³, b₂ = x₂² D_M³, b₆ = D_M³, D_M = 1 + e_M cosn, x₁ = cosg_E cos a_E, x₂ = cosg_E sin a_E, g_E, a_E – углы, определяющие ориентацию радиус-вектора Земли относительно орбиты КА, T – период орбиты КА. Полагаем приближенно, что эксцентриситет Лунной орбиты e_M = 0, e = 1, dg_E / dt = 0. Тогда

D₁₂a » (15/2) (m_E / m_M) (a / r_M)³ a T b₃^' n_M,                                                 (4.15)

где b₃^’ дано в (4.2), n_M = 2p / T_M – угловая скорость орбитального движения Луны. Переходим от большой полуоси a к энергии E, тогда получим

DE @ (15/2) p m_E (-m_M /2 E r_M)³ n_M b₃^' / (-2 E)^-1/2.                                       (4.16)

Положим, что среднее значение энергии E » DE / 2. Тогда

DE @ sign (b₃^') ((15/2) p m_E (m_M / r_M)³ n_M  |b₃^'|)^2/9.                                      (4.17)

Оценка (4.17) дает DE = - 0,11 км² / с² при b₃^’ = -1, DE = - 0,09 км² / с² при b₃^’ = -0,5. Следовательно, оценка возмущения Землей показывает, что необходимое уменьшение энергии до E_f = - 0,07 км² / с² может быть обеспечено с помощью гравитации Земли примерно за виток орбиты спутника при ее хорошей ориентации. Отметим, что в данном случае, в отличие от п. 4.2, где поднимался перигей и требовалось положительное значение b₃^’, теперь необходима отрицательная величина b₃^’, что отмечено в [12].

Суммируя результаты анализа в п. 4.2, 4.3, 4.4, видим, что есть принципиальная возможность реализации «обходного» полета от Земли к Луне с захватом КА Луной на указанную эллиптическую орбиту. Численные расчеты траекторий подтверждают этот вывод. Сначала, почти одновременно, американские и японские ученые получили «обходные» траектории полета к Луне [10, 11], причем японские ученые реализовали такой полет в проекте Hiten. Недавно европейские ученые также получили траектории данной схемы [12]. Ниже приводятся характеристики некоторых полученных автором траекторий.

5. Численные характеристики траекториЙ «обходного» полета к Луне

5.1. Модель расчета траектории

Траектория полета КА определяется численным интегрированием системы дифференциальных уравнений движения точки в невращающейся геоэкваториальной геоцентрической системе прямоугольных координат в поле притяжения Земли, Луны и Солнца с учетом главной гармоники С₂₀. Эта система уравнений имеет вид:

dr / dt = V, dV / dt = - m_E r / r³ + a_E + a_M + a_S.                                             (5.1)

Здесь r, v – радиус-вектор и вектор скорости КА, r = |r|, a_E , a_M , a_S – возмущающие ускорения, вызванные нецентральностью поля тяготения Земли, притяжением Луны, Солнца [19]. Интегрирование системы (5.1) производится методом ИПМ им. М.В.Келдыша РАН [20], с определением координат Луны и Солнца по JPL – эфемеридам DE403. При этом используется среднее равнодействие и средний геоэкватор стандартной эпохи J2000.0. Расчет ведется с двойной точностью.

5.2. Характеристики «обходного» полета к Луне

Приведем некоторые характеристики двух из полученных автором траекторий «обходного» полета от Земли к Луне. Для одной траектории (T₁) осуществляется довольно быстрый захват КА Луной, для другой (T₂) сравнительно долгий.

На рис. 6 - 10 приведены характеристики первого варианта. Рис. 6 дает проекцию геоцентрической траектории на плоскость XY. Отлет от Земли (r_p0 » 6578 км) происходит 1.1.1997 г. КА отлетает от Земли на расстояние r_max » r_amax » 1,54×10⁶ км. После этого КА летит к Луне, при этом перигейное расстояние под влиянием Солнечной гравитации увеличивается до r_p » 480 тыс. км, большая полуось достигает значения a » 890 тыс.км, затем уменьшается. После подлета КА к Луне на расстояние r » 182 тыс. км в течение ~ 2,6 сут его селеноцентрическая скорость V_¥ уменьшается от ~ 0,4 км/с до 0 при r » 105 тыс. км. Далее, в течение ~ 14 сут. происходит эволюция окололунной эллиптической орбиты. Через ~ 130,5 сут. полета КА приходит в периселений конечной орбиты, для которой r_pf = 1838 км, r_af = 75072 км, а_f = 38455 км, i = 90⁰, W = -45⁰ (относительно геоэкваториальной селеноцентрической системы координат). Полет является полностью пассивным после отлета от Земли. Основные характеристики траекторий этого типа приведены в третьем столбце таблицы 1 (стр. 9).

          Рис. 6. Геоцентрическая траектория Т₁ полета от Земли к Луне

                     в проекции на плоскость XY.

На рис. 6 отмечены точки достижения максимального расстояния от Земли r_max , направление на Солнце в этот момент, моменты достижения максимального перигейного расстояния r_pmax = r_pmax и максимальной (до сближения с Луной) большой полуоси a_max. Здесь же отмечены точки, соответствующие определенной скорости на бесконечности V_¥=V_h.

На рис. 7 приведены для данной траектории Т₁ зависимости от времени расстояния до Земли r(t), большой полуоси a(t), апогейного расстояния r_a=r_a (t), перигейного расстояния r_p=r_p(t) для начальной части геоцентрического полета до сближения с Луной. На кривых даны точки, соответствующие началу захвата, V_¥ = 0.

                На рис. 8 приведена селеноцентрическая траектория полета для второй, окололунной фазы полета, в проекциях на плоскости XY, XZ. Здесь отмечены у Луны точка, соответствующая концу полета, и точки, соответствующие определенной величине скорости на бесконечности V_¥ = V_h. Отмечено также направление на Землю в начале захвата, при V_¥=0.

На рис. 9 дана зависимость от времени константы селеноцентрической энергии h = 2E для этой части траектории. Отмечена точка начала захвата, E = 0. Конечная энергия соответствует приходу в периселений конечной эллиптической окололунной орбиты.

На рис. 10 дано продолжение селеноцентрической траектории полета КА за конечную точку, если бы там не произошло торможение КА. Через ~ 37 сут. КА опять  приобретет положительную энергию, захват будет закончен.

На рис. 11 – 15 приведены некоторые характеристики для другой траектории, T₂, с более длительным захватом. В данном случае отлет происходит 20.1.2001. Максимальное расстояние отлета от Земли r_max » r_amax » 1,3 млн. км, максимальное перигейное расстояние r_pmax » 510 тыс. км, максимальная большая полуось до сближения с Луной a_max » 876 тыс. км, время полета ~ 111 сут., время захвата ~ 28 сут. При уменьшении скорости на бесконечности с ~0.4 км/с до 0 расстояние до Луны уменьшается с ~ 161 тыс. км. до ~ 105 тыс.км.

На рис. 11 дана проекция на плоскость XY геоцентрической траектории T₂ полета КА. Здесь отмечены точки достижения максимального расстояния от Земли r_max , направление на Солнце в этот момент, моменты достижения максимального перигейного расстояния r_{p max}   и максимальной (до сближения с Луной) большой полуоси a_max. Здесь же отмечены точки, соответствующие определенной скорости на бесконечности V_¥=V_h.

Рис. 12. Заключительная часть селеноцентрической траектории Т₂ полета

             к Луне в проекции на плоскость XY.

Рис. 12 представляет селеноцентрический полет для заключительной части траектории, фазы захвата, в проекции на плоскость XY. Точкой отмечено направление на Землю в момент начала захвата КА Луной.

Рис. 13. Заключительная часть селеноцентрической траектории Т₂ полета

             к Луне в проекции на плоскость XZ.

На рис. 13 для данной траектории представлена проекция на плоскость XZ этой заключительной фазы полета к Луне – фазы захвата КА Луной. В конце движения фактически видна половина витка конечной орбиты, ортогональной плоскости XY.

     Рис. 14. Зависимость селеноцентрической константы энергии h =2E от

                    времени для заключительной фазы полета при захвате Луной.

          Рис. 14 представляет зависимость от времени механической энергии селеноцентрического движения КА на заключительной фазе полета при захвате Луной до достижения периселения конечной орбиты. Точкой отмечен момент начала захвата, E(t)=0.

                Рис. 15. Геоцентрическое продолжение траектории Т₂ за конечный

                                периселений

На рис. 15 дано геоцентрическое продолжение этой траектории за конечную точку (при условии отсутствия в ней торможения). Следовательно, данный захват тоже временный, он продолжается ~ 42 сут.

Анализ численных результатов полного моделирования движения КА в системе Земля –Луна – Солнце показывает, что они хорошо соответствуют качественным аналитическим оценкам, выполненным в п. 4.

6. Выводы

Анализ выполненных современной наукой исследований позволяет сделать вывод о существовании траекторий полета от Земли к Луне нового, «обходного» типа. В отличие от обычных траекторий прямого полета с малым временем перелета, они используют отлет от Земли на большое расстояние, ~ 1,5 млн. км, имеют большое время полета, ~ 100 сут, и временный захват Луной КА на эллиптическую орбиту, что заметно сокращает энергетику торможения для перехода на орбиту искусственного спутника Луны или для посадки на ее поверхность. В работе дан анализ особенностей траекторий этого типа, даны оценки влияния гравитационных возмущений, раскрывающие качественную картину эволюции орбиты точки, приводящей к захвату ее Луной. Выполнен численный анализ, позволивший получить ряд траекторий этого типа. Приведены характеристики таких полетов, находящиеся в хорошем соответствии с качественными результатами.

В заключение автор считает своим приятным долгом отметить поддержку, которую оказали автору Ф.И.Баум, С.В.Белоусов, В.Г.Жаров, В.В.Сазонов, В.А.Степаньянц, А.В.Чернов в разработке численного алгоритма. А.В.Чернов оказал также большую помощь в подготовке препринта к печати.

Автор благодарен д-ру M. Bello Mora и безвременно скончавшемуся проф. J.J.Martinez Garcia, которые обратили внимание автора на данную проблему и оказали поддержку в ее анализе.

7. ЛИТЕРАТУРА

1.     Егоров В.А. О некоторых задачах динамики полета к Луне. – Успехи физических наук, т. 63, вып. 1а, 1957, 73 – 117.

2.     Егоров В.А. К вопросу о захвате в ограниченной круговой проблеме трех точек. – ИСЗ, 1959, вып. 3.

3.     Егоров В.А. Пространственная задача достижения Луны. – М.: Наука, 1965.

4.     Егоров В.А., Гусев Л.И. Динамика перелетов между Землей и Луной. – М.: Наука, 1980.

5.     Фесенков В.Г. О возможности захвата при близком прохождении. – Астрономический журнал, N 23, вып. 1, 1946, 45-58.

6.     Davidson M.C. Numerical examples of transition orbits in the restricted three body problem. – Astronaut. Acta, 1964, 10, 308-313.

7.     Себехей, В. Теория орбит. Ограниченная задача трех тел. – Пер. с англ. М.: Наука, 1982.

8.     Belbruno E. Lunar capture orbits, a method of constructing Earth – Moon trajectories and the Lunar GAS Mission. – AIAA Paper No 87-1054. International Electric Propulsion Conference, 1987.

9.     Космонавтика, энциклопедия. Главный редактор В.П. Глушко. М.: Советская энциклопедия, 1985.

10.  Miller, J.K. and Belbruno, E.A. A method for the construction of Lunar transfer trajectory, using ballistic capture. – AAS/AIAA Spaceflight Mechanics Meeting, AAS Paper 91-100, 1991.

11.  Belbruno, E.A., and Miller, J.K. Sun-perturbed Earth – to – Moon transfer with ballistic capture. – Journal of Guidance, Control and Dynamics, vol. 16, No. 4,  July – August 1993, 770-775.

12.  Biesbroek R., Janin G. Ways to the Moon? - ESA Bulletin, 103, August 2000, 92-99.

13.  Hohmann W. Die Erreichbarkeit der Himmelskoerper. – R. Oldenbourg, Munich, 1925.

14.  Цандер Ф.А. Полеты на другие планеты. (Теория межпланетных путешествий), 1924-1925 гг.

 а) Проблема полета при помощи реактивных аппаратов. Межпланетные полеты. Сборник статей». - М.: Оборонгиз, 1961, 285-360.

б) Пионеры ракетной техники. Кибальчич, Циолковский, Цандер, Кондратюк. Избранные труды. - М.: Наука, 1964, 277 – 359.

15.  Штернфельд А.А. Введение в космонавтику. – а) М. – Л. – ОНТИ. 1937. б) М.: Наука, 1974.

16.  Штернфельд А.А. Искусственные спутники. - М.: Гоcтехиздат, 1958.

17.  Лидов М.Л. Эволюция орбит искусственных спутников планет под действием гравитационных возмущения внешних тел. – Искусственные спутники Земли, вып. 8, 1961, 5 –45.

18.  Эльясберг П.Е. Введение в теорию полета искусственных спутников Земли. – М.: Наука, 1965.

19.  Ивашкин В.В., Баум Ф.И. Использование гравитационного маневра у Луны для полета космического аппарата к сближающемуся с Землей астероиду. – Препринт ИПМ им. Келдыша РАН, 2000, N 67.

20.  Степаньянц В.А. , Львов Д.В. Эффективный алгоритм решения системы дифференциальных уравнений движения. - Математическое моделирование. 2000, т.12, вып. 6, 9-14.