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

u_lectures

.pdf
Скачиваний:
9
Добавлен:
29.02.2016
Размер:
1.62 Mб
Скачать

Тогда обобщенная потенциальная сила равна

 

П

 

 

 

1

S

S

 

 

 

 

 

QmП

=

 

cij

qi qj

=

 

 

qm

 

 

 

 

 

 

qm

2 i 1 j 1

 

 

 

 

 

 

S

S

 

 

 

 

S S

 

 

 

S

(24.12)

 

 

 

 

 

 

 

 

 

=

1 cijδmi qj

 

1 cijδmjqi

cmjqj

(m 1, 2, ..., s),

 

2 i 1 j 1

 

 

 

 

2 i 1 j 1

 

 

j 1

 

где использовалось (24.7) и свойство симметрии обобщенных коэффициентов жесткости cij cji .

Рассмотрим определение обобщенной силы сопротивления QmR . Пусть

на точки механической системы действуют силы линейно-вязкого сопротивления, т. е. сопротивления пропорционального первым степеням их скоростей:

Rк bкVк

(к 1, 2, ..., n) ,

(24.13)

что соответствует малым скоростям.

Согласно определению (23.4) для обобщенной силы найдем выражение для обобщенной силы сопротивления QmR :

n

r

 

r

 

 

 

n

r

 

r

 

 

QmR FкR

 

к

bкVк

 

 

 

к

=

qm

qm

к 1

 

 

 

к 1

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

V

2

 

 

 

 

r

 

 

 

 

 

 

 

к

 

 

n

 

 

V

 

 

n

 

 

 

 

2

 

= bкVк

 

 

к

=

bк

 

 

 

=

 

 

q&

&

 

 

к 1

 

 

 

 

 

к 1

 

 

 

 

m

(24.14)

 

 

 

 

 

qm

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

2

 

 

R

 

 

 

=

 

 

 

bкVк

 

 

 

 

 

 

,

 

 

q&

2

q&

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

к

1

 

 

 

m

 

 

 

где введена диссипативная функция Рэлея R 1 n bкVк2.

2 к 1

Для склярономной системы в потенциальном стационарном поле диссипативная функция Рэлея может быть выражена через полную энергию Е = Т + П механической системы в виде равенства

dEdt 2R ,

которое означает, что удвоенная функция Рэлея равна скорости убывания полной энергии в диссипативной системе.

По аналогии с кинетической энергией Т системы выразив функцию Рэлея R через обобщенные скорости, разложив ее в ряд в окрестности положения равновесия и отбросив члены выше второго порядка малости, получим

 

R

1 S

S

b

q&&q

 

 

 

(24.15)

 

 

ij 0

 

i

j .

 

 

2 i 1

j 1

 

 

 

 

 

Постоянные

коэффициенты

 

bij

 

0

называются

приведенными

 

 

 

 

 

 

 

 

коэффициентами сопротивления, которые ниже будем обозначать bij .

Подставляя (24.15) в (24.14), получим выражение для обобщенной силы сопротивления QmR , входящее в линеаризованное уравнение Лагранжа

(24.1), аналогичное (24.6)

QR R b q&

 

 

 

 

S

 

m

 

 

mj j .

(24.16)

 

&

 

 

qm

j 1

 

Важную категорию

образуют возмущающие силы

внешнего

происхождения, описываемые заданными функциями времени и не зависящие от движения системы.

Обобщенная возмущающая сила Qm (t) определяется по формуле

Q (t)

δAm

(t)

(24.17)

m

δq

,

 

m

 

где δAm (t) возможная работа возмущающих сил системы на возможном перемещении системы при δqm 0 .

В этом параграфе при определении обобщенных возмущающих сил Qm (t) ограничимся случаем, когда на механическую систему действуют

только гармонические возмущения. Тогда обобщенные возмущающие силы изменяются во времени по закону

Qm (t) Hm sin( pt δm )

(m 1, 2, ..., s)

(24.18)

т. е. силы Qm (t) имеют одинаковые частоты р, но различные достаточно малые амплитуды Hm и начальные фазы δm .

С учетом (24.8), (24.9), (24.12), (24.16) и (24.18) линеаризованное уравнение Лагранжа (24.1) принимает вид

S

a

 

q&

S

 

S

 

sin( pt δ

 

)

