Все научные статьи

Ферронский В.И. Вириальный подход к решению задачи о глобальной динамике Земли

Научная статья

 

Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 12и'а Ферронский В.И. (ferron@aqua.laser.ru)

Институт водных проблем Российской академии наук

С времен Ньютона и Клеро считают, что Земля является инертным телом, ее вращательное движение происходит под действием сил инерции, а динамические эффекты определяет лунно-солнечный потенциал сил [5]. Идея инерционного вращения Земли появилась при рассмотрении задачи о фигуре планеты как о вращающемся сферическом теле, заполненном жидкостью и находящемся в гидростатическом равновесии в однородном силовом поле Солнца. Эта идея основывалась на представлении о том, что сумма внутренних сил взаимодействующих масс и вращающих моментов планеты равна нулю. Такое представление пришло из задачи двух тел, где последние в рамках оговоренных допущений были приняты за точечные массы, из которых, по Ньютону, исходили силы их притяжения, образуя центральное поле. В результате, в геодинамике утвердилось и до сих пор остается умозрительное представление о гидростатическом равновесии масс планеты и ее инерционном вращении. Что касается наблюдающихся эффектов прецессии и нутации оси вращения при движении Земли, то их объясняют возмущением Луны и Солнца, связанным с возможным избытком массы в зоне экватора из-за эллиптичности планеты.

Современные наблюдательные факты свидетельствуют о том, что идея гидростатического равновесия Земли не подтверждается. По результатам анализа большого числа измерений зональных и тессеральных гравитационных моментов, выполненных в последние десятилетия с помощью геодезических искусственных спутников при исследовании гравитационного поля Земли, установлено, что планета не находится в состоянии гидростатического равновесия. Ее фигура отклоняется от нормального эллипсоида вращения на величину квадрата сжатия, т.е. Ч(1/300) [1, 3, 6]. Как отмечает Мельхиор [6], теперь есть доказательство того, что Земля не находится в гидростатическом равновесии и это вызывает трудности в интерпретации распределения плотности на основе сейсмических измерений. Однако, добавляет он, мы вынуждены использовать это условия, поскольку нет ничего лучшего. Следует отметить, что в модели инерционного вращения Земли значение ее потенциальной энергии оказывается на три порядка выше кинетической энергии, что противоречит


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 11та Проблема вращения Земли обсуждалась на недавнем семинаре НАТО [14]. Было отмечено, что две задачи этой проблемы остаются неразрешенными. Ими являются вариации продолжительности суток и наблюдаемое Чандлеровское движение полюсов с периодом 14 месяцев против 10 месяцев, которые дает модель твердого тела Эйлера. Отмеченные факты, а также наблюдаемое изменение гравитационного поля, неравномерности угловой скорости вращения Земли, тектоника литосферных блоков и плит и глубинные геотектонические процессы свидетельствуют о необходимости развития новых физических подходов в решении задач динамики планеты.

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

Приведение сил тяжести и инерции к равнодействующему сфероиду (эллипсоиду) силового давления

Будем исследовать задачу о динамике Земли как самогравитирующего одномерного шара с однородным и неоднородным распределением плотности массы, непрерывно распределенной по его объему. Движение шара будет относительным и происходит в собственном силовом поле и в силовом поле Солнца.

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


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 12иуа Для рассмотрения движения одномерного шара в собственном силовом поле его поступательное (орбитальное) движение относительно фиксированной точки (Солнца) нужно отделить от двух других составляющих движения. После чего можно рассматривать как вращение тела относительно центра масс под действием собственного поля сил, так и движение, связанное с изменением структуры и формы. Такого отделения требует лишь момент инерции шара, который зависит от выбора системы координат. Силовая же функция, определяемая как эффект взаимодействия всех пар частиц массы шара, не зависит от ее выбора [2]. Момент инерции шара относительно солнечной системы отсчета необходимо разложить на момент инерции его центра масс относительно той же системы отсчета и на момент инерции планеты, взятый в собственной системе отсчета. Чтобы сохранить условия инерциальности собственной системы отсчета, совместим ее с геометрическим центром масс.

Итак, примем абсолютную декартову систему координат Ос^т]^ с началом в геометрическом центре Солнца и перенесем ее параллельно осям в геометрический центр симметрии масс шара, обозначив эту систему через Oxyz(рис.1). Момент инерции шара как инертного тела относительно солнечной системы отсчета будет

где nij- инертная частица массы шара в солнечной системе отсчета; Rj- ее расстояние от начала солнечной системы координат.

Для разделения момента инерции (1) воспользуемся методом Лагранжа, который основан на его же алгебраическом тождестве вида

М<г'<иа 'аа М<г<иа 'а М<г<иа 'а ^ \<i<n\<j<n

где atи bj- какие угодно величины; п - юбое целое положительное число.

