прибору данной конструкции. Точность определения относительного ускоре
ния силы тяжести в этом случае зависит только от случайных ошибок в опреде
лении периода качания маятников и изменений систематических ошибок
и постоянных прибора, которые могут произойти после наблюдений на исход
ном пункте.
Применением соответствующей методики и конструкции прибора ослаб
ляется влияние и этих ошибок.
Точность окончательного вывода периода качания маятника после учета всех поправок характеризуется ошибкой около ±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 |
равно |
|
|
|
|
1·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 |
dт • |
|
Тогда будем иметь |
|
|
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 - |
у2- z2) dm. |
(5/2.18) |
Аналогично две другие |
части представятся |
интегралами: |
|
|
~~ S(2у2-х2- 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