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

shimbirev_b_p_teoriya_figury_zemli

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

Если в подынтегральной функции (1.20) вместо г подставить гга1п или г„ (рис. 13), то, очевидно, будут иметь место следующие неравенства:

гшах

т

Н'К'Ш' б йт

Но гт а х и гт ! п , как постоянные величины, могут быть вынесены из-под знака интеграла, а интеграл

где М — масса тела.

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате получаются неравенства

 

 

 

 

 

' тах

 

 

' т т

 

 

 

 

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

 

 

 

 

 

и ^ < г У ( Р ) < 1 т . .

 

 

 

 

гтах

 

 

 

' т!п

 

 

 

 

При удалении точки Р в бесконечность

 

 

 

 

 

Нш —=Нш

 

—=

1.

 

 

 

 

гтах

 

 

гщ1п

 

 

 

 

Отсюда получаем, что

 

 

 

 

 

 

 

 

 

Н т г У =

/М .

 

 

 

(1.74)

 

г->-со

 

 

 

 

 

 

Аналогично можно доказать,

что

 

 

 

 

 

 

 

Н т

Г2

- -

--1М.

 

 

 

 

г-> оэ

дУ

 

 

 

 

 

 

 

дг

 

 

 

 

 

Функция,

удовлетворяющая

полученным

трем условиям

(1.73),

(1.74)

и (1.75), называется регулярной

на

бесконечности.

Потенциал

объемных

масс

регулярен на

бесконечности.

 

 

 

 

 

 

 

 

 

П о т е н ц и а л п р о с т о г о

с л о я

 

 

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

Производная потенциала V по произвольному направлению I

Т Г - ' И

а

