- •ПРЕДИСЛОВИЕ
- •ЧАСТЬ I. ПРОСТРАНСТВА И КООРДИНАТЫ
- •Аффинные пространства
- •Аффинные евклидовы пространства
- •Аффинные координаты и преобразования
- •Криволинейные координаты
- •Криволинейные системы координат
- •Локальные базисы
- •ЧАСТЬ II. КИНЕМАТИКА
- •Кинематика точки
- •Коэффициенты Ламе. Проекции скорости точки на оси криволинейной системы координат
- •Проекции ускорения точки на оси ортогональной криволинейной системы координат
- •Описание движения точки в естественных координатах
- •Определение кривизны траектории точки по движению
- •Два примера движения точки
- •Кинематика твердого тела
- •Движение механической системы. Твердое тело. Число степеней свободы положения. Аффинное пространство и координаты, связанные с твердым телом
- •Группа движений аффинного евклидова пространства
- •Поступательное движение твердого тела
- •Вращение твердого тела вокруг неподвижной оси
- •Плоское движение твердого тела
- •Движение твердого тела вокруг неподвижной точки
- •Угловая скорость, формула Эйлера и движение твердого тела в общем случае
- •Сложное движение
- •Сложное движение точки. Основные понятия
- •Относительная производная
- •Теорема сложения скоростей в сложном движении точки
- •Теорема сложения ускорений в сложном движении точки
- •Теорема о сложении угловых скоростей в сложном движении твердого тела
- •ЧАСТЬ III. ДИНАМИКА
- •Уравнения движения и основные законы динамики механической системы
- •Принцип детерминированности и уравнение Ньютона
- •Инерциальные системы координат
- •Сила и масса. Второй и третий законы Ньютона
- •Законы сил
- •Две задачи динамики
- •Уравнения движения механической системы
- •Теорема об изменении главного вектора количества движения
- •Уравнение движения центра инерции
- •Кинетический момент относительно неподвижной точки и теорема о его изменении
- •Движение точки в центральном поле сил
- •Изменение кинетического момента, вычисляемого относительно подвижного полюса
- •Работа силы и изменение кинетической энергии материальной точки
- •Условия потенциальности силового поля
- •Кинетическая энергия системы и теорема Кенига
- •Теорема об изменении кинетической энергии системы
- •Движение точки в центральном поле сил (продолжение)
- •Динамика твердого тела
- •Масса и плотность. Геометрия масс
- •Основные законы динамики твердого тела
- •Вращение твердого тела вокруг неподвижной точки
- •Уравнения движения свободного твердого тела
- •Динамика точки с переменной массой
- •Уравнение Мещерского
- •Две задачи Циолковского
- •ЧАСТЬ IV. АНАЛИТИЧЕСКАЯ ДИНАМИКА
- •Общее уравнение механики
- •Cвязи, реакции. Обобщенные кооpдинаты
- •Изохронные ваpиации
- •Идеальные связи. Общее уравнение механики и принцип возможных перемещений
- •Общее уравнение механики в лагранжевых координатах
- •Уравнения Лагранжа
- •Уpавнения Лагpанжа II pода, их инваpиантность
- •Разрешимость уравнений Лагранжа II рода относительно старших производных
- •Обобщенный потенциал и уравнения Лагранжа II рода
- •Уравнения Лагранжа I рода и реакции идеальных связей
- •Канонические уравнения механики
- •Вывод канонических уравнений
- •Первые интегралы канонических уравнений
- •Метод Якоби решения канонических уравнений
- •Решение задачи о движении точки в центральном поле методом Якоби
- •Вариационные принципы механики
- •Дифференциальный принцип Даламбера-Лагранжа в декартовых переменных
- •Дифференциальный принцип Даламбера-Лагранжа в канонических переменных
- •Функционал и функция действия
- •Интегральный принцип наименьшего действия при изохронном варьировании (Принцип Гамильтона)
- •ЧАСТЬ V. ДОПОЛНИТЕЛЬНЫЕ СВЕДЕНИЯ
- •Структуры и пространства
- •Группы, кольца, поля
- •Векторные пространства
- •Метрические пространства
- •Банаховы, гильбертовы и евклидовы пространства
- •Пространства Rn
- •Тензоры
- •Сопряженные пространства
- •Два определения тензора и их эквивалентность
- •Примеры тензоров
- •Алгебраические операции над тензорами
- •ЧАСТЬ VI. МАТЕРИАЛЫ ДЛЯ ЗАНЯТИЙ
- •Проверочные вопросы
- •Экзаменационные вопросы
- •Тесты
- •Упражнения
- •Литература
- •Предметный указатель
- •Оглавление
Перейти к оглавлению на странице: 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) можно |
разрешить |
относительно |
величин |
||||||
q˙1, . . . , q˙s : |
|
|
|
|
|
|
|
|
|
q˙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 q˙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 ∂∂q˙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
Преобразование Дирака
Рассмотрим систему обыкновенных дифференциальных уравнений
q˙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 |
|
|
|||||
|
|
|
|
|
|
|
||
q˙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 |
q˙i + ∂ pi |
p˙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 |
piq˙i(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 |
|
q˙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 |
q˙k + ∂ pk p˙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