Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

shimbirev_b_p_teoriya_figury_zemli

.pdf
Скачиваний:
100
Добавлен:
04.06.2015
Размер:
14.99 Mб
Скачать

Рассмотрим решение внутренней задачи Дирихле на конкретном примереПусть заданная на сфере функция имеет вид

/(0, Я) = С + с о з 2 0 .

Представим / (0, Я) как сумму сферических функций. Постоянная С является сферической функцией нулевого порядка, ибо любая постоянная величина может рассматриваться в качестве сферической функции нулевого порядка, поскольку при п= 0 по формуле (111.29) получаем У 0 (0, Я) = А 0 0 . Функцию соз2 0 можно выразить через полином Лежандра второй степени

 

соз 2 0= - | - Р я

(0080) + -^-.

 

Таким образом, заданная на сфере функция может быть следующим

обра-

зом представлена суммой сферических функций

 

 

/(0, х) = С + 4 - + -|Р2 (соз0).

 

С другой стороны, на основании (111.25) можно написать

 

/(0, Я) =

2 г „ ( 0 , Я) = У0 (0,

Я) + 1М0, Я) + У2 (0, Я) + . • .

 

 

п-0

 

 

Следовательно,

должно соблюдаться условие

 

п=О

 

 

 

В силу теоремы о единственности

разложения функции в ряд по

сфери-

ческим функциям должны соблюдаться

условия:

 

У1 = 0,

У2 = - § - ^(соз0),

Уз = У4 = . . . = Уп = . . . = 0.

Подставляя далее найденные значения Уп (0, Я) в (111.22), получим решение внутренней задачи Дирихле

^(Р, 0> я) = 2 - & Г » ( 0 . Я) = С + 4 + -|-т^(соз0).

п=0

Читателю представляется

возможность самостоятельно убедиться в

том,

чт о полученное решение задачи Дирихле значительно проще, чем решение

при

по мощи интеграла Пуассона (П.42)

 

У,.(р, 0, я) =

4 И ( С + С О 8 2 0 ) - ^ 7 Г - Й 0 -

 

 

с

 

90

При решении внешней задачи Неймана для сферы краевое условие имеет

вид

дУе(р, 9, X)

др

Р ^ - т * ) .

(П1.37)

где / (0, X) — заданная на сфере функция. Она раскладывается в ряд (111.25); будем считать сферические функции Уп (0, X) известными. Искомая функция Ге (р, 0, X) — гармоническая во внешнем пространстве и регулярная на бесконечности — должна иметь вид (111.24)

СХ1

п=Ор"

где Хп (0, X) — сферические функции, подлежащие определению. Выполняя дифференцирование функции Уе (р, 0, X) по р, получаем предельное равенство

дУе

др

Р-+В

11=О

 

 

 

 

 

 

 

 

 

 

Краевое условие (111.37) примет

вид

 

 

 

 

2 ^ ± 1 Х „ ( 0 , Л) = 2 ^ „ ( 8 , Я).

 

 

 

 

 

 

п=0

Исходя из теоремы о единственности разложения функции в ряд по сфери-

ческим функциям, получаем

 

 

 

 

 

 

п +я1 - Х „ ( 0 ,

Л) = Г„(0,

X).

(п — 0, 1, 2, . . .)

Отсюда

неизвестные

функции

Х п

(0, X)

определяются через известные

функции Уп

(0, X)

 

 

 

 

 

 

 

 

 

хп(в,

х) =

- - ^ г У п

( в , X)

л решение внешней задачи Неймана будет иметь вид

 

 

 

 

 

со

Н у+1 Уп(9, X)

 

 

у д р ,

0, я,) = - д

2 (

 

 

- г )

и + 1

 

 

 

 

 

п-0

 

 

 

 

 

 

 

 

Пусть, например, требуется решить внешнюю задачу Неймана, если краевое условие задано

 

/(0, Х) = а-\-Ь 8 Ш 2 0 + С 81П 0 соз 0.

Выразим з т 2 0

з т 0 соз 0 через сферические

функции

 

з т 2 8 = - у Р22 (0), з т

0 соз 0 =

Р21 (0);

таким образом,

 

 

 

 

 

 

00

/(0,

Я) = а + - у Р22 (б) +

- у Рц (0)

= 2

 

 

 

п=О

91

На основании сказанного выше находим

