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

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

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

 

 

x1

 

 

 

x1

...

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

q

 

 

 

 

 

 

 

 

q

 

 

 

1

2

 

 

 

 

n

 

 

 

 

x2

 

 

x2

...

 

x2

 

 

 

 

 

 

 

 

qn

 

 

J =

q1

 

 

 

 

 

q2

 

,

(1.7)

 

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

 

 

 

 

xn

 

xn

...

xn

 

 

 

 

 

 

q

 

 

 

 

q

 

 

 

 

 

q

 

 

 

 

1

2

 

 

 

n

 

 

известную как матрица Якоби. Невырожденность этой матрицы обеспечивает взаимную однозначность преобразований (1.3), (1.4). При этом определитель матрицы Якоби служит коэффициентом искажения элементарных объемов, выраженных через координаты {xi} и обобщенные координаты {qi}, т.е.

dx1 dx2 ... dxn = det J dq1 dq2 ...dqn.

(1.8)

Преобразования (1.5) позволяют записать такие тождества:

xi

=

xi

.

(1.9)

qj

 

 

qj

 

Нетрудно также убедиться, что полная производная по времени может быть внесена под знак частной производной по координате. Действительно:

 

 

d

x

 

 

 

n

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

=

 

 

 

 

 

 

 

 

i

q

 

=

 

 

 

 

 

dt q

 

 

q

 

 

q

 

 

 

 

 

 

 

j

 

k=1

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

j

 

 

 

 

 

 

n

x

 

 

 

 

 

 

 

 

 

n

x

 

 

 

 

x

=

 

 

 

 

i

q

=

 

 

 

 

 

 

i

 

q

=

i

.

 

 

 

 

 

 

 

 

 

qk

 

 

k=1

qj

qk k

 

 

 

qj k=1

 

 

k

 

qj

Итак:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

xi

 

=

xi

.

 

 

 

 

 

 

(1.10)

 

 

 

 

 

 

 

 

dt qj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qj

 

 

 

 

 

 

 

 

Два тождества – (1.9) и (1.10) помогут нам составить еще одну важную систему тождеств с участием вторых производ-

11

ных от координат по времени, т.е. с участием компонент ускорения. А именно рассмотрим полную производную

по времени от произведения x

 

xi

. Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i qj

 

 

 

 

 

 

 

 

 

 

 

 

 

d

x

 

 

 

 

x

 

 

 

 

d

 

 

x

 

 

 

 

 

x

 

 

i

 

=

 

 

 

i

 

x

+ x

 

 

 

 

 

i

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

i

dt qj

 

 

dt

qj

 

 

qj

 

 

 

 

 

Вносим сюда равенства (1.9), (1.10):

 

 

 

 

 

 

 

 

 

d

 

x

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

i

 

 

=

 

 

 

 

i

 

x

+ x

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

qj

 

 

 

 

 

qj

 

i

 

 

i

qj

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

x2i

 

 

 

 

x2i

=

xi

x .

