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

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

ГЛАВА 11. КАНОНИЧЕСКИЕ УРАВНЕНИЯ МЕХАНИКИ

§1. Вывод канонических уравнений

Преобразование Лежандра

В связи с уравнениями Лагранжа (1.4) главы 10 введем в рассмотрение величины

pi =

∂L

, i [1 : s],

(1.1)

∂ q˙i

которые называют импульсами или обобщенными импульсами. Как мы знаем, система уравнений Лагранжа состоит из s

обыкновенных дифференциальных уравнений второго порядка относительно обобщенных координат q1, . . . , qs . Существует бесконечно много способов сведения ее к системе из 2s обыкновенных дифференциальных уравнений первого порядка. Для этого вместо s переменных вводят 2s каких-то новых переменных. В частности, к s переменным q1, . . . , qs можно добавить еще s дополнительных переменных. Мы так и поступили, введя дополнительные переменные p1, . . . , ps . Как мы увидим, уравнения относительно переменных q1, . . . , qs , p1, . . . , ps имеют удобную симметричную форму так называемых канонических или гамильтоновых уравнений, для которых развита содержательная математическая теория.

Так как L = T − Π, а потенциальная энергия Π не зависит от обобщенных скоростей q˙1, . . . , q˙s , то истинно равенство:

 

2L

=

 

2T

, i, k

[1 : s].

(1.2)

 

∂ q˙i∂ q˙k

∂ q˙i∂ q˙k

Так как

 

 

2T

 

 

 

 

det

 

 

 

6= 0,

(1.3)

 

∂ q˙i∂ q˙k i,k=1

 

 

 

 

 

 

 

s

 

 

то равенства (1.1) можно

разрешить

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

величин

1, . . . , q˙s :

 

 

 

 

 

 

 

 

i = q˙i(q, p, t), i [1 : s],

q = (q1, . . . , qs), p = (p1, . . . , ps).

 

 

 

 

 

 

 

 

 

(1.4)

154

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

Введем в рассмотрение функцию

s

 

H(q, p, t) = Xk=1 pk k(q, p, t) − L(q, q˙(q, p, t), t),

(1.5)

которую называют функцией Гамильтона или гамильтонианом. Используя формулы (1.1), (1.5), получаем равенства:

 

 

∂ pi

 

 

 

s

pk

∂ pi

 

∂ q˙k

∂ pi = q˙i,

(1.6)

 

 

= q˙i + k=1

 

 

 

∂H

 

X

 

 

 

∂ q˙k

 

 

∂L ∂ q˙k

 

 

 

 

 

∂ qi

 

s

 

 

 

∂ q˙k ∂ qi

∂ qi

= − ∂ qi ,

(1.7)

= k=1 pk ∂ qi

∂H

 

X

 

∂q˙k

 

 

∂ L ∂ q˙k

 

∂L

 

∂ L

 

 

∂ H

 

s

pk tk

∂L ∂ q˙

 

 

∂L

 

∂L

(1.8)

 

∂ t

= k=1

∂ q˙k

 

∂ tk ∂ t

= − ∂ t

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при i [1 : s].

Из уравнений (1.4) главы 10 и формул (1.1) и (1.7) выводим,

что

 

 

 

 

 

 

d

∂L

 

∂H

 

 

 

 

pi =

 

= −

 

, i [1 : s],

(1.9)

dt

∂ qi

∂ qi

и, присоединяя сюда равенства (1.6), получаем искомые канонические или гамильтоновы уравнения:

 

d

∂H

 

d

∂H

 

 

 

 

qi =

 

,

 

pi = −

 

, i [1 : s].

(1.10)

dt

∂ pi

dt

∂ qi

Переход от уравнений Лагранжа II рода к уравнениям (1.10) при помощи введения обобщенных импульсов (1.1) и функции Гамильтона (1.5) называют преобразованием Лежандра.

