那么:
Hjr
、2
2-3
2-4
2、模态检测
模态检测是根据频域中的模态模型对复模态
〔或实模态〕向量进行局
部估计的一种单自由度方法。
在hj
Qr
j
中略去剩余
项那么单个频响函数在r处的值近似为:
Qr1rjr
Htjjr-——
jrrjr
Qr1rjrA1jr
2-5
由此式可见,频响函数在r处的值乘以模态阻尼因
r,就是留数〔的
估计值如图1。
利用这种模态检测方法之前,先要估计出
图1对频响应函数的幅值进行峰值和模态检测
3、圆拟合
圆拟合是一种单自由度方法,用频域中的模态模型对系统极点和复模
态〔或实模态〕向量进行局部估计。
此方法依据事实是:
单自由度系统的速度频响函数〔速度对力〕在奈奎斯特图〔即实部对虚部〕上呈现为一个圆。
如果把其他模态的影响近似为一个复常数,那么在共振频率r附近,频响函
数的根本公式为:
Htjj——-―———Rj12-6
jr
因此,首先要选择共振频率附近的一组频率响应点,通过这些点拟合
成一个圆。
阻尼固有频率r可以看成是复平面上数据点之间角度变化率最大〔角间隔最大〕的那个点的频率,也可以看成是相位角与圆心的相位角最为接近的那个数据点的频率。
对于分得幵的模态而言,二者的差异是很小。
阻尼比「估计如下:
2-7
式中!
,2:
分居在r两侧的两个频率点:
1,2:
分别为频率点在
r得半径之间的夹角
圆的直径和阻尼固有频率点的角位置含有复留数
H,tan
r
U+jV的信息:
2-8
式中:
圆的直径
:
园心与固有频率点的连线跟虚轴之间的夹角.
圆拟合法速度也很快,但为防止结果出错,特别是在模态节点附近,
需要操作者参与
〔二〕单自由度与多自由度系统
粘性阻尼单自由度SDO系统如图2的力平衡方程式表示惯性力、阻尼
力、弹性力与外力之间的平衡
MxtCxtKxtft2-9
其中M:
质量C:
阻尼K:
xxx:
加速度,速度,位移f:
外力t时间变量,把结构中所呈现出来的全部阻尼都近似为一般的粘性阻尼。
把上面的时间域方程变换到拉氏域复变量P,并假设初始位移和初始
速度为零,那么得到拉氏域方程:
Mp2CpKFp,或ZpXpFp
乙动刚度经过变换可得传递函数的定义,HpZ1p即XpHpFp
1/M
p2C/MpK/M
2-10
上式右端的分母叫做系统特征方程,它的根即是系统的极点是:
12C/2MVC/2M2K/M2-11
如果没有阻尼c=o,那么所论系统是保守系统。
我们定义系统的无阻尼固有频率为:
2-4
临界阻尼Cc的定义为使〔〕式中根式项等于零的阻尼值:
2-5
Cc
而临界阻尼分数或阻尼比Zi为:
Zi二CG,阻尼有时也有用品质因数即
Q因数表示:
Q1/2i2-6
系统按阻尼值的大小可以分成过阻尼系统〔Z1>1〕、临界阻尼系统〔Z
1=1〕和欠阻尼系统〔Z1<1〕。
过阻尼系统的响应只含有衰减成分、没有振荡趋势。
欠阻尼系统的响应时一种衰减振动,而临界阻尼系统那么是过阻尼系统与欠阻尼系统之间的一种分界。
实际系统的阻尼比很少有大于10%的,除非这些系统含有很强的阻尼机
制,因此我们只研究欠阻尼的情形。
11j1,11j1
2-7
在欠阻尼的情况下式2-11两个共轭复根:
其中1为阻尼因子1为阻尼固有频率。
有关系统极点的另外一些关系式
2-8
有:
1j.112
2-9
2-10
2-11
2-2式写成如下形式:
1/M
2-12
在展幵成局部分式形式,那么有:
-AU,这里A
P1
1/M
72~1
这里的A,和A是留数。
多自由度系统
多自由度系统可以用简单的力平衡代数方程演化成形式相似的一个
矩阵的方程。
下面是以而自由度系统为例。
如图:
―
品曲
K,
治
一J
"T7T7CT
图3多自由度系统
该系统的运动方程如下:
M1x1C1C2x1t
C2X2t
K1
K2x-!
t
K2x2tf
『1t
2-14
M2x2C2C3x2t
C2X1t
K2
K3X2t
K2x1t
f2t
写成矩阵形式是
M10为C1C2
C2
x
K1K
2K2
X1
f12-15
0M2x2C2
C2C3
x>
K2
K2K3
X2
f2
或者M
xCx
K
xf
2-16
其中[M、[q、[的、{f〔t〕}和{x〔t〕}分别为质量矩阵、阻尼矩阵、刚度
矩阵、方向量和响应向量。
把这个时间域的矩阵方程变换到拉氏域〔变量为p〕且假定初始位移和初始速度为零,那么得:
p2MpCKXpFp2-17
或者是ZpXpFp
式中:
[Z(p)]动刚度矩阵2-18
可以得到传递函数矩阵为:
HpZp12-19
式中adjZp:
Z叫的伴随矩阵,等于jZjT;
Zj:
Zp去掉第行第列后的行列式
1如果ij等于偶数
1如果ij等于奇数
传递函数矩阵含有幅值函数。
2-19式中的分母,即是Zp的韩烈士,叫做系统的特征方程。
与单
自由度情况一样,系统特征方程的根,即系统极点,决定系统的共振频率。
根据特征值问题,可以求出系统特征方恒的根。
为了把系统方程2-17转
化为
般的特征值问题公式,参加下面的恒等式:
pMpMX
0
2-20
将此式与2-17式结合在-
起得:
pABY
F'
2-21
其中
0
A
M
,B
M0
M
C
0K
pX
0
Y
,F
0
X
F
如果力函数等于零,那么式2-19就成了关于实值矩阵的一般特征值问题,其特征值马祖以下方程的p值:
pAB02-22
它的根就是特征方程Zp0的根。
对于N各自由度系统,此方程有
2N个呈复共轭对出现的特征根:
N
1
0
N
NJN
1J1
0
NJN
同单自由度系统一样,多自由度系统的极点的实部r是阻尼因子,虚
部r是阻尼固有频率。
〔三〕实模态和复模态
按照模态参数〔主要指模态频率及模态向量〕是实数还是复数,模态
可以分为实模态和复模态。
对于无阻尼或比例阻尼振动系统,其各点的振
动相位差为零或180度,其模态系数是实数,此时为实模态;对于非比例阻尼振动系统,各点除了振幅不同外相位差也不一定为零或180度,这样
模态系数就是复数,即形成复模态。
1复模态与实模态理论在拟合频段,实模态理论中传递函数在k点鼓励Z点响应的留数表达式为
Hki
rRklJr
■e
22
r12irVr
rarctan—r
Vr
k,l1,2,,n
(1)
rrJVr
其中,rRkl为留数;r和Vr构成的复数为系统的复特征值
拟合频段复模态理论中传递函数在k点鼓励f点响应的留数表达式为
n
Hki()
rRkleJr
r12jrj
rRkleJr
=eJr
Vrr12i%2
arctan
由〔1〕、〔2〕式中可以看出,传递函数共振峰处复模态的相位与实模态相
位的差异在于多出的复留数相位r,由传递函数的逆变换可以得到脉冲响
应函数,由此可以得到物理坐标系中结构的自由响应表达式
Xkr
krYrsinrtr
粘性比例阻尼:
t时刻第
r阶模态k点的振动为
Xkr
krYrertSindrt「
一般粘性阻尼:
t时刻第
r阶模态k点的振动为
对于无阻尼结构,t时刻第r阶模态k点的振动为
Xkr
2Trkremrtcosdrt
r
(3)
(4)
kr
式中,©kr表示振型幅值;Q表示模态频率;。
表示相位角
可以看出,无阻尼和比例阻尼系统的初相位与初始条件有关,与物
理坐标无关,具有模态〔振型〕保持性;而一般粘性阻尼系统的初相位
还与物理坐标k有关,每个物理坐标振动时并不同时到达平衡位置和最
不变,各点的振动周期、衰减率仍保持相同J.从物理坐标点的自由
响应公式还可看出,即使各测点留数为复数,但如果留数的相位差,即
振型的幅角相同,那么还是可以得到振动周期内形状不变且节点固定的
振型.这样模态虽是复模态,但表现出实模态的性质.因此实模态理论的实振型与复模态理论中复模态的差异在于各测点峰值相位差的大小.
2实模态提取方法
复模态理论中模态参数〔特征值和特征向量〕均为复数,在进行结构
模型修正时大量采用复数矩阵和复数迭代运算,计算工作量大,效率低;实模态理论中模态参数为实数,物理概念明确,后续结构模型修正计算公式简单,计算工作量小又节约空间,故实模态得到广泛的应用,实际测试得到的传递函数留数一般都为复数,要由复模态经过实模态提取技术才能得到实模态。
复模态提取实模态的方法主要有:
根据复模态的实部、虚部或相位确定实模态的传统方法;Ibrahim的扩大模型法;Chen的传递函数提取法等。
目前的模态分析软件中普遍使用的为传统方法。
由复模态实部或虚部获得实模态向量的方法为:
直接取复留数的实部或虚部作为实模态理论中的留数,进行规格化得到实模态振型.
由复模态相位获得实模态向量的方法为:
取复留数的幅值作为实模态理论中的留数,根据sinr的数值接近1或-1,将留数相位归为90?
°或-90?
°,然后尽享振型规格化,得到实模态振型,此振型中各测点相位差即为0?
°或180?
°。
用复模态理论获得的复模态向量,由复振型的周期变化中t=0即振动到达最大幅度时的振幅之比表示。
三、模态分析的应用与开展模态分析技术的应用可归结为以下几个方面:
1)评价现有结构系统的动态特性;
2)在新产品设计中进行结构动态特性的预估和优化设计;
3)诊断及预报结构系统的故障;
4)控制结构的辐射噪声;
5)识别结构系统的载荷。
对于实际的工程,用有限元软件分析需要的频率段,可查找振动原因,或校核。
模态分析可以看出在那些频率段需要防止或防止共振时很有用。
首先,频率和振型是结构的固有特性,任何结构都可以进行模态分析;其次,结构的功能是不同的,不同结构对应的模态分析的用途是有差异的。
对建筑结构,模态分析可以知道结构的避频设计、用于抗震设计计算以及考虑动力荷载的放大作用等。
另外,还可以挖掘振型有关的信息。
机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动千姿百态、瞬息变化。
模态分析提供了研究各种实际结构振动的一条有效途径。
首先,将结构物在静止状态下进行人为激振,通过测量激振力与胯动
响应并进行双通道快速傅里叶变换(FFT)分析,得到任意两点之间的机
械导纳函数(传递函数)。
用模态分析理论通过对试验导纳函数的曲线拟合,识别出结构物的模态参数,从而建立起结构物的模态模型。
根据模态叠加原理,在各种载荷时间历程的情况下,就可以预言结构物的实际振动的响应历程或响应谱。
模态分析软件以美国的MEScopeVES的功能最为全面。
MEScopeVES软件的功能包括信号处理(signalProcessing)、运行挠曲振型
(OperatingDeflectionShapes)、模态分析(ModalAnalysis)、结构改正(SDM)和声学分析(AcousticsAnalysiS)等,解决和分析机器与结构的
振动噪声问题。
主要可用于:
1、可以显示被测物体的实际工作形态(0DS)、模态、声学分布形态和
工程数据的形态等;
2、模块化结构便于用户根据自己的需要选择适宜的产品;
3、强大的图形显示、结构编辑、数据处理及动画显示功能;
4、软件的开放性好,能够与全球的多家厂商的硬件兼容;
5、主要应用的领域:
航空航天、建筑桥梁、汽车制造、钢铁冶金、军工装备等。
模态分析与参数辨识作为结构动力学中的一种逆问题分析方法并在工程实践中应用是从60年代中、后期开始,至今已有近四十年的历史了。
这一技术首先在航空、宇航及汽车工业中开始开展。
由于电子技术、信号处理技术与设备的开展,到80年代末这项技术已成为工程中解决结构动态性能分析、振动与噪声控制、故障诊断等问题的重要工具。
目前这一技术已渐趋成熟。
经过二十余年的研究开展,到目前为止模态分析技术已在我国各个工程领域中广泛应用,成为一种解决工程问题的重要手段。
在工程应用方面模态分析已渗透到我国各个工程领域,并取得了不少成就。
例如,某型火箭全装置的实物模态试验保证了火箭的准确发射与导航,防止了发射的失败;模态分析与参数识别技术曾被成功地用于解决某型航空发动机的严重振动故障,取得重大经济及社会效益;某型鱼雷全装置实物水下模态试验为鱼雷的振动与噪声控制确保导航性能提供了技术依据;远东第一高塔的上海东方明珠电视塔的振动模态试验,为高塔的抗风抗地震平安性设计提供了技术依据;目前世界上跨度第一的斜拉索杨浦大桥的振动试验对大桥抗风振动的平安性分析与故障诊断提供了技术依据;建立在模态分析技术上的桩基断裂检测技术已在高层建筑施工中广泛应用,提咼了桩基的质量,确保咼层建筑的平安;等等,这些成就不胜枚举。
总之,二十余年的开展是迅速的,成就是显着的,回忆这一开展过程和取得的成就,可更鼓励我们朝着新的目标发奋前进。
模态分析技术开展到今天已趋成熟,特别是线性模态理论方面的研究已日臻完善,但在工程应用方面还有不少工作可做。
首先是如何提高模态分析的精度,扩大应用范围。
增加模态分析的信息量是提高分析精度的关键,单靠增加传感器的测点数目很难实现,目前提出的一种激光扫描方法是大大增加测点数的有效方法,测点数目的增加随之而来的是增大数据采集与分析系统的容量及提高分析处理速度,在测试方法、数据采集与分析方面还有不少研究工作可做。
对复杂结构空间模态的测量分析、频响函数的耦合、高频模态检测、抗噪声干扰……等等方面的研究尚需进一步幵展。
模态分析当前的一个重要开展趋势是由线性向非线性问题方向开展。
非线性模态的概念早在1960年就由Rosenberg提出,虽有不少学者对非线性模态理论进行了研究,但由于非线性问题本身的复杂性及当时工程实践中的非线性问题并引引起重视,非线性模态分析的开展受到限制。
近年来在工程中的非线性问题日益突出,因此非线性模态分析亦日益受到人们的重视。
最近已逐步形成了所谓非线性模态动力学。
关于非线性模态的正交性、解耦性、稳定性、模态的分叉、渗透等问题是当前研究的重点。
在非线性建模理论与参数辨识方面的研究工作亦是当今研究的热点。
非线性系统物理参数的识别、载荷识别方面的研究亦已开始。
展望未来,模态分析与试验技术仍将以新的速度,新的内容向前发展。
姓名:
徐海滨
学号: