Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Тер. Мех. (Леушин,Нигматуллин,Прошин).pdf
Скачиваний:
856
Добавлен:
10.02.2015
Размер:
2.45 Mб
Скачать

Теоретическая физика. Механика (практический курс)

215

 

 

Примеры решения задач

Задача 3. Записать уравнение Гамильтона-Якоби для простейших механических систем:

a) свободной частицы; б) гармонического осциллятора; в) частицы в центральном поле.

Решение. Функции Гамильтона этих систем найдены в задаче 1 раздела 9. а) Функция Гамильтона свободной частицы (число степеней свободы n = 3)

 

 

 

 

 

 

 

 

 

 

 

 

p2

 

+ p2

+ p2

 

H (x, y, z, px , py , pz ) =

 

x

 

 

 

y

 

z

.

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходное уравнение Гамильтона-Якоби для полной функции действия S –

 

1

S 2

S

2

S

 

2

 

S

= 0 .

 

 

 

 

 

 

 

+

y

+

 

z

 

 

 

+

t

 

 

 

 

 

2m

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пользуемся цикличностью координат y, z и времени t, получаем связь

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S = S(x) + α2y + α3z Et

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

и упрощенное уравнение для сокращенного действия S(x)

 

 

 

 

1

 

%

2

2

 

 

 

2

 

 

 

 

 

 

 

 

 

dS

 

+ α2

+ α3 = E .

 

 

 

 

 

 

 

2m

 

 

 

 

 

 

 

dx

 

2m

 

 

 

 

 

 

 

 

б) Функция Гамильтона гармонического осциллятора (n = 1)

 

 

 

 

H(x,p) = p2/2m + mω2 x2/2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

Исходное уравнение Гамильтона-Якоби для полной функции действия S –

 

 

 

 

1 S 2

 

mω2

x

2

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

+

0

 

+

t

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m

x

 

2

 

 

 

 

 

 

 

 

Пользуемся консервативностью системы (цикличностью времени t) и переходим к укороченному действию S0

S = S0(x) Et.

Для S0 получаем упрощенное уравнение

Канонические преобразования. Уравнение Гамильтона-Якоби

216

 

 

 

 

dS0

2

= 2mE m2ω2 x2 .

 

 

 

0

 

dx

 

 

 

в) Функция Гамильтона частицы в центральном поле (n = 2)

 

p2

 

p2

 

H =

ρ

+

ϕ

+U (ρ) .

2m

2mρ2

 

 

 

Исходное уравнение Гамильтона-Якоби для полной функции действия S –

1

S 2

+

1

 

S 2

+U (ρ) +

S

= 0 .

 

2mρ

2

t

2m ∂ρ

 

 

∂ϕ

 

 

Пользуемся цикличностью ϕ и t, выражаем полный интеграл S через со-

~

кращенное действие S

= ~ ρ + α ϕ −

S S( ) ϕ Et.

~ ρ

Упрощенное уравнение для S( ) выглядит следующим образом

dS% 2

αϕ2

 

 

= 2mE 2mU (ρ)

 

.

ρ2

dρ

 

 

Задача 4. Найти закон движения тела массы m брошенного с начальной скорость v0 под углом ϕ к горизонту в поле силы тяжести.

Замечание: Перед решением этой известной школьной задачи с помощью уравнения Гамильтона-Якоби вспомните ее решение "обычным" способом с помощью второго закона Ньютона.

Решение. Введем декартовую систему координат. Для описания движения тела достаточно 2 координаты x и y. Ось y направим вертикально вверх, ось x – горизонтально. Начальные условия: x0 = 0, y0 = 0, vx0 = v0cosϕ, vy0 = v0sinϕ. Далее действуем по алгоритму.

1.Гамильтониан. Задача имеет две степени свободы: q1 = x, q2 = y. Кинетическая энергия – T = (m/2)(x2 + y2) = T (2) – квадратична по скоростям, потенциальная – U = mgy – имеет обычный вид. Следовательно, функция Лагранжа

L(q1,q2,q1,q2,t) = L(x,y,x,y) = T U = (m/2)( x2 + y2) mgy.

Теоретическая физика. Механика (практический курс)

217

 

 

Обобщенные импульсы (в данном случае они совпадают с обычными!)

px =

L

&

py

=

L

&

 

 

&

= mx;

