振动力学第四章
大学物理 第4章机械振动(完全版)
注意:研究简谐振动时,坐标原点只能取在平衡位置。
平衡位置: k
F
外
0, 或 M外 0
o (原长) (平衡位置)
m
x
21
４．4简谐振动的能量 x =Acos( t+ ) =- Asin( t+ )
振动势能： E p 振动动能：Ek
1 2 1 kx
2
k
o (原长) (平衡位置)
m
o (原长) (平衡位置) x
o
x
xo (平衡位置)
(原长)
m
1 2 kx
2
Ep
弹
E p振
E p振 E p弹 1 2
1 2
kx
2
k ( xo x)
2
24
例题如图，有一光滑水平面上的弹簧振子，弹簧的 倔强系数k=24N/m, 物体的质量m=6kg, 静止在平衡位 置。设以一水平恒力F=10N向左作用于物体，使之由 平衡位置向左运动了s=0.05m, 此时撤去外力F。取物 体运动到左方最远处开始计时，求：(1)物体的运动方 程; (2)何处Ek=Ep? 解 (1) k
质点的简谐振动状态由下面两个物理量确定： x =Acos( t+ )
dx dt
A sin( t )
显然，它们由相位唯一确定。
5
五 .振动的超前与落后
设有两个同频率的谐振动： x1=A1cos( t+1) x2=A2cos( t+2) 振动x2超前x1(2 -1) ； >0, =0, 振动x2和x1同相 ； 相差 =2 -1 <0, 振动x2落后x1(2 -1) ； =, 振动x2和x1反相 。
振动力学［PDF］
第四章多自由度系统的振动4.1多自由度系统运动方程的建立4.2 耦合与坐标变换4.3 固有频率和主振型4.4振型矩阵、主坐标和正则坐标4.5 固有频率相等的情况4.6 固有频率为零的情况4.7 无阻尼系统对初始条件的响应4.8 无阻尼系统对任意激励的响应4.9 多自由度系统的阻尼4.10 有阻尼系统的响应4.11 一般粘性阻尼系统的响应一般粘性阻尼系统的响应i nj nj j j i j j i nj j j i ==•=••111i n j n j j j i j j i n j j j i Q q k q c q m =++∑∑∑==•=••111•••[][]{}[]{}{}Q q k q c q m =++⎭⎬⎫⎩⎨⎧•••nn n n n n 212222111211212222111211nn n n n n 212222111211212222111211nn n n n n 212222111211212222111211••••••n 11•••n 11n 11n 11inj nj j j i j j i nj j j i ==•=••111i n j n j j j i j j i n j j j i P x k x c x m =++∑∑∑==•=••111•••[][]{}[]{}{}P x k x c x m =++⎭⎬⎫⎩⎨⎧•••n 2121•i ••i1m 12m 23m 3ii •i i Q q Dq T =∂∂+∂∂•1m 12m 23m 32222)2221k +2222•2221⎟⎠⎞+•x c jjj j q W δδ11x δ11P 111x P δ1x δ22P 33P 2⎟⎠⎞21221212212111••••122121221211123323212332321222•••••2332321233232122233323332333••••33323332333•••••••••321333322221321321321333322221321321000000003213213333222213213213333222210•1•1θv ••2θv •1•=1θ•2•=2θ22θ−+mg l 22θl +k Oθ222yk+=•••[][]{}[]{}{}P x k x c x m =++⎭⎬⎫⎩⎨⎧•••••••••n n i j j i i 1111•••nn i j j i i 1111in n i j j i i 1111i j刚度影响系数k i j 若系统各自由度的广义速度和广义加速度为零，除j i i j i 。
大学物理 新 4 振动
v0 M
v m
x 解：(1).x处物块动力学方程
(m M ) d 2x dt 2
x (m M )g k( Mg x)
o
k
o ＝ mg kx
正确解：(m
M)
d2x dt 2
(m
M
)g
k（ m
M）g k
x＝
kx
k T 2 M m
22
mM
t 0 : A Acos
2
3
v0 A sin 0
3
66
振动（Vibration）
t 1: x1 0, v1 0
x
A
3
2
0
t 1: 0 Acos( )
3
32
v1
A
sin(
3ቤተ መጻሕፍቲ ባይዱ
)
0
(2). 2
t T
t
2
t 0.5 5 5
6 12
A 1.0
dt
1。运动方程中各物理量
(1)周期、频率、角频率
周期T：完成一全振动所需的时间
x = Acos(ω t + ) = A cos ω (t + T )+
一个周期后位移相等，所以
ω T = 2π 11
振动（Vibration）
数学式
T 2
弹簧振子T 2 m /k
单 摆T 2 l /g
复摆T 2 J / mgh
P
x
33
振动（Vibration）
M
P
x
34
振动（Vibration）
M
P
x
35
机械振动基础
第4章 机械振动基础4-1 图示两个弹簧的刚性系数分别为k 1 = 5 kN/m ，k 2 = 3 kN/m 。
物块重量m = 4 kg 。
求物体自由振动的周期。
解：根据单自由度系统自由振动的固有频率公式 mk =n ω 解出周期 nπ2ω=T图（a ）为两弹簧串联，其等效刚度 2121eq k k k k k +=所以 )(2121n k k m k k +=ω2121n)(π2π2k k k k m T +==ω代入数据得s 290.0300050003000)4(5000π2=⨯+=T图（b ）为两弹簧串联（情况同a ） 所以 T = 0.290 s图（c ）为两弹簧并联。
等效刚度 k eq = k 1 + k 2 所以 mk k 21n +=ω21nπ2π2k k mT +==ω代入数据得 T = 0.140 s图（d ）为两弹簧并联（情况实质上同（c ））。
所以 T = 0.140 s4-3 如图所示，质量m = 200 kg 的重物在吊索上以等速度v = 5 m/s 下降。
当下降时，由于吊索嵌入滑轮的夹子内，吊索的上端突然被夹住，吊索的刚度系数k = 400 kN/m 。
如不计吊索的重量，求此后重物振动时吊索中的最大张力。
解：依题意，吊索夹住后，重物作单自由度自由振动，设振幅为A ，刚夹住时，吊索处于平衡位置，以平衡位置为零势能点，当重物达到最低点时其速度v = 0。
根据机械能守恒，系统在平衡位置的动能与最低点的势能相等。
即 T max = V max 其中 2max 2v m T = ， 2max 21kA V =v km A =吊索中的最大张力 mk v mg kA mg F +=+=max 代入数据得 kN 7.461040020058.92003max =⋅⋅+⋅=F4-5 质量为m 的小车在斜面上自高度h 处滑下，而与缓冲器相碰，如图所示。
缓冲弹簧的刚性系数为k ，斜面倾角为θ。
理论力学 第四章
(m
(m
2J r
2
) x kxx 0 x
) kx 0 x
2J r
2
－－自由振动微分方程
系统的固有频率为
0
k r2 mr 2 2 J
§ 4－2 计算固有频率的能量法
如图所示无阻尼振动系统
当系统作自由振动时，运动规律为
x A sin(0t )
2
2
当圆柱体作微振动时， 可认为 sin
2
2
1 V mg ( R r ) 2 2
设系统作自由振动时θ的变化规律为 A sin(0 t )
3m 2 ( R r ) 2 0 A 2 则系统的最大动能 Tmax 4 1 2 系统的最大势能 Vmax mg ( R r ) A 2 由机械守恒定律 有 Tmax Vmax
2 0
k m
0
k m
0 只与表征系统本身特性的质量m和刚度k有关
而与运动的初始条件无关
它是振动系统固有的特性
所以称为固有角（圆）频率（一般也称固有频率） m=P/g
k P / st
0
g
0
k m
st
（2）振幅与初相角
x A sin( 0 t )
速度为
dx v 0 A cos(0t ) dt
在瞬时t 物块的动能为
1 2 1 2 T mv m0 A2 cos2 (0t ) 2 2
若选平衡位置为零势能点，有
1 2 V k[( x st ) 2 st ] Px 2
k st P
1 2 1 2 V kx kA sin 2 (0 t ) 2 2
振动力学基础
k1 = k2 = 2k0
k0 4k0 ϖ= =ϖ0 ϖ = = 2ϖ 0 m m
例. 质量为m的比重计，放在密度为 ρ 的液体中。已知比 重计圆管的直径为d。试证明，比重计推动后，在竖直方向 的运动为简谐振动，并计算周期。 解： 取平衡位置为坐标原点 平衡时： mg 浮力：
−F =0
F
F = ρ Vg
o
2
2
d πρg ω= 2 m
x
x
k的轻弹簧、一 例.如图所示，振动系统由一倔强系数为 如图所示，振动系统由一倔强系数为k 半径为 R、转动惯量为 J的定滑轮和一质量为 m的物体所 半径为R 转动惯量为J 的定滑轮和一质量为m 组成。使物体略偏离平衡位置后放手，任其振动，试证 物体作简谐振动 . 物体作简谐振动. 解：取位移轴 ox ， 解：取位移轴ox ox， m的平衡位 原点在 原点在m m在平衡位置 置。 置。m 时，设弹簧伸长量 : 为∆l，则有 ，则有:
x1
x2
1 1 1 2 2 2 m1v1 + m2 v2 + kx = c 2 2 2
m1v1 + m2v2 = 0 x = x2 − x1 − l
dv1 dv2 dx m1v1 + m2 v2 + kx =0 dt dt dt
d ( x2 − x1 −l) dv1 dv2 mv +mv +k ( x2 − x1 −l) =0 11 11 dt dt dt dv1 dv2 k − + ( x2 −x1 −l)( v2 −v1) =0 dt dt m v1
o
其中 V 为比重计的排水体积
mg
2 ⎡ ⎤ d x ⎛d ⎞ mg − ⎢V + π ⎜ ⎟ x ⎥ ρ g = m 2 2 d t ⎝ ⎠ ⎢ ⎥ ⎣ ⎦ 2 2 d x ⎛d ⎞ mg − ρ Vg − ρ gπ ⎜ ⎟ x = m 2 dt ⎝2⎠
《机械振动》张义民—第4章第1、2节ppt
◆当振动系统需要两个独立坐标描述其运动时， 那么这个系统就是两个自由度系统。
◆两自由度系统是最简单的多自由度系统。 ◆两自由度系统的振动微分方程一般由两个联立 的微分方程组成。 ◆两自由度系统有两个固有频率及固有振型。
◆在任意初始条件下的自由振动一般由这两个固 有振型叠加，只有在特殊的初始条件下系统才按某 一个固有频率作固有振动。
大象体积庞大，走起路来 更是别具一格，四只脚移动 时分别各自相差90度的位移 差。没有一只脚做的是相同 位移的移动。
◆四只脚动物可以看作是“四个振动体耦合在一起的 系统”吗？事实上，四个振动体组成的系统的基本运动 模式，确实与所提到的那四种走路方式一模一样。
◆可是动物们为什么会按照耦合振动体的方式来行走 呢？虽说现在关于这个问题还没有定论。生物学家们认 为，掌管运动的脑神经网(由数突连接起来的神经细胞) 看起来更接近“耦合振动体”一些。有推测认为，正是 脑神经网的动力学特性，使得动物走起路来才会表现出 振动体的特点。
1998年匈牙利的物理学家塔 马斯·维塞克在布达佩斯音乐学 院举行的一场音乐会上意外地发 现了同步化的现象。
演出相当成功，落幕后观众们热烈的掌声长达 3分钟之久，而维塞克博士便在这里发现了有趣 的东西。音乐会刚一结束，观众们雷鸣暴雨般的 掌声响起，然而过了一段时间之后，观众们的热 烈的掌声显然同步化了，变成了同一种节奏的拍 手。为了答谢观众们的热情，演奏者重新走上台 来谢幕，这时的掌声又突然之间失去了刚才的节 奏，雨点般疯狂地响起。在最后长达3分钟的鼓 掌声中，狂热的掌声和同步的掌声依次交替出现。
◆强迫简谐振动发生在激励频率，而这两个坐标 的振幅将在这两个固有频率下趋向最大值。共振时 的振型就是与固有频率相应的固有振型。
振动力学4.4
n值取决于精度要求，n越多精度越高，但同时也带来计算量 越大，i (x)
i 1,2,...,n 是系统的实际模态函数，计算时以假设
模态近似，满足部分或全部边界条件，n越大， i (x) 越接近真实
的模态，解的精度越高。 用假设模态法可以建立由有限个广义坐标表示的动力学方程。
二、广义坐标的动力学方程
d L L Qi dt qi qi
i 1,2,...,n
1 1 T l 1 T 1 T q Mq qT Kq q Mq q Kq 0 f ( x, t )i ( x)dx F (t )i ( xc ) 2 2 2 qi 2
对于有集中质量和弹性支承的情况有
1 l 1 T l ( x) 2 ( x)dx m 2 ( xb ) 2 0 2
*
Vmax
2 1 l 1 1 2 ' ( xb ) k2 2 ( xa ) EI ( x) "( x) dx k1 2 2 0 2
j
k ij q j Qi
i 1,2,...,n
矩阵形式
Mq Kq Q
方程显然与集中质量法求法（解法）不同，
结果形式相同，即都离散为一个有限自由度系统。
例题：
变截面圆轴一端自由，一端固定，如图截面的极惯性
矩 I ( x) I 0 (1 x ) ，求轴扭转振动的前两阶固有频率。
1 a (1) 0.0681
0.1995 a ( 2) 1
( x) a T
x
2l 3x 2l
1 2
3x 1 ( x) sin 0.0681 sin 2l 2l
振动力学2.4
it
Fp T F0 Fp1
Fp 2 ... Fpn
T
M pi pi K pi x pi Fpi eit x
pi i2 xpi Bii2eit x
任意 i 有标准式
Fpi M pi
Bi i2
Bi
F pi M pi
2 i
F pi K pi
1 x p1 1 x x p 1.366 0.366 x p 2
F0 F0 x1 0.4553 (1 cos 1t ) 0.122 (1 cos 2t ) k k
F0 F0 x2 0.6219 (1 cos 1t ) 0.0447 (1 cos 2t ) k k
T
主质量：
1 0 1 1.366 1 0 1 4.7319 M p T M m m 1.2679 1 0.366 0 2 1.366 0.366 0
主刚度：
1 0 1 1.366 2 1 1 3k K p T K k k 1 0.366 1 2 1.366 0.366 0 2.9999k
T F0
K pi
si
主坐标运动微分方程的解为：
Bi it x pi e 2 1 si
i
变换为物理坐标
x x p (i ) x pi
i 1
n
x
i 1
n
(i ) (i )T F0 eit
K pi 1 si2
计算简谐激励响应可先计算 i， 观察
x H () F0 eit
x X sin t
2
大学物理第4章机械振动 机械波课件讲义
则系统受到的合力为
F mg FS mgi k(x l0 )i
Fx mg k(x l0 ) max
mgBiblioteka k(xl0
)
m
d2x dt 2
k
x
m
d2 dt
x
2
2 k
m
d2x dt 2
2
x
0
动力学方程
l
0A
x
F
A
x
mg
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
§4.1 简谐振动的动力学特征
振动中最简单最基本的是简谐振动 简谐振动:一个做往复运动的物体，如果其偏离平
衡位置的位移x(或角位移)随时间t按余弦(或正弦)
Fx F ,x , ax a
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
谐振动的位移、速度、加速度之间的相位关系
x A cos( t 0 ) A cos
A
sin(
t
0 )
第四章振动优秀课件
xo Acos
vo Asin
xo2
vo2
2
A2 (sin2 cos2 )
A2
A
x02
v0
2
tg vo xo
x0 0, v0 0 :在第I象限
x0 0, v0 0 :在第II象限
x0 0, v0 0 :在第III象限 x0 0, v0 0 :在第IV象限
（5）两个物体做简谐运动位相差
速度相位比位移相位超前/2。
a
dv dt
2
A cos(
t
)
am
cos(
t
)
am 2 A 称为加速度幅。
加速度与位移反相位。
4-1-3 简谐运动的旋转矢量表示法
旋转矢量A在 x 轴上的 投影点 M 的运动规律：
x Acos( t )
结论：
投影点M的运动 为简谐振动。
y
PAt源自ooMx
• 旋转矢量的模A：振幅
第四章振动
机械振动：
物体在一定的位置附近作来回往复的运动。
振动：任何一个物理量在某个确定的数值附近作 周期性的变化。
波动：振动状态在空间的传播。
任何复杂的振动都可 以看作是由若干个简 单而又基本振动的合 成。这种简单而又基 本的振动形式称为简 谐运动。
§4-1 简谐运动
4-1-1 简谐运动的基本特征
(1)弹簧振子的运动 弹簧振子： 一根轻弹簧和一个刚体构成的一个
振动系统。
根据胡可定律：
F
kx
（k为劲度系数）
（a） 在弹性限度内，弹性力F 和位移x 成正比。
（b） 弹性力F 和位移x 恒反向，始终指向平衡位置。
F
o
x
振动力学(倪振华)讲解学习
第1 章 导 论
21
第2章 单自由度系统 自由振动
第2章 单自由度系统自由振动
22
2.1 引 言
➢ 单自由度系统: 可以用一个独立坐标来确 定系统的位置及其运动规律的振动系统;
➢ 单自由度线性系统的振动是最简单的振动 系统;
确定性激励下的响应不一定是确定的 ，但随机激励下的响应一定是随机的。
第1 章 导 论
10
2. 工程振动分析的类别 振动分析：研究振动系统、激励(输入)和
响应(输出)三者之间的关系。
理论上讲，只要知道两者就可以确定第 三 者 。 这样，工程振动分析所要解决的问题 可以归纳为下面几类。
第1 章 导 论
通常系统设计要依赖于响应分析， 所以在实际工作中，响应分析和系统设 计这两个问题是交替进行的。
第1 章 导 论
13
（3）系统识别 已知振动系统的激励(输入)和响应(输出)求系
统参数,以便了解系统的特性。 系统识别包括物理参数识别(确定系统的物理
参数:质量、刚度、阻尼等)和模态参数识别(确定 或估计系统的固有特性:固有频率、振型等)。 （4）环境预测
方程解中的A称为振幅，是质量偏离静平
衡位置的最大距离; f 称为初相位。
第2章 单自由度系统自由振动
2.2 自由振动系统
34
从方程的解中还可以看出，系统属于周 期振动，振动的周期为
T 2 n
周期是系统振动一次所需要的时间，单位 为秒（s）。
周期的倒数称为频率，是系统每秒钟振动 的次数,单位为1/秒(1/s)或赫兹(Hz)。记作 f
第1 章 导 论
振动力学基础
= 1 kA2 2
结论： � 振子在振动过程中，动能和势能分别随时间变化，但
任一时刻总机械能保持不变。 � 动能和势能的变化频率是弹簧振子振动频率的两倍。
� 频率一定时，谐振动的总能量与振幅的平方成正 比。（适合于任何谐振系统）
第一个重点：如何判断物体是否作简谐振动？ 简谐振动的特征方程
（1） f = −kx
E
=
1 2
J
( dθ dt
)2
+
mgrc
(1 −
cosθ
)
o
两边对时间t求一阶导数：
J
dθ dt
d2θ ⋅ dt2
+
mgrc
dθ dt
⋅ sin θ
=
0
d2θ + mgrc θ = 0 dt2 J
θ rc c
x
解为： θ = θ0 cos(ωt + ϕ)
例. 小球在半径很大的光滑凹球面底部来回滚动， 试分析小球的运动是否简谐振动。
力学量（如位移）
机械振动
电磁量（如I 、V、 E、 B） 电磁振动 最基本、 最简单、最重要的振动是简谐振动。
§4.1 简谐振动动力学
简谐振动：凡质点的运动遵从余弦（或正弦） 规律时，其运动形式为简谐振动。
x = Acos(ω t + ϕ )
x
o
t
1. 弹簧振子模型
弹簧振子： 一根轻弹簧和一个质点构成的一 个振动系统
其中 V 为比重计的排水体积
F
o mg
mg
−
⎡ ⎢V ⎢⎣
+π
⎛ ⎜ ⎝
d 2
2
⎞ ⎟ ⎠
⎤ x⎥
⎥⎦
分析力学第四章
设
是系统的平衡位置，在
附近展开L ，则
注：a(q)展开到零阶小量
略去 m: 等效质量 k: 等效倔强系数
代入拉格朗日方程，得
二、自由振动方程的解
振动方程
解 积分常数：A—振幅； —初相位
由初始条件确定。
保守系能量守恒
解的复数形式（指数形式）
运用指数解的进行运算： 如果 那么
能量
？
条件：线性运算才可以先用指数解运算，最后再取实部。
2.实际运动存在阻尼，振幅不会随时间无限增大。
二、阻尼振动方程的解
阻尼的作用：使机械运动的能量耗散，转化为热能，使
机械运动停止(无外力时)。 1.振动系统，不再是保守系，不能引入势能函数； 2.办法：在运动方程中加进阻力项： 振子的弹力 受力 分析： 周期性策动力 介质阻尼 运动方程
求通解 齐次方程：
而这种组合一定是符合分子所属的对称性群的一个对称类的。
画出一个分子可能的结构，就能够根据这个结构求算出分子的 可以预测分子在红外光谱和拉曼光谱中的特征吸收峰。 /wiki/Normal_mode
简正坐标，通过考查分子的简正坐标可以了解分子内部运动的能量，
1 s 1 s 2 L m x k x x 2 1 2 1
一般性思路： 拉格朗日方程 m x
k x 0
1
s
( 1, 2, s)
设 x Re[C eit ] ( 1, 2, s) 代入运动方程得：
能量
三、受迫振动方程的解
系统处在随时间变化的外场中：
在平衡位置附近展开外场： 忽略
运动方程
运动方程
则运动方程降阶为 齐次方程 的通解
(1)
振动理论习题答案汇总
《振动力学》——习题第二章 单自由度系统的自由振动2-1 如图2-1 所示，重物1W 悬挂在刚度为k 的弹簧上并处于静止平衡位置，另一重物2W 从高度为h 处自由下落到1W 上且无弹跳。
试求2W 下降的最大距离和两物体碰撞后的运动规律。
解：222221v gW h W =，gh v 22=动量守恒：122122v gW W v g W +=，gh W W W v 221212+=平衡位置：11kx W =，kW x 11=1221kx W W =+，kW W x 2112+=故：kW x x x 21120=-= ()2121W W kgg W W k n +=+=ω故：tv t x txt x x n nn n nn ωωωωωωsin cos sin cos 12000+-=+-=xx 0x 1x 12平衡位置2-2 一均质等直杆，长为l ，重量为w ，用两根长h 的相同的铅垂线悬挂成水平位置，如图2-2所示。
试写出此杆绕通过重心的铅垂轴做微摆动的振动微分方程，并求出振动固有周期。
解：给杆一个微转角θ2aθ＝h α2F ＝mg由动量矩定理：ah a mg a mg Fa M ml I M I 822cos sin 12122-=-≈⋅-====αθαθ其中12cossin ≈≈θααh l ga p ha mg ml n 22222304121==⋅+θθ g h a l ga h l p T n 3π23π2π222===2-3 一半圆薄壁筒，平均半径为R , 置于粗糙平面上做微幅摆动，如图2-3所示。
试求其摆动的固有频率。
图2-3 图2-42-4 如图2-4 所示，一质量m连接在一刚性杆上，杆的质量忽略不计，试求下列情况系统作垂直振动的固有频率：（1）振动过程中杆被约束保持水平位置；（2）杆可以在铅垂平面内微幅转动；（3）比较上述两种情况中哪种的固有频率较高，并说明理由。
图T 2-9 答案图T 2-9解：（1）保持水平位置：m kk n 21+=ω（2）微幅转动：mglllF2112+=mgl1l2xx2xx'mglll2121+=k2k1ml1l2()()()()()()()()()mgk k l l k l k l mgk k l l k l l k l l l k l mg k k l l k l k l l l l k l l mg l mgk l l l k l l l l l l k l l mg l l l l x x k F x x x 2122122212121221221121212221212211211121212122211211121221112111 ++=+-++=+-⋅+++=⎥⎦⎤⎢⎣⎡+-++++=+-+='+=故：()22212121221k l k l k k l l k e++=mk en =ω 2-5 试求图2-5所示系统中均质刚性杆AB 在A 点的等效质量。
清华大学物理学概论第4章振动
x
0
振动方向旋转 2
圆
线偏振光 椭圆偏振光 圆偏振光－－二(四)分之一波片
33
例2. 振动方向垂直的两个频率比是简单 的正整数谐振动的合成 1 m ，m , n 2 n 合成轨迹为稳定的闭合曲线—李萨如图 y A x x x达到最大的次数 2 y y y达到最大的次数
36
所用的最短时间
14 受益匪浅
解：设 t 时刻到达末态
由已知画出t = 0 时刻的旋矢图
再画出末态的旋矢图
由题意选蓝实线所示的位矢 设始末态位矢夹角为 因为 t 得
o
t0
x
繁复的三角函数的运算用匀速 圆周运动的一个运动关系求得
t
7π 7π k 6 6 m
15
波
动
1
第4章 振动 §4.1 简谐振动的描述
§4.2 同方向简谐振动的合成
§4.3 垂直方向简谐振动的合成
2
§4.1 简谐振动的描述
机械振动：
物体位置在某一值附近来回往复的变化 广义振动： 一个物理量在某一定值附近往复变化 该物理量的运动形式称振动
物理量：r E H Q i 等等
1 2
则
2 A0 cos
1 2
2
t
较
2 1 2 cos t 2
随时间变化缓慢 将合成式写成谐振动形式 x A(t ) cos t
A(t ) 2 A0 cos
1 2
2
t27
合振动可看做是振幅缓变的谐振动 合成振动如图示
x1
1
t
x2
2
t
x
= 1 - 2
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
]
[u]T
[m][u]
2mL2 0
0
2mL2
[K
]
[u]T
[k
][u]
2mgL 0 0 2mgL
运动方程已解耦。
2020年8月6日 20
《振动力学》
多自由度系统振动的基本知识
教学内容
4.1 广义坐标 4.2 线性系统的运动方程及其矩阵表达式 4.3 线性变换和坐标耦合 4.4 无阻尼自由振动，特征值问题 4.5 模态向量的正交性和展开定理 4.6 系统对初始激励的响应
{x(t)} [u]{(t)}
代入运动方程，则有:
[m][u 2020年8月6日 ]{} [k ][u]{} {P(t)} 17 《振动力学》
4.3 线性变换与坐标耦合
左乘 [u]T
[u]T [m][u]{}[u]T [k][u]{} [u]T {P(t)}
记 [M ] [u]T [m][u] [K ] [u]T [k][u] {N (t)} [u]T {P(t)}
其中，[m]，[c]，[k]为n维矩阵，分别为质量矩阵，阻尼矩
阵和刚度矩阵。
其中，刚度矩阵为：
2020年8月6日 《振动力学》
k11...k1 j ...k1n
[k
]
k21...k2
j
...k2
n
.....................
kn1...knj ...knn
5
4.2 线性系统的运动方程及其矩阵表达式
k3
x2
P2
(t)
0 0 m3 x3 0
k3
k3 k4 x3 P3(t)
矩阵形式为：
2020年8月6日 [m]{x} [k]{x} {P(t)} 15 《振动力学》
多自由度系统振动的基本知识
教学内容
4.1 广义坐标 4.2 线性系统的运动方程及其矩阵表达式 4.3 线性变换和坐标耦合 4.4 无阻尼自由振动，特征值问题 4.5 模态向量的正交性和展开定理 4.6 系统对初始激励的响应
因此，只有两个坐标独立。
L2
2
m2
y
(x2 , y2 )
能完备的描述系统运动的一组独立的坐标叫广义坐标。
本例202中0年8，月6可日 选(x1, x2 ) 作为广义坐标。 3
本例《振中动力，学》也可选(θ1,θ2 ) 作为广义坐标。
多自由度系统振动的基本知识
教学内容
4.1 广义坐标 4.2 线性系统的运动方程及其矩阵表达式 4.3 线性变换和坐标耦合 4.4 无阻尼自由振动，特征值问题 4.5 模态向量的正交性和展开定理 4.6 系统对初始激励的响应
例如： k12 k2
k13 0
3.刚度矩阵为对称矩阵。
例如： k12 k21
k1
k2
k3
kn
m1
m2
m3
mn
c12020年8月6日 c2
c3
cn
7
《振动力学》
4.2 线性系统的运动方程及其矩阵表达式
4.阻尼矩阵和刚度矩阵规律相同。
例如： c22 c2 c3
c23 c3
5.取系统质心为坐标原点，则质量矩阵为对角矩阵。
(2 ) kij 2mij 0
本式称为系统的频率方程，该行列式叫特征行列式。
2020年8月6日 25
《振动力学》
4.4 无阻尼自由振动，特征值问题
展开的频率方程为:
k11 2m11 k21 2m21
k12 2m12 k22 2m22
k1n 2m1n k2n 2m2n 0
kn1 2mn1 kn2 2mn2
第四章
多自由度系统振动的基本知识
多自由度系统振动的基本知识
教学内容
4.1 广义坐标 4.2 线性系统的运动方程及其矩阵表达式 4.3 线性变换和坐标耦合 4.4 无阻尼自由振动，特征值问题 4.5 模态向量的正交性和展开定理 4.6 系统对初始激励的响应
2020年8月6日 2
《振动力学》
4.1 广义坐标
4.1 广义坐标
先看一个例子
图示双摆，质量m1， m2在平面摆动。
可以取四个直角坐标 (x1, y1) (x2 , y2 ) 来描述系统的运动。 但这四个直角坐标不独立，有：
x 1
L1
m1 (x1, y1 )
x12 + y12 = L12 (x2 - x1)2 + ( y2 - y1)2 = L22
n
n
{x(t)} Cr{x(t)(r)} Cr{u(r)}cos(rt r )
r 1
r 1
2020年8月6日 27
《振动力学》
4.4 无阻尼自由振动，特征值问题
例: 求图示系统的自然频率。
x1
x2
2k m
k
k
m
x3
2k m
解：系统的质量矩阵为
m 0 0
[m]
0
m
0
系统的刚度矩阵为 0 0 m
10
4.2 线性系统的运动方程及其矩阵表达式
P1(t)
x1 P2(t)
x2
k1
k2
k3
m1
m2
刚度矩阵为：
[k
]
k1 k2
k2
k2
k2
k3
激励向量为：
2020年8月6日
{P(t
)}
P1(t) P2 (t)
11
《振动力学》
4.2 线性系统的运动方程及其矩阵表达式
例2： 直接写出图示系统的刚度矩阵。
刚度矩阵 [k] 中的元素 kij 是使系统仅在第 j 个坐标上产 生
单位位移而相应于第 i 个坐标上所需施加的力
kij
Qi
qj qr
1 0(r
1, 2..., n, r
j)
例如
Q1
k11 Q1 k1 k2
k1 1
m1
k21 Q2 k2
k1
k2
m1
m2
c12020年8月6日 c2
《振动力学》
knn 2mnn
n 阶方阵 A 正定 A >0
是指对于任意的 n 维列向量 y，总有 yT Ay 0 成立
并且等号仅在 y 0 时才成立
若系统的质量矩阵和刚度矩阵正定，则有频率方程可解得 n个
正实根：
1 2 ... n
2最020小年8的月6日 称为基频。
《振动力学》 1
i 称为第i阶自然频率。
j 1, 2..., n
其中，u j ( j 1, 2..., n) 是一组常数。
即为同步振动: xj (t) u j =const xi (t) ui
i, j 1, 2..., n
同20步20年振8月动6日：系统在各个坐标上除了运动幅值不相同外，随时
间《变振动化力学的》规律都相同的运动
22
4.4无阻尼自由振动，特征值问题
三自由度系统
振动形式1
振动形式2
振动形式3
同步振动：系统在各个坐标上除了运动幅值不相同外，随时 间变化的规律都相同的运动
2020年8月6日 23
《振动力学》
4.4 无阻尼自由振动，特征值问题
将同步解代入运动方程:
n
n
f (t) miju j f (t) kiju j 0
j 1
j 1
i 1, 2..., n
k3
0
k3
k3 k4
2020年8月6日 14
《振动力学》
4.2 线性系统的运动方程及其矩阵表达式
k5
P2(t)
k6
k1
P1(t) k2 m2
k3
P3(t) k4
m1
m3
系统的运动方程为：
m1 0 0 x1 k1 k2
k2
0 x1 P1(t)
ቤተ መጻሕፍቲ ባይዱ
0
m2
0
x2
k2
k2 k3 k5 k6
2020年8月6日 21
《振动力学》
4.4 无阻尼自由振动，特征值问题
4.4 阻尼自由振动，特征值问题
无阻尼系统的运动方程：
[m]{x(t)}[k]{x(t)} {0}
一般形式为：
n
n
mij x j (t) kij x j (t) 0
j 1
j 1
i 1, 2..., n
设其有同步解: x j (t) u j f (t)
12
0 0
若有线性变换:
{
(t)}
1 1
11{(t)}
计算变换以后的质量矩阵和刚度矩阵。
解：
mL2 0
[m] 0
mL2
mgL ka2 ka2
[k
]
ka
2
mgL
ka2
2020年8月6日 19
《振动力学》
4.3 线性变换与坐标耦合
[u]
1 1
1 1
[u]T
1 1
1 1
则:
[M
2020年8月6日 9
《振动力学》
4.2 线性系统的运动方程及其矩阵表达式
P1(t)
x1 P2(t)
x2
k1
k2
k3
m1
m2
直接写出矩阵形式的运动微分方程：
[m]{x}[k]{x} {P(t)}
其中，质量矩阵为：
[m]
m1 0
0
m2
位移向量为：
2020年8月6日 《振动力学》
{x}
x1(t) x2 (t)
3k k 0
[k] k
2k
k
2020年8月6日
0 k 3k
28