Серия «Естественные науки»
Гравитационный аналог Земли и ее фигура
к.т.н. доц. Русанов П.Г.
МГТУ им. Н.Э. Баумана 8 (499) 263-63-75, pa4rus@gmail.com
Аннотация. Методом численного моделирования фигуры Земли установлено, что ее реальная фигура не соответствует центральному гравитационному полю. Предложен двухмассовый аналог источника нецентрального поля гравитации Земли в виде ядра и экваториального кольца. Установлены численные значения масс двухэлементой модели и радиуса кольца, при которых расчетная фигура Земли и зависимость ускорения силы тяжести на ее поверхности от широты соответствуют экспериментальным данным и международным соглашениям.
Ключевые слова: фигура Земли, масса, гравитация, вращение Земли, ядро, кольцо, равновесная форма поверхности, ускорение поля силы тяжести.
Введение
Фигурой Земли принято считать равновесную форму поверхности спокойного мирового океана, покрывающего, условно, всю планету. На форму этой поверхности влияют, главным образом, силы гравитации Земли и ее суточное вращение. Очевидно, что при однородном строении планеты и в отсутствии суточного вращения фигура Земли была бы близка к сферической. Экспериментальные данные свидетельствуют о том, что реальное внешнее гравитационное поле Земли не является центральным [1 - 4]. Главной причиной отклонений гравитационного поля Земли от его расчетного центрального образа является неоднородность распределения масс внутри планеты.
Одна из целей исследователей фигуры Земли - приблизить расчетные модели к ее реальному образу. Существующие аналитические модели фигуры Земли построены на аппроксимации формы всей планеты в целом или ее локальных участков подходящей поверхностью эллипсоида по экспериментальным данным, не используя уравнений равновесного состояния водной поверхности мирового океана.
Анализ фигуры Земли и ее гравитационного поля выполнен в рамках твердотельных моделей на основе конструктивного моделирования масс - источника гравитационного поля и численного решения дифференциального уравнения (ДУ) относительного равновесия частицы жидкости на поверхности мирового океана с учетом физических сил и сил инерции.
В ходе исследований установлено, что результаты численных расчетов характерных размеров фигуры Земли для модели № 1 с центральным гравитационным полем не отвечают экспериментально измеренным величинам Re, Rp - экваториального и полярных радиусов
планеты.
В связи с этим, для уточнения внутреннего строения Земли предложена комбинированная модель источника гравитационных сил планеты № 2, состоящая из двух тел - центрального ядра и экваториального кольца, которая позволила существенно повысить точность представления фигуры Земли и параметров поля результатов. Для двух-массовой модели источника гравитации найдены оптимальное соотношение масс ядра и кольца, а также величина радиуса кольца, при которых получаемые результаты приближаются к экспериментальным данным и к критерию международных соглашений в отношении: 1) расчетного образа фигуры Земли и 2) величины ускорения поля силы тяжести на ее поверхности [4].
Постановка двух задач и допущения
Задача оценки фигуры Земли рассматривается ниже в двух вариантах ее постановки. В первом варианте предполагаем, что внешнее гравитационное поле планеты вблизи ее поверхности является центральным, сферическим, при этом основным источником гравитации является однородное сферическое ядро планеты. Во втором варианте гравитационное поле
считаем симметричным относительно ОС - оси суточного вращения планеты и плоскости экватора. Здесь источником гравитации служит система из двух однородных симметричных тел с общей осью симметрии, состоящая из ядра и экваториального кольца (пояса).
В обоих случаях неизвестную фигуру Земли считаем поверхностью тела вращения, ось которой совпадает с осью Ог. Дополнительным элементом симметрии этой поверхности служит ее плоскость экватора. Вследствие принятых допущений о принадлежности фигуры Земли к числу поверхностей тел вращения с двумя элементами симметрии, задача расчета ее формы сводится к расчету плоской формы ее меридиана от экватора до полюса. Анализ формы меридиана в обоих вариантах задачи проведем с применением цилиндрической и декартовой систем координат (СК), жестко связанных с Землей, с общим началом координат в т. О, центре ядра, и общей осью Ог (рисунок 1).
Искомую линию меридиана фигуры Земли в плоскости Оху рассматриваем как геометрическое место равновесных положений на поверхности океана частиц воды, принимаемых за материальные точки некоторой массы т. Неизвестную линию меридиана представим в цилиндрической СК функцией г = г(ф), где ф - географическая широта места. В каждом из двух модельных вариантов источника гравитационного поля Земли, рассматриваемых ниже, эту функцию будем рассчитывать исходя из условия покоя частиц воды по отношению к Земле.
Считаем известными величины О - гравитационная постоянная в законе всемирного тяготения, М - масса Земли, О - постоянная угловая скорость суточного вращения Земли, Яе, Яр - экваториальный и полярный радиусы Земли, gО - ускорение силы тяжести на экваторе, gе. - ускорение гравитационного поля Земли на экваторе. При этом gО = ge — О2Яе и
gA2 = ом .
Для частицы воды, расположенной на равновесной поверхности мирового океана Земли, пренебрегаем силами поверхностного натяжения, силами гравитации со стороны соседних частиц воды, движением отдельных масс Земли относительно ядра и влиянием космических тел.
Математическая модель № 1 источника гравитации Земли в виде одиночного ядра
В этой модели основным источником гравитационных сил Земли считаем однородное сферическое ядро. В число сил, действующих на частицу воды массы т, включаем (рисунок 2) Г - силу гравитации со стороны ядра Земли:
mM R2e
F = G-^ = mge—e ; (1)
2 o e 2 '
r r
Р - силу воздействия атмосферного воздуха и А - силу Архимеда со стороны соседних частиц воды.
Полагаем, что в состоянии относительного равновесия масс океана по отношению к ядру Земли силы А и Р перпендикулярны к свободной поверхности воды.
В этих условиях относительному равновесию точечной массы m в равномерно вращающейся системе отсчета движения, жестко связанной с Землей, отвечает уравнение динамики в проекции на местное касательное направление к меридиану:
(F — Ф cos j)sin a-Ф sin j cos a = 0, (2)
где: r [м], j [рад] - полярные координаты точечной массы m;
Ф = mW2r cos j - величина переносной силы инерции Эйлера в неинерциальной системе
отсчета, вращающейся вместе с Землей с постоянной угловой скоростью W; a [рад] - угол между касательным направлением к меридиану и окружным направлением в плоскости меридиана (рисунок 2), при этом по геометрическим соображениям < = r ■ tga •
Рисунок 1 Рисунок 2
С учетом развернутых выражений для Ф, rj представим (2) после сокращения на множитель m как ДУ первого порядка относительно r = г(ф) :
dr _ W2r4 sin ф cos ф
~J = Ü2 r\2 3 2 . (3)
dф geRe -W r cos ф Введением величин р и e уравнение (3) можно было бы представить в безразмерном
виде:
r = ReP ,
W2 Re dp sp4sin ф cos ф
ge dj 1 -sp cos j Тем не менее, для удобства физического восприятия результатов, получим численное решение r = r(j) непосредственно по уравнению (3) с исходными данными [4]: W =7.292115
5 1 2
10- [с- ], ge = 9.814234780 [мс- ], при очевидном начальном условии j = 0, r = Re= 6378245 [м]. Отметим, что на нулевой широте j = 0 значение параметра gQ отвечает общепринятой
2 2
величине ускорения поля силы тяжести на экваторе gW= ge - W2Re = 9.7803185 [мс- ].
Рисунок 3. Зависимость^^ =АД(ф>)
Для наглядности полученные результаты численного решения (3) г = г(ф) представим зависимостью АЯ =АЯ(ф) (рисунок 3), где ф - широта в [град], АЯ = Яе — Г [м] - величина
снижения расстояния от текущей точки меридиана до центра ядра при возрастании угла ф от 0 до 90°.
e
Согласно расчетам величина полярного радиуса Земли оказалась на 11.002 [км] меньше экваториального радиуса, что почти в 2 раза меньше реальной разности указанных радиусов по данным космических измерений:
АЯ0 = Я — Яр = 21.382 [км].
Поэтому делаем вывод, что в рамках гравитационной модели № 1 точность определения формы фигуры Земли в полярных областях планеты недостаточна, и в связи с этим необходима дополнительная коррекция модели источника гравитации Земли.
Математическая модель № 2 двухмассового источника гравитации
Не имея достаточной точности, результат, полученный выше для АЯ величины снижения расстояния от текущей точки меридиана до центра ядра фигуры Земли с центральным полем гравитации, тем не менее, является качественно верным в отношении условия Яе > Яр. Он указывает на то, что у планеты, целиком покрытой слоем жидкости (воды) и обладающей центральным полем гравитации, вследствие суточного вращения должно наблюдаться частичное перераспределение массы жидкости в область экваториального пояса. В свою очередь, изменённая структура масс планеты должна нарушить центральность ее внешнего поля гравитации.
Кроме того, по мнению ученых [1, 3], некогда ранее поверхностный слой Земли был горячей расплавленной массой и Земля имела в 2 - 3 раза большую угловую скорость суточного вращения. Согласно результатам, полученным для модели № 1, вследствие суточного вращения у фигуры еще горячей Земли размер ее экваториального радиуса Яе также превышал Яр. Общая конфигурация когда-то расплавленной планеты должна была естестсвенным образом сохраниться по мере ее остывания и затвердевания части жидкой массы. В результате прошедших геофизических явлений современная фигура Земли стала близкой к эллипсоиду вращения, приплюснутому со стороны полюсов.
Принципиальная «приплюснутость» формы Земли к плоскости экватора настолько сильно сказалась на ее внешнем гравитационном поле, что для ряда современных технических приложений это поле уже не допустимо считать центральным. Факт нецентральности поля гравитации Земли уже давно считается неоспоримым и экспериментально проверенным. Тем более, что к настоящему времени экспериментально зафиксированы более «слабые» локальные отступления (аномалии) гравитационного поля Земли от принятых стандартов.
Воспользуемся этими доводами для организации модифицированной модели № 2 масс - источника поля гравитации Земли. С учетом вышесказанного, в новой модели представим источник гравитации Земли в виде комбинации из двух тел, состоящей из центрального сферического ядра массы Мя и однородного экваториального кольца массы Мк, радиуса Як (рисунок 4).
Полагаем, что кольцо и ядро расположены концентрично, Мя + Мк =М, при этом Мя =^М, Мк = тМ л+т=1. При заданных величинах Яе, О, ge значения двух других параметров, т, Як, далее подлежат выбору из условий наилучшего соответствия двух расчетных величин:
1) полярного радиуса фигуры Земли - значению Яр = 6356863 [м],
2) функции gО(ф) - ускорения свободного падения зависимости go = go(ф) для «нормальной Земли», принятой в [4]:
g0 = 9.7803185(1 + 0.0053024 мп2ф - 0.0000059 Бш22ф) [м с-2]. (4)
Модель № 2 имеет иной источник гравитации. Это отразится лишь на физической картине гравитационных сил, действующих на частицу воды массы т. А именно, гравитационное воздействие на массу т далее будем оценивать двумя силами (рисунок 4).
Одна сила - от ядра:
т • М Я2е
Гя =Л-О--— = Л-т • ge —г, (5)
г г
а вторая сила - от кольца.
Оценку силы выполним, следуя методу твердых тел [5]. Представим кольцо массы Мк системой из 2п одинаковых масс - однородных шаров массыМ0 =0.5Мк/п = 0.5 цМП, центры которых расположенны равномерно по периметру кольца и симметрично по отношению к плоскости исследуемого меридиана (рисунок 5).
P
Рисунок 4 Рисунок 5
Силу Ек считаем равнодействующей сходящейся системы элементарных сил, действующих на массу т со стороны каждого шара. Для упрощения расчетов воспользуемся свойствами симметрии в расположении дискретных масс шаров по отношению к плоскости меридиана Оху.
Очевидно, что т/к - равнодействующая двух сил гравитационного притяжения точечной массы воды т к к-той паре симметрично расположенных шаров (к =1, п , отсчет к начинается с пары шаров, ближайших к меридиану) - будет лежать в плоскости меридиана, то есть в плоскости Оху.
Поэтому проекции силы /к на оси Ох и Оу будут равны удвоенным значениям соответствующих проекций сил, ее составляющих.
FKxy = 2шЁ fk,y ,
(6)
k-1
, „Mlx R2J, где: fx = —f - = -mge
Ml_
e'kx . Г _
k- ; Jky
= -mge
Relky .
nlk3 nlk'' ' nlk3 6 nlk3
hz = yj^L + lfy +- расстояние между центром шара k -той пары и точкой массы m;
lx =r cos ф-rk cos Pk; =r sin ф ; = r sin Pk;
bk = pk/n - угловой шаг расположения центров шаров по осевой окружности кольца;
fkr - fx cos j + fy sinф ; fp - fy cosФ-fx sinф - проекции вектора Jk на радиальную и
поперечную оси в плоскости меридиана. На этот раз, с учетом (6), относительному равновесию точечной массы m будет отвечать уравнение динамики (7) в проекции на местное касательное направление к меридиану:
-2m^i fkr -Фcos j)sina + (2m^i fp -Фsin j)cosa - 0 .
(7)
k-1
k-1
Подобно уравнению (2), (7) тоже допускает сокращение на множитель m и представле-
ние его в виде ДУ первого порядка относительно г = г(ф) - функции, определяющей форму меридиана:
42,
q2г sin ф cos ф- fp
dr „3 k=i p
dj= r -n-"-• (8)
dj hgR - 2r2 -Q2r3cos2 j
k=1
Напомним, что в ДУ (8) содержатся два свободных параметра т, ЯК, величины которых еще предстоит выбрать так, чтобы численное решение (8) имело результаты, близкие к экспериментальным значениям Яе, Яр и принятому соглашению [4]. Проблему подбора параметров т, ЯК можно было бы решить методами минимизации двух невязок 81 и 82:
e = r.
RP и s2=ж! g (r, Ф) - ф
1 j=p/2 p °2 _ 0<ф<п/2
где: r = r(j, m, RK) - численное решение ДУ (8) с начальным условием ф = 0, r = Re,
I R2 " n
g(r, ф) =, [hge Rr - 2^ fkr - r • (Qcos ф)2 ]2 + [2^ fkp - r • Q2 cos jsin j)2 ]2 •
V r k=1 k=1
Тем не менее, в связи с качественным отношением к получаемым результатам данного исследования, а также ввиду громоздкости представления функций, входящих в невязки, оценку оптимальных значений параметров h, Rk выполним методом сведения данной задачи оптимизации к двум одномерным случаям.
Численные результаты
На первом этапе анализа модели № 2 для фиксированных значений h, начиная с h=0.1=10%, при 2и=100, численными итерационными методами были подобраны соответствующие h значения параметра Rk из условия близости решения ДУ (8) в конце участка интегрирования по ф, то есть значения r(p/2,h, Rk) к величине Rp = 6356863 [м].
В результате расчетов установлено, что граничному условию на конце участка интегрирования удовлетворяют следующие пары параметров h, RK (таблицы 1а, 1b, 1с).
Таблица 1а
h % 10 20 30 40 50 60 70 80 90
Rk [км] 314.4 333.5 356.6 385.2 422.0 471.9 545.0 667.9 946.1
Таблица 1 b
h % 91 92 93 94 95 96 97 98 99
Rk [км] 997.7 1059 1132 1224 1342 1503 1740 2140 3054
Таблица 1с
h % 99.1 99.2 99.3 99.4 99.5 99.6 99.7 99.8
Rk [км] 3223 3422 3659 3948 4307 4759 5320 5942
По таблице 1 видно, что с уменьшением массы кольца (л®1, М®0) резко возрастает значение его радиуса ЯК (ЯК® Яе), при котором г = г(ф, т, ЯК), решение ДУ (8) может удовлетворять граничным условиям на экваторе ф=0 и на полюсе ф=р/2.
Отсюда следует вывод, что экваториальный поверхностный слой планеты, несмотря на его сравнительно малую массу, оказывает существенное влияние на ее гравитационное поле.
В каждом случае подбора величины ЯК, обеспечивающей схождение численного решения ДУ (8) к значению г(р /2, ЯК) = КР при фиксированном значении производим повторное численное интегрирование (8). Но на этот раз с целью оценки А? = оф/2|g (Л ф) - ?о(ф)|
- максимальной погрешности разности, ?(г, ф) - расчетного значения ускорения свободного падения на поверхности фигуры Земли в широтах из промежутка 0 < ф < 90 [град], отличающихся на 1 [град], и значениями ?0 = ?0(ф) для тех же широт.
Расчеты показали, что величина Д^ имеет глобальный минимум в диапазоне 80% < л £ 90%. Выборочные результаты подобных расчетов приведены в таблице 2. _Таблица 2
Л % 82 83 84 85 86 87 88 89 90 91
Як [км] 682.3 701.1 721.5 743.8 768.3 795.4 825.6 859.5 897.9 945.3
Дg [мм с-2] 6.941 6.928 6.914 6.899 6.883 6.867 6.853 6.850 6.851 6.854
Согласно таблице 2, предложенная математическая модель №2 дает наиболее близкие результаты к [4] при л » 89% и ЯК » 860 [км] с погрешностью оценки величины g(r, ф) -ускорения поля силы тяжести менее 0.07% на всем протяжении меридиана фигуры Земли.
Об отличии форм расчетной линии меридиана при оптимальных значениях л и ЯК от меридиана фигуры Земли в виде эллипсоида вращения [1 - 3]:
2 2 X У 1
--ь — = 1
22 а Ь
с полуосями а= Ке и Ъ= Кр можно судить по рисунку 6.
Рисунок 6 Рисунок 7
На нем изображена зависимость от широты ф [град] величины Дг(ф) - расстояния
между этими кривыми, измеряемого в [м] по радиальному направлению из центра О: Дг(ф) = г (ф, л, ЯК) - гЭ (ф), где гЭ (ф) = ^ (а соб у)2 + (Ь Бт у)2; у = аг^^ф • а / Ь). На широте ф =350 функция Дг(ф) имеет максимальное значение, равное 90.9 [м]. При
этом все промежуточные точки меридиана расчетной фигуры Земли для модели №2 отстоят дальше от центра О, чем точки меридиана у эллипсоида. Это означает, что эллипсоид Земли [1-3] расположен внутри поверхности фигуры Земли, определяемой по модели №2, и касается ее лишь в точках экватора и на полюсах.
На рисунке 7 отражена зависимость от широты ф [град] а - угла отклонения в [дуг. минутах] направления местного вектора ускорения поля силы тяжести g(r, ф) в плоскости меридиана, или местной вертикали от направления на центр О.
tga =
X1¥ - г °2 соб ф ^ ф
к=1
/
Я
л^-г1кг- г • соб Ф)2 .
В районе широты ф =440 угол а достигает наибольшего величины, равной 11.544 [дуг.мин.], что соответствует общепринятой норме отклонения данного направления от радиального.
Заключение
Численными расчетами подтверждено, что математическая модель центрального поля
п
гравитации Земли с центральным ядром гравитации массой М , равной массе Земли, по точности результатов не обеспечивает экспериментально определенной разности Яе - Яр = 21382 [м] радиусов осесимметричной фигуры Земли на экваторе и на ее полюсах. Ошибка составляет примерно 50%.
Модифицированная расчетная схема расположения основных гравитирующих масс Земли в виде центрального сферического ядра и однородного экваториального пояса (центральное ядро массой МЯ » 0.89 М и экваториальное кольцо радиусом ЯК » 860 [км], массой МК » 0.11 М) позволила получить существенно более точные результаты уже не только в отношении (Яе - Яр) - величины разности характерных размеров фигуры Земли, но также и в отношении параметров нецентрального гравитационного поля Земли, в частности величины и направления локального вектора ускорения силы тяжести на ее поверхности. Полученная расчетная форма Земли внешним образом охватывает эллипсоидальную форму Земли. Причем обе поверхности контактируют между собой лишь на экваторе и на полюсах, а в точках параллели с широтой ф=35 имеют в радиальном направлении максимальное расстояние между собой, равное 90.9 [м].
Кроме геометрических методов воспроизведения фигуры Земли, основанных на аппроксимации экспериментальных данных, ее математические модели, базирующиеся на физических законах гравитационного взаимодействия, также позволяют уточнять внутреннее строение Земли и оценивать положение и размеры локальных источников гравитации.
Литература
1. Пантелеев В.Л. Курс лекций «Теория фигуры Земли» // М.: МГУ им. М.В.Ломоносова, 2000, 98 с.
2. Молоденский М.С. Гравитационное поле, фигура и внутреннее строение Земли // М.: Наука, 2001, 569 с.
3. Кузнецов В.В. Физика Земли // Новосиб.: 2011, 840 с.
4. Параметры общего земного эллипсоида и гравитационного поля Земли // Москва: РИО ТС ВС РФ, 1991, 68 с.
5. Русанов П.Г. Алгоритмические концепции метода твердых тел // Известия МГТУ «МАМИ». Научный рецензируемый журнал. Серия «Естественные науки». - М., МГТУ «МАМИ» № 3(17), 2013, т. 1. С. 124-136.