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

Литература / Закатов П.С. - Курс высшей геодезии (1976)

.pdf
Скачиваний:
733
Добавлен:
10.06.2017
Размер:
34.58 Mб
Скачать

прибору данной конструкции. Точность определения относительного ускоре­

ния силы тяжести в этом случае зависит только от случайных ошибок в опреде­

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

и постоянных прибора, которые могут произойти после наблюдений на исход­

ном пункте.

Применением соответствующей методики и конструкции прибора ослаб­

ляется влияние и этих ошибок.

Точность окончательного вывода периода качания маятника после учета всех поправок характеризуется ошибкой около ±0,00000004 с; соответству­

ющая ошибка в определении ускорения силы тяжести равна приблизительно

±0,2 мгал.

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

пружины.

Пусть мы имеем тело с некоторой постоянной массой; вес этого тела опре­ деляется силой тяжести. Определим вес этого тела р в некотором исходном

пункте, для которого известно значение силы тяжести g 0

Вес этого тела в другой точке, для которой требуется определить силу

тяжести, будет иной; назовем его через р 1 , а искомую силу тяжести в опреде­

ляемой точке - через g 1. Очевидно,

изменение веса тела

1 - р) = Лр яв­

ляется следствием изменения силы

тяжести (Лg) = g 1 -

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

будем иметь:

(Лg) = f (Лр), g1 = go + (Лg).

Простейшая возможная зависимость между (Лg) и (Лр) имеет вид

(Лg) = а Лр+ Ь.

Определение постоянных а и Ь может быть произведено посредством изме­ рения веса теJш в трех точках А, В и С сила тяжести для которых известна. Если обозначить:

gв -gл = (Лg1) и РвРА = Лр1, gc - gв = (Лg2) и РеРв = Лр2,

то

(Лg1) = а Лр1+ Ь, (Лg2) = а Лр2 + Ь,

из которых легко находятся неизвестные коэффициенты а и Ь.

:Измерение веса тела должно быть произведено при помощи прибора, работа которого не зависит от силы тяжести; таким прибором могут служить

пружинные весы, так как упругая сила пружины не зависит от силы тяжести.

Определение силы тяжести при помощи гравиметров занимает мало вре­

мени и производится с высокой точностью.

§ 53. Потенциал силы притяжения

Пусть в пространстве имеем две материальные точки А и В с координатами

а, Ь, с их, у, z (рис. 103). Обозначим массу точки А через т, массу точки В

примем равной единице. Расстояние между точками обозначим через r.

230

Согласно закону всемирного тяготения, сила взаимного притяжения точек А и В :выразится формулой

 

 

F=f.!!!._.

 

(53.1)

 

 

. r2

 

 

Расстояние

r

равно

 

 

 

 

2 = (а-х)2 + (Ь-у)2 + (c-z)2

 

(53.2)

Обозначим

составляющие силы F по осям координат

через Fx, Fy, F2 ,

а углы, образуемые вектором ВА с осями координат, -

сх.,

~' у.

Очевидно,

Fx,

Fy, Fz будут проекциями вектора F

на оси координат ох,

оу, oz.

 

 

 

 

Из рис. 103

 

 

 

 

 

Fx=F cosa}

 

 

 

 

Fy =F cos ~

 

(53.3)

Fz=F cosy

(на рис. 103 для простоты построения начало координат совмещено с одним

из концов вектора F без изменения направления осей). Но

а-х )

cosa=-r- (

 

cos~=Ь--rу- 1·

(53.4)

c-z

 

соsу=-r-

 

Подставляя (53.4) в (53.3), получаем

 

а-х

 

F х = F cos а= fт ----,-:з-

 

Ь-у

(53.5)

F у= F cos ~ = f т ---,:з--

c-z

Fz = F cos у= fт --,:з

Сила притяжения и создаваемое ею напряжение является в е к т о р о м,

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

Поэтому характеристика силового поля в пространстве выражается тремя

уравнениями (53.5).

Однако при определенных условиях поле сил может быть выражено одной

функцией от координат х, у, z, как независимых переменных.

Возьмем функцию

V = f !!!._

(53.6)

r

 

и найдем частные производные ее по координатам х, у и z, явно входящим в нее

через r,

 

 

 

 

 

дV

 

т dr

 

 

дх= -! -;:i dx.

 

Производная

дr

 

 

 

дх вычисляется путем дифференцирования выражения (53.2)

 

дr

=

а-х

(53.7)

 

дх

-- r- ,

 

 

 

 

231

О··.

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

 

 

дV

 

а--х

- =fm -- ·

дх

'

1·3 '

аналогично этому

 

 

~=fm Ь-у }

 

