Способ приведения геосинхронного космического аппарата к рабочей долготе двигателями малой тяги

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

 

Предлагаемое изобретение относится к области космической техники. Оно может быть использовано для приведения космического аппарата (КА) с электрореактивной двигательной установкой (ЭРДУ) или другими двигателями малой тяги к рабочей геофизической долготе на геосинхронной или геостационарной орбите (ГСО) при больших начальных отклонениях КА от рабочей долготы и от номинального периода, а также для перевода КА на другую долготу.

Известен способ приведения КА к рабочей долготе на ГСО, описанный в статье [1] А.Е. Назаров «Использование комбинированного метода расчета программы коррекций приведения геостационарного КА для минимизации эксцентриситета орбиты», Вестник ФГУП НПО им. С.А. Лавочкина, 2009, №2, стр. 23-33. В нем приведение выполняется тремя импульсными коррекциями, рассчитываемыми с учетом ошибок измерений элементов орбиты КА и ошибок исполнения коррекций двигательной установкой. Этот способ принципиально не может быть применен, если используются двигатели малой тяги.

Известен способ, описанный в книге [2] Г.М. Чернявский, В.А. Бартенев, В.А. Малышев «Управление орбитой стационарного спутника», М., Машиностроение, 1984 г., стр. 120-126. В нем выполняют многошаговый процесс приведения двигателями малой тяги. На каждом шаге во время полета КА выполняют измерения орбитальных параметров и по результатам измерений вычисляют время начала и длительность воздействия коррекций на этом шаге. При вычислении прогнозируют результат путем решения дифференциальных уравнений движения КА с учетом корректирующего воздействия.

Недостатки этого способа в следующем. Способ не предусматривает ограничений на длительность работы двигателей при очередном включении. Точные измерения на каждом шаге приведения требуют около суток и более времени между коррекциями. Поскольку при применении двигателей со сверхмалой тягой порядка 0,1 Н и ограниченной длительностью работы могут потребоваться десятки шагов приведения, то его общая длительность будет недопустимо большой, особенно при больших начальных отклонениях периода и, соответственно, при большой начальной скорости движения КА по долготе.

В качестве прототипа принят способ, описанный в книге [3] «Космический мусор». В 2 книгах, Кн. 2. «Предупреждение образования космического мусора», под научной редакцией Г.Г. Райкунова, М.: ФИЗМАТЛИТ, 2014 г., стр. 50-52.

В этом способе после испытаний КА с пробными включениями и уточнением тяги двигателей последовательно выполняют: измерение и по их результатам расчет элементов орбиты КА; расчет коррекций, которым уточняют план приведения; коррекции периода орбиты для придания КА заданной скорости движения по долготе; измерение и расчет элементов орбиты для уточнения достигнутой скорости пассивного перемещения; уточнение плана торможения КА на заданной долготе; коррекции периода с целью торможения в окрестности заданной долготы; измерение и расчет элементов орбиты после коррекции периода; уточнение тяги и последнее включение двигателей для точной постановки спутника в заданную орбитальную позицию; выполнение последней коррекции.

Такой способ обычно используют для приведения при нормальной, штатной работе КА и его системы управления. При этом скорость пассивного перемещения, т.е. дрейфа КА по долготе от 0,4 градуса до 1,0 градуса в сутки. Это соответствует отклонению периода обращения КА от 95 до 240 секунд от номинального периода ГСО, составляющего 86164,09 секунды. При необходимости уменьшения длительности приведения скорость дрейфа и отклонение периода могут быть больше, платой за это является риск столкновений с отработавшими КА.

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

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

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

На фиг. 1 приведен пример схемы фазовой траектории приведения в изобретении, в координатах долготы и отклонения периода от номинального периода ГСО. На фиг. 1 и в описании введены обозначения:

L₀ - целевая рабочая долгота;

L_н - начальная долгота приведения;

L_к, ΔT_к - долгота и отклонение периода в начале коррекций приведения;

L₁₂, ΔT₁ - промежуточные расчетная долгота и положительное отклонение периода;

L₁₁, ΔT₁ - долгота начала завершающих коррекций и положительное отклонение периода;

L₂₂, ΔT₂ - промежуточные расчетная долгота и отрицательное отклонение периода;

L₂₁, ΔT₂ - долгота начала завершающих коррекций и отрицательное отклонение периода;

Начальной долготой приведения является L_н. От нее КА пассивно, без коррекций, перемещается к долготе L_к. За время этого перемещения выполняют измерения параметров орбиты КА. После этого рассчитывают план коррекций, начинающихся при достижении долготы L_к и отклонения периода ΔT_к.

