一种具有频偏修正效果的时间-距离域映射方法技术

技术编号:8078187 阅读:218 留言:0更新日期:2012-12-13 20:13
本发明专利技术提供了一种具有频偏修正效果的时间-距离域映射方法,包括下列步骤:1)求取频散波数和非频散波数关系;2)计算原始频散信号频谱:3)求取信号波数谱;4)修正激励信号波数谱;5)计算得到频偏修正的非频散距离域信号。本发明专利技术的方法可补偿Lamb波信号中的频散效应,再压缩信号中因频散发生扩展和变形的各个波包,提高信号分辨率和信噪比。而且,在处理过程中对激励信号波数谱进行了修正,从而保持频散补偿前后信号波数成份不变,便于信号解释。

【技术实现步骤摘要】

本专利技术属于Lamb波信号处理领域,尤其涉及Lamb波频散补偿方法。
技术介绍
诸如载人飞行器、核反应堆和桥梁等重要结构在其服役过程中很容易产生各种形式的损伤,为了避免结构损伤带来的灾难或损失,必须对这些结构进行长期有效的监测,而传统的无损检测技术已不能很好满足该监测要求。为此,研究者们提出了结构健康监测的概念。结构健康监测技术是一种在线、动态、实时的监测技术,其近几年得到高速发展并在工程结构的安全和可靠性评估中发挥着日益重要的作用。Lamb波是一种板类结构中传播的超声导波,因其能够传播较远距离并且对结构表面和内部损 伤均敏感的优点,在结构健康监测领域得到日益广泛的应用,Lamb波监测技术已成为一个研究热点。Lamb波具有多模和频散特性,即使对于单模式的Lamb波信号,频散效应也使其波包发生扩展和变形,幅值随之减小,这种现象会随着传播距离的增大表现得越明显,严重降低了传感信号的分辨率和信噪比。现有技术中,时间反转方法无需Lamb波在结构中传播的先验知识便可自动补偿频散特性,但同时也消除了 Lamb波的传播时间,为后续损伤识别增加了难度。Alleyne通过信号再次激励来抑制传播距离已知的特定模态的频散现象。Sicard等基于后向传播函数提出一种频散补偿方法,该方法需要积分运算,计算量较大。Wilcox通过把信号从时域变换到空间域来消除频散特性,由于未修正处理过程中对激励信号波数谱的影响,使处理结果发生频偏现象,为信号解释造成困难。 Sicard R,Goyette J,Zellouf D. A numerical dispersioncompensation technique for time recompression of Lamb wave signals. Ultrasonics,2002,40(1-8) :727-732.Wilcox PD. A rapid signal processing technique to remove theeffect of dispersion from guided wave signals. IEEE Transactions on Ultrasonics,Ferroelectrics, and Frequency Control,2003,50 (4) :419-427. Xu B,Yu L,Giurgiutiu V. Lamb Wave Dispersion Compensationin Piezoelectric Wafer Active Sensor Phased-Array Applications In KunduTribikram(ed. ). Health Monitoring of Structural and Biological Systems 2009.Proceedings ofthe SPIE,2009,7295.。
技术实现思路
本专利技术的目的是克服现有技术的不足之处,提供一种具有频偏修正效果的时间-距离域映射方法。本专利技术的具有频偏修正效果的时间-距离域映射方法,包括下列步骤 I)求取频散波数和非频散波数关系数值求解或实验测得Lamb波信号的频散波数关系Κ(ω),并计算非频散波数关系Kmax (ω):Kmax( ) = /cgo 其中,ω为解频率,Cg0为信号在中心频率处的群速度; 2)计算原始频散信号频谱;对原始频散信号v(t)进行傅立叶变换得到其频谱为ν(ω); 3)求取信号波数谱 首先,根据Κ(ω),计算出插值映射序列Ω(10 ; Ω (k) = Γ1 ( ω ) 其中,IT1(Co)为Κ(ω)的逆函数。 其次,对ν(ω)在ω = Ω (k)处进行插值处理,得到信号波数谱V' (k)V' (k) = ν 其中,ν为对ν(ω)在ω = Ω (k)处的插值处理结果 4)修正激励信号波数谱 首先,计算插值映射序列ΩηΧη(10 il . - K K+ :. 其中,λ. 〃+> 为ΚηΧη(ω)的逆函数; 其次,计算信号波数谱的修正因子Cmax (k);C-(k) =V0/Va 其中,为对Κω)在ω = ΩηΧη(10处的插值处理结果=VJQ (k)]为对Vq(Co)在ω = Ω (k)处的插值处理结果; 再其次,计算考虑激励信号修正的信号波数谱V (k) v(k) =r (k)Cmax(k); 5)计算得到频偏修正的非频散距离域信号;对¥(10进行逆傅立叶变换得到非频散的距离域信号V (r)。 本专利技术的方法可补偿Lamb波信号中的频散效应,再压缩信号中因频散发生扩展和变形的各个波包,提高信号分辨率和信噪比。而且,在处理过程中对激励信号波数谱进行了修正,从而保持频散补偿前后信号波数成份不变,便于信号解释。附图说明 图I是具有频偏修正效果的时间-距离域映射方法的流程 图2是压电片P1和P2在铝板中的位置和Atl模式传感信号的传播路径示意 图3是激励信号; 图4是原始Atl模式传感信号; 图5是频散波数关系Κ(ω)和非频散波数关系曲线Κ_(ω); 图6是插值映射序列ΩηΧη(10和Ω (k); 图7是经具有须偏修正效果的时间-距离域映射方法处理后的A0模式传感信号; 图8是经未考虚频偏修正的时间-距离域映射方法处理后的A0模式传感信号。具体实施例方式下面结合附图和具体实施方式来对本专利技术作进一步地详细说明。实施例I 本实施例采用2024铝板作为待测结构,尺寸为ImX I. 5mX I. 5mm。该结构材料参数见表I。在结构中布置有2个压电片P1和P2,在铝板结构中的位置如图2所示。下面将以Lamb波Atl模式传感信号为例给出本专利技术的频散补偿方法的实现过程。为了在结构中激发Atl模式为主的Lamb波监测信号,选择中心频率为IOOkHz的三波峰正弦调制信号作为激励波形,如图3所示。以P1作为激励,P2作为传感器,采集到的传感信号如图4,传感信号中每个波包的传播路径如图2所示。在图4显示的原始传感信号中,第一个波包为幅值较大的A0模式直达信号,其时域宽度从原来的30 μ s (如图3所示)增大到60 μ S,波包形状不再对称,幅值也由于波包拉伸而下降,随着传播距离的增大这种频散效应表现得更严重,使得Atl模式边界反射信号发生了严重混叠,给信号分析增加了难度。 表I 2024铝板材料参数权利要求1.一种具有频偏修正效果的时间一 距离域映射方法,其特征在于,包括下列步骤I)求取频散波数和非频散波数关系 数值求解或实验测得Lamb波信号的频散波数关系10 ),并计算非频散波数关系 ΚΜ ) = ω ^ 其中,为角频率,Cfi0为信号在中心频率处的群速度;2)计算原始频散信号频谱对原始频散信号@进行傅立叶变换得到其频谱为 3)求取信号波数谱 首先,根据,计算出插值映射序列_ O0L)=IT1Cffl) 其中,Jr1(O)为jr(<p)的逆函数, 其次,对FXfiPj在I = 处进行插值处理,得到信号波数谱#( :) Z(Jt)=F 为对在β = Ω( )处的插值处理结果;4)修正激励信号波数谱 首先,计算插值映射序列《^(*5 本文档来自技高网
...