Математическое преобразование, связанное с разделением момента инерции п взаимодействующих материальных частиц относительно любой системы координат на две алгебраические суммы впервые было выполнено Якоби в его "Лекциях по динамике" [2, 10, 15]. Якоби было показано, что если ввести обозначения (рис. 1)

гi = Xi + A; л. = У + В;аа Ci = z + C;

2>=M; ^m^=MA;аа ?вд = MB;а щ? = МС,(3)

где А, В, С - координаты центра масс в абсолютной системе отсчета,


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 121иа

то на основе тождества (2) получим

+Z ^ + 2В - ШгУг + В' Z Щ1 +S Ш? + 2С - Щ^ + ^ S ^ '

Поскольку

М А = ?>?. = ?>*,. + YjnA= 1>Л + Ш > то

2] дагхг = 0,аа а также ? "У.- = >а S "^ = Х

Теперь момент инерции (1) принимает форму

2>Д2 = м(А2 + В2 + C2) + ^ml(xf +yf+zf),(4)

где

М(А2 + В2 + С2) = MR2m,(5)

2]Щ ( х2 + у2 + z2) =Мгт2 ,(6)

М - масса шара; Rmи гт - радиусы инерции шара, описывающие сферические поверхности с сосредоточенной и распределенной массой, которые определяют момент ее инерции в солнечной и собственной инерциальной системе отсчета.

Таким образом, момент инерции вращающейся вокруг Солнца массы шара в инерциальной системе отсчета мы разложили на два алгебраических слагаемых. Первое (5) представляет момент инерции шара в солнечной системе отсчета Ос^т]^. Второе слагаемое (6) представляет момент инертной массы шара в собственной системе отсчета Oxyz. Этот момент инерции может определяться в любой системе отсчета с началом в центре масс О. Учитывая симметрию одномерного шара, примем полярную систему отсчета с началом в центре О. Тогда выражение (6) для полярного момента инерции шара примет вид