ду

 

r3

(53.8)

дV _ f

c-z ·

--т--

 

дz

,

r3

 

Сравнивая (53.8) и (53.5), получаем

 

 

!!]!_ = F

х

)

 

дх

/

 

!]!_=F

у

\

(53.9)

ду

(.

!!I_=F

z

1

 

дz

J

 

Ф у н к ц и я V, ч а с т н ы е п р о и з в о д н ы е к о т о р о й п о

прямоугольным координатам притягиваемой точ­

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

к о о р д и н а т, н а з ы в а е t;q, с я п о т е н ц и а л ь н о й ф у н к ц и е й,

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

/

, - - - -

- -

-;i

'', /

 

 

/

/

\

/

f-/ - -

'

А/(а.ь.с) 1

Fz

 

 

 

1

1

 

\

 

1

1

 

 

 

 

\

1

 

 

 

1

 

 

 

\

 

 

1

/

/

 

 

 

 

'

1

-~'_J/

/

 

1

 

у

 

Рис. 103

 

Рассмотрим более общий случай,

когда точка притягивается некоторым

телом. На рис. 104 изображено тело 't,

о

!J

Рис. 104

создающее вокруг себя силовое поле притяжения. Определим напряжение этого поля в точке А. Разобъем объем тела 't на элементарные объемы.

dт: = da db dc.

(53.10)

Обозначим через б плотность в единице массы текущей точки М, т. е.

s:.=

dm

(53.11)

u

Тогда будем иметь

 

 

dm = б dadb dc.

(53.12)

232

Аналогично с (53.5) проекциями силы притяжения элементарной массы в точке Jvl на точку А будут:

dF

= f ~- x) dm = f

6 (а:..._х)

da db dc 1

 

х

r3

 

 

r3

 

1

 

dF

= f (ь-

У)

dm = f O(Ь-

У)

da db dc }

(53.13)

У

r3

 

 

r3

 

,

 

dF

= f (с-

z)

dm = f

6 (с-

z)

da db dc 1

 

z

r3

 

 

r3

 

J

 

Суммируя действие элементарных масс по всему объему тела, получаем

da db dc

1

 

dadb dc

1

(53.14)

 

dadb dc

J'

 

Потенциал притяжения тела М на точку А

выразится формулой

V = f SIS~ da db dc.

(53.15)

 

дV

дV дV

В этом легко убедиться, если вычислить производные дх'

ду и дz' кото-

рые будут равны соответственно Fx, F_y, FZ'

 

 

Вводя в последнее выражение массу, можем

написать

 

 

 

(53.16)

В (53.16) интегрирование выполняется по всему объему тела 't. Следова­

тельно, пределы интегрирования определяются в зависимости от формы тела. Если под телом понимать Землю, то выражение (53.16) будет представлять потенциал притяжения Земли на внешнюю точку; в этом случае интегрирова­

ние должно выполняться по всему объему Земли.

При бесконечном удалении тела от притягиваемой точки, т. е. когда рас­ стояние r н'еограниченно велико по сравнению с размерами тела, выражения для силы притяжения ·и потенциальной функции напишутся:

Р= JMr2

}

(53.17)

JM

.

V=r-

 

 

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

Из последних выражений имеем:

lim V = О;

IimrV = fM;

lim r2 /

~~ 1 = fM.

(53.18)

r-oo

r-oo

 

 

 

Функции, удовлетворяющие равенствам

(53.17),

называют

ре гул яр -

ными на бесконечности.

До сих пор мы полагали, что притягиваемая точка находится вне притяги­ tшющего тела. Допустим: теперь, что притягиваемая точка расположена внутри

233

притягивающего тела. Тогда величина r в выражении (53.15) может при­

нять малые значения и стремиться к нулю. При r, стремящемся к нулю, его

значение будет становиться веJiичиной первого порядка малости. Следова­