Сначала рассчитывают план коррекций приведения непосредственно к целевой начальной долготе и к номинальному периоду геосинхронной орбиты. Затем для такого плана рассчитывают диапазон конечных отклонений долготы и периода после выполнения коррекций, возможные из-за ошибок измерений и исполнения коррекций. Если рассчитанные конечные отклонения не превышают допустимые значения, то выполняют коррекции и завершают приведение. Такая ситуация является редкой, особой и возможна в случае малых начальных отклонений.

Если рассчитанные конечные отклонения превышают заданные допустимые значения, то для приведения предварительно находят два варианта промежуточных величин. С этой целью выбирают две величины отклонений периода ΔT₁, ΔT₂ с разным знаком, ΔT₁ от 95 до 240 секунд и ΔT₂ от минус 240 до минус 95 секунд. Для каждого из этих выбранных отклонений ΔT₁, ΔT₂ вычисляют долготу соответственно L₁₁, L₂₁ такую, что для данного КА минимально время его приведения к заданной долготе стояния L₀ и к номинальному периоду из начальных положений L₁₁, ΔT₁ и L₂₁, ΔT₂. Для отклонений периода ΔT₁, ΔT₂ вычисляют соответственно долготы L₁₂, L₂₂ такие, что при отклонениях периода ΔT₁, ΔT₂ длительность некорректируемого движения КА от долготы L₁₂ к долготе L₁₁ и от долготы L₂₂ к долготе L₂₁ равна сумме длительности измерения орбиты и длительности смещения КА по долготе на величину возможной ошибки приведения к долготе соответственно L₁₂ с отклонением периода ΔT₁ и к долготе L₂₂ с отклонением периода ΔT₂.

Для выбора одного из двух вариантов рассчитывают планы коррекций приведения КА от долготы начала коррекций L_к к долготе L₁₂ с отклонением периода ΔT₁ и к долготе L₂₂ с отклонением периода ΔT₂. Рассчитывают полную длительность приведения к заданной целевой долготе L₀ и к номинальному периоду, а также требуемый для приведения расход характеристической скорости для двух вариантов приведения, через долготы L₁₂, L₁₁ с отклонением периода ΔT₁ или через долготы L₂₂, L₂₁ с отклонением периода ΔT₂. Исходя из рассчитанной длительности приведения и расхода скорости, выбирают вариант приведения через долготы L₁₂, L₁₁ с отклонением периода ΔT₁ или через долготы L₂₂, L₂₁ с отклонением периода ΔT₂.

На фиг. 1 показан вариант приведения через долготы L₁₂, L₁₁ с отклонением периода ΔT₁. Долгота L₁₂ с отклонениями периода ΔT₁ является той самой, после достижения которой при приведении в соответствии с сутью изобретения не существенно влияние ошибок измерения орбиты и исполнения оставшихся коррекций. Далее на схеме фиг. 1 от долготы L_к коррекциями по рассчитанному и выбранному варианту плана приводят КА с ошибкой на долготу L₁₂ с отклонением периода ΔT₁. Выполняют измерения элементов орбиты КА при его некорректируемом движении от долготы L₁₂ к L₁₁. По данным этих измерений вычисляют длительность достижения долготы L₁₁ и рассчитывают план коррекций приведения от долготы L₁₁ к заданной рабочей долготе L₀ и периоду Т₀. От момента достижения КА долготы L₁₁ выполняют коррекции приведения к заданной долготе стояния L₀ и периоду Т₀ и завершают приведение, с возможным выполнением контрольных измерений.

