Нецентральность гравитационного поля Земли - Возмущенное движение космического аппарата

Возмущенный движение гравитационный орбита

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

Потенциальная функция поля тяготения такого тела приближенно описывается формулой, рекомендованной Международным астрономическим союзом для практического использования, т. е. стандартной формой записи потенциала сил притяжения Земли:

карта гравитационного поля земли

Рисунок 2 - Карта гравитационного поля Земли

типы гармоник потенциала поля тяготения земли и расположение на сфере оластей изменения их знаков

Рисунок 3 - Типы гармоник потенциала поля тяготения Земли и расположение на сфере оластей изменения их знаков

Если рассматривать фигуру Земли как сжатый сфероид, что оказывается достаточным для решения большинства практических задач, то можно использовать формулу:

Первое слагаемое в последнем выражении есть ньютоновский потенциал. Второе слагаемое учитывает сжатие Земли по полюсам.

Используя формулы сферической тригонометрии для геоцентрической широты ц' получаем выражение:

Где U - аргумент широты, а I - наклонение орбиты. Обозначая, получаем для потенциала от эффекта сжатия земного сфероида следующее выражение:

Для того чтобы найти проекцию гравитационного ускорения от сжатия земного эллипсоида на какое-либо направление достаточно найти производную от потенциала по этому направлению. Нас в дальнейшем будут интересовать направления, задаваемые орбитальной системой координат S, T,W. Поэтому находим:

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

Для исследования влияния возмущений, задаваемых в форме приведенной выше (S1, T1, W1), используют систему дифференциальных уравнений для оскулирующих элементов, где независимой переменной является U. Подстановка возмущающих ускорений S1, T1, W1 в систему приводит к системе:

Для любых орбит множитель г никогда не отклоняется от единицы более чем на ?0,003. Такого же порядка ошибки могут произойти вследствие осреднения значений оскулирующих элементов Щ...ф. Поэтому не имеет смысла усложнять задачу учетом г, несмотря на то, что в первом приближении этот множитель нетрудно учесть. Примем г =1 за исключением случаев, оговоренных особо.

При интегрировании 1-го, 3-го, 5-го уравнений за один оборот ИСЗ (в пределах от "u0" до u0+2p ) получаем, что вековые уходы элементов P, i, e в первом приближении отсутствуют. Указанные элементы подвержены только периодическим возмущениям.

Таким образом, Фокальный параметр орбиты спутника не имеет вековых возмущений из-за нецентральности гравитационного поля Земли, принимаемой сжатым сфероидом. В течение одного оборота вокруг Земли он (с точностью до эксцентриситета "е") испытывает два полных колебания с амплитудой:

Периодическое возмущение Эксцентриситета (с точностью до "е") происходят с периодом, равным периоду обращения спутника вокруг Земли (две гармоники), и с периодом втрое меньшим. А вековые отклонения отсутствуют. Амплитуда колебаний зависит от i и от щ. Чем меньше значение эксцентриситета, тем относительно больше он колеблется.

Наклонение орбиты в первом приближении, как и параметр орбиты "p" не имеет вековых возмущений. За один оборот параметр "i" так же, как и "p", совершает (с точностью до "e") два полных колебания с амплитудой:

Возмущения Долготы восходящего узла Щ и Аргумента перигея щ кроме периодических имеют вековые уходы.

Периодические возмущения линии узлов имеют второй порядок малости по отношению к вековым возмущениям, поэтому график функции dW = f (u) (по крайней мере, за один виток) представляет собой почти прямую линию.

Для прямых спутников, имеющих узел двигается в сторону обратную возрастанию Щ (к западу), для обратных спутников - (сторону возрастания Щ (к востоку). Это движение определяет поворот плоскости орбиты вокруг оси Земли и называется Прецессией узла орбиты.

Изменение оскулирующих элементов в возмущенном движении приводит к возмущению радиуса орбиты, следовательно, и высоты полета. Эти возмущения таковы, что происходит как бы частичное "отслеживание" поверхности Земли высотой полета, при этом аппарат "поднимается" над экваториальными областями и "проседает" над полюсами.

Пространственный поворот плоскости орбиты при возмущенном движении, связанный с изменениями элементов I И Щ, Приводит к появлению бокового "ухода" возмущенной орбиты (по отношению к невозмущенной).

Таким образом, получены формулы учитывающие сжатие земного сфероида для элементов эллиптических орбит P, Щ, i, щ, e в зависимости от аргумента широты U.

В этих формулах в скобках - невозмущенные элементы орбиты, а значения Dp, dе, dw, dЩ, di как функции "u" определены ниже:

Для данного значения U и соответствующих ему элементов орбиты, рассчитанных по этим формулам, можно вычислить координаты x, y, z и скорости.

Похожие статьи




Нецентральность гравитационного поля Земли - Возмущенное движение космического аппарата

Предыдущая | Следующая