Баев Теория колебаниы 2015
.pdfОбобщенные координаты и импульсы уравнений Гамильтона (1.129), (1.130) называют канонически сопряженными.
Отметим два важных свойства функции Гамильтона, следующих непосредственно из одноименных уравнений. Первое свойство состоит в том, что ее полная и частная производные по времени
совпадают. Действительно, |
|
o |
|
o |
|
|
o |
|
||||
&o |
|
o |
|
|
|
|
|
|||||
& |
|
|
|
m ! |
|
|||||||
|
|
o |
m |
|
o |
|
||||||
|
|
o |
o |
|
o |
|
. |
1.131 |
||||
|
|
m |
m |
! |
|
|||||||
Второе свойство связано со склерономными системами, для которых справедливы равенства
|
|
9 |
1.132 |
o m |
9 9 2( 9 ( 6. |
||
|
|
|
|
Из них следует, что у склерономной системы функция Гамильтона является ее полной энергией. Эти свойства функции Гамильтона делают понятной принятую в аналитической механике терминологию. Система называется обобщенно-консервативной, если ее функция Гамильтона не зависит явно от времени.
Функция Гамильтона еще называется обобщенной полной энергией системы, а равенство H(q,p)=const, справедливое для обобщенно-консервативных систем, обобщенным интегралом энергии.
Наряду с уравнениями Лагранжа и Аппеля канонические уравнения Гамильтона представляют третью разновидность математических моделей систем. Они получили столь же широкое распространение в физике, как и уравнения Лагранжа, а в некоторых ее разделах им даже отдается предпочтение. Это, прежде всего, объясняется большой информативностью понятия фазового пространства, а также функции Гамильтона благодаря ее физическому
41
смыслу. Кроме того, уравнения Гамильтона удобны для их численного интегрирования и математического моделирования движения систем, поскольку имеют вид, к которому, как правило, должны приводиться дифференциальные уравнения в численных методах их интегрирования.
Недостатком канонических уравнений Гамильтона является то, что они справедливы только для натуральных систем, тогда как уравнения Лагранжа распространяются и на системы с не потенциальными силами.
В качестве примера рассмотрим простой маятник, функция Гамильтона которого имеет вид
o |
m |
|
#DBcosφ; |
1.133 |
2#B |
|
здесь # масса маятника; B – его длина; φ и p – обобщенные координата и импульс соответственно.
Поскольку функция Гамильтона маятника не зависит явно от времени, она является интегралом энергии H(q,p)=const, который, в свою очередь, определяет фазовые траектории маятника, представленные на рис. 1.5.
Каждая фазовая траектория определяется значением const. Совокупность всех фазовых траекторий системы называется ее фазовым портретом, так что изображенная на рис. 1.5 картинка является фазовым портретом маятника. Как видно из рисунка, фазовые траектории можно разбить на две группы. Первую группу образуют траектории, занимающие ограниченную область фазового пространства, а это означает, что система не может сколь угодно далеко удаляться от стартовой точки, колеблясь около положения, которое называют положением устойчивого равновесия, а само движение системы называют устойчивым. Причем если, как в данном случае, система является консервативной, она будет периодически возвращаться к своему стартовому положению, и поэтому фазовые траектории должны быть замкнутыми. Такое движение системы называют устойчивым.
42
p
|
0 |
φ |
Рис. 1.5. Фазовый траектории маятника
Вторая группа состоит из неограниченных в фазовом пространстве траекторий, когда система может сколь угодно далеко удаляться от исходного положения. В случае с маятником они соответствуют его вращению вокруг точки подвеса. Подобные траектории в фазовом пространстве соответствуют так называемому неустойчивому движению системы.
Толстой линией на рис. 1.5 выделена фазовая траектория, которая служит границей между областями устойчивого и неустойчивого движения. Подобная траектория проходит через точки неустойчивого равновесия, которыми в случае простого маятника являются точки, кратные π. Эта граничная фазовая траектория играет важную роль в теории колебаний и в некоторых работах, как по теории колебаний, так и в ее приложениях, получила название сепаратриса [4,5].
Если ограничиться малыми отклонениями маятника от его
положения равновесия, когда |
|
|
|
|
, то функцию |
Гамильтона |
||||
можно аппроксимировать |
выражениями вида |
|
|
|
||||||
|
m |
|\φ| q 1 |
|
|
|
|
||||
o r |
|
s |
#DB |
\φ |
|
, |
1.134 |
|||
2#B |
|
2 |
|
|||||||
|
|
|
|
43 |
|
|
|
|
|
|
точек |
φ 2%π, % 0, s1, s2, …, |
|
где знак «плюс» соответствует окрестностям точек устойчивого |
||
Фазовые портреты маятника вφ 2% 1 e, % 0, s1, s2, … |
||
равновесия, |
|
а знак «минус» – окрестностям |
неустойчивого равновесия,
окрестностях этих точек показаны на рис. 1.6. Сами эти точки являются особыми точками фазовой плоскости, так как в них не определено направление скорости
системы.
типа «центр». О том, что особая точка типа «центр» это точка устойчивого равновесия, говорит вид фазовых траекторий в ее окрестности, которыми в приближении (1.134) являются эллипсы.
Особая точка, изображенная на рис 1.6,а называется особой точкой
Особая точка на рис. 1.6,б называется особой точкой типа
«седло».Вблизи данной точки фазовые траектории в приближении (1.134) представляют гиперболы, а их вид свидетельствует о том, что особая точка типа «седло» это точка неустойчивого равновесия. На фазовой плоскости встречаются и другие разновидности особых точек.
|
p |
p |
|
φ |
|
а |
б |
Рис. 1.6. Фазовые портреты движения маятника около точек его устойчивого (а) и неустойчивого (б) равновесия
Из уже сказанного о фазовой плоскости следует, что она содержит важную информацию о системе, а именно: данные об областях устойчивости, их форме и размерах, координаты точек устойчивого и неустойчивого равновесия. По виду фазовых траекторий в области
44
устойчивости можно судить о характере сил действующих на систему:
если фазовые траектории замкнутые, то система |
|
консервативная, ее |
|||
полная энергия сохраняется, сворачивающиеся |
|
фазовые траектории |
|||
|
|
|
|
|
|
свидетельствуют о действии на систему диссипативных сил и т.п. |
|||||
В частности, поскольку особая точка типа «седло» |
|
это точка |
|||
равновесия системы, в которой она покоится, а |
сепаратриса проходит |
||||
|
|
|
|
|
|
через эту точку, то логично предположить, что и в любой другой точке сепаратрисы система должна покоиться. Из данного утверждения следует, что по мере удаления от центра области устойчивости и приближения к ее границе (сепаратрисе), т. е. по мере роста амплитуды колебаний маятника частота его колебаний будет уменьшаться, и на самой сепаратрисе она обратится в нуль.
Как уже об этом говорилось в предисловии, концепция фазового пространства представляет мощный инструмент для анализа математических моделей систем, что будет проиллюстрировано ниже на примерах решения сравнительно сложных задач. Здесь же ограничимся ссылкой на очень полезную монографию [6], посвященную приложению понятия фазового пространства к решению задач по динамике частиц.
1.10.Уравнения Рауса
Уравнения Рауса получаются, если в качестве основных переменных, определяющих пространство, в котором строится математическая модель системы, взять часть переменных Лагранжа и
часть переменных Гамильтона.
Получающийся при этом набор переменных t, qi, qα, , pα ; i = =1,…,m; α = m+1,…,n (здесь m произвольное фиксированное число,
меньшее n), называют переменными Рауса. Чтобы перейти к этим
переменным, необходимо выразить обобщенные скорости |
через |
||||
обобщенные импульсы p |
|
|
|
используются |
|
|
α |
|
(α = m+1,…,n), для чего |
|
|
выражения (1.122), определяющие эти импульсы. |
|
|
|||
Обобщенные скорости |
следует выбирать таким образом, чтобы |
||||
гессиан лагранжиана был |
отличен от нуля: |
|
|
||
|
|
|
|
|
|
45
det , 9 - 1 0. 1.135
Если условие (1.123) выполнено, то в соответствие с упоминавшейся выше теоремы Донкина вводится функция переменных (t, qi, qα, , pα), которая называется функцией Рауса и определяется выражением
|
uuu |
u |
1.136 |
|
|
m |
9; |
||
|
|
|
|
|
здесь черточка над символами означает, что отмеченные ею
величины выражены через pα. |
|
|
При этом переменные t, qi, qα, |
(i = 1,…,m; α |
= m+1,…,n) |
рассматриваются как параметры, |
для которых |
справедливы |
выражения (1.121): |
|
|
|
9 |
, |
|
9 $ 1, … , # ; |
|
|
|
|
|
|
|
|
|
|
9 |
α # 1, … , % ; |
|
|
|
|
|||
|
|
|
|
|
9 |
|
|
|
. |
||
1.137 1.138 1.139
|
|
|
|
|
|
|
В итоге для координат |
|
уравнения Лагранжа сохраняют свой вид, |
||||
а для координат |
преобразуются в уравнения Гамильтона: |
|
||||
|
|
|
0 $ 1, … , # ; |
1.140 |
||
|
& |
|
||||
|
|
|
|
|
||
|
& |
|
|
|
|
|
46
& |
|
|
|
&m |
|
α # 1, … , % . 1.141 |
& |
|
m |
, |
& |
|
Уравнения (1.140), (1.141) и представляют уравнения# Рауса, состоящие из m уравнений Лагранжа и 2(n уравнений Гамильтона, так что функция Рауса в первых уравнениях играет роль функции Лагранжа, а во второй группе уравнений – роль функции Гамильтона. Заметим, что уравнения Рауса – уже третья разновидность математических моделей систем.
1.11.Циклические координаты
Обобщенная координата qα называется циклической, если она не входит явно в функцию Лагранжа, так что
9 0. 1.142
Из уравнений Гамильтона, учитывая выражения (1.127), следует, что обобщенный импульс pα, соответствующий циклической координате qα, в процессе движения системы сохраняет постоянное значение pα= const= cα.
Допустим теперь, что у системы r циклических координат. Соответствующие им обобщенные импульсы имеют постоянные значения, определяемые начальными условиями, и функцию Гамильтона можно представить в следующем. виде: H = H(t,q, p1, p2,…, pm, cm+1, cm+2,…, cn), здесь n = m
Циклические обобщенные координаты как функции времени получаются непосредственно из уравнений Гамильтона (1.130) в квадратурах
|
|
o |
& v |
, α # 1, # 2, … , %. 1.143 |
; v |
||||
|
|
|
! |
|
|
|
|
|
|
47
Таким образом, при наличии у системы r циклических координат ее математическая модель упрощается, поскольку фактически она будет состоять из 2(n– уравнений Гамильтона. В частности, если удастся найти такую систему обобщенных координат, в которой все координаты будут циклическими, то описание движения системы получается автоматически в виде постоянных значений обобщенных импульсов и обобщенных координат в виде квадратур (1.143).
Поэтому при выборе системы обобщенных координат необходимо стремиться к максимально большему количеству циклических координат. Обычно циклические координаты появляются в связи с разного рода симметриями в пространственных структурах силовых полей, в которых движется система. Например, если в цилиндрических координатах силовое поле не зависит от азимута (так называемые азимутально-симметричные поля, часто встречающиеся в инженерной практике), то, выбрав его в качестве одной из обобщенных координат, получим постоянный азимутальный импульс, а количество канонических уравнений Гамильтона уменьшится на две
единицы. |
|
|
|
|
Как уже отмечалось выше, если у системы r циклических |
||
координат, то ее движение описывается |
|
уравнениями Гамильтона, |
|
где |
– С помощью уравнений Рауса математическую модель |
||
|
|
2# |
|
# % .
системы можно представить в виде системы из дифференциальных
#
уравнений второго порядка типа уравнений Лагранжа. Удобство этих уравнений состоит в том, что, если не включить циклические координаты в функцию Рауса, то для их решения отпадает необходимость предварительно определять функции, описывающие зависимость циклических координат от времени.
Действительно, заменив обобщенные импульсы pα циклических |
|||||||||
координат |
α # 1, # 2, … , % |
|
определяемыми |
начальными |
|||||
условиями |
|
|
$ 1, … , # |
, |
ведем |
функцию Рауса в |
виде |
||
( |
|
|
α |
||||||
R=R(t, qi, |
, |
сα) ( |
|
|
. |
Тогда |
уравнения |
Рауса |
для |
нециклических координат |
|
0 $ 1, … , # |
1.144 |
||||||
|
|
|
& |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
& |
|
|
|
|
|
|
48
не содержат циклические координаты, являясь в этом смысле по отношению к ним автономными, что позволяет, сначала решив эти уравнения, найти функцию Рауса, а уже потом – циклические координаты через квадратуры
|
|
& v |
, α # 1, # 2, … , %, 1.145 |
; v |
|||
|
|
! |
|
|
|
|
|
|
|
|
естественно, предварительно заменив в функции Рауса все величины qi и на найденные в результатеv $ 1,решения… , # . уравнений Рауса функции от 2m + 1 аргументов t, ci, !
В качестве примера рассмотрим движение заряженной частицы в электромагнитном поле, обладающем в цилиндрической системе координат азимутальной симметрией.
Азимутально-симметричное магнитное поле может быть описано двумя компонентами векторного потенциала Aφ и Az [1], так что лагранжиан заряженной частицы в соответствии с выражением (1.74) будет иметь вид
|
9 2 |
|
φ ] Q] 3 |
Q φ 3 Q6; |
1.146 |
||||||||
|
# |
|
|
|
|
|
|
|
" |
|
|||
здесь U(r,z,t) – потенциал электростатического поля. |
|
||||||||||||
Координата |
|
|
циклическая, поэтому сохраняется азимутальный |
||||||||||
импульс |
|
|
" |
9 |
|
|
φ 3 |
" |
const. |
1.147 |
|||
|
|
m |
φ # |
|
|||||||||
|
φ |
|
|
|
|
|
|
|
|
|
|
||
Введем функцию Рауса, которая в данном случае будет иметь вид |
|||||||||||||
|
|
|
|
|
# |
|
|
|
|
1.148 |
|||
где |
9 φ φ 2 |
|
|
] Q] 3 Q6 f, |
|||||||||
|
|
|
|
|
|
1 g |
m" |
|
|
|
|
||
|
|
|
|
f |
|
Q3"h . |
1.149 |
||||||
|
|
|
|
|
|
2# |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
49 |
|
|
|
|
|
Ограничимся случаем, когда |
|
|
что соответствует отсутствию |
|||||||
|
азимутальной компоненты, |
|
|
|
В этом случае |
|||||
у магнитного поля |
|
3 |
0, |
|
|
|
" |
|
|
|
математической моделью |
движения |
заряженной частицы будут |
||||||||
|
|
В |
|
0. |
|
|||||
следующие уравнения Рауса: |
|
|
|
|
|
|
|
|
||
|
|
|
|
; |
|
|
|
1.150 |
||
|
|
|
|
|
|
|
|
|||
|
|
|
|
, |
|
|
|
1.151 |
||
где f Q6. |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
1.12.Интегралы уравнений движения
Функция f(t,q,p) называется интегралом уравнений движения системы, если она сохраняет постоянное значение на фазовых траекториях системы или, что то же самое, на решениях уравнений движения.
Например, у обобщенно-консервативной системы интегралом является функция Гамильтона, а у системы с циклическими координатами – соответствующие им обобщенные импульсы.
Хотя произвольная функция от интеграла уравнения движения также является интегралом, интерес представляют только независимые интегралы движения, т.е. такой набор функций, являющихся интегралами уравнений движения, в котором ни одна из функций не может быть выражена через другие функции из того же
набора. Если известна полная система |
интегралов, состоящая |
из |
|||||||||||||
|
, , m |
v |
интегралов |
движения, |
и m , |
|
|
выражения |
|||||||
2nнезависимых |
то, |
разрешив |
|||||||||||||
|
|
|
|
|
относительно величин |
|
|
получим |
|
функции, |
|||||
полностью |
|
|
, |
|
. |
|
|
, v , v , … , v , |
|||||||
содержащие |
2n |
констант, |
определяемых начальными |
условиями, |
|||||||||||
m |
m , v , v , … , v $ 1, … , % |
|
|
|
|
|
# |
|
|
||||||
|
|
|
описывающих движение системы |
|
|
|
|
|
|
||||||
|
Если известны только l независимых интегралов |
|
|
, l |
, то |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
системы, и чем |
|||
они дают лишь частичное представление о движении , … , |
|
N 2% |
|
||||||||||||
больше |
l, |
тем |
полнее это |
представление. |
Поэтому при анализе |
||||||||||
|
|
|
|
|
|
|
50 |
|
|
|
|
|
|
|
|