У0 = а,

П= о,

У 2 = 4 - 1 ^ 8 2 ( 6 ) + <*а(в)1,

 

 

 

у3 4 =

 

. . . = о.

 

 

При этом решение внешней задачи

 

Неймана получим в виде

 

У Д р , 0, Я) =

а

 

 

[6Р22 (0) +

сР21 (0)].

 

Решение третьей краевой задачи более сложно.

 

Найдем решение внешней задачи

 

 

 

 

 

Уе(р,

В, Я) =

2

(

 

т

 

 

 

 

 

п=0

р

 

 

где сферические функции Хп

(9, Я) неизвестны и подлежат определению. Крае-

вое

условие запишем

 

 

 

 

 

 

 

 

аУДр,

0,

%)+ а т Г '( р 'р е '

к),

где /

(0, к) — заданная на сфере функция и ее разложение в ряд (111.25) сфери-

ческих функций Уп (0, к) считаем

известным.

 

 

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

ных

сферических функций

Хп (0, к) через известные функции У„ (0, к).

 

Для этого получим предельное равенство

 

п=О

и представим заданную на сфере функцию / (0, X) двумя рядами

оо

со

 

/(0, =

Я) = 2

хп (0, к).

Л=0

п=О

 

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

у Л 0 , Я ) =

Х Д О , К),

откуда определяются все Хп (0, к) при условии, что аВ—п—1 Ф 0. Запишем решение внешней третьей краевой задачи

удр, 0, =

Я

Уп (9, к)

л=о

аЯ — п — 1

 

 

 

92

Таким образом, если аВ—п—1 Ф 0, то решение содержит сферические функции всех порядков. Если же при некотором п имеет место равенство, аВп—1 = 0, то для существования решения соответствующая сферическая функция порядка п в краевом условии должна отсутствовать.

Рассмотрим, наконец, внутреннюю третью краевую задачу. Ее решение примет вид

