<< 4.1. Основные параметры Земли | Оглавление | 4.3. Геоцентрическая и геодезическая >>
4.2. Уравнение геоида
Для вычисления уравнения геоида используем формулу (4.1) и вычислим потенциал притяжения тела произвольной формы. Полученные выражения будут использоваться при вычислении прецессионного и нутационного движения осей Земли.
Потенциал точки с массой на расстоянии равен:
Сила притяжения этой точки другой точки с массой равна:
Подставляя в (4.5) выражение (4.4) и используя (4.3), получим хорошо известную формулу Ньютона:
Используем выражение для потенциала точки (4.4) и найдем потенциал притяжения элемента массы тела произвольной формы в точке ( - плотность, зависящая от координат элемента объема ). Для этого определим систему координат с началом в центре масс тела (рис. 4.1).
Будем считать, что точка расположена вне тела. Расстояние от точки до элемента массы , находящегося в точке с координатами , равно , а до точки равно . Тогда потенциал в точке равен
где - расстояние от точки до точки , - угол между отрезками и . Найдем разложение знаменателя по степеням до членов третьего порядка, предполагая, что :
Традиционное представление потенциала имеет вид:
Чтобы найти потенциал в точке от всех точек тела, проинтегрируем (4.6) во всему объему тела, используя разложение (4.7):
Первый член соответствует потенциалу точки с массой , расположенной в начале координат . Второй член равен нулю, так как начало координат совпадает с центром масс тела. Это легко доказать. В самом деле, центр масс тела - это точка , радиус-вектор которой относительно некоторой системы координат равен:
Определим осевые моменты инерции Земли относительно осей , соответственно, следующим образом:
Тогда третий член в (4.8) можно представить в виде:
Так как произведение равно перпендикуляру , опущенному из точки на прямую , то четвертый интеграл в (4.8) равен моменту инерции тела относительно оси . Следовательно, гравитационный потенциал потенциал тела произвольной формы в точке , расположенной вне тела, равен:
Выразим теперь момент инерции тела относительно оси через направляющие косинусы этой прямой относительно осей . Пусть . Тогда
Так как вектор имеет компоненты , то имеем согласно (2.17):
Отсюда
Интегралы
Таким образом, если координатные оси системы совпадают с главными осями тензора инерции, то
Обозначим угол между и осью через . Тогда
Для Земли главные моменты инерции равны кг$·$ м$^2$, кг$·$ м$^2$, кг$·$ м$^2$, кг$·$ м$^2$. Относительная разность экваториальных моментов и равна . Поэтому часто считают, что и тем самым полагают, что Земля - двухосный эллипсоид или эллипсоид вращения. Гравитационный потенциал Земли называется геопотенциалом.
Перепишем (4.10) в следующем виде:
где называется динамическим форм-фактором Земли, - экваториальный радиус Земли ( м). Приведенные значения взяты из стандартов Международной службы вращения Земли (IERS2003).
Уравнение (4.11) в пределе должно выполняться на поверхности Земли. Если ограничиться разложением до , то полный геопотенциал на поверхности Земли равен:
где - центробежный потенциал. Запишем его в виде:
Формула верна с точностью до .
Чтобы найти уравнение поверхности геоида, перепишем формулу (4.12) в следующем виде:
Так как с ошибкой , то перемножая скобки и пренебрегая членами, содержащими произведения малых величин, запишем уравнение геоида в виде
В первом приближении уравнение геоида представляет фигуру, близкую к эллипсоиду с геометрическим сжатием .
В самом деле, рассмотрим меридиональное сечение двухосного эллипсоида, которое является эллипсом с большой и малой полуосями:
Так как малая полуось равна и , , то из (4.15) получим:
С точностью до первого порядка сжатие равно: .
В настоящее время коэффициент определяется по наблюдениям искусственных спутников Земли. Забегая вперед, скажем, что по скорости прецессии определяется динамическое сжатие Земли
Найдем теперь ускорение силы тяжести на геоиде. Для этого достаточно продифференцировать полный геопотенциал (4.12). Так как , то с учетом (4.3) найдем, что модуль равен:
Отсюда ускорение силы тяжести на экваторе равно:
Знак минус в формуле (4.17) подчеркивает тот факт, что радиальная компонента направлена в противоположную сторону, чем единичный вектор . В соответствии с принятыми значениями постоянных равно (см. табл. 9.2): . Используя определение (4.17) и (4.16), найдем ускорение силы тяжести на широте :
Более точное выражение для ускорения силы тяжести , измеряемого в точке с геодезической широтой на высоте над эллипсоидом с полуосями и сжатием , имеет вид:
где
где высота измеряется в метрах.
<< 4.1. Основные параметры Земли | Оглавление | 4.3. Геоцентрическая и геодезическая >>
Публикации с ключевыми словами:
астрометрия - сферическая астрономия - системы координат - шкалы времени
Публикации со словами: астрометрия - сферическая астрономия - системы координат - шкалы времени | |
См. также:
Все публикации на ту же тему >> |