!р = Z тг ( Х>2 + У'2 + Z>2 )=Цтг Г^ = М ГЩ'

Откуда радиус инерции гт , описывающий сферическую поверхность, будет

r*=J^7J--(8)

М

Здесь М = У^ да. - масса шара относительно собственной системы отсчета. При сферической симметрии шара выражение (8) можно записать в виде

Га 2'а '

MJ0' v 'а MR2

Ч f г24г 2p(r)dr = ЧЧ\rAp(r)dr,аа (9)

или


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 1111а

4^ = ^_____________ = Р2МЯЪ

AttR2MR2MR2

4ж\г4 p{r)dr

Р(Ю)


откуда

rm2=jB2R2,

где р(г)- закон радиального распределения плотности массы шара; R- его радиус; р -безразмерный численный коэффициент, представляющий отношение поверхностей оболочек с приведенным радиусом сфероида инерции гт и шара радиусом R.

Значение р в зависимости от закона распределения плотности р(г) изменяется в пределах 1 > р >0. Этот коэффициент мы ранее назвали структурным форм-фактором момента инерции [15].

Аналогичным образом получим выражение для приведенного радиуса сфероида сил тяжести rgкак отношение момента сил гравитационного взаимодействия масс оболочек шара с радиальной плотностью р(г), к моменту сил взаимодействия массы шара, распределенной по его внешней оболочке с радиусом i?, т.е.


47ir2g

-------- \rp{r)m{r)dr

Го

GM2 R2

2GM2 а2ЧЧ R2

- п'1

4ж R2

GM2 R2


(11)


Выражение для радиуса тяготения, записанное через силовую функцию, будет


R2

.2

AnG\rp(r)m(r)drа а^ GM

4^ =_________ о____________ =______ R_ = a2(и)

4л-i?2а GM2GM2


RR

r

где в соотношениях (11) и (12) m(r) = 4ж\ г2p{r)dr.

о Безразмерный коэффициент = rgIRесть отношение площадей поверхности

сфероида сил тяжести с приведенным радиусом rgи радиусом шара R. Его значение

зависит от закона распределения плотности р(г) и изменяется в пределахаа 1 > аа >0.

Ранее [15] коэффициент а был назван структурным форм-фактором силовой функции.

Численныеа значенияа безразмерныхаа структурныхаа коэффициентоваа аираа для

некоторых законов радиального распределенияа плотности р(г),аа полученныеа путем

интегрирования выражений (10) и (12) для полярного момента инерции и силовой


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 121iа Таблицаа 1. Численные значенияа форм-факторов аа иа |3аа для радиального распределения плотности массы и политропных моделей шара


Закон распределения Индекс политропы


Радиальное распределение плотности массы

0.6аа 0,4

0.74 0.27

0.71аа 0.29

0.16k

Р(г) = Ро

p(r) = pc(l-r/R)

р(г) = рс(1 - r2/R2),

p(r) =рс ехр(1 - kr/R) р(г) = рсехр(1 - kr2/R2) р(г) = рс 5(1 - r/R)

0 1

1,5

2 3 3,5


8/kz

1/k

0.5аа 0.67

0.6

0.75

0.87

1.0

1.5

2.0

Политропные модели

0.4

0.26

0.20

0.15

0.08

0.045


0.6

0.4

0.42

12/к2

1.5/к

1.0

0.6

0.38

0.30

0.23

0.12

0.07


Из таблицы видно, что для однородного шара, где p(r)=const., приведенные радиусы

инерции и тяготения совпадают. Их безразмерные структурные коэффициенты аир

численно равны 3/5,а вращательныеа моменты сил тяжести и инерции

уравновешиваются и поэтому вращение масс отсутствует. Так что

г'а г~а 3

.2а .2а ^

RlRl

(13)

откуда

rm = rg = л/з/5Я2 = 0,7745966R.а (14)

Для неоднородного по плотности шара при p(r) ^ constиз (10)-(12) имеем


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа lIli 0<-^-<-<--------- *-<\,(15)

4Я2а 4Я2

Из неравенства (15) и табл.1 видно, что у неоднородного шара по сравнению с

однородным при возрастании плотности к центру радиус инерции уменьшается, а

радиус тяготения увеличивается. Так как rm ^rguгт< 0,77R< rg, то между объемными

силамиа взаимодействияаа масса оболочека шараа иа силамиа иха инерцииа появляется

неуравновешенный вращающий момент. Теперь из соотношения (15) следует, что

rm=rmo-Srmtаа иа rg = rgo+8rgtаа ,(16)

где индексы 0 и tотносятся к однородному и неоднородному шару.

Согласноаа (15)аа иаа (16)аа вращениеаа оболочекаа одномерногоа шар будетаа не

твердотельным, а оболочечным и асинхронным.аа При возрастании плотности масс к

поверхности знаки в выражениях (15) и (16) изменятся на обратные. Это замечание

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

определяет прямое и обратное направление вращение тела.

Основной вывод, который следует из приведенного выше рассмотрения состоит в

том, что поле сил самогравитирующего тела приводится не к равнодействующей силе,

проходящейа череза геометрическийа центра симметрииа масс,а аа ка давлениюа сила по

замкнутой поверхности сфероида (эллипсоида). В случае однородного тела

приведенныеаа радиусыаа гравитацииаа иаа инерцииаа совпадают,аа моментыаа вращения

гравитационныхаа иаа инерционныхаа силаа уравновешиваются. Уаа неоднородногоаа по

плотности тела приведенные радиусы инерции и гравитации не совпадает. Между

силами взаимодействия и инерции образуются неуравновешенные момент вращения и

возникает сжатие оболочек. В приводимом ниже аналитическом рассмотрении задачи

мы определим некоторые параметры той и другой системы.

Динамическое равновесие движения и уравнения вращения и колебания

В небесной механике и геофизике для выражения силовой функции и момента инерции неоднородного шара и эллипсоида вращения обычно прибегают к разложения этих величин в ряд по сферическим и эллиптическим функциям. Анализ результатов измерения гравитационных моментов Земли с помощью искусственных спутников показал, что все четные моменты в разложении потенциала по сферическим функциям, кроме второго, содержат член, равный квадрату сжатия планеты [3]. Соотношения (15)-(16) выражают ту же идею. Больше того, они говорят о том, что для решения задачи о динамикеаа неоднородногоаа самогравитирующегоаа шар егоаа силовуюаа функциюаа и


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 1214а 4>(S)=fk^U&,

оа Ро где s= r/Rаа - отношение текущего к полному радиусу; р0 - средняя плотность шара

радиусом г ; рг - радиальная плотность неоднородного шара;а х - текущая координата;

R

значение (рг- р0) удовлетворяет условию I (рг - p0)r2dr= 0, а функцияаа Щ1)=0.


Видно, что функция *F( s) выражает радиальное изменение плотности массы неоднородного шара относительно ее среднего значения на расстоянии r/R. После замены переменныха с помощью этой функции выраженияа потенциальнойа Uи

жу

(17) (18)

(19)

(20)

кинетической энергии K=Jpa>аа неоднородного самогравитирующего шара разлагаются в форме [16, 17]

GM2

GM2

+ 3 Г щек + Ч\\ Ч dx

U=a

3 аа ,а 9 р( у/

R

2-чх

Цщск

2 Д2

MRla

RK=p2MR2co2

или после соответствующего приведения

GM2

K=(j302-2j3t2)MR2a)2,

R

гдеа а о = Роаа и 2at= Д , а индексы о, t, у означают радиальную, тангенциальную и диссипативную компоненты рассматриваемых величин.

Поскольку потенциальная и кинетическая энергии однородного шара равны между собой (<%2=/?о2=3/5), то

Uo = K0,(21)

