Перелеты со встречей продолжительностью не более витка между близкими околокруговыми компланарными орбитами
Автор: Брагазин Александр Федорович, Усков Алексей Викторович
Журнал: Космическая техника и технологии @ktt-energia
Рубрика: Динамика, баллистика, управление движением летательных аппаратов
Статья в выпуске: 3 (30), 2020 года.
Бесплатный доступ
Рассматриваются перелеты со встречей космических аппаратов в классе компланарных непересекающихся околокруговых орбит корабля и орбитальной станции. Продолжительность перелета предполагается ограниченной одним витком. Рассмотрены возможности реализации встречи с использованием оптимального двухимпульсного межорбитального перелета. Для определения единственного свободного параметра перелета - момента его начала, обеспечивающего встречу в заданный момент времени или с заданной скоростью в конце перелета, - получены соответствующие уравнения. Для реализации в алгоритмах наведения предложены оптимальные трехимпульсные программы коррекций для достижения встречи в заданный момент времени с заданной относительной скоростью в момент контакта космических аппаратов. Определен диапазон разностей фаз на начало маневрирования, в котором характеристическая скорость встречи совпадает с минимальной характеристической скоростью межорбитального перелета. Приведены результаты моделирования схем «быстрой» встречи с использованием предложенных программ.
Космический корабль, орбитальная станция, "быстрая" встреча, перелет, программа встречи
Короткий адрес: https://sciup.org/143177936
IDR: 143177936 | DOI: 10.33950/spacetech-2308-7625-2020-3-82-93
Текст научной статьи Перелеты со встречей продолжительностью не более витка между близкими околокруговыми компланарными орбитами
брагазин а.ф.

