模态分析报告.docx
《模态分析报告.docx》由会员分享,可在线阅读,更多相关《模态分析报告.docx(13页珍藏版)》请在冰豆网上搜索。
模态分析报告
汽车挡风玻璃模态试验报告
(模态分析理论与试验第三小组)
姓名:
张朋
小组:
三组
学号:
1133036
1试验目的
1.学习模态分析原理和模态测试方法;
2.试验分析得到汽车挡风玻璃的前15阶模态的模态参数;
3.试验分析汽车挡风玻璃的动态振动特性;
4.为汽车挡风玻璃的有限元分析计算模型的修改提供可靠依据。
2试验仪器
试验仪器如表1所示:
表1试验仪器列表
编号
仪器名称
数量
仪器型号
仪器编号
厂家
1
力锤
1
LC-02系列
75704
江苏联能
2
加速度传感器
10
CA-YD107压电式
江苏联能
3
电荷放大器
2
YE5853A
江苏联能
3
数据采集系统
1
INV306DASP
DASP
4
数据分析处理软件
1
CoinvDASP2005
DASP
3试验对象
试验对象:
POLO三厢车前挡风玻璃。
实验对象附件描述见表2
表2实验对象描述
序号
名称
有
无
1
密封胶条条
√
2
后视镜安装塑料座
√
3
粘接胶水残留
√
4试验测量和分析系统
试验测量分析系统由三大部分组成:
试验试验激振系统,响应采集系统,模态分析和处理系统。
其中,
(1)试验激振系统包括:
江苏联能LC系列力锤;
(2)响应采集系统包括加速度传感器、和DASP信号采集系统;
(3)模态分析和处理系统主要是DASP和Matlab软件。
具体的组成方式如图1和图2所示。
图1模态试验测量分析系统模型示意图
图2模态试验测量分析系统—电荷放大器
图3模态试验测量分析系统—INV306UDASP数采系统
5实验原理
5.1传递函数
试验模态分析是基于系统响应和激振力的动态测试,即通过振动测试,经信号处理和参数识别确定系统的模态参数,建立以模态参数表示的运动方程。
从模态分析理论可知,这些参数可以通过传递函数或频响函数曲线进行分析求得。
在时域内,系统的极点、固有频率和阻尼比由系统的特征方程决定,即由传递函数的分母决定;在频域内,以输入-输出信号在频域中的关系,可以采用模态参数识别法求得动态参数。
试验模态分析必须借助于模态分析理论的指导,才能够正确的求出系统的动态参数。
系统的传递函数是输出的拉氏变换与输入的拉氏变换之比,用频率f代替s即可得到频响函数
但由于单次数据测取并不可靠,常采取平均技术来减少干扰和噪声带来的影响。
所以有
从上式可以看出,频响函数与外界输入无关,它是系统的固有属性,可以避免输入含有噪声而影响输出结果分析,从而影响模态参数的正确性。
为了判断输入输出线性关系的可靠性,运用相干函数来进行检验,即相干函数是判断试验可靠性的有力工具。
5.2相干函数
一般情况,不知道系统的响应信号x(t)中有多少成分是由激励信号f(t)产生的,这时用相关函数和传递函数的关系来衡量这一问题,可以得出一种有效的方法,这即是相干函数法。
对于单输入单输出,可以用下图简单描述:
为了反映第i个输入所提供的分量在第j个输出总量中所占的权重以及检查系统有否噪声干扰和非线性的程度,引入相干函数通常用球(f)表示,其定义:
如果相干函数为零,表示输出信号与输入信号完全不相干,所得传递函数是无效的。
那么,当相干函数为1时,表示输出信号与输入信号完全相干。
若相干函数在(0一l)之间,则表明有如下三种可能:
(l)测试中有外界噪声干扰;
(2)输出x(t)是输入f(t)和其它输入的综合输出;
(3)联系x(t)和f(t)的线性系统是非线性的。
相干函数的一个重要意义是检验系统非线性程度。
将F(t)作为系统输入的随机过程,X(t)作为系统输出的随机过程,由系统是在假设的线性条件下得到的:
所以,谱相关函数值的大小可以用来检查系统的非线性程度。
5.3误差控制
本试验采用提高激振频率的方法减小迭混误差。
采用加窗与猝发激励减小泄漏误差。
本试验使用的冲击锤是产生脉冲信号最常用的激励装置,理想脉冲信号即δ(t)函数,其傅氏谱为一水平直线,频率成分宽。
现实中的脉冲信号是一有限宽度和有限高度的信号,其时间历程和自谱见图(a),(b)。
从图(a),(b)可以看出,在低频能量近似均匀分布,而在高频段能量逐渐衰减。
脉冲信号的频率成分和能量大致可控,试验周期短,无泄漏,但是信噪比差,特别是对大型机构,冲击锤产生的激励能量往往不足以激起足够大的响应信号,所以冲击信号的高频响应较差。
图(a)冲击锤产生的激振信号时间历程
图(b)冲击锤产生激振信号自谱图
6模态分析方法和测试过程
6.1激励方法
从频响函数的物理意义可知,若知道激励和响应,就可推知系统的特性。
从这个意义上来说,有两种激励方法可供选择,其一是对结构上某点激励,测得所有点响应,即单点激励的方法。
其二是对结构某些点同时激励,测得各点的响应,即通常所说的多点激励方法。
单点激励的优点是在模态试验中,只要同时测量记录激励和响应的信号,再经数字信号处理,可获得与响应激励自由度对应的频响函数,简单易行。
单点激励的缺点是能量有限,在大型复杂结构的模态试验中,响应信号信噪比差,有时难于激励出所关心的整体模态,当试验结构的模态比较密集且阻尼又比较大时,难于激出结构的纯模态。
多点激励的优点是可以激发出试验结构的纯模态,对大型结构尤其有效。
而且多点激振可以解决系统存在重根的问题,但是多点激励也存在不足,其缺点是需要采用几套激振器激励。
模态试验中,各激振器的激励信号应互不相关,输入力的大小和相位要反复调试,比较费时,加上设备昂贵,目前,尚未普遍推广。
在进行模态试验时,试件采用单点激励还是多点激励方法取决于试件被整体激振的难度。
如果单点激励就可以测得试件上任意点的响应,且响应幅度足够大,则采用单点激振即可,否则需要对试件进行多点激振。
由于本试验研究的是挡风玻璃的模态参数,其刚度比较小,总体比较容易被激振。
而且由于低频段是我们所感兴趣的频带,不考虑重根的存在,所以这里采用单点激振。
6 .2结构安装方式
在测试中使结构系统处于什么状态,是试验准备工作的一个重要方面。
一种经常采用的自由状态。
即使试验对象在任一坐标上都不与地面相连接,自由地悬浮在空中。
如放在很软的泡沫塑料上;或用很长的柔索将结构吊起而在水平方向激振,可认为在水平方面处于自由状态。
另一种是地面支承状态,结构上有一点或若干点与地面固结。
如果在我们所关心的实际情况支承条件下的模态,这时,可在实际支承条件下进行试验。
但最好还是自由支承为佳。
因为自由状态具有更多的自由度。
本试验采用轮胎简支的方式固定试件,如下图所示。
图(c)试件支撑方式
7、实验步骤
7.1测点的确定
对于需要得到模态振型的测试件,测点选取原则如下:
1.在基本反应测试件轮廓的前提下,测点尽量少;2.避开各阶振型节点;3.能明确显示模态振型特征。
对于仅需得到模态频率的测试件,测点选取原则如下:
1.在保证数据采集全面的前提下,测点尽量少;2.避开各阶振型节点。
试验时,布置了28个测点,对这些测点只进行Z向的振动加速度信号的采集。
分别完成测点模型的建立(见图5),以及实物机构的测点布置(图2)。
图5测点布置示意图
表3测点坐标
7.2仪器连接
仪器连接如图1所示,其中力锤上的力传感器和压电加速度传感器接电荷放大器的电荷输入端。
经过放大器的低通滤波,输出通道分别接到INV306U接口箱(盒)上的第输入通道。
7.3结构生成及约束
仪器安装好后,启动DASP2005专业版,选择模态分析按钮,进入模态分析界面。
输入测点号及其坐标,生成结构的节点及线,面等几何元素。
7.4参数设置与采样
在DASP信息窗口中选择随机采样,进入采样界面。
在测量设置中设置传感器类型、总测点数和原点导纳位置,总测点数根据结构自动读取,不可更改。
用力锤击激振点测点,观察有无波形,如果有一个或两个通道无波形或波形不正常,就要检查仪器是否连接正确、导线是否接通、传感器、仪器的工作是否正常等等,直至示波波形正确为止。
在采样参数设置中选定采样频率1200Hz,采样长度2K,程控倍数为1,使用适当的敲击力敲击各测点,调节放大器的放大倍数和INV的程控倍数,直到力的波形和响应的波形即不过载也不过小。
选定采样时自动增加测点号,准备采样。
采样类型设为:
变时基;单位类型设为:
第1-9通道工程单位设为:
m/ss(加速度),,第10通道工程单位设为:
N(牛顿);其余通道本试验不用。
输入标定值和工程单位。
通过示波已经定好了放大器的档位,传感器灵敏度为KCH(PC/U)表示每个工程单位输出多少PC的电荷,如是力,而且参数表中工程单位设为牛顿N,则此处为PC/N;如是加速度,而且参数表中工程单位设为m/s2,则此处为PC/m/s2;试验采用的是江苏联能生产的型号为YE5853A的电荷放大器,在试验前必须先用示波器调节波形的幅值,以确定放大器的增益及灵敏度的设定,为采样程序参数设置表中输入各通道的工程单位和标定值。
在试验时,我们把与力锤相连的电荷放大器的低通拨到10K挡,高通拨到加速度传感器挡,根据计算其标定值为100mv/EU;与加速度传感器相连的电荷放大器的低通拨到也拨为I0K挡,高通拨到加速度传感器挡,根据计算其标定值为I0mv/EU;按下式设定DASP采样参数表中的标定值:
参数设置完以后,选择自动增加测点号,按左窗下面的开始采样按钮,进入触发变时基采样状态,等待触发,并提示当前采样的点号和触发次数。
根据提示从第一点设定的触发次数测试到最后一个测点。
自动记录下每次测试结果。
测试过程中尽量避免连击现象,如果有连击现象,按中止采样按钮,改变测点号重新开始采样,将覆盖原来数据。
7.5实验数据分析处理
7.5.1调采样数据
采样完成后,进入DASP信号分析功能模块,对采样数据进行传函分析。
首先选择要调入的测点号,按调入波形按钮,右面窗口中显示该测点的波形。
给每一通道的力信号加力窗,对响应信号加指数窗,选择系数1。
选择通用传函分析,具体设置如图6所示。
图6传递函数设置图7信号加窗
设置完成后,按进行传递函数计算,完成测点的传函分析,显示分析结果。
按自动计算全部传函按钮,可以分析完全部采样点的传函分析。
在实验进行过程中,检查信号的相干性,只有相干系数(除去节点或反节点外)在0.8以上的信号才有效,如图6所示,可见本次试验的信号在研究频段内相干系数在研究频率范围内均大于0.8,信号符合要求。
图8测点的传递函数以及相干系数
7.5.2模态定阶与拟合
采用集总平均的方法进行模态定阶,按开始模态定阶,显示集总平均后的结果,用鼠标分别点峰值点,收取该阶频率,依次收取各阶频率,按保存按钮存盘。
如果收取有误可按清除按钮清除当前结果。
(1)集总平均:
将所有频响函数进行平均,得到一条平均后的曲线,由此来确定模态阶数,其好处是不容易丢失某阶模态。
(2)一点频响函数:
选所有频响函数中的一条来确定模态阶数,往往和多自由度模态拟合方法配合使用,所选的一条频响函数必须具有代表性。
(3)集总显示:
同时显示所有频响函数曲线,来确定模态阶数。
其好处是不容易丢失某阶模态,但当频响函数曲线条数很多时,画面有时较凌乱。
本试验采用集总平均的方法确定试件各阶模态频率。
如图9所示
图9集总平均法模态定阶
模态拟合采用复模态单自由度拟合方法,按开始模态拟合得到拟合,得到拟合结果。
如图10所示。
图10拟合结果
3振型编辑
质量归一和振型归一两种方式随各自需要任选,本试验选择质量归一,模态振型编辑完毕后,进行模态振型相关矩阵校验,如图11所示。
振型相关矩阵校验可用来校验各阶模态振型之间的正交性,矩阵关于主对角线对称,主对角线的元素都为1。
矩阵元素的行号和列号分别代表了两阶模态,其大小表示了这两阶模态振型的正交性,为归一化后的两阶模态振型标量乘积,值越小表示正交性越好。
理想的模态分析,其结果的振型相关矩阵除主对角线元素外,其它元素的值都很小。
图11为试验的振型相关矩阵校验图,从中可以看出,虽然存在一些误差,但仍符合主振型的特点,说明试验是正确可行的。
图11模态振型相关矩阵校验
至此,模态分析完毕,以后可以观察、打印和保存分析结果,也可以观察模态振型的动画显示。
4动画显示
按进入动画显示对话框,根据每个对话框中的相应按钮可以动画进行控制,如更换显示阶数、显示轨迹;在视图选择中选取显示方式:
单视图、多模态和三视图;改变显示色彩方式;振幅、速度和大小,以及几何位置。
对于当前动画可输出AVI格式的动画文档,可直接用媒体播放器播放。
按模态输出为AVI文件,弹出保存文件对话框,模态动画视频压缩对话中的压缩程序选择MicrosoftVideo1.0方式,确定后即可生成动画文件。
在保存的目录下调用文件,显示动画。
8、实验结果和分析
8.1模态频率和阻尼
通过以上步骤,得到挡风玻璃的前13阶模态,其模态参数如表4所示。
表4:
模态频率和阻尼
阶数
频率(Hz)
阻尼(%)
模态质量
模态刚度
模态阻尼
1
9.701
8.661
1.0000e+000
3.7155e+003
1.0559e+001
2
19.962
4.584
1.0000e+000
1.5731e+004
1.1499e+001
3
22.481
3.712
1.0000e+000
1.9952e+004
1.0486e+001
4
40.400
1.961
1.0000e+000
6.4435e+004
9.9562e+000
5
46.576
1.113
1.0000e+000
8.5641e+004
6.5133e+000
6
67.950
0.729
1.0000e+000
1.8228e+005
6.2211e+000
7
96.827
0.850
1.0000e+000
3.7012e+005
1.0339e+001
8
100.904
0.522
1.0000e+000
4.0196e+005
6.6174e+000
9
138.252
1.624
1.0000e+000
7.5458e+005
2.8209e+001
10
145.658
1.271
1.0000e+000
8.3759e+005
2.3255e+001
11
164.649
0.943
1.0000e+000
1.0702e+006
1.9504e+001
12
206.187
1.002
1.0000e+000
1.6783e+006
2.5960e+001
13
226.917
1.064
1.0000e+000
2.0328e+006
3.0344e+001
2、打印出各阶模态振型投影图
图12一阶模态振型(9.701Hz)
图13二阶模态振型(19.962Hz)
图14三阶模态振型(22.481Hz)
图15四阶模态振型(40.4Hz)
图16五阶模态振型(46.576Hz)
图17六阶模态振型(67.950Hz)
图18七阶模态振型(96.827Hz)
图19八阶模态振型(100.904Hz)
图20九阶模态振型(138.252Hz)
图21十阶模态振型(145.658Hz)
图22十一阶模态振型(164.649Hz)
图23十二阶模态振型(206.187Hz)
图24十三阶模态振型(226.917Hz)
表5挡风玻璃Z向前10阶振动模态阵型描述
阶数
模态频率(Hz)
模态振型描述
1
9.701
一阶弯
2
19.962
一阶扭转
3
22.481
局部模态
4
40.400
弯扭组合
5
46.576
二阶弯曲
6
67.950
三阶弯曲
7
96.827
四阶弯曲
8
100.904
弯扭组合
9
138.252
五阶弯曲
10
145.658
二阶扭转
注:
振阵型描述引人而异,具体以实际振动动画为准
8.2试验与仿真对比
表6试验与仿真所得固有频率对比
阶数
试验值(Hz)
仿真值(Hz)
误差(%)
1
9.701
7.38
-23.9254
2
19.962
18.269
-8.48111
3
22.481
21.169
-5.83604
4
40.4
34.651
-14.2302
5
46.576
47.608
2.215733
6
67.95
67.686
-0.38852
7
96.827
90.311
-6.72953
8
100.904
104.98
4.039483
9
138.252
133.19
-3.66143
10
145.658
142.41
-2.22988
11
164.649
165.19
0.328578
12
206.187
208.97
1.349746
表7仿真各阶模态振型(前10阶)
阶数
频率(Hz)
振型
1
7.38
2
18.269
3
21.169
4
34.651
5
47.608
6
67.686
7
90.311
8
104.98
9
133.19
10
142.41
8.3分析结论
含集中质量的弹性板有限元计算所得的固有频率与试验值的相对误差绝对值大体上小于6%。
两种分析结果存在误差的原因是多方面的,主要有:
(1)有限元建模时,模型结构的简化,不可避免要产生一些误差;结构本身可能存在非线性因素,而在计算时,假定其为线性系统进行线性求解,必然造成一些误差;
(2)有限元方法本身也产生离散误差。
(3)试验条件所限,试件的边界条件和受力产生时变。
总体上,有限元计算结果与试验值误差较小,表明有限元模型能较好地反映结构的动态特性,模型简化较为合理。