Eo = Uo + K0 = 2U0.(22)

Для выражения динамического равновесия между потенциальной и кинетической энергией взаимодействия неоднородностей с однородной массой из (17)-(18) имеем


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 1215а Et=Ut + Kt=3Ut,(24)

где Ео, Et ,Uo, Ко, Ut, К - полная потенциальная и кинетическая энергия колебания и вращения соответственно.

Уравнения (21)-(24) представляют выражения усредненной теоремы вир нала для самогравитирующей однородной и неоднородной систем, которые определяют условия их динамического равновесия [4]. Потенциальная энергия Urвзаимодействия самих неоднородностей теряется с граничной поверхности тела в виде излучения и является механизмом эволюционных процессов, которые непрерывно происходят.

Из эффектов взаимодействия однородной и неоднородной по плотности массы находим, что, как и следует из классической механики для диссипативных систем, вращательный момент сил N неоднородной гравитирующей системы относительно ее центра масс не равен нулю, угловой момент Lсистемы не является постоянной во времени величиной, а энергия непрерывно расходуется при движении системы во внешнем поле на преодоление сопротивления трения и на поддержание равновесия в пространстве, т.е.

N = Ч Ф0,аа L Ф const.ЕФ const. >0.

dt

Система физически не может быть консервативной, если в ней присутствует трение или иные диссипативные процессы, поскольку из-за них величина Fdsбудет всегда положительной, а интеграл не может исчезнуть, т.е. [19]

зF-ds>0.

Теперь, после того как найдено, что результирующая внутреннего гравитационного поля не равна нулю и что динамическое равновесие системы определяется вириальным соотношением между потенциальной и кинетической энергией, могут быть записаны уравнения движения самогравитирующего тела.

Ранее для описания и исследования движения однородного и неоднородного самогравитирующего шара нами использовалось вириальное уравнение Якоби [10, 15, 18]. Якоби (1884) вывел его из уравнений движения Ньютона для системы п взаимодействующих точечных масс и свел задачу многих тел к ее частному случаю - к задаче одного тела с двумя независимыми переменными вида [2,15]

Ф = 2Е-и,(25)

где Ф=1/21 - функция Якоби; / - полярный момент инерции; Е = U+ К - полная энергия; UиК - потенциальная и кинетическая энергия системы.


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 111ьа

теоремы о вириале сил Клаузиуса. Как известно, Клаузиус выводил эту теорему в приложении к термодинамической задаче кинетической теории газов и, в частности, применительно к условиям появившихся в то время машин Карно, которые реально работали во внешнем для этих машин поле сил тяжести Земли. В этой связи, для кинетической энергии системы газ - камера сгорания машины в выражение для кинетической энергии был введен коэффициент 1/2. Так что при обозначении суммарной кинетической энергии, которую тогда называли "живой силой", появился коэффициент 1/2, т.е.

K = -Ymv2.

Как отмечал сам Якоби [10], смысл этого коэффициента состоит в том. что у машины Карно учитывалась только та работа кинетической энергии, которая оплачивалась, т.е. работа, совершаемая машиной, а не силой тяжести Земли. Например, при работе парового молота для забивки свай в машине учитывалась и оплачивалась только кинетическая энергия, необходимая для подъема молота, а энергия его падение происходила за счет силы тяжести в поле Земли и она не оплачивалась. Выше в (21)-(22) было показано, что в случае движения тела в собственном силовом поле коэффициент 1/2 исчезает, поскольку в этом случае тело движется только за счет собственной энергии.

Нами было найдено приближенное решение уравнения (25) для неоднородных гравитирующих систем с высокой симметрией распределения плотности массы, для

которых С/л/ф = const. [15]. Теперь, после разложения силовой функции и момента инерции тела, при Ur= 0 мы получим строгое решение уравнения движения (25), записав с учетом (22) и (24) два уравнения для радиальной и тангенциальной составляющей в виде


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 121'а

(26)

<Ьй=-Ей-иД

(27)

<bt=\Et-Ut.

Учитывая найденную строгую связь (21) и (23) между потенциальной энергией и моментом инерции через структурные коэффициенты а о = Ро и 2at = Pt, уравнения (26) и (27) приводятся к уравнению с одной неизвестной вида

В

(28)

Ф = -А +

где А и В - постоянные величины.

л/ф = Ч Г1 - scos(^- ф ) ] А1

Общим решением уравнения (28) в полярной системе отсчета будет [17] В<

(29)

АВ

ж[г-esin(г-q>)],