тельно, выражение

dп db dc r

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

.содержащие интеграл вида

(53.19)

13 этом случае будут конечными и непрерывными, так как, если (а - х) будет

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

азнаменатель - третьего. Поэтому отношение их будет первого порядка

малости.

Отсюда вытекает существенный вывод: потенциал пр и тяже ни я

Земли и его первые производные всюду конечны

и однозначны; можно доказать, что они будут также непрерывными.

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

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

§ 54. Разложение потенциала земного притяжения в ряд

Использование выражения для потенциала jсилы земного притяжения

ввиде

встречает известные трудности. Более удобное выражение для потенциала V

можно получить путем разложения 1 /r в ряд.

 

z

 

Примем систему пространственных прямо-

 

 

угольных координат

с началом в центре Земли

 

 

и с

осью z,

совпадающей с осью вращения

 

 

Земли

(рис.

105). Тогда плоскость

ху совпа­

 

 

дает с

плоскостью земного экватора.

 

 

 

 

Напишем

выражение для потенциала при­

 

 

тяжения

 

 

 

 

 

 

 

 

 

 

V=tS~

(54.1)

 

 

 

 

 

 

 

р' '

 

!/

 

где

р'

-

расстояние

притягиваемой

точки А

 

Рис. 105

(а,

Ь,

с)

от текущей

точки М (х, у,

z) с эле­

 

ментарной массой dm.

 

 

 

 

Обозначим расстояния от точек А и М до начала координат О соответст-

венно r

и р.

 

 

 

 

 

 

 

234

С введенными обозначениями имеем:

р•1 =(а-х)2 + (Ь-у)2+ (c-z) 2 = а2 + ь2 + с2 -

- 2 (ах+ Ъу+ cz)+x2 + у2 + z2 , r2 = а2+ь2 +с2 } •

 

 

 

 

 

р2 = х2+ у2+ z2

 

 

 

 

 

Тогда

 

 

 

 

ах+Ьу+сz + р2

-

 

 

 

2р ах+Ьу+сz + р2)

 

р

~2

r:: r

2

-

2

(

1 -

.

 

 

2rp ----

--' r

 

 

-

r

--------'-

-

 

 

 

 

 

rp

 

 

 

 

 

rp

r2

 

Из треугольника ОМА (см. рис. 105)

р'2 = r2 - 2rp cos \jJ + р2 .

Сравнивая (54.4) с (54.5), находим

ах+Ьу+сz

1 ,~

---- =COS't" .

rp

 

Подставляя (54.6) в (54.4), получаем,

 

р~2 = r2 ( 1 - 2 ~ cos \jJ + ~:) .

Для получения выражения i1Т' входящего

в (54.1), возведем (54.7)

(54.2)

(54.3)

(54.4)

(54,5)

(54.6)

(54. 7)

в сте-

пень - ~ и разложим правую часть по биному Ньютона, тогда :получим ряд

по возрастающим степеням ...rЕ. 1

Ограничиваясь членом с ~:, получаем

1

1

(1

р

.

р2 (3

cos

2

}

(54.8)

-, = - {

+ -

cos 'Ф

+ - 2

 

\j)- 1) + ... .

р

r

t

r

 

r

 

 

 

 

Попутно укажем, что коэффициенты при степенях ..,е._ являются так павы­

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

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

значение --;- в

(54.1), находим

V =; Sdт+ ; 2

р

 

SР cos ~'dm+ 2; 3

SР2 (3 cos2 \j)-1) dm+.. .. (54.9)

Первый интеграл, распространенный па весь объем притягивающего тела

будет равен его массе М.

Второй интеграл, принимая во внимание (54.6), можно представить таю

--f

spcos 1 ,~dm=f

- sр -ахdm+!- sрЬу-'dт+-f

spcz-dm.

(54.10)

r2

r2 rp

r2

rp

r2

rp

 

Но координаты центра массы тела (центра инерции) определяются фор­

мулой

 

 

 

r xdm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x•=_.J__

·

 

 

(54.11)

 

 

 

 

Sam

 

 

 

 

 

 

 

 

 

235

Аналогичные выражения получаются для других осей координат.

Так как начало координат нами было принято в

центре массы Земли,

то х' = О. Принимая во внимание, что

 