Технический результат изобретения достигается тем, что способ приведения геосинхронного космического аппарата к рабочей долготе двигателями малой тяги, включающий измерения орбиты КА, расчет плана коррекций приведения к заданной рабочей долготе и к номинальному периоду геосинхронной орбиты, выполнение коррекций орбитального движения КА по рассчитанному плану, отличается тем, что выполняют измерения параметров орбиты КА при его некорректируемом движении от начальной долготы приведения L_н к начальной долготе коррекций приведения L_к, рассчитывают план коррекций приведения КА от долготы L_к к заданной долготе L₀ и к номинальному периоду геосинхронной орбиты, рассчитывают отклонения долготы КА от долготы L₀ и отклонение периода от номинального периода, возможные после приведения по рассчитанному плану из-за ошибок измерений орбиты и ошибок исполнения коррекций, выбирают две величины отклонения периода ΔT₁, ΔT₂ с разным знаком, ΔT₁ от 95 до 240 секунд и ΔT₂ от минус 240 до минус 95 секунд, для каждого из двух выбранных отклонений периода ΔT₁, ΔT₂ вычисляют долготу соответственно L₁₁, L₂₁ такую, что для данного КА минимально время приведения КА к долготе L₀ и к номинальному периоду из начальных положений L₁₁, ΔT₁ и L₂₁, ΔT₂, для отклонений периода ΔT₁, ΔT₂ вычисляют соответственно долготы L₁₂, L₂₂, такие, что при отклонениях периода ΔT₁, ΔT₂ длительность некорректируемого движения КА от долготы L₁₂ к долготе L₁₁ и от долготы L₂₂ к долготе L₂₁ равна сумме длительности измерения орбиты и длительности смещения КА по долготе на величину возможной ошибки приведения к долготе соответственно L₁₂ с отклонением периода ΔT₁ и к долготе L₂₂ с отклонением периода ΔT₂, рассчитывают планы коррекций приведения КА от долготы начала коррекций L_к к долготе L₁₂ с отклонением периода ΔT₁ и к долготе L₂₂ с отклонением периода ΔT₂, рассчитывают полную длительность приведения к заданной долготе L₀ и к номинальному периоду, а также требуемый для приведения расход характеристической скорости для двух вариантов приведения, через долготы L₁₂, L₁₁ с отклонением периода ΔT₁ или через долготы L₂₂, L₂₁ с отклонением периода ΔT₂, исходя из рассчитанной длительности приведения и расхода скорости, выбирают вариант приведения через долготы L₁₂, L₁₁ с отклонением периода ΔT₁ или через долготы L₂₂, L₂₁ с отклонением периода ΔT₂, от долготы L_к коррекциями по рассчитанному и выбранному варианту плана приводят КА на долготу L₁₂ с отклонением периода ΔT₁ или на долготу L₂₂ с отклонением периода ΔT₂., при этом с ошибками в этих величинах, вызванными ошибками измерений и исполнения коррекций, выполняют измерения параметров орбиты КА при его некорректируемом движении от долготы L₁₂ к L₁₁ или от долготы L₂₂ к L₂₁, по данным измерений вычисляют длительность достижения долготы L₁₁ или L₂₁ и рассчитывают план коррекций приведения от долготы L₁₁ или L₂₁ к заданной рабочей долготе L₀ и к номинальному периоду, от момента достижения КА долготы L₁₁ или L₂₁ выполняют коррекции приведения к долготе L₀ и к номинальному периоду.

Предлагаемый способ реализуется с применением известных методов расчета [4-7].

Вначале методами [4] или [5, 6] рассчитывают план коррекций приведения за минимальное время от начальной долготы L_к к целевой долготе L₀ и к номинальному периоду. Затем, исходя из рассчитанного плана коррекций и известных для данной космической системы ошибок измерений и исполнения коррекций, методом [7] вычисляют возможные конечные отклонения долготы и периода от целевых значений. Если эти отклонения не превышают допустимые значения, то выполняют коррекции и завершают приведение. Такая ситуация является редкой, особой и возможна в случае малых начальных отклонений.

Затем выбирают две величины отклонения периода, ΔT₁ от 95 до 240 секунд и ΔT₂ от минус 240 до минус 95 секунд, в зависимости от характеристик КА и требований по точности его приведения.

Далее рассчитывают долготу L₁₁, такую, что для данного КА минимально время приведения КА от этой долготы L₁₁ к целевой рабочей долготе L₀ и к номинальному периоду. Для этого планируют коррекции максимального уменьшения периода от его отклонения ΔT₁ при долготе L₀ до номинальной величины периода, затем прогнозируют корректируемое движение КА на этом интервале. Долгота, получаемая по результату прогнозирования, принимается равной величине L₀-(L₀+L₁₁) и из нее получают величину L₁₁, с несущественной ошибкой. Далее получают предварительную величину L₁₂ как сумму долготы L₁₁ и разницы долгот, которую преодолевает КА при отклонении периода ΔT₁ и некорректируемом движении за время, достаточное для измерений в данной космической системе. Полученную предварительную величину L₁₂ уточняют несколькими итерациями следующим образом. Методом [6] планируют коррекции приведения к долготе L₁₂ и к отклонению периода ΔT₁ от начальной долготы коррекций L_к и отклонения периода ΔT_к. Методом [7] вычисляют утроенное среднеквадратическое отклонение долготы из-за ошибок измерений и исполнения коррекций. Долготу L₁₂ увеличивают на вычисленное отклонение и повторяют такие итерации планирования приведения к долготе L₁₂ и ее вычисление до тех пор, пока изменение отклонения в очередной итерации не станет достаточно малым, например, менее 1/10 его полной величины. Запоминают требуемый расход скорости и длительность приведения в последней итерации приведения к долготе L₁₂ и к отклонению периода ΔT₁.

