Проектно-баллистический анализ развертывания группировки спутников связи на окололунных орбитах

клонением 58 ° и высотой 1500 км. Третий вариант группировки - 6 КА в двух плоскостях на ОИСЛ с наклонением 58 ° и высотой 4260 км.

Рассматривалась следующая схема выведения КА на рабочие ОИСЛ:

• 1.    Ракета-носитель (РН) Союз 2.1б выводит головной блок на круговую опорную орбиту искусственного спутника Земли (ОИСЗ) с наклонением 51.8 ° и высотой 210 км. В состав головного блока входят межорбитальный транспортный аппарат (МТА), выполненный на основе разгонного блока (РБ) «Фрегат», и блок КА. Под МТА понимается модифицированный вариант разгонного блока с увеличенным сроком активного существования, достаточным для обеспечения перелета Земля-Луна, и спо-

• собный осуществлять коррекцию перелетной траектории.

• 2.    МТА обеспечивает разгонный импульс скорости Δ V ₁ и переводит транспортную систему на траекторию полета к Луне.

• 3.    Астроориентация и коррекция перелетной траектории с помощью межорбитального транспортного аппарата. Предполагается одна коррекция перелетной траектории с выдачей корректирующего импульса до 70 м/с.

• 4.    Выдача тормозного импульса скорости в перицентре пролетной гиперболической траектории в окрестности Луны А V 2 с помощью двигательной установки (ДУ) МТА. В результате этого транспортная система переходит на окололунную орбиту фазирования. Наклонение этой орбиты 58 ° , высота перицентра равна высоте рабочей орбиты, а эксцентриситет (вернее, отношение периодов обращений по рабочей и фазирующей орбите) выбирается в соответствии с располагаемым временем на разведение всех КА в орбитальной плоскости.

• 5.    Последовательное отделение КА от межорбитального транспортного аппарата с последующим переводом его на рабочую орбиту. Импульс скорости, который реализуется ДУ КА, на приведение его в рабочую точку обозначим Δ V ₃.

В результате удается разместить космические аппараты связи в одной из рабочих плоскостей. Соответственно, для развертывания всей группировки, расположенной в двух плоскостях, потребуется два пуска РН.

Оценка энергетических затрат для перелета «опорная околоземная орбита - окололунная орбита фазирования» осуществлялась в рамках импульсной аппроксимации активных участков МТА. Первый разгонный импульс скорости А V 1 прикладывается коллинеарно вектору скорости в некоторой точке (характеризующейся аргументом широты u ₀) опорной орбиты c долготой восходящего узла Ω₀. Далее, численно интегрируя уравнения пассивного движения (1), строится траектория, обеспечивающая пролет Луны на некоторой высоте H . При этом наклонение полученной пролетной гиперболической (относительно Луны) траектории обозначим i .

Анализ пассивного участка движения предлагается проводить в геоцентрической экваториальной системе координат.

где

Здесь U – гравитационный потенциал 33 Земли как сжатого сфероида [1];   – радиус- вектор КА относительно Земли; pj - гравитационный параметр j-го небесного тела (индексом 1 обозначено Солнце, индексом 2 – Луна, индексом 3 – Земля); j – радиус-вектор j-го небесного тела; ф - геоцентрическая широта КА.

Таким образом, для некоторого момента времени выдачи первого импульса 1 1 необходимо найти такие значения Q 0 , и 0 и модуля импульса скорости Δ V ₁, чтобы высота пролета Луны H и наклонение пролетной гиперболы i были бы равны высоте и наклонению рабочей ОИСЛ соответственно.

В случае если данные условия выполнены, в перицентре пролетной гиперболы прикладывается тормозной импульс скорости А V 2 , обеспечивающий выход на окололунную орбиту фазирования, эксцентриситет которой e_f выбирается из условия разведения всех КА в одной орбитальной плоскости за время не больше заданного. Расчетную последовательность нахождения e_f а также величины третьего импульса Δ V ₃ представим ниже.

Окололунная орбита фазирования – это эллиптическая орбита, период обращения по которой в целое число раз ( n ) превышает период обращения по рабочей круговой ОИСЛ. Высота апоцентра орбиты фазирования больше, чем высота рабочей круговой ОИСЛ, высота перицентра равна высоте рабочей ОИСЛ.