(1.76,

а

 

40

является проекцией силы притяжения простого слоя на направление I. В том случае, когда притягиваемая точка пересекает поверхность простого слоя, первые производные потенциала терпят разрыв. Мы уже убедились в этом при рассмотрении частного случая простого слоя — бесконечно тонкого однородного сферического слоя (см. § 4). Теперь покажем, что это имеет место и в более общем случае.

Обозначим внешнюю точку относительно слоя через А, внутреннюю — через В, а точку на самой поверхности слоя — через С (рис. 14). Пусть точки А ж В неограниченно приближаются к точке С, в этом случае производные потенциала У по произвольному напра-

влению

I

будут иметь

различные

 

 

пределы.

 

 

 

 

 

 

 

 

Будем называть

внешней произ-

 

 

водной

потенциала

простого

 

слоя

Н46

 

в точке С величину

 

 

 

 

 

 

 

А1

= Н т

ауа

 

 

 

 

 

 

А-уС

А1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= — Н т

/ СС [ХС08(Га,

I) йа

(1.77)

 

 

 

А-УС

 

«>

л

 

 

 

Рис. 14

 

внутренней производной потенциала простого слоя в точке С величину

 

 

 

аУс

=

Н т 1

-

Н т / СС (х соз (гв , I) йа.

(1.78)

 

 

 

61

в-+с

А1

 

 

 

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

внешнюю

и

внутреннюю производные можно

определить

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

Прямое значение производной —йУ^/й/ определим, как составляющую силы притяжения в самой точке слоя. йУ^'/й/ определяется без предельных переходов непосредственно в точке С, расположенной на слое

 

 

АУ$

-'И

Р СОВ (Гс,

йа.

(1.79)

 

 

 

 

 

6,1

 

 

 

 

 

 

 

 

Для

того

чтобы показать,

что при приближении

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

кам С, А и В

(а точки А и В

неограниченно приближаются к точке С), т. е.

когда г

будет

стремиться к нулю, интегралы в (1.77)

(1.78) и (1.79) сохранят

конечное значение, поступим следующим образом. Выделим вокруг точки С

бесконечно малый круг

а 0 с постоянной

плотностью

рс-

Тогда

 

все три

произ-

водные можно будет представить в виде сумм

 

 

 

 

 

<11

 

А-+С

 

г \

 

 

А-^-С

 

 

 

 

 

ЛУе

 

- Н т

с-СГо

|ЛС03(ГА, I)

йа- Н т

цсоз(гА ,

I) йа,

(1.80)

АУУ

 

 

(X СОЗ (Гв,1)

 

Г р,С08(гв, I)

 

 

=

— Н т / СГ

 

йа,

 

61

 

г2

 

йа — Н т / Г С

 

 

 

( 1 . 8 1)

 

В-+С

^ ^

 

 

 

Ов

 

 

 

 

 

 

 

 

О-Оо

с

 

 

 

 

 

 

 

 

 

 

 

 

, 1

 

•Н'

 

г)

 

 

 

 

 

 

 

(X соз с

 

ЦС08(гс,

 

 

 

 

 

 

 

 

г2

)

йа-

 

 

йа.

 

(1.82)

41

Первые слагаемые в формулах (1.80), (1.81)

и (1.82)

всегда

ограничены,

так как в области о — сг0 г нигде не обращается в нуль.

 

равенства

При совпадении точек А и В с точкой

С будут иметь место

А^Н-тС ои-а.

цсо8(гЛ,

I) йа =

в^сН т

ои-о„

ЦС08(/-В, I)-йа -

и[.1С08(гс,

I) йа. (1.83'

Вторые слагаемые в (1.80), (1.81) и (1.82) можно рассматривать как соста-

вляющие притяжения плоского круга сг0 по направлению I.

 

(1.63), следует

Значения

этих составляющих

определяются

по формулам

только заменить в них направление

оси

т,

направлением

внешней нормали п

к поверхности а в точке С. Подставляя (1.63) в

(1.80),

(1.81) и (1.82) и считая,

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

 

 

 

 

 

41

^

цсоз(гс ,

I) йа — 2я/р с соз (п,

I)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4У^

А \а

ЦС08(гс,

I)

йа-|-2я/рс соз (п.

I)

 

(1.84)

 

си

 

 

 

 

 

 

 

 

 

 

 

 

 

-'И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

йУк

 

(Д,С08(/'С,

I)

йа

 

 

 

 

 

 

 

 

41

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или опуская индексы

С,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е

ауо

 

•2я/рсоз (п,

I)

 

 

 

 

 

 

41

41

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.85)

 

 

йУ{

4У(1

-р2я/рсоз (п,

I)

 

 

 

 

 

 

 

 

 

Отсюда получим

41

41

 

 

 

 

0

1

VГ4Уее .

4У14УГ1

 

 

 

 

 

 

 

 

 

(1.86)

 

 

41

2

[_

41

'

41 _Г

 

 

 

 

 

 

 

 

 

Эта формула называется формулой

Плёмели. Пользуясь формулой Плёмели,

мояшо получить прямое значение производной на слое йУ0/й1 через предельные значения.

Вычитая из первого соотношения в (1.85)

второе, получим формулу Пуас-

сона

е

АУ,

 

 

 

 

 

 

-4я/рсоз(и, I).

(1.87)

 

41

41

 

 

 

 

 

 

Формула Пуассона

позволяет

вычислить

величину скачка производной

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

поверх-

ности слоя.

 

 

 

 

 

 

Если производную

потенциала

брать по

направлению внешней

нормали

к поверхности слоя, то соз (п, I) =

+ 1 и получим

 

 

 

е

4УI

—4я/р.

(1.88)

 

 

4п

4п

 

 

 

 

 

Итак, первые производные потенциала силы притяжения простого слоя терпят разрыв при пересечении слоя. Поэтому па простом слое всегда следует различать три значения производных потенциала: два предельных и одно прямое.

42

П о т е н ц и а л д в о й н о г о с л о я

Потенциал двойного слоя можно представить как сумму производных иенциалов простых слоев (1.33). В самом деле, интегралы, стоящие справа

1.33),

являются потенциалами простых

слоев с плотностями

—V сов (п, х),

-V соз (п, у) и —V соз (п, г).

 

 

Во

внешнем пространстве производные

потенциала простого

слоя непре-

^вны, удовлетворяют уравнению Лапласа (1.70) и регулярны на бесконеч- :сти. Следовательно, во внешнем пространстве этими свойствами будет область и потенциал двойного слоя.

Однако при пересечении притягиваемой точкой поверхности слоя первые :оизводные потенциала простого слоя терпят разрыв, поэтому и потенциал ЕОЙНОГО слоя также будет терпеть разрыв. Так же, как и для первых произлных потенциала простого слоя, следует различать два предельных е и И7,- олно прямое значение Ш0 потенциала двойного слоя. Величина разрыва : генциала двойного слоя в точках самого слоя определяется из соотношений

е = 1У0 + 2я/у

(1.89)

§ 6. УГОЛ ВИДИМОСТИ И ИНТЕГРАЛЫ ГАУССА

Если двойной слой имеет постоянную плотность /V = 1, то его потенциал :>значим через Й. В этом случае формула (1.30) примет вид

„ .

ч

=

Г Г соз (г, п)

Я (г,

 

йа.

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

Рис. 15

круг точки Р (х, у,

2.)

радиусом, равным единице длины. Обозначим

пло-

.^дку участка, вырезываемого из единичной сферы конусом, вершиной

кото-

го является точка Р

(х,

у, г), а основанием — элемент площади йа, через й<л.

43

Величина йа называется элементарным телесным углом или же углом видимости, под которым элемент йа виден из точки Р. Для вычисления йен необходимо элемент йа поверхности а в точке М (рис. 15) спроектировать на сферу с центром в точке Р и радиусом г, что даст йа соз (г, п), и поделить на г2, так что

,

=

Ла соз (г,

п)

а &

— ^ — >

- .

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

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

а.

 

Вычислим значение потенциала

& в трех точках: внутренней относительно

замкнутой поверхности а, внешней и на

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

На основании предыдущего можно написать

& (ж, у, =

а.

 

 

сг

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

с/0 соз (/•, п)

—йа.

Г 2

 

Сумма углов видимости, под которыми площадки йа усматриваются из

точки Р, равна 4л и, следовательно,

 

й(ж, у, г) = — 4я.

Если потенциал вычисляется во внешней относительно поверхности а

точке Р

(х, у, т), то каждый элементарный конус, вершиной которого

является

точка Р,

вырезает на этой поверхности два элемента: йаг и йа2 (см.

рис. 15),

причем угол между вектором г и внешней нормалью п будет острым в точке Мг

и тупым в точке М 2 . Следовательно, проекции элементов йах,

йа2 на единичную

сферу

будут равными по величине, но противоположными

по знаку. Поэтому

 

-^р-соз^!, т^) -(-

соз (гс2, г2) = йсо — йа = О,

где гг

= МгР, г2 = М2Р (см. рис.

15).

 

 

Таким образом, для внешней точки получаем

 

 

Й(ж, у, 2) =

0.

 

При вычислении потенциала в точке Р, находящейся на самой поверх-

ности а, будем иметь (см. рис. 15)

 

 

 

 

с!а соз (г, п)

,

 

 

г2

—=

— йа,

 

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

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

 

О (х, у, г) =

—2я.

Полученные

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

О (х, у, г) в

различных точках, называются

интегралами Гаусса. Как видно из

предыдущего, интеграл Гаусса О (х, у, г) испытывает скачок, равный -)- 4л, при перемещении точки Р (х, у, г) изнутри наружу любой замкнутой поверхности а.

Глава II

КРАЕВЫЕ ЗАДАЧИ ТЕОРИИ ПОТЕНЦИАЛА

§ 7. ФОРМУЛЫ ГРИНА

Прямая задача теории потенциала заключается в определении потенциала тяготения V, если известны строение тела и его форма.

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

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

В прямой задаче теории потенциала эта зависимость определяется формулой (1.20). Эта формула не позволяет дать положительный ответ на поста-

вленный

вопрос.

Но

рассмотрим обратную задачу: выразим плотность через потенциал.

Эта задача может быть решена с помощью уравнения Пуассона

 

ДУ = —4я/6.

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

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

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

Рассмотрим в качестве примера однородный шар. Как видно из формулы (1.54), его внешний потенциал определяется массой М и, следовательно, при любых значениях плотности б и радиуса Я, при условии, что общая масса остается неизменной, внешний потенциал V будет иметь одинаковые значения.

Потенциал вне притягивающих масс удовлетворяет уравнению Лапласа, которое не содержит плотности. Для его решения необходимо знать область, в которой потенциал удовлетворяет уравнению Лапласа.

45

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

В основе современных методов решения уравнения Лапласа лежат формулы преобразования некоторых объемных интегралов в интегралы по поверх-

ности, полученные Грином. Дадим их вывод.

 

Пусть имеются непрерывные функции Т]1 (х, у, г) и

(х, у, г), имеющие

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

Рассмотрим тождество

 

 

д

( и

дУ1\_дЦ1

 

дУ1

д*У[

 

 

 

 

дх

\

'дх у

дх

дх

'

1 дх%

 

 

Проинтегрировав по переменной х, получим

 

 

 

 

 

Х\

 

 

 

XI

Я1

 

 

 

Интегрируя по двум другим переменным, будем иметь

х2

Уг

г2

 

7/2 г2

 

хг

 

х2 у2

 

г,

111

 

Ио

. ^

I

 

11

$*

XI

г/1

21

 

У1 21

 

XI

 

Ж1 1/1

21

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

нием области, т на плоскость уг (рис. 16).

 

Допустим, что объем т пересекается любой прямой, параллельной оси

х,

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

2,

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

х,

пересекает поверхность а. Элемент йу йг. плоской области

2 можно

рассматри-

вать как проекцию элементов йах и йо2

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

которые

вырезаются

на этой поверхности бесконечно тонкой

призмой, имеющей основание йу йг.

46

В этом случае йу йг =

—Й0Х

соз (пх, ж) =

йа>2 соз (га2>

Угол

(гах, ж) —

>тйл между внешней нормалью п1

к элементу (1а1 и осью х

всегда будет тупым,

а угол (пг, х) — всегда острым.

 

 

 

 

 

 

 

 

 

 

 

2

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

: =>?сти к интегрированию по поверхности сг

 

 

 

 

 

 

 

 

 

 

 

 

 

V I 21

 

 

 

 

 

 

У1 21

 

 

 

 

 

 

 

 

 

 

 

 

• ^

111

- соз (п,

х)йа.

 

 

 

 

 

 

Таким

образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тт д2У1 Л

 

 

 

Напишем еще две аналогичные формулы,

заменяя в

них

х

на

у и на

г

и все три формулы сложим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ас/г

дУ{

,

91/;

дУI

,

дцI

 

97/ •

 

 

д*У1

азV,

, 927/

йт =

дх

дх

1

%

 

 

 

 

 

 

 

 

 

 

91/2

 

322

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЗУ/

СОЗ (п,

у)-{--~С08(п,

 

2)

 

 

 

 

 

 

 

 

 

 

 

 

 

ду

 

 

 

 

 

 

 

 

 

 

На основании

(1.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЛУ{

дУ{

 

 

ж) +

ЗУ/

соз (и,

 

ЗУ/

 

г),

 

 

 

 

 

 

 

^ - с о з ( п ,

 

у) + -5 -г-сов(п,

 

 

 

 

поэтому

 

 

Ж

 

 

 

 

дХ]ь

ду(

,

ди{

дуг

 

 

 

 

 

 

 

 

 

 

дУг

ЗУ/

 

 

 

 

 

 

 

 

 

 

 

 

дх

 

дх

 

ду

ду

1

Зг

Зг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

йп йо.

 

 

 

 

Обозначим

3[7/

ЗУ/ .

3^7/

ЗУ;

, 9(7/

ЗУ/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх

дх

 

ду

 

ду

Зг

 

 

 

 

 

 

 

 

Из этого определения

следует

 

 

 

 

 

 

 

 

 

 

 

 

 

так как перестановка функций

 

и У,- сводится здесь к перестановке

сомножи-

телей в трех слагаемых. При [ / ; =

сонзЪ

или У / =

сопзЪ Б

(?7/, У;) =

0. Также

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

щи»

В

оба эти выражения представляют собой сумму квадратов первых производных от и У,-.

Далее, сумму вторых производных от Уг

по х, у, г обозначим по-прежнему

через ДУ,-, так что

С)2VI

д*Ус д*У1

д.

 

 

дх* '

ду2 "Т"

I-

Тогда первая или, как говорят иногда, предварительная формула Грина

будет иметь вид

 

 

 

х

 

х

а

Отметим, что в последнем интеграле производная берется по направлению внешней нормали к поверхности о, восстановленной на элементе йсг.

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

Положим теперь в (II.1) 17г = У,-, это дает

 

 

X

X

 

 

о

 

Второй интеграл слева имеет важное значение в теории потенциала; он

носит название интеграла Дирихле.

Обозначим его через 1

 

 

 

 

 

 

 

йх.

 

Все элементы I, а следовательно, интеграл в целом, суть величины суще-

ственно

положительные, /

может

быть равен нулю, в случае, если

каждый

из его

элементов обратится

в нуль, но это возможно только при 17{

= сопзЪ

во всех точках области т.

 

 

 

 

Таким образом, /

^ О, причем

знак равенства относится только к случаю,

когда 171 =

сонзЪ. Вследствие этого для

всякой функции 171 имеет место не-

равенство

 

 

 

 

 

 

 

 

 

1

 

X

 

 

Положив

в (11.1)

27,- =

сопзЪ, найдем,

сократив на 171

 

X а

Наконец, поменяем в (II.1) функции 17{ и У,- местами, получим

 

 

+

 

Щ Ъ - Ц у

^

А а .

 

 

(П.5)

Т

 

Т

 

 

 

 

 

 

Вычтем из (11.1)

(II.5)

 

 

 

 

 

 

 

 

1 1 $ Ш

~

*

= $ $

О

(

П

.

6

)

т.

 

 

о

 

 

 

 

 

 

Полученная формула называется второй формулой

Грина.

 

 

 

48

Выведем теперь формулы Грина для функций IIе и Vе, непрерывных вне поверхности сг, имеющих вне этой поверхности ограниченные вторые производные. Кроме того, для функций Vе и Уе ставится дополнительное (и обязательное) условие: они должны быть регулярны на бесконечности, т. е. удовлетворять условиям (1.73), (1.74) и (1.75). Проведем радиусом В сферу 2 так, чтобы поверхность а была бы целиком заключена внутри сферы 2 (рис. 17), и к области т, ограниченной двумя поверхностями сг и 2, применим формулу (II.1)

х

о

2

Элемент поверхности сферы й2 = В2 з т 0Й0ЙА, и потому подынтегральное выражение во втором интеграле правой части можно представить в виде

^е Ф -

У = ч д * 4 т г З1п е ю ак,

е йп

ал

так как направление внешней нормали к поверхности 2 совпадает с направле-

нием радиуса сферы В. При В

на основании свойства

регулярности на

бесконечности,

имеем

 

 

 

 

 

В 2

сШе =

сопз!,

а Н т V'

= 0.

 

 

 

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

при В,

стремящемся

к беско-

 

нечности, предел интеграла

 

^

равен

ну-

 

лю. Поэтому

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

И.7)

Рис. 17

Отрицательный знак перед интегралом в правой части поставлен потому, что за положительное направление нормали по-прежнему будем считать направление во внешнее пространство (внутрь т). Меняем Ле и Уе местами

е йп

(".8)

 

х

а

Вычитая из (II.7) (II.8), получим

V.

йУе

•V.

йЦе

Ло.

(П.9)

 

 

 

 

йп

 

 

Получим формулы Грина для функции 111 = 1/г, понимая под г расстояние между текущей точкой М и данной точкой Р. В этом случае формула (11.6) примет вид

 

 

йп

 

 

а

Для точки Р, лежащей вне объема т,

 

 

 

( П . Ю )

х

о

а

4 Заказ 1379

 

49

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