(m 1, 2, ..., s)

 

mj

c q

 

b q& H

m

m

 

j

mj j

 

mj j

 

 

 

j 1

 

 

 

j 1

 

j 1

 

 

 

 

 

или

S

a

 

q& b

q& c

 

q

 

H

 

sin( pt δ

 

)

(m 1, 2, ..., s)

 

 

mj

mj

j

m

m

. (24.19)

 

j mj

j

 

 

 

 

 

j 1

Линейные дифференциальные уравнения с постоянными коэффициентами (24.19) можно представить в векторно-матричной форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ &

ˆ &

ˆ

ˆ

 

(24.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aq

Bq

Cq

H ,

 

где введены обозначения: q

 

 

 

qj

 

 

 

 

вектор-столбец обобщенных координат,

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

aij

 

 

 

 

квадратная

матрица

обобщенных

коэффициентов

инерции,

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

bij

 

квадратная матрица приведенных коэффициентов сопротивления,

 

 

 

 

 

B

 

 

 

 

 

 

 

 

ˆ

 

 

 

cij

 

 

квадратная

матрица

обобщенных

коэффициентов

жесткости,

 

 

 

C

 

 

 

 

ˆ

Hm sin( pt δm )

вектор-столбец возмущающих сил.

H

В зависимости от действующих на механическую систему сил различают три вида колебательных движений:

1)собственные (гармонические) колебания, происходящие под действием только потенциальных (восстанавливающих) сил,

2)затухающие колебания, происходящие под действием потенциальных сил и сил линейно-вязкого сопротивления среды,

3)вынужденные колебания, когда кроме потенциальных сил и сил сопротивления среды, действуют возмущающие силы, зависящие от времени.

Свободные колебания механической системы с одной степенью свободы

При изучении данного вида колебаний будем исходить из уравнения Лагранжа второго рода (24.1), которое при S = 1 имеет вид

d

(

T )

T

 

П

QR .

(24.21)

dt

q

q

 

q&

 

 

 

В этом случае кинетическая энергия Т (24.5) системы определяется зависимостью

T12 аq&2 ,

апотенциальная энергия (24.11)

П12 сq2 ,

где а и с соответственно обобщенная масса системы и обобщенный коэффициент жесткости.

Обобщенная сила вязкого трения

QR Rq&;

где диссипативная функция рассеяния Рэлея R 12 bq&2 ; b обобщенный

коэффициент вязкости.

Вычислив соответствующие производные, входящие в уравнение (24.1), получим уравнение (24.19) для S = 1 и в отсутствии возмущающих сил:

&

&

cq

 

0 .

(24.22)

aq

bq

 

Уравнение (24.22) называется дифференциальным уравнением свободных колебаний при линейно-вязком сопротивлении.

Если на механическую систему не действуют силы линейно-вязкого сопротивления, то уравнение (24.22) соответствует дифференциальному уравнению свободных колебаний системы:

 

&

cq

 

 

0

 

 

 

(24.23)

или

aq

 

 

 

 

 

 

 

 

&

k

2

q

 

 

 

 

 

 

 

 

0 ,

 

 

 

 

 

q

 

 

 

 

 

 

где k

c / a собственная

(круговая)

частота

свободных

колебаний

системы.

Решение уравнения

 

 

(24.23)

при

начальных

условиях

q(0) q0 ; q&(0) q&0 имеет вид:

q

q&

sin kt q

cos kt.

 

0

(24.24)

 

k

0

 

 

 

 

 

После несложных преобразований, решение (24.24) можно записать более наглядной форме:

q A sin(kt ).

(24.25)

Здесь амплитуда свободных колебаний А определяется выражением

A q02 ( q&k0 )2 ;

а начальная фаза колебаний

arctg kqq&0 .

0

Зависимости (24.24) и (24.25) описывают свободные гармонические колебания механической системы около положения устойчивого равновесия.

Величина Т, определяемая выражением

T

2

2

a

,

 

k

 

c

 

является периодом свободных колебаний.

Теперь вернемся к уравнению (24.22). Разделим его на коэффициент а,

обозначим отношения b / a 2h,

c / a k 2 и запишем уравнение (24.22) в

форме

 

&

&

k

2

x

 

0 .

(24.26)

x

2hx

 

 