Тормозной импульс скорости А V 3 , реализуемый ДУ единичного КА для приведения его в рабочую точку ОИСЛ, прикладывается в момент прохождения всего блока КА перицентра орбиты фазирования. Он может быть найден следующим образом.

Орбитальная скорость КА, находящегося на рабочей круговой ОИСЛ:

U ИССЛЕДОВАНИЯ

Havko-

■ ГРАДА

где µ – гравитационный параметр Луны, r – радиус рабочей ОИСЛ.

Большая полуось орбиты фазирования может быть найдена из интеграла энергии:

Тогда период обращения по фазирующей орбите

Угловая дальность между КА, равномерно распределенными в плоскости рабочей орбиты, может быть найдена следующим об-

Задав n и подставляя (2) - (5) в (6), получим уравнение относительно Δ V ₃. Решив его, можно с помощью (3) и (4) также найти большую полуось орбиты фазирования а_f и период обращения по фазирующей орбите Т . Тогда время разведения всех аппаратов в орбитальной плоскости найдем следующим образом:

/_ =7„|\lll I). (7)

Число витков n необходимо выбирать таким, чтобы суммарное время разведения КА в орбитальной плоскости не превышало допустимого.

Эксцентриситет орбиты фазирования

разом:

где N _КА - количество КА в орбитальной плоскости.

В табл. 1 представлены результаты расчета орбит фазирования для трех рассматриваемых высот рабочих ОИСЛ. В расчетах

Тогда для обеспечения фазирования КА должно выполняться условие где n – число витков, совершаемых по орбите фазирования блоком КА вместе с МТА между поочередным отделением КА.

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

Для примера в табл. 2 приведены результаты расчета импульса А V 3 , а также суммарного времени разведения КА для различных отношений периодов обращения КА по фазирующей и рабочей орбитам.

Итак, краевая задача перелета «опорная околоземная орбита - окололунная орбита фа-

Таблица 1

Результаты расчета орбит фазирования для трех рассматриваемых высот рабочих ОИСЛ ( Н _ОИСЛ)

Н ОИСЛ , км	N КА	j	n	А ' , м/с	Т ,, мин	Т Σ	e f
1000	5	72	16	5,507	216,89	9,64	0,008247
1500	4	90	17	5,945	279,55	9,901	0,009685
4260	3	120	10	9,723	717,80	9,969	0,021623

Таблица 2

Импульс ΔV3, и суммарное время разведения КА ТS

n Вариант группировки КА НОИСЛ 1000 км, NКА 5 НОИСЛ 1500 км, NКА 4 НОИСЛ = 4260 км, NКА = 3 А V3, м/с ТΣ А V, м/с ТΣ А V3, м/с ТΣ 1 74,509 0,714 82,298 0,717 75,723 1,286 2 40,579 1,309 45,621 1,291 43,122 2,251 4 21,246 2,499 24,136 2,439 23,193 4,181 8 10,881 4,879 12,431 4,735 12,057 8,04 16 5,507 9,64 6,311 9,327 6,151 15,758 32 2,771 19,16 3,18 18,51 3,107 31,195 зирования» имеет три граничных условия: высота перицентра орбиты фазирования должна быть равна высоте рабочей ОИСЛ, наклонение орбиты фазирования должно быть равно наклонению рабочей ОИСЛ, эксцентриситет орбиты фазирования должен быть равен ef .

Выбираемыми параметрами краевой задачи являются: момент времени подачи первого разгонного импульса 1 1 , величина первого разгонного импульса А V 1 , долгота восходящего узла опорной околоземной орбиты Q 0 , аргумент широты КА в момент подачи первого импульса u ₀, величина второго тормозного импульса скорости, прикладываемого в перицентре пролетной (относительно Луны) гиперболической траектории А V 2 - всего пять параметров. Очевидно, что существует множество решений краевой задачи перелета, и из них необходимо выбрать наиболее оптимальное. В качестве минимизируемого функционала задачи будем рассматривать сумму импульсов скорости:

• -    - ^А V , -^А V 1 + ^А V 2 .

Таким образом, задача перелета сведена к задаче поиска условного экстремума (минимума).

