Lamb(兰姆)波是二维波,与三维体波相比具有衰减速度慢,传播距离远的特点,因此常被用于大型板材的长距离及快速无损检测中。板材中兰姆波与管中、变截面波导介质中的导波一样,具有频散性与多模态性。加上环境噪声等多方面因素的影响,导波检测时传感器接收到的Lamb波信号非常复杂,属于非平稳随机信号,需要利用有效的信号处理技术提取有用的信息成分才能确定合适的激励方式,获得更好的检测成像效果。传统的处理Lamb波信号的方法包括反射系数法、傅里叶变换法、小波变换法、动态光弹法等,但是这些方法都有各自的不足。反射系数法是通过测量漏兰姆波的频散曲线来确定材料的性质,但测量难度较大。傅里叶变换只能处理线性非平稳的信号。小波变换法虽然在理论上能处理非线性非平稳信号,但是同傅里叶变换、短时傅里叶变换法一样,都受Heisenberg测不准原理制约,即时间窗口与频率窗口的乘积为一个常数,这就意味着如果要提高时间精度就得牺牲频率精度,反之亦然。当兰姆波中不同模态的频率比较接近时,不适用小波变换处理信号。动态光弹法能从Lamb波的应力分布观察到传播和频散,但是在实际检测中对硬件要求较高。
HilbertGHuang变换(HHT)是一种近几年发展起来的一种自适应信号处理方法,不受Heisenberg测不准原理制约,可以在时间和频率上同时达到很高的精度,非常适用于分析突变信号。笔者以薄壁铝板为研究对象,利用双重时间尺度的方法,即采用二维傅里叶变换法整体传播时间尺度,HilbertGHuang变换从单一信号时间尺度,将二者相结合对在铝板中不同位置采集到的Lamb波信号作数据处理与分析,与半解析有限元法得到兰姆波的频散曲线相对照,进而识别与分析铝板中兰姆波模态,获得较高的时间分辨率。
薄板中超声导波分离模态信号分析方法
一、传播特性的半解析有限元分析
在边界自由的固体板中,板材厚度与激励声波波长数量级相当时,在板中产生的应力波就是Lamb波.根据薄板两表面质点振动相位关系,Lamb波分为对称型Lamb波和反对称型Lamb波,分别用Sn、An 表示,n 为模态的阶次(n=0,1,2,3)。对称模态和反对称模态示意如图1所示,对称模态中沿着板厚方向中心点对称的各节点在厚度方向的位移方向相反;而反对称模态中,沿着板厚方向中心点对称的各节点在厚度方向的位移方向相同。
以铝板中Lamb波的传播为例,其频散曲线可采用半解析有限元法求得,只需要在波导介质的截面上作有限元离散,而沿波导介质传播方向的位移则以简谐波的振动方式表示,在对介质截面进行有限元离散后,根据哈密顿原理可以推导出导波在介质中的波动方程,求解特征值可以得到波数和频率的关系,进而绘制出频散曲线。
通过半解析有限元法,以铝板厚度d=0.8mm绘制铝板中导波的相速度和群速度频散曲线,得到的模态分布分别如图2,3所示。
二、二维傅里叶变换
Lamb波在时间和空间上都可以通过二维傅里叶变换转换为二维各个离散频率点的频率G波数能量谱,从而分解出单个Lamb波,并可对其幅值进行测量。单个波动组分在时间上的频度称为频率,而在空间(距离)上的频度称为波数.由频率波数谱中某个波动组分的频率和波数,可以确定周期和波长。通过对接收信号的二维傅里叶变换,与理论计算得到的波数G频率的频散曲线进行对比,从而确定检测信号中包含的Lamb波模态。
三、HilbertGHuang变换