уСков а.в.
БРАГАЗИН Александр Федорович — кандидат технических наук,
Аппаратура спутниковой навигации, установленная на современных космических аппаратах (КА), значительно повысила оперативность и точность навигационного обеспечения системы управления сближением и обеспечила необходимые условия для реализации так называемых схем «быстрой» встречи КА [1, 2]. Расчетная продолжительность встречи российских кораблей серий «Союз» и «Прогресс» для исходных некомпланарных орбит с момента начала маневрирования до момента встречи в алгоритмах наведения [1, 3] ограничена снизу двумя витками. В работе [4] предложены программы встречи КА для околокруговых компланарных непересекающихся орбит, дополняющие алгоритм наведения [1, 3] для продолжительностей встречи от одного до двух витков. Программы содержат на конечном участке перелета трехимпульсную программу встречи продолжительностью виток, определенную по методу распределения импульсов оптимального межорбитального перелета. Фазирование для достижения необходимой разности фаз на начало трехимпульс-ной программы встречи обеспечивается маневрированием на предшествующем участке полета и достигается разделением импульсов оптимального межорбитального перелета, не приводящим к увеличению расхода топлива.
Программы для продолжительностей встречи не более витка рассматриваются в предлагаемой статье. Назначение продолжительности встречи в пределах витка позволяет рассмотреть возможности использования оптимального двухимпульсного перелета между околокруговыми компланарными орбитами для достижения встречи. Для продолжительности встречи, превышающей продолжительность оптимального двухим-пульсного межорбитального перелета, определены условия существования трех-импульсного перелета с фиксированным значением конечного импульса, обеспечивающего заданные априори конечные условия встречи.
Статья дополняет исследование алгоритмов встречи [4] для околокруговых компланарных непересекающихся орбит. Предложенные оптимальные решения задачи встречи продолжительностью менее двух витков содержат только трансверсальные импульсы одного направления и относятся к так называемым компланарным особым решениям в соответствии с терминологией, предложенной в работе [5].
Представление относительногодвижения ка для компланарных орбит [3]
Вектор относительного положения двух КА в проекциях на оси опорной орбитальной системы координат (ОСК), вращающейся с угловой скоростью ω0 = const, определяется как разность ∆r→ ≡ (x, y) = r→П – r→A, где r→A — вектор положения активного корабля (АК), выполняющего маневры; r→ П — вектор положения пассивного корабля, например, орбитальной станции (ОС). Аналогично определяется вектор ∆v→ относительной скорости КА. Безразмерный вектор относительного положения определяется по формуле ∆r , где R0 — радиус опорной орбиты; безразмерный вектор скорости — ∆v→ по формуле ω . В качестве независимой переменной в линеаризованных дифференциальных уравнениях движения [6] принимается центральный угол ϑ = ω0t, определяющий положение ОСК, являющийся безразмерным аналогом времени, и по отношению к ϑ в статье исполь- зуется терминология, принятая для реального времени. Параметры относительного движения АК и ОС определим по формулам:
c 1 = x – 2 y ′ ; c 2 = 2 y + x ′ ;
c 3 = –3 y – 2 x ′ ; c 4 = y ′
(штрихами помечены производные координат по аргументу ϑ ) и представим решение уравнений относительного движения в виде:
c 1 ( ϑ ) = c 1 ( ϑ 0 ) – 3 c 2 ( ϑ – ϑ 0 );
c 2 ( ϑ ) = c 2 ( ϑ 0 );
c 3 ( ϑ ) = c 3 ( ϑ 0 )cos( ϑ – ϑ 0 ) + c 4 ( ϑ 0 )sin( ϑ – ϑ 0 );
c 4 ( ϑ ) = c 4 ( ϑ 0 )cos( ϑ – ϑ 0 ) – c 3 ( ϑ 0 )sin( ϑ – ϑ 0 ).
Физический смысл параметров относительного движения поясняется в работах [3, 7].
реализация встречи на двухимпульсном оптимальном межорбитальном перелете между близкими околокруговыми компланарными непересекающимися орбитами
Для непересекающихся орбит АК и ОС параметры относительного движения удовлетворяют ограничению:
J = 4 c 2 2 – c 3 2 – c 4 2 ≥ 0, (3) где J — инвариант невозмущенного относительного движения двух КА. Без ограничения общности далее будем полагать, что c 2 > 0.
Оптимальный межорбитальный перелет продолжительностью менее витка содержит два трансверсальных импульса
J ∆ u 0 = ;
1 4(2 С 2 + C8 >
1 (4)
∆u20 = 4(2c2 + c3) [(2c2 + c3)2 + c42] = с2 – ∆u10, если предположить дополнительно, что 2c2 + c3 ≠ 0.
Начало перелета ϑ 1 — произвольный момент времени; момент окончания перелета ϑ 2 = ϑ 1 + τ 0, где
c
τ 0 = π + 2arctg 4 (5)
2 c + c
— продолжительность оптимального перелета, 0 ≤ τ ≤ 2 π .
При 2 c 2 + c 3 = 0 единственный импульс, прикладываемый в начале перелета в точке касания орбит, определяется по формуле
∆ u 1 0 = 0,25(2 c 2 + c 3 2 + c 4 2 ) = c 2 . (6)
Графическое представление характеристик перелета содержится в работе [3].
Характеристики оптимального межорбитального перелета зависят только от параметров c 2 , c 3 , c 4 , определяющих размеры и форму исходных орбит. Встреча по завершении оптимального двухимпульсного перелета достигается при условии, что на момент ϑ 1 разность фаз встречающихся КА удовлетворяет уравнению [3]
c1(ϑ1) – Σn 3(ϑi – ϑ1)∆ui = i = 1 (7)
= c 1 ( ϑ 1 ) – 3 τ 0( ϑ 1 ) ∆ u 2 0 ( ϑ 1 ) = 0.
Связь (7) между параметрами относительного движения определяет положение соединения , когда одновременно с завершением межорбитального перелета достигается встреча КА.
Представим параметры относительного движения на момент ϑ1 как функции аргумента 0 ≤ ϕ ≤ 2π, отсчитываемого от момента прохождения максимального по модулю отклонения по высоте y* траектории относительного движения, и параметра формы траектории e*, определяемого как отношение минимального y** и максимального отклонения по высоте (e* = y**/y*):
c 2 ( ϕ ) = 0,25( y * + y **) = 0,25 y *(1 + e *);
c 3 ( ϕ ) = 0,5( y * – y **)cos ϕ =
= 0,5 y *(1 – e *)cos ϕ ; (8)
c 4 ( ϕ ) = –0,5( y * – y **)sin ϕ =
= –0,5 y *(1 – e *)sin ϕ .
Без ограничения общности положим y * = 1.
Представление о совокупности положений соединения в декартовых координатах ОСК можно получить, выразив координаты как функции аргумента ϕ :
x ( ϕ ) = c 1 *( ϕ ) + 2 c 4 ( ϕ ) = c 1 *( ϕ ) + 2 y ′ ( ϕ ) = c 1 *( ϕ ) – (1 – e *)sin ϕ ;
y(ϕ) = 2c2 + c3(ϕ) = 0,5[1 + e* + (1 – e*)cosϕ], где c1*(ϕ) определяется связью (7).
Общий вид положений соединения в ОСК с началом в центре масс орбитальной станции (ОС) представлен на рис. 1. Параметром кривых является e *. Отметим, что при значениях ϕ = 0 и ϕ = π значение τ 0 = π, а положения соединения принадлежат выделенной синим пунк-
3 π
тиром прямой
x = y , называемой линией 4
хомановских переходов [7].

Рис. 1. Положения соединения на плоскости X0Y ОСК орбитальной станции
Возможности практического использования двухимпульсного оптимального перелета для достижения встречи ограничены, так как единственным свободным параметром является момент ϑ 1 — начало перелета. В известных алгоритмах управления сближением [1, 3, 8] назначаются момент встречи ϑ пр и так называемая прицельная скорость v → пр = ( V пр , 0, 0) на момент контакта, выбираемая исходя из требований навигационного обеспечения и безопасного сближения при подлете АК к ОС. Выбором начала перелета можно обеспечить либо встречу КА в заданный момент времени, либо величину скорости в момент контакта КА в диапазоне, ограниченном апсидальными импульсами оптимального межорбитального перелета.
Будем полагать, что на момент ϑ N разрешенного начала маневрирования заданы параметры относительного движения c 2 , c 3 , c 4 .
Оптимальный двухимпульсный перелет продолжительностью τ 0, заканчивающийся встречей КА, начинается в момент ϑ 1 = ϑ пр – τ 0. Значение ϑ 1 является решением трансцендентного уравнения на интервале ( ϑ N , ϑ пр ):
ϑ пр – ϑ = τ 0( ϑ ). (10)
Рис. 2 иллюстрирует нахождение решения уравнения (10). Аргумент ϑ и функции f ( ϑ ) = ϑ пр – ϑ ; τ 0( ϑ ) выражены в градусах. Решение ϑ 1 для начала оптимального двухимпульсного перелета далее будем обозначать ϑ 1 0 по аналогии с другими характеристиками такого перелета. Оно определяется как абсцисса точки пересечения прямой f ( ϑ ) = ϑ пр – ϑ , выделенной синим цветом, и функции τ 0( ϑ ) на интервале ( ϑ N , ϑ пр ).
как функции начала перелета ϕ , удовлетворяют очевидному ограничению
0,25 y ** ≤ ∆ u 0 ( ϕ ) ≤ 0,25 y *;
1 (12)
0,25 y ** ≤ ∆ u 0 2 ( ϕ ) ≤ 0,25 y *,

Рис. 2. Иллюстрация к решению трансцендентного уравнения (10)
где y ** и y * — отклонения вдоль радиус-вектора в апсидальных точках невозмущенного относительного движения. Если заданная величина прицельной скорости удовлетворяет неравенству
y ** ≤ 4 V пр ≤ y *, (13)
то момент начала перехода с прицельной скоростью ∆ u 2 0 = V пр находится из уравнения:
∆ u 2 0 = c 2 – ∆ u 1 0 = V пр .
С помощью формул (4), (8), получим
уравнение относительно ϕ :
2 e * – (1 + e *) 1 + e *
= (1 – e *)
1 + e *
V пр y * _
- 4 2 nL У * .
cos ϕ .
На рисунке представлены решение ϑ 1 0 уравнения (10) для заданного ϑ пр при фиксированном значении e * = 0,4 и интервал допустимых значений ( ϑ m пр in, ϑ пр ), а также максимальный в предположениях задачи интервал ( ϑ N , ϑ m пр ax) продолжительностью виток. Будем предполагать, что заданное значение ϑ пр ≥ ϑ m пр in. Очевидно, граница ϑ m пр in определяется по формуле
ϑmin = ϑ + τ0(ϑ ), (11) пр N N , где τ0(ϑN) — продолжительность оптимального перелета, начинающегося в момент ϑN, рассчитывается по формуле (5). Разность фаз на момент начала перелета определяется связью (7). Разность фаз на момент ϑN определяется по формуле прогнозирования переменной c1(ϑ).
Параметры, определяющие размер и форму исходных орбит КА, ограничивают возможность назначения произвольного значения V пр при реализации оптимального двухимпульсного межорбитального перелета. В самом деле, значения импульсов, рассматриваемые
Напомним, что ϕ отсчитывается от большей по величине апсиды траектории относительного движения. В соответствии с формулами (8) аргумент ϕ N этой апсиды на момент ϑ N определяется через параметры c 3 ( ϑ N ) и c 4 ( ϑ N ). Начало перелета ϑ 1 нетрудно привязать к ϑ N :
ϑ 1 = ϑ N + ( ϕ – ϕ N ).
Очевидно, момент встречи не может задаваться независимо от V пр и определяется по формуле ϑ пр = ϑ 1 + τ 0( ϑ 1 ), где τ 0( ϑ 1 ) — продолжительность перелета. Разность фаз на момент начала перелета определяется связью (7). Значение c 1 ( ϑ N ) на разрешенное начало маневрирования определяется по формулам прогнозирования невозмущенного движения.
оптимальная трехимпульсная встреча для непересекающихся орбит
Для обеспечения заданной прицельной скорости V пр в заданный момент встречи ϑ пр рассмотрим трехимпульсную программу с трансверсальными импульсами. Для фиксированного начала перелета ϑ 1 всегда существует оптимальный двухимпульсный межорбитальный перелет продолжительностью τ 0( ϑ 1 ), рассчитываемый по формулам (4)–(6). Будем
далее полагать, что заданный момент встречи ϑ пр > ϑ 1 + τ 0( ϑ 1 ), а программа встречи имеет вид
ϑ 1 |
ϑ 2 |
ϑ 3 = ϑ пр |
∆ u 1 |
∆ u 2 |
∆ u 3 = V пр . |
Представим уравнения связи между параметрами движения и импульсами [3] в следующем виде:
c 1 ( ϑ 1 ) – 3( ϑ 2 – ϑ 1 ) ∆ u 2 – 3( ϑ пр – ϑ 1 ) V пр = 0;
c 2 ( ϑ 1 ) – ∆ u 1 – ∆ u 2 – V пр = 0;
c 3 ( ϑ 1 ) + 2 ∆ u 1 + 2 ∆ u 2 cos( ϑ 2 – ϑ 1 ) +
+ 2 V пр cos( ϑ пр – ϑ 1 ) = 0; (15)
c 4 ( ϑ 1 ) + 2 ∆ u 2 sin( ϑ 2 – ϑ 1 ) +
+ 2 V пр sin( ϑ пр – ϑ 1 ) = 0.
Последние три уравнения системы (15) определяют характеристики ∆ u 1 , ∆ u 2 , ϑ 2 межорбитального перелета, начинающегося в момент ϑ = ϑ 1 . Выделим их как отдельную систему:
c 2 ′ – ∆ u 1 – ∆ u 2 = 0;
c 3 ′ + 2 ∆ u 1 + 2 ∆ u 2 cos τ = 0; (16)
c4′ + 2∆u1sinτ = 0, где c2′ = c2′ (ϑ1) = c2(ϑ1) – Vпр;
c 3 ′ = c 3 ′ ( ϑ 1 ) = c 3 ( ϑ 1 ) + 2 V пр cos( ϑ пр – ϑ 1 ); (17)
c 4 ′ = c 4 ′ ( ϑ 1 ) = c 4 ( ϑ 1 ) + 2 V пр sin( ϑ пр – ϑ 1 );
τ = ϑ 2 – ϑ 1 .
Система (16)–(17) определяет оптимальный двухимпульсный перелет с началом в произвольный момент ϑ = ϑ 1 только для непересекающихся компланарных орбит. Представим инвариант относительного движения с параметрами c i ′ ( i = 2, 3, 4) в виде:
J = 4 c 2 ′ 2 – c 3 ′ 2 – c 4 ′ 2 = 4 c 2 2 – c 3 2 –
– c 4 2 – 4 V пр y 0( ϑ пр ) = J 0 – 4 V пр y 0( ϑ пр ), где
J 0 = 4 c 2 2 – c 3 2 – c 4 2 =
= (2c2 + c32 + c42)(2c2 – c32 + c42) = y*y**, y0(ϑпр) = 2c2 + c3cos(ϑпр – ϑ1) +
+ c 4 sin( ϑ пр – ϑ 1 ).
Условие непересечения орбит принимает вид:
J 0 – 4 V пр y 0( ϑ пр ) ≥ 0. (18)
Отсюда, в частности, следует, что оптимальное решение существует даже для исходных орбит, касающихся в момент ϑ = ϑ пр , если c 2 – V пр ≥ 0.
Решением системы (16)–(17) является оптимальный двухимпульсный перелет с начальной орбиты, определяемой параметрами c i ( ϑ 1 ) ( i = 2, 3, 4) движения относительно круговой опорной орбиты, на некоторую промежуточную орбиту, при движении по которой в момент ϑ 3 = ϑ пр достигается цель. На момент ϑ = ϑ 2 параметры относительного движения, определяющие промежуточную орбиту, суть:
c 2 ( ϑ 2 ) = V пр ;
c 3 ( ϑ 2 ) = –2 V пр cos( ϑ пр – ϑ 2 ); (19)
c4(ϑ2) = –2Vпрsin(ϑпр – ϑ2), и значения ординат апсидальных точек траектории относительного движения в орбитальной системе координат:
y = 2 c + c 2 + c 2 = 4 V ;
α 2 3 4 пр yπ = 2c2 – c32 + c42 = 0, т. е. переходная орбита касается орбиты ОС. Траектория перелета с момента ϑ2 до момента ϑ3 = ϑпр является частью траектории хомановского перелета между близкими околокруговыми компланарными орбитами. В самом деле, представим декартовы координаты в ОСК корабля на момент ϑ2 через параметры относительного движения:
x ( ϑ 2 ) = c 1 ( ϑ 2 ) + 2 c 4 ( ϑ 2 ) =
= V пр [3( ϑ пр – ϑ 2 ) – 4sin( ϑ пр – ϑ 2 )]; (20)
y ( ϑ 2 ) = 2 c 2 ( ϑ 2 ) + с 3 ( ϑ 2 ) = 2 V пр [1 – cos( ϑ пр – ϑ 2 )].
Формулы (20) являются параметрическим представлением траектории хомановского перелета на отрезке [ ϑ 2 , ϑ пр ]. Для представления траектории в ОСК ОС в формулах (20) следует изменить знаки x ( ϑ 2 ) и y ( ϑ 2 ) на противоположные.
В уравнениях (15)–(17) варьируемые моменты начала и конца перелета удовлетворяют ограничениям ϑ N ≤ ϑ 1 < ϑ пр и ϑ пр – ϑ N ≤ 2 π . Обозначим ∆ϑ = ϑ пр – ϑ 1 , где ∆ϑ ≥ τ 0( ϑ 1 ).
Если ∆ϑ = τ 0( ϑ 1 ), то трехимпульсный перелет в точку встречи вырождается в двухимпульсный перелет.
Если ∆ϑ = 2 π , то трехимпульсный перелет для встречи, начинающийся в момент ϑ 1 = ϑ N и рассчитанный по уравнениям (16)–(17), иллюстрирует метод распределения (расщепления) импульсов двухимпульсного межорбитального перелета [7] между исходными орбитами АК и ОС. В этом случае первый импульс ∆ u 1 0 межорбитального перелета делится на два импульса программы встречи следующим образом:
∆ u 1 ( ϑ 1 ) = ∆ u 1 0 – V пр ;
∆ u 3 ( ϑ пр ) = ∆ u 3 ( ϑ 1 + 2π) = V пр .
Если τ 0( ϑ 1 ) < ∆ϑ ≤ 2 π , то уравнения (15)–(17) позволяют определить множество трехимпульсных перелетов с заданной прицельной скоростью V пр . Далее рассматриваются два типа перелетов, представляющих практический интерес.
-
1. Трехимпульсные перелеты первого типа. Начало перелета ϑ 1 — варьируемый параметр, конец перелета ϑ пр — фиксированный параметр. Фиксация конца перелета используется, например, для обеспечения приемлемых условий освещенности на случай перехода к ручному управлению, для обеспечения резервирования связи между кораблем и ЦУП на конечном участке встречи.
-
2. Трехимпульсные перелеты второго типа. Начало перелета ϑ 1 — фиксированный параметр, конец перелета ϑ пр — варьируемый параметр. Перелет используется для обеспечения минимальной продолжительности операции с момента выведения.
Будем полагать, что для заданных на момент ϑ = ϑ N параметров относительного движения неравенство (18) справедливо.
Рассмотрим трехимпульсные перелеты первого типа. Найдем решение уравнения (10), определяющего момент ϑ 1 0 начала оптимального двухимпульсного перелета, заканчивающегося встречей в момент ϑ пр . Начало ϑ 1 трехимпульсного перелета назначим на интервале ( ϑ N , ϑ 1 0). Тогда назначенная продолжительность трех-импульсного перелета ϑ пр – ϑ 1 > τ 0( ϑ 1 0).
Первые два импульса трехимпульс-ного перелета составляют оптимальный перелет с исходной орбиты корабля на промежуточную хомановскую орбиту, при движении по которой в момент ϑ = ϑ пр обеспечивается заданная прицельная скорость V пр . Продолжительность τ двухимпульсного перелета на промежуточную орбиту рассчитывается по формуле (5) для параметров относительного движения c ′ i ( ϑ 1 ) ( i = 2, 3, 4).
Для существования трехимпульсного перелета необходимо, чтобы было удовлетворено неравенство:
ϑ 1 + τ ( ϑ 1 ) < ϑ пр . (21)
Покажем, что, если ϑ 1 < ϑ 1 0, то (21) справедливо. Рассмотрим функцию
f(ϑ) = ϑ + τ0(ϑ) = ϑ + π + 2arctg c 4(S)
2 c 2 + c3(9)
где τ 0( ϑ ) — продолжительность оптимального двухимпульсного перелета для исходных орбит. Производная этой функции с помощью (8) представляется в виде:
f ′ ( ϑ )
2 e*
(1 - e *2)cos 9 + (1 + e *2)
и положительна при e * > 0. Поскольку ϑ 1 0 — решение уравнения (10), то
ϑ10 + τ(ϑ10) = ϑпр, и для ϑ1<ϑ10 справедливо неравенство
ϑ 1 + τ 0( ϑ 1 ) < ϑ пр . (22)
Для справедливости неравенства (21) достаточно показать, что
τ ( ϑ 1 ) ≤ τ 0( ϑ 1 ). (23)
Поскольку arctgx — монотонно возрастающая функция при –∞ ≤ x ≤ ∞, преобразуем неравенство (23) с помощью формулы (5) к виду c′4 ≤ c4 2c′2 + c′3 2c2 + c3 , и подставим переменные ci (i = 2, 3, 4), определенные по формулам (17). Используя уравнения (9) и (17), перепишем неравенство (23) в виде у ‘(9,) + 2 KpSWnp — 91) у '(9J
≤ , y ( ϑ 1 ) + 2 V пр [cos( ϑ пр – ϑ 1 ) – 1] y ( ϑ 1 )
после несложных преобразований приведем к неравенству y′(ϑ1) sin(ϑпр – ϑ1) ,
y ( ϑ 1 ) ≤ cos( ϑ пр – ϑ 1 ) – 1
которое равносильно неравенству ctg
τ 0 ( ϑ 1 )
ϑ – ϑ пр 1
≥ ctg 2
.
Последнее неравенство равносильно неравенству (22), и существование трехимпульсных перелетов для 9 1 < 9 0 доказано.
Разность фаз на фиксированное начало 9 = 9 1 трехимпульсного маневра встречи, определяющая единственное положение соединения для заданных параметров относительного движения с ( 9 1 ) ( i = 2, 3, 4) и параметров встречи 9 пр , V , рассчитывается с помощью первого уравнения системы (15) по формуле:
с 1 *( 9 1 ) = 3 тА u 2 + 3 А9 V пр . (24)
Параметры относительного движения на начало маневрирования 9 = 9 N определяются по формулам прогнозирования невозмущенного движения (2). Положение соединения достигается при пассивном полете с момента 9 = 9
N до момента 9 = 9 1 .
На рис. 3 представлены границы оптимального диапазона разности фаз с 1 min( 9 ) = с 1 *( 9 ) и с m ax( 9 ) для трех-импульсных перелетов с варьируемым моментом 9 ( 9 N < 9 < 9 0 ) приложения первого импульса. Расчет произведен для параметра e * = 0,6.
На момент 9 N = 0 значение параметра, определяющего ориентацию разности векторов эксцентриситетов, ф = 90 ° . Фиксированный параметр 9 п р = 360 ° . Значения 9 N и 9 пр представлены на рисунке вертикальными прямыми.
Продолжительность оптимального двухимпульсного перелета между исходными орбитами т 0 = 208,07 ° , отсюда следует, что 9 0* = 151,93 ° . Нижняя граница диапазона определена по формуле (24).
Верхняя граница диапазона разности фаз определена по формуле:
с 1 тах( 9 ) = с 1 *( 9 0) + 3 с 2 ( 9 0 - 9 ). (25)
Пунктирные линии на рисунке соответствуют пассивному полету корабля с момента 9N до момента 91 начала трехимпульсного перелета со встречей

