WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

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

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

 

Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     12и'     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

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

Ферронский В.И. (ferron@aqua.laser.ru)

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

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

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     11т     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

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

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

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

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

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

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     12иу     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

изменением его формы и структуры [2]. В задаче двух тел об орбитальном движении тела последними двумя эффектами пренебрегают из-за их малости.

Для рассмотрения движения одномерного шара в собственном силовом поле его поступательное (орбитальное) движение относительно фиксированной точки (Солнца) нужно отделить от двух других составляющих движения. После чего можно рассматривать как вращение тела относительно центра масс под действием собственного поля сил, так и движение, связанное с изменением структуры и формы. Такого отделения требует лишь момент инерции шара, который зависит от выбора системы координат. Силовая же функция, определяемая как эффект взаимодействия всех пар частиц массы шара, не зависит от ее выбора [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и     http://zhumal.ape.relarn.ru/articles/2005/120.pdf


то на основе тождества (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     http://zhumal.ape.relarn.ru/articles/2005/120.pdf


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     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

функции, представлены в табл. 1 [15]. Заметим, что значения полярного и осевого моментов инерции одномерного шара связаны соотношением 1р=3/21а .

Таблица  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) имеем


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     lIlihttp://zhumal.ape.relarn.ru/articles/2005/120.pdf

4ж2      3     4яг2

0<-^-<-<--------- *-<\,(15)

4лЯ2     5     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     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

полярный момент инерции можно разложить на составляющие, которые соответствуют однородной по массе плотности и ее неоднородностям. Эти составляющие согласно принципу суперпозиции определяют нормальные и тангенциальные динамические эффекты неоднородного тела. Такое разложение для безразмерных структурных коэффициентов а и /? было выполнено Гарсия Ламбасом с соавторами [17] в нашей физической интерпретации [15]. Для разложения использована вспомогательная функция относительного радиального изменения плотности шара вида

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     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

2Ut=Kt,(23)

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ь     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

Уравнение (25) выражает закон сохранения энергии системы и описывает ее поведение в скалярных интегральных характеристиках Uи Ф. Было показано, что кроме уравнений Ньютона, оно также выводится из уравнений движения Эйлера для сплошной среды, а также из уравнений Гамильтона, Эйнштейна и квантовой механики [15]. В русской литературе уравнение (25) известно как уравнение Лагранжа-Якоби, поскольку при его выводе Якоби использовал тождество Лагранжа (2) для разделения движения системы п материальных точек на движение их центра инерции и относительного движения тех же точек вокруг него. При усреднении движения по времени,  когда   Ф=0,  уравнение  (25)  принимает форму  и  содержание уравнения

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

K = -Ymv2.

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

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

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     121'     http://zhumal.ape.relarn.ru/articles/2005/120.pdf


(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


1                       S2(2                 2      \

Ф

' sB

sinZ + — Јsin2ZH------ sinZI2cos Z-sin  Zl-

2                   2          v                          ;

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


(Ос


4Вп


ип

Л


\<GMr

1 PIS


j-xGpor


(31)


(24)

3/2

(32)

лО р0гке

COtr

У.

1РУе

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     121«     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

р0г = —   p{r)dVr; р(г) - закон радиального распределения плотности; рог - среднее

V

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

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

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

2                          2

_ cotat 0     -лСро

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

Найдено, что безразмерный коэффициент ке е[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у     http://zhumal.ape.relarn.ru/articles/2005/120.pdf


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

IX. _U0 _а] _       0.6

Uа2

Cos©-

(33)

0.66-сс,

IX

0.918,   0=23.5°

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

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

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

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

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

Два тела с массами М и да взаимодействуют по закону притяжения Ньютона (рис.2а). Пусть М»т и R»r, где г - радиус тела т, 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и     http://zhumal.ape.relarn.ru/articles/2005/120.pdf


4     з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.(R»r)

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

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     12il     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

данный эффект был одним из определяющих при образовании Земли и Солнечной системы в целом. Например, теперь значение средней плотности Луны менее 2/3 средней плотности Земли, т.е. рм < 2/3рт. Если предположить, что такое соотношение плотностей сохранялось на стадии формирования двух тел, то, согласно уравнению (35), Луна стала самогравитирующим телом на ранней стадии дифференциации исходной массы Земли. Формирования самого тела из выделявшейся оболочки могло происходить с помощью механизма циклонических вихрей, который в свое время предлагал Декарт и который, как теперь становится ясно, был незаслуженно отвергнут. Если учесть появление тангенциальных сил взаимодействия, возникающих в неоднородных по плотности массах, то такой механизм является вполне реальным [16]. На основе анализа орбит ИСЗ найдено значение структурного коэффициента осевого момента инерции, равное Р j_ = 0,3315 [3]. Тогда из выражений (19) и (20)

полярный радиус инерции равен 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 г/см со скачками значений на границе Мохоровичича, между верхней и нижней мантией и на границах внешнего и внутреннего ядра. Буллен ввел скачки плотности после неудачной, как он решил, аппроксимации сейсмических данных параболической кривой, при которой получалось, что плотность ядра к центру должна падать. Буллен, естественно, не предполагал, что радиус инерции и радиус тяготения тела не находится в его центре масс, а поэтому там не находится и максимум плотности. Теперь, когда наши представления об этих величинах изменились, задача о радиальном распределении плотности требует нового рассмотрения.


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     lIZIhttp://zhumal.ape.relarn.ru/articles/2005/120.pdf

Для поиска возможных путей решения такой задачи нами рассчитано и проанализировано ряд формальных кривых радиального распределения плотности, изменяющейся по параболическому закону вида рг = ро(ах + Ъх + с) ('где x=r/R; а, Ъ, с - численные коэффициенты; ро - средняя плотность тела). Численные коэффициенты подбирались при различных заданных значениях плотности поверхностной оболочки psи таким образом, чтобы масса планеты оставалась постоянной, т.е.

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     lIIihttp://zhumal.ape.relarn.ru/articles/2005/120.pdf

Для построения графиков в безразмерной форме значения потенциала Ur отнесены к Ur=GM/R, а значения силы тяжести qrк qR=GM/R.

Как видно из рисунков 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     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

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

Прецессия и нутация оси вращения

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

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     lIIihttp://zhumal.ape.relarn.ru/articles/2005/120.pdf

каждом обороте оболочки ее угловая скорость пол-оборота ускоряется и вторые полоборота замедляется. При этом кора и океаны за счет инерции массы, будучи во взвешенном состоянии, при ускорении оболочки тормозят, а при замедлении -ускоряют ее движение. Эти эффекты инерции движения коры и океанов мы наблюдаем в форме суточных и полусуточных нутаций. Аналогичный эффект имеем при движении Луны. Силовое поле Земли ускоряет и замедляет движение Луна по орбите и вызывает месячные и полумесячные нутации верхней оболочки (оси) Земли. Тот же эффект наблюдается при движении Земли по орбите вокруг Солнца. Он дает не только годовую и полугодовую нутацию, но и Чандлеровскую составляющую качания оси с периодом 420 дней. Последний эффект пропорционален 27 дневному периоду вращения Луны и Солнца и составляет 365/27 = 14 мес.

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

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

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

Заключение

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


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     111ь     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

[11-13] для потока Хаббла при взаимодействии однородной среды и для получения углового момента при взаимодействии неоднородностеи с однородной массой находят успешное приложение в динамике небесного тела и геофизике в рамках дальнейшего развития и решения задач методами классической механики.

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

Литература

  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. Classical Mechanics, 2nd Ed. Reading-Massachusetts: Addison-Wesley, 1980. 642 P.
  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'     http://zhumal.ape.relarn.ru/articles/2005/120.pdf

Подписи к рисункам

Рис. 1. Отделение момента инерции Земли в собственной инерциальной системе отсчета при ее относительном движении в силовом поле Солнца.

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

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

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

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

(а)


Электронный научный журнал «ИССЛЕДОВАНО В РОССИИ»     122«     http://zhumal.ape.relam.ru/articles/2005/120.pdf


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


 



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.