第三章多自由度体系的振动2
第六讲--多自由度系统振动-2
解: 1)求柔度系数
m
31
k/5
m
21
k/3
P=1
2m k
11
32 4
P=1
22 4 12
P=1
33 9
23 4 13
11 1/ k 21 31 11
22
1 k
1 k /3
4
22
1 k
1 k/3
1 9
k /5
3.3.1 柔度法
1 1 1
柔度矩阵: [ ] 1 4 4
1 4 9
2)求频率
2 0 0
质量矩阵: [M] m 0 1 0
0 0 1
由频率方程: M I 0
2 1 1 m 2 4 4 0 ,
2 4 9
展开式为: 3 15 2 42 30 0
1 m m2
方程三个根为: 1 11.601 2 2.246 3 1.151
三个频率为:
1 0.2936
k m
4Y
4 4
3.4.1 主振型矩阵与正则坐标
(2)正则坐标 任意一个质点的位移 y 都可按主振型来组合:
y1 1Y11 2Y12 3Y13 y2 1Y21 2Y22 3Y23
yi 1Yi1 2Yi2 3Yi3
yn 1Yn1 2Yn2 3Yn3
nY1n nY2n
y1
y2
Y1 Y121
Y YYY132111
Y2 1
Y2 2
Y32
Y3 1
Y3 2
Y33
Y14 Y4
2
Y34
Y41
Y2 4
Y3 4
Y44
主 振
型 矩 阵
第一振型
1
第三章 多自由度系统
例题
例题
再将初始条件(2)代入式,得
A1(1) 0,
1
2
,
A1(2) 1,
2
π 2
x1 (t) cos p2t cos3
kt m
(cm), x2 (t) cos p2t cos3
kt m
(cm)
这表明,由于初始位移之比等于该系统的第二振幅比,因 此,系统按第二主振型以频率p2作谐振动。
将第一固有频率p1代入 x1 A1 sin( pt )
x2 A2 sin( pt )
normal mode 第一主振动
x11 x21
A11 A21
sin( sin(
p1t p1t
1 1
) )
第二主振动
x12 x22
解:系统的质量矩阵和刚度矩阵为
M
m
0
m0 ,
K
k1 k2
k2
k2 k2 k3
5k 4k
4k
5k
将M、K代入频率方程,得
p1
k, m
p2 3
k m
对应的两个主振型和振幅比为
1
A2(1) A1(1)
1,
2
A2(2) A1(2)
代入上式得到
1
2
2
(1)
(2)
1 2
0
因此得到双摆作自由振动的规律
位置的距离x1、x2为广义坐标,两物体沿x方向的受力如图示, 它们的运动微分方程分别为
机械振动学(第三章)-多自由度振动系统
装备制造学院
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
机械振动二自由度
m mL1
mL1 mL12
I
C
y"
q "
k1 k2
k2L
k2L k2 L2
yA
q
{0}
这个方程存在弹性耦合和惯性耦合
3. 取广义坐标为yA,yB
系统 的动
ET
1 2
y
A
',
yB '
mL22 I L2
mL1L2
c
I
c
mL1L2 L2
mL12 I
I
c
c
yA' yB '
k2
k2 k2 k3
x1 x2
1 {x}T 2
[ K ]{x}
系统的能量 耗散函数
D
1 2
c1x'12
1 2
c2
(x'2
x'1
)2
1 2
c3 x'22
1 2
{x'1
,
x'2
}c1c2c2
c2 c2 c3
xx''12
1 2
{x'}T
[C]{x'}
利用这三个函数可以分别求出三 个矩阵的各个元素
工程振动问题中遇到的质量矩阵一般都是正定矩阵。对于静定和超静定结构,刚度矩
阵也是正定矩阵。
上面关于质量矩阵、阻尼矩阵和刚度矩阵情况的讨论完全可以推广到任意的二自
由度系统和n自由度系统。
[M
]
m1
0
0
m2
[C
]
c1 c2
c2
c2
c2
c3
[K
]
k1 k2
《多自由度系统振动》课件
课程目的
理解多自由度系统振动的 特性,包括固有频率、模 态振型等。
掌握多自由度系统振动的 基本原理和数学模型。
学习多自由度系统振动的 分析方法,包括直接法、 模态法和传递矩阵法等。
控制算法则是实现控制策略的具体计算方法。常见的控制算法包 括PID控制、状态反馈控制、最优反馈控制等。这些算法可以根 据系统的特性和要求进行选择和优化。
05
多自由度系统振动应用
机械系统振动控制
机械系统中的多自由度振动问题广泛存在,如旋转机械、往复机械和柔性机械等 。控制这些振动可以提高机械系统的稳定性和可靠性,减少磨损和疲劳,延长使 用寿命。
多自由度系统振动
CONTENTS
• 引言 • 多自由度系统振动基础 • 多自由度系统振动特性 • 多自由度系统振动控制 • 多自由度系统振动应用 • 课程总结与展望
01
引言
课程背景
机械系统振动是工程领域中常见的问题,多自由度系统振动 更是其中的重要分支。随着科技的发展,多自由度系统在许 多领域如航空航天、交通运输、能源等都得到了广泛应用, 因此对多自由度系统振动的研究具有重要意义。
多自由度系统振动与多个学科领域密切相关,如结构力学、流体力学 和声学等,需要加强这些交叉学科领域的应用研究。
多自由度系统振动实验平台的搭建与验证
为了验证多自由度系统振动理论和方法的有效性,需要搭建更加先进 的实验平台,并开展更加系统的实验研究。
谢谢您的聆听
THANKS
被动控制技术
被动控制技术是通过改变系 统的刚度、阻尼和/或质量分 布来减小系统的振动。被动 控制技术不需要外部能源, 而是利用自然现象或物理效 应来减小系统的振动。
结构动力学-多自由度系统的振动
sin(1t sin(1t
1) 1)
A2Y1(2) A2Y2(2)
sin(2t sin17
m1 y1 m2 y2
(k1 k2
y1
k2 ) y1 (k2
k2 y2 k3 ) y2
0 0
方程的全解:
y1(t) y2 (t)
A1Y1(1) A1Y2(1)
2
1 2
k11 m1
k22 m2
1
2
k11 m1
k22 m2
2
k11k22 k12k21 m1m2
正实根,仅依赖于结构体系的物理性质,
即质量和弹簧刚度。
2021/6/24
13
2
1 2
k11 m1
k22 m2
1
2
k11 m1
k22 m2
2
k11k22 k12k21 m1m2
具有两个自由度的体系共有两个自振频率, 1 表示其中最小的圆频率,称为第一圆频率或 基本圆频率(fundamental frequency); 2 称为第二圆频率。
y1 (t) y2 (t)
YY12ssiinn((tt))
1)、在振动过程中,两个质点具有相同的频
率 和相同的相位角 。
2)、在振动过程中,两个质点的位移在数值上 随时间而变化,但两者的比值始终保持不变:
y1(t) Y1 常数 y2 (t) Y2
2021/6/24
10
主振型:结构位移形状保持不变的振动形式称
设方程的解为:
y1(t) Y1 sin(t ) y2 (t) Y2 sin(t )
2k m 2
k
2k
k
m
2
YY12
0
第三章(多自由度系统的振动)
固有振型的正交性
加权正交性的简洁表示
T r M s 0, r s
M s M r , r s
T r
rT M s M r rs
rs
def
1, r s 0, r s
rT K s 0, r s
rT K s K r , r s
x
x1 1
节点
x3 1
3 2
k m
x2 1
理解固有振型
理解固有振型
理解固有振型
返回
固有振型的正交性
1.固有振型的归一化
2 r 1 3 2 r 1 3
都是固有振型向量 ① 按某一自由度的幅值归一化
k m
理解固有振型
3k k 0 m 0 0 1 0 k 2k k 2 0 m 0 2 0 0 k 3k 0 0 m 3 0
u(t ) sin(t )
对任意时间都成立
( M ) 0, 2
特征方程 特征值
det( K M ) 0
r (r 1, 2, N )
有非零
r (r 1, 2, N )
特征向量
u(t ) ψa sin(t ) φ sin(t )
结论: 系统存在形如 形式的同步振动。
u(t ) φ sin(t )
多自由度系统的固有振动
2.多自由度系统的固有振动
Mu(t ) Ku(t ) 0
( K 2 M ) sin(t ) 0
第r阶模态质量
固有振型关于刚度矩阵加权正交性 T 当 rs 时 r K s 0 T r K s K r 当 rs时
第三章多自由度体系的振动2
10
§ 3.3 主振型的正交性和正则坐标
如果两个主振型 {Y (1) }和 {Y ( 2) }彼此不正交,即
{Y (1) }T [M ]{Y ( 2) } 0
取一个由 {Y (1) } 和 {Y ( 2) } 组成的新的主振型,即
{Y
(1, 2)
} {Y } c{Y }
主振型正交的物理意义:
1)每一主振型相应的惯性力在其他主振型上不做功,其数学表 达式为 :
s2 {Y (t ) }T [M ]{Y ( s ) } 0
2)当一体系只按某一主振型振动时,不会激起其他主振型的 振动。 3)各个主振型都能够单独出现,彼此间线性无关。
9
§ 3.3 主振型的正交性和正则坐标
(t ) T (t ) T (t ) T 2 (s)
{Y } [ M ]s {Y }sin(s t )
2 (s)
s {Y } [ M ]{Y }sin(s t ) 0
2 (t ) T (s)
s2 {Y (t ) }T [M ]{Y ( s ) } 0
8
§ 3.3 主振型的正交性和正则坐标
(1) T ( 2)
T
同理:
{Y (1) }T [ K ]{Y (3) } 0.001 k 0 {Y ( 2) }T [ K ]{Y (3) } 0.00001 k 0
6
§ 3.3 主振型的正交性和正则坐标
对任意一个位移向量{y} ,将其写成主振型的线性组合:
{ y} 1{Y (1) } 2 {Y ( 2) } n {Y ( n ) } i {Y (i ) }
Y
1
第三章二自由度系统
二自由度系统振动 / 不同坐标系的运动微分方程
以汽车的二自由度振动模型为例
汽车板簧以上部分被简化成为一根刚性杆,具有质量m和绕质心 的转动惯量Ic。质心位于C 点。分别在A点和B点与杆相联的弹性 元件k1、k2为汽车的前,后板簧。
若系统有 n 个自由度,则各项皆为 n 维矩阵或列向量
二自由度系统振动 / 运动微分方程
式中:
[M
]
m11 m21
m12
m22
m1
0
0
m2
[K
]
k11 k 21
[C]
c11 c21
k12
k
22
k1 k2
k2
c12
c22
2 ET x1x1
2 ET x12
m1
m12
2 ET x1x2
2 ET x2x1
m21
0
m22
2ET x2x2
2 ET x22
m2
[M
]
m11 m21
m12
m22
m1
0
0
m2
二自由度系统振动 / 能量法
(t ) (t)
如同在单自由度系统中所定义的,在多自由度系统中 也将质量、刚度、位移、加速度及力都理解为广义的。
第三章-多自由度系统振动6.19
第三章 多自由度系统振动多自由度系统和单自由度系统的振动特性是有区别的。
单自由度系统受初始扰动后,按系统的固有频率作简谐振动。
多自由度系统有多个固有频率,当系统按某一个固有频率作自由振动时,各独立坐标在振动过程中相互关系是固定的,这个关系叫振幅比,也叫作主振型或模态。
主振型是多自由度系统以及弹性体振动的重要特征。
多自由度系统的振动方程是多个二阶微分方程组,这些方程一般是耦合的。
多自由度振动的求解有两种方法:直接积分法和振型叠加法。
直接积分法可直接根据微分方程求出响应,涉及的概念不多且有应用软件,本章不做介绍。
振形叠加法要先求出系统的固有频率和振型,在此基础用叠加法求响应,物理概念清楚、并且是模态分析与参数识别的理论基础。
因此本章将先用较多的篇幅介绍多自由度系统的固有振动特性、振型叠加法和传递函数。
3.1 振动微分方程虽然一些多自由度系统数目较多,有些相当复杂,但建立多自由度系统振动微分方程并没有新理论和方法,都是动力学基本理论和方法,本节只通过例题介绍多自由度系统振动微分方程基本形式。
[例一] 试建立图3-1所示3自由度系统的运动微分方程。
三个质量只作水平方向的运动,并分别受到激振力()t P 1,()t P 2和()t P 3的作用,质量块的质量分别为1m ,2m 和3m ,弹簧刚度分别为1k ,2k 3k 和4k ,阻尼分别为1c ,2c 3c 和4c 。
图3-1 3自由度系统解:分别用三个独立坐标1x ,2x 和3x 描述三个质量块的运动,坐标原点分别取在1m ,2m 和3m 的静平衡位置。
质量块的速度分别为1x,2x 和3x ,加速度分别为1x,2x 和3x 。
每个质量块的受力图如3-2(a 、b 、c )所示,则由受力图根据牛顿第二定律,得系统的运动方程为:图3-2 (a) 图3-2(b)图3-2(c))()()(1212112121111t P x x c x c x x k x k xm +------= )()()()()(232321232321222t P x x c x x c x x k x x k x m +---+---= )()()(3343233432333t P x c x x c x k x x k xm +--+--= 或)()()(1221212212111t P x k x k k x c x c c xm =-++-++ )()()(23323212332321222t P x k x k k x k x c x c c x c x m =-++--++- )()()(3343233432333t P x k k x k x c c x c xm =++-++- 上述方程组可以用矩阵表示为:⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+--+--++⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+--+--++⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡)()()(000032132143333222213214333322221321321t P t P t P x x x k k k k k k k k k k x x x c c c c c c c c c c x x x m m m三个二阶微分方程是耦合的,这是因为矩阵中有非零的非对角元素。
第三章二自由度系统振动的理论及工程应用
� �
x1 x2
�0� � 0��
�a�
令
b
�
�12m2 �11m1
�
c
�
�21m1 �22m2
�
d
�
1�
�11m1
e
�
1
�22m2
�b�
则方程�a�可改写为
�x�1 � b�x�2 � dx1 � 0� c�x�1 � �x�2 � ex2 � 0
�c�
设上述方程解的形式为
x1 � Asin(� t � � )� x2 � B sin(� t � � )
整理后得
(b � � 2 )A � cB � 0� � dA � (d � � 2 )B � 0
系统发生振动时�方程具有非零解
b��2
�d d
�c
��2
�0
� 4 � (b � d )� 2 � d(b � c) � 0
�2 1,2
�
b
�d 2
�
(b � d )2 � d(b � c) 2
整理得
A2 B2
�
c
b � �22
�
d � �22
d
�
1
�2
A1 B1
�
c
b � �12
�
d � �12
d
�
1
�1
�1
�
B1 A1
�
b � �12
c
�
1[b � d c2
�
( b � d )2 � cd ] � 0 2
�2
�
B2 A2
�
b � �22
c
�
1[b � d c2
多自由度体系自由振动讲解
代入振动方程可得
K 2 M A 0 -----------振型方程
K11 m1 2
K 21 K n1
K12
K22 m2 2
K1n
K2n
0
Kn2
Knn mn 2
-----频率方程
3. 柔度矩阵与刚度矩阵的关系
FEK2 K21 y1(t) K22 y2 (t)
振动方程------受力平衡方程
m1
mm21yy21
(t (t
) )
FEK1 FEK 2
0 0
m1y1(t)
mm21yy21
(t) (t)
K11 y1 (t) K21 y1 (t)
方程(1)的解设为 : y(t) Asin(t ) 式中, A A1 A2 An T
把 y(t) Asin(t ) 代入(1)
A 2 M A
M EA 0
记
1
2
----------------(2)
3)振型图
画振型图时, 完全按照2个振型中的量值,与假定的2 个位移方向相协同。
A1 1.0 1.0T A2 1.0 1.0T
1
1
1
1
第一振型
第二振型
多自由度振动
重 点:频率、振型 难 点:建立方程、求刚度
系数、柔度系数
多自由度体系的自由振动
主要内容:振动方程、振型方程、频率方程及振型图
一、柔度法建立振动方程 1. 两个质点的振动
m2
y2 (t)
y1(t) y2(t) 由质点1与质点2的惯性力共同产生 m1
第三章 多自由度系统振动
U = U ( q1 , q2 ,..., qn )
通常将静平衡位置作为势能零点, 并且以静平衡 通常将静平衡位置作为势能零点, 位置为坐标原点。 位置为坐标原点。 我们研究的是在静平衡位置附 近的微振动, 近的微振动,则将 U 在静平衡位置作泰勒展开有
∂U U = U0 + ∑ i =1 ∂qi
0
q
对应的广义力,阻尼力,耗散力。 对应的广义力,阻尼力,耗散力。系统的第 k 个 质点受到的阻尼力
& Rk = − β k ⋅ rk
与势能形式上对应存在一个耗散函数
m n 1 ∂rk dqi n ∂rk dq j 1 & & Φ = ∑ β k ⋅ rk ⋅ rk = ∑ β k ⋅ ∑ ⋅ ⋅∑ ⋅ dt j =1 ∂q j dt k =1 2 k =1 2 i =1 ∂qi
kn 2 − mn 2ωi2 ) ⋅ ϕ 2i + ... + ( knn − mnnωi2 ) ⋅ ϕ ni = ( mn1ωi2 − kn1 ) ϕ1i (
n − 1 个方程,n − 1 未知数, 个方程, 未知数, 最终可求出 ϕ2i ,..., ϕni 用 ϕ1i
表示,其余都与其成一定比例。 表示,其余都与其成一定比例。 与其成一定比例
系统的能量等于各阶主振动的能量之和不同阶之间能量不发生变换每一阶主振动的动能和势能在内部交换总和保持常数34多自由度系统的受迫振动mxcxkx1特征值分析求出无阻尼的各阶固有频率和各阶主振型2模态叠加方法分解解耦期望阻尼阵也和mk一样具有正交性即如果这样就可以使用模态叠加法进行解耦分析求解
结 构 动 力 学
1 n n ∂ 2U U = ∑∑ 0 qi q j 2 i =1 j =1 ∂qi ∂q j , 令
《多自由度系统振动》PPT课件教案资料
2022/2/12 《振动力学》
代入,得: (FM I)φ 0 特征方程: FMI 0 18
多自由度系统振动 / 多自由度系统的自由振动/模态
• 多自由度系统的模态(主振型)
正定系统: M X KX 0
主振动: X φ asi nt ()
XRn M、 KRnn
0 φRn
特征值问题: (K2M)φ0
7
《振动力学》
多自由度系统振动 / 多自由度系统的自由振动
• 多自由度系统的固有频率
作用力方程: MXKXP(t) XRn
自由振动方程: MXKX0
和单自由度系统一样,自 由振动时系统将以固有频 率为振动频率。
在考虑系统的固有振动时,最感兴趣的是系统的同 步振动,即系统在各个坐标上除了运动幅值不相同外, 随时间变化的规律都相同的运动 。
则有:(TTAT)T=TTAT(TT)T=TTAT 正定性质:若原来的刚度矩阵K 正定,则(TTKT)仍正定。
因此坐标变换X =TY 不改变系统的正定性质。 对于质量矩阵也如此。
2022/2/12
5
《振动力学》
回顾:单自由度系统自由振动-无阻尼自由振动
小结:
单自由度系统自由振动分析的一般过程:
1、由工程装置建立自由振动的一般方程,并写出振动的标准方程; 2、根据标准方程,建立本征方程并计算得到本征值; 3、根据本征值,写出标准方程的通解; 4、根据初始条件,计算标准方程的特解。
f(t)asint(), 0
f(t)a tb,
0
主振动
(1)正定系统 0
只可能出现形如 X φ asi nt ()的同步运动。
系统在各个坐标上都是按相同频率及初相位作简谐振动。
(2)半正定系统 0
结构动力学之多自由度体系的振动问题ppt课件
448 (1 536)2
m1m2l 6 (EI )2
0
解得
21
23l3 (m1 m2 2 1 536EI
)
529(m1 m2 )2l6 41 5362 (EI )2
448m1m2l 6 1 5362 (EI )2
从而得第一和第二阶自振频率
1
1
1
2
1
2
为了确定第一阶振型,可将1代入平衡方程。
其展开式是关于λ的n次代数方程,先求出λi再求 出频率ωi
柔度法
(11m1 )
12m2
...
21m1 ( 22m2 ) ...
...
...
...
1n mn 2nmn 0
...
n1m1
n2m2 ... ( nnmn )
将λi代入 ( [δ] [M] - λi [I ] ){Y(i)}={0} 可求出n个主振型。
多个自由度体系的自由振动
结构在受迫振动时的动力响应与结构的动力特性 密切相关;另外,当用振型叠加法计算任意干扰力 作用下结构的动力响应时,往往要用到自由振动的 频率(frequency)和振型(mode)。
为此,要需要首先分析自由振动。
自振频率和振型的计算
m1
m2
mi
mn
y1(t) y2(t)
yi(t)
刚度法
其中最小的频率1 称为最低自振频率,或称
基本频率。 通常将上述每一个频率所对应的振动都称为
主振动,对应于每一个主振动的形状称为主振 型。
1)如果各质体的初速度为零,而初位移和某 一振型成比例,然后任其自然,则系统就按 这个振型作简谐自由振动,此解答就相应于 该振动的一组特解;
结构力学课件之多自由度体系的振动
y y y i点位移: yi = − fi1m1 ɺɺ1 − fi2m2 ɺɺ2 − ⋯ − fin mn ɺɺn 点位移: y y y 即: fi1m1 ɺɺ1 + fi2m2 ɺɺ2 + ⋯ + finmn ɺɺn + yi = 0
二、柔度法
同理, 同理,体系中的每一个质点都可以列出相应的动力位移方程 即用柔度法多自由度体系自由振动时的运动微分方程式。 式,即用柔度法多自由度体系自由振动时的运动微分方程式。
m2
ɺ MYɺ + KY = 0
方程特解: 方程特解:
X (1) y1 y X ( 2 ) 2 ω sin( t + ϕ ) = ⋮ ⋮ yn X ( n )
ω 即 :Y = X sin( t + ϕ )
二、柔度法
柔度法:由各质点运动的位移协调条件建立微分方程 柔度法:由各质点运动的位移协调条件建立微分方程; 按照力法的概念求解: 按照力法的概念求解: n 1 3 2 a. 确定体系的振动自 由度; 由度; b. 依次给予质点施加 f n1 f 11 f 31 f 21 一单位力。 一单位力。在此力 1 作用下, 作用下,各质点产 生的位移。 生的位移。 f n2 f 12 f 32 f 22 如质点受力: 如质点受力: y 惯性力: 惯性力: − mi ɺɺi 1
用柔度法。可分别在 、 、 用柔度法。可分别在1、2、 3点作用单位力,画出弯 点作用单位力, 点作用单位力 矩图, 矩图,利用图乘法就可以 求出各柔度系数值fij。
4m m
1
m
2
0.5m
3
(a)
4m
4m
4
多自由度振动
2
§3.1二自由度系统的自由振动
p1 p2
振幅比 振幅比
γ γ
1
A 2(1 ) a − p 12 f = = (1 ) = A1 b e − p 12
2 A 2( 2 ) a − p2 f = = (2 ) = 2 A1 b e − p2
第3章 多自由度系统振动 章
3.1 二自由度系统的自由振动
一. 二自由度系统振动研究的意义
工程中的大量实际问题是不能简化为单自由度系统, 而只能简化为多自由度系统才能充分描述其振动特性。 二自由度系统振动问题具有一定的代表性,我们通过处 理二自由度系统振动问题及实际应用来熟悉多自由度振 动系统。另一方面,二自由度振动理论在实际中广泛的 应用,因此讨论二自由度系统的振动,具有特别重要的 意义。
2 1 2 2 2 1 2 1 2 2 1 2 2
( ) ( ) x (t ) = x ( ) + x ( ) = A ( ) sin ( p t + ψ ( ) ) + A ( ) sin ( p t + ψ ( ) )
注意到实际振动形式 是两种主振动的叠加
二自由度汽车自由振动分析
二自由度汽车力学模型
2
2
= f2
2 1 2
B1 = B2 =
(e − ω ) f + bf ( p − ω )( p − ω ) (a − ω ) f + ff ( p − ω )( p − ω )
2 1 2 2 2 2 2 2 2 2 1 2 2 2 2
§3.2二自由度系统的强迫振动 二自由度系统的强迫振动
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
设频率方程具有一个二重根,即两个主振型 {Y (1) } 和{Y ( 2) }对应 的固有频率彼此相等,记为 1 2 1, 2 ,而其他频率都彼 此不同。 [K ]{Y ( S ) } 2 [M ]{Y ( s ) }
s
[K ]{Y (1) } 12,2 [M ]{Y (1) } (a)
§ 3.3 主振型的正交性和正则坐标
{Y } [K ]{Y } {Y } [M ]{Y }
(s) T (t ) 2 t (s) T (t )
转置
[ K ]T [ K ], [M ]T [M ]
{Y (t ) }T [K ]{Y (s) } s2 {Y (t ) }T [M ]{Y ( s) }
是一个与频率 1, 2 对应的主振型向量。
10
§ 3.3 主振型的正交性和正则坐标
如果两个主振型 {Y (1) }和 {Y ( 2) }彼此不正交,即
{Y (1) }T [M ]{Y ( 2) } 0
取一个由 {Y (1) } 和 {Y ( 2) } 组成的新的主振型,即
{Y
(1, 2)
} {Y } c{Y }
{Y } [ K ]{Y } 0
(t ) T (s)
相对于刚度矩阵 [K]来说,不同频率相应的主振型是彼此正 交的。 主振型第二正交条件 两个正交关系式是建立在s≠t 基础的。
2
§ 3.3 主振型的正交性和正则坐标
{Y (t ) }T [M ]{Y ( s ) } 0
对于s=t的情形,令:
{Y ( j ) }T [M ]{y} j Mj
可将任一位移按主振型展开。
7
§ 3.3 主振型的正交性和正则坐标
主振型正交的物理意义: 1)每一主振型相应的惯性力在其他主振型上不做功,其数学表 达式为 :
ys {Y ( s ) }sin(s t ), ys s 2 {Y ( s ) }sin(s t ) f s [ M ] ys [ M ]s2{Y ( s ) }sin(s t ) {Y } f s {Y } [ M ]s {Y }sin(s t ) 0
Y
1
0.163 0.569, 1
Y
2
0.924 1.227, 1
Y
3
2.760 3.342 1
解:(1)演算第一正交性。
0.163 {Y (1) }T [ M ]{Y (2) } 0.569 1
主振型正交的物理意义:
1)每一主振型相应的惯性力在其他主振型上不做功,其数学表 达式为 :
s2 {Y (t ) }T [M ]{Y ( s ) } 0
2)当一体系只按某一主振型振动时,不会激起其他主振型的 振动。 3)各个主振型都能够单独出现,彼此间线性无关。
9
§ 3.3 主振型的正交性和正则坐标
(1) ( 2)
{Y (1) }T [M ]{Y (1, 2) } {Y (1) }T [M ]{Y (1) } c{Y (1) }T [M ]{Y ( 2) }
{Y
(1,2)
}=0
{Y (1) }T [ M ]{Y (1) } c (1) T {Y } [M ]{Y ( 2) }
11
[M ]{ y(t )} [ K ]{ y(t )} 0
[Y ] [M ][Y ]{ (t )} [Y ] [ K ][Y ]{ (t )} 0
T T
17
§ 3.3 主振型的正交性和正则坐标
利用广义质量矩阵和广义刚度矩阵的定义,有
(t )} [ K * ]{ [M * ]{ (t )} 0
3
§ 3.3 主振型的正交性和正则坐标
[K ]{Y
(S )
(s) T { Y } } [M ]{Y }
2 s
(s)
左乘
{Y ( s ) }T [K ]{Y ( s ) } s2 {Y ( s ) }T [M ]{Y ( s ) }
K s {Y ( s ) }T [ K ]{Y ( s ) }
T
2 0 0 0.924 0 1 0 1.227 m 0.0006m 0 1 0 0 1 5
§ 3.3 主振型的正交性和正则坐标
同理:
{Y (1) }T [ M ]{Y (3) } 0.002m 0 {Y ( 2) }T [ M ]{Y (3) } 0.0002m 0
将{Y ( j ) }T [M ]左乘方程的两边:
i 1
n
{Y
( j) T
} [ M ]{y} i {Y ( j ) }T [ M ]{Y (i ) }
i 1
n
由主振型的正交性:
{Y } [M ]{y} j {Y } [M ]{Y } j M j
( j) T ( j) T 0
(t ) T (s)
M s {Y ( s ) }T [M ]{Y ( s ) }
K s {Y } [ K ]{Y }
(s) T (s)
Ms和Ks分别称为第s个主振型相应的广义质量 (generalized mass)和广义刚度(generalized stiffness) 每个主振型都有相应的广义质量 和广义刚度。
(1) T ( 2)
T
同理:
{Y (1) }T [ K ]{Y (3) } 0.001 k 0 {Y ( 2) }T [ K ]{Y (3) } 0.00001 k 0
6
§ 3.3 主振型的正交性和正则坐标
对任意一个位移向量{y} ,将其写成主振型的线性组合:
{ y} 1{Y (1) } 2 {Y ( 2) } n {Y ( n ) } i {Y (i ) }
M1 0 0
0 M2 0
0 0 [M * ] Mn
14
§ 3.3 主振型的正交性和正则坐标
矩阵 [ M * ] 中的非对角元素全为零,对角线的元素就是广义质 量 M 1 , M 2 , M n 对角矩阵 [ M * ] 称为广义质量矩阵。
K1 0 [Y ]T [ K ][Y ] 0
[K ]{Y (2) } 12,2 [M ]{Y ( 2) } (b)
(a) a (b) b
[K ](a{Y (1) } b{Y (2) }) 12,2 [M ](a{Y (1) } b{Y (2) })
{Y (1, 2) } a{Y (1) } b{Y ( 2) }
i (t ) K ii (t ) 0 M i
(i 1,2,, n)
利用正则变换,可以把一个n元联立方程组简化为n个独立的一 元方程,将一个具有n个自由度的结构体系的耦合振动问题简化 为n个独立的单自由度体系的振动问题,计算工作大为简化。 解耦条件: (1)线性结构 (2)[M]、[K]具有正交性
{Y (1) }和 {Y (1,2) } 就是两个彼此正交的主振型。
§ 3.3 主振型的正交性和正则坐标
由于 1, 2与其余 i (i 3,4,, n) 不相等,与 1, 2 对应的
任意一个主振型{Y (1, 2) } 都与其余频率的主振型{Y ( i ) } (i=3,4,
…,n) 彼此正交。 在具有n个自由度的体系中,即使在频率方程中出现两重根, 仍然可以选到n个主振型,使它们彼此正交。 n个自由度的体系一定有n个彼此正交的主振型。
2 s
2 t
{Y } [M ]{Y
(t ) T (s)
(t ) T
(s)
} 0
1
{Y } [K ]{Y } 0
§ 3.3 主振型的正交性和正则坐标
{Y } [M ]{Y
(t ) T (s)
} 0
振型的正交关系式(orthogonality relation) 相对于质量矩阵 [M]来说,不同频率相应的主振型是彼此正交 的。 主振型第一正交条件 另一个正交关系式:
M s {Y ( s ) }T [M ]{Y ( s ) }
K s s2 M s
可以利用广义质量Ms和广义刚度Ks计算多自由度体系的第s个自 由振频率 s 。
Ks s Ms
由广义刚度和广义质量求频率的公式。 是单自由度体系频率公式的推广。
4
§ 3.3 主振型的正交性和正则坐标
0 K2 0
0 0 [K * ] Kn
* [ K ]称为广义刚度矩阵。 为广义刚度,对角矩阵 Ki
15
§ 3.3 主振型的正交性和正则坐标
n 个自由度体系的振动方程:
[M ]{ y(t )} [ K ]{ y(t )} 0
质量矩阵[M]和刚度矩阵[K]都是对角矩阵, 方程组就是n个独立的方程,每个方程只有一个未知量。 相当于求解n个单自由度体系的振动问题。 质量矩阵[M]和刚度矩阵[K]不是对角矩阵。方程是一个耦合 方程。
(t ) T (t ) T (t ) T 2 (s)
{Y } [ M ]s {Y }sin(s t )
2 (s)
s {Y } [ M ]{Y }sin(s t ) 0
2 (t ) T (s)
s2 {Y (t ) }T [M ]{Y ( s ) } 0
8
§ 3.3 主振型的正交性和正则坐标
例:图示体系的刚度矩阵[K]和质量矩阵[M]为: m
m
2m
k 5 k 3
k
三个主振型分别如下,演算正交性。
20 5 0 k [K ] 5 8 3 15 0 3 3