模态分析理论基础
模态分析中的几个基本概念模态分析中的几个基本概念分析
模态分析中的几个基本概念物体按照某一阶固有频率振动时，物体上各个点偏离平衡位置的位移是满足一定的比例关系的，可以用一个向量表示，这个就称之为模态。
模态这个概念一般是在振动领域所用，你可以初步的理解为振动状态，我们都知道每个物体都具有自己的固有频率，在外力的激励作用下，物体会表现出不同的振动特性。
一阶模态是外力的激励频率与物体固有频率相等的时候出现的，此时物体的振动形态叫做一阶振型或主振型；二阶模态是外力的激励频率是物体固有频率的两倍时候出现，此时的振动外形叫做二阶振型，以依次类推。
一般来讲，外界激励的频率非常复杂，物体在这种复杂的外界激励下的振动反应是各阶振型的复合。
模态是结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。
这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。
有限元中模态分析的本质是求矩阵的特征值问题，所以“阶数”就是指特征值的个数。
将特征值从小到大排列就是阶次。
实际的分析对象是无限维的，所以其模态具有无穷阶。
但是对于运动起主导作用的只是前面的几阶模态，所以计算时根据需要计算前几阶的。
一个物体有很多个固有振动频率（理论上无穷多个），按照从小到大顺序，第一个就叫第一阶固有频率，依次类推。
所以模态的阶数就是对应的固有频率的阶数。
振型是指体系的一种固有的特性。
它与固有频率相对应，即为对应固有频率体系自身振动的形态。
每一阶固有频率都对应一种振型。
振型与体系实际的振动形态不一定相同。
振型对应于频率而言，一个固有频率对应于一个振型。
按照频率从低到高的排列，来说第一振型，第二振型等等。
此处的振型就是指在该固有频率下结构的振动形态，频率越高则振动周期越小。
在实验中，我们就是通过用一定的频率对结构进行激振，观测相应点的位移状况，当观测点的位移达到最大时，此时频率即为固有频率。
实际结构的振动形态并不是一个规则的形状，而是各阶振型相叠加的结果。
固有频率也称为自然频率( natural frequency)。
模态分析的基础理论
模态分析的基础理论模态分析是一种研究系统中不同模式的分布、生成和演化规律的方法。
在这个理论中，模态是指系统中不同状态或形式的存在形式，例如质量分数、温度、湿度等。
模态分析的基础理论包括概率论、统计学和模态分析技术等。
概率论是模态分析的基础之一、它研究随机事件的发生概率和规律。
在模态分析中，我们可以利用概率论来描述不同模态出现的概率分布，并通过分析系统中的模式，得出不同模态的生成规律。
通过概率论的方法，我们可以预测不同模态的变化趋势，从而指导系统的优化设计和运行管理。
统计学也是模态分析的基础理论之一、统计学研究如何收集、处理、分析和解释数据，通过对大量数据的统计分析，揭示数据背后的规律和趋势。
模态分析中，统计学的方法可以用于分析模态数据的分布情况，寻找模态之间的相关性和影响因素，并建立相应的模型来预测和优化系统的运行情况。
在模态分析技术方面，主要包括聚类分析、主成分分析和模态分析方法等。
聚类分析是一种将相似的对象分组的方法，通过对模态数据进行聚类分析，我们可以将相似的模态归为一类，从而描述系统中的不同模态分布情况。
主成分分析是一种降维技术，它可以将高维的模态数据降低到低维，并保留大部分信息。
这可以帮助我们更好地理解系统模态之间的关系和重要性。
模态分析方法包括有限元模态分析、频响函数法和模态参数识别等。
通过这些方法，我们可以对系统的模态进行分析，包括振型、频率和阻尼等，并找出模态的摄动源和分布规律。
模态分析的基础理论对于理解和优化系统具有重要意义。
通过对模态的分析和研究，我们可以了解系统的特性和不同模态之间的关系，从而指导系统的设计和运行。
同时，模态分析也可以帮助我们发现和解决系统中存在的问题，提高系统的稳定性和可靠性。
因此，深入理解和应用模态分析的基础理论对于各个领域的研究和实践具有重要价值。
第3章 实验模态分析的基本理论
实验模态分析第三章：实验模态分析的基本理论振动系统的特性可以用模态来描述:固有频率、固有振型(主振型)、模态质量、模态刚度和模态阻尼等。
建立用模态参数表示的振动系统的运动方程并确定其模态参数的过程使称为模态分析。
—种理解可以认为，振动系统的物理模型、物理参数和以物理参数表示的运动方程都是已知的，引入模态参数、建立模态方程的目的是为了简化计算，解除方程耦合，缩减自由度。
另一种理解可以认为，通过对实际结构的振动测试，识别振动系统的模态参数，从而建立起系统的以模态参数表示的运动方程，供各种工程计算应用。
试验模态分析指的是后一种过程,即通过振动测试(称模态试验),识别模态参数,建立以模态参数表示的运动方程这样一个过程。
1 多自由度系统振动基础回顾&&&++=M x C x K x f t []{}[]{}[]{}{()} 2实模态理论一个n 自由度线性定常振动系统，其运动方程可以如下表示:现对两端作付氏变换得:[]{}[]{}[]{}{()}M x C xK x f t ++=&&&2([][][]){()}{()}M j C K X F ωωωω−++=式中和分别是x(t)和F(t)的付氏变换,并有()X ω()F ω()()j t X x t e dt ωω+∞−−∞=∫()()j t F f t e dtωω+∞−−∞=∫(){()}{()}Z X F ωωω=111212122212()()()()()()()()()()n n n n nn Z Z Z Z Z Z Z Z Z Z ωωωωωωωωωω⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦L L L L L L L 1()[()]{()}{()}{()}X Z F H F ωωωωω−==2[][][]K M j C ωω=−+阻抗矩阵中各元素值无法在实际振动测试中获得，因为人们不可能在实际结构上固定其它坐标，令其不动，仪留下J坐标，待其作出响应；也不可能仅使某个坐标运动，在其余坐标上测量力。
_模态分析理论基础
IVE
Institute of Vibration Engineering, Northwestern Polytechnical University, China
Iration Engineering, Northwestern Polytechnical University, China
有限元分析软件(如ANSYS、NASTRAN、SAP、MAC等)在结 构设计中被普遍采用,但在设计中,由于计算模型和实
际结构的误差,而且受到边界条件很难准 确确定的影响,特别是结构的形状和动态特性很复杂时，
IVE
Institute of Vibration Engineering, Northwestern Polytechnical University, China
e. 为结构动力学优化设计提供目标函数或约束条件
动力学设计，即对主要承受动载荷而动特性又至关重要的结构，以 动态特性指标作为设计准则，对结构进行优化设计。它既可在常规静力 设计的结构上，运用优化技术，对结构的元件进行结构动力修改；也可 从满足结构动态性能指标出发，综合考虑其它因素来确定结构的形状， 乃至结构的拓扑（布局设计、开孔、增删元件）。动力学优化设计就是 在结构总体设计阶段就应对结构的模态参数提出要求，避免事后修补影 响全局。
•解的形式（s为复数）及拉氏 变换： x Xest (ms2 cs k ) x(s) f (s)
IVE
Institute of Vibration Engineering, Northwestern Polytechnical University, China
实验模态分析基本理论
机械式（Mechanical）：频带窄（10-100Hz）、行程一般数毫米、 噪音大、位移和波形控制不精确；
3.频响函数的测量
(2) 冲击激励（Impact Excitation）： 力锤（Hammer），适用于小阻尼线性结构。还有夯锤、 落锤、摆式冲击锤、小火箭等；
(3) 脉动等（Enviromental Excitation）： 利用大地、地震、人工爆炸模拟地震风等条件引起的振动。
3.频响函数的测量
1.3 测量系统(Measurement System)： 传感器+放大器（Transducer+Amplifier） ICP传感器（内装IC放大电路） （Integrated Circuit Piezoelectric）
第8章 实验模态分析初步
1.结构动特性的建模方法
理论建模：动力学仿真分析
耦合的力平衡方程 Mx(t) Cx(t) Kx(t) f (t)
通过有限元离散模型，采用模态分析理论、特征值求解技术和矩阵变换 技术获得结构的动力学特性：频率、阻尼和振型；建立结构完整模态模 型（一组解耦的方程组）、系统的传递特性（传递函数或频响函数）， 进一步建立结构的频域响应计算模型和时域响应计算模型。
3.频响函数的测量
4 测点布置与激振点的选择 测点布置 1.能够较好地反映结构物的构型 2.能够充分显示结构的模态振型 例：一个梁单元无法求解简支梁的10个模态。计算上一般要求至少 20个单元，计算出的20阶模态，只有前10阶准确。
激励点 应避开节点节线。多点激励进行校核。
激励力的选择 在不破坏试件的情况下，尽可能大的激励力，有助于提高信噪比。 不同大小的激励力，可以定性考查结构非线性的程度
模态分析理论基础
有限元分析软件(如ANSYS、NASTRAN、SAP、MAC等)在结
构设计中被普遍采用,但在设计中,由于计算模型和实
际结构的误差,而且受到边界条件很难准
确确定的影响,特别是结构的形状和动态特性很复杂时，
有限元简化模型和计算的误差较大。通过对结构进行实验模态分 析，可以正确确定其动态特性，并利用动态实验结果修改有限元 模型，从而保证了在结构响应、寿命预计、可靠性分析、振动与 噪声控制分析与预估以及优化设计时获得有效而正确的结果。
•传递函数和频率响应函数
H(s)m2s(11jg)k
H()m21(1jg)k
（1＋jg）k — 复刚度
–用实部和虚部表示
H ()1 k (1 1 2 )22 g2j(1 2)g 2g2
与粘性阻尼系统相比频响函数形式相同 g和2 相互置换即可得各自表达式
位移、速度和加速度传递函数
Hd (s)
e. 为结构动力学优化设计提供目标函数或约束条件
动力学设计，即对主要承受动载荷而动特性又至关重要的结构，以 动态特性指标作为设计准则，对结构进行优化设计。它既可在常规静力 设计的结构上，运用优化技术，对结构的元件进行结构动力修改；也可 从满足结构动态性能指标出发，综合考虑其它因素来确定结构的形状， 乃至结构的拓扑（布局设计、开孔、增删元件）。动力学优化设计就是 在结构总体设计阶段就应对结构的模态参数提出要求，避免事后修补影 响全局。
x(s) f (s)
Hv(s)
v(s) f (s)
Ha(s)
a(s) f (s)
• 位移、速度和加速度频率响应函数
()
x() f ()
Hv()
v() f ()
• 三者之间的关系
Ha()
模态分析理论基础
点，有图可知节点并不唯一，而且修改前后节点的位置未变。
对应尽可能避开结构振动的节点，以免给测量带来误差。
４．４试验模态分析试验模态分析的目的是为了验证理论模态分析的正确性的基础上进行深入研究奠定基础。
４．４．１试验模态分析的理论基础阻１所以在进行模态实验为在理论模态分析在物理坐标下，描述Ｎ自由度离散振动系统的运动微分方程为阻】耕＋【ｃ】扛｝＋医】Ｍ＝沙｝（４．２）式中：【Ｍ］——质量矩阵（对称且正定），Ｍ∈Ｒ～，【Ｃ】——阻尼矩阵，Ｃ∈Ｒ“”，晖】——刚度矩阵（对称且正定或半正定），Ｋ∈Ｒ“”，｛ｘ），｛卦，｛封——Ｎ维位移、速度和加速度响应向量，｛厂（ｒ））——＿Ｎ维激振力向量。
设系统的初始状态为零，对式（４．２）两边进行拉普拉斯变换可得（［Ｍｌｓ２“Ｃ］ｓ＋【Ｋ］）｛Ｘ０））＝【Ｚ（ｓ）］｛工０））＝｛Ｆ０））式中的矩阵【Ｚ（ｓ）］－（［Ｍ］ｓ２＋［ｃ］ｓ＋［Ｋ】）反映了系统的动态特性，称为系统动态矩阵或广义阻抗矩阵，其逆阵［日（５）】＝［Ｚ（ｓ）】～＝（【Ｍ］ｓ２＋【Ｃ］ｓ＋［Ｋ］）。
１称为广义导纳矩阵，也就是传递函数矩阵。
由式（２．２）可知｛ｘ（Ｊ））＿【日０）】（Ｆ（Ｊ）｝在上式中．令Ｓ＝ｊｏＪ，即可得到系统在频域内输出和输入的关系式｛并（国）｝＝【日（脚）】（Ｆ（国））（４．３）（４．４）（４．５）（４．６）（４．７）式中［Ｈ（ｃｏ）】为频率响应函数矩阵。
［Ｈ（∞）】矩阵中第ｆ行＿，列的元素％（叻２篇（４８）表示仅在』坐标激振（其余坐标激振力为零）时，ｉ坐标的响应与激振力之比。
在式（４．４）中令Ｓ＝＿，∞，可得阻抗矩阵［ｚ（∞）】＝（［Ｋ］一曲２【吖］）＋ｊｃｏ［Ｃ］（４．９）它和导纳矩阵有类似式（４．５）的关系［日（珊）］＝［ｚ（国）】～＝｛（【置卜。
２［＾卅）＋ｊｃｏ［Ｃ】｝１（４．１０）对于一般机械、结构，假设矩阵［ｃ］也对称，这样矩阵【ｚ（∞）】对称，频率响应函数矩阵［日＠）］也对称，故有ｑ（脚）＝ＨⅣ（０３）（４．１１）上式反映了机械、结构频率响应有互易性，可作为频率响应测试精度的一项重要检验手段。
模态分析的基础理论
运动微分方程
单自由度系统无阻尼自由振动是简谐振动
2π
m
T 2π
n
k
fn
1 T
n
2π
1 2π
k m
能量关系
mx dx kx dx 0 dt dt
意义：惯性力的功率Fm与弹性力的功率Fs之和为零
d dt
1 2
mx2
1 2
kx 2
0
ET
1 mx2 2
单自由度系统
自由振动 简谐振动 非周期强迫振动
自由振动
振动系统在初始激励下或外加激励消失后的 运动状态。
自由振动时系统不受外界激励的影响，其运 动时的能量来自于初始时刻弹性元件和惯性 元件中存储的能量。
振动规律完全取决于初始时刻存储的能量和 系统本身的性质。
运动微分方程
•使该矢量以等角速度在复平面内旋转（复数旋转矢量）
虚轴
ei x cos i sin
P A
t
z Acost i sint Aeit
实轴
y Asint Im z Im Aeit
运动学
速度、加速度的复数表示
位移 x Aeit
速度 x d Aeit iAAeeiitt / 2
2.0
0.5 和 0.7 临 界 阻 尼 比 无
c/cc=0
抛物线
阻尼曲线更接近理想加
1.5
速度计曲线
c/cc=0.5
1.0
c/cc=0.7
0.5
0 0 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00
模态分析基本理论
+ +
(C1 (C2
+ +
C2 C3
)x&1(t) - C2x&2 (t) )x&2 (t) - C2x&1(t)
+ +
(K1 (K2
+ K2 )x1(t) - K2x2 (t) + K3 )x2 (t) - K2x1(t)
= =
f1 (t) f2 (t)
第三节 多自由度振动系统举例
一 系统方程
写出矩阵形式：
\
eλ1 t
\
e∧t
\
=
O
0
eλN t
eλ*1 t
0
O
e
λ*N
t
第四节 多自由度系统相关模态概念
一 无阻尼系统
阻尼矩阵[C]为零矩阵的系统
系统阻尼因子σ r = 0 ，全为纯虚数极点 λ1 = jω,L , λ*N = − jωN
系统方程：
P2 + Pα Pβ+1
[M]
+
[K
]{x}
=
{0}
比例阻尼系统频响函数
\
或 [H( jω)] = [ψ]
[ ] ∑ { } { } H(jω)
=
N
j2ωr Qr
ψ
r
ψ
T
r
r =1
(σ
2
r
+
ω
2
r
-
ω
2
)
-
2σ
r
模态分析理论
e t
sin dt
就是脉冲响应函数。
很容易证明频响函数和脉冲响应函数是一对傅氏变换对：
H () Fh(t)
（1） 简谐激励
结构在简谐激励下的稳态响应也是同频率的简谐振动。但有相位差。
f (t) Fe j(t ) x(t) Xe j(t )
H（） X e j( )
F
工程中，应变常常是非常重要的，而且易于测量。应变片体积小、质量小、成分低，对试验结
结构动力修改
模态分析的目的是了解系统的动态特性。在已知结构动态特性参数后，我们应该寻求改进系统动态 特性的方法。 有两种情况： 1） 由于制造和设计原因，不得不对现有结构进行局部修改。
word 格式-可编辑-感谢下载支持
机械模态分析理论基础
假设：系统是线性、定常与稳定的线性时不变系统
线性：描述系统振动的微分方程为线性方程，其响应对激励具有叠加性;
定常：振动系统的动态特性（如质量、阻尼、刚度等）不随时间变化，即具有频率保持性;如系统受简谐 激励-响应的频率必定与激励一致。 稳定：系统对有限激励必将产生一个有限响应，即系统满足傅氏变换和拉氏变换的条件。 振动系统分类：
word 格式-可编辑-感谢下载支持
ˆ
2 fx
()
1
GMM G ff ()
1 1
GNN Gxx ( )
输入存在噪声，会使估计的频响函数偏小；
输出存在噪声，会使估计的频响函数偏大；
还可用下面一些估计方法：
Hˆ 3 ()
Hˆ 1 ( )
2
Hˆ 2 ()
Hˆ 4 () Hˆ1() Hˆ 2 ()
K s2M φs 0
右乘 φs ，得到：
φsT KT s2MT φr 0
模态分析理论基础
损耗因子
E 2U
U——最大势能
根据能量等效原则（一个周期内等效粘性阻尼与结构阻尼耗散能量相等 ΔW=ΔE），求得结构阻尼的等效粘性阻尼系数ce ：
ce g rk 2
或
ce
k
注意与频率有 关，非常数
其中 g
rk k 称为结构阻尼系数，具有刚度量纲。 2
损耗因子常称为结构阻尼比。
5
2014/12/18
表1.3-1 单自由度结构阻尼系统频响函数各种图形及数字特征（续）
(b) 相频特性曲线
半功率带宽
O 4 2 4
A B
① 拐点M （位移谐振点）：
A M B
=1 ( =D =0 )
固有频率
0 D
-
② 半功率点A、B： A ，B 结构阻尼比
(1.3-1)
其中频率比（无量纲频率） B. 频响函数的极坐标表达式：
0
H H e j
(1.3-2) (1.3-3) (1.3-4)
8
其中：幅频特性 相频特性
H
1 k
1
2
1
2
tg 1
2 1
i. ii.
各种特性曲线不象结构阻尼系统那样具有较简单的特征； 粘性阻尼系统具有三种不相等的谐振频率：位移谐振频率D、阻尼谐 振频率d和无阻尼谐振频率0，它们出现在各种曲线的不同特征点上， 具有如下关系：
D d 0
iii.
粘性阻尼系统的Nyquist图也不在是一个圆，而是一个近似桃形的图形。 不过，在小阻尼情形下，使用Nyquist图作参数识别时仍可将其视为圆 来处理。
模态分析的基础理论-PPT精品文档109页
k
m
c x
kx c·x
m F0 cos t
简谐强迫振动
系数
B
2
x
2 0
x0
n d
x0
tan 1 x0 n x0 d x0
X
A
1
(
n
)
2
2
2
n
2
ET
U1kA2 2
12(x02x02n2)
ET UE
Rayleigh商 动能系数
能量关系
T1mA2 2
12mxm 2ax
n2
k m
Umax T
阻尼自由振动
方程
mxcxkx 0 x(0) x0, x0(0) 0
x2nxn2x 0
自激振动：输电线的舞动 1940年美国塔可马(Tacoma Narrows)吊桥在中速
风载作用下，因桥身发生扭转振动和上下振动造 成坍塌事故 1972年日本海南的一台66×104kW汽轮发电机组， 在试车过程中发生异常振动而全机毁坏； 步兵在操练时，不能正步通过桥梁，以防发生共 振现象造成桥梁坍塌
x ( t) e n t( c 1 c o sd t c 2 s ind t)
x (t)X e n tco s(dt)
c 1 x 0 ,c 2 (xn x 0)/ d
阻尼自由振动
对数衰减率
x1 x2
X Xeenntt12ccooss((ddtt11)
单自由度系统
自由振动 简谐振动 非周期强迫振动
模态分析的基础理论
通常用物理量：
质量M，刚度K，阻尼C，和外界激励F表示。
x
1
x k
2
x
1
x c
2
振动分类
按系统分： 线性系统和非线性系统 离散系统和连续系统 确定性系统和随机系统 按激励分： 自由振动 受迫振动 自激振动 参数共振
振动分类
按响应分： 简谐振动 周期振动 非周期振动 随机振动
按自由度分： 单自由度振动 多自由度振动 连续体振动
运动学
一、简谐运动
按时间的正弦函数(或余弦函数)所作的振动
x Asin t
振幅 相位 圆频率 初相位
运动学 简谐振动的速度和加速度
位移
x Asin t x A cos t x A sin t
颤振：大气紊流和其他振源都会使飞机等飞行器 产生振动（舒适性，机载仪表） 自激振动：输电线的舞动 1940年美国塔可马(Tacoma Narrows)吊桥在中速 风载作用下，因桥身发生扭转振动和上下振动造 成坍塌事故 1972年日本海南的一台66×104kW汽轮发电机组， 在试车过程中发生异常振动而全机毁坏； 步兵在操练时，不能正步通过桥梁，以防发生共 振现象造成桥梁坍塌
机械振动的积极作用
共振放大 利用颗粒的振动进行清洗，抛光，零件去毛刺； 利用振动减小零部件之间的摩擦阻力和间隙
阀体 阀芯 电磁铁
学习机械振动的意义
1. 2. 3.
4.
进行结构动强度设计的需要 消除有害的振动 利用振动有利的一面 是学好相关知识的基础
离散系统的基本元件
第一章模态分析理论基础
共振频率点
ds max d 1
• 粘滞阻尼系统
– Nyquist图
2
2
[H
R
( )]2
(H
I
( ))2
1
4k
1
4k
» 特点
»桃子形，阻尼比越小
轨迹圆越大
» ( 是变的，所以不是圆 )
在固有频率附近，曲线 接近圆，仍可利用圆
的特性
第20页/共60页
速度与加速度频响函数特性曲线
• 关系回顾
HR 1, 2
(
)
4k
1 (1
)
2
1
g
2
半功率带宽反映阻尼大小 阻尼越大，半功率带宽
越大，反之亦然
第17页/共60页
• 虚频图
• •
H
I
( )
g
k[(1 2 )2
（结构阻尼） （g粘2 ] 性阻尼）
• 以H结I构(阻) 尼k[为(1例：2 )22(2 )2 ]
– 系统共振时虚部达到最大值
– 系统共振时实部为零
m1
机架线
第30页/共60页
• 一般多自由度约束系统
机架线
– N自由度约束系统有N个共振频率,（N－1）个反共振频率 – 对原点函数共振反共振交替出现 – 对跨点频响函数无此规律 – 一般两个距离远的跨点出现反共振的机会比较近的跨点少
第31页/共60页
– 自由系统
• 两自由度系统运动方程（无阻尼）
第7页/共60页
单自由度系统频响函数分析
粘性阻尼系统
•阻尼力（与振动速度成正比）：
•强迫fd振动方c程x 及其解
..
.
m x•解c的x形式k（xs为复f 数）及拉氏变换：
第10讲：模态分析
在这种状态下，系统具有六个刚体模态， 即三个平移模态和三个转动模态。前者由结 构的质量所确定，后者由三个转动惯量所确 定。刚体模态所对应的固有频率为零。这种 状态下，可获得系统无边界约束时的模态。 试验结构若与其他结构联接时，则可以通过 增加联接约束，通过模态综合分析获得总体 结构的模态。
实际应用中，除太空零重力环境下和飞 行器俯冲失重情况下可实现被测结构的自由 状态外，地面试验中常采用的气悬浮和磁悬 浮方法都只能实现平面三个自由度的自由状 态，而且这种自由支承的制造精度和费用很 高，目前仍未广泛采用。对于大多数工程实 际情形，所谓自由状态还是要通过某种支承 来实现的。
H11 H12
H
21
H 22
H
n1
Hn2
H1n
1r
H
2
n
H
nn
m r 1
kr
1
2mr
jcr
r
2nrr
[1r
,
2r ,
nr ]
式中n为测点数，m为所取模态数。
一般情况下n>m，高阶模态作剩余影响 处理。由上式可见，为了取得全部模态信息， 仅需测量频响函数矩阵中的一行或一列就够 了。为了测一列，可采用在结构上选定适当 的一点，作单点固定激振，而在所有测量点 (含激振点)依次测量响应，称为SISO测量法; 或在多点同时测量响应，称为SIMO法。
例如，把一个齿轮放在一个很软的泡沫 塑料上，进行锤击试验时，可以认为齿轮处 于自由状态；用长柔性索悬挂结构时，在水 平方向激振也可以认为在水平方面处于自由 状态；而用很柔的橡胶绳或弹簧将结构吊起， 则可以认为结构在垂直方向处于自由状态。
尽管悬索很长或橡胶的柔性很好，但既 有支承就只能是自由状态的一种近似。其结 果是：刚体运动模态所对应的固有频率不再 是零。悬挂较好的情况应能保证最高的“刚 体模态”(结构不产生任何弹性变形的运动模 态)频率比最低的结构变形模态要低得多，至 少应能保证小于0.1~0.2。
模态分析理论
精心整理模态分析指的是以振动理论为基础、以模态参数为目标的分析方法。
首先建立结构的物理参数模型，即以质量、阻尼、刚度为参数的关于位移的振动微分方程；其次是研究其特征值问题，求得特征对（特征值和特征矢量），进而得到模态参数模型，即系统的模态频率、模态22¨330m 0z k 2k k z 000m 0k k z 0z +--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦（9） 定义主振型由于是无阻尼系统，因此系统守恒，系统存在振动主振型。
主振型意味着各物理坐标振动的相位角不是同相（相差0o ）就是反相位（相差180o ），即同时达到平衡位置和最大位置。
主振型定义如下：()i i j ωt+i i sin ωt+=Im(e )φφi mi mi z =z z （10）其中为第i 阶频率下，各自有度的位移矢量，为第i 个特征矢量，表示第i 阶固有频率下的振型，i ω为第i 阶频率下的第i 个特征值，i φ为（去除项化简得以矩阵的形式展开得：2i 2i mi 2i k-ωm -k 0-k 2k-ωm -k z =00-k k-ωm ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦（15）有非零解，则2i 2i 2i k-ωm -k 0-k 2k-ωm -k =00-k k-ωm ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦（16）即()234222ω-m ω+4km ω-3k m =0（17）阶固有频率，每一个特征根对应一个特征矢量，表示对应模态下该由式3i i 21=z k 如果设定了1z 值，则就可以求出三个特征根值下，2z 和3z 相对于1z 的位移。
假设m=k=1， 一阶模态，1ω=0：21z =1z ，31z =1z ，即；二阶模态，223kω=m ：21z=0z，31z=-1z，即；三阶模态，23kω=m ：21z=-2z，31z=1z，即。
运动方程的解耦图错误！未指定顺序。
运动方程解耦过程在进行坐标变换之前需对刚度矩阵和质量矩阵进行归一化。
模态分析理论
模态分析指的是以振动理论为基础、以模态参数为目标的分析方法。
首先建立结构的物理参数模型，即以质量、阻尼、刚度为参数的关于位移的振动微分方程；其次是研究其特征值问题，求得特征对（特征值和特征矢量），进而得到模态参数模型，即系统的模态频率、模态矢量、模态阻尼比、模态质量、模态刚度等参数。
特征根问题以图3所示的三自由度无阻尼系统为例，设123m =m =m =m ，123k =k =k =k ，图 1 三自由度系统其齐次运动方程为：（8）其中分别为系统的质量矩阵和刚度矩阵，123m 00m 00m=0m 0=0m 000m 00m ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦，11212221k -k 0k -k 0k=-k k +k -k =-k 2k -k 0-k k 0-k k ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦，则运动方程展开式为： ¨11¨22¨33z m 00k k 0z 00m 0z k 2k k z 000m 0k k z 0z ⎡⎤⎢⎥-⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥+--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦（9）定义主振型由于是无阻尼系统，因此系统守恒，系统存在振动主振型。
主振型意味着各物理坐标振动的相位角不是同相（相差0o ）就是反相位（相差180o ），即同时达到平衡位置和最大位置。
主振型定义如下：()i ij ωt+i i sin ωt+=Im(e)φφi mi mi z =z z （10）其中为第i 阶频率下，各自有度的位移矢量，为第i 个特征矢量，表示第i 阶固有频率下的振型，i ω为第i 阶频率下的第i 个特征值，i φ为初始相位。
对于三自由度系统，在第i 阶频率下，等式可以写成1m1i 2m2i i i 3m3i z z z =z sin(ωt+)z z φ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦（11）mki z 表示第k 个自由度在第i 阶模态下的模态矩阵。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
点，有图可知节点并不唯一，而且修改前后节点的位置未变。
对应尽可能避开结构振动的节点，以免给测量带来误差。
４．４试验模态分析
试验模态分析的目的是为了验证理论模态分析的正确性的基础上进行深入研究奠定基础。
４．４．１试验模态分析的理论基础阻１所以在进行模态实验为在理论模态分析
在物理坐标下，描述Ｎ自由度离散振动系统的运动微分方程为
阻】耕＋【ｃ】扛｝＋医】Ｍ＝沙｝（４．２）式中：【Ｍ］——质量矩阵（对称且正定），Ｍ∈Ｒ～，
【Ｃ】——阻尼矩阵，Ｃ∈Ｒ“”，
晖】——刚度矩阵（对称且正定或半正定），Ｋ∈Ｒ“”，
｛ｘ），｛卦，｛封——Ｎ维位移、速度和加速度响应向量，
｛厂（ｒ））——＿Ｎ维激振力向量。
设系统的初始状态为零，对式（４．２）两边进行拉普拉斯变换可得
（［Ｍｌｓ２“Ｃ］ｓ＋【Ｋ］）｛Ｘ０））＝【Ｚ（ｓ）］｛工０））＝｛Ｆ０））式中的矩阵
【Ｚ（ｓ）］－（［Ｍ］ｓ２＋［ｃ］ｓ＋［Ｋ】）
反映了系统的动态特性，称为系统动态矩阵或广义阻抗矩阵，其逆阵
［日（５）】＝［Ｚ（ｓ）】～＝（【Ｍ］ｓ２＋【Ｃ］ｓ＋［Ｋ］）。
１称为广义导纳矩阵，也就是传递函数矩阵。
由式（２．２）可知
｛ｘ（Ｊ））＿【日０）】（Ｆ（Ｊ）｝
在上式中．令Ｓ＝ｊｏＪ，即可得到系统在频域内输出和输入的关系式
｛并（国）｝＝【日（脚）】（Ｆ（国））（４．３）（４．４）（４．５）（４．６）（４．７）
式中［Ｈ（ｃｏ）】为频率响应函数矩阵。
［Ｈ（∞）】矩阵中第ｆ行＿，列的元素
％（叻２篇（４８）表示仅在』坐标激振（其余坐标激振力为零）时，ｉ坐标的响应与激振力之比。
在式（４．４）中令Ｓ＝＿，∞，可得阻抗矩阵
［ｚ（∞）】＝（［Ｋ］一曲２【吖］）＋ｊｃｏ［Ｃ］（４．９）它和导纳矩阵有类似式（４．５）的关系
［日（珊）］＝［ｚ（国）】～＝｛（【置卜。
２［＾卅）＋ｊｃｏ［Ｃ】｝１（４．１０）对于一般机械、结构，假设矩阵［ｃ］也对称，这样矩阵【ｚ（∞）】对称，频率响应函数矩阵［日＠）］也对称，故有
ｑ（脚）＝ＨⅣ（０３）（４．１１）上式反映了机械、结构频率响应有互易性，可作为频率响应测试精度的一项重要检验手段。
利用振型矩阵的加权正交条件，即
蛾ｎ圳¨：｛ｏ，Ｈ９蛾ｎ州。
，）－｛？’Ｈ９（４．１２）１２
｛中。
｝７【＾卅｛中，）＝｛’１｛ｍ。
｝。
【置】｛ｏ，）＝｛，’１（４．。
【“，，，２ｑ１Ｉ疗，，ｒ２９
式中：｛ｏ。
）、砷，）＿一第ｇ、，阶特征向量（固有振型），
ｍ．、ｋｒ——第，阶模态质量和模态刚度。
并设矩阵【Ｃ】也可由振型矩阵【叫对角化，那么可以对阻抗矩阵进行如下变换
＝［中］。
【ｏ］７“［置卜出２［盯】）＋＿，珊［Ｃ］）【∞］【中】－１
＝ｃ。
，４Ⅱ、七一、］一∞２［、ｍ，、］＋－，国［、ｃ一、］｝ｃ西，～。
．。
，，
『＼］
＝【中】。
ｌｏ｜【叫“
式中：ｚ，＝（七，一∞２ｍ，）＋，珊ｃ，，ｃ，是第ｒ阶模态阻尼。
由（４．９）、（４．１０）式可以得到用模态参数表示的频率响应函数矩阵为
江苏大学硕士学位论文
［Ｈ（珊）】矩阵中第ｉ行，列的元素为（４．１４）
吲妒喜而暑‰㈣式中：功：：生——第，阶模态频率，
ｍ．
善，＝＝＿三Ｌ——第ｒ阶模态阻尼比
２ｍ，∞，
｛中，｝——第ｒ阶模态振型。
不难发现，Ｎ自由度系统的频率响应等于Ｎ个单自由度系统的频率响应的线性迭加。
为了确定全部模态参数曲，、善，、｛中，）ｐ＝ｌ…２．．Ｎ），实际上只需测量频率响应函数矩阵的一列或一行就行了。
４．４．２模态参数识别
由ｆ４．１５）式可知，当珊趋近于某阶模态的固有频率时，该模态起主导作用，称为主模态。
在主模态附近，其它模态影响较小。
若模态密度不很大，各阶模态比较远离，其余模态的频率响应函数值在该主模态附近很小，且曲线比较平坦，即几乎不随频率而变化，因此其余模态的影响可用一复常数，ｔ来表示，对第，阶模态（４．１５）式可近似表示成
日”（珊）。
而弭ｉ１再丽丽＋日ｃ
＝去［矿豢‰一，矿考‰Ｊ（４峋
２而Ｉ五【矿万话历Ｆ叫瓦‘再丐蕊万矛ｊ¨’”’＋Ｈ：＋ｊＨ２
Ｈ，＠）的实部和虚部可表示如下
蹦咖去ｌ矿≤岳］＋彬
㈣。
去［矿耘卜根据上面两式可以作出实频图４．１９和虚频图４．２０。
Ｏ
Ｏ
图４．１９第，阶模态频率响应函数的实频图
图４．２０第ｒ阶模态频率响应函数的复频图下面就可进行模态参数的确定口］，［３Ｈ。
１）固有频率的确定（４．１７）（４．１８）
从实频图可以看出，固有频率就是实频曲线与剩余惯性直线日？＠）＝Ｈｆ的交点所对应的０９值。
另外一种方法是从虚频图来确定，此时珊＝∞，，正好对应虚频曲线的峰值，因为峰值较尖，所以更容易确定固有频率。
２）阻尼比的确定
阻尼比系数可以由半功率带宽Ａｃｏ来确定，由图４．１９的实频曲线对６９求导可得
ＡＣａ＝钆一０９。
（４．１９）对于模态阻尼系统，其阻尼比系数为
善，：．竺：掣（４．２０）
¨
２珊，２∞，、
３）确定模态振型
由（４．１８）式知，对主模态而言（不计剩余模态），当∞＝珊，时，
Ｈ：‘甜２ｑ’２霸－－１‘４’２１’
分别测出车身上各点的日：∞）ｔ『－１…２．．１６１）值，则可得到第，阶摸态的振型系数向量：
｛中，）＿扣：洄＝（Ｏｒ）｝
＝归ｊ（出＝ｑ）日二（珊＝（－Ｏｒ）’・・日ｊｓ，（ｃ０２啡）ｊ
一南‰％…‰｝
Ｈ．２２）对于第，阶模态，当采用逐点激励，单点响应时，一—乓为常数。
因此
２ｍ，鲁山，
ｐ，｝＝溉（ｃｏ＝（Ｏｒ）｝即可代表模态振型。
因为振型只反映振动形式，与振动大小无关。
故常取归一化后的振型向量，在此归一化的方法是对响应点归一化，即取中。
＝１・那么
ｐ，）＝砩洄＝ｃｏｔ）｝
＝一去‰中ｎ…叫＠２３’
４１模态质量和刚度的确定
当０９＝ｑ时，考虑，点激励，ｉ点响应的频率响应函数，其虚部可以由（４．２１）式近似表示
３７—————————＿－＿＿ｌ－———————＿———＿＿——＿——————＿＿－＿＿－－————一一
咸‘ｍ２ｑ’２丽ｍｌ㈣１
中。
中。
‘
２ｍ，毒ｑ２
若取响应点ｉ的频率响应函数，且对响应点归一化，则车架结构阻尼系统的第ｒ阶模态质量ｍ，可以表示为
卅ｒ２币丽耘一２５’应该指出，模态质量的数值与归一化的方法有关，不同的归一化方法将得到不同的模态质量。
第，阶模态刚度可由下式求出
ｋ，＝ｍ，ｃｏ；（４．２６）至此已求出全部的模态参数。
４．４．３白车身模态试验
由上可知，传递函数矩阵的一列或一行都含有模态向量．通过传递函数的一列或一行识别复模态矢量和极点位置，进而可以确定模态参数。
本文采用锤击法依次单点激励、单点测响应的试验方案．对ＮＪ６４００轿车的白车身进行了模态试验分析。
在试验过程中，输入信号由型号为５２１０的力锤产生，力传感器的型号为ＹＤＬ．１２．５，响应信号由型号为ＢＫ４３２１的三向加速度传感器拾取，激励和响应信号经电荷放大器放大后输入ＳＤ３８０动态分析仪，输出频率响应函数和相干函数曲线。
这些曲线可以硬拷贝输出，也可经计算机系统中的ＳＴＡＲ５．７结构分析软件分析后，得出车身结构的各阶固有频率和振型。
由于汽车所受边界非完全自由而受弱约束，但它所受的约束相对结构自身的刚度来说仍小得多，因此采用自由支撑也是适当的。
本试验采用气垫支撑，相当于给结构增加了柔软约束，刚体模态频率不再是零。
经测量，气垫支撑的固有频率为５Ｈｚ左右，远低于白车身的第一阶固有频率。
所以它对车身模态分析的影响很小，
可以忽略不计。
１．测试系统框图如图４．２ｌ所示。
在测试过程中，由于车身上有些点不方便固定三向加速度传感器，因此采用拾振的加速度传感器固定不动．激励的力锤移动的方法进行测量。