РАЗРАБОТКА МЕТОДИКИ ОПРЕДЕЛЕНИЯ И КОРРЕКЦИИ ПАРАМЕТРОВ РАБОЧЕЙ ОРБИТЫ КОСМИЧЕСКОГО АППАРАТА ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ ЗЕМЛИ

Полный текст

Введение. Одной из наиболее важных задач Федеральной космической программы России на период до 2025 г. является создание космических аппаратов дистанционного зондирования Земли. Неотъемлемой частью проектирования системы ДЗЗ считается определение параметров орбит, которые наиболее эффективны с позиций информативности, энергообеспечения и длительности активного существования КА. Разнообразные и многочисленные исследования в данной области проводятся в нашей стране и за рубежом; их результаты обобщены в ряде изданий [1-10]. Так в работе [6] рассмотрен подход к выбору технических параметров корректирующей двигательной установки КА ДЗЗ, исходя из характеристик съемочной аппаратуры, системы ориентации и требуемого срока активного существования КА, но не решается проблема видимости пункта приема информации (ППИ) и целевых задач использования КА ДЗЗ. В труде [7] представлен анализ параметров орбит КА ДЗЗ для обеспечения хорошего покрытия ближневосточных территорий, однако нет практических рекомендаций для анализа орбитальных параметров. В издании [8] показаны возможности расчета движения КА с использованием специализированного программного обеспечения моделирования траекторных параметров. В работах [9; 10] приведены современные подходы к более точному моделированию отдельных траекторных параметров движения КА. Стоит отметить, что в области расчета характеристик орбит КА ДЗЗ еще имеются пробелы в части комплексного анализа характеристик орбит, актуализирующих разработку соответствующей методики. Цель настоящей работы заключается в формировании универсальной методики, позволяющей на ранних стадиях проектирования достаточно быстро и одновременно точно проводить расчетные оценки параметров орбит КА ДЗЗ и анализировать условия их коррекции для поддержания этих параметров. Методика включает в себя: оценку максимального времени нахождения КА в зоне видимости ППИ и оценку зон видимости; обеспечение постоянства солнечной освещенности и возможности проведения наблюдений путем использования солнечно-синхронной орбиты (ССО); обеспечение видимости земной поверхности в заданное время полета КА, путем выполнения условий кратности орбит; оценку углового размера теневого участка орбиты; построение трассы полета КА и определение максимального просвета между географическими долготами восходящих узлов орбиты в зависимости от высоты; анализ влияния земной атмосферы на орбитальное движение КА; моделирование коррекций поддержания высоты орбиты с учетом атмосферного торможения КА. Для решения этих задач в рамках методики построена математическая модель условий функционирования КА ДЗЗ в целях анализа освещенности и видимости, кратности, покрытия земной поверхности трассами полета. При расчетах учтено воздействие атмосферы на орбитальное движение и выполнена коррекция поддержания высоты орбиты. Расчет выполнен с помощью программного обеспечения MatLab. Апробация методики проведена для КА «БелКА» [11]. Задача определения наклонения орбиты в зависимости от высоты. По замыслу КА «БелКА» выводится на ССО. Далее представлено решение задачи определения зависимости наклонения орбиты (i, град) от высоты (Н, км) и определения значения i для заданной высоты. Рис. 1. Зависимость наклонения от высоты орбиты где ω s = 0,2·10 -6 - средняя угловая скорость движения Земли вокруг Солнца, рад/с; R e = 6378,14 - экваториальный радиус Земли, км; μ = 3,99·10 5 - гравитационный параметр Земли, км 3 /с 2 ; J 2 = 1,08·10 -3 - безразмерный коэффициент, характеризующий форму Земли. Результаты расчета угла i (зависимость наклонения ССО от средней высоты полета КА) показывают, что свойство постоянства условий освещенности КА может быть достигнуто только для орбит, превосходящих по наклонению полярные орбиты (рис. 1). В пределах высот полета от 200 км до 2000 км изменение наклонения i составляет от 97 до 105°. Построение трассы полета и определение максимального просвета между географическими долготами восходящих узлов ССО. На примере КА «БелКА» была построена трасса полета, в географических координатах, в виде зависимости φ(λ). Вычисления проводились с помощью системы уравнений: где φ - географическая широта, град; λ - географическая долгота, град; ω 3 = 7,3·10 - 5 - средняя угловая скорость суточного вращения Земли вокруг оси, рад/с; r, x, y, z - вектор положения в орбитальной геоцентрической системе координат. Зависимость φ 1 (λ 1 ) приведена на рисунке 2. Далее решается задача по графическому построению зависимости максимального просвета между географическими долготами восходящих узлов ССО (λ) от высоты орбиты на временном интервале n суток (за целое число витков). Выбраны значения высот в интервале от 500 до 550 км. Строится зависимость максимального просвета от высоты орбиты (рис. 3). Определение размера тени. Зависимость углового размера тени (θ, град) определяется с учетом плоскости орбиты (β, град) для заданной высоты. Изменение половинного размера теневого участка: На рисунке 4 представлена кривая зависимости углового размера тени от угла наклона солнечных лучей к плоскости орбиты (β). Значения угла β для орбиты с теневым участком лежат в пределах: Орбиты с углами arcsin e e R RH β> + являются чисто солнечными, т.е. с отсутствием теневого участка. Для высоты орбиты с тенью Н = 511 км предельное значение угла β max = 67,8°. Оценка максимального времени нахождения КА на круговой орбите в зоне видимости ППИ и определения размера полосы обзора бортовой аппаратуры. Считается, что в некоторой точке сферической поверхности Земли расположен ППИ. Плоскостью местного горизонта является плоскость, касательная сферической поверхности в данной точке. Зона видимости (обзора) ППИ - это область пространства, в пределах которой с данного пункта виден КА. Практически область «радиовидимости» охватывает не все полусферическое пространство, а некоторую его часть, возвышающуюся над местным горизонтом на минимальный угол возвышения δ = δ min = 5-10°. Минимальный угол «радиовидимости» определяется исходя из условий распространения радиоволн и требований приема радиосигналов. Исходные данные: высота орбиты H = 511 км; минимальный угол места δ min = 5°, угол съемочной аппаратуры КА составляет ±40°. Размером поля зрения аппаратуры можно пренебречь. Геоцентрический радиус зоны обзора ППИ (φ δ , град) определяется по формуле Граница зоны видимости КА с ППИ определяется также (помимо минимального угла радиовидимости) максимальной дальностью связи (D max , км), обеспечиваемой аппаратурой ППИ и КА (рис. 5): Минимальный угол места, при котором работает антенна ППИ, составляет δ min = 5°. Максимальная дальность, на которой работает антенна, составляет D max = 2100 км и достигается при значении угла места δ D  5° (рис. 6). Из рисунке 5 находится соответствующий геоцентрический угол φ D  17°. Угловая скорость спутника равна 360°/98 мин, где за примерный период обращения КА берется 98 мин. Следовательно, спутник, проходящий через зенит ППИ, будет находиться в зоне радиовидимости примерно 2·17·98/360  9 мин. Если спутник проходит не через зенит ППИ, то время его нахождения в зоне видимости будет меньше. Рассчитано значение углового радиуса поля обзора по формуле (φ ε , град) Линейный размер (L, км) на поверхности Земли связан с центральным углом φ° соотношением L  111φ° = ±433 км. Проиллюстрированы границы зоны видимости на примере 3-х пунктов наблюдения, находящихся в городах: Отрадное, Берлин и Хабаровск. Границы зоны видимости для каждого пункта определяются (по времени) с помощью зависимостей: λ 1 (j), Δt 1 (j), λ 2 (j), Δt 2 (j) (на восходящей и нисходящей частях витка), j - количество точек (рис. 7). Определяется время полета (Δt, c) от начала витка до входа в зону видимости и выхода из нее: где u 1 , u 2 - аргументы широты прилета в точку на границе зоны видимости на восходящей и нисходящей части витка соответственно, град; ω КА - угловая скорость КА, град. Географическая долгота восходящего узла витка, проходящего через точку на границе зоны видимости: где Δλ - расстояние от географической долготы восходящего узла до географической долготы пункта по азимуту, град; dλ - смещение долготы восходящего узла за счет наклонения, град; ω 3 Δt - смещение долготы восходящего узла за счет поворота Земли за время полета от восходящего узла до точки, град. Моделирование коррекций. Во время движения по орбите КА подвергается воздействию: гравитационных полей Земли, Солнца, Луны, аэродинамическому напору, возмущениям от включений двигательной установки КА. Все перечисленные факторы, поддаются точному аналитическому расчету. Воздействие хотя и составляет малую величину (при плотности 10 -14 - 10 -13 кг/м 3 на высоте 600 км), но при постоянном действии приводит к существенному снижению высоты орбиты. Торможение КА и темп снижения высоты орбиты в значительной мере определяются уровнем солнечной активности, носящем случайный характер и оцениваемом по средней величине. Текущие же (случайные) значения высоты и темпа ее падения определяются по данным измерений орбиты. При этом темп падения высоты незначителен, что избавляет от необходимости принимать срочные решения по подъему высоты орбиты. Для поддержания свойств ССО достаточно отслеживать и регулировать величину смещения момента прихода Δt н в восходящий узел фактического движения от опорного. Задача заключается в том, чтобы определить, до какого уровня необходимо поднять высоту реальной орбиты (в соответствии с эквивалентными затратами ΔV на коррекцию), чтобы после подъема максимальное значение Δt не вышло за установленный предел Δt м . Ускорение с учетом торможения определяется по формуле где c б = cS m /2m - баллистический коэффициент КА, м 2 /кг; S m - площадь миделева сечения, м 2 ; m - масса КА, кг; c = 2,0-2,2 - коэффициент аэродинамического сопротивления; ρ - средняя плотность атмосферы на высоте H 0 , кг/м 3 . Процесс падения высоты и ее поднятие с помощью коррекций рассматривается за период времени: 0,5 года (180 суток). Число витков за 180 суток полета составляет: n = 2700. Падение высоты орбиты за сутки составляет 130 м. Допустимое значение падения высоты составляет 2000 м. Продолжительность падения на допустимую величину составляет 15 суток. Импульс необходимый для подъема на высоту Δy составляет 1,1 м/c. Необходимое число коррекций (подъемов) высоты орбиты 12 раз. Суммарный импульс: ΔV Σ = ΔVk = 1,1 · 12 = 13,2 м/c. Расход топлива: Далее строится зависимость периодической функции H(n), где диапазон изменения высоты орбиты лежит в пределах [-1-1] км (рис. 8). Основные параметры для построения: допустимый диапазон падения высоты 2000 м; продолжительность падения на диапазон высоты 15 суток; число коррекций 12. Расчетная формула: Заключение. Разработана комплексная методика расчета параметров орбиты для КА ДЗЗ с учетом анализа условий освещенности, кратности, покрытия и видимости поверхности Земли, позволяющая сократить временные затраты на этапе формирования технических предложений. В методике предложена система математических моделей, реализованная в виде алгоритмов в программном пакете MatLab, с помощью которой выполнена оценка параметров рабочей орбиты КА «БелКА» с учетом: постоянства солнечной освещенности; видимости земной поверхности в заданное время полета КА над одними и теми же участками обследуемой местности; глобального обзора Земли; наилучшей оперативности наблюдения районов земной поверхности; требуемых коррекций высоты орбиты, понижающейся вследствие торможения.

Об авторах

К В Михайловский

Московский государственный технический университет им. Н.Э. Баумана (национальный исследовательский университет)

Автор, ответственный за переписку.
Email: konst_mi@mail.ru

Михайловский Константин Валерьевич, кандидат технических наук, доцент кафедры «Ракетно-космические композитные конструкции» МГТУ им. Н.Э. Баумана.

2-я Бауманская ул., 5, стр. 1, Москва, Россия, 105005

М А Городецкий

Московский государственный технический университет им. Н.Э. Баумана (национальный исследовательский университет)

Email: mikgord@gmail.com

Городецкий Михаил Алексеевич, магистрант 2-го года обучения кафедры «Ракетно-космические композитные конструкции» МГТУ им. Н.Э. Баумана.

2-я Бауманская ул., 5, стр. 1, Москва, Россия, 105005

Список литературы