В зависимости от соотношений между k и h общее решение уравнения (24.26) имеет разный характер. Если сопротивление мало (k > h), т. е.

значение коэффициента вязкости b будет удовлетворять условию b 2 ac,

то механическая система будет совершать затухающие колебания около положения устойчивого равновесия q 0 . В этом случае решение уравнения

(24.26) следует искать в виде

 

 

q e ht (c1 sin k1t c2 cos k1t),

(24.27)

где k1

k2 h2

круговая

частота

затухающих колебаний.

Постоянные

интегрирования с1 и с2 с учетом начальных

условий движения q(0) q0 ,

q&(0) q&

определяются выражениями

 

 

 

0

 

 

 

 

 

c

q& hq

c

q .

 

 

 

0

0

;

 

 

 

 

 

 

 

 

1

 

k1

2

0

 

 

 

 

 

 

 

 

После преобразования решения (24.27) в форму одночлена, получим

q Be ht sin(k1t ).

(24.28)

Колебания, происходящие по закону (24.8), называются затухающими колебаниями механической системы. Коэффициент В в (24.28) равен

&

 

hq0 )

2

 

 

 

B

(q0

 

 

q2 ,

 

 

 

 

 

k2 h2

 

 

0

 

 

 

 

 

а начальная фаза

 

 

 

 

 

 

 

 

q

 

k2

h2

 

 

arctg

0

 

 

 

 

.

 

hq

 

 

 

 

q&2

 

 

 

 

 

0

 

0

 

 

Уравнение (24.28) показывает, что в рассматриваемом случае движение системы имеет гармонический характер с убывающей во времени по экспоненциальному закону величиной отклонения от положения равновесия.

Эту величину Be ht называют амплитудой затухающих колебаний. Количественной мерой быстроты затухания (убывания амплитуды)

колебаний механической системы является логарифмический декремент, который определяется выражением

hT

h 2

 

b 2

.

(24.29)

k2 h2

4ac b2

 

 

 

 

Для наглядного представления о характере колебаний механической системы удобно использовать понятие фазовой плоскости. Фазовая плоскость представляет собой плоскость, образованную прямоугольной системой координат, в которой по осям откладываются обобщенная координата и обобщенная скорость. В этом случае любое состояние механической системы определяется точкой фазовой плоскости с координатами (q, q ), которая называется изображающей точкой. При

движении механической системы эта точка описывает на фазовой плоскости кривую. Геометрическое место положений изображающей точки на фазовой плоскости называется фазовой траекторией.

Совокупность фазовых траекторий для всех возможных случаев движения системы называется фазовой диаграммой. Для случая свободных гармонических колебаний системы, ее фазовая траектория имеет форму эллипса; для затухающих колебаний фазовая траектория представляет собой логарифмическую спираль, накручивающуюся на начало координат, которое называется устойчивым фокусом.

ЛЕКЦИЯ 25 Свободные колебания механической системы

с двумя степенями свободы

Уравнения Лагранжа второго рода (24.1) для системы с двумя степенями свободы, находящейся под действием только потенциальных (восстанавливающих) сил, запишутся в виде

d

 

T

T

 

П

 

 

 

 

 

(

 

)

 

 

 

q

;

 

dt

q&

q

 

 

 

 

1

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(25.1)

 

d

(

T

)

 

T

 

 

П

.

 

 

 

 

 

 

q

 

 

dt

 

q&

 

 

 

q

 

 

 

 

 

2

 

 

2

 

2

 

 

Согласно (24.5) и (24.11) кинетическая и потенциальная энергии системы являются квадратичными формами обобщенных скоростей и обобщенных координат соответственно и определяются выражениями:

T a q& 2a q& q& a q&,

 

2

 

 

 

 

2

11

1

12

1

2

22

2

П c q2

2c q q c q2.

11

1

12

1

2

22

2

Здесь коэффициенты а11, а12, а22 и с11, с12, с22 по-прежнему имеют смысл обобщенных масс и обобщенных коэффициентов жесткости.

Вычислим производные от кинетической энергии и потенциальной энергии, входящие в (25.1), и подставим найденные выражения в эти уравнения. После не сложных преобразований получим уравнения (24.19) в отсутствие сопротивления и возмущающих сил:

a

 

q& a

 

q& c

 

q

c q

0;

 

 

11

1

12

2 11

 

1

12

2

 

(25.2)

a

21

q& a

22

q& c

21

q

c

q

0.

 

 

1

 

2

1

22

2

 

 

Уравнения (25.2) называются дифференциальными уравнениями свободных колебаний механической системы с двумя степенями свободы.

В случае малых колебаний механической системы около положения её устойчивого состояния равновесия ( q1 q2 0 ) решение уравнений (25.2) будем искать в виде

q1

Asin(kt ),

(25.3)

q2 Bsin(kt ),

 

где А, В амплитуды колебаний, k круговая частота колебаний, α начальная фаза колебаний являются неизвестными постоянными величинами.

Подставим решение (25.3) в исходные дифференциальные уравнения движения (25.2) и сгруппируем слагаемые:

 

a11k

2

)A (с12

 

a12k

2

)B

 

 

 

 

