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

Основы релятивистской механики. Монография

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
841 Кб
Скачать

Из определения (П2.1) функции Tn+1 находим

 

Tn+1

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2a

+

q

+

2a + +

q

+ ,

 

 

 

 

 

 

 

qn+1

 

 

k,n 1

k

 

 

n 1,n 1

n 1

 

 

 

 

 

k=1

 

 

 

 

 

 

 

 

 

 

 

 

Tn+1

 

 

 

n

 

 

 

 

 

 

 

 

 

q +

 

= 2q +

a

+

q

+ 2a +

+

q

2 + .

(П2.5)

q +

 

n 1

 

n

1

k=1

k,n 1

k

n

1,n 1

 

n 1

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

Сложение левых и правых частей равенств (П2.4), (П2.5) дает

n+1

Tn+1

 

n

 

qk

= 2 Tn

+ 2qn+1 ak,n+1qk + an+1,n+1q2n+1

.

q

k=1

 

k=1

 

k

Отсюда, учитывая формулу (П2.3) приходим к искомому выражению

n+1

Tn+1

 

q

= 2T + ,

q

k

n 1

k=1

k

 

которое и доказывает утверждение леммы Эйлера.

131

ПРИЛОЖЕНИЕ 3.

ОСТАТЬЕ А. ЭЙНШТЕЙНА 1915 г. «ОБЪЯСНЕНИЕ ДВИЖЕНИЯ

ПЕРИГЕЛИЯ МЕРКУРИЯ В ОБЩЕЙ ТЕОРИИ ОТНОСИТЕЛЬНОСТИ»

Параграф 2 «Движение планет» начинается так: «Уравнения движения материальной точки в поле тяже-

сти, вытекающие из общей теории относительности, имеют вид

d2xν

= Гν

dxσ

 

dxτ

,

(7)

 

 

 

ds2

στ

στ

ds ds

 

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

................................................................................

Тогда первые три уравнения будут в точности ньютоновскими. Если ввести в плоскости орбиты полярные координаты r, ϕ , то сохранение энергии и момента количества движе-

ния дает известные уравнения

 

 

 

1

u2

 

 

 

 

 

 

 

+ Φ = A

 

2

 

 

 

 

 

 

 

,

(8)

 

 

 

2

 

 

r

= B

 

 

 

 

 

ds

 

 

 

 

 

 

 

 

 

 

 

где А и В – постоянные законов сохранения энергии и момента количества движения; здесь, кроме того, использованы сокращенные обозначения

Φ = − α ,

2r

u2 = dr2 + r22 .

(8a)»

ds2

(конец цитаты).

132

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

s = ct, α =

2GM

,

(8б)

c2

 

 

 

где с – скорость света в вакууме. Тогда формулы (8), (8а) записываются так

1

v2

 

 

 

 

 

+ Φ0 = E

 

2

 

 

 

 

 

,

(8)*

 

 

 

 

r

2

= l

 

 

 

dt

 

 

 

 

 

 

 

 

Φ0 = − GM , r

v2 = dr2 + r22 .

(8а)*

dt2

Здесь постоянные E, l связаны с постоянными А, В равенствами

E = c2 A , l = cB.

(8в)

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

2

x

α

x

 

1 + α + 2u2

 

2

 

 

 

d

 

 

 

 

 

 

ν

= −

 

ν

 

− 3

dr

 

 

(ν = 1,2,3) .

(7б)

 

 

 

 

 

 

 

2

 

2 r

3

 

r

ds

 

 

 

ds

 

 

 

 

 

Оговоримся здесь, что формула (7б) взята из оригинального немецкого издания работы Эйнштейна 1915 г. [22]. Русский перевод этой работы, представленный в академическом четырехтомнике 1965 г. [16], дает эту ключевую формулу с грубой опечаткой (!?), полностью исказившей ее содержание. Вот как она там напечатана:

133

d2x

αx

 

 

3 dr 2

 

ν

= −

 

ν

u2

 

 

 

 

 

,

2

r

3

 

 

ds

 

 

 

2 ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и без нумерации.

