一种基于测量波数曲线的Lamb波频散补偿方法技术

本发明专利技术提供了一种基于测量波数曲线的Lamb波频散补偿方法，其包括下列步骤：（1）测量波数曲线；（2）对测得波数曲线进行线性化；（3）求取Lamb波脉冲响应频谱；（4）求取线性化波数下的脉冲响应；（5）构建频散补偿的Lamb波信号。本发明专利技术具有如下技术效果：（1）借助于快速傅立叶和反傅立叶变换算法，便于快速实现；（2）可再压缩信号中每个频散波包，提高了信号的分辨率和信噪比；（3）维持信号中每个波包的时域位置不变，便于后续损伤定位或成像；（4）直接在结构中测出波数曲线，适用于结构材料参数未知无法求取理论波数曲线的场合。

【技术实现步骤摘要】

本专利技术涉及一种Lamb波信号处理方法，尤其涉及一种Lamb波频散补偿方法。
技术介绍
诸如载人飞行器、核反应堆和桥梁等重要结构在其服役过程中难免产生各种形式的损伤，为了避免结构损伤带来的灾难或损失，必须对这些结构进行长期有效的监测，而传统的无损检测技术已不太能满足该监测要求。为此，研究者们提出了结构健康监测的概念。结构健康监测技术是一种在线、动态、实时的监测技术，其近几年得到高速发展并在工程结构的安全和可靠性评估中发挥着日益重要的作用。Lamb波作为一种板类结构中传播的超声导波，因其能够进行长距离传播并且对结构表面和内部损伤均敏感，在结构健康监测领域得到日益广泛的关注，Lamb波监测技术已成为一个研究热点。Lamb波具有多模特性，在实际应用中通常采用合适中心频率的窄带激励信号向被测结构中选择性地激发单一的基本对称(Stl)或反对称(Atl)模式的Lamb波信号。但由于Lamb波频散特性的存在，使得监测信号中的波包在传播过程中发生扩展，波包形状发生畸变，幅值也随之降低，影响了信号的信噪比和分辨率，为后续信号分析和损伤识别增加了难度。现有技术中，时间反转方法无需Lamb波在结构中传播的先验知识便可自动补偿频散特性，但同时也消除了 Lamb波的传播时间，为后续损伤识别增加了难度。Alleyne通过信号再次激励来抑制传播距离已知的特定模态的频散现象。Sicard等基于后向传播函数提出一种频散补偿方法，该方法需要积分运算，计算量较大。Wilcox通过把信号从时域变换到空间域来消除频散特性，由于信号域的变换，处理过程较为复杂。Liu和Yuan对传感信号频谱进行插值处理来补偿铝板中频散的A0模式信号，该方法没有考虑插值处理对原始激励信号的影响，而且需要已知结构材料参数来计算理论波数曲线，应用场合受到限制。
技术实现思路
本专利技术的目的是克服现有技术的不足之处，提供。本方法无需知道结构材料参数，而是根据传感信号直接测量出波数曲线，然后基于该曲线和结构脉冲响应构建出非频散的传感信号，实现频散补偿。本专利技术的基于测量波数曲线的Lamb波频散补偿方法，包括下列步骤(1)测量波数曲线选择具有合适中心频率的窄带激励信号向结构中激发Lamb波信号并采集相应的传感信号，从传感信号中提取出所选Lamb波模式的直达波，分别求取激励信号和直达波的相位变化曲线，则波数曲线可计算为Κ(ω) = /1其中Φ1(ω)和Φ2(ω)分别为激励信号和直达波的相位曲线，L为激励和传感器之间的距离；(2)对波数曲线进行线性化把κ(ω)在激励信号的中心频率ω。处进行泰勒级数展开并只保留其中第0，1阶次项，得到线性化的波数曲线Κ1 η(ω) = Κ(ωε)+1/ο8(ωε)δ(ω-ωε)其中cg( 。)为ω。处的群速度；(3)求取Lamb波脉冲响应频谱在结构中的激励上加载阶跃激励信号，利用传感器采集结构的阶跃响应，对该响应进行求导运算可获得对应于该激励-传感器对的结构脉冲响应h (t)，最后对脉冲响应进行傅立叶变换得到Lamb波脉冲响应频谱Η(ω)；(4)求取线性化波数下的脉冲响应根据线性化前后的波数，可把每个频率值ω改变为ω1 η(ω) = Π — )]其中IT1M为测量波数曲线Κ(ω)的逆函数，对Η(ω)在ω = ω1 η(ω)处进行插值处理得到新的传感信号频谱Hlin( )，对Η1 η(ω)进行逆傅立叶变换可得到线性化波数下的脉冲响应hlin(t)；(5)构建频散补偿的Lamb波信号对窄带激励信号a(t)和线性化波数下的脉冲响应hlin(t)进行卷积运算便能得到频散补偿后的传感信号^in (t)。优选地，所说的步骤(1)中测得的波数曲线可以为绝对波数曲线，也可以为相对波数曲线，即只要求测得波数曲线与实际波数曲线具有相同的变化趋势，以便于该方法在工程实际中的应用。本专利技术的基本原理是从信号与系统的角度可把Lamb波单模式传感信号在频域中简单表示为激励信号频谱与相位延迟因子的乘积，即ν(ω) = Va( )不同频散情况下的相位延迟因子不同，可把它写成复函数形式，E (ω) = E 。母函数E(K)中K E之间为指数函数关系，与频散特性无关，而是各频散特性下的子函数Κ(ω)不同。非线性频散情况下，Κ(ω)中ω K之间的非线性变化关系使激励信号中不同频率成份Ι(ω)在频域中的相位延迟不同，所以对应于时域中的波达时间也不一致，造成传感信号ν(ω)中波包在时域中扩展，幅值降低，波包形状也发生畸变，表现出频散效应。如果对非线性的波数曲线Κ(ω)线性化，得到Klin(Co)。新的传感信号可表示为= 。由于Klin(Co)中线性的ω K变化关系，使不同频率下的Va(co)具有相同的相位延迟，在时域中会同时延迟到某一位置发生同相叠加，保持了波包形状和时域宽度。如果在ν(ω)的基础上构建出Vlin(co)，就能实现频散补偿。信号构建中最关键的是如何求取相位延迟因子=由于不同频散特性下的相位延迟因子Ε(ω)具有相同的母函数E (k)，可以对E (ω)进行插值得到Elinl (ω),从而最终得到Vlin (ω)。由于本专利技术所公开的方法只对脉冲响应频谱进行插值，能够避免Liu和Yuan所提方法中因直接对传感信号频谱进行插值改变原始激励信号频谱的问题。而且，本专利技术中只是根据原始波数曲线Κ(ω)和线性化波数曲线Κ1 η(ω)之间的相对关系来求取ω1 η(ω)的，所以只需测出的相对波数曲线与实际波数曲线具有相同的变化趋势即可，提高了所提频散补偿方法的适用性。本专利技术的基于测量波数曲线的Lamb波频散补偿方法，具有如下技术效果(1)借助于快速傅立叶和反傅立叶算法，便于快速实现；(2)可再压缩信号中每个频散波包，提高了信号的分辨率；(3)维持信号中每个波包的时域位置不变，便于后续损伤定位或成像；(4)直接在结构中测出波数曲线，适用于结构材料参数未知无法求取理论波数曲线的场合。附图说明图1是基于测量波数曲线的Lamb波频散消除方法的流程图图2是压电片和损伤的分布情况；图3a是激励信号；图3b是提取的A0模式直达波；图3c是得到的相位差曲线；图3d是测得的波数曲线；图3e是线性化后的波数曲线；图4 是 P3--P5中频散的Atl模式传感信号；图5 是 P3--P5中的脉冲响应时域波形；图6 是 P3--P5中线性化波数下的脉冲响应时域波形；图7 是 P3--P5中经频散补偿的传感信号。具体实施例方式实施例1本实施例采用玻璃纤维环氧树脂复合材料板作为待测结构，尺寸为IOOOmmX IOOOmmX 2mm0该结构材料参数未知，无法求出理论波数曲线。在结构中布置有8个压电片P1 P8组成的压电阵列，以结构中心为坐标原点建立直角坐标系，则压电片在坐标系中的分布情况及其位置分别见图2和表1。监测设备由波形发生与数据采集系统、功率放大器和矩阵开关控制器组成。其中，波形发生与数据采集系统中的LAI200-ISA任意波形发生卡、信号放大卡和PCI-9812数据采集卡分别实现Lamb波激励波形的产生、传感信号的放大以及信号的采集功能，功率放大器增强激励信号以扩大Lamb波在结构中的监测范围，矩阵开关控制器则根据事先确定的扫查策略选择相应的监测路径。表1压电片的坐标(单位mm)权利要求1.，其特征在于，包括下本文档来自技高网...

【技术保护点】

【技术特征摘要】

【专利技术属性】
技术研发人员：石立华，蔡建，周颖慧，李炎新，付尚琛，张兵，马丁，
申请(专利权)人：中国人民解放军理工大学，
类型：发明
国别省市：