Равенства (1.10) задают канонические уравнения алгоритмически (вспомним, что так же алгоритмически задавались и уравнения Лагранжа второго рода равенствами (1.1) и (1.4) главы 10). Чтобы выписать канонические уравнения для данной механической системы, необходимо составить для нее по формуле (1.5) гамильтониан H(q, p, t) и подставить эту функцию в правую часть уравнений (1.10), произведя там необходимые дифференцирования.

155

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

Преобразование Дирака

Рассмотрим систему обыкновенных дифференциальных уравнений

i = fi (q1, . . . , qn, t) , i = 1, . . . , n,

(1.11)

или, в векторной форме:

 

q˙ = f(q, t),

(1.12)

при f = (f1, . . . , fn), q = (q1, . . . , qn).

Преобразование Дирака произвольной системы обыкновенных дифференциальных уравнений вида (1.11) состоит в том, что в дополнение к “координатам” q = (q1, . . . , qn) вводятся в рассмотрение

“импульсы” p = (p1, . . . , pn) и гамильтониан H

по формуле:

 

 

 

 

 

n

 

 

 

 

 

H(q, p, t) = f · p = Xi=1 fi(q, t) · pi.

 

(1.13)

Это приводит к гамильтоновой системе:

 

 

 

∂ H

∂ H

 

n

 

 

 

Xj

 

 

 

 

 

 

 

 

 

i = ∂ pi

= fi(q, t), p˙i = ∂ qi

= −fq0 · p = −

· pj.

(1.14)

(fj)q0 i

 

 

 

 

 

=1

 

 

§2. Первые интегралы канонических уравнений

Интеграл механической энергии

Используя канонические уравнения (1.10), получаем:

 

dt

= ∂ t

+ i=1

∂ qi

i + ∂ pi

i =

 

dH

 

∂H

 

s

∂H

 

 

∂H

 

 

 

= ∂ t

 

 

 

P

 

 

 

 

 

(2.1)

+ i=1 ∂ qi ∂ pi

∂ pi ∂ qi = ∂ t .

 

∂H

 

s

 

∂H ∂H

 

∂H ∂H

 

∂H

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

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

H(q, p) = h,

(2.2)

156

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

на любом решении q, p этих уравнений, причем величина h не меняется со временем, ее называют постоянной механической энергии.

Используя формулы (см. (1.1), (1.5) и (1.3), (2.2) главы 10)

 

 

s

 

 

 

 

 

Xi

pii(q, p, t) − L(q, q˙(q, p, t), t),

(2.3)

 

H(q, p, t) =

 

 

=1

 

 

 

pi =

∂L

, i [1 : s], L = T − Π = T2 + T1 + T0 − Π

(2.4)

 

∂ q˙i

и легко проверяемое равенство

 

 

 

s

∂L

 

 

 

Xi

 

i = 2T2 + T1,

(2.5)

 

 

 

∂ q˙

 

 

=1

i

 

 

 

 

 

 

получаем интеграл механической энергии в форме ЯкобиОстроградского:

H = T2 − T0 + Π.

(2.6)

Если соотношения (1.1) главы 9, при помощи которых введены лагранжевы координаты q = (q1, . . . , qs), стационарны, то из формул (2.3), (2.4) главы 10 получаем, что T1 = T0 = 0, T2 = T , и тогда интеграл механической энергии принимает вид:

H = T + Π.

(2.7)

Циклические координаты и отвечающие им первые интегралы

Пусть функция Лагранжа L = L(q, q,˙ t) не зависит явно от координат q1, . . . , qσ при σ < s. Такие координаты называются циклическими, а остальные s − σ координат qσ+1, . . . , qs называются тогда позиционными.

Из формулы (1.7) следует, что

∂H

qν = 0, ν [1 : σ],

атогда из канонических уравнений (1.10), получаем:

pν =

∂L(q, q,˙ t)

= cν, ν [1 : σ],

∂ q˙ν

(2.8)

(2.9)

157

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

где все величины cν не меняются со временем, а это означает, что функции pν = ∂L(q, q,˙ t)/∂ q˙ν , ν [1 : σ] являются первыми интегралами канонической системы. Говорят, что pν первый интеграл, отвечающий циклической координате qν .