После вывода уравнений движения в форме (7б) Эйнштейн с помощью равенств (8) и замен переменных преобразует эти уравнения к виду

 

d2x

 

∂Φ

 

 

 

 

 

ν

= −

x

,

 

 

 

 

 

 

 

 

ds2

 

 

 

 

 

 

 

 

 

 

ν

 

 

(7в)

 

 

α 1

 

 

B2

 

 

 

 

 

Φ = −

 

 

1 +

 

 

.

 

2 r

r

2

 

 

 

 

 

 

 

 

Оговоримся еще раз, что в последней записи потенциала Ф допущена опечатка как в оригинальном немецком издании, так и в перепечатанном русском переводе. Именно пропущен множитель 1/r. Очевидно, что потенциал не может

быть таким:

α

 

B2

 

+

 

1

 

 

.

2

r

2

 

 

 

 

 

Уравнения (7б) вполне определены и разрешимы без привлечения классических законов сохранения (8), (8а). Причем они допускают вывод не только угла смещения перигелия, но и вывод явного уравнения релятивистской орбиты (т.е. явной зависимости полярного радиуса от угловой координаты), чего в статье Эйнштейна не сделано. Покажем, как осуществить эту процедуру.

Прежде всего отмечаем, что из уравнений (7б) следует центральность вектора ускорения (коллинеарность его ради- ус-вектору). Поэтому имеет место тот же закон сохранения количества движения (вторая формула (8)). Соответственно из представления квадрата скорости в полярных координатах (вторая формула (8а)) имеем

dr 2

2

 

2

dϕ 2

2

 

B2

 

 

 

 

 

= u

r

 

 

 

 

= u

 

.

(8г)

 

 

 

r2

ds

 

 

 

ds

 

 

 

 

134

С учетом этого равенства уравнения (7б) принимают вид

d2xν

 

α xν

 

 

α

 

2

 

B2

 

 

= −

 

 

 

1

+

 

u

 

+ 3

 

 

(ν = 1,2,3) .

(7г)

2

2 r

3

r

 

r

2

ds

 

 

 

 

 

 

 

 

 

 

Умножаем обе части последнего уравнения на компоненту

скорости uν = dxν и суммируем его по индексу . Получаем ds

 

du2

 

 

 

 

α dr2

 

1

 

 

 

 

 

 

α

 

 

 

2

 

 

B2

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

 

u +

3

 

 

,

 

 

 

 

 

 

2 ds

 

r

3

 

r

r

2

 

 

ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du2

 

 

 

 

 

α

 

 

 

 

α

 

 

2

 

 

B2

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

1 +

 

 

 

u

 

 

+

3

 

 

,

 

 

 

 

 

dr

r

2

 

r

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du2

α

2

 

 

 

 

 

 

α

 

 

 

 

α

 

B2

 

 

 

 

 

 

 

 

 

 

u

 

 

= −

 

 

 

 

1 +

 

+ 3

 

 

.

(9)

 

 

 

 

 

 

 

r

2

 

 

r

2

 

r

r

2

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это линейное уравнение относительно квадрата скоро-

сти. Видим, что в нем удобно ввести замену

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ =

1

,

 

 

= −

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

Теперь уравнение (9) приобретает вид

 

 

 

 

 

 

 

 

du2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ αu2 = α + α2λ + 3αB2λ2 .

(9а)

 

 

 

 

Его точное решение записывается так:

 

 

 

 

 

 

 

6B2

 

 

 

 

 

 

 

 

 

6B2

 

 

 

 

 

 

 

 

 

 

 

u2 =

 

 

 

 

 

+ α 1 −

 

 

 

 

 

λ + 3B2λ2 + Ce−αλ ,

(10)

 

 

α

2

 

 

α

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где С – постоянная интегрирования.

 

 

 

 

 

 

 

Для учета

 

только

первой

 

релятивистской

поправки

к классическому закону

 

 

сохранения

энергии

нам здесь

в экспоненте необходимо удержать только члены до порядка α3λ 3 включительно. Тогда выражение (10) запишется так:

135

 

 

 

 

 

 

 

 

 

 

u2 = C +

6B2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α2

 

 

 

 

 

 

 

 

 

 

 

 