{2АГ- "................ ^ Ч (30)

где sи сра - постоянные интегрирования, зависящие от начальных значений функции

Якоби Ф и ее первой производнойа Фа в момент времени to;аа <^

вспомогательная независимая переменная; А = Ао = -1/2 Ео >0,а В = Во = UoJ(г> для

радиальных колебаний; и А = At=а 1/3 Et>0 , В =Bt =Utд/ф"а для вращения.

Выражения для функции Якоби и ее первой производной в явном виде получим после соответствующих преобразований в форме рядов Лагранжа [15]


Ф


В2


\ + -s2 + 2


2s+ -


,зЛ

J


cos L------ cos 2L----- cos 3Z +.

2а 4


S2(2а 2аа \

Ф

' sB

sinZ + Ч гsin2ZH------ sinZI2cos Z-sinа Zl-

2 2 vа ;

Частота колебания a>orи угловая скорость cotrоболочки на радиусе г будут [15]


(Ос


4Вп


ип


\<GMr

1 PIS


j-xGpor


(31)


2f/.

(24)

3/2

(32)

О р0гке

COtr

У.

1РУе

гдеаа Uorиаа Utrаа -аа радиальнаяаа иаа тангенциальнаяаа компонентыаа силовойаа функции (потенциальной энергии); Jorи Jtr= 2/3Jor- полярный и осевой момент инерции;


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 121ла V

значение плотности шара радиусом г; Vr- объем шара радиусом г;аа 2atr=ptr;а кег -безразмерный коэффициент, учитывающий эффект динамического сжатия оболочек.

Выражения (29)-(32) представляют законы движения Кеплера. При однородном распределении плотности массы частота колебания всех оболочек шара будет единой.

Вращение каждой оболочки зависит от ее плотности и приливного трения со стороны внутренней массы, определяемого коэффициентом кег . Его значение для внешней оболочки определится из выражений (31) - (32) и будет равно отношению частоты радиального колебания к угловой скорости, т.е.

2а 2

_ cotat -Сро

По наблюдательным данным не трудно найти, что значения ке для тел Солнечной системы соответствуют величинам динамического сжатия этих тел.

Найдено, что безразмерный коэффициент ке е[0,1] в случае трехосного эллипсоида с осями а, Ъ, с для эллипсоидального закона распределения плотности равен [15]

к= FjcpJ)1а2+Ь2+с2

sirup/ За2

жа а2-г fа К"*2аа Wfi

где ф = arcsinJ------ Ч , 1= , Ч----- -,аа a t(ф,1) - неполный эллиптический интеграл

V аV а - с

первого рода в нормальной форме Лежандра.

Колебание и вращение Земли

Итак, в дополнение к уже полученному ранее решению о радиальных колебаниях Земли [15], теперь мы имеем строгое решение о ее вращении. Из формулы (31) видно, что радиальное колебание оболочек тела не зависит от фазового состояния массы и определяется ее средней плотностью. Корректность выражения (31) подтверждается результатами наблюдений. Так, период радиальных колебаний внешней оболочки Земли по формуле (31) и по нашим измерениями равен ~1.4часа [8, 9, 15], а Солнца по формуле (31) и по измерениям - -2.8 часа [7].

Период вращения внешней оболочки Земли по формуле (32) равен ~ 24 часам. Здесь коэффициент геодинамического сжатия кег = 1/289.37 принят по данным измерения гравитационных моментов с ИЗС [6].


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 121уа

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

IX. _U0 _а] _ 0.6

Uа2

Cosй-

(33)

0.66-сс,

IX

0.918,аа 0=23.5

гдеаа ара -аа поправк н солнечно-лунноеаа возмущениеаа ваа третьемаа знаке,аа которую принимаем раной 0.006.

Дляа выявленияа природыа прецессииа иа нутацииа оси Землиа необходимоа иметь сведения об условиях формирования ее оболочек, которые рассмотрим ниже.

Дифференциация массы Земли по плотности

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

Как известно из теоремы Ньютона о гравитационном взаимодействии материальной точки и шарового слоя, последний в силу условий симметрии не оказывает никакого силового воздействия на точку, находящуюся внутри этого слоя. Напротив, материальная точка, находящаяся за пределами шарового слоя, подвергается с его стороны силовому воздействию. На этой теореме основана приливная динамика известного французского астронома Роша. Его подход состоит в следующем [16].

Два тела с массами М и да взаимодействуют по закону притяжения Ньютона (рис.2а). Пусть Мт и Rr, где г - радиус тела т, a R- расстояние между телами Ми т. Полагая, что масса тела М равномерно распределена в пределах сферы радиусом R, запишем выражения для ускорений точек А и Ва тел т в виде


GMОт


GMGm


1Аа (R-rfг2 'аа "а (R+rfг2 Х

Относительное приливное ускорение точек А и В будет выражаться


М

М

2/V

1АВ


G

2т

2

(R-r)аа (R + r)аа r

4ж Д ^аа 4Rr

G

pMR'

(i?2-r2)2


8л-


Gr(2pM-pm).


(34)


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 122иа

з4аа 3

Здесьаа рм=М/ЧЕаа иаа рт=т1Чжга выражают средние значенияаа плотности

распределения массы сфер радиусом Rи г. Критерий Роша говорит, что тело массой т устойчиво против приливной силы тела М, если средняя плотность тела т вдвое превосходит среднюю плотность сферы радиусом R.

Теперь с помощью динамикиаа Рош аоценим приливную устойчивость сферического слоя радиусом Rи толщинойа г = RB- RA(рис.2b). На слой массой т и

средней плотностью рт = т / 4nRа г в точке А действует приливная сила от сферы

4 радиусомаа RA. Масса сферы равн М, а средняя плотностьа рм = М / ЧnR3Aаа .

Приливную силу в точке В генерирует сфера радиусом R+ г и массой М+ т. Тогда ускорения в точках А и В будут

GMG(M + m)

ra(RA +r)

Относительное приливное ускорение в точках А и Ва запишется как

GM

Gm

(RA+r)2 =

К

&+')

(35)

nGpM- 4Срт )r= 4Сг[ -рм - pm),(R г)

3-----лаа -----;аа ж--аа уаа P"~oJ.(Rr)

Выражения (34) и (35) позволяют понять природу геодинамических эффектов реальной Земли. Из приведенного выше рассмотрения задачи о приливном ускорении внешнего неоднородного сферического слоя шара следует, что при рм ^ рт соотношение (35) определяет эффект гравитационной дифференциации масс шара по плотности. В частности, при рм<Рт оболочка погружается (притягивается) до уровня, где рм=Рт ж Если рм >рт , то оболочка всплывает (отталкивается) до уровня рм=Рт-, При рм >2/Зрт оболочка становится самогравитирующей. Отсюда находим, что в случае возрастания плотности от поверхности к центру, как это наблюдается у Земли, каждый вышележащий слой находится во взвешенном состоянии и "плывет" на подстилающем слое за счет "выталкивания" силами Архимеда, которые физически представляют радиальную компоненту сил гравитационного взаимодействия масс.

Эффект гравитационной дифференциации масс объясняет природу формирования земной коры и океанов, геотектонических, орогенных и сейсмических процессов, включая землетрясения. Все эти процессы являются следствием непрерывно продолжающегося процесса разделения планетной массы по плотности.а Очевидно,


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 12ilа полярный радиус инерции равен rm=3/2 r"Lm=Vl.5- 0.3315R2 = 0,70516R=4.493'106м,аа а

радиус тяготения будет rg = ]а2К2 = 0,8164R= 5,20110 м. В этой связи известная интерпретация радиального распределения плотности по данным сейсмического зондирования, которую дал Буллен [6], требует пересмотра. Оценим возможный характер усредненного непрерывного радиального распределения плотности, опираясь на спутниковые и сейсмические данные.

Радиальное распределение плотности Земли

Как известно, основой современного представления о радиальном распределении плотности являются опытные данные о скорости распространения продольных и поперечных сейсмических волн. В интерпретации Буллена [3, 6] сейсмические данные дают следующую картину. Плотность земной коры равна 2.7-2.8 г/см и возрастает по некоторой кривой к центру планеты до -13.0 г/см со скачками значений на границе Мохоровичича, между верхней и нижней мантией и на границах внешнего и внутреннего ядра. Буллен ввел скачки плотности после неудачной, как он решил, аппроксимации сейсмических данных параболической кривой, при которой получалось, что плотность ядра к центру должна падать. Буллен, естественно, не предполагал, что радиус инерции и радиус тяготения тела не находится в его центре масс, а поэтому там не находится и максимум плотности. Теперь, когда наши представления об этих величинах изменились, задача о радиальном распределении плотности требует нового рассмотрения.


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа lIZI Rrг г2а г4 3 3

М=4ж\г2р(г)Ф = 4ж\г2р0(-а Ч + ЪЧ + c)dr= Чжр0Я3(Чал- Ч Ь + с),

оаа о

3 3 где члена Ч алЧ Ь + с = 1аа позволяет вычислять и строить кривые в безразмерной

форме. На рис.3 показан спектр этих кривых от линейной зависимости 1 с максимумом в центре шара, до линейной зависимости 7 с максимумом на его поверхности. Кривые 1-7 пересекают оболочку 10 со средней плотностью и радиусом инерции rm=0.775R, который здесь совпадает с радиусом тяготения rg. Спектр кривых распределения плотности дает принципиальную картину ее перераспределения в результате гравитационной дифференциации за историю Земли. Находим, что начальная кривая распределения плотности находилась в районе кривых 7-8 со значением плотности ~7-8 г/см . Кривая плотности современного этапа эволюции находится около кривых 5-6 с плотностью ~ 2-3 г/см (или меньше) на поверхности, ~1-4 г/см в геометрическом центре и 7-8 г/см между оболочками rmи rg. Эти величины кривой плотности соответствуют найденным по орбитам ИСЗ значениям структурных форм-факторов |3 = 0.49725 и p2j_ = 0.3315 и форм-фактора а=0.6601 силовой функции.

Распределение гравитационного потенциала и силы тяжести Земли

На основе принятого параболического закона радиального распределения плотности массы для приведенных на рис.3 кривых были рассчитаны и кривые радиального распределения гравитационного потенциала и силы тяжести (рис.4а, 4b) для пробной массы т=], взаимодействующей с массой шараМ. Расчеты выполнены по известным уравнениям теории притяжения [15]

