СПОСОБ УДЕРЖАНИЯ КОСМИЧЕСКОГО АППАРАТА НА ГЕОСИНХРОННОЙ 24-ЧАСОВОЙ ОРБИТЕ
METHOD OF HOLDING SPACECRAFT IN GEOSYNCHRONOUS 24-HOUR ORBIT

Предлагаемое изобретение относится к области космической техники и может быть использовано для удержания космического аппарата (КА) по долготе восходящего узла геосинхронной 24-часовой орбиты на условиях удержания КА на геостационарной орбите, когда регламент удержания по широте и гринвичской долготе не превышает ±0,05° относительно рабочей точки стояния. При этом считается, что поскольку допускается активное существование КА на геосинхронных (наклонных) 24-часовых орбитах, то для таких КА колебания по долготе вне области ±0,05° по широте, которые укладываются в регламент, и уходы по широте за пределы ±0,05° не противоречат международным требованиям.
В книге [1] Г.М.Чернявский, В.А.Бартенев, В.А.Малышев «Управление орбитой стационарного спутника», М., Машиностроение, 1984 г., на стр.42, 43, 134-136 приведен метод использования гравитационных полей Луны и Солнца для целенаправленного изменения наклонения плоскости орбиты спутника (КА), взятый за прототип. Суть метода сводится к тому, что осуществляя выведение КА на геостационарную околокруговую орбиту, определенным образом расположенную относительно Луны и Солнца, можно обеспечить практически вдвое меньший интервал изменения наклонения орбиты КА, чем при выведении его на строго экваториальную орбиту. Это же относится и к выведению на орбиты с наклонениями, измеряемыми градусами, которые правильнее считать геосинхронными 24-часовыми орбитами.
Согласно данному методу:
0
0
0
a.с.
1. Определяют начальные параметры наклонения i и долготы восходящего узла орбиты выведения Ω с учетом эпохи запуска КА D на орбиту и срока его активного существования (САС) t
0
0
0
0
a.c.
0
0
a.c.
0
<mrow><msub><mi>i</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>Δ</mi><mi>i</mi></mrow>
Выбранный с учетом конкретных расположений орбит Луны и Солнца на эпоху D начальный параметр Ω обеспечивает для
(Δi - максимальное изменение наклона плоскости орбиты КА в конце САС к плоскости орбиты КА в начале САС) сначала уменьшение наклонения вплоть до нуля, а затем его монотонное возрастание вновь до значения i. Зависимости Δi=ƒ(D, t) и Ω=ƒ(D, t) известны расчетным путем, пример таких зависимостей приведен в [1] на стр.136. Читая [1], следует понимать, что изменение за САС наклонения при i=0 равносильно максимальному монотонному изменению наклона Δi плоскости орбиты КА в конце САС к плоскости орбиты КА в начале САС.
0
0
2. Определяют новое значение i, уменьшая прежнее значение i на величину максимальной ошибки выведения по наклонению δi.
3. Прогнозируют расчетным путем значение наклонения орбиты КА на конец САС при максимальных ошибках выведения по наклонению и долготе восходящего узла орбиты ±δi, ±δΩ.
0
0
4. Из полученных вариантов расчета, с учетом варианта [δi=0; δΩ=0], определяют такие начальные значения i и Ω, при которых максимальное реализованное начальное наклонение будет равно максимальному конечному наклонению в конце САС.
0
Метод предполагает отсутствие активных коррекций наклонения. Оптимальное значение Ω порядка 270°, конечное значение всегда порядка 90°.
ε
опт
ε
опт
Метод использования гравитационных полей Луны и Солнца для целенаправленного изменения наклонения плоскости орбиты КА следует считать оптимальным вспомогательным методом, который вкупе с тактикой практического удержания геостационарного КА по долготе и широте на основе данных, полученных по траекторным измерениям, являет стратегию удержания КА в заданной области стояния (области рабочей точки стояния). При расчетах коррекций удержания геостационарных КА вектор изменения наклонения Δi направлен в центр малой области ε радиуса r, если вектор наклонения i находится вне этой области, и направлен в противоположную сторону вековым изменениям наклонения, если i находится внутри этой области. В координатах [i;Ω] i=i; Ω=270°, где i - оптимальное значение «точки прицеливания», например равное половине области стояния по наклонению. В случаях перерыва в коррекциях наклонения такая стратегия гарантирует пребывание КА в заданной узкой области по широте [±(0,05 - 0,1)°] относительно точки стояния в течение 30 суток и более. Далее по тексту, для краткости, выражение «относительно точки стояния (орбитальной позиции)» опущено. Метод не имеет недостатков.
В данной работе рассматривается способ, использующий цикличность изменения эксцентриситета в течение САС и использующий в полном объеме вышеприведенный метод пассивного управления наклонением орбиты КА.
В настоящее время, кроме собственно удержания геосинхронного КА в заданной области долгот, требуются работы по коллокации (безопасному сосуществованию) космических аппаратов, находящихся в одной области стояния. Коллокацию следует рассматривать как способ удержания. Известен способ управления кластером находящихся на геостационарной орбите спутников (RU 2284950 С2). Согласно данному способу, включающему измерение параметров орбиты каждого КА, определение по ним текущих значений орбитальных элементов каждого КА, сравнение их с требуемыми и проведение коррекций периода обращения, наклонения и эксцентриситета орбиты, маневры на каждом из КА проводят с помощью двигателей малой тяги, переводя векторы наклонения КА в разнесенные относительно друг друга [кольцевые] области их допустимого изменения так, чтобы угол между линией, соединяющей текущее положение конца каждого вектора с центром соответствующей ему кольцевой области, и направлением на Солнце был равен увеличенной на 180° величине прямого восхождения Солнца, одновременно проводят коррекции векторов эксцентриситета с целью перевода этих векторов в разнесенные относительно друг друга кольцевые области их допустимого изменения так, чтобы линия, соединяющая текущее положение каждого вектора с центром соответствующей ему кольцевой области (далее варианты):
1- отставала от направления на Солнце на половину углового расстояния при движении вектора эксцентриситета по окружности естественного дрейфа в пределах кольцевой области, далее на протяжении всего полета производят изменение относительного расстояния между КА в требуемых пределах за счет компенсации квазивекового приращения вектора наклонения каждого КА в сочетании с коррекцией вектора эксцентриситета, при которой в момент прохождения вектором эксцентриситета середины интервала между точкой входа окружности естественного дрейфа в кольцевую область допустимого изменения вектора эксцентриситета и точкой выхода из этой области линия, соединяющая центр окружности естественного дрейфа и центр соответствующей ему кольцевой области допустимого изменения вектора эксцентриситета, совпадала с направлением на Солнце, приводя тем самым к постоянству относительные векторы наклонения и эксцентриситета между КА;
2- совпадала с направлением на Солнце, далее на протяжении всего полета производят изменение относительного расстояния между КА в требуемых пределах за счет компенсации квазивекового приращения вектора наклонения каждого КА без коррекции вектора эксцентриситета, приводя тем самым к постоянству относительные векторы наклонения и эксцентриситета между КА.
Здесь «окружность естественного дрейфа» - окружность радиуса устойчивого эксцентриситета (определение и вывод формулы (14') см. ниже).
x
y
x
y
x
y
x
y
Суть данного способа сводится к синхронизации движения КА в фазовых плоскостях [i; i] и [e; e], где i=i·cosΩ; i=i·sinΩ; e=e·cos(Ω+ω); e=e·sin(Ω+ω); ω - аргумент широты перигея, причем синхронизации движения КА в обеих плоскостях, считая, что взаимная ориентация относительных векторов эксцентриситета и наклонения ΔЕ и ΔI 2-3 КА сохраняется, если использовать годовую цикличность естественного дрейфа вектора эксцентриситета и полугодовую цикличность вектора наклонения при условии компенсации вековой составляющей естественного дрейфа вектора наклонения. Обе цикличности обязаны своим существованием исключительно Солнцу, потому начальная и текущая ориентации векторов наклонения и эксцентриситета относительно Солнца являются необходимым условием достижения технического результата.
Предлагается способ, который не требует синхронизации изменения векторов наклонения (что исключительно важно для геосинхронного КА без соответствующего энергетического обеспечения) и не требует согласованных действий по управлению движением КА.
Целью изобретения является:
а) в отсутствие коллокации - гарантированное удержание геосинхронного КА в области ±(0,1-0,05)° по долготе восходящего узла орбиты без ограничений на наклонение (колебания по широте) со стороны баллистики;
б) гарантированная коллокация КА с другими КА в узкой области;
в) многократное сокращение расхода топлива на удержание по сравнению с расходами топлива на геостационарных КА, без потерь качества решения целевых задач на орбите.
Последнее возможно, когда на борту КА есть антенно-приводная система. А она всегда есть даже на геостационарных КА для узконаправленных антенн. Такая же в принципе система есть и на средствах наземных пунктов управления и связи. Требуется только определенным образом запрограммировать ее работу на витке движения КА. Что касается глобальных антенн, то для них антенно-приводная система не нужна в принципе.
Топлива на борту никогда много не бывает, стартовая масса КА на высотах ~36000 км имеет тенденцию к увеличению, поэтому заявленный способ удержания КА является актуальным и даже альтернативным в выборе существования КА - или на геосинхронной 24-часовой орбите, или геостационарной орбите.
уст
min
min
уст
min
0
min
0
пред
пред
пред
min
пред
тек
пред
уст
Поставленная цель достигается тем, что в способе удержания КА на геосинхронной 24-часовой орбите, включающем определение начальных параметров наклонения и долготы восходящего узла орбиты выведения с учетом эпохи запуска КА на орбиту и срока его активного существования, определение нового значения начального наклонения уменьшением прежнего значения на величину максимальной ошибки выведения по наклонению, прогнозирование расчетным путем значения наклонения орбиты КА на конец срока активного существования при максимальных ошибках выведения по наклонению и долготе восходящего узла орбиты, определение начальных значений наклонения и долготы восходящего узла, при которых максимальное реализованное начальное наклонение будет равно максимальному конечному наклонению в конце срока активного существования, введены новые операции, заключающиеся в том, что определяют величину устойчивого эксцентриситета е для выводимого КА; определяют предварительную величину минимального эксцентриситета е как сумму максимально допустимого эксцентриситета при удержании КА по долготе при нулевом наклонении (на геостационарной орбите) и отношения радиуса зоны отчуждения между КА, занимающими одну и ту же рабочую позицию на орбите, к радиусу геостационарной орбиты; определяют годограф изменения вектора эксцентриситета под действием светового давления в полярной системе координат [е;θ], где θ - угол между направлениями на Солнце в начальную эпоху (дату окончания приведения в область стояния) и на перигей орбиты КА, полярная ось направлена из центра Земли на Солнце и включает радиус-вектор е, при котором θ равно нулю, полярное расстояние равно эксцентриситету, вектор полярного расстояния в течение года описывает окружность радиуса е, и, зная амплитуду колебаний угла θ, уточняют величину е и годограф в целом; после выведения КА на геосинхронную 24-часовую орбиту и определения фактических параметров движения определяют срок «0», в течение которого фактическое начальное наклонение уменьшится до нуля; по годографу определяют оптимальное значение вектора начального эксцентриситета е на начальную эпоху, при котором линия апсид совмещена с радиус-вектором Солнца, и циклические годовые эволюции вектора е [е,θ], в отсутствие негативного влияния на них активных сил, гарантируют в конце срока «0» минимальное расхождение между векторами е и е; на этапе приведения КА в заданную область стояния одновременными коррекциями периода обращения и эксцентриситета доводят фактическое текущее значение вектора эксцентриситета до расчетного е; на этапе удержания, ближе к окончанию срока «0» проводят уточнение времени перехода в режим функционирования на геостационарной орбите, когда наклонение орбиты КА достигнет предельно допустимого наклонения i, при котором выход по широте за предельное значение для геостационарных КА осуществляется на границе номинальной области стояния по долготе; проводят коррекцию вектора эксцентриситета с целью, чтобы он равнялся номинальному для коллокации КА в заданной области стояния, а линия апсид орбиты геосинхронного КА совпала с линией узлов орбиты; в течение времени, пока наклонение производит эволюции в диапазоне [0-i], проводят активную коллокацию КА без взаимодействия с центрами управления смежными КА; при наклонении большем i увеличивают эксцентриситет до е с установкой вектора Лапласа в направлении от Солнца и в оставшееся время до окончания срока активного существования коррекции вектора эксцентриситета не проводят; в отсутствие коллокации на этапе приведения проводят одновременные коррекции периода обращения и эксцентриситета с тем, чтобы текущий эксцентриситет стал равен устойчивому; при уменьшении наклонения до значения i проводят коррекции эксцентриситета с тем, чтобы е был не более допустимого эксцентриситета при удержании КА по долготе (≡2угл.мин по долготе для областей стояния ±0,05°), при увеличении наклонения до значения i восстанавливают эксцентриситет до е с установкой вектора Лапласа в направлении на Солнце и в оставшееся время до окончания срока активного существования коррекции вектора эксцентриситета не проводят.
Суть способа заключаются в следующем.
1. Коллокация с геостационарными КА
z
z
ст
min
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
2
доп
<mrow><mi>Δ</mi><msub><mi>e</mi><mrow><msub><mi>z</mi><mn>0</mn></msub></mrow></msub></mrow>
<mrow><mo stretchy="false">(</mo><mi>Δ</mi><msub><mi>e</mi><mrow><msub><mi>z</mi><mn>0</mn></msub></mrow></msub><mo stretchy="false">)</mo></mrow>
Если учесть, что на интервале между траекторными измерениями 7 суток погрешность определения и прогнозирования положения вдоль орбиты δl составляет 4,5 км ([2] ОАО «ЭКА». Научно-технический отчет. Отработка технологии и оценка характеристик навигационно-баллистического обеспечения полетом КА 17Ф15М на этапе летных испытаний, М.,2010 г., стр.82), умножить ее на два (имеем ошибку определения положения КА пары) и прибавить 1 км (гарантия исключения всяких случайностей), получим гарантированную зону отчуждения обоих КА S = 10 км по всем опорным направлениям. Величина
отношения S к номинальному радиусу (r) геостационарной орбиты 42164 км есть допустимая минимальная разность эксцентриситетов при совпадающих направлениях на перигеи (векторов Лапласа) и минимальная сумма при диаметрально противоположных направлениях на перигеи - 0,00024 - 0,0003, которая не позволяет геосинхронному и геостационарному КА сближаться на расстояние меньшее S = 10 км в областях стояния по долготе ±0,05° и более при любых значениях i, Ω, при условии, что аргумент широты перигея (ω) орбиты геосинхронного КА не выходит за пределы ±(20-30)° относительно нуля или 180°, т.е. линия апсид орбиты геосинхронного КА находится в совпадающем положении относительно линии узлов, что подтверждено расчетами межспутниковых расстояний, отображенными на фиг.1-3 (на фиг.1: i=5°, i=0; Ω=Ω=270°, ω=ω=0, е=0,0006, e=0,0003; на фиг.2: i=4'30″, i=4'30″, Ω=0, Ω=90°, ω=0, ω=135°, e=0,0004, e=0,0001; на фиг.3: i=30″, i=30″, Ω=Ω=0, ω=180°, ω=0, е=0,00015, e=0,00015; индекс «1» относится к геосинхронному КА). Если прибавить к этой величине
максимально возможный эксцентриситет е при удержании геостационарных КА, рассчитанный по формуле (все составляющие формулы в радианах):
<mrow><msub><mi>e</mi><mrow><mi>д</mi><mi>о</mi><mi>п</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mrow><mi>Δ</mi><msub><mi>λ</mi><mrow><mi>н</mi><mi>о</mi><mi>м</mi></mrow></msub><mo>−</mo><mi>δ</mi><msub><mi>λ</mi><mrow><mi>о</mi><mi>п</mi><mi>р</mi></mrow></msub><mo>−</mo><mi>δ</mi><msub><mi>λ</mi><mrow><mi>у</mi><mi>п</mi><mi>р</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>
ном
где Δλ - половина ширины номинальной (данной) области стояния по долготе;
опр
ст
δλ=δl/r - погрешность определения и прогнозирования движения КА вдоль
орбиты;
упр
min
доп
min
0
δλ - степень свободы при управлении удержанием по долготе, порядка ±(0,01-0,015), получим востребованный позднее минимально допустимый эксцентриситет е для КА на геосинхронной 24-часовой орбите. Для областей стояния ±0,05° е≤0,0003, и e≤0,0006. Важная величина, дающая возможность рассчитать начальный вектор эксцентриситета е. До тех пор, пока наклонение орбиты геосинхронного КА не упадет с начальных (5-7,5)° (в зависимости от САС КА) до нескольких минут (например - 4,5 угл.мин для областей стояния по долготе ±(0,1-0,05°)), можно не учитывать колебания по долготе за счет наклонения и эксцентриситета. Последнее возможно ввиду того, что выход из номинальной области стояния по долготе на широтах, превышающих установленную норму для геостационарных КА, не следует считать нарушением международного регламента, поскольку в определенное время на витке геосинхронный КА находится в неоговоренном регламентом пространстве и реально никому не мешает.
min
уст
На этапе приведения КА в заданную область стояния в самом щадящем режиме, когда только лишь гасится пассивный дрейф в эту область (~2°/сутки), вырабатывается импульс приращения скорости порядка 17 м/с. На создание эксцентриситета с нуля до (e+2e)≈0,00150 требуется, исходя из формулы [3], К.Эрике «Космический полет», т.II, часть 1, стр.388, приращение скорости
<mrow><mi>Δ</mi><msub><mi>V</mi><mi>T</mi></msub><mo>=</mo><mfrac><mrow><mi>Δ</mi><mi>e</mi><mo>⋅</mo><mi>V</mi></mrow><mrow><mn>2</mn><mo>+</mo><mn>2</mn><mi>Δ</mi><mi>e</mi></mrow></mfrac><mo>=</mo><mn>2,38</mn><mtext> </mtext><mi>м</mi><mo>/</mo><mi>с</mi><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>
max
min
уст
где Δе=е=(е+2е)=0,0015;
V - скорость движения КА по круговой орбите, 3074 м/с.
Коррекции эксцентриситета, по необходимости, совмещаются с коррекциями периода обращения, поэтому энергозатраты на реализацию начального эксцентриситета самостоятельными не являются и общие энергозатраты на приведение не увеличивают.
2. Удержание геосинхронного КА относительно заданной точки стояния на тех же условиях, на которых удерживаются геостационарные КА, без коррекций наклонения вообще и с минимальными затратами на коррекции эксцентриситета в течение срока активного существования. За счет наклонения, уже большем 5 угл.мин, можно иметь большой эксцентриситет орбиты, что дает возможность вообще не проводить коррекции эксцентриситета (см.п.10, фиг.8).
Энергозатраты на коррекции эксцентриситета не являются существенными в общем топливном бюджете и их всегда, при условии не срочности исполнения этих коррекций, можно совмещать по времени с коррекциями периода обращения КА.
Введем понятие устойчивого эксцентриситета.
<mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow>
уст
Влияние малого импульса скорости (приращения скорости за секунду)
на эксцентриситет e и аргумент широты перигея ω описывается соотношениями [3], стр.388 (см. фиг.6, цифрами на которой обозначены: 1 - Солнце; 2 - Земля; 3 - КА; 4- орбита КА на текущую эпоху; 5 - перигей орбиты КА; 6 - орбита КА в эпоху, когда центр Земли, перигей и Солнце находятся на одной линии; 7 - радиус большого круга е), учитывающими тангенциальное и нормальное возмущения в плоскости орбиты:
<mrow><mfrac><mrow><mi>d</mi><mi>e</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mn>2</mn><mo>⋅</mo><mrow><mo>(</mo><mrow><mi>e</mi><mo>+</mo><mi>cos</mi><mi>η</mi></mrow><mo>)</mo></mrow><mo>⋅</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mi>sin</mi><mi>α</mi><mo>−</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mi>sin</mi><mi>η</mi><mo>⋅</mo><mi>cos</mi><mi>α</mi><mo>;</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>
<mrow><mfrac><mrow><mi>d</mi><mi>ω</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>2</mn><mi>sin</mi><mi>η</mi></mrow><mi>e</mi></mfrac><mo>⋅</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mi>sin</mi><mi>α</mi><mo>+</mo><mfrac><mrow><mn>2</mn><mi>e</mi><mo>+</mo><mi>cos</mi><mi>η</mi></mrow><mi>e</mi></mfrac><mo>⋅</mo><mi>cos</mi><mi>α</mi><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>
где ϑ - средняя скорость движения КА, 3074 м/с;
η = α - θ - истинная аномалия;
θ - угол между направлением на Солнце и на перигей орбиты КА.
Тогда, подставляя выражение η в (3) и используя формулы разности двух углов, получаем:
<mtable columnalign="left"><mtr><mtd><mfrac><mrow><mi>d</mi><mi>e</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mrow><mo>(</mo><mrow><mn>2</mn><mi>e</mi><mo>⋅</mo><mi>sin</mi><mi>α</mi><mo>+</mo><mn>2</mn><mi>cos</mi><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><mo>⋅</mo><mi>sin</mi><mi>α</mi></mrow></mrow><mo>+</mo></mtd></mtr><mtr><mtd><mrow><mrow><mo>+</mo><mn>2</mn><mi>sin</mi><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo>⋅</mo><mi>sin</mi><mi>α</mi><mo>−</mo><mi>sin</mi><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><mo>⋅</mo><mi>cos</mi><mi>α</mi><mo>+</mo><mi>cos</mi><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo>⋅</mo><mi>cos</mi><mi>α</mi></mrow><mo>)</mo></mrow><mo>=</mo></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mo stretchy="false">(</mo><mn>2</mn><mi>e</mi><mo>⋅</mo><mi>sin</mi><mi>α</mi><mo>+</mo><mi>sin</mi><mn>2</mn><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><mo>−</mo><mn>0,5</mn><mo>⋅</mo><mi>sin</mi><mn>2</mn><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><munder><munder><mrow><mo>+</mo><mn>2</mn><msup><mrow><mi>sin</mi></mrow><mn>2</mn></msup><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo>+</mo><msup><mrow><mi>cos</mi></mrow><mn>2</mn></msup><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo></mrow><mo stretchy="true">︸</mo></munder><mrow><mi>sin</mi><mi>θ</mi><mo>⋅</mo><mrow><mo>(</mo><mrow><msup><mrow><mi>sin</mi></mrow><mn>2</mn></msup><mi>α</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></munder></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mrow><mo>(</mo><mrow><mn>2</mn><mi>e</mi><mo>⋅</mo><mi>sin</mi><mi>α</mi><mo>+</mo><mn>0,5</mn><mo>⋅</mo><mi>sin</mi><mn>2</mn><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><mo>+</mo><msup><mrow><mi>sin</mi></mrow><mn>2</mn></msup><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo>+</mo><mi>sin</mi><mi>θ</mi></mrow><mo>)</mo></mrow><mtext> </mtext><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mtd></mtr></mtable>
Первое слагаемое, как минимум, на два порядка меньше остальных и не является постоянным членом, тогда
<mrow><mfrac><mrow><mi>d</mi><mi>e</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>⋅</mo><mi>sin</mi><mn>2</mn><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><mo>+</mo><msup><mrow><mi>sin</mi></mrow><mn>2</mn></msup><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo>+</mo><mi>sin</mi><mi>θ</mi></mrow><mo>)</mo></mrow><mtext> </mtext><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>
При θ = 0
<mrow><mfrac><mrow><mi>d</mi><mi>e</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>⋅</mo><mi>sin</mi><mn>2</mn><mi>η</mi><mo>⋅</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mtext> </mtext><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>
Аналогично рассуждая, будем иметь для скорости движения (аргумента широты) перигея:
<mrow><mfrac><mrow><mi>d</mi><mi>ω</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mi>e</mi></mfrac><mo>⋅</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mi>ϑ</mi></mfrac><mo>⋅</mo><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>⋅</mo><mi>sin</mi><mn>2</mn><mi>α</mi><mo>⋅</mo><mi>sin</mi><mi>θ</mi><mo>+</mo><msup><mrow><mi>sin</mi></mrow><mn>2</mn></msup><mi>α</mi><mo>⋅</mo><mi>cos</mi><mi>θ</mi><mo>+</mo><mi>cos</mi><mi>θ</mi></mrow><mo>)</mo></mrow><mtext> </mtext><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>
При θ = 0
<mrow><mfrac><mrow><mi>d</mi><mi>ω</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow><mrow><mi>e</mi><mo>⋅</mo><mi>ϑ</mi></mrow></mfrac><mo>⋅</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msup><mrow><mi>sin</mi></mrow><mn>2</mn></msup><mi>η</mi></mrow><mo>)</mo></mrow><mtext> </mtext><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>
Далее,
<mrow><mi>d</mi><mi>t</mi><mo>=</mo><mfrac><mrow><mi>d</mi><mi>α</mi></mrow><mrow><mi>n</mi><mo>−</mo><mi>n</mi></mrow></mfrac><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>10</mn></mrow><mo>)</mo></mrow></mrow>
<mrow><mi>n</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>T</mi></mfrac><mo>=</mo><mn>0,000073</mn></mrow>
-1
где
- среднее движение КА, с;
тогда
Интегрируем на сутках (на витке):
<mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow>
Дадим оценку
. Световое давление описывается формулой
<mrow><mi>P</mi><mo>=</mo><mfrac><mrow><mi>S</mi><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mi>A</mi></mrow><mo>)</mo></mrow></mrow><mi>c</mi></mfrac><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>15</mn></mrow><mo>)</mo></mrow></mrow>
2
2
где S - мощность световой волны, падающей на 1 м поверхности тела, вт/м;
А - коэффициент отражения (А=0 для абсолютно черного тела);
с - скорость света, км/с.
3
2
S=1,4·10вт/м.
-6
2
Значение А зависит от отражающей способности деталей конструкции КА и в контексте данного технического решения должно включать (условно) гравитационное воздействие Солнца как «±» относительно положения, когда вектор Лапласа направлен на Солнце. Для реальных КА на высоте стационарной орбиты значения А находятся в пределах [0,28-0,44]. Исходя из А = 0,44, будем иметь Р, равное 6,72·10н/м. Поскольку сила светового давления F=S'·P, где S' - площадь миделевого сечения, то
<mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub><mo>=</mo><mi>a</mi><mo>=</mo><mfrac><mrow><mi>S</mi><mo>'</mo></mrow><mi>m</mi></mfrac><mo>⋅</mo><mi>P</mi><mtext> </mtext><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>16</mn></mrow><mo>)</mo></mrow></mrow>
<mrow><mi>k</mi><mo>=</mo><mfrac><mrow><mi>S</mi><mo>'</mo></mrow><mi>m</mi></mfrac></mrow>
<mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub><mn>0,174</mn><mo>⋅</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>−</mo><mn>6</mn></mrow></msup><mi>м</mi><mo>/</mo><msup><mi>с</mi><mn>2</mn></msup><mo>.</mo></mrow>
<mrow><msub><mover accent="true"><mi>ϑ</mi><mo>˙</mo></mover><mn>0</mn></msub></mrow>
<mrow><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mrow><mn>390</mn></mrow></mfrac></mrow>
-2
уст
Отношение
для современных отечественных геостационарных КА более или менее постоянно и равно порядка (2,3-2,6)·10. Тогда, к примеру, при k=0,0259
Как показывает численное интегрирование, период цикличности для эксцентриситета составляет несколько больше года - порядка 390 суток. Это происходит из-за того, что возмущения движения перигея при устойчивом эксцентриситете от гравитационного поля Солнца имеют не годовой, а полугодовой период с амплитудой колебания, как показывает раздельное интегрирование, порядка 0,00005. Подстановка
в уравнение (14') дает при
значение е = е порядка 0,00045, т.е. для того, чтобы скорость движения перигея равнялась скорости движения Солнца, необходимо иметь устойчивый эксцентриситет. Оригинальный вывод формул (13-14') приведен для раскрытия сущности понятия устойчивого эксцентриситета.
уст
Энергозатраты на поддержание е практически отсутствуют: каждый год необходимо проводить коррекцию эксцентриситета на величину его векового ухода Δe = 0,000052 за счет нецентральности гравитационного поля Земли. Это очень малая величина, соизмеримая с суточными уходами эксцентриситета за счет совокупного влияния всех пассивных возмущающих факторов, однако если вековой уход не компенсировать, текущий эксцентриситет к концу САС будет далек от устойчивого.
Реализация предлагаемого способа предполагает выполнение следующей последовательности операций:
1. Определяют начальные значения наклонения и долготы восходящего узла орбиты выведения.
Данная операция аналогична совокупной последовательности операций 1-4 прототипа.
2. Определяют величину устойчивого эксцентриситета для выводимого на орбиту геосинхронного КА.
уст
Операция заключается в расчете е на основе формул (14', 15, 16) с учетом реальных данных по отражающей способности деталей конструкции КА (в основном панелей солнечных батарей), средней площади миделевого сечения и массы КА. Влияние гравитационного воздействия на КА Солнца можно не учитывать на годовом интервале, но по окончании годового интервала скорректировать положение перигея в направлении на Солнце. Скорректировать придется не более (4-5)°, что потребует не более 0,055 м/с, которые не превышают 2% топливного годового бюджета удержания по долготе при функционировании КА по целевому назначению и 0,09% топливного годового бюджета при сравнении с геостационарными КА, удерживаемых по широте (наклонению).
3. Определяют первое приближение величины минимального эксцентриситета, реализуемого на годовом интервале.
min
Величину минимального эксцентриситета е определяют из соотношения:
<mrow><msub><mi>e</mi><mrow><mi>min</mi></mrow></msub><mo>=</mo><msub><mi>e</mi><mrow><mi>д</mi><mi>о</mi><mi>п</mi></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>S</mi><mi>Z</mi></msub></mrow><mrow><msub><mi>r</mi><mrow><mi>c</mi><mi>m</mi></mrow></msub></mrow></mfrac><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>17</mn></mrow><mo>)</mo></mrow></mrow>
доп
min
уст
как сумму максимально допустимого эксцентриситета е при удержании КА по долготе при наклонении орбиты близком к нулю, рассчитанного по формуле (1), и отношения радиуса зоны отчуждения между КА в состоянии коллокации к радиусу геостационарной орбиты. Поскольку e>e, центр полярной системы координат вынесен за большой круг годографа (см. п.5).
4. Проводят траекторные измерения.
Траекторные измерения проводят после выведения КА на геосинхронную орбиту. Эта операция - типовая на всем САС.
5. Определяют (строят) годограф изменения вектора эксцентриситета под действием светового давления.
Принципиальная схема годографа приведена на фиг.7. Введены следующие обозначения:
1 - вектор е;
min
2 - вектор минимального эксцентриситета е, направлен на Солнце;
mах
3 - вектор максимального эксцентриситета е, направлен на Солнце;
уст
4 - указатель направления движения против часовой стрелки вектора е;
уст
5 - большой круг, описываемый концами векторов е и е;
уст
6 - текущий вектор е.
Большой круг годографа разделен на равные участки, представляющие временные интервалы (сутки, недели). Зная текущую дату и значение е, можно всегда сказать, когда будем иметь те или иные реальные значения эксцентриситета и угла θ отклонения вектора Лапласа (направления на перигей) от направления на Солнце.
min
max
θ
z
1
2
1
2
1
2
1
2
1
2
1
2
1
2
z
1
2
1
2
min
max
уст
z
min
max
<mrow><mi>Δ</mi><msub><mi>e</mi><mrow><msub><mi>z</mi><mn>0</mn></msub></mrow></msub></mrow>
Соотношение e/e годографа гарантирует амплитуду A = 48° колебаний угла θ (фиг.7) (колебания в пределах ±24° относительно нуля - начального направления на Солнце), что соответствует тем же колебаниям линии апсид относительно линии узлов, восходящий узел которой фиксирован в пространстве, согласно прототипу, сначала в районе 270°, а с середины САС - в районе 90°, благодаря выбору начальных параметров выведения (операция 1). Эти колебания соизмеримы с допуском на уход линии апсид от линии узлов орбиты КА, находящегося в состоянии коллокации с другими КА с момента приведения в точку стояния. В связи с этим уточним величину Δe для расхождения по аргументу широты перигея на 24°. На фиг.4, фиг.5 приведены наихудшие варианты расчета межспутниковых расстояний: фиг.4 - i=30″; i=30″ (величина наклонения роли не играет); Ω=Ω=0; ω=336°; ω=24°; e=0,00068; е=0,00030; фиг.5 - Ω=Ω=0; ω=336°; ω=90°; e=0,00068; е=0,00030, при которых, однако, межспутниковое расстояние Sz больше 10 км. Окончательно будем иметь Δе = 0,00038 (что оптимально для функционирования при диаметрально разнесенных перигеях орбит: e=е=0,00019 и (Ω+ω)=(Ω+ω)±180°), e=0,00068, e≈0,00158 (в зависимости от конкретного значения е), или, если быть более точным (см. стр.7,
=0,00024), Δе=0,00032, e=0,00061, e≈0,00151. Угол θ (фиг.7) изменяется во все тех же пределах: ±24° относительно нуля.
6. Определяют срок «0», в течение которого реализованное начальное наклонение уменьшится до нуля.
пред
Это большой срок, исчисляемый годами. Хотя нас на самом деле интересует срок, в течение которого наклонение уменьшится до известного значения i, составляющего несколько угловых минут, нет необходимости его определять на начальном этапе функционирования КА, поскольку его все равно приходится уточнять.
0
7. Определяют оптимальное значение вектора начального эксцентриситета е на начальную эпоху.
min
max
Годографом на момент начала функционирования КА по целевому назначению определены два значения: е, е, расходящихся на 180° и удовлетворяющих условию совмещения линии апсид с радиус-вектором Солнца.
0
min
ном
Из этих двух векторов эксцентриситета выбирают один - е, эволюция которого в конце срока «0» гарантирует минимальное расхождение между текущим вектором е и е и потребует проведение коррекции текущего вектора эксцентриситета до номинального для коллокации в заданной области стояния значения (е) с наименьшими энергозатратами.
0
Пункты 5-7 вероятно потребуют ограничений на даты старта. Запуски КА возможно будут проводиться в «окна» зимнего и летнего солнцестояний. Однако выставление ограничений на даты старта по тем или иным соображениям - обычная практика. Итак, значение е снимают с годографа.
0
8. На этапе приведения КА в заданную область стояния одновременными коррекциями периода обращения и эксцентриситета доводят фактическое текущее значение эксцентриситета до расчетного е, тем самым совмещая линию апсид с линией узлов.
9. На этапе удержания проводят уточнение времени перехода КА в режим функционирования на геостационарной орбите.
Уточнение проводят ближе к окончанию срока «0».
пред
10. Определяют предельное значение наклонения i, при котором выход по широте за предельное значение для геостационарных КА осуществляется на границе номинальной области стояния по долготе.
пред
пред
Наклонение i необходимо знать и при достижении наклонения значения i переходить в режим удержания (возможно коллокации) на условиях активного существования геостационарных КА, поскольку, если не в заданной области стояния по широте и долготе, то в смежных областях находятся (могут находиться) геостационарные КА, и существует норма (регламент) ширины области стояния на геостационарной орбите. Для произвольных областей стояния относительно центральной точки
<mrow><msub><mi>i</mi><mrow><mi>п</mi><mi>р</mi><mi>е</mi><mi>д</mi></mrow></msub><mo>=</mo><mn>2</mn><msub><mi>e</mi><mrow><mi>min</mi></mrow></msub><mo>⋅</mo><mfrac><mrow><mi>Δ</mi><msub><mi>ϕ</mi><mrow><mi>н</mi><mi>о</mi><mi>м</mi></mrow></msub></mrow><mrow><mi>Δ</mi><msub><mi>λ</mi><mrow><mi>н</mi><mi>о</mi><mi>м</mi></mrow></msub></mrow></mfrac><mo>+</mo><mi>a</mi><mi>r</mi><mi>c</mi><mi>t</mi><mi>g</mi><mrow><mo>(</mo><mrow><mfrac><mrow><mn>1</mn><mo>−</mo><mi>cos</mi><mi>Δ</mi><msub><mi>ϕ</mi><mrow><mi>н</mi><mi>о</mi><mi>м</mi></mrow></msub></mrow><mrow><mn>2</mn><msqrt><mrow><mi>cos</mi><mi>Δ</mi><msub><mi>ϕ</mi><mrow><mi>н</mi><mi>о</mi><mi>м</mi></mrow></msub></mrow></msqrt></mrow></mfrac></mrow><mo>)</mo></mrow><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>18</mn></mrow><mo>)</mo></mrow></mrow>
ном
ном
пред
где ±Δφ и ±Δλ - соответственно номинальные области стояния КА по широте и долготе. Первое слагаемое прямо следует из геометрии общей области стояния, второе учитывает колебания по долготе за счет наличия наклонения орбиты. Для области ±(0,1-0,05)° i=4,5 угл.мин (см. фиг.8).
Коррекции гринвичской долготы восходящего узла (то же - коррекции периода обращения КА) планируют так, чтобы график изменения вектора эксцентриситета согласно годографу строго соблюдался.
11. При уменьшении текущего наклонения до предельного значения проводят коррекцию вектора эксцентриситета.
ном
Текущий вектор эксцентриситета переводят в номинальный вектор е.
12. Проводят активную коллокацию КА (в отличие от пассивной коллоакации, когда совмещение линии узлов и линии апсид в пределах ±24° происходило автоматически) без взаимодействия с центрами управления смежными КА.
n
n
n
n
n
n
n
n
n
n
Как правило, активную коллокацию КА проводят по согласованным схемам. Все схемы - аналоги сводятся к равноудалению точек прицеливания векторов е [e,(Ω+ω)] (n=1,2,…) и i [i,Ω] (n=1,2,…) в соответствующих фазовых плоскостях КА и поддержанию концов векторов е и i внутри соответствующих областей выбранных радиусов, центрами которых являются соответствующие точки прицеливания. Идеальным вариантом считается для двух КА разнесение долгот восходящих узлов (Ω) и прямых восхождений перигеев (Ω+ω) точек прицеливания на 180°, причем аргументы широты перигеев КА должны быть близки нулю. Для трех КА цифру 180 заменяют на 120.
Операция 12 является существенным и отличительным признаком изобретения. Суть ее сводится к следующему:
а) для цели активной коллокации нет необходимости заниматься коррекциями векторов наклонения: ни на геостационарных КА, ни на геосинхронном КА, достаточно поддерживать линию апсид орбиты только геосинхронного КА в заданных пределах относительно линии узлов его орбиты проведением коррекций вектора эксцентриситета, т.е.
<mrow><mi>ω</mi><mo>=</mo><mo>±</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>A</mi><mi>θ</mi></msub><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>19</mn></mrow><mo>)</mo></mrow></mrow>
либо
<mrow><mi>ω</mi><mo>=</mo><mi>π</mi><mo>±</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>A</mi><mi>θ</mi></msub><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>20</mn></mrow><mo>)</mo></mrow></mrow>
z
доп
ном
а, поскольку величина Δe≥е и больше е, то оптимальным вариантом активной коллокации будет разнесение перигеев орбит двух КА на 180° (фиг.3, стр.7), когда
<mrow><mi>ω</mi><mo>=</mo><mi>π</mi><mo>±</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>A</mi><mi>θ</mi></msub><mtext> </mtext><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mrow><mo>(</mo><mrow><mn>21</mn></mrow><mo>)</mo></mrow></mrow>
1
2
1
2
1
2
1
2
(на фиг.9: i=i=1', Ω=90°, Ω=0°, ω=45°, ω=315°, e=0,00019, е=0,00019),
притом что активная коллокация актуальна в течение не более двух месяцев;
б) нет необходимости во взаимодействии с центрами управления смежными КА, поскольку относительно простую задачу разнесения векторов эксцентриситета (или векторов Лапласа) можно решить мониторингом данных по тактике и стратегии удержания смежных КА (например, с помощью международной системы слежения NORAD) и адаптацией управления геосинхронным КА к этим данным, что может оказаться более эффективным, чем пытаться наладить тесный контакт с центрами управления смежными КА.
пред
min
13. Как только наклонение станет больше i, увеличивают эксцентриситет до значения е с установкой вектора Лапласа в направлении от Солнца и в оставшееся до окончания САС время коррекции вектора эксцентриситета не проводят.
14. В отсутствие коллокации на этапе приведения проводят одновременные коррекции периода обращения и эксцентриситета с тем, чтобы текущий эксцентриситет стал равен устойчивому.
пред
15. При уменьшении наклонения до значения i проводят коррекции эксцентриситета с тем, чтобы текущий эксцентриситет был не более допустимого при удержании КА по долготе.
пред
16. При увеличении наклонения до значения i восстанавливают эксцентриситет до устойчивого значения с установкой вектора Лапласа в направлении на Солнце и в оставшееся до окончания САС время коррекции вектора эксцентриситета не проводят.
0
к
Отметим, что согласно прототипу начальное и конечное значения долготы восходящего узла орбиты КА составляют соответственно 270° и 90°, более того, переход от Ω к Ω происходит скачкообразно, поэтому совмещение на всем САС линии узлов и линии апсид геосинхронного КА (действие по п.13 соответствует данной установке) гарантирует успешное сосуществование нескольких КА в заданной области стояния.
Отметим еще раз, что энергозатраты на коррекции эксцентриситета самостоятельными не являются и общие энергозатраты не увеличивают.
Предлагаемый способ удержания КА по долготе на геосинхронной 24-часовой орбите обеспечивает:
1 - удержание по долготе восходящего узла орбиты без энергозатрат на коррекции эксцентриситета;
2 - гарантированную коллокацию КА с другими КА в области рабочей точки стояния по долготе;
3 - сокращение на порядок энергозатрат на удержание против геостационарных КА за счет полного исключения коррекций наклонения;
4 - принципиальную альтернативу геостационарной орбите.