Sdm=M,

(54.12)

получаем

 

Jxdm=O.

(54.13)

Аналогичные формулы получаются по другим осям

координат. Тогда пер­

..вый интеграл правой части (54.10) будет

 

!2 5Р ;; dm = 5ахdm = :.: 5хdm = О.

По тем же соображениям будут равны нулю и остальные два интеграла

(54.10).

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

Jpcos-фdm=O. (54.14)

Таким образом, второй член разложения потенциала силы притяжения

вряд равен нулю.

Рассматривая третий интеграл выражения (54. 9)

5p2 (3cos2 ч,-1)dm= 5{r~ (ax+by+cz)2-p 2 }dm,

(54.15)

замечаем, что при возведении в квадрат трехчлена (ах+ Ьу + cz) появится,

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

 

Sах Ьу dm = аЬ Sху dm.

(54.16)

Интегралы вида Sxydm называются произведениями инерции. При рас­

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

равны нулю. Поэтому (54.15) можно представить в виде

5{r~ (ах+Ъу+cz)2 - р2} dm = /2 5{3 (q,2x 2 +Ь2у2+c2z2) - p2r 2 } dm =

 

= /2 S{3 (а2х2 +Ь2у2 +c2z2 ) - р2 2 +Ь2+ с2)} dm.

(54.17)

Соединяя вместе члены, содержащие а2 , затем Ь2 и с2, разбиваем интеграл

на три части, из которых первую напишем так:

r\ Sа2 (3х2-

р2)dm = ;: S(2 -

у2- z2) dm.

(5/2.18)

Аналогично две другие

части представятся

интегралами:

 

 

~~ S(2у22- z2) dm,

 

;: S(2z2 - x2 -y2 ) dm.

Из механики известно, что главные моменты инерции А, В, С относительно

осей координат ~выражаются формулами:

А= j'(y2 +z2) dm

С=J(х2+у2)dm } ·

В= S(x2 +z2 )dm;

 

(54.19)

236

Легко проверить, что

f:::=::=::;::::::=:~ \·

(54.20)

 

J(2z 2 - x2 - y2 )dm=A +В-2С

Поэтому, принимая во внимание (54.18) и (54.20), интеграл (54.15) при­

мет вид

 

} 2 2 (В+ С-2А)+ Ь2 (С+А-2В)+ с2 (А+ В-2С)].

(54.21)

Введем теперь геоцентрические координаты - широту Ф и

долготу L,

()тсчитываемую от плоскости xz. Тогда на основании (4.15) и (4.32) будем иметь:

а = r cos Ф cos L ]

 

 

 

 

Ь = r cos Ф sin L .

 

 

 

(54.22)

 

 

 

 

c=rsin Ф

 

 

 

 

 

Далее, заменяя

cos2L

и sin2L через

косинусы кратных дуг,

получаем:

 

 

а2 =:1 r2 cos2 Ф(1+cos2L) 1

 

 

 

 

 

Ь2 =

2

r 2 cos2

Ф (1- cos 2L)

·

 

 

(54.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2 =r2 sin2 Ф

 

 

 

 

 

 

 

 

Заменяя в (54.21) величины а2 ,

Ь2 и с2 через их выражения (54.23), группи­

руя члены,

содержащие cos 2L, после некоторых

преобразований

получаем

Еыражения для третьего интеграла в (54.9)

 

 

 

 

 

'

А +в)