Для позиционных координат можно выписать следующую замкнутую систему канонических уравнений:

 

d

 

qi =

 

H(qσ+1, . . . , qs, c1, . . . , cσ, pσ+1, . . . , ps),

 

 

 

 

 

 

dt

∂ pi

 

d

 

 

 

(2.10)

pi = −

H(qσ+1, . . . , qs, c1, . . . , cσ, pσ+1, . . . , ps),

 

 

 

dt

∂ qi

i [σ + 1 : s].

Если позиционные координаты qσ+1, . . . , qs найдены, то, подставив их в равенства (2.9), можно найти циклические координаты. Действительно, из критерия Сильвестра положительной определенности квадратичной формы T2(u) (см. (2.7), (2.8) главы 10) следует, в частности, что истинно неравенство

 

 

2L

σ

 

 

det

 

i,k=1

> 0,

(2.11)

∂ q˙i∂ q˙k

которое является достаточным условием разрешимости равенств (2.9) относительно величин q˙1, . . . , q˙σ .

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

Метод Пуассона построения первых интегралов

Скобки Пуассона

Пусть f, g скалярные дважды непрерывно дифференцируемые функции канонических аргументов q = (q1, . . . , qn), p = (p1, . . . , pn) и времени t.

Положим

n

∂ qi

 

∂ pi

∂ pi

 

∂ qi .

(2.12)

{f, g} = i=1

 

 

X

 

∂ f

 

∂ g

 

∂ f

 

∂ g

 

Величина {f, g} носит название скобки Пуассона

функций

f, g . Она удовлетворяет следующим свойствам:

158

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

(а) {f, g} = −{g, f}, {f, f} = 0;

(б) если c = const, то {f, c} = {c, f} = 0;

(в)

{f1 + f2, g} = {f1, g}+ {f2, g}, {f1f2, g} = f1{f2, g}+ f2{f1, g};

(г)

 

∂ t{ϕ, ψ} =

∂ t

, ψ +

ϕ, ∂ t ;

 

 

∂ ϕ

 

 

∂ ψ

 

{f, qk} = −

 

∂ f

 

=

∂ f

(д)

 

,

{f, pk}

 

;

∂ pk

∂ qk

(е)

{qi, qk} = 0, {pi, pk} = 0, {qi, pk} = ±1;

(ж)

{f, {g, h}} + {g, {h, f}} + {h, {f, g}} = 0 тождество Пуас-

 

сона.

 

 

 

 

 

 

Упражнение 2.1. Доказать свойства (а)–(е).

Доказательство. Докажем свойство (ж). Как следует из опреде-

ления, скобка Пуассона {f, g} есть билинейная однородная функция производных первого порядка от f и g . Тогда ясно, что левая часть тождества Пуассона есть линейная однородная функция производных второго порядка от f , g , h, то есть каждый член тождества должен содержать одну из этих производных второго порядка. Тождество будет доказано, если мы покажем, что оно не содержит ни одной производной от f, g, h второго порядка. В силу симметрии тождества по f, g, h, достаточно доказать, что оно не содержит вторые производные от одной из функций f, g или h, например от f .

Соберем в тождестве все члены, содержащие вторые производные от f . Они содержатся только во второй и третьей скобках. Введем обозначения

D1(ϕ) = {g, ϕ}, D2(ϕ) = {h, ϕ},

(2.13)

тем самым мы вводим в рассмотрение два линейных дифференциальных оператора

2n

 

 

 

2n

 

 

 

 

X

 

, D2 =

X

 

 

 

D1 = ξk

 

 

ηk

 

 

,

(2.14)

k=1

∂ x

k

k=1

∂ x

k

 

 

 

 

 

 

 

 

 

 

 

159

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

где x = (x1, . . . , x2n) = (q1, . . . , qn, p1, . . . , pn), а ξk = ξk(x), ηk =

ηk(x) некоторые функции своих аргументов. Используя обозначения (2.13) получаем

{g, {h, f}} + {h, {f, g}} = {g, {h, f}} − {h, {g, f}} =

