振动及冲击JOURNALOFVIBRATIONANDSHOCK.docVIP

振动与冲击 JOURNAL OF VIBRATION AND SHOCK 运用小波分析方法进行结构模态参数识别 朱宏平，翁顺 (1.华中科技大学土木工程与力学学院，武汉430074) 摘要结构的模态参数反映了结构自身特性，是基于动态特性的结构损伤识别和健康评估的重要因子。本文首 先介绍了环境激励下基于小波分析的模态参数识别方法，针对土木工程结构的前几阶自振频率处于低频区域以及环境激 励下结构响应信号信噪比很低的特点，着重论述了采用小波方法抑制原始测量信号中的高频成分（即噪音），从而突出结 构低频特性的降噪处理方法的基本原理。通过比较传统傅里叶变换、短时傅里叶变换和小波变换三种方法对一实际高层 建筑结构现场测试信号的处理结果以及有限元分析结果，认为小波分析方法可以更精确、更有效地识别工程结构的模态 参数。 关键词：傅里叶变换，短时傅里叶变换，小波变换，降噪，模态参数 中图分类号：TN911.6; TU311.4 文献标识码：A 在高层建筑抗震、抗风、健康监测及损伤诊断等研 究中，结构模态参数是非常重要的设计参数之一，基于 环境激励的模态参数识别方法越来越受到人们的重 视[1]。目前国内外在结构模态参数识别方面的研究方 法有很多，主要可以分为：①频域方法：它是建立在频 响函数的理论基础上的，频域法的最大优点是利用频 域平均技术，最大限度地抑制了噪声影响，使模态定阶 问题容易解决，但也存在着如功率泄露、频率混叠、离 线分析等问题)②时域方法:是直接利用响应的时域信 号进行模态参数识别。与频域法相比，时域法对于分 离密集模态有更好的效果;③小波分析法[2(:比短时傅 里叶变换具有更好的时频窗口特性，克服了傅里叶变 换中时-频分辨率恒定的弱点，因此它能在具有足够 时间分辨率的前提下分析信号中的短时高频成分，又 能在很好的频率分辨率下估计信号中的低频。④基于 H HT变换的非平稳信号的处理方法：它以瞬时频率 为基本量，以固有模态信号为基本信号，用于非平稳 信号处理;⑤基于模拟进化的模态参数识别的方法： 该方法实现了基于达尔文进化理论的整体优化算法 用于识别线性振动结构的模态参数。基于模拟进化 的模态参数识别方法用于测试噪音是很可靠的，但该 方法用于识别更复杂模态的现实问题上，还需要更进 一步研究。 本文针对实际工程，对环境激励下的高层建筑振 动响应信号采用小波方法进行分析[3 (，有效地识别高 层建筑固有模态参数，并同传统傅里叶变换、短时傅里 叶变换的结果相比较，证实了小波分析方法在处理随 基金项目：国家自然科学基金和教育部新世纪优秀人才基金 (2004年）资助 收稿日期：2006 -06 -21修改稿收到日期:2006 -07 -27 第一作者朱宏平博士，教授，博士生导师,1965年11月生 机信号方面的优越性。由于环境激励下的振动测试信 号信噪比低[4(，对信号分析造成一定的干扰，本文用小 波分析方法对测试信号进行降噪处理，结果表明小波 方法能有效抑制噪音，还原真实信号，提取更多有用 信息。 1小波分析方法基本理论 傅里叶变换的实质是把波形分解成许多不同频率 正弦波的叠加，是傅里叶级数在连续情况下的推广，函 数!) L ($)的傅里叶变换为： %(!) = fV!/( ( D ⑴ J — A 短时傅里叶变换在傅里叶分析基础上引入时域信 息的最初尝试，它的基本思想是:把信号划分成许多小 的时间间隔，用傅里叶变换分析每一个时间间隔，以便 确定该时间间隔存在的频率。以高斯函数ga (/)=从频域上看，用不同尺度作小波变换相当于用一 组带通滤波器对信号进行带通滤波处理。在表达式 (6 )中.求得信号在固定小波函数l#,b( t)上的分量.对 参数a和b进行展开后，就得到任何时刻，任意精度的 频谱了。但是对实际计算来讲,这样的代价太高,类似 傅里叶级数的想法，把参数a或6或（a, b)同时做离 散,得到离散小波变换。离散小波函数可表示为： #3 = a#( a - 40 )，j,4 6 (7 ) 离散小波变换的系数可表示为： 1,4 ( t) = \_Af( t) #3,4 ( t) dt (8 ) 设函数##( t) L2 ($)，如果存在两个常数A.B. 且04B^使得稳定条件： A % 36 #(2!) I2 % B (9) 几乎处处立，如果A=B则称为最稳定条件。那 么函数序列丨12,(4)丨叫做函数/的二进小波变换： 1/(4) = [/(t)，#2,(4)] 广 =去f/(t)#$ (2-t - 4)dt (10) J-A 由小波的定义可知，小波函数和其傅里叶变换 #o,b( !)满足窗函数的条件，小波函数的时间窗和频率 窗可以表示为： #a，b(t) ： [ b + at 0 - a$#，b +