Рис. 3. Диапазон оптимальных фаз для трехимпульсных перелетов первого типа
На рис. 4 представлены траектории относительного движения с момента 9 N до момента 9 пр . На момент 9 N все траектории начинаются на линии с разностью высот y = 0,8, траектории отличаются значением момента первой коррекции 9 1 , координаты первой коррекции трехимпульсного перелета отмечены на рисунке красным кружком.

Рис. 4. Траектории относительного движения с мо мента 9N до момента 9пр для перелетов первого типа
Рассмотрим перелеты второго типа. Будем считать, что оптимальный двух-импульсный перелет начинается в момент 3 = 3 N а его продолжительность — т ° ( 3 N ). Трехимпульсные перелеты с закрепленным началом 3 1 = 3 N существуют, если назначаемый момент окончания перелета 3 п варьируется на интервале 3 N + т ° ( 3 N ) < Р3 пр < 3 N + 2 п . В самом деле, продолжительность двухимпульсного перелета, определяемого системой уравнений (16)–(17), не превышает продолжительности оптимального двухимпульс-ного перелета между исходными орбитами АК и ОС, если оба перелета начинаются в один и тот же момент времени. Очевидно, что минимальное значение 3 пр для трехимпульсного перелета второго типа
3min = 3„ + т°(3,) + е, пр N N , где е > ° — сколь угодно малая величина.
Разность фаз на фиксированное начало 3 1 = 3 N трехимпульсного маневра встречи, определяющая единственное положение соединения, рассчитывается по формуле (24) и представлена на рис. 5. Исходные параметры относительного движения на момент 3 1 = 3 N = ° суть e * = 0,6; ф (°) = 90 ° . Назначаемый момент встречи 3 п варьируется на интервале т ° (°) < 3 пр < 36 п ° ° , где т ° (°) = 151,93 ° .