TTf^ 4О} 2, . ,а _Л , ..аа GM\ 3 г4а 1 г3а 1аа г2а 3 , 3а п

и (г) =------ \г. руглаг +АЬ\г.р{г.)аг. =------------ аЧ----- оЧ----- сЧг---- а + о + Чс\,

гоа Jа R


riR


3 г3 Ъ г2 г ЧаЧг + ЧЬЧг + сЧ 5аа R34аа R2R


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа lIIi Как видно из рисунков 4а и 4b, при любом законе распределения плотности массы гравитационный потенциал в центре тела имеет положительный знак и его численное значение выше, чем на поверхности тела. Это означает, что в центре тела имеет место не гидростатическое давление вышележащих масс, как это считается в механике однородных консервативных систем, а давление излучения.

Для дальнейшего анализа возможного закона распределения плотности нами выбраны четыре из 10 приведенных выше парабол с уточненными численными коэффициентами (рис.5). Они удовлетворяют условию равенства относительного радиуса осевого момента инерции значению Р _i_ = 0,3315, найденному по данным анализа орбит ИСЗ, и соответствуют значению форм-фактора полярного момента инерции Р =3/2|3 j_ = 0,49725. Численные значения параметров этих кривых представлены в табл. 2.

Таблица 2. Физические иаа динамическиеаа характеристикиаа Землиаа дляаа кривых радиального распределения плотности, представленных на рис.П.5.