На рис. 1 представлена зависимость отлетного импульса Δ V ₁ от момента времени t ₁ – фактически от даты старта. Зависимость построена для эпохи начала 2020 года в диапазоне одного лунного месяца для трех случаев высоты перицентра орбиты фазирования.

Из данных, представленных на рис. 1, видим, что первый отлетный импульс для полета к Луне крайне слабо зависит от высоты перицентра орбиты фазирования (высоты рабочей ОИСЗ).

На рис. 2 представлена зависимость тормозного импульса А V ₂, прикладываемого в перицентре пролетной гиперболы в окрестности Луны, от t ₁ – фактически от даты старта РН.

На рис. 3 представлена зависимость характеристической скорости маневра перелета на окололунную орбиту фазирования Δ V _Σ от даты старта.

Данные рис. 3 показывают, что наилучшей датой старта с точки зрения минимума характеристической скорости перелета является 12 февраля 2020 года. Также можно видеть, что высота рабочей ОИСЛ достаточно сильно влияет на характеристическую скорость.

А V _Е для оптимальной даты старта составляет:

• -    для Н = 1000 км - 3,867656 км/с;

• -    для Н = 1500 км - 3,837065 км/с;

• -    для Н = 4260 км - 3,761112 км/с.

Для оценки основных массовых характеристик космической транспортной системы были использованы следующие исходные данные:

• -    масса головного блока, выводимого РН на опорную околоземную орбиту: 8200 кг;

• -    конечная масса МТА, выполненного на основе РБ «Фрегат»: 1275 кг;

  
  
  

  Рис. 1. Зависимость отлетного импульса Δ V ₁ от даты старта

  
  

U ИССЛЕДОВАНИЯ

Havko-

Ограда

• -    максимальная масса заправляемого топлива МТА: 6550 кг;

• -    удельный импульс маршевой ДУ МТА: 333,2 с;

• -    удельный импульс ДУ КА: 255 с.

В табл. 3 представлены основные массовые характеристики рассматриваемой транспортной системы при использовании ОИСЛ различной высоты (НОИСЛ). В таблице исполь- зуются следующие обозначения: VХАР – запас характеристической скорости МТА с учетом затрат на коррекцию перелетной траектории; МК - масса транспортной системы, доставленной на орбиту фазирования; МТ – требуемая масса рабочего топлива МТА; МБ - масса полезной нагрузки МТА (масса блока КА, доставленного на орбиту фазирования).

■ Н=1500 км

▲ 11 Н=4260 км со си си си си си

0J

о:

сг>

Ю си

-е- 00

Н=1000км

0J

-е-

Дата старта

Рис. 2. Зависимость тормозного импульса А V 2 прикладываемого в перицентре пролетной гиперболы в окрестности Луны, от даты старта

Рис. 3. Зависимость характеристической скорости маневра перелета на окололунную орбиту фазирования Δ V _Σ от даты старта

Проектно-баллистический анализ развертывания группировки спутников

Таблица 3

Массовые характеристики космической транспортной системы

НО ОИСЛ , км	V ХАР м/с	М Р кг	М к , кг	М Б , кг
1000	3937,656	5742,7	2457,3	1182,3
1500	3912,065	5723,3	2476,6	1201,6
4260	3836,112	5665,1	2534,8	1259,8

Итак, проанализировав три варианта развертывания группировки КА связи на ОИСЛ, оценив затраты характеристической скорости перелета ОИСЗ-ОИСЛ и затраты на перевод единичного КА в рабочую точку ОИСЛ, можно заключить, что суммарная характеристическая скорость перелета ОИСЗ-ОИСЛ для большинства дат старта не превышает 4 км/с, дата старта не сильно влияет на величину характеристической скорости перелета ОИСЗ-ОИСЛ. Более значительно на характеристическую скорость перелета влияет высота рабочей круговой ОИСЛ. В варианте с наибольшей высотой ОИСЗ (4260 км) выигрыш по массе блока составляет около 80 кг.

Затраты характеристической скорости единичного КА на приведение его в рабочую точку в случае, если располагаемое время фазирования достаточно велико (около 10 суток), незначительны – не превышают 10 м/с.