(10а)

 

 

 

 

 

 

 

 

2

λ + α

2

 

 

 

 

 

 

 

2

 

 

 

 

α

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ α 1

 

C

+

6B

 

 

 

 

C

+

6B

 

 

λ2

Cλ3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

2

 

 

2

 

 

 

 

α

2

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

закону сохранения энергии:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 A = C +

6B2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α2

 

 

 

 

 

 

 

 

 

 

 

 

Теперь (10а) приобретает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2

 

 

 

α2

2

 

α

 

 

 

 

 

 

 

 

αB2

 

Aα2

3

 

 

 

 

= A

1

+

 

λ

 

+

(1 − 2 A)λ +

 

 

 

 

1 −

 

 

 

λ

.

 

2

2

 

 

2

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

3B

 

 

 

 

Учитывая, что α

1 2 , A 1

2 , B2

1

2

 

вторые слага-

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

c

 

 

 

 

 

c

 

 

 

1

 

емые в первой и третьей скобках имеют порядок

. По-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c4

 

этому без потери членов второго порядка малости они могут быть опущены. Тогда имеем

 

 

u2

 

= A + α (1 − 2 A)λ + αB2

λ3 .

 

 

(10б)

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

С другой стороны, с учетом замены

r = 1

λ

из формулы

(8г) следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

2

 

2

 

 

 

 

 

 

 

 

 

 

u

 

= B

 

 

 

 

+ λ

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Внося это выражение квадрата скорости в (10б), получа-

ем уравнение траектории

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

2 A

 

 

1 − 2 A

 

3

 

 

 

 

 

 

+ λ

 

=

 

 

 

+ α

 

 

 

 

λ + αλ

.

(11)

 

 

B2

 

 

B2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

136

В статье Эйнштейна это уравнение записано так:

dx 2

 

2 A

 

α

 

2

3

 

 

 

 

 

=

 

+

 

x x

 

+ αx

.

(11)

 

B2

B2

 

 

 

 

 

 

 

 

Здесь использована другая буква: x = λ = 1r . И говорится:

«Это уравнение отличается от соответствующего уравнения теории Ньютона только последним членом в правой части».

Видим, что эйнштейновское уравнение траектории (11) отличается от выведенного нами из эйнштейновского же уравнения движения (7б) видом константы при первой сте-

пени переменной x = λ = 1r . Однако путем переобозначения

постоянных они переходят одно в другое.

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

 

d2λ

+ λ = α

1 − 2 A

+

λ2 .

(12)

 

2

 

 

 

 

 

 

 

 

 

2 B2

2

 

 

 

Введем две новые константы:

 

 

 

 

 

 

 

 

1

 

= α

1 − 2 A

,

ε =

.

(13)

 

 

 

 

 

 

 

 

 

h

2

 

B2

 

 

2h

 

Теперь уравнение (12) предстает в виде

 

 

 

 

 

 

d2λ

 

+ λ =

1

+ εhλ2 .

(12а)

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

Постоянная h выполняет роль фокального параметра классической эллиптической орбиты. Вводя еще переменную ξ = hλ , записываем уравнение (12а) в безразмерной форме:

d2ξ

2

 

 

+ ξ = 1 + εξ .

(12б)

2

 

 

Ищем решение в виде разложения по малому параметру:

ξ = ξ0 + εξ1 + O2 ).

137

Внося это представление в уравнение (12б) и приравнивая члены при одинаковых степенях , найдем

 

 

e

2

 

e

2

 

+ O2 ) , (13)

ξ = 1 + e cos ϕ + ε 1

+

 

+ eϕ sin ϕ −

 

cos2ϕ

 

 

 

 

 

2

6

 

 

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

Точки афелия и перигелия определяются

= 0. Имеем dϕ

 

e

 

 

= −e sin ϕ + εe sin ϕ + ϕ cos ϕ +

 

sin 2ϕ .

 

 

3

 

условием

(14)

Пусть первый перигелий приходится на точку ϕ = 0 . Ищем следующий перигелий через оборот в 2 радиан смещенным на угол порядка :

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

= −e sin δ + εe sin δ + (2π + δ) cos δ +

 

