Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
АДУС.pdf
Скачиваний:
119
Добавлен:
21.03.2016
Размер:
2.37 Mб
Скачать

Перейти к оглавлению на странице: 256

§4. Решение задачи о движении точки в центральном поле методом Якоби

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

Уравнения Лагранжа второго рода

Различные аспекты задачи о движении точки в центральном поле мы рассмотрели в §§10, 13, 16 главы 6.

Сила, действующая на материальную точку M массы m в центральном поле сил в системе координат с началом в центре сил O задается формулой

~

= δ · Φ(r) ·

~r

 

(4.1)

F

r

,

а соответствующее уравнение Ньютона имеет вид (см. (10.2) гл. 6)

¨

~r

 

 

m~r = δ · Φ(r) ·

 

,

(4.2)

r

−−→

где ~r = OM , r = |~r|, δ = ±1, а Φ(r) модуль силы, действующей на рассматриваемую материальную точку.

Напомним, что законы сил, дающих примеры центральных силовых полей, мы рассмотрели в §4 главы 6, а в §10 главы 6 мы выписали для них соответствующие величины δ · Φ(r).

Центральное поле сил является потенциальным, причем потенциал и потенциальная энергия задаются формулами (§13 гл. 6)

Z

U(x, y, z) = (r) = δ · Φ(r)dr, Π(x, y, z) = − (r), (4.3)

167

Перейти к оглавлению на странице: 256

а тогда функция Лагранжа вычисляется по формуле:

L = T − Π =

m v2

+ (r).

(4.4)

2

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

Уравнения Лагранжа в декартовых координатах

Если в качестве обобщенных координат использовать декартовы координаты q = (q1, q2, q3) = (x, y, z), то

p

r = |q| = q12 + q22 + q32,

(4.5)

v = |q˙| = p12 + q˙22 + q˙32, L = m |q˙|2 /2 + (|q|).

Если в уравнения Лагранжа

 

d ∂L

∂L

= 0

(4.6)

 

 

 

 

 

dt

∂ q˙

∂ q

подставить выражение L по формуле (4.5) и произвести все дифференцирования слева, то в результате получим уравнения Ньютона (4.2).

Уравнения Лагранжа в сферических координатах

Введем в рассмотрение в качестве обобщенных координат точки M ее сферические координаты q = (q1, q2, q3) = (r, ϕ, ϑ), задаваемые формулами x = r sin ϑ cos ϕ, y = r sin ϑ sin ϕ, z = r cos ϑ. Пользуясь формулами (1.2), (1.5) главы 3, последовательно получаем, что

Hr = 1, Hϕ = r sin ϑ, Hϑ = r,

(4.7)

vr = r,˙ vϕ = r sin ϑ ϕ,˙

 

˙

vϑ = rϑ,

 

L = T − Π = m (r˙

2

+ r

2

2

2

2 ˙2

)/2 + (r).

(4.8)

 

 

sin ϑ ϕ˙

 

+ r ϑ

Упражнение 4.1. Выпишите уравнение Лагранжа в сферических координатах для рассматриваемой задачи, подставив функцию Лагранжа в (4.6).

168

Перейти к оглавлению на странице: 256

Канонические уравнения

Канонические уравнения относительно декартовых переменных

Для того, чтобы получить канонические уравнения для рассматриваемой задачи, введем импульсы p = (p1, p2, p3) и гамильтониан H по формулам (см. (1.1), (1.4), (1.5) и (4.4), (4.5)):

p =

∂L

= mq,˙

q˙(q, p, t) =

p

,

(4.9)

 

 

 

∂ q˙

m

H(q, p, t) = pq˙(q, p, t) − L(q, q˙(q, p, t), t) =

(4.10)

= p2/m − L(q, p/m, t) = p2/2m − (|q|).

Если в канонические уравнения

q˙ = ∂H/∂p, p˙ = −∂H/∂q

(4.11)

подставить H по формуле (4.10) и произвести справа дифференцирования, то в результате получим уравнения:

q˙ =

p

 

p˙ = δ · Φ (|q|) ·

q

(4.12)

 

,

 

.

m

|q|

Эти же уравнения получаются непосредственно из уравнений Ньютона (см. (4.2)), если использовать замену переменных q = (x, y, z), p = (mx,˙ my,˙ mz˙).

Канонические уравнения относительно сферических переменных

Используя формулу (4.8) для функции L в сферических переменных, введем импульсы p = (pr, pϕ, pϑ), соответствующие координатам q = (r, ϕ, ϑ), и гамильтониан H(q, p, t):

 

∂L

 

 

 

