DEVELOPMENT OF METHOD FOR DETERMINING AND CORRECTING PARAMETERS OF THE WORKING ORBIT OF THE EARTH REMOTE SENSING SATELLITE

Cover Page

Abstract


A technique for the complex analysis of the characteristics of orbits used for space vehicles for remote sensing of the Earth (RSE), taking into account their trajectory, is proposed. The results of mathematical modeling of the operating conditions of such devices are presented. The technique is tested taking into account the input parameters of the orbit of the BelKA spacecraft. The technique will be useful in choosing of complex orbital characteristics for RSE spacecraft at the stage of technical proposals.


Введение. Одной из наиболее важных задач Федеральной космической программы России на период до 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, с помощью которой выполнена оценка параметров рабочей орбиты КА «БелКА» с учетом: постоянства солнечной освещенности; видимости земной поверхности в заданное время полета КА над одними и теми же участками обследуемой местности; глобального обзора Земли; наилучшей оперативности наблюдения районов земной поверхности; требуемых коррекций высоты орбиты, понижающейся вследствие торможения.

K V Mikhaylovskiy

Bauman Moscow State Technical University (National Research University)

Author for correspondence.
Email: konst_mi@mail.ru
2-nd Baumanskaya str., 5/1, Moscow, Russia, 105005

Mikhailovskiy Konstantin Valerievich, candidate of technical sciences, associate professor of the department “Rocket and space composite construction”. Bauman MSTU.

M A Gorodetsky

Bauman Moscow State Technical University (National Research University)

Email: mikgord@gmail.com
2-nd Baumanskaya str., 5/1, Moscow, Russia, 105005

Gorodetsky Mikhail Alekseevich, master student of the 2nd year of the department “Rocket and space composite construction”. Bauman MSTU.

  • Montenbruck O. and Gill E. Satellite orbits: models, methods, applications. Springer — Verlag Berlin Heidelberg, 2000. 369 p.
  • Hastings D. and Garrett H. Spacecraft-Environment interactions. Cambridge University press, 2004. 320 p.
  • Davis V.A. and Gordon L.W. Spacecraft surface charging handbook. NASA, technical report, 1992. 343 p.
  • Hughes P.C. Spacecraft attitude dynamics. Dover publications, Mineola, 2004. 592 p.
  • Schaub H. and Junkins J.L. Analytical mechanics of space systems, Second edition. AIAA, Reston, 2009. 794 p.
  • Khodnenko V.P., Khromov A.V. Selection of the design parameters of the orbit correction system for the space remote sensing satellite // Questions of electromechanics. 2011. Vol. 121. P. 15—22.
  • Asghar E., Mehran M., Ali A. Orbit analysis of a remote sensing satellite for local observation of the earth surface // Applied Science and Research Association (ASRA) Department of Mechanical Engineering, Toosi University of technology, 2005. 5 p.
  • Lovera M. Control-oriented modelling and simulation of spacecraft attitude and orbit dynamics // Journal of Mathematical and Computer Modelling of Dynamical Systems, Special issue on Modular Physical Modelling. 2006. Vol. 12. No. 1. P. 73—88.
  • Sidorenko V.V. and Neishtadt A.I. Investigation of the stability of long-periodic planar motion of a satellite in a circular orbit // Cosmic Research. 2000. Vol. 38. No. 3. P. 289—303.
  • Casella F., Lovera M. High-Accuracy orbital dynamics simulation through keplerian and equinoctial parameters / The Modelica Association, Modelica, March 3rd—4th. 2008. P. 505—514.
  • Site: https://www.energia.ru/rus/news/news-2006/press_release-07-26.html (Data: 26.07.2006).

Views

Abstract - 147

PDF (Russian) - 86


Copyright (c) 2017 Mikhaylovskiy K.V., Gorodetsky M.A.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.