sin 2δ .

 

3

 

 

2π+δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Удерживая здесь лишь первые порядки по и , получим

 

 

 

 

 

 

= −eδ + 2πεe = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

2π+δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ = 2πε,

 

 

(15)

или, вспоминая значения параметров (6б), (13):

 

 

 

 

 

 

δ =

3πα

=

GM

.

 

 

(15а)

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

c2h

 

 

 

Фокальный параметр выражается через большую полуось и эксцентриситет эллипса так: h = a(1 − e2 ). Поэтому формулу (15а) обычно записывают в виде

δ =

GM

 

.

(15б)

2

e

2

)

 

c a(1

 

 

 

138

Для наглядности выпишем явное выражение зависимости полярного радиуса от угловой координаты (т.е. уравнение орбиты). Переписываем функцию (13) так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

2

 

 

 

e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

 

+ eϕ sin ϕ −

 

cos2ϕ

 

(ε

 

).

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

6

 

 

 

 

 

2

ξ = r

= (1

+ e cos ϕ)

1

+ ε

 

 

 

 

 

 

 

 

 

+ O

 

 

 

 

 

 

1 + e cos ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

2

 

 

 

 

e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

 

 

+ eϕ sin ϕ −

 

 

cos2ϕ

+ O (ε

 

) ,

 

 

r

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

1 − ε

2

 

 

 

6

 

 

 

 

 

2

 

 

h

1 + e cos ϕ

 

 

 

 

 

 

 

1 + e cos ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

3GM

1 +

e2

+ eϕ sin ϕ −

e2

cos2ϕ

 

1

 

 

 

 

 

 

 

 

r =

 

 

 

 

2

6

 

+ O

 

.

(16)

 

+ e cos ϕ

c2

 

(1 + e cos ϕ)2

 

1

 

 

 

 

c4

 

 

Посмотрим теперь на эйнштейновские уравнения движения (7б) с точки зрения описания прямолинейного движения. Другими словами, рассмотрим случай падения материальной точки на точечную массу М, создающую гра-

витационное поле. В этом случае u = dr и мы вместо трех ds

уравнений (7б) имеем только одно уравнение движения. А именно:

 

d2r

 

α 1

 

 

 

α

2

 

 

 

 

 

= −

 

 

 

 

 

1

+

 

u

,

2

2 r

2

r

 

ds

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du

 

α 1

 

 

 

α

2

 

 

 

 

 

= −

 

 

 

 

 

1

+

 

u

.

 

 

ds

2 r2

r

 

 

 

 

 

 

 

 

139

Так как

du = du dr = u du = 1 du2 , ds dr ds dr 2 dr

то мы приходим к линейному уравнению относительно квадрата скорости, которое в точности совпадает с уравнением (9), если в последним положить В = 0. То есть

du2

α

u2 = −

α

 

+

α

 

 

 

 

 

1

.

dr

r

2

r

2

 

 

 

 

r

Его решение, очевидно, записывается так:

u2 =

α

α

 

+ Ce r .

r

Положим s = ct и перейдем от безразмерного времени к обычному:

dr 2

1 dr 2

 

 

2GM 1

 

 

=

 

 

 

 

 

=

 

 

 

 

 

 

 

2

 

 

 

c

2

 

r

ds

с

dt

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

2GM

 

 

 

 

 

 

 

 

= −

 

 

1

+

 

 

 

 

 

 

 

Are

 

 

 

 

 

 

 

 

 

dt

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2GM 1

+Ce c2 r .

2GM 1

c2 r .

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

по степеням 1 и строить различные приближения к точной r

зависимости r(t) . Но вот при малых r ( r 0 ), как следует из последней формулы, релятивистский член уничтожается. Однако в этой области сильного гравитационного поля исходные уравнения (7б), по-видимому, не применимы вообще. Какими тогда должны быть исходные уравнения движения вместо уравнений (7б), пригодные вблизи точечной массы, создающей гравитационное поле? Каким должен быть тензор Гστν в уравнениях движения материальной точки (7), вытекающих из общей теории относительности, чтобы их можно было считать точными во всем пространстве

140

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