多自由度振动系统微分方程的建立
第3章 多自由度机械振动系统 作业答案
⎤ ⎡ x1 ⎤ ⎡ p1 ( t ) ⎤ ⎢x ⎥ = ⎢ p t ⎥ − k3 ⎥ ⎥ ⎢ 2 ⎥ ⎢ 2 ( )⎥ k3 + k 4 ⎥ ⎦⎢ ⎣ x3 ⎥ ⎦ ⎢ ⎣ p3 ( t ) ⎥ ⎦ 0
d ∂T ∂T ∂U ∂D ( )− + + = Qi i ∂qi ∂qi ∂q i dt ∂q
2、拉格朗日法：
1 1 2 12 + m2 x 2 T = m1 x 2 2
U=
1 2 1 1 2 ⎤ k1 x1 + k2 (2 x2 − x1 ) 2 = ⎡ (k1 + k2 ) x12 + 4k2 x1 x2 + 4k2 x2 ⎣ ⎦ 2 2 2
Dr. Rong Guo
School of automotive studies, tongji university
⎡ k1r 2 K =⎢ 2 ⎣ − k1r
⎡3 2 ⎢ 2 Mr ⎢ ⎢ 0 ⎢ ⎣ 0
⎤ ⎥ ( k1 + k2 ) r 2 ⎦ − k1r 2
− k1r 2 ⎤ ⎡θ1 ⎤ ⎡0 ⎤ ⎥⎢ ⎥ = ⎢ ⎥ θ 2 ⎦ ⎣0 ⎦ ( k1 + k2 ) r 2 ⎦ ⎣
⎤ ⎤ ⎡ k1r 2 ⎥ ⎡θ ⎥ ⎢ 1 ⎥ + ⎢ 3 −k r 2 θ Mr 2 ⎥ ⎣ 2 ⎦ ⎣ 1 ⎥ ⎦ 2
x1 2l + k1 x1 2l + m2 x2l = 0 ⎧m1 ⎨ ⎩m2 x2l + k2 ( 2 x2 − x1 ) 2l = 0 x1 + m2 x2l + 2k1 x1 = 0 ⎧2m1 ⎨ x2 − 2k2 x1 + 4k2 x2 = 0 ⎩ m2 ⎡ 2m1 ⎢ 0 ⎣ m2 ⎤ ⎡ x1 ⎤ ⎡ 2k1 ⎢ ⎥ + ⎢ −2 k m2 ⎥ x 2 ⎦⎣ 2⎦ ⎣ 0 ⎤ ⎡ x1 ⎤ ⎡0 ⎤ ⎢ x ⎥ = ⎢0 ⎥ 4k 2 ⎥ ⎦⎣ 2⎦ ⎣ ⎦
第4章 振动系统的运动微分方程
(d)
分析杆 AB ，列写 AB 的运动微分方程，如图（c）
m2 &x&C = − X A
(e)
m2 &y&C = −YA − m2 g
(f)
1 12
m2l 2ϕ&&
=
X
A
l 2
cosϕ
+ YA
l 2
sin ϕ
(g)
运动学方程
xC
=
xA
+
l 2
sin
ϕ
,
x&C
=
x& A
+
l ϕ& cosϕ 2
yC
=
−
l cosϕ , 2
y& C
=
l ϕ& sinϕ 2
&x&C
=
&x&A
−
l ϕ& 2 2
sin ϕ
+
l ϕ&& cosϕ 2
(h)
&y&C
=
l ϕ&& sin ϕ 2
+
l ϕ& 2 2
cos ϕ
(i)
上述 9 个方程包含 &x&A ,ε , &x&C , &y&C ,ϕ&&, X A ,YA , F, N 等 9 个未知量，由上述 9 个方程消去
解：系统具有两个自由度，选图示 AB 与铅垂线的夹角ϕ 及圆轮中心 A 的位移 xA 为广
义坐标。
分析圆轮 A ，受力图如图（b）所示。列写圆轮 A 的运动微分方程：
汽车振动基础第4章-多自由度(定稿)
k11 k1 x1 k2 x1 k1 k2
k21 k12 k2 x1 k2
k22 k2 x2 k3 x2 k2 k3
j2
k31 k13 0
k32 k23 k3 x2 k3
0 k1 k 2 k 2 K k 2 k 2 k3 k3 0 k3 k3
– 拉格朗日法
• 方程的形式
广义坐标
qi (i 1, 2,3,, n)
T：系统的总动能
d T T ( ) Qi 0 dt qi qi
i 1, 2,3, , n
对应于第i个广义 坐标的广义力
– 保守系统
» 系统作用的主动力仅为势力 Qi
d T T U ( ) 0 dt qi qi qi
m2 m22 m3 4
④柔度矩阵的影响系数法
F ij
柔度影响系数 ij 的意义是在第j个坐标上施加单位力作用时，在第i个坐 标上引起的位移。 例题4-8 用影响系数法求图示系统的柔度矩阵
11 F 21 31
12 22 32
13 23 33
也可写成 其中
或
或
MX KX 0
力方程 位移方程
K 1MX X 0
m x 0 或 x
称为柔度，而
FMX X 0
1 称为柔度矩阵
1 k
FK
②刚度矩阵的影响系数法
K kij
刚度影响系数 k 的意义是使系统的第j个坐标产生单位位移，而其它的 ij 坐标位移为零时，在第i个坐标上所施加的作用力的大小。
仅代表外部激励 广义力
机械振动学(第三章)-多自由度振动系统
装备制造学院
College of Equipment Manufacture
利用直接法，对下图所示的三自由度振动系统建立微分方程。。
装备制造学院
College of Equipment Manufacture
解：1）受力分析 选取 m1， m2和m3离开平衡位置的坐标x1， x2和 x3 为3 个独立 坐标。受力分析如图所示 2）建立振动微分方程 (c c ) x c x ( k k ) x k x p (t ) x m1： m 2 2 2 2 2 ( c 2 c 3 ) x 2 c2 x 1 c 3 x 3 ( k 2 k 3 ) x 2 k 2 x1 k 3 x 3 p 2 ( t ) x m2: m 2 2 2 2 3 c 3 x 3 c3 x 2 k 3 x3 k 3 x 2 p 3 (t ) x m3: m 3
装备制造学院
College of Equipment Manufacture
本章结束
装备制造学院 College of Equipment Manufacture
3 ）如果将应为能量耗散函数 D 引起的阻尼力也从其他的非势 力的广义力中分离出来，并使Qi仅代表外部作用的广义激振力， 则可将非保守系统的拉格朗日方程改为：
d dt ( T i q ) T i q U qi D i q Q i ( i 1, 2 , 3 ,...., n )
车 身 车 轮 二 自 由 度 振 动 问 题
装备制造学院
College of Equipment Manufacture
装备制造学院
College of Equipment Manufacture
结构动力学之多自由度体系的振动问题
2.760 3.342 1
0.163
0.924
2.76
柔度法
利用刚度法的方程间接导出柔度法方程：
由刚度法振幅方程：
令λ=1/ω2 得频率方程：
( [K]－ω2 [M] ){Y}={0}
前乘[K]－1=[δ]后得： ( [I ]－ω2 [δ] [M] ){Y}={0} ( [δ] [M] － λ [I ] ){Y}={0} ┃ [δ] [M] － λ [I ] ┃=0
刚度法
2）如果初始条件是任意的，则任其自然 后， 系统所发生的振动就不是按主振型的简谐自由 振动，而是复杂的周期振动，这时可以用各阶 主振动的线性组合来描述它，也就是说其通解 表为各个特解之和，即
y j sin( j t v j )
j 1 n
所以系统的任意振动可以表示为各个主振动 的叠加。
Yij为正时表示质1 1.293 5Y11 6.70Y21 3 0 量mi的运动方向与单 3Y 1.707 0
21
Y
(1)
0.163 0.569 1
0.569
5Y13 5.027Y23 3 0 (1) Y 3Y21 10.027 0 3.342 1.227
1 1 4 0 , m m 2 9
展开得: 解之:
3 15 2 42 30 0
ξ1=11.601，ξ2=2.246，ξ3=1.151
1 m
三个频率为：
1 0.2936
1 1 3 0.9319 m m 3）求主振型： （令Y3i=1）将λ1代入振型方程： （[δ] [M ]－λ1[I]）{Y}＝0的前两式：
第二章(多自由度系统的运动微分方程)详解
k11 k 21 kN1
k1 j k2 j k Nj
k1N k2 N k NN
刚度影响系数 kij ：第 j 个自由度产生单位位移，其他自由度位移为零时， 需要在第i 自由度处沿着位移方向施加的力。
用影响系数法建立系统的运动微分方程
【例】用影响系数法写出图示系统的刚度矩阵。
多自由度振动系统
Piezoelectric actuator
基于压电作动器的垂尾抖振主动抑制 （此系统有一、两千个自由度（3D实体单元） ）
Z Y
X
第二章： 多自由度系统的运动 微分方程
第二章:多自由度系统的运动微分方程
第一讲：
1.建立多自由度系统运动微分方程的 各种方法的概述 2.用牛顿第二定律列写系统的运动微 分方程 3.用影响系数法建立系统的运动微分 方程
F1 1
k3
m2
k2 (d11 d21 )
m1
k2 (d11 d21 ) k1d11 1
d 21 k2 (d11 d21 )
F2 0
d11
k3d21
k2 k3 k1k2 k1k3 k2 k3 k2 k1k2 k1k3 k2 k3
m2
d 21
k2 (d11 d21 ) k3d21 0
上次课内容回顾
3.刚度影响系数
刚度影响系数 kij ：第 j 个自由度产生单位位移，其他自由度位移为零时， 需要在第 i 自由度处沿着位移方向施加的力。
4.柔度影响系数
柔度影响系数 dij ：第 j 个自由度上作用单位力，其他自由度作用力为零时，
在第 自由度上产生的位移。 i
5.刚度矩阵和柔度矩阵的关系
第2章——多自由度系统的振动——运动方程建立方法0425
船体振动基础1第章多自由度系统的振第2章多自由度系统的振动一、引言二、两自由度系统的振动三、多自由度系统的振动四、振动方程建立的其他方法2有阻尼的多自由度系统振动1、拉格朗日方程式1、拉格朗日方程式P38拉格朗日法是建立微分方程一种简单的方法：先求出系统的动能、势能，进而得出质量矩阵和刚度矩阵.优点：系统的动能和势能都是标量，无需考虑力的方向。
141、拉格朗日方程式P38拉格朗日第二类方程式适用于完整约束的系统。
完整约束完整约束：当约束方程本身或约束方程通过积分后可以下式所示的形式表示时，称为完整约束。
不完整约束：当约束方程本含有不能积分的速度项时，系统的约束称为不完整约束。
具有不完整约束的系统，系统的自由度不等于广义坐标数自由度数小于广义坐标数于广义坐标数，自由度数小于广义坐标数。
151、拉格朗日方程式P3811•位移方程和柔度矩阵P40对于静定结构，有时通过柔度矩阵建立位移方程比通过对于静定结构有时通过m1x1x2以准静态方式作用在梁上。
梁只产生位移（即挠度），不产生加速度。
的静平衡位置为坐标P1=1 f11 f21 f12P2=1 f22（1）P1 = 1、P2 = 0 时 m1 位移：x1 = f11 m2 位移：x2 = f 21 （3）P1、P2 同时作用 m1 位移： 位移 x1 = f11 P 1 + f12 P 2 m2 位移：x2 = f 21 P 1 + f 22 P 2（2）P1 = 0、P2 = 1 时 m1 位移：x1 = f12 m2 位移：x2 = f 22P1 m1 x1 x2 P2 m2P1=1 f11 f21 f12 P1 m1 x1P2=1 f22 P2 m2 x2P 同时作用时 1、P 2 同时作用时：x1 = f11P 1 + f12 P 2 x2 = f 21P 1 + f 22 P 2矩阵形式 X = FP 矩阵形式：⎡ x1 ⎤ X =⎢ ⎥ ⎣ x2 ⎦f ij 柔度影响系数f12 ⎤ f 22 ⎥ ⎦⎡ f11 F=⎢ ⎣ f 21⎡P 1⎤ P=⎢ ⎥ ⎣ P2 ⎦物理意义： 系统仅在第 j 个坐标受到 单位力作用时相应于第 i 个坐标上产生的位移柔度矩阵P1 m1 x1P2 m2 x2P1(t) m1 m2P2(t)&1 m1 & x&2 m2 & xX = FP⎡ x1 ⎤ ⎡ f11 ⎢x ⎥ = ⎢ f ⎣ 2 ⎦ ⎣ 21 f12 ⎤ ⎡ P 1⎤ ⎢P ⎥ f 22 ⎥ ⎦⎣ 2 ⎦当P 1、P 2 是动载荷时 集中质量上有惯性力存在⎡ x1 ⎤ ⎡ f11 ⎢x ⎥ = ⎢ f ⎣ 2 ⎦ ⎣ 21 f12 ⎤ ⎡ P && 1 (t ) − m1 x1 ⎤ ⎢ P (t ) − m & ⎥ f 22 ⎥ & x 2 2⎦ ⎦⎣ 2⎡ x1 ⎤ ⎡ f 11 ⎢x ⎥ = ⎢ f ⎣ 2 ⎦ ⎣ 21位移方程:f 12 ⎤⎛ ⎡ P1 (t ) ⎤ ⎡m1 ⎜⎢ −⎢ ⎥ ⎥ ⎜ f 22 ⎦⎝ ⎣ P2 (t ) ⎦ ⎣ 0&1 ⎤ ⎞ 0 ⎤⎡ & x ⎟ ⎥ ⎢ ⎥ &2 ⎦ ⎟ m2 ⎦ ⎣ & x ⎠&& ) X = F ( P − MXP1(t) m1 m2P2(t)⎡ x1 ⎤ X =⎢ ⎥ ⎣ x2 ⎦⎡P 1 (t ) ⎤ P=⎢ ⎥ P ( t ) ⎣ 2 ⎦&1 m1 & x&2 m2 & x位移方程 位移方程：&& ) X = F ( P − MX也可按作用力方程建立方程：&& + KX = P MX刚度矩阵&& + X = FP FMX柔度矩阵与刚度矩阵的关系 柔度矩阵与刚度矩阵的关系：&& KX = P − MX若K非奇异F=K−1FK = I&& ) X = K −1 ( P − MX应当注意：对于允许刚体运动产生的系统（即具有刚体自由度的系统） ， 柔度矩阵不存在。
3振动系统的运动微分方程
机械与结构振动
Mechanical and Structural Vibration Mechanical and Structural Vibration
制作与设计 贾启芬
返回总目录
第3章 振动系统的运动微分方程
3.1 牛顿定律和普遍定理 3.2 拉格朗日运动方程 3.3 刚度影响系数 作用力方程 3.4 柔度影响系数 位移方程
1 δ 13 k1 1 δ 23 = k δ 33 11 k1
1 k1 1 1 + k1 k 2 1 1 + k1 k 2
1 1 + k1 k 2 1 1 1 + + k1 k 2 k 3 1 k1
Mechanical and Structural Vibration
K =K
Mechanical and Structural Vibration
T
第3章 振动系统的运动微分方程
3.4 柔度影响系数 位移方程
Mechanical and Structural Vibration
3.4 柔度影响系数 位移方程
在单自由度的弹簧—质量系统中，若弹簧常数是 ， 在单自由度的弹簧 质量系统中，若弹簧常数是k，则 质量系统中
画出各物块的受力图根据平衡条件， 画出各物块的受力图根据平衡条件，有
k11 = k1 + k 2，k 21 = −k 2，k 31 = 0
Mechanical and Structural Vibration
3.3 刚度影响系数 作用力方程
同理， 同理，令 x1 = 0，x 2 = 1，x 3 = 0
&& && && m11 x1 + m12 x 2 + ⋅L+ m1n x n + k 11 x1 + k 12 x 2 +L+ k 1n x n = 0 m x + m x +L+ m x + k x + k x +L+ k x = 0 21 &&1 22 &&2 2 n &&n 21 1 22 2 2n n LL mn1 x1 + mn 2 x 2 +L+ mnn x n + k n1 x1 + k n 2 x 2 +L+ k nn x n = 0 && && &&
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
D x1
D
.
.
c1 x1 c2 ( x2 x1 ) (c1 c2 ) x1 c2 x2
c2 ( x2 x1 ) c3 ( x3 x2 ) c2 x1 (c2 c3 ) x2 c3 x3
. . . . . . .
.
.
.
.
.
x2
D x3
解毕
拉格朗日方程法
系统存在粘性阻尼时：
• 其中：qi为广义坐标；T为体系的总动能；V为体系的势能；Qi为 广义坐标qi对应的非保守主动力；D为耗能函数
弹性势能 耗能函数
体系动能
分别计算偏导
பைடு நூலகம்
T 0 xi
T xi
.
mi xi
.
.. d T ( . ) m xi dt x i
写成矩阵形式得
m1 0 0 m 2 0 0 .. . 0 x1 (c1 c2 ) c2 0 x1 (k1 k2 ) k2 0 x1 F1 (t ) .. . x k x F (t ) 0 x c ( c c ) c ( k k ) k 2 2 3 3 2 2 2 3 3 2 2 ..2 . m3 c3 c3 k3 k3 x3 0 x3 0 x3 F3 (t )
V k1 x1 k2 ( x2 x1 ) (k1 k2 ) x1 k2 x2 x1
V k2 ( x2 x1 ) k3 ( x3 x2 ) k2 x1 (k2 k3 ) x2 k3 x3 x2
V k3 x2 k3 x3 x3
多自由度系统振动微分方 程的建立
例：如下图所示三个自由度系统的有阻尼强迫振动，分别用牛顿第二定律、
拉格朗日方程建立其振动微分方程。
牛顿第二定律法
由牛顿第二定律建立上述系统振动微分方程，分析各质量单元受力图，此时 不加惯性力。 牛顿第二定律表达式: F=ma 其中F为物体所受所有外力之和
• 分别以m1，m2，m3为研究对象列方程
.
c （ c3 x2 c3 x3 3 x3 x2 ）
.
.
.
.
以上式中i=1,2,3
带入拉格朗日方程，得出
写成矩阵形式得：
• 整理得：
. . .. m1 x1 (c1 c2 ) x1 c2 x2 (k1 k2 ) x1 k 2 x2 F1 (t ) . . . .. m2 x2 c2 x1 (c2 c3 ) x2 c3 x3 k2 x1 (k2 k3 ) x2 k3 x3 F2 (t ) .. . . m3 x3 c3 x2 c3 x3 k3 x2 k3 x3 F3 (t )