∂L

 

 

2

 

2

 

 

∂L

 

˙

2

 

 

pr =

 

= mr,˙

pϕ

=

 

= mϕr˙

 

sin

 

ϑ, pϑ =

 

= mϑr

,

(4.13)

 

∂ ϕ˙

 

 

˙

 

∂ r˙

 

 

 

 

 

 

 

 

 

 

∂ ϑ

 

 

 

 

 

 

 

 

 

pr

 

 

 

 

pϕ

 

 

˙

 

 

pϑ

 

 

 

 

r˙(q, p, t) =

m

, ϕ˙(q, p, t) =

mr2 sin2 ϑ

, ϑ(q, p, t) =

mr2

,

 

(4.14)

169

Перейти к оглавлению на странице: 256

H(q, p, t) = pq˙(q, p, t) − L(q, q˙(q, p, t), t) =

 

 

 

1

2

 

pϕ2

 

pϑ2

! − (r).

(4.15)

=

 

 

pr

+

 

+

 

 

2m

r2 sin2 ϑ

r2

 

 

 

 

Выпишите

канонические

уравнения в

 

Упражнение 4.2.

 

 

 

 

 

 

 

 

 

 

 

 

 

сферических переменных для рассматриваемой задачи, подставив функцию Гамильтона в равенства (4.11).

Уравнение Гамильтона-Якоби в декартовых и сферических переменных

Для рассматриваемой задачи о движении материальной точки в центральном поле уравнение Гамильтона-Якоби в декартовых переменных (см. (3.12) и (4.10)) имеет следующий вид:

 

 

 

∂ t + 2m

k=1 (∂S/∂ qk)2

!

 

q12 + q22 + q32

= 0, (4.16)

 

 

 

∂S

 

 

1

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

X

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ(r) dr (см. (4.3)).

 

 

 

 

 

 

 

 

 

где (r) = δ ·

 

 

 

 

 

 

 

 

 

 

 

 

Используя формулу (4.15), в случае сферических переменных

получаем другое уравнение Гамильтона-Якоби:

 

! − (r) = 0.

 

∂ t +

2m

 

∂ r

 

 

+ r2 sin2

ϑ

∂ϕ

 

 

+ r2

∂ϑ

∂S

1

 

 

 

∂S

 

2

1

 

 

 

∂S

 

2

1

∂S

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.17)

Общий интеграл уравнений движения

Так как уравнение Гамильтона - Якоби (4.17) явно от t и ϕ не зависит, то его полный интеграл будем искать в виде:

 

 

 

 

 

 

 

S = a1t + a2ϕ + Z(r, ϑ).

 

(4.18)

 

Подставляя (4.18) в (4.17), получаем:

 

!

− (r) = 0.

 

a1

+ 2m

∂ r

 

 

+ r2 sin2 ϑ

+ r2

∂ϑ

(4.19)

 

1

 

 

∂Z

 

2

 

a22

1

∂Z

 

2

 

 

170

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Перейти к оглавлению на странице: 256

 

Функцию Z будем искать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z(r, ϑ) = Z1(r) + Z2(ϑ),

 

 

(4.20)

что, после подстановки (4.20) в (4.19) дает равенство:

 

 

 

a1

+ 2m

dr

 

 

 

 

+ r2

 

 

 

sin2

ϑ +

!!

− (r) = 0. (4.21)

 

1

 

dZ

1

 

2

 

 

 

 

1

 

 

 

 

a2

 

 

 

 

 

 

dZ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

Полагая здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

+

dZ

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= a3,

 

 

(4.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin2 ϑ

 

 

 

получаем:

 

a1

 

+ 2m

dr

 

+ r23

! − (r) = 0.

(4.23)

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

dZ

2

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система (4.22), (4.23) решается в квадратурах:

 

 

 

 

 

Z1 = e1

Z r

 

 

 

 

 

 

dr,

e1 = ±1,

(4.24)

 

 

2m ( (r) − a1) − r23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

Z2 = e2

Z s

 

 

 

 

 

 

e2 = ±1.

(4.25)

 

 

 

 

 

a32 sin22

ϑ dϑ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно (4.18), (4.20), мы получили:

 

 

 

 

 

 

S = a1t+a2ϕ+e1

 

Z r

 

 

dr+e2

Z s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a32

sin22 ϑ dϑ.

 

2m ( (r) − a1) − r23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.26)

Упражнение 4.3. Покажите, что функция (4.26) удовлетво-

ряет условиям теоремы Якоби и примените эту теорему для нахождения общего интеграла уравнений движения материальной точки в центральном поле.

171