Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Баев Теория колебаниы 2015

.pdf
Скачиваний:
22
Добавлен:
12.11.2022
Размер:
2.85 Mб
Скачать

этом случае величины

,

 

смогут принимать

произвольные значения,

иначе говоря, они будут независимыми.

 

π

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

φ gα

 

γ

 

D

cosφh cosφ ε;

 

I gα

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, α12, …, αm), то функция Y будет содержать те же параметры Y(y1,y2,…,yn1, α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