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

Закатов Вища геодезія 1

.pdf
Скачиваний:
12
Добавлен:
28.06.2022
Размер:
23.3 Mб
Скачать

Позтому

 

dV2

4

,с /

.

 

 

 

- ^ - = r n f S (a -x ).

 

 

Дифференцируя еще раз по х, долучаєм

 

 

 

 

^

= -

4

"/б-

 

(57.12).

Так как правая часть полученного вьірашения не зависит от

х %то по ана-

логии можем написать

 

 

 

 

 

 

 

dW

 

d*V

dW

 

4

 

дх2

dy*

~dz^~

 

 

 

Суммируя все три вторьіе

производньїе, долучаєм

 

 

 

A2F2 =

-4:rt/6.

 

(57.13)'

Теперь, принимая во внимание (57.9) и (57.10), находим искомое уравне-

ниє Пуассона

dW

ow

dW

 

 

A2F

—4jt/6,

(57.14)

д х 2

д у *

 

dz%

 

 

 

 

где под б следует понимать плотность тела в малой окрестности вокруг притягиваемой точки. Если в исследуемой точке нет притягивающих масо, т. е. 6 = 0, уравнение Пуассона обращается в уравнение Лапласа. Иначе говоря, уравнение Лапласа можно рассматривать как частнмй случай уравнения Пуассона, когда 6 = 0 .

Уравнения Лапласа и Пуассона являются фундаментальними в теории потенциала.

Определение внешнего потенциала Земли и изучение ее фигурьі основьіваются на интегрировании уравнения Лапласа при дополнительннх условиях, вмтекающих из существа конкретной задачи.

Функции, непрернвньїе во всех точках данной области вместе с их первьіми и вторшми производньїми и удовлетворяющие уравнению Лапласа (57.6), назьіваются г а р м о н и ч е с к и м и ф у н к ц и я м и в зтой области.

Потенциал притяжения для точек в области, не занятой притягивающими массами, будет всегда гармонической функцией; иначе говоря, в любой точке пространства вне притягивающего тела потенциальная функция V всегда будет удовлетворять уравнению Лапласа A2F = 0.

В каждой точке внутри притягивающего тела потенциал сильї притяже­ ния будет удовлетворять уравнению Пуассона

A2F = —4rt/6.

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

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

250

§ 58. Потенциал сильї тяжести, его основньїе свойства

Напишем известное вьіражение для сили тяжести

g = F + Q .

Потенциал сили притяжения найден в § 53; для получения сили тяжести найдем потенциал центробежной сили Q. Центробежная сила при т = 1 внражается уравнением

Q = со2р.

Если во взятой системе прямоугольннх коор­ динат ось z совместить с осью вращения, то цен­ тробежная сила будет параллельна плоскости ху и составляющая зтой сили по оси z будет равна нулю. Определим составляющие ее по осям х и у. Центробежная сила направлена по радиусу р (рис. 111). Имеем:

х = р cos а

у р cos (5

 

 

 

(58.1)

 

 

 

 

 

 

Проекции центробежной сили Q на оси х и у будут равньї:

 

Qx=

Q cos а =

(о2р cos а,

(58.2)

 

Qy =

Q cos {3 =

оо2р cos р.

(58.3)

Определяя из (58.1) значення cos а и cos р и подставляя их в (58.2) и (58.3),

получаем составляющие

центробежной

сили:

 

 

 

Q x—

со2,£

і

 

 

 

Qu = <**y

*

(58.4)

 

 

Q г ~

^

 

J

 

Зти составляющие являются частньши производннми по прямоугольннм

координатам функции

 

 

 

 

 

 

 

U = ^ - ^ - + y%

(58.5)

которая, следовательно, является п о т е н ц и а л ь н о й

ф у н к ц и е й ц е н -

т р о б е ж н о й с и л и .

Действительно:

 

 

ди

=

=

Qx

 

 

дх

 

 

ди

= to2y = Qy

(58.6)

 

ду

 

 

 

 

 

251

Далее, проекция сили тяжести на оси координат будет:

£х

F

4 - 0

 

dx ‘

 

 

х '

 

 

дх

 

 

 

г

,

ri

 

dV .

dU

8»~ Fy'r

 

 

+

 

 

gz

 

: F

—1—О = = —

4-

 

 

r 2-TVz dz 1-

dz

 

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

 

 

 

 

 

 

 

 

w

= У +

V = f J ~

+

- f - + y*).

Очевидно,

 

dW _

 

<W_JL_ dV_ _

 

 

 

 

g'

 

 

dx

 

fi'r*

I

 

 

 

 

dx

 

 

 

 

 

dW

 

dV

dU

_

 

 

 

dy

 

dy

■"ТьГ -

gy

 

 

 

 

 

 

dy

 

 

 

 

dW _

dV _ ,d dU _

 

 

 

dz

 

dz

dzя*

 

 

(58.7)

(58.8)

(58.9)

Отсюда заключаем, что функция W, определяемая уравнением (58.8), является п о т е н д и а л о м с и л ь ї т я ж е с т и . П о т е н ц и а л с и л ь ї т я ж е с т и р а в е н с у м м е п о т е н ц и а л о в с и л ь ї з е м н о г о п р и т я ж е н и я и ц е н т р о б е ж н о й с и л ь ї .

Обозначая (g ,

х), (g, у), (g,z) углм, которьіе образует

направление сильї

тяжести с осями

координат,

находим:

 

 

 

 

 

gx = gcos(g,

х)

)

 

 

 

gy = gcos(g,

у)

(58.10)

 

 

g2 = gcos(g,

z)

)

 

