Баев Теория колебаниы 2015
.pdfэтом случае величины |
, |
|
смогут принимать |
|
произвольные значения, |
иначе говоря, они будут независимыми. |
|||
|
π |
a 1,2, … , %, |
|
Тогда, рассматривая линейные формы (1.78), (1.79) как систему линейных алгебраических уравнений относительно обобщенных
скоростей и решая ее, получим следующие выражения: |
|
|
1.80 |
||||||||
|
|
b |
|
π b |
, |
$ 1,2, … , #. |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
|
|
|
|
|
|
π |
|
|
рассматривать |
a 1,2, … , %, |
|
|
|
|
|
|||||
Величины |
, |
|
|
|
|
и |
интегралы от |
них |
|
|
, будем |
|
как |
новые |
обобщенные скорости |
и |
координаты |
соответственно, только теперь, в отличие от исходных обобщенных скоростей и координат, число независимых обобщенных скоростей равно числу независимых обобщенных координат.
Запишем выражение для работы системы на виртуальных перемещениях:
|
|
|
|
|
|
|
|
|
|
||
δ3 5 δ @ |
4 |
A δ . |
|||
|
|
|
|
|
|
|
|
|
|
|
С другой стороны, как следует из выражений (1.80),
δ |
b |
|
δ π |
, $ 1,2, … , #. |
|
|
|
|
|
|
|
|
|
|
Объединив выражения (1.81) и (1.82), найдем, что
δ3 Π δπ ,
где
1.81
1.82
1.83
31
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Π |
|
|
|
1.84 |
||
|
b 5 b |
4 |
. |
Величины Π , a 1,2, … , %, есть не что πиное,. как обобщенные силы в новой системе обобщенных координат
Введя выражения (1.80) в выражения (1.75), преобразуем последние к виду
|
|
|
Q , π Q , ; |
1.85 |
||||||||
|
|
|
|
|
|
|
|
|
|
|||
из которых, в свою очередь, следует, что |
|
1.86 |
||||||||||
|
|
|
δ Q , δπ ; |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
… , |
ν 1,2, … , ; |
1.87 |
|||||
Q , πL |
||||||||||||
L |
|
|
|
|
|
|
|
|
|
|
|
|
Далее будем |
|
|
πL a 1,2, … , %. |
|
|
|
||||||
в правых частях выражений (1.87) выделены лишь члены, содержащие |
||||||||||||
вторые производные |
, |
|
|
|
|
|
|
|
|
|||
|
|
исходить из следующего очевидного выражения |
||||||||||
|
|
|
|
|
|
|
|
|
ν 1,2, … , . |
1.88 |
||
|
|
4 # δ 0, |
||||||||||
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Если в него ввести выражения (1.86), то оно преобразуется к виду |
||||||||||||
|
|
|
@Π |
|
|
# Q A δπ |
0. |
1.89 |
||||
|
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
32
Так как вариации δπ независимы, то равенства (1.89) возможны, если
|
# Q Π , a 1,2, … , %. |
1.90 |
|
L |
|
|
|
Оригинальность идеи Аппеля состоит в том, что он ввел понятие энергии ускорений
6 |
2 |
|
# |
6 , , π, πL. |
|
1 |
|
|
|
|
|
|
L |
|
Поскольку, как следует из выражений (1.87), |
||||
Q , |
1,2, … , ; a 1,2, … , %, |
|||
|
eL |
|
|
|
|
L |
|
|
то выражения (1.90) можно переписать в следующем виде: |
||
6 |
|
, a 1,2, … , %, |
πL Π |
||
|
|
|
1.91
1.92
1.93
которые и являются уравнениями Аппеля.
Решения уравнений (1.93) совместно с условиями (1.78) дают полное описание движения неголономной системы. В частности, при переходе к новым переменным, если это удобно, частично могут быть сохранены исходные обобщенные координаты вместе с их обобщенными скоростями.
Уравнения Аппеля являются второй разновидностью математических моделей, и если первая разновидность уравнения Лагранжа описывает движения голономных систем, то вторая разновидность математических моделей уравнения Аппеля описывает движения неголономных систем.
В качестве примера применения уравнений Аппеля решим следующую задачу: построить математическую модель системы, изображенной на рис. 1.4 и представляющей две точечные единичные массы М1 и М2, соединенные жестким стержнем длиной l c
33
пренебрежимо малой массой; система может двигаться только в вертикальной плоскости и только так, что скорость середины стержня всегда направлена вдоль стержня.
Решение задачи. У плоской системы (системы, движущейся в плоскости) с двумя точечными массами и одной геометрической связью должно быть три независимых обобщенных координаты.
y |
|
|
2М2
2
,
φ
М1
x
Рис. 1.4. Система, описываемая уравнениями Аппеля
В качестве обобщенных возьмем координаты центра стержня (x,y)
и угол его наклона φ к оси абсцисс, тогда |
|
||
|
B |
1.94 |
|
H |
|
H 2 cosφ; |
|
|
|
B |
1.95 |
I |
I 2 sinφ; |
||
|
B |
1.96 |
|
H |
|
H 2 cosφ; |
|
|
|
B |
1.97 |
I |
I 2 sinφ. |
Связь, наложенная на скорость центра стержня с координатами (x,y), очевидно является дифференциальной.
34
Выведем уравнение этой связи. |
Согласно условию задачи, если |
P |
||
скорость середины стержня, то |
|
|||
|
H Pcosφ, |
1.98 |
||
или |
I Psinφ |
1.99 |
||
H |
I |
|
|
|
|
1.100 |
|||
|
cosφ |
sinφ. |
Выражение (1.100) представляет записанную в аналитическом виде дифференциальную связь, которая не интегрируется (не сводится к конечной или геометрической связи), поскольку не является полным дифференциалом. Следовательно, эта система неголономная, и ее движение не может быть описано уравнениями Лагранжа.
Перейдем к выводу уравнений Аппеля, для которых требуется
построить энергию ускорений: |
|
|
|
|
1.101 |
|||
6 |
2 |
HL IL HL IL . |
||||||
|
1 |
|
|
|
|
|
|
|
Объединив это выражение с выражениями (1.94) – (1.97), получим |
||||||||
следующее выражение для энергии ускорений: |
1.102 |
|||||||
6 |
HL IL |
4 |
φL |
φ . |
||||
|
|
|
|
B |
|
|
|
|
В соответствии с изложенной выше методикой, теперь требуется
ввести две новые обобщенные скорости, в качестве которых выберем |
|
величины π и φ, положив |
1.103 |
H πcosφ; |
|
I πsinφ, |
1.104 |
и тогда выражение для энергии ускорений преобразуется к виду
35
6 πL |
4 |
φL … , |
1.105 |
|
B |
|
|
где точками в скобках, как и выше, обозначена та часть выражения, которая не содержит ускорений и поэтому не представляет интереса.
выбрана исходная обобщенная скорость φ. О таком варианте вывода уравнений Аппеля с частичным сохранением исходных обобщенных
Как видим, в качестве одной из двух новых обобщенных скоростей
переменных выше говорилось, и в решении данной задачи этот
вариант реализуется. Из выражений (1.98), (1.99) следует, что роль |
|||||||||
|
π P. |
|
|
|
|
|
|
|
|
второй обобщенной скорости играет скорость движения середины |
|||||||||
стержня |
|
|
Как увидим ниже, такой выбор новых обобщенных |
||||||
скоростей оказался очень удачным. |
|
|
|
|
|||||
Запишем выражение для работы системы на виртуальных |
|||||||||
перемещениях: |
|
δ3 Πδπ fδφ. |
|
|
1.106) |
||||
|
|
|
|
|
|
||||
С другой стороны, учитывая, что |
|
|
1.107 |
||||||
имеем |
|
|
|
δI δπsinφ, |
|
|
|||
|
|
δ3 2DδI 2Dδπsinφ. |
|
1.108 |
|||||
|
|
π, φ |
(1.106) и (1.108) следует, что в новой системе |
||||||
Из выражений |
|||||||||
координат |
|
|
обобщенные |
силы, |
действующие |
на |
|||
рассматриваемую систему, описываются выражениями |
1.109 |
||||||||
|
|
|
Π 2Dsinφ; f 0, |
|
|
||||
и, таким образом, уравнения Аппеля для данной системы имеют вид |
|
||||||||
|
|
|
|
6 |
Π; |
|
|
1.110 |
|
|
|
|
|
πL |
|
|
|||
|
|
|
|
6 |
f. |
|
|
1.111 |
|
|
|
|
|
φL |
|
|
36
Если в них подставить найденные выше выражения для энергии ускорений и обобщенных сил, то окончательно они принимают
следующий вид: |
πL Dsinφ; |
1.112 |
|
||
|
φL 0. |
1.113 |
Выбор новых обобщенных координат оказался настолько удачным, что уравнения движения построенной математической модели оказались интегрируемыми.
Опуская промежуточные несложные выкладки, приведем
найденные в аналитическом виде решения уравнений движения |
||||||||
D |
|
γ |
D |
cosφh sinφ δ; |
1.114 |
|||
H 2α |
φ gα |
2α |
||||||
|
γ |
|
D |
cosφh cosφ ε; |
|
|||
I gα |
2α |
1.115 |
||||||
|
|
φ α β, |
|
|
1.116 |
|||
в которых константы |
|
интегрирования |
|
определяются |
||||
начальными условиями. Задача полностью |
решена. |
|
||||||
|
α, β, γ, δ, ε |
|
1.9.Канонические уравнения Гамильтона
В этом параграфе будет рассмотрена третья разновидность математических моделей систем – уравнения Гамильтона. Для вывода этих уравнений потребуется понятие гессиана. Гессианом функции X(x1,x2,...,xn) называется определитель det(aik), i,k = 1,2,…,n,
элементами которого являются вторые производные этой функции:
j . 1.117
H H
Воспользуемся теоремой Донкина, которую приведем здесь без доказательства (с доказательством теоремы можно ознакомиться в книге [3]). Хотя в данном случае теорема Донкина имеет и прикладное
37
значение, она представляет самостоятельный интерес в вопросах, связанных с преобразование координат.
Теорема Донкина. Пусть дана некоторая функция X(x1,x2,...,xn), гессиан которой отличен от нуля, и пусть имеется преобразование переменных, порождаемое этой функцией:
|
j |
, $ 1,2, … , %; |
1.118 |
I |
H |
тогда существует обратное преобразование, порождаемое функцией
Y(y1,y2,…,yn): |
k |
|
|
|
, $ 1,2, … , %; |
1.119 |
|
H |
I |
и при этом порождающая функция Y обратного преобразования связана с порождающей функцией прямого преобразования выражением
k H |
I j; |
1.120 |
|
|
|
если функция X содержит параметры α1, α2, …, αm X= X(x1,x2,...,xn, α1,α2, …, αm), то функция Y будет содержать те же параметры Y(y1,y2,…,yn,α1, α2, …, αm ) и
k |
j |
, l 1,2, … |
1.121 |
α |
α |
Как видно из формулировки теоремы, она важна сама по себе, так как дает простую связь (1.120) между функциями прямых и обратных
преобразований переменных. |
|
|
|
|
||
Переменные |
|
, i= |
1,2,…,n, |
через |
которые выражаются |
|
|
называют переменными Лагранжа. В 1834 году У. |
|||||
функции Лагранжа, , , |
|
|
|
|
||
Гамильтон предложил |
в |
качестве |
переменных, определяющих |
|||
|
|
|
величины |
, , m |
,i= 1,2,…,n, где |
|
состояние системы, |
использовать38 |
|
|
переменные |
|
получившие |
название |
обобщенных импульсов, |
||
выражениями |
|
|
|
|||
определяютсяm |
, |
|
9 |
$ 1,2, … , %. |
1.122 |
|
|
|
|
||||
|
|
m |
, |
|||
Переменные |
|
, , m , i= |
1,2,…,n, |
называются |
переменными |
Гамильтона, пространство обобщенных координат и обобщенных импульсов получило, , название фазового пространства, а пространство переменных i= 1,2,…,n, называют расширенным фазовым пространством. Система в фазовом пространстве изображается точкой, а движение системы изображается линией, описываемой этой точкой, которую называют фазовой траекторией системы.
Забегая вперед, отметим, что понятие фазового пространства играет важную роль в теории колебаний и в ее многочисленных приложениях, так как форма фазовых траекторий системы дает обширную информацию о ее свойствах.
Перейдем теперь непосредственно к выводу канонических уравнений Гамильтона, для чего перепишем уравнения Лагранжа в следующем виде:
|
|
|
|
|
|
|
|
9 |
|
|
|
|
|
|
|
L … 0, $ 1,2, … , %; |
1.123 |
|
|
|
|
здесь точками обозначены члены уравнений, не содержащие обобщенных ускорений.
Поскольку уравнения движения в механике должны разрешаться относительно ускорений, т.е. должны приводиться к виду
L n , , , $ 1,2, … , %, |
|
|
1.124 |
то должен быть отличен от нуля определитель det( |
|
9/ 1 0, |
|
|
|
|
|
являющийся гессианом функции Лагранжа. |
|
|
|
39 |
|
|
|
Воспользуемся теоремой Донкина, а именно: будем считать обобщенные скорости и импульсы исходными и новыми переменными соответственно, а функцию Лагранжа – порождающей функцией прямых преобразований (1.118). Тогда, согласно теореме,
должна |
существовать |
порождающая |
функция |
обратных |
||||||||
преобразований, которую обозначим символом Н, |
|
|||||||||||
|
|
|
o |
|
$ 1,2, … , %, |
|
1.125 |
|||||
|
|
m , |
|
|
||||||||
и которая определяется выражением |
|
|
|
1.126 |
||||||||
|
|
|
|
o |
m 9. |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Функция |
получила название функции Гамильтона. Поскольку |
|||||||||||
|
координаты играют роль параметров, то |
|
||||||||||
обобщенные o |
|
9 |
o |
|
$ 1,2, … , %. |
1.127 |
||||||
|
|
|
|
|
, |
|
||||||
С помощью уравнений Лагранжа и равенств (1.127) получаем |
||||||||||||
следующую цепочку равенств: |
|
|
o |
|
|
|
||||||
&m |
|
|
& 9 |
|
9 |
, |
$ 1,2, … , %. |
1.128 |
||||
& |
m |
! |
|
|
|
|||||||
|
|
& |
|
|
|
|
|
|
||||
В итоге приходим к следующей системе из 2% уравнений вида |
||||||||||||
|
|
|
|
|
o |
; |
|
|
|
1.129 |
||
|
|
|
|
m |
|
|
|
|
||||
|
|
|
|
|
o |
, |
$ 1,2, … , %, |
1.130 |
||||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
m |
|
|
|
|
|
|
которые называются каноническими уравнениями Гамильтона.
40