【技术保护点】
一种具有频偏修正效果的时间-距离域映射方法,其特征在于,包括下列步骤:????1)求取频散波数和非频散波数关系数值求解或实验测得Lamb波信号的频散波数关系????????????????????????????????????????????????,并计算非频散波数关系:????其中,为角频率,为信号在中心频率处的群速度;????2)计算原始频散信号频谱:对原始频散信号进行傅立叶变换得到其频谱为;3)求取信号波数谱:首先,根据,计算出插值映射序列:????其中,为的逆函数,其次,对在处进行插值处理,得到信号波数谱:????其中,为对在处的插值处理结果;????4)修正激励信号波数谱:首先,计算插值映射序列:????其中,为的逆函数;其次,计算信号波数谱的修正因子:????其中,为对在处的插值处理结果;为对在处的插值处理结果;再其次,计算考虑激励信号修正的信号波数谱:????;????5)计算得到频偏修正的非频散距离域信号:对进行逆傅立叶变换得到非频散的距离域信号。106773dest_path_image001.jpg,989278dest_path_image002.jpg,341762dest_path_image003.jpg,471261dest_path_image004.jpg,909196dest_path_image005.jpg,329813dest_path_image006.jpg,802382dest_path_image007.jpg,853515dest_path_image001.jpg,778746dest_path_image008.jpg,3054dest_path_image009.jpg,64551dest_path_image010.jpg,34387dest_path_image001.jpg,446914dest_path_image007.jpg,412596dest_path_image011.jpg,594179dest_path_image012.jpg,970802dest_path_image013.jpg,870625dest_path_image014.jpg,436736dest_path_image007.jpg,738404dest_path_image011.jpg,36661dest_path_image015.jpg,423780dest_path_image016.jpg,793582dest_path_image017.jpg,372593dest_path_image002.jpg,169648dest_path_image018.jpg,778483dest_path_image019.jpg,951976dest_path_image020.jpg,900340dest_path_image021.jpg,868296dest_path_image022.jpg,964428dest_path_image023.jpg,862983dest_path_image021.jpg,728171dest_path_image011.jpg,867028dest_path_image024.jpg,450456dest_path_image025.jpg,637855dest_path_image024.jpg,623129dest_path_image026.jpg...

【技术特征摘要】

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

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1