Аналогично L₁₁, L₁₂ рассчитывают долготы L₂₁, L₂₂, планируют коррекции приведения к долготе L₂₂ и к отклонению периода ΔT₂, определяют требуемый расход скорости и длительность этого приведения.

Для каждого из двух вариантов приведения определяют расход скорости и длительность от начальной долготы коррекций L_к и отклонения периода ΔT_к до целевой рабочей долготы L₀ и номинального периода. Их получают добавлением к найденным величинам соответствующих величин для интервалов некорректируемого движения при отклонениях периода ΔT₁, ΔT₂, а также интервалов приведения от этих отклонений к номинальному периоду долготе L₀. На интервалах некорректируемого движения расход отсутствует, а длительность вычисляют по разницам крайних долгот этого движения и отклонениям периода. На интервалах приведения эти величины известны из прогнозирования при вычислении долгот L₁₁, L₂₁.

По результатам сравнения длительностей приведения и расходов скорости выбирают вариант для осуществления. В примере на фиг 1 приведение выполняют через промежуточную расчетную долготу L₁₂ и положительное отклонение периода ΔT₁.

По рассчитанному плану коррекций приводят КА к долготе L₁₂ и к отклонению периода ΔT₁ с ошибками. С учетом этих ошибок уточняют долготу L₁₁, вычисляют время прибытия КА на эту долготу и уточняют план коррекций окончательного приведения к долготе L₀ и к номинальному периоду.

При использовании описанного выше способа приведения конечные отклонения от долготы L₀ и от номинального периода определяются интервалом приведения от долготы L₁₁ и отклонения периода ΔT₁, с учетом ошибок предшествующего интервала приведения.

В случае наличия эксцентриситета и наклонения орбиты КА совершает суточные колебания по долготе. В изобретении долготой считается ее среднее значение в течение звездных суток. Эксцентриситет при необходимости корректируют одновременно с периодом известными способами [2, 4, 5].

Для оценки конечных отклонений используем метод [7]. При расчете в качестве примера примем обозначения и значения величин:

ΔT₁=200 с;

m=2500 кг - масса КА;

Δt=8 ч - длительность максимального воздействия на витке;

τ₀=0,36 с - среднеквадратическая ошибка измерения периода;

δ_n=0,02 - относительная среднеквадратическая ошибка трансверсальной коррекции;

Т=86164,09 с - номинальный период ГСО;

k=0,32 с/км - гравитационный параметр;

n - длительность приведения в целых сутках;

F - тяга ЭРДУ.

Длительность приведения в целых сутках:

Ошибка приведения на долготу L₀ и период Т, по периоду:

Ошибка приведения на долготу L₀ и период Т, по долготе:

Из формул (1-3) ошибки приведения:

ΔT=12 с, ΔL=0,05° при тяге F=0,2 Н;

ΔT=8,2 с, ΔL=0,08° при тяге F=0,1 Н.

Такие ошибки достаточно малы и позволяют перейти к удержанию КА на целевой рабочей долготе L₀.

Оценки ошибок приведения на долготу начала завершающих коррекций L₁₁ и отклонение периода ΔT₁:ΔT=15 с, ΔL=0,25° при тяге F=0,1 Н и примерном суммарном изменении периода на 600 с. Эти ошибки также достаточно малы для завершения приведения к долготе L₀ и периоду Т с необходимой точностью согласно изобретению.

Источники информации.

1. А.Е. Назаров Использование комбинированного метода расчета программы коррекций приведения геостационарного КА для минимизации эксцентриситета орбиты. Вестник ФГУП НПО им. С.А. Лавочкина, 2009, №2, стр. 23-33.

2. Г.М. Чернявский, В.А. Бартенев, В.А. Малышев «Управление орбитой стационарного спутника», М., Машиностроение, 1984 г., стр. 120-126.

3. «Космический мусор». В 2 книгах, Кн. 2. «Предупреждение образования космического мусора», под научной редакцией Г.Г. Райкунова, М.: ФИЗМАТЛИТ, 2014 г., стр. 50-52.