Вьіражение для потенциала W может бьіть написано и в виде

 

w = V +

V = f § M L + -^-(*2 + уї),

(58.11)

где dx — злемент

обгема.

т

 

 

 

 

 

 

 

б — обгемная плотность.

Остановимся подробнее на свойствах потенциальной функции, изложенннх в § 55, применительно к потенциалу сили тяжести.

1. Согласно (55.6), можем написать, что

dW = gcos(g, s) = gSi

(58.12)

т. е. п р о и з в о д н а я п о т е н ц и а л а с и л ь ї т я ж е с т и

п о л ю ­

б о м у н а п р а в л е н н ю р а в н а п р о е к ц и и с и л ь ї н а

з т о н а ­

пр а в л е н и е .

2.На оснований (55.12) при cos (g-, s) = 0, т. е. при перемещении материальной точки в направлений, перпендикулярном силе тяжести, будем иметь

W — пост. = С.

(58.13)

Полученное внражение (58.13) является

общим уравнением у р о вен н ьїх

поверхностей сильї тяжести. В каждой точке такой поверхности с и л а т я ж е -

252

с т и н а п р а в л е н а п о н о р м а л и к з т о и п о в е р х н о с т и , а с о с т а в л я ю щ и е с и л ь ї т я ж е с т и п о к а с а т е л ь н о й к п о ­ в е р х н о с т и в л ю б о й т о ч н е р а в н ь ї н у л ю .

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

венную поверхность, п о в е р х н о с т ь

р а в н о в е с и я , или, как при-

нято говорить, г о р и з о н т а л ь н у ю

п о в е р х н о с т ь , слагающая

сила тяжести по касательной к которой в любой точне и в любом направлений равна нулю.

Меняя в (58.13) значение С, получаем разнне уровеннне поверхности. У р о в е н н а я п о в е р х н о с т ь , с о в п а д а ю щ а я с н е в о з - м у щ е н н о й п о в е р х н о с т ь ю о к е а н о в , н а з и в а е т с я п о -

в е р х н о с т ь ю г е о и д а .

3. Вьіражения (55.17) для потендиала силн тяжести примут вид:

(58.14)

Зти зависимости соответствуют случаю, когда в виражений (58.12) cos (g, s) = —1, т. е. когда перемещение точки происходит в направлений, противоположном направленню сили тяжести, a dh представляет собой расстояния между уровенньїми поверхностями W и W + dW.

Последнее уравнение из (58.14) показьівает, что расстояния между двумя близкими уровенннми поверхностями не равнн в разньїх точках, а обратно пропорциональньї силе тяжести, действующей в зтих точках. На полюсе, где сила тяжести имеет максимальнеє значение, уровеннне поверхности сближаются, а на зкваторе расходятся.

Из последних уравнений также следует, что dh — величина одного по­ рядна с dW и ни при каких условиях dh не может обратиться в нуль, если dW Ф 0 (так как g — конечная величина), позтому уровеннне поверхности между собой не пересекаются.

Из тех же формул следует, что уровеннне поверхности не «параллельнн» между собой.

Непараллельность уровенннх поверхностей влияет на определение висот точек земной поверхности, получаемнх из геометрического нивелирования.

Интегрируя второе из уравнений (58.14), получаем

(58.15)

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