(с11

 

 

 

 

sin(kt ) 0;

 

 

a21k

2

)A (с22

a22k

2

)B

 

 

) 0.

 

(с21

 

 

 

 

sin(kt

 

В общем случае

sin(kt ) 0 .

Следовательно,

должны

выполняться

равенства:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

(с11 a11k

)A (с12 a12k

)B 0;

 

 

 

 

 

 

 

a k2 )A

(c

a

 

k2 )B 0.

(25.4)

 

(с

22

 

 

21

 

 

21

 

 

 

22

 

 

 

 

 

 

 

 

 

 

Эти алгебраические линейные однородные уравнения относительно А и В должны иметь ненулевое решение ( q1 q2 0 соответствует покою, а не движению). Поэтому определитель системы (25.4) должен равняться нулю:

 

c

 

a

k2

c

 

a

 

k2

0 .

 

 

 

11

11

2

12

12

 

k

2

 

 

 

c

21

a

k

 

c

22

a

22

 

 

 

 

 

 

 

21

 

 

 

 

 

 

 

 

Раскрывая его, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c11 a11k2 )(c22

a22k2 ) (c12

a12k2 )2

0.

(25.5)

Уравнение (25.5) называется уравнением частот, или вековым уравнением.

Корни k12 и k22 этого уравнения вещественные и положительные. Это доказывается математически, но может быть обосновано и тем, что если

k1 k12 и k2 k22 не будут вещественными, то система дифференциальных

уравнений (25.2) не будет иметь решения вида (25.3), что для механической системы, находящейся в положении устойчивого равновесия не возможно, так как после малых возмущений она должна двигаться около положения q1 q2 0 .

Каждому корню k1 и k2 будет отвечать свое частное решение (25.3), причем каждой частоте k1 и k2 соответствуют свои значения А, В и .

Поэтому, определив из (25.5) k1 и k2 , получим две совокупности частных решений вида (25.3):

q(1)

A sin(k t ),

q(1)

B sin(k t ),

 

1

1

1

1

 

 

2

1

1

1

 

 

(25.6)

q(2)

A

sin(k

t

2

),

q(2)

B

sin(k

t

2

).

1

2

2

 

 

2

2

2

 

 

 

Из уравнений (25.4) следует, что для каждой частоты ki амплитуды Ai иBi ( i 1, 2 ) связаны друг с другом. Для записи этой связи можно воспользоваться любым из уравнений (25.4):

 

 

B

 

 

c a k2

 

 

c

21

a k2

 

,

 

 

 

 

1

 

11

 

11 1

 

 

 

 

 

 

21 1

 

 

 

 

A

 

 

a k2

c

 

a

 

 

 

 

 

 

 

c

 

 

 

22

22

k2

 

 

1

 

 

 

 

1

 

 

12

 

12

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

B

 

 

 

c

a k2

 

 

 

c

 

a k2

 

 

 

 

 

(25.7)

 

 

 

 

 

 

 

 

21

 

 

.

 

 

 

 

2

 

 

 

11

 

11 2

 

 

 

 

 

 

 

 

21 2

 

 

 

 

A

 

 

c

a k2

 

 

c

 

a

 

 

 

 

 

 

 

 

 

 

 

 

22

22

k2

 

 

2

 

 

 

 

2

 

 

12

 

12

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

т. е. B1 1A1 , а B2 2 A2 . Тогда частные решения (25.6) принимают вид:

q(1)

A sin(k t ),

 

 

q(1)

 

A sin(k t ),

1

1

 

 

1

 

1

 

 

 

2

 

 

 

 

1

1

 

 

1

 

1

 

 

q(2)

A sin(k

t

 

 

 

 

q(2)

 

 

 

 

 

A sin(k

t

 

(25.8)

2

),

 

 

 

 

2

).

1

2

 

 

2

 

 

 

 

2

 

 

 

 

 

2

 

2

 

 

2

 

 

 

Колебания, определяемые уравнениями (25.8), называются главными

колебаниями, а их частоты k1

и k2 собственными частотами системы.

Колебание с частотой k1 (всегда меньшей) называют первым главным колебанием, а с частотой k2 вторым главным колебанием.

Числа 1 и 2 , определяющие отношения амплитуд (или обобщенных

координат q2 / q1 ) в каждом из этих колебаний, называются коэффициентами

формы. Эти коэффициенты показывают во сколько раз амплитуда соответствующего главного колебания одной из обобщенных координат больше (или меньше) амплитуды другой обобщенной координаты.

Поскольку частные решения (1.37) являются линейно независимыми, то общее решение системы дифференциальных уравнений (25.2) будет равно их линейной комбинацией:

q1 q11 q12 A1 sin(k1t 1) A2 sin(k2t 2 ),

q2 q12 q22 1A1 sin(k1t 1) 2 A2 sin(k2t 2 ).

Теперь найдем постоянные интегрирования А1, А2, 1 , и в общем виде запишем начальные условия движения системы

t 0,

q (0) q

,

q&(0) q&

,

 

1

 

10

 

 

1

10

. ;

 

q

(0)

q

 

,

q&(0)

q&

 

2

 

20

 

2

20

 

(25.9)

2 . Для этого

(25.10)

Вычислим обобщенные скорости, соответствующие обобщенным координатам (25.9):

q& k

A cos(k t ) k A

cos(k

t

),

 

 

1

1

1

1

1

2

2

 

2

 

2

 

 

(25.11)

q& k

 

A cos(k t ) k

 

A cos(k

t

 

).

2

1

1

1

1

1

2

2

2

2

 

2

 

 

Подставляя начальные условия (25.10) в (25.9) и (25.11), получаем:

q10 A1 sin 1 A2 sin 2 ,

q20 1A1 sin 1

2 A2 sin 2 ,

 

 

 

 

 

 

q&

k A cos

k

A

 

cos

,

 

 

 

 

 

 

10

 

1

1

 

 

 

1

 

2

2

 

 

 

 

 

2

 

 

 

 

 

 

 

q&

k

 

A cos k

 

2

A cos

.

 

 

 

 

20

 

1

1

1

 

 

 

1

2

 

 

 

2

 

 

 

2

 

 

 

 

 

Решая эти уравнения, найдем амплитуды и начальные фазы А1, А2, α1, α2:

A

 

 

 

 

1

 

 

 

(q

 

q

 

)2

 

(q&

 

q& )2

,

 

 

 

 

 

 

 

 

 

 

20

 

2 10

 

 

 

 

 

 

 

 

k2

 

 

1

 

 

2

 

 

20

 

2

10

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

A

 

 

 

 

1

 

 

(q q )2

 

(q&

q& )2

,

 

 

 

 

 

 

 

 

20

1 10

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

20

 

1

10

 

 

 

 

 

 

k2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(25.12)

 

 

 

 

 

 

 

 

 

 

 

 

q20

2q10

 

 

 

 

 

 

 

 

 

1

arctg

k1

 

 

,

 

 

 

 

 

 

q&

q&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

2

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q20

1q10

 

 

 

 

 

 

 

2

 

arctg k2

 

.

 

 

 

 

 

 

 

 

 

q&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

1 10

 

 

 

 

 

 

 

Проведенный анализ показывает, что движения системы с двумя степенями свободы около положения устойчивого равновесия слагаются из двух независимых колебаний, которые называются собственными формами

колебаний с соответствующей собственной частотой (k1 или k2). Коэффициенты форм 1 и 2 в совокупности описывают

конфигурацию системы при наибольшем ее отклонении от положения равновесия в процессе свободных колебаний с соответствующей собственной частотой.

Из (25.8) видно, что координаты q1 и q2 в каждом главном колебании изменяются по гармоническому закону, имея одинаковые частоты и фазы:

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]