У({р, е,

X) =

 

п=О

Дифференцируя V (р, 6, X) по

р и складывая с а У,- (р, 0, Я), при р В полу-

чим предельное равенство

 

п=О

тогда краевое условие выразится формулой

аЛ + тг

л

п-О

Неизвестные Х„ (0, X) через известные Уп (0, Я) определятся следующим образом:

и окончательно получим

7,<р,

Глава IV

ПОТЕНЦИАЛ СИЛЫ ТЯЖЕСТИ И МЕТОДЫ ЕГО ОПРЕДЕЛЕНИЯ

§ 20. СИЛА ТЯЖЕСТИ И ЕЕ ПОТЕНЦИАЛ. ХАРАКТЕРИСТИКИ ГРАВИТАЦИОННОГО ПОЛЯ ЗЕМЛИ

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

->-

-V

где § — сила тяжести; Р — сила тяготения,

К — центробежная сила.

При

этом будем полагать, что Земля

вращается с постоянной угловой

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

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

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

Потенциал силы тяжести всегда имеет смысл использовать при анализе наблюдений, производимых на поверхности Земли, на самолете, дирижабле

94

п т. п., поскольку при этом приходится иметь дело с характеристиками гравитационного поля Земли. Так, в одних случаях в результате измерений мы получаем значение силы тяжести (гравиметрические измерения), в других же1 приходится использовать в качестве основных ориентирных линий направления силовых линий поля силы тяжести (астрономические и геодезическиеизмерения).

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

Для определения потенциала силы тяжести необходимо знать потенциал центробежной силы, который можно найти, если известны проекции центробежной силы на оси координат. Найдем эти проекции.

Предположим,

что Земля вращается с постоянной угловой скоростью со.

В этом случае на

каждую единичную массу, находящуюся на Земле, будет

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

 

К = со2р,

 

 

где р — радиус параллели

текущей точки М.

 

 

Выберем прямоугольную систему координат,

начало которой

совместим

с центром параллели. Ось

г направим вдоль оси

вращения Земли,

оси ж и у

расположим произвольно в плоскости малого круга, по которому совершается вращение точки М.

Так как направление центробежной силы будет совпадать с направлением радиуса р параллели точки М, то, очевидно, что проекции центробежной силы

на оси координат можно вычислить по формулам

 

Кх

= Ксов(К,

х) = К соз (р,

х);

Ку

= Ксо&(К,

у) = К соз(р,

у);:

К г = Ксо&(К, г) — К соз(р, г).

соз (р,

=

соз (р,

= у

соз(р, г) = О,

и проекции центробежной силы на оси координат будут

Кх = а2х

К у = со2^ к г = 0.

Нетрудно прийти к выводу, что потенциал () центробежной силы будет"

<? = ^-(^2+г/2);

(IV.1)

95

в самом деле, взяв от

функции () частные производные по

осям координат,

найдем, что функция

удовлетворяет

условиям Лежандра

(1.14)

 

дОдх =

со2а; = Кх;

 

4

=

 

 

 

 

да

=

к2г.

 

 

Из формулы (IV. 1) видно, что

на

оси

вращения

потенциал

равен нулю,

а на экваторе имеет максимальное значение.

Оператор

Лапласа

от функции ()

будет

 

 

 

 

 

Д<? =

2.

 

 

 

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

— + = +

8у = ру+Ку=-щ Ь - ^ - ;

Эти соотношения позволяют определить потенциал силы тяжести IV как

1У = У + <?.

(IV.2)

Проверим это утверждение. Дифференцируя (1У.2) почленно по координатам, убеждаемся, что функция IV удовлетворяет условиям (1.14)

д\у _ дх

дШ _

ду

д\У дт, =

Формулу (1У.2) можно представить в ином виде

IV-- ' 1 1 1 1

+ "у-О^2 + 2/2). ( \

'

(1У.З)

Для вычисления потенциала силы тяжести по этой формуле необходимо знать не только форму Земли, но и закон распределения плотностей внутри ее. Поскольку фигура Земли и распределение масс внутри ее неизвестны, точное вычисление потенциала IV по формуле (IV.3) невозможно. Однако если первый

*Т.7ЮГЕ этой формулы

— хготеицшал я-яготежия

V — разложить

в ряд по

сферш-

96

ческим функциям, можно вывести приближенное

значение потенциала ТУ

В следующих разделах будет дан этот вывод.

 

Определим вектор силы тяжести

 

ЗРР т , д\у - , ди> г

/ТЛГ

Таким образом, физический смысл первых производных потенциала IV вполне ясен. Что касается физического смысла вторых производных

32\у

дш

дЧУ

дЧУ

дЧУ

дШ

дх% '

дуг '

дгг '

дхду '

дхдъ '

ду дг '

то о нем будет сказано

ниже.

 

 

 

 

В потенциал силы тяжести в качестве слагаемого входит потенциал притяжения V, поэтому первые производные потенциала IV непрерывны, а вторые производные IV терпят разрыв там, где плотность притягивающих масс ме-

няется

скачком.

 

 

 

 

Из

(IV.2) следует,

что

 

 

 

 

 

ДТУ = Д У + Д ( ? .

 

Во

внешнем пространстве потенциал тяготения удовлетворяет

уравне-

нию Лапласа (1.70), следовательно, во

внешнем пространстве

 

 

 

 

Д1У =

2со2,

(1У.5)

т. е. потенциал силы тяжести не

является гармонической функцией.

 

Во

внутреннем пространстве

V удовлетворяет уравнению Пуассона (1.72)

п поэтому справедливо

соотношение

 

 

 

 

АИ/ = - 4 я / б + 2ю2.

(1У.6)

Для потенциала силы тяжести будут также справедливы формулы (1.35), (1.40), (1.41) и все следствия из них. На основании (1.35) определим дифференциал потенциала силы тяжести

Ш = %С05(8, «)<&. (IV.7)

Из этой формулы вытекает, что приращение потенциала между двумя бесконечно близкими точками может быть определено, если известны: сила тяжести расстояние с1я между точками и угол (§, в), являющийся дополнением зенитного расстояния до 180°. Зенитные расстояния измеряются со сравнительно малой точностью. Поэтому имеет смысл выразить приращение потенциала с11У через соответствующий отрезок йп нормали к уровенной поверхности потенциала силы тяжести.

Итак, проектируя отрезок Аз на направление внешней

относительно уро-

венной поверхности нормали

 

Ап = йзсо5 [180 — (§, в)1 = —соз(^, в)йз,

 

получим

(1У.8)

№ = —§йп.

Величину Лп можно рассматривать как элементарное превышение йк между двумя бесконечно близкими уровенными поверхностями и потому в формуле (1У.8) величину йп можно заменить величиной йк.

7 Заказ 1379

97

Из (1У.8) следует, что •

Ш'

 

 

г — ъ г -

.

( 1 У - 9 )

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

Разность потенциалов между двумя.точками земной поверхности О и А, находящимися на конечном друг от друга расстоянии, на основании (1У.8) может быть получена как сумма элементарных работ, т. е.

 

А

 

Ж А - № 0

= - \ § й к .

(IV.10)

 

о

 

Разность потенциалов Ц'л IV0

*

производимых

находим из измерений,

на физической поверхности Земли. Для этого нужно, чтобы между точками О и А был проложен нивелирный ход, который позволит получить превышения к, и чтобы в точках этого хода были измерены значения силы тяжести При этих условиях интеграл (IV. 10) может быть вычислен.

Важной характеристикой поля силы тяжести являются силовые линии, которые получили название отвесных линий.

Получим в общем виде систему дифференциальных уравнений для всего семейства отвесных линий.

Для этого возьмем на какой-либо отвесной линии элементарный отрезок йк, его проекции на оси координат обозначим через йх, йу и йг,. Согласно опре-

делению вектор силы тяжести § должен касаться отрезка йк. В этом случае

проекции §х , §у ,

вектора д должны быть пропорциональны проекциям йк

 

йх

йу

йг

ИЛИ

ёх

ёу

ёг

 

 

 

 

й х _

йУ

/IV 1 <П

 

дУУ

дУУ

дУУ '

 

дх

ду

дг

Из системы уравнений (IV. 11) следует, что отвесные линии определяются первыми производными потенциала ТУ. Поскольку первые производные потенциала обладают свойством непрерывности, то и отвесные линии также будут выражаться непрерывными функциями. Определим вектор кривизны отвесной линии. Вектором кривизны кривой называют предел Отношения приращения единичного вектора касательной т к приращению йк длины кривой, когда последнее стремится к нулю, т. е.

где V — единичный вектор главной

нормали

кривой; р — радиус кривизны

кривой.

 

 

Единичный вектор касательной к

отвесной

линии

98

отсюда

йх

 

~йк

йк

йк "

Будем считать, что направление касательной к противоположно направле-

нию силы тяжести

В этом

случае,

согласно (1У.9)

 

 

 

д\У

и, следовательно,

 

 

дк

 

 

дт

 

 

 

 

 

дк

дк2

Учитывая, что вектор силы тяжести § определяется формулой (IV.4), получим

йк ~

е

дх дк

ду дк

-

дг

дк

0+М'дх

I

ду

дг

Л

д№

йт _

1

Г дт

г . дт

дт

 

 

 

д\У

т

д]У .

дУУ уН

дЩ

Если ось 2 совместить с направлением силы тяжести, то

 

 

 

дУУ

ду

=0,

дШ

дт

 

 

дт

'

дт

 

дт

дт

дт

дх

 

дг

дх дк

 

 

дх

дг

ду дк

 

ду дг '

дг

дк

дка '

После

сокращений

 

1

_ Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дт

 

-

дт

-

 

 

 

(IV. 13)

 

 

 

Р

~"

8

1_ дх

дг

I •

ду дг

]

 

 

 

отсюда следует, что отвесная линия является пространственной кривои двоякой кривизны.

Поскольку вектор кривизны отвесной линии определяется через вторые производные потенциала И7, то при пересечении отвесной линией масс с раз-

личными

плотностями,

ее вектор кривизны

Касательная плоскость

будет меняться скачком. Однако вне масс вто-

 

рые производные потенциала непрерывны и от-

 

весные линии изменяют свою кривизну также

 

непрерывно.

 

 

 

 

 

Получим уравнение касательной пло-

 

скости к уровенной поверхности. Обозначим

 

через х0 , у0,

20 — координаты точки, нахо-

 

дящейся

на

уровенной

поверхности,

а

че-

 

рез X, У, 2 — координаты текущей точки пло-

 

скости

(рис.

30). Плоскость касается уровенной поверхности в точке

0, у0,

2о).

Рассмотрим

два вектора:

а

и

Очевидно, что

 

 

 

а = ах1 +

ау]

+

 

Эти два вектора взаимно перпендикулярны, поэтому должно соблюдаться условие

ах§х + ау8у + агёг = 0.

7*

99

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]