Влияние возмущающих гравитационных сил, связанных с нецентральностью гравитационного поля Земли, на эволюцию орбиты космического аппарата

Размещено на http://www.allbest.ru/

Влияние возмущающих гравитационных сил, связанных с нецентральностью гравитационного поля Земли, на эволюцию орбиты космического аппарата

А.В. Ващенко

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

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

При баллистических расчетах исследуется эволюция околокруговых орбит, учитывается влияние возмущающих гравитационных сил, связанных с «неправильностью» формы Земли и возмущающие аэродинамические силы от сопротивления земной атмосферы. Исследуются орбиты высотой 600 км и наклонением орбиты i от 0 до 180°. Расчет производится с помощью программы баллистических расчетов (на базе математической системы компьютерной алгебры Maple 7).

В расчетах используется система дифференциальных уравнений в оскулирующих элементах с модификацией - вместо времени прохождения перицентра используется аргумент широты. Данная система может использоваться для анализа движения любого типа (эллипс, парабола, гипербола), но она имеет особенности при е=0 и i=0. Для исключения особенности е=0 используются две новых формальных переменных (составляющие вектора Лапласа): l(t)=e(t)sin(щ(t)) и q(t)=e(t)cos(щ(t)).

Метод оскулирующих элементов сводится к тому, что исследование возмущенной траектории КА может быть сведено к анализу совокупности невозмущенных траекторий, соответствующих каждому моменту времени. Т.е., при анализе возмущенного движения мы можем считать, что в любой момент времени КА находится на дуге конического сечения с определенными значениями элементов орбиты, а коническое сечение при движении КА изменяется. При этом изменяются элементы орбиты [1]:

космический земля возмущенный гравитационный

p=p(t); e=e(t); щ=щ(t); Щ=Щ(t); i=i(t); u=u(t).

Элементы орбиты, рассматриваемые как функции времени, через которые координаты и составляющие скорости в возмущенном движении выражаются теми же формулами, что и в невозмущенном движении, называются оскулирующими элементами [1].

Рисунок 1 Элементы возмущенного движения

Обозначения: p(t) - параметр, e(t) - эксцентриситет, щ(t) - агрумент широты перицентра, ?(t) - долгота восходящего узла, i(t) - наклонение, u(t) - аргумент широты, A(t) - большая полуось, r(t) - радиус, f(t) - истинная аномалия, _S(t), _T(t), W(t) проекции возмущающих ускорений на оси орбитальной системы координат (по радиус-вектору, по нормали к нему в плоскости орбиты и по нормали к плоскости орбиты).

После всех преобразований система дифференциальных уравнений движения в оскулирующих элементах при исключении особенностей е>0 принимает вид:

Данная система дифференциальных уравнений используется при исследовании эволюции орбиты КА в программе баллистических расчетов.

Возмущения от гравитационного поля Земли.

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

в = (rc-Rз.п.)/Rз.э.?0,0033523,

где rc - средний экваториальный радиус, равный 6378,245 км.;

Rз.п. - полярный радиус, равный 6356,863 км.

Потенциал гравитационного поля Земли как сжатого сфероида равен:

,

где fm - гравитационный параметр Земли, r - расстояние от центра Земли до КА, - широта точки, = 66,07*10^3 км^2 [1].

Однако для точных расчетов орбиты спутника необходимо учитывать аномалии земного потенциала. Соответственно потенциал гравитационного поля Земли может быть записан, как [1]:

где rc - средний экваториальный радиус Земли, r,,л - геоцентрический радиус, геоцентрическая широта и долгота, - полином Лежандра i-го порядка, - присоединенная функция nk-го порядка, выраженная через k-ю производную полинома Лежандра n-го порядка, , , - постоянные, характеризующие фигуру Земли.

Зональные составляющие:

Потенциал гравитационного поля Земли для зональных составляющих принимается в форме [1]:

Рисунок 2 Положительные и отрицательные значения зональной гармоники для i=4

Тессериальные и секториальные составляющие:

Потенциал гравитационного поля Земли для тессериальных и секториальных составляющих принимается в форме [1]:

Подставляя соотношение для долготы л=l+Щ-щ3t в выражение для тессериальных и секториальных составляющих получим:

Рисунок 3 Положительные и отрицательные значения тессеральной гармоники для n=10 и k=4.

Рисунок 4 Положительные и отрицательные значения секториальной гармоники для n=10

Используя геометрические соотношения sin(ц)=sin(u)sin(i)

,

представим силовую функцию в виде:

UE(r,u,i,Щ,t)=Uno(r,u,i)+Unk(r,u,I,Щ,t)

И проекции ускорений на соответствующие направления через элементы орбиты:

, ,

Постоянные силовой функции (потенциала) принимаются для fm=398603,2 rc=6378,165.

По данным из [1] и [2]:

Коэффициенты In0 зональных гармоник потенциала

Коэффициенты Ink тессериальных и зональных гармоник потенциала

Фазы лnk тессериальных и зональных гармоник потенциала

Результаты баллистических расчетов

На рисунке 5 показана зависимость изменения долготы восходящего узла в зависимости от наклонения орбиты КА.

Характеристики космического аппарата следующие:

Масса КА на опорной орбите 1950 кг.

Наклонение орбиты от 0 до 180°

Высота орбиты 600 км

Характерная площадь аппарата 20 м2

Рисунок 5 Зависимость изменения долготы восходящего узла от наклонения орбиты КА

Как видно из графика угол долготы восходящего узла изменяется, то есть происходит поворот плоскости орбиты вокруг оси Земли, называемый прецессией узла орбиты. С наклонением орбиты близким к 90° прецессия орбиты уменьшается. Этот свойство орбиты необходимо учитывать при анализе обеспечения КА энергией, получаемой от солнечных батарей. При больших изменениях угла долготы восходящего узла КА будет чаще оказываться в тени Земли. Соответственно можно выбрать орбиту так, чтобы КА всегда был освещен Солнцем. В таком случае орбита КА будет поворачиваться каждые сутки вместе с Землей на ?1°.

На рисунках 6 и 7 показаны зависимости изменения большой полуоси и долготы восходящего узла от массы аппарата и наклонения орбиты.

Рисунок 6 Зависимость изменения большой полуоси от массы КА и наклонения орбиты

Рисунок 7 Зависимость изменения долготы восходящего узла от массы КА и наклонения орбиты

Список использованной литературы

1. Дубошин Т.Н. Справочное руководство по небесной механике. М.: Изд. «Наука», 1971. 431 с.

2. Дубошин Т.Н. Небесная механика (основные задачи и методы). М.: Изд. «Наука», 1975. 385 с.

3. Константинов М.С., Каменков Е.Ф., Перелыгин Б.П., Безвербый В.К. Механика космического полета. М.: Машиностроение, 1989. 408 с.

4. Гушин В.Н. Основы устройства космических аппаратов. М.: Машиностроение, 2003. 272 с.

Размещено на Allbest.ru