4 кривой

1

2

3

4

ps, г/см3

2.76

2.08

1.65

1.03224

Рс г/см

13.8

10.455

6.315

1.6284

Ртах, Г/СМ3/ КМ

13.8/0

10.455/0

8.26/2096

8.57/3122

р\

0.33(3)

0.3315

0.3315

0.3315238

р2

0.50

0.49725

0.49725

0.49725858

p2t

0.10

0.10275

0.102752

0.10 2714

а

0.6607142

0.6607374

0.6607374

0.660143

at

0.05

0.05

0.0513714

0.0513571

ау

0.0107142

0.009366

0.009366

0.0087859

гё,км

5178.6

5178.7

5178.6

5176.4

гт, км

4504.9

4492.6

4492.6

4492.7

Обозначения: ps, рс, ртах-плотность поверхностная, в центре планеты и максимальная ; Р _ц Р , Р t -форм-факторы осевой, полярной и тангенциальной компонент радиуса инерции; a, at , ау - форм-факторы радиальной, тангенциальной и диссипативной составляющих силовой функции; rg , rm - радиусы тяготения и инерции.


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 1224а Прецессия и нутация оси вращения

По сейсмическим данным скачки плотности планеты наблюдаются на границах между литосферой и верхней мантией (-350-400 км), верхней и нижней мантией (-1000 км), нижней мантией и внешним ядром (-2700-2900 км) и между внешним и внутренним ядром (-5400 км). Из условия, определяемого выражением (32), эти границы можно рассматривать как поверхности изменения угловых скоростей вращения оболочек Земли и их динамического сжатия. Согласно выражению (33), интегральный эффект вращения всех оболочек и изменение формы приведенных эллипсоидов сил инерции и тяжести демонстрирует прецессия земной оси. Но вклад в этот эффект каждой из оболочек будет разный. Наблюдаемое суточное вращение Земли очевидно относится в к ее верхней оболочке, ограниченной поверхностью Мохоровичича (350 км), где обнаружен первый разрыв плотности. Верхняя и нижняя мантия составляют большую часть массы планеты. Они вносят основной вклад в этот эффект. Учитывая известные из наблюдений значения периодов прецессии и климатических ледниковых и межледниковых эпох, можно предположить, что средний период вращения верхней мантии может быть примерно вдвое меньше (-13000 тыс. лет), чем период нижней (-26000 тыс. лет). Внешнее ядро должно иметь весьма малый коэффициент динамического сжатия и его период вращения велик, а угловая скорость мала. Внутреннее ядро с нулевой скоростью поперечных волн не имеет вращения и скорее всего представляет жидкое или газообразное образование с низкой плотностью и низким давлением (-1.5 бар).