Рис. 5. Оптимальная разность фаз на начало трехим-пульсных перелетов второго типа
На рис. 6 представлены траектории относительного движения с момента 3N= ° до назначенного (варьируемого) момента 3пр. На момент 3N все траектории начинаются на линии с разностью высот y = 0,8. Красным цветом выделена траектория двухимпульсного перелета со встречей.

Рис. 6. Траектории относительного движения с момента 3 N = 0 до варьируемого момента 3 пр для перелетов второго типа
Построение трехимпульсных перелетов для заданных параметров относительного движения
Перелеты первого типа ( с фиксацией 3 пр ). Первый этап построения перелета — решение трансцендентного уравнения (10) относительно 3 1 °, определяющего верхнюю границу значений 3 1 для начала трехимпульсных перелетов, 3 N < 3 1 < 3 ° . Верхняя граница диапазона оптимальных фаз на момент 3 N определяется по формулам (24) и (25). Нижняя граница — c mi n( 3 N ) = c *( 3 N ).
Трехимпульсный перелет для заданных значений ci(3N) (i = 1, 2, 3, 4) существует при условии, что c min(3 n) < c 1(3 n ) < c max(3 n ). (26)
Момент 31 начала перелета находится как решение трансцендентного уравнения на интервале 3N < 31 < 3° c 1(3) = c min(3), (27)
где нижняя граница оптимальных фаз c1min(3) определяется по формуле (24) и достигается после пассивного полета с момента 3N до момента 3. Разность фаз на момент 3 определяется по формуле c 1(3) = c 1(3n) - 3c2(3 - 3n). (28)
На рис. 3 моменты 3 1 определяются как точки пересечения пунктирных прямых (28) с функцией c 1 *( 3 ) = c “n ( 3 ).
При подготовке полета АК достаточно хорошо известны характеристики ракетоносителя, определяющие размер и форму орбиты АК, т. е. параметры c 2 , c 3 , c 4 относительного движения АК и ОС на момент 3 N разрешенного начала маневрирования. Разность фаз назначается, исходя из неравенства (26).
Перелеты второго типа (с фиксацией 31 = 3N). Продолжительность т0(3N) оптимального двухимпульсного перелета определяется по формуле (5). Разность фаз на фиксированное начало 31 = 3N трехимпульсного маневра встречи рассчитывается по формуле (24). Нижняя граница диапазона разностей фаз соответствует 3п = т0(3 N); верхняя граница 3 = ^max, р е ^max ^ 3 , 2^ — априори пр пр пр N заданная продолжительность трехим-пульсного перелета. Значение параметра c 1(3 N) должно удовлетворять неравенству вида (26). Значение 3пр находится как решение уравнения c 1(3N) = c*(3N, 3пр), (29)
где c *( 3 n 3 пр ) = 3 тА и 2 + Жр - 3 N ) V _,, в соответствии с формулой (24).
Список литературы Перелеты со встречей продолжительностью не более витка между близкими околокруговыми компланарными орбитами
- Микрин Е.А., Орловский И.В., Брагазин А.Ф., Усков А.В. Новые возможности автономной системы управления модернизированных кораблей "Союз" и "Прогресс" для реализации "быстрой" встречи с МКС // Космическая техника и технологии. 2015. № 4(11). С. 58-67.
- Муртазин Р.Ф. Схемы ускоренного доступа к орбитальной станции для современных космических кораблей // Космические исследования. 2014. Т. 52. № 2. С. 162-175.
- Брагазин А.Ф. Управление сближением космических аппаратов (навигация, наведение, коррекция движения). Королёв: РКК "Энергия", 2018. 472 с.
- Брагазин А.Ф., Усков А.В. Перелеты со встречей продолжительностью от одного до двух витков между околокруговыми компланарными орбитами // Космическая техника и технологии. 2020. № 1(28). С. 85-97. 10.33950/ spacetech-2308-7625-2020-1-85-97. DOI: 10.33950/spacetech-2308-7625-2020-1-85-97
- Баранов А.А. Маневрирование космических аппаратов в окрестности круговой орбиты. М.: Спутник, 2016. 512 с.
- Clohessy W.H., Wiltshire R.S. Terminal guidance system for satellite rendevous // Journal of Aerospace Sciences. September 1960. P. 653-658.
- Ермилов Ю.А., Иванова Е.Е., Пантюшин С.В. Управление сближением космических аппаратов. М.: Наука, 1977. 448 с.
- Бажинов И.К., Гаврилов В.П., Ястребов В.Д. и др. Навигационное обеспечение полета орбитального комплекса "Салют-6" - "Союз" - "Прогресс". М.: Наука, 1985. 376 с.