& = my .

 

 

 

x

 

 

 

 

y

 

 

 

Начальные условия в гамильтоновом формализме при t = t0 = 0

 

x = x0 = 0, y = y0 = 0;

px = px0 = mv0cosϕ, py = py0 = mv0sinϕ.

(10.25)

Запишем функцию Гамильтона

 

 

 

 

 

 

 

 

 

(2)

 

2

2

/2m + mgy

 

H(x, y, px, py, t) = pxx + pyy L = Т

 

+ U = px

/2m + py

 

2. Уравнение Гамильтона-Якоби для полной функции действия S(x,y,t)

1

S 2

1

 

S 2

S

 

 

 

 

+

 

 

 

+ mgy +

t

= 0 .

 

 

2m

x

 

2m

y

 

 

3.Упрощение (интегралы движения). Функция Гамильтона не зависит от времени t и от x. Следовательно, сохраняются обобщенная энергия H (свойство а) и импульс px (свойство б). Для их определения воспользуемся начальными условиями (10.25)

H = p2

/2m + p2

/2m + mgy = p2

/2m + p2

/2m = mv2

/2 = E,

(10.26)

x

y

x0

y0

0

 

 

 

 

px = px0 = mv0cosϕ = αx.

 

 

10.27)

Полное действие S(x,y,t) в соответствии со свойствами (10.14) и (10.16) имеет вид

 

 

 

~

 

αx, E) + αxxEt,

(10.28)

 

S = S(y,

 

 

 

 

 

 

~

 

а упрощенное уравнение для сокращенного действия S(q2) –

 

 

α2x

+

1

 

S% 2

+ mgy = E .

(10.29)

 

 

 

 

2m

2m

y

 

 

~

4. Решение уравнения Гамильтона-Якоби. Поскольку S функция только одной переменной y, то, заменив частную производную в (10.29) обычной

~

производной dS/dy и разрешив относительно нее уравнение (10.29), имеем dSdy% = 2mE 2m2 gy − α2x .

Находим решение, в котором опускаем аддитивную константу

Канонические преобразования. Уравнение Гамильтона-Якоби

218

 

 

 

 

 

%

= −

(2mE 2m2 gy − α2x )3/ 2

.

(10.30)

S

3m2 g

 

 

 

 

Подставляя (10.30) в (10.28), получаем полный интеграл уравнения Га- мильтона-Якоби – функцию, зависящую от 2 координат (x и y), 2 констант (αx и E) и времени t

S(x, y,αx , E,t) = −

(2mE 2m2 gy − α2x )3/ 2

+ αx x Et .

(10.31)

3m2 g

 

 

 

5.Закон движения – "производная по константе дает константу". Под-

ставляем найденный полный интеграл (10.31) в первые из соотношений (10.23) и получаем

S

=

 

S

=

αx (2mE 2m2 gy − α2x )1/ 2

+ x

,

(10.32а)

 

 

 

 

 

 

 

 

 

 

∂α1

 

 

 

∂αx

 

 

m2 g

 

1

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

=

S

= −

(2mE 2m2 gy − α2x )1/ 2

t 2 .

 

(10.32б)

 

∂α2

E

mg

 

 

 

 

 

 

 

 

 

 

Константы αx = α1 и E = α2 определены выше в (10.26), (10.27). Используем начальные условия для нахождения β1 и β2, для этого подставим в

(10.32) t = t0 = 0 и применим (10.25) – (10.27)

 

 

 

αx (2mE − α2x )1/ 2

v2 cosϕsin ϕ

 

v2 sin 2ϕ

 

 

β =

 

 

 

=

 

0

 

=

0

 

 

 

 

 

 

 

 

 

1

 

 

m2 g

 

 

 

 

g

 

2g

 

 

 

 

 

 

 

 

 

 

 

.

(10.33)

 

 

 

 

(2mE − α2x )1/ 2

 

 

 

 

 

 

 

β

 

= −

= −

v

sin ϕ

 

 

 

 

2

mg

0

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь все константы определены, решаем систему алгебраических уравнений (10.32) относительно неизвестных x и y

αx

2 2

 

m

g β2 2gy + gx = gβ1.

(10.34)

g2β22 2gy + gt = gβ2

 

Из второго уравнения определяем значение

g2β22 2gy , которое под-

ставляем в первое уравнение, и находим