(1.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt qj 2 ∂qj 2 ∂qj

 

 

i

 

 

 

 

 

 

 

 

 

 

Напомним, что здесь индексы i и j пробегают значения от 1 до n. Суммирование равенств (1.11) по индексу i дает

 

 

=

i

x ,

d

 

T

 

T

n

x

 

dt qj

 

qj

i=1 qj

i

 

 

где

1 n

T = x2i .

2 i=1

(1.12)

(1.13)

Заимствуя термин из механики, назовем функцию Т кинетической энергией.

Чтобы кинетическую энергию выразить в обобщенных координатах и скоростях, нужно равенства (1.5) подставить в выражение (1.13). Тогда получим

 

 

1

 

n

n

 

 

 

 

 

 

 

T =

∑∑a q

q

 

,

 

 

j

 

 

2 i=1 j=1

ij

i

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

x

 

x

 

 

a

= a

ji

=

k

 

 

 

k

.

 

 

 

 

 

ij

 

 

 

 

qi

 

 

qj

 

 

 

 

 

 

k=1

 

 

(1.14)

(1.15)

12

Таким образом, равенства (1.12) представляют собой выражения компонент вектора ускорения, преобразованные транспонированной матрицей Якоби.

Возьмем теперь произвольную дважды дифференцируемую функцию координат и скоростей:

Φ = Φ(x1,x2 ,...xn,x1,x2 ,...xn ),

(1.16)

образуем из нее вектор с компонентами

f =

d

∂Φ

∂Φ

(1.17)

 

i

dt xi

xi

 

 

 

и посмотрим, как преобразуются его компоненты при преобразовании координат согласно формулам (1.4). Имеем

∂Φ = ∂Φ xi

+ ∂Φ xi ,

 

 

n

 

 

 

n

 

 

 

 

q

 

 

 

 

 

 

x

 

 

 

j

i=1

x q

j

i=1

q

j

 

i

 

i

 

∂Φ = ∂Φ x

,

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

q

 

 

x

 

 

 

 

 

 

 

j

i=1

q

j

 

 

 

 

 

 

i

 

 

 

 

 

 

или с учетом (1.9), (1.10)

 

 

∂Φ = ∂Φ xi + ∂Φ d ∂xi ,

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

j

i=1

x q

j

i=1

dt q

j

 

 

 

 

 

 

i

 

 

i

 

 

 

 

 

 

 

∂Φ = ∂Φ xi .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

i=1

q

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

∂Φ

∂Φ

 

n

x

 

 

d

∂Φ

∂Φ

 

 

 

 

=

 

i

 

 

 

.

 

dt qj

 

 

 

 

 

 

 

qj

 

i=1

qj dt xi

 

 

xi

И соответственно, обратные преобразования:

d ∂Φ

 

∂Φ

n

qj

 

d

∂Φ

 

∂Φ

 

=

 

.

 

 

 

 

 

 

 

 

 

dt x

x

 

 

 

 

q

 

 

j=1

x dt q

j

 

j

 

 

i

 

i

i

 

 

 

 

 

(1.18)

(1.19)

13

1.2. Уравнения движения

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

xi = fi (x1,x2 ,...,xn,x1,x2 ,...,xn ),

(1.20)

i = 1, 2, …, n.

 

В механике такая ситуация интерпретируется как движение системы с n степенями свободы в силовом поле с компонентами

{f1,f2,...,fn }.

(1.21)

При этом весьма часто встречаются системы, в которых компоненты силового поля представимы в виде формул (1.17). Тогда функция Φ(x1,x2 ,...xn,x1,x2 ,...xn ) называется потенциалом силового поля {f1,f2,...,fn }. Соответственно си-

стема (1.20) приобретает вид

 

 

 

 

x

=

d

∂Φ

∂Φ .

(1.22)

 

 

 

 

 

i

 

dt xi

xi

 

 

 

 

 

 

 

Записываем координаты вектора ускорения через функ-

цию кинетической энергии (1.13):

 

 

 

d ∂T

=

d

∂Φ

∂Φ .

(1.23)

 

 

 

 

 

 

 

 

dt x

 

 

 

dt x

x

 

 

 

 

 

i

 

 

 

 

i

i

 

Теперь, вводя функцию

 

 

 

 

 

 

 

 

 

 

 

L = T – Ф,

 

(1.24)

известную как функция Лагранжа, приводим уравнения движения к следующей стандартной форме записи:

d ∂L

L

= 0.

(1.25)

 

 

 

 

dt x

 

 

x

 

 

 

i

 

i

 

Благодаря формулам преобразований (1.18) эта система

вобобщенных координатах qj(j = 1, 2, …, n) предстает точно

втаком же виде:

14

d

 

L

L

= 0.

(1.26)

 

 

 

dt qj

qj

 

Итак, мы пришли к инвариантной (неизменной по отношению к любым координатным преобразованиям) форме записи системы из n дифференциальных уравнений второго порядка. Эти уравнения, т.е. уравнения, записанные в форме (1.26), называются уравнениями движения Лагранжа.

Система уравнений Лагранжа (1.26) допускает сразу, в общем случае, вывод одного из ее первых интегралов. Записываем полную производную функции Лагранжа по времени:

 

 

 

 

dL

 

n

 

L

 

 

 

n

L

 

 

 

 

 

 

 

 

=

q

 

 

+

q .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

j=1

qj

j

j=1

qj

j

 

 

В силу системы (1.26)

 

 

 

 

 

 

 

 

 

 

 

 

 

dL

n

d ∂L

 

 

n

 

 

L

 

 

n

 

L

 

=

q

 

+

 

q

=

d

 

 

 

 

 

 

 

 

j

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

qj

 

j=1 dt qj

 

 

j=1 qj

j=1 dt

или

d

 

n

 

 

 

 

L

q

 

L

= 0.

 

 

j

 

 

 

 

 

 

 

dt j=1 qj

 

 

 

q ,

j

(1.27)

Отсюда, очевидно, следует

n

L

 

 

 

 

 

 

 

 

q

 

L = E = const.

(1.28)

 

j

j=1

qj

 

 

 

 

 

 

Функция

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

H =

L

q

 

L

(1.29)

 

 

j

 

 

 

 

j=1 qj

 

 

называется функцией Гамильтона.

Внесем в (1.29) выражение (1.24), определяющее функцию Лагранжа:

H = ∂(T − Φ) q

 

− (T − Φ) = T q

 

T + Φ − ∂Φ q .

n

 

n

 

n

 

 

j

 

 

 

j

j=1 qj

j

j=1

qj

 

j=1

qj

 

15

Исходя из представления кинетической энергии (1.14), можно показать, что имеет место равенство

n

T

 

 

 

 

q

 

= 2T,

(1.30)

 

j

j=1

qj

 

 

выражающее так называемую лемму Эйлера (см. вывод, приложение 2). Поэтому функция Гамильтона может быть записана так:

H = T + Φ − ∂Φ q

 

.

(1.31)

n

 

 

 

j=1 qj

j

 

 

Тогда первый интеграл (1.28), который принято называть интегралом энергии, приобретает вид

T + Φ − ∂Φ q

 

= E.

(1.32)

n

 

 

 

j=1 qj

j

 

 

Системы, в которых потенциал Ф не зависит от скоростей, называются консервативными. В этом случае равенство (1.32) уменьшается на одно слагаемое, и оно становится таким:

T + Φ = E.

(1.33)

В физике это соотношение известно как закон сохранения энергии.

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

Этот переход осуществляется следующим образом. Определим с помощью равенств

p

 

=

L

, j = 1,2,...,n

(1.34)

j

 

 

 

qj

 

n новых неизвестных функций времени – так называемых обобщенных импульсов. Тогда, согласно определению (1.29), функция Гамильтона записывается так:

16

 

n

 

 

H = pjqj L.

 

j=1

 

 

Отсюда

 

 

 

H = qj

, H = −

L

.

 

pj

qj

qj

Но из системы уравнений Лагранжа следует, что

L = d ∂L = pj ,

qj dt qj

и мы приходим к системе, известной как система уравнений движения Гамильтона:

qj

=

H ,

pj

= −

H .

(1.35)

 

 

pj

 

 

qj

 

1.3.Движение в трехмерном пространстве

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

q1 = r, q2 = , q3 = .

(1.36)

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

x = rsin cos , y = rsin cos , z = rcos .

Записываем три уравнения Лагранжа:

d

L

L = 0,

 

dt r

r

d

L

L = 0,

 

dt ∂ϕ

∂ϕ

d

L

L = 0.

 

dt ∂θ

∂θ

(1.37)

(1.38)

(1.39)

(1.40)

Рассмотрим сначала наиболее очевидный и часто встречающийся случай, когда в функции Лагранжа L=Т – Ф кинетическая энергия пропорциональна квадрату скорости,

17

а потенциальная энергия не зависит от угловых координат и угловых скоростей. Такой потенциал называют радиальным. Итак, пусть

 

 

 

 

 

 

 

 

 

 

 

 

T =

v2

, Φ = Φ(r,r ),

 

 

 

 

 

 

 

 

 

(1.41)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

2

= x

2

+ y

2

 

2

= r

2

+ r

2 2

+ r

2

sin

2

 

2

.

(1.42)

 

 

 

 

 

 

 

+ z

 

 

 

θ

 

 

 

 

θ ϕ

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d T

 

 

T

 

 

 

 

2

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= r rθ

 

r sin

 

 

θ ϕ

 

,

 

 

(1.43)

 

 

 

 

 

 

dt r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d T

 

 

T

 

 

 

2

 

 

 

2

 

θ ϕ + 2rr sin

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

∂ϕ

 

= r

 

sin

 

 

 

 

θ ϕr

 

 

sin2θ θϕ,

(1.44)

 

dt

∂ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d T

 

T

 

 

 

2

 

 

 

 

r2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂θ

= r

 

θ +

2rrθ −

 

 

 

sin2θ ϕ

,

 

(1.45)

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

dt ∂θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и соответственно уравнения (1.38)–(1.40) в данном случае предстают в виде

r rθ2 r sin2 θ ϕ2 = d ∂Φ ∂Φ , dt r r

r2 sin2 θ ϕ + 2rr sin2 θ ϕr2 sin2θ θϕ = 0,

r2θ + 2rrθ − r2 sin2θ ϕ2 = 0.

2

(1.46)

(1.47)

(1.48)

Видим, что эти уравнения могут описывать два характерных движения – движение по прямой, соответствующее нулевому значению угловой долготы = 0, и плоское движе-

ние, соответствующее значению = π . В первом случае оста-

ется только одно уравнение:

 

2

 

 

 

 

r = z =

d

∂Φ

∂Φ .

(1.49)

 

 

dt z

 

z

 

18

Во втором случае имеем два уравнения:

2

 

d ∂Φ

 

∂Φ

 

 

r rϕ

=

 

r

r

,

(1.50)

dt

 

d

(r2ϕ) = 0.

 

 

(1.51)

 

 

 

 

 

dt

 

 

 

 

 

 

 

На них и остановимся.

Из уравнения (1.51) получаем еще один первый интеграл

r2ϕ = l = const,

(1.52)

выражающий закон сохранения момента вектора скорости. Внося (1.52) в уравнение (1.50), получим уравнение,

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

= π :

 

 

 

 

 

 

 

2

l2

 

d ∂Φ

 

∂Φ

 

r

 

=

 

 

.

(1.53)

r3

dt r

 

 

 

r

 

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

Первый интеграл (1.52) предоставляет нам правило перехода от дифференцирования по времени к дифференцированию по угловой координате, которое выглядит как

d = l d . Согласно этому правилу переходим от уравнения dt r2

движения (1.53) к уравнению его траектории:

 

2

2

r4

 

d ∂Φ

 

∂Φ

 

r′′ −

 

r′ − r =

 

 

 

 

,

(1.54)

 

 

 

 

 

r

 

l2 dt r

 

r

 

где штрихами обозначено дифференцирование по угловой координате .

В уравнении (1.54) удобно применить замену

r =

1

, r′ = −

1

λ′,

(1.55)

λ

λ2

 

 

 

 

19

приводящую его к виду

 

 

 

 

 

 

 

 

 

 

 

λ′′ + λ = −

1 1

 

d ∂Φ

∂Φ

 

 

 

 

 

 

 

 

 

.

(1.56)

l

2

 

λ

2

 

 

 

 

 

 

dt r

 

r

 

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

1.4. Потенциал Ньютона

Это самый известный и самый простой потенциал. Он имеет вид

Φ = −

k

(k = const).

(1.57)

 

 

r

 

Вносим его в уравнение (1.56):

 

λ′′ + λ =

 

 

k

.

(1.58)

 

 

 

 

 

 

 

 

l2

 

Обозначим

 

 

h =

l2

.

 

(1.59)

 

 

 

 

 

k

 

Тогда

 

λ′′ + λ =

1

.

(1.60)

 

 

 

 

 

 

 

h

 

Решением этого уравнения является одна из трех кривых второго порядка. В случае эллипса (рис. 1.1) в полярной системе координат решение записывается так:

λ =

1 + e cosϕ

,

 

(1.61)

 

 

 

 

 

 

 

 

h

 

 

r =

1

 

=

 

h

.

(1.62)

λ

1 + e cosϕ

 

 

 

 

Мы выбрали решение с одной постоянной интегрирования е (0 < e < 1), выполняющей роль эксцентриситета, так, чтобы полярная ось была совмещена с большой полуосью эл-

20

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