(1-3sin

2

 

3

 

2

Ф cos 2L.

(54.24)

 

( С--2-

 

Ф)+ 2(B-A)cos

 

Принимая во внимание (54.24), а также (54.14), получаем разложение

потенциала V в ряд по (54.9) в окончательном виде

 

 

 

V= fM +~ (с- А+В) (1-3sin2 Ф)+~ (B-A)cos2 Фcos2L.

(54.25)

r

2r,З

2

 

 

 

 

 

4r3

 

 

 

 

Первый член полученного выражения

(54.25)

представляет собой потен­

циал шара с массой, равной массе Земли.

Второй член, зависящий от широты, представляет собой влияние сжатия

:Земли или, иначе говоря, влияние возрастания масс от полюса :к экватору.

Третий член представляет собой влияние неравномерного распределения масс

no долготе.

Таким образом, полученная формула (54.25) справедлива и для трехос-

ной Земли. '

Если рассматривать Землю как тело вращения и моменты инерции относи­

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

-собой, т. е. А = В, то третий член в формуле (54.25) пропадает; соответственно

изменится в этой формуле и второй член, и для данного случая формула (54.25)

шерепишется

 

V= f: +(C-A)(1-3sin2 Ф).

(54.25')

237

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

Перепишем (54.8) так:

1

1 [

1

р

cos '1' +

( р )2 (3

 

 

1)

+ ...

р' = r

+r

r

2 cos

2

'1'- 2

 

 

 

 

 

 

 

 

 

 

 

 

 

+ (~урп<cos '1') +...J.

 

(54.26)

Ряд (54.26) сходится

при

р ~r. Рп (cos ,J,) -

 

выражение, являющееся

функцией cos ,J, порядка п, представляет собой один из видов сферических функций; это многочлены степени п от cos ,t,. Они называются многочленами Лежандра.

В более общем ви,тQ формула (54.26) перепишется

 

n=-oo

 

 

1

~ pn

(54.27)

-, =

r

п+1 Рп (cos -ф).

р

 

 

n=O

Выражение для потенциала притяжения, принимая во внимание (54.1)

и (54.27), в общем виде напишется:

n=oo

 

V=f Sd~ =f s~ ~:1 Pп(cos,t,)dm

(54.28)

n=O

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

(54.8):

V=f sа; +f s~; cos-фdm+fs~: ;(3cos2 ,t,-1)dm+

(54.29)

Интегрируя первые три члена, получим формулу (54.25), конечно, без

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

fM

f

(

А +в)

2

Ф)+

V=-+-

 

С--- (1-3соs

 

r 2r3

 

2

 

 

+ ::3

(В-А)sin2 Фcos 2L.

 

(54.30)

Последовательное изучение

вопросов программы

 

данного курса дается

без использования сферических функций. Приведенные самые начальные и эле­

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

для общего ознакомления, которое может быть полезно при более детальном

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

трудам (например, [1], [12]; [39] и многие другие).

238

§ 55. Основные свойства потенциала притяжения

Пусть задана материальная точка В с координатами х, у, z; на бесконечно

малом расстоянии ds от В возьмем другую материальную точку В 1

с координа­

тами х + dx, у + dy и z + dz.

 

 

Очевидно,

 

 

 

 

(55.1)

но

 

 

dx = ds cos (s,

х) }

 

dy .= ds cos (s,

у) ,

(55.2)

dz = ds cos (s,

z)

 

где (s, х), (s, у) и (s, z) - углы, образуемые направлением ds с осями коорди­

нат х, у, z.

Тогда потенциал V силы притяжения, как функция координат х, у, z, получит приращение при перемещении из точки В в точку В 1 , равное полному дифференци.алу

дV

дV

дV

(55.3)

dV =дх dx+75y dy+ дz dz.

Или, имея в виду (53.9),

 

 

 

dV =Fxdx+F0 dy+F2 dz.

(55.4)

Принимая во внимание (55.2), (55.4) и (53.3),

 

dV =F ds {cos (s, х) cos (F,

х) +cos (s, у) cos (F,

у)+

+ cos (s, z) cos (F, z)}.

(55.5)

Но выражение, стоящее в фигурных скобках, есть косинус угла между

направлением силы тяжести F и направлением элемента ds, т. е. cos (F, s).

Поэтому

 

dV = F ds cos (F, s).

(55.6)

Так как F cos (F, s) = F s - составляющая силы F по направлению эле­

мента ds, то

 

dV =Fs ds.

( 55. 7)

Из равенств (55.6) и (55. 7) вытекает ряд важных свойств потенциальной функции.

Rак известно из механики, элементарная работа силы F при перемещении точки на расстояние ds будет Fds; отсюда следует, что бесконечно малое при­

ращение потенциала есть работа, которую совершает сила F при перемещении

единицы массы на расстояние ds.

При конечном перемещении единицы массы между некоторыми точками

М и N работа R, совершаемая силой F, будет равна разности потенциалов

в этих точках, т. е.

мм

R= 5Fdscos(F, ds)= 5dV=Vм-VN=ЛV.

(55.8)

N

N

 

239