4. А.В. Соколов, Ю.П. Улыбышев Многовитковые маневры с малой тягой в окрестности геостационарной орбиты. Известия Академии наук. Теория и системы управления, 1999, №2, стр. 95-100.

5. А.Р. Агишев, В.Г. Касаткин Планирование коррекций для приведения КА на заданную позицию ГСО двигателями малой тяги. Труды МАИ, №102, 2018, Режим доступа: http://trudymai.ru/upload/iblock/782/Agishev-_Kasatkin_rus.pdf (Дата обращения: 15.11.2018).

6. А.Р. Агишев Расчет коррекций орбиты квазистационарного КА двигателями малой тяги для приведения к заданным долготе и периоду. Труды 61-й Всероссийской научной конференции МФТИ. Аэрокосмические технологии, 2018, стр. 98-99.

7. А.Р. Агишев Расчет среднеквадратических отклонений долготы и периода после завершения серии коррекций вблизи геостационарной орбиты. Труды МАИ, №100, 2018, Режим доступа: http://trudymai.ru/upload/iblock/04c/Agishev_rus.pdf (Дата обращения 04.10.2018).

Способ приведения геосинхронного космического аппарата (КА) к рабочей долготе двигателями малой тяги, включающий измерения орбиты КА, расчет плана коррекций приведения к заданной рабочей долготе и к номинальному периоду геосинхронной орбиты, выполнение коррекций орбитального движения КА по рассчитанному плану, отличающийся тем, что выполняют измерения параметров орбиты КА при его некорректируемом движении от начальной долготы приведения L_н к начальной долготе коррекций приведения L_к, рассчитывают план коррекций приведения КА от долготы L_к к заданной долготе L₀ и к номинальному периоду геосинхронной орбиты, рассчитывают отклонения долготы КА от долготы L₀ и отклонение периода от номинального периода, возможные после приведения по рассчитанному плану из-за ошибок измерений орбиты и ошибок исполнения коррекций, выбирают две величины отклонения периода ΔТ₁, ΔТ₂ с разным знаком: ΔТ₁ от 95 с до 240 с и ΔТ₂ от минус 240 с до минус 95 с, для каждого из двух выбранных отклонений периода ΔТ₁, ΔТ₂ вычисляют долготу соответственно L₁₁, L₂₁ такую, что для данного КА минимально время приведения КА к долготе L₀ и к номинальному периоду из начальных положений L₁₁, ΔТ₁ и L₂₁, ΔТ₂, для отклонений периода ΔТ₁, ΔТ₂ вычисляют соответственно долготы L₁₂, L₂₂ такие, что при отклонениях периода ΔТ₁, ΔТ₂ длительность некорректируемого движения КА от долготы L₁₂ к долготе L₁₁ и от долготы L₂₂ к долготе L₂₁ равна сумме длительности измерения орбиты и длительности смещения КА по долготе на величину возможной ошибки приведения к долготе соответственно L₁₂ с отклонением периода ΔТ₁ и к долготе L₂₂ с отклонением периода ΔТ₂, рассчитывают планы коррекций приведения КА от долготы начала коррекций L_к к долготе L₁₂ с отклонением периода ΔТ₁ и к долготе L₂₂ с отклонением периода ΔТ₂, рассчитывают полную длительность приведения к заданной долготе L₀ и к номинальному периоду, а также требуемый для приведения расход характеристической скорости для двух вариантов приведения: через долготы L₁₂, L₁₁ с отклонением периода ΔТ₁ или через долготы L₂₂, L₂₁ с отклонением периода ΔТ₂, исходя из рассчитанной длительности приведения и расхода скорости, выбирают вариант приведения через долготы L₁₂, L₁₁ с отклонением периода ΔТ₁ или через долготы L₂₂, L₂₁ с отклонением периода ΔТ₂, от долготы L_к коррекциями по рассчитанному и выбранному варианту плана приводят КА на долготу L₁₂ с отклонением периода ΔТ₁ или на долготу L₂₂ с отклонением периода ΔТ₁, при этом с ошибками в этих величинах, вызванными ошибками измерений и исполнения коррекций, выполняют измерения параметров орбиты КА при его некорректируемом движении от долготы L₁₂ к L₁₁ или от долготы L₂₂ к L₂₁, по данным измерений вычисляют длительность достижения долготы L₁₁ или L₂₁ и рассчитывают план коррекций приведения от долготы L₁₁ или L₂₁ к заданной рабочей долготе L₀ и к номинальному периоду, от момента достижения КА долготы L₁₁ или L₂₁ выполняют коррекции приведения к долготе L₀ и к номинальному периоду.