
книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений
.pdf1.3. Формализм Гамильтона |
|
17 |
|
Уравнения преобразования имеют вид |
|
|
|
z = b sin Ѳ, |
|
|
|
X = а cos гр + Д sin ф, |
ф = |
соt, |
(1.13) |
у = a sin ф — А cos тр, |
А = |
Ъcos Ѳ. |
|
Существует только одна обобщенная координата, Ѳ. Остальные переменные удовлетворяют записанным выше уравнениям связи. Кривая, изображающая путь частицы в трехмерном простран стве, задается параметрическими уравнениями:
z = Ъsin Ѳ,
X = а cos соt Ъcos Ѳ sin соt, |
(1.14) |
у = a sin co£ — Ъcos Ѳ cos cof,
где решение Ѳ = Ѳ(t) завершает параметризацию. Уравнения (1.14) являются уравнениями преобразования особого вида — неста ционарного преобразования. С важными свойствами таких пре образований мы вскоре встретимся.
Используя уравнения (1.14) и проводя ряд длинных вычисле ний, получаем выражение для Т, а следовательно, и для L (= Т). Теперь можем записать уравнения Лагранжа. Остается, однако, неисследованной интересная проблема устойчивости. Мы рассмот рим ее в задаче 1.8.
Задача 1.8. Из рис. 1.3 видно, что существуют четыре «равно весных» решения. Исследуя уравнения движения, убедимся
в этом. Затем подставим решения в форме Ѳ(1) = Ѳ4 + ееѴ в урав |
||
нения движения и пренебрежем членами порядка е2. Получим |
||
условие для X (в виде алгебраического уравнения), которое долж |
||
но удовлетворяться, чтобы записанные выше выражения для Ѳ |
||
действительно были решениями. Исследуя корни характеристиче |
||
ского |
уравнения^!, |
. . ., Х4, можно судить об устойчивости равно |
весных решений. Например, если ^ > 0, то любое бесконечно |
||
малое |
возмущение |
равновесного решения Ѳ экспоненциально |
растет со временем, т. е. равновесие не является устойчивым.
Задача 1.9. Пренебрегая силой тяжести, решить задачу 1.5, когда ось винтовой линии расположена в горизонтальной пло скости и вращается относительно фиксированной оси с угловой скоростью со.
1.3. Формализм Гамильтона
Уравнения Лагранжа являются лучшим средством для иссле дования многих задач, возникающих в классической механике. При подробном рассмотрении уравнений обнаруживается инте
ресный факт. Если L не зависит от qi, то величина дЫдді цостояаг,
2 -0124 3
«I
ч1
»V
18 Гл. I. Элементы классической механики
на по времени. Такие координаты дг называются циклическими. В механике существует метод, который не только использует этот факт, но и позволяет представить обширное поле механики как науку с принципами симметрии. Этот метод воплощен в фор мализме Гамильтона.
Первой задачей нашего нового подхода является переход
от переменных (q, q) к (q, dL/dq). Такой переход достигается при использовании преобразования Лежандра, которое переводит
L — функцию g, q и t — в гамильтониан Н , являющийся функ
цией g, dL/dq, t. Преобразование имеет вид
N
(1-15)
і=і Hi
Образуя дифференциал от Н, легко проверить, что Н действитель
но является функцией g, dL/dqi и t. Если мы введем
^ r ^ P i |
(1-16) |
Hi |
|
— канонические импульсы, или импульсы, сопряженные координа там qh то Н = II (д, р, t).
Чтобы найти уравнения движения системы через эту новую ди намическую функцию, мы сначала запишем дифференциал Н как функцию д, р (и t, если Н явная функция времени). Это дает
N |
N |
|
dH = 2 |
9z dpt — 2 Гг dqt— dt, |
(1.17) |
i=i |
г=і |
|
где при построении дифференциала мы воспользовались уравне ниями Лагранжа. Отсюда следует, что
|
дН |
(1.18а) |
|
орі |
|
|
|
|
Рі = |
dH |
(1.186) |
dqi ' |
и dH/dt = —dL/dt. Если L не зависит явно от времени, то имеют смысл только уравнения (1.18). Они называются уравнениями Гамильтона и, как мы видим, представляют собой 2N уравнений первого порядка в противоположность уравнениям Лагранжа, которые являются N уравнениями второго порядка.
Из уравнения (1.186) видно, что если гамильтониан не зависит от координаты qi (т. е. ф циклическая), то импульс рі, сопряжен ный этой координате, сохраняется. Этот факт, связанный с неотъ-
1.3. Формализм Гамильтона |
19 |
емлемыми свойствами симметрии пространства и времени, приво дит к основным теоремам сохранения в механике. Сначала мы исследуем свойство однородности пространства.
Рассмотрим изолированную систему из N частиц, координаты которых относительно некоторой фиксированной системы коорди нат равны г;. Пусть система может иметь поступательное виртуаль ное перемещение А, так что координаты гг -у г\ = гг + А. Урав нения движения системы после такого перемещения не отличают ся от уравнений, которые описывали систему в первоначальной ориентации. Простой эксперимент покажет, что принцип одно родности пространства будет нарушен, если эти уравнения ока жутся различными. Если уравнения движения инвариантны при указанных перемещениях, то этим же свойством должен обладать и гамильтониан, т. е.
e * = s - § £ - e * = A - s 4 “ °- |
(I-«) |
Согласно определению виртуальных перемещений, ни импульсы, ни время не варьируются. Так как приращение А — произвольно, то наш эксперимент означает, что гамильтониан изолированной системы подчиняется условию
N
Объединяя это свойство с уравнениями Гамильтона, получаем желаемый результат:
• s r ( 2 » = ° -
Что представляют собой г; и рг? Ясно, чтогг являются радиус-век торами частиц, составляющих систему, так что кинетическая энергия Т равна
N |
(1.21) |
т=4 2 "^ ? . |
Поскольку потенциал взаимодействия частиц от скорости не зави сит, импульс
Pi = - ^ - = m«rj, |
(1-22) |
дТі |
|
т. е. это обычный линейный импульс.
Таким путем устанавливается, что однородность пространства означает постоянство по времени полного линейного импульса изолированной системы.
2*
20 |
Гл. I. Элементы классической механики |
Другим исследуемым свойством пространства является его изотропность. Это свойство означает, что уравнения движения изолированной системы инвариантны относительно виртуального бесконечно малого поворот системы вокруг некоторой фиксированной оси. Такой поворот показан на рис. 1.4. Каждая частица испытывает перемещение, при котором гг ->■ ѵ\ —
=гг + бгI, где бгі = бѳ X гI. Снова гамильтониан должен быть
инвариантным; таким образом,
(1.23)
Поскольку 6Ѳ произвольна, мы полу чаем второе важное свойство гамиль тониана изолированной системы:
N
У,чх дЩдті 0. (1.24)
Это соотношение совместно с уравнением (1.22) и кинематиче
ским равенством d (гг х r t)/dt = rt |
х (drt/dt) дает |
желаемый |
результат: |
• |
|
N |
|
|
г*X Рі ) == 0, |
(1.25а) |
|
N |
|
|
З г гХр г = |
Ж. |
(1.256) |
Полный угловой момент X (относительно произвольного фикси рованного начала координат) постоянен во времени.
Задача 1.10. Как изменяется полный угловой момент систе мы при преобразовании координат гг К + а, где а — постоян ный вектор?
Задача 1.11. Построить в декартовых координатах (исходя из лагранжиана) гамильтониан для свободной частицы, движу щейся в трехмерном пространстве. Какие величины сохраняются при таком движении?
Задача 1.12. Решить задачу 1.11, выражая*# в сферических координатах (г, Ѳ, ф). Построить уравнения движения (т. е. уравнения Гамильтона) для координат, не являющихся цикли ческими.
1.3. Формализм Гамильтона |
21 |
Задача 1.13. Построить гамильтониан для частицы, движу щейся по поверхности цилиндра радиуса а, если ось цилиндра горизонтальна. Сила тяжести направлена вниз. Найти движение при следующих начальных условиях (рис. 1.5):
z (0) = 0, Ѳ (0) = я/2, z (0) = V, è (0) = 0.
Мы подошли теперь к важной теореме о гамильтониане: для консервативных систем, связи которых не зависят от времени, гамильтониан равен полной энергии Т + V и постоянен по вре мени. Точнее, условие независимости связей от времени означает,
что уравнения преобразования (ср. с (1.10)), общий вид которых есть гг = Ti (g1? . . ., qN, t) (индекс i связан с полным числом частиц в системе), явно от времени не зависят (заметим, что урав нения (1.14) являются хорошим примером связи, зависящей от времени).
Так как силовое поле консервативно и связи не зависят от вре мени, то лагранжан L не содержит время явно. Следовательно,
* dL _
dt
II
N
= 2
/ dL |
■fl« |
Г“ |
|
\ dqi |
|||
|
|
d / dL \ dt \ dqi .
dL dqi \ = dqi dt j
’ . |
dL |
dqj |
N |
dL |
|
|
d I |
(1.26) |
|||||
4r+ |
dqi |
dt |
г |
* |
||
|
При получении этого результата мы использовали уравнения Лагранжа. Результат можно переписать в виде
(1.27)
Таким образом,
Я = const. |
(1.28) |
Кроме того, если потенциал V не является функцией скорости,
то рг = dLldqi = dTldqi и сумма в уравнении (1.27) имеет вид
2<?г (dTldqi). Если уравнения преобразования не содержат время
22 |
Гл. I. Элементы классической механики |
|
явно, то |
Т будет однородной квадратичной формой |
скоростей: |
|
Т = 2 |
(1.29) |
|
I, k |
|
Коэффициенты aih не содержат обобщенных скоростей.
Задача 1.14. Исходя из
N
т=\ 2
г=і
(N — число частиц в системе) и уравнений преобразования в фор ме Г; = гг (qi, . . ., qN), доказать, что коэффициенты aik не содер жат обобщенных скоростей.
Из теоремы Эйлера, которая утверждает, что если g является
однородной функцией порядка ѵ от переменных хи то ^jX^dg/dzi) —
= vg (для нашего случая ѵ = 2), следует, что
= |
(1.30) |
dqi |
|
и |
|
H = ' Z i lpl - L = 2 T - ( T - V ) = T + V |
(1.31) |
— полная энергия системы. Для консервативной системы с неза висящими от времени связями гамильтониан равен полной энергии и постоянен во времени. Отметим, что изолированная система не может описываться уравнениями преобразований, явно содер жащими время. В таком случае гамильтониан также не будет явной функцией времени, и Н сохраняется. Тот факт, что гамиль тониан изолированной системы не может содержать явно время, выражает однородность времени. В данном случае, как мы уже видели, Н представляет собой константу движения, а именно энергию.
Время, также как и пространство, обладает еще одним фунда ментальным свойством — свойством изотропности. Этого рода симметрия ведет к принципу динамической обратимости. Если одномерный объект, такой, как время, является изотропным, то не существует различия между временем прогрессирующим и вре менем регрессирующим. При обращении времени законы движе ния остаются инвариантными, так что Н (t) = Н (—t). В свете уравнений Гамильтона это означает, что если [q (t), р (t)] являет ся динамическим решением, то [q (—t),— р (—£)] также будет решением.
Чтобы сделать доказательство более убедительным, рассмот рим сначала начальные и конечные условия:
1.3. Формализм Гамильтона |
23 |
Далее найдем динамический путь, соответствующий новым началь ным условиям р * (0) = —pi, q* (0) = qx. Имеем решение
Р* = —Р i—t* +'h), Я* = q (—t* + О 7
где t* есть время в задаче с переменными, отмеченными звездочкой (заметим, что р и q — те же самые функции, что и в первом слу чае). Таким образом, мы нахо дим, что
Р* (t* = 0) = —р (h) = —ри q* (0) = qu
так что начальные |
условия вы |
|
||||
полняются. |
|
|
|
|
||
При t* |
= ti имеем р* (£t) = |
|
||||
= — ро |
и |
q* iti) = q0. Следова |
|
|||
тельно, полученное решение яв |
|
|||||
ляется |
точным обращением ре |
|
||||
шения [q (і), р (£)]. Уравнения, |
|
|||||
которые порождают |
решения с |
|
||||
такими свойствами, |
называются |
|
||||
динамически обратимыми. |
На |
|
||||
рис. 1.6 |
изображены два таких |
|
||||
движения. |
Каждому |
решению |
|
|||
динамических уравнений |
соот |
|
||||
ветствует |
другое |
обращенное |
|
|||
решение. |
|
|
|
|
Р и с . 1.6. Движение а и динами |
|
Задача |
1.15. Как в задаче 1.6 |
чески обращенное движение б. |
изменяется по времени энергия шарика? Если эта энергия возрастает (убывает), то откуда (куда) идет прирост (убыль) энергии?
Утверждение, что некая динамическая переменная постоянна по времени, означает, что существует функция координат и им пульсов, которая не изменяется при изменении qi и рі по времени. Такая функция
^ (?І7 • • •? qNi Pii ■ • ч Р Nt
постоянная вдоль «динамического пути системы», называется константой движения. Эта концепция, в сущности, является
•одной из наиболее важных во всей классической динамике. Сколь ко независимых констант движения существует для системы с N обобщенными координатами? Из второго закона Ньютона либо из уравнений Лагранжа или Гамильтона мы видим, что решение включает 2N постоянных интегрирования. Обращая эти решения так, что постоянные интегрирования становятся функциями дина мических переменных, получаем совокупность 2N независимых констант движения. Разрешенные относительно индивидуальных
24 Гл. I. Элементы классической механики
координат и импульсов, эти 2N первых интеграла уравнений движения, как мы видим, будут эквивалентны решению задачи. Существование 2N независимых констант движения можно также усмотреть из геометрического построения.
2іѴ-мерное декартово пространство с координатами (д1? . . ., qN; Pi, . . ., p N) называется фазовым пространством (а также Г-про- странством). В этом пространстве вся система изображается одной
|
|
|
|
точкой. Состояние системы в любой |
||||
|
|
|
|
момент есть точка Г-пространства. |
||||
|
|
|
|
Так как система изменяется по |
||||
|
|
|
|
времени, то точка системы прочер |
||||
|
|
|
|
чивает путь (динамическую траек |
||||
|
|
|
|
торию системы) в Г-пространстве. |
||||
|
|
|
|
Однако в этом пространстве вре |
||||
|
|
|
|
менной параметр отсутствует. Та |
||||
|
|
|
|
ким |
образом, |
кривая |
в Г-про |
|
|
|
|
|
странстве не показывает временную |
||||
|
|
|
|
эволюцию точки системы. Она дает |
||||
|
|
|
|
соотношение, |
которому |
должны |
||
|
|
|
|
удовлетворять координаты и им |
||||
|
|
|
|
пульсы системы. |
|
|||
|
|
|
|
Решение задачи для одномер |
||||
Р и с . |
1.7. |
Динамическая |
траек |
ного гармонического осциллятора |
||||
можно записать в виде p = acosa>t, |
||||||||
тория |
точки системы для |
просто |
||||||
го гармонического осциллятора в |
X = bsinco^. В |
соответствующем |
||||||
трехмерном Г-пространстве. |
двумерном Г-пространстве дан |
|||||||
рическим |
представлением |
ные |
уравнения будут |
парамет |
||||
эллипса |
(р/а)2 +, (x/b)2 = |
1. Хотя |
уравнение этой динамической кривой является важным соот ношением, которому должны удовлетворять р и х , она не дает информации о развитии системы во времени. С другой стороны,
в (2N + 1)-мерном расширенном фазовом пространстве (Г-про странстве), которое включает ось времени, кривая, прочерчивае мая точкой системы, дает полное динамическое решение задачи. Для рассмотренной задачи с гармоническим осциллятором эта кривая является пересечением двух поверхностей а = р/cos соt, Ъ — xlsmuit. Динамический путь изображен на рис. 1.7.
Для одномерного движения свободной частицы кривая системы
в трехмерном Г-пространстве представляет собой пересечение
двух плоскостей р |
= р 0 и х 0 — х — v0t, как показано на рис. 1.8. |
||
Обычно решением для системы с N обобщенными координатами |
|||
служит |
кривая в |
(2N + 1)-мерном |
Г-пространстве. «Кривая» |
в таком |
пространстве соответствует |
пересечению 2N «поверхно |
стей» в этом же пространстве. Кривая в пространстве с любым числом измерений является одномерным геометрическим местом
1.3. Формализм Гамильтона |
25 |
точек. В трехмерном пространстве два независимых |
уравнения |
/і (X, г/, z) = 0 и / 2 (х, у, z) = 0 определяют кривую, так как только одна координата точки этого многообразия может быть выбрана произвольно. Геометрически кривая является пересече нием двух соответствующих поверхностей. Поверхность в іѴ-мер- ном пространстве есть (N — 1)-мерное многообразие. Таким обра-
Р и с. 1.8. Динамическая траектория свободной частицы в трехмерном Г- пространстве.
зом, в трехмерном пространстве уравнение g (х, у, z) = 0 опре деляет поверхность, так как только две координаты могут быть выбраны произвольно.
В (2N + 1)-мерном Г-пространстве уравнение
h (qi, . . ., qд.; pi, . . ., p N, t) — C1= 0
представляет собой поверхность. Пересечение 2N таких различ ных поверхностей образует одномерное многообразие в Г-про
странстве, т. е. кривую. Но кривая в Г-пространстве может быть записана в параметрической форме дг = qt {х); p t = p t (х); t = х, что дает динамическое решение задачи. С другой стороны, 2N — 1 констант определяют двумерное многообразие, т. е. две пере менных могут быть выбраны произвольно. Обозначим их через Рі и t. Если рі и t независимы, то задача не решена. Итак, мы заключаем, что 2N независимых констант движения полностью определяют решение.
Хотя концепция констант движения является полезным геомет рическим понятием в теории неравновесной статистической меха ники, в любой задаче классической механики только часть этих констант имеет значение. Например, насколько важны начальные значения координат? После утомительного решения системы 2N
26 Гл. I. Элементы классической механики
дифференциальных уравнений такие константы, конечно, необхо димы, чтобы единственным образом определить движение. Но до этого они, несомненно, были не нужны. С другой стороны, как мы видели, существуют константы, которые следуют из неотъем лемых свойств симметрии системы. Эти функции координат и импульсов, постоянные вдоль динамической траектории систе мы, являются крайне важными. Во многих случаях с помощью этих констант решение можно представить в виде квадратур (мы вернемся к этой теме в нашем последнем обсуждении эргодической теоремы в гл. V).
Задача 1.16. а) Рассмотреть две частицы в трехмерном про странстве с потенциалом взаимодействия V (1 rt — г2 |). Выбрать обобщенные координаты таким образом, чтобы они включали координаты центра масс системы. Построить гамильтониан и пока зать, что задачу можно описать через движение относительно центра масс. Доказать также, что это движение происходит в пло скости, вектор нормали которой постоянен по времени.
б) Построить гамильтониан для плоского движения, выбирая за обобщенные координаты скаляр г = | г4 — г2 | и угол Ѳ между векторами г4 (0) — г2 (0) и rj (t) — г2 (t). Получить константы движения и свести задачу к квадратурам.
1.4. Скобки Пуассона и канонические преобразования
Рассмотрим некоторую динамическую характеристику систе мы, т. е. функцию и от координат, импульсов и времени. Скорость изменения этой функции со временем задается равенством
N
(1.32)
Подставляя уравнения Гамильтона (1.18) в это выражение, полу чим:
du |
N |
I |
ди |
dH |
ди |
дІІ |
|
ди |
|
|
__•sri |
\ |
(1.33) |
||||||||
dt |
г=іI |
\ |
dqi |
dpi |
dpi |
dqi |
) ' |
dt ' |
||
|
или, что эквивалентно,
du
dt = [u,H] + ^ . (1.34)
Скобки Пуассона для любых двух динамических переменных А и В определяются как
N |
дА |
дБ |
dB |
дА \ |
|
|
[А, В] = 2 ( |
(1.35) |
|||||
dqi |
dpi |
dqi |
dpI ) |
і= і