Если зто суммирование превншений исполнено между некоторнми двумя точками физической земной поверхности А и В (рис. 112), расположенннми на конечном расстоянии и на разннх уровенннх поверхностях, то будем иметь

в

(58.16)

А

ЙЛИч,

(58.17>

g

где g — некоторое значение сили тяжести:

Нв Н А = АН — превншение точки В над точкой А.

253

Разность потенциалов W А W B постоянна. Отсюда следует, что разность внсот точек А и В будет долучать различнне значення при разном ви­ боре величини g.

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

лирования.

уровне, на котором значение

Бели принять, что точка А находитея на

потенциала сили тяжести равно W Q, то

 

 

в

 

 

W B = W 0 - ^ g d h .

(58.18)

 

о

 

Из (58.17) можно также написать

 

 

w0- w B

(58.19)

 

Нв =

 

8

 

Для

последнего принято, что НА = 0, a

g — по-прежнему некоторое

значение

сили тяжести.

 

Из (58.18) следует, что, зная значение потенциала сили тяжести в исходном пункте, можно получить из результатов нивелирования значение потен­ циала для любой точки земной поверхности. При зтом, учитнвая современную високую точность нивелирования и измерений сили тяжести, ошибка определения потенциала WB будет зависеть от ошибки, с которой известна постоянная W 0. Заметим также, что разность потенциалов W А W B, при внчислении ее по формуле (58.16), не зависит от пути нивелирования.

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

Рассмотрим две уровеннне поверхности ab и а 1 Ь1 (см. рис. 112) и допустим,

что они весьма близки одна к другой. Ускорение сили тяжести в точках а и b обозначим через g t и g 2, а отрезки нормали к уровенной поверхности ааг и bbx

— через ЛАХи АА2. Из механики известно, что материальная точка, перемещаясь из одной уровенной поверхности на другую, производит механическую работу, внражаемую произведением gh, причем зта работа не зависит от пути перемещения, а только от положення конечних точек. Следовательно, мате-

254

риальная точка, перемещаясь по нормали аах и ЬЬ1 произведет одинаковую

работу. Поатому можно написать равенство

gi Ahx= g2 Ah2.

Так как в разннх точках ускорения сили тяжести в общем

случае

раз-

личнн, то для сохранения написанного равенства должно бнть

lS.h1 Ф

Ай2,

что и

определяет непараллельность уровенньїх поверхностей.

 

V + U =

4.

Фигура Земли, определяемая в общем виде уравнением W =

= С, зависит главньїм образом от силн притяжения, так как влияние центробежной сили мало. Действительно, взяв наибольшее значение центробежной

силн, которое

она достигает

на зкваторе со2а,

вичислим отношение

 

 

 

 

q = —

,

(58.20)

где g о — значение

силн тяжести

на зкваторе.

и © — согласно (52.5), по-

Принимая

g 0

978 049

мгл,

а = 6 378 245

лучаем

 

 

 

 

 

 

 

 

 

 

Ї - - Щ

Г -

(58'21)

Следовательно, величина q одного порядка со сжатием Земли а.

При последующих виводах формул величини а и q будут приниматься малими величинами первого порядка. Если удерживать члени только с а или q и пренебрегать членами с а 2 и q2, то будем допускать ошибку порядка

("Шо)* =0.00001.

Такая точность формул даст ошибку в ускорении силн тяжести около 10 мгл. Поскольку ускорение силн тяжести внчисляют с ошибкой до сотих частей миллигала, то в соответствующих формулах необходимо удерживать величини до второго порядка малости включительно, т. е. до а 2 или q2.

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

ности, атмосфера в пределах некоторого расстояния от 3

мли будет передавать

на материальнне тела действие центробежной силн.

силн тяжести. Для зтого

5.

Найдем оператор

Лапласа для

потенциала

предварительно определим оператор Лапласа для центробежной силн. Напишем

вторне

производнне из (58.6):

 

 

 

 

d W

дЮ

8 W

_

(58.22)

 

дх2

ду* = о 2;

dz2

= 0.

 

A =

dW . dW .

д ю

D

(58.23)

 

дх2 ' Зг/2 1

dz%

= 2( 2.

 

Оператор Лапласа для потенциала силн тяжести тогда внразится фор­

мулами:

 

 

 

 

для внзшней точки

 

 

 

(58.24)

 

Д2И' = A2F + A%U = 2со2;

255

для внутренней точки

A2W = ~4л/б + 2w%.

(58.25)

 

 

Поскольку A 2U = 2ю2 постоянно, то все соображения, изложенньїе в § 57

относительно

функций А 2V,

остаются в силе и для Д2ТК.

характеризуют

Как уже

отмечалось,

вторьіе производньїе потенциала

кривизну уровенной поверхности, и там, гдеплотность меняется скачкообразно, кривизна уровенной поверхности также претерпевает скачкообразньїе изме­ нения. Зтот внвод имеет важное значение для изучения фигурьі геоида. Его иоверхность пересекает массьі разной плотности; в зтих местах кривизна по­ верхности геоида также меняется скачкообразно. Резкие изменения плотностей происходят и на самих материках в зависимости от структурьі и состава земной кори и форми земной поверхности. В атом случае они визивают скачкообразнне изменения кривизни геоида.

Кривизна геоида меняется скачком прежде всего на берегах морей и океанов, а также там,1где геоид пересекает горньїе породи разной плотности.

Вместе с тем все уровеннне поверхности и геоид, как одна из атих поверх-

ностей нигде не

терпят никаких разрнвов; ато витекает из непрернвности

потенциала

силя

тяжести.

6.

Первьіе производньїе потенциала сили тяжести так же, как и сила при-

тяжения, всюду непрерьівнн. Поскольку первне производньїе определяют направление векторов сили, т. е. силових линий, то и последние также непреривнн.

§ 59. Теорема Клеро

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

Приведем один из внводов теореми Клеро, основанний на теории потен­ циала и, в частности, на использовании вираження для потенциала сили тя­

жести в виде

ряда.

 

 

 

Напишем

на оснований (54.25)

 

 

V = Щ- + -JL- ( с

( 1 - 3 Sin2 0) + J i L ( B - A ) cos2 Ф cos 21.

(59.1)

Так как для аллипсоида вращения моменти инерции в плоскости акватора

между собой равнн, т. е. А

= В, то (59.1) примет вид

 

 

V = JT - + 4 > ( C - A) <~1 - З а іп ‘ Ф).

(59.2)

Разность С — А имеет размерность массьі, помноженной на квадрат рас-

стояния.

 

 

 

 

Положим

 

С — А==ца2,

(59.3)

 

 

где а — большая полуось

земного

зллипсоида,

вдоль

ц — некоторая воображаемая

добавочная масса, расположенная

зкватора, обусловливающая различие С от А и В.

 

256

Тогда

 

 

 

 

l/ = i ^ - + - ^ - ( l - 3 s i n >Ф).

(59.4)

Теперь получим внражение для потенциала сильї тяжести W.

 

Так как W =^= V +

Q,

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

 

W

+

( 1 - 3 sin2 Ф) + -^ -г2 cos2 Ф.

(59.5)

Преобразуем внражение для потенциала центробежной силн, т. е, последний член формульї (59.5).

Возьмем отношение центробежной силн к силе тяжести на екваторе, т. е.

q = ^ ~ .

(59.6)

Так как q — малая величина первого порядка, можем заменить g Q черев

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

со2а3

(59.7)

я = • f M *

 

откуда

 

CD2-

f M q

вз *

Тогда потенциал Q примет ВИД

< ? = — г COS2 Ф = ~ ^ - q r 2 co s2 0 *

Подставляя (59.9) в (59.5), находим

II І*-** Ч* -----і і

+ гої

Г2

v

3

sin2 Ф) -f-

соз2Ф^|.

 

 

К -

(1

(59.8)

(59.9)

(59.10)

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

Зто уравнение имеет вид

 

 

 

 

 

W = C0.

 

 

 

г- Для определения

искомой

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

поставим

условие,

чтобн на

екваторе, т. е. при Ф =

0°, радиус поверхности равнялся большой полуоси а.

Подставив зти значення

для

г

и Ф в (59.10),

получим внражение

для С0

 

 

“ Г

( 1+

2м + 1 ) -

Мо­

 

(59.11)

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

 

 

 

 

 

 

 

 

 

 

 

 

Г - [ 1+ W

-

3 sin2 ф) +

1

cosSф]

= - т - ( 1

2М

(59.12)

 

>рли

 

 

 

 

 

 

 

 

 

,

 

i + w

^

- (f-

3si”s<p) + l 5 - Ms20

 

(59.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2М

 

 

 

п , С, Закатов

257

Сохраняя лишь мальїе величиньї порядка

и

q,

последнее

внражение

можем переписать

 

 

 

 

 

г

1 + '2 ¥ ( 1 “ 3sin2a>)+ T COs20~

^

Ч

 

а

2

 

или

 

 

 

 

 

 

1 “ ( 4 і г + у ) 8 і п 3 ф -

 

 

<5 9 Л 4 >

Но на оснований (4.16) с удержанием членов порядка е2 и а

для зллип-

соида вращения

 

 

 

 

 

г

У і —е2

Ф = 1— a sin2 Ф.

(59.15)

а

= 1 — 4гsin2

У і — Є2'cos2 Ф

 

 

 

 

Сравнивая (59.14) и (59.15), приходим к заключению, что уравнение (59.14)

является з л л и п с о и д о м в р а щ е н и я

с о с ж а т и е м

 

З р

(59.16)

а = - 2 ¥

 

Зтот внвод получен на основе теории потенциала сильї тяжести на Земле; в процессе вьівода принята только одна гипотеза, что Земля — тело вращения, по форме незначительно отступающее от шара.

Существенно заметить, что полученннй вьівод справедлив только в том случае, если все массн находятся внутри геоида, так как только при зтом условии можно разложить потенциал сильї тяжести в ряд. Зто понятно, так как ми уже установили, что в пересечении геоидом материков, когда притягивающие массн оказнваются вне геоида, его кривизна меняется скачком, а поверхность геоида под материками внражается другой аналитической функцией, чем на море. Отступления от допущенного условия распределения масс в теле Земли но долготе также вьізьівают соответствующие отклонения фигурн геоида от поверхности зллипсоида вращения. Зти замечания должнн приниматься во внимание при детальном изучении фигурьі Земли; в то же время влияние зтих отступлений от реальних условий вьіражается величинами, хотя далеко не пренебрегаемнми, но не могущими изменить сделанньїй внвод о том, что фигура геоида близка к фигуре зллипсоида вращения.

Зллипсоид вращения, поверхность которого внражается уравнением (59.14),

иногда назьівают « и д е а л ь н н м

г е о и д о м » .

Теперь перейдем к получению

вираження для величини сильї тяжести

на уровенной поверхности, определяемой уравнением (59.14), т. е. на «идеальном геоиде». Из (58.14) имеем

dW

(59.17)

dk

 

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

В (59.17) значение силн тяжести внражается как производная от потенциала по нормали «к идеальному геоиду» с точностью до малих величин первого порядка, имеющему форму зллипсоида вращения. Однако в (59.5), внражающем потенциал силн тяжести, направление нормали не входит; определим, к а к у ю

-

dW

производную

dW

ошибку мн допустим, взяв вместо

 

 

258

Различие направлений нормали h и радиуса-вектора гвнражается разностью геодезической и геоцентрической широт, которая максимально достигает значення 11,8'. На оснований (58.12) можно сказать, что, дифференцируя по г, а не по h, мн вместо g получим g cos (ф — Ф), т. е.

g cos ( ф - Ф ) = -

где ф — широта, образованная н о р м а л ь ю к поверхности земного зллипсоида.

Так как cos (ф — Ф) =

cos 11,8' = 0,999995, то с точностью 5-10“ 6 можем

написать

 

dW

 

 

 

 

 

£~~

дг

 

Дифференцируя (59.10)

по г, получаем

 

g = -71г +

(1 -

з sin* Ф ) - fMq — cos2 Ф.

(59.18)

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

Чтобьі получить g для точек рассматриваемой поверхности, заменим г в главном члене его вьіражением согласно (59.14), а в остальншх членах, вследствие их малости, положим г = а, тогда

 

sin2 ФТ ['1 +4 Ж ^- 3зіп2 ф)- qcosa ф] •

<59л9>

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

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

 

 

*=4Н 1+Ч 4 £ + т )зіп2ф] [ 1+4тг<1- 3з“2ф)-«соз2ф]

<59-2°)

или

 

 

 

 

г = - ^

[ і + г - ж - * +

( 2* - ? ж ) * і»г ф ] -

<59-21>

Из (59.16) имеем

 

 

 

 

После подстановки

в (59.2) а ---- 1

вместо ~ ^ получим

 

г = ^ [ 1 + а - т « + ( т ® - “ ) зіп 2ф ] -

(59.22)

 

Обозначая через g 0 ту часть (59.22), которая не зависит от широти, т. е.

Волучаем

 

УЇГ ( і + “’— § •« )= в»

(59.23)

 

 

 

 

 

 

 

 

-jq-CC

 

 

8 =

1

sin2 Ф 1.

(59.24)

17*

 

1+ <Х|” т 9

 

 

 

 

259

Соседние файлы в предмете Вища геодезія