На внешней оболочке Земли "плавают" во взвешенном состоянии и имеют иную собственную угловую скорость твердая кора и океаны. Согласно уравнениям (29)-(32) вращение оболочек происходит по третьему закону Кеплера. Это означает, что при


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа lIIi Инерционное качания земной коры и океанов с их реверсивным ускорением и торможением вращательного движения верхней оболочки известно в геофизике как эффект Кориолиса. Кроме нутации оси он объясняет природу таких важных геодинамических процессов, как тектоника плит, океанические течения, неравномерность вращения Земли и продолжительности суток, короткопериодные изменения погоды и климата и синхронность качания с лунными приливами.

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

Такова природа прецессии и нутации оси Земли, вытекающая из рассмотрения ее динамики как самогравитирующего тела.

Заключение

Мы получили принципиально новое решение задачи о вращении и колебании небесного тела в собственном силовом поле, которое генерируется самим телом, представляющем как консервативную, так и диссипативную систему. Из рассмотрения этой задачи следует, что двухкомпонентная модель плотности массы Мак-Миллана [20] для объяснения природы потенциальной энергии взаимодействия и подход Коимми


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 111ьа Решение задачи о динамике самогравитирующего тела может найти множество новых приложений в астрономии, науках о Земле и о природных изменениях окружающей среды обитания.

итература

  1. Грушинский Н.П. Теория фигуры Земли. М.: Наука, 1976. 512 с.
  2. Дубошин Г.Н. Небесная механика. Основные задачи и методы. М.: Наука, 1975. 800 с.
  3. Жарков В.Н. Внутреннее строение Земли и планет. М.: Наука, 1978. 192 с.
  4. Киттель Ч., Найт У., Рудерман М. Курс физики, Т.1. (Перев. с англ). М.: Наука, 1971. 480 с.
  5. Клеро А.К. Теория фигуры Земли, основанная на началах гидростатики. (Перевод с франц.). М.-Л.: Изд. АН СССР, 1947. 288 с.
  6. Мельхиор П. Физика и динамика планет. (Перевод с франц.). М: Мир. 1976. 411 с.
  7. Северный А.Б. Некоторые проблемы физики Солнца. М.: Наука, 1988. 222 с.
  8. Ферронский СВ. // Физика атмосферы и океана. 1984. Т. 20. С. 922.-928.

9.а Ферронский В.И., Денисик С.А.,Ферронский СВ. // Физика атмосферы и океана.

1984. Т. 20. С. 802-809

  1. Якоби К. Лекции по динамике. (Перевод с нем.). М.-Л.: Техлитиздат, 1936. 252 с.
  2. Caimmi R. //Astron. Nadir. 1992. V. 313. Р.165-182.
  3. Caimmi R. Private communication. 1997.
  4. Caimmi R, Secco L. //Astron., Astrophys. 1990. V. 237. P. 336-344.

14.а Cazenave A. (Ed.). Earth Rotation: Solved and Unsolved Problems. Proc. NATO Advanced

Research Workshop. Dordrecht: Reidel, 1986. 320 p.

  1. Ferronsky V.I.,Denisik SA.,Ferronsky S.V. Jacobi Dynamics. Dordrecht: Reidel, 1987.366 p.
  2. Ferronsky V.I., Denisik S.A., Ferronsky S.V. // Celest. Mech. & Dynam. Astron. 1996. V.64. P.

167-183.

  1. Garcia Lambas D., Mosconi M.B., Sersic J.L. // Astroph. & Space Sci. 1985.V. 113. P. 89-98.
  2. Giordano СМ., Plastino A.R // Celest. Mech. & Dynam. Astron. 1999. V. 75. P.165-183.
  3. Goldstein H.>
  4. Mac Millan W.D. The theory of the Potential. New York: Dover, 1930. 260 P.

21.а Whitteker E. Т., A Treatise on the Analytical Dynamics ofаа Particles and Rigid Bodies.

Cambridge: Cambridge Univ. Press, 1937. 420 P.


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 12i'а Рис. 1. Отделение момента инерции Земли в собственной инерциальной системе отсчета при ее относительном движении в силовом поле Солнца.

Рис. 2. Приливная гравитационная устойчивость тела (а) и шарового слоя (Ь) по Рошу.

Рис.3. Спектр кривых радиального распределения плотности массы Земли по параболическому закону в безразмерной форме.

Рис. 4. Спектр кривых радиального распределения силовой функции (а) и силы тяжести (Ь) Земли для пробной массы при параболическом законе распределения плотности массы в безразмерной форм.

Рис.5. Возможный спектр непрерывного распределения плотности массы Земли при значении структурного коэффициента J3_i_ = 0.3315, найденного на основе анализа траекторий ИСЗ.

(а)


Электронный научный журнал ИССЛЕДОВАНО В РОССИИа 122ла

AV ^mаа 1В

ftlli


М


R



1.0 r/R

0.5 Рис.3


1.0аа r/R


(b)

Ur/Uo


Рис.2


qr/qo


1.0 r/R


1.0аа r/R



Рис. 4b


     Все научные статьи