(2.15)

= D1(D2(f)) − D2(D1(f)) = (D1D2 − D2D1) (f)

Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D1D2

= ξkηl

 

 

 

 

2

 

 

+ ξk

∂ ηl

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ x

 

∂ x

∂ x

 

 

 

 

 

 

 

k,l

 

k

 

k,l

 

∂ x

 

 

 

 

 

 

 

l

 

k

 

 

 

l

 

 

P

 

 

 

2

 

P

 

 

 

 

 

 

 

 

 

(2.16)

 

P

 

 

 

 

P

∂ ξl

 

 

 

 

 

 

D2D1 = ηkξl

 

 

 

 

 

 

+ ηk

 

 

 

 

 

,

 

 

∂ x

k

∂ x

∂ x ∂ x

 

 

 

k,l

 

 

 

 

l

 

k,l

k

 

 

 

l

 

то

 

 

 

ξk ∂ xk

− ηk

 

 

∂ xl .

 

D = D1D2 D2D1 = k,l

∂ xk

(2.17)

 

X

 

 

 

 

∂ ηl

 

∂ ξl

 

 

 

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

Что и требовалось.

Теорема Пуассона

Скобки Пуассона мы ввели формально. Выясним, откуда они берутся, как связаны с системой Гамильтона. Рассмотрим уравнения Гамильтона:

q˙ =

∂ H(q, p, t)

, p˙ = −

∂ H(q, p, t)

,

(2.18)

∂ p

 

∂ q

где q = (q1, . . . , qn), p = (p1, . . . , pn).

Пусть f = f(q, p, t) скалярная функция от q, p и t. Вычислим ее полную производную по t:

dt =

∂ t

n

∂ qk

k + ∂ pk k .

(2.19)

+ k=1

df

∂ f

X

 

∂ f

 

∂ f

 

160

 

 

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

Используя здесь (2.18), получаем

 

df

∂ f

 

 

 

=

 

+ {f, H}.

(2.20)

 

dt

∂ t

Отсюда следует необходимое и достаточное условие того, чтобы функция f(q, p, t) была первым интегралом уравнений (2.18):

∂ f

+ {f, H} ≡ 0.

(2.21)

∂ t

Пуассон предложил метод, позволяющий (к сожалению не всегда) по двум известным интегралам гамильтоновых уравнений строить третий.

Теорема 2.1. (Теорема Пуассона) Пусть ϕ(q, p, t), ψ(q, p, t)

первые интегралы канонической системы (2.18). Тогда функция f = {ϕ, ψ} также будет интегралом этой системы или f ≡ const.

Доказательство.

Если ϕ, ψ первые интегралы

уравнений

 

 

 

 

 

 

 

 

 

 

 

(2.18), то

 

 

 

 

 

 

 

 

∂ϕ

∂ψ

 

 

 

 

 

 

 

+ {ϕ, H} = 0,

 

 

+ {ψ, H} = 0.

(2.22)

 

 

∂ t

∂ t

 

Выпишем тождество Пуассона для H, ϕ, ψ :

 

 

{H, {ϕ, ψ}} + {ϕ, {ψ, H}} + {ψ, {H, ϕ}} = 0,

(2.23)

 

Из (2.22) следует, что

 

 

 

 

 

 

 

 

 

 

 

∂ϕ

∂ψ

 

 

{H, ϕ} = −{ϕ, H} =

 

 

, {ψ, H} = −

 

.

(2.24)

 

∂ t

∂ t

 

Подставляя эти скобки в тождество (2.23), получаем последо-

вательно:

 

 

 

 

 

 

 

{H, {ϕ, ψ}} + {ϕ, −∂ψ∂ t } + nψ, ∂ϕ∂ t o = 0,

 

 

{{ϕ, ψ} , H} + nϕ, ∂ψ∂ t o + n∂ϕ∂ t , ψo = 0,

(2.25)

{ϕ,ψ} + {{ϕ, ψ} , H} = 0,

t

ft + {f, H} = 0,

161