u_lectures
.pdfТогда обобщенная потенциальная сила равна
|
П |
|
|
|
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) в форму одночлена, получим
Подставим решение (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.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 в каждом главном колебании изменяются по гармоническому закону, имея одинаковые частоты и фазы: