CN107884058B - 一种复杂结构Lamb波群速度频域测量方法 - Google Patents

一种复杂结构Lamb波群速度频域测量方法 Download PDF

Info

Publication number
CN107884058B
CN107884058B CN201710951005.6A CN201710951005A CN107884058B CN 107884058 B CN107884058 B CN 107884058B CN 201710951005 A CN201710951005 A CN 201710951005A CN 107884058 B CN107884058 B CN 107884058B
Authority
CN
China
Prior art keywords
lamb wave
signal
wave
scattered
lamb
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710951005.6A
Other languages
English (en)
Other versions
CN107884058A (zh
Inventor
蔡建
袁慎芳
汪晓鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710951005.6A priority Critical patent/CN107884058B/zh
Publication of CN107884058A publication Critical patent/CN107884058A/zh
Application granted granted Critical
Publication of CN107884058B publication Critical patent/CN107884058B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H5/00Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
    • G01H11/06Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
    • G01H11/08Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means using piezoelectric devices

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提供了一种复杂结构Lamb波群速度频域测量方法,属于Lamb波结构健康监测领域。该方法包括下列步骤:(1)在结构中布置压电片并采集结构Lamb波参考信号;(2)在结构中引入散射体并采集当前结构状态下的Lamb波传感信号;(3)根据Lamb波传感信号及其参考信号求取Lamb波散射信号;(4)从Lamb波散射信号中提取所选模式的散射波包;(5)计算Lamb波激励信号和散射波包的频谱;(6)求取Lamb波激励信号和散射波包的相位差谱;(7)由相位差谱计算得到Lamb波所选模式的群速度。本发明避免了传统时域测量方法中Lamb波频散的影响,提高了测量的准确性。

Description

一种复杂结构Lamb波群速度频域测量方法
技术领域
本发明涉及一种复杂结构Lamb波群速度频域测量方法,属于Lamb波结构健康监测领域。
背景技术
结构健康监测技术对于提高工程结构安全性、降低维护费用具有重要意义。Lamb波作为薄板类结构中的一种超声导波,其传播距离远并对结构内部和表面损伤均敏感,所以广泛应用于结构健康监测中。Lamb波具有多模和频散特性。多模特性是指结构中至少存在两种Lamb波模式。按照Lamb波传播时结构质点振动相对于结构中界面是否为对称和反对称的特点,将Lamb波模式分为对称模式和反对称模式两类,而且每类模式又分为不同阶数。不同阶数的对称和反对称模式分别表示为S0,S1,S2...Sn和A0,A1,A2...An,其中n为阶数。Lamb波频散特性是指其相速度和群速度随频率变化而变化的特性。
在Lamb波结构健康监测方法的实际应用中,Lamb波的传播速度,尤其是群速度作为一种重要参数常需要测量。目前常用的群速度测量方法为时域测量方法,该方法一般需要提取Lamb波信号中直达波包的波达时间。然而,Lamb波频散容易造成信号波包发生时域扩展和变形,使Lamb波直达波包波达时间的准确提取变得困难,并最终降低群速度测量结果的准确性。Lamb波频散越严重,基于上述传统时域测量方法测得的群速度值越不准确。除此之外,实际工程结构较为复杂,边界反射严重,Lamb波传感信号中所选模式的直达波包很容易与边界反射波包发生混叠,难以提取纯净的直达波包进行该模式的群速度测量。
综述所述,复杂结构中Lamb波的群速度测量成为Lamb波工程应用中急需解决的一个实际问题。
发明内容
本发明提出了一种复杂结构Lamb波群速度频域测量方法,该方法通过频域求取Lamb波激励信号和散射波包的相位差谱,避免了传统时域测量方法中Lamb波频散的影响,提高了测量的准确性。另一方面,该方法基于Lamb波散射信号进行,而在通过比较Lamb波传感信号及其参考信号求取散射信号的过程中已大部分消除了边界反射信号,从而降低了复杂结构边界反射对群速度测量的影响。
本发明为解决其技术问题采用如下技术方案:
一种复杂结构Lamb波群速度频域测量方法,包括下列步骤:
(1)在结构中布置压电片并采集结构Lamb波参考信号
在结构中布置两个压电片分别作为激励和传感器,将窄带激励信号va(t)加载在激励上向结构中激发出Lamb波信号,并利用传感器采集结构中没有引入散射体情况下的参考信号v0(t);
(2)在结构中引入散射体并采集当前结构状态下的Lamb波传感信号
在结构中引入单个散射体,并利用传感器采集结构当前状态下的Lamb波传感信号v1(t);
(3)根据Lamb波传感信号及其参考信号求取Lamb波散射信号
对所述传感信号v1(t)和参考信号v0(t)进行信号差运算,得到对应于该散射体的Lamb波散射信号s(t);
(4)从Lamb波散射信号中提取所选模式的散射波包
利用窗函数从所述Lamb波散射信号中提取出散射体产生的所选模式的散射波包;
(5)分别计算Lamb波激励信号和散射波包的频谱
分别对所述Lamb波激励信号va(t)和散射波包进行傅里叶变换得到两者的频谱Va(ω)和SW(ω),其中Va(ω)和SW(ω)分别为va(t)和散射波包的频谱,ω表示角频率;
(6)求取Lamb波激励信号和散射波包的相位差谱
利用Va(ω)和SW(ω),将所述Lamb波激励信号va(t)和散射波包的相位差谱先计算为
其中Φ1(ω)为未进行相位展开的Lamb波激励信号va(t)和散射波包的相位差谱,arctan[]表示反正切计算。然后对Φ1(ω)进行相位展开处理得到绝对或相对相位差谱Φ(ω)。
(7)根据相位差谱计算得到Lamb波所选模式的群速度
利用Φ(ω)将Lamb波所选模式的群速度计算为
其中cg(ω)为Lamb波所选模式在角频率ω下的群速度,L为该模式信号从激励经散射体传播至传感器整个散射路径的距离。
所述的步骤(2)中在结构中引入散射体的位置应该在所述激励和传感器的直达路径上。
本发明的有益效果如下:
(1)本发明群速度测量中通过计算相位差谱来求取所选模式的群速度,避免了传统时域测量方法中计算信号飞行时间时,Lamb波频散对确定信号波达时间的影响,提高了测量的准确性。
(2)本发明利用Lamb波散射波包计算群速度,可减小实际复杂多反射结构中边界反射信号对测量的影响,提高了测量的准确性。
(3)本发明可根据某一Lamb波传感信号一次性计算出信号中心频率附近多个频率点的群速度值,相比于传统时域测量方法每次只能测量信号中心频率处的群速度值,测量效率高。
附图说明
图1是复杂结构Lamb波群速度频域测量方法的实施流程图。
图2是实施例提供的铝板结构中压电片和螺母的分布情况图。
图3是实施例提供的激励信号波形图。
图4是实施例提供的铝板结构中未粘接螺母时的Lamb波参考信号波形图。
图5是实施例提供的铝板结构中粘接有螺母时的Lamb波传感信号波形图。
图6是实施例提供的Lamb波散射信号波形图。
图7是实施例提供的从Lamb波散射信号中提取的A0模式散射波包图。
图8是实施例提供的相位展开后的激励信号和A0模式散射波包的相对相位差谱图。
图9是实施例提供的A0模式群速度频域测量值、时域测量值和理论值。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下将参照本发明实施例中的附图,通过实施方式清楚、完整地描述本发明的技术方案,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了便于验证,本实施例采用LY21铝板结构,尺寸为800mm×800mm×4mm,材料参数见表1。
表1.铝板材料参数
厚度d(mm) 密度ρ(kg/m<sup>3</sup>) 泊松比μ 弹性模量E(GPa)
4 2780 0.33 73.1
本实施例的复杂结构中Lamb波的群速度测量方法如图1所示,包括下列步骤:
(1)在结构中布置压电片并采集结构Lamb波参考信号
为了考察复杂多反射结构情况,在铝板的左下角,并靠近左边界和下边界处布置有2个压电片P1和P2分别作为激励和传感器,两个压电片在铝板结构中的具体位置如图2所示,其中两个圆表示分别表示相距200mm的压电片P1和P2。加载到P1上的激励信号为1.5波峰正弦调制信号va(t),如图3所示,中心频率选为100kHz以向结构中激发出A0模式为主的Lamb波信号,同时利用P2采集到结构中没有引入散射体情况下的参考信号v0(t),如图4所示,可在60-100μs时间段内看到,A0模式直达波包和边界反射波包发生严重混叠。
(2)在结构中引入散射体并采集当前结构状态下的Lamb波传感信号
为了产生能量尽量大的散射波包,以及尽可能降低二次边界反射波包的能量。在P1和P2直达路径的中间用502胶粘剂一个边长为1cm的六边形螺母D作为散射体,如图2中的六边形所示。向P1加载图3所示的激励信号va(t),并利用P2采集结构当前状态下的Lamb波传感信号v1(t),如图5所示。
(3)根据Lamb波传感信号及其参考信号求取Lamb波散射信号
对Lamb波传感信号v1(t)和参考信号v0(t)进行信号差运算,得到对应于螺母的Lamb波散射信号s(t),如图6所示,可在60-100μs时间段内看到边界反射信号被有效抑制,使A0模式散射波包变得更为纯净。
(4)从Lamb波散射信号中提取所选模式的散射波包
根据Lamb波A0模式散射波包的时域位置,确定合适的窗函数,从图6中的Lamb波散射信号s(t)中提取出由散射体产生的A0模式散射波包,如图7所示。
(5)分别计算Lamb波激励信号和散射波包的频谱
基于快速傅里叶变换,分别对计算出Lamb波激励信号va(t)和所提取的A0模式散射波包的频谱Va(ω)和SW(ω)。
(6)求取Lamb波激励信号和散射波包的相位差谱
将Lamb波激励信号va(t)和A0模式散射波包的相位差谱计算为再经过相位展开处理得到相对相位差谱,如图8所示。
(7)根据相位差谱计算得到Lamb波所选模式的群速度
L=20cm,根据计算得到Lamb波A0模式在中心频率100kHz附近的群速度,如图9所示,图中的理论值根据表1中的材料参数数值计算得到,时域测量值则依次利用中心频率分别为80~130kHz(间隔为1KHz)的散射信号进行了51次时域群速度测量得到。
从图9可看到,在80-130kHz频率范围内,频域测量结果在A0模式的理论值附近波动,两者非常吻合。而对于时域测量结果,虽然所利用的散射信号中边界反射已得到抑制,但由于A0模式频散对其散射波包的影响,使得最终时域测量结果与理论值偏差很大。而且,相比于通过1次频域测量就得得到80~130kHz频率范围内的群速度曲线,时域测量方法则需要进行51次测量,测量效率较低。
本发明的基本原理是:传统Lamb波群速度测量方法中将群速度计算为cg(ω)=R/Ttof,其中Ttof为Lamb波模式信号从激励传播至传感器的飞行时间(Time of Flight,TOF),R为传播距离。Ttof为Lamb波传感信号中所选模式直达波包和激励信号波包波达时间(Timeof Arrival,TOA)之差。由于信号波包的TOA一般定义为信号包络的最大值,而频散传感信号中的直达波包由于波包时域扩展和变形,使其包络顶点的时域位置相比于理想非频散下的情况发生变化,从而影响了最终群速度测量结果的准确性。
本发明则是将传统时域测量方法中时域确定Ttof替换为频域求取Lamb波激励信号和散射波包之间的绝对或相对相位差谱,故能避免频散影响信号时域波形引起的问题。
另一方面,本发明没有直接利用Lamb波传感信号,而是利用散射信号进行群速度测量。散射信号是通过对比结构引入散射体前后的传感信号和参考信号得到。由于两者均包含共同的边界反射信号,所以通过对比可以大大降低边界反射的影响。

Claims (3)

1.一种复杂结构Lamb波群速度频域测量方法,其特征在于,包括下列步骤:
(1)在结构中布置压电片并采集结构Lamb波参考信号;具体过程如下:在结构中布置两个压电片分别作为激励和传感器,将窄带激励信号va(t)加载在激励上向结构中激发出Lamb波信号,并利用传感器采集结构中没有引入散射体情况下的参考信号v0(t);
(2)在结构中引入散射体并采集当前结构状态下的Lamb波传感信号;所述在结构中引入散射体的位置应该在所述激励和传感器的直达路径上;
(3)根据Lamb波传感信号及其参考信号进行信号差运算,求取得到对应于该散射体的Lamb波散射信号;
(4)利用窗函数从Lamb波散射信号中提取出散射体产生的所选模式的散射波包;
(5)分别计算Lamb波激励信号和散射波包的频谱;
(6)求取Lamb波激励信号和散射波包的相位差谱;具体过程如下:
利用Va(ω)和SW(ω),将所述Lamb波激励信号va(t)和散射波包的相位差谱先计算为
其中Φ1(ω)为未进行相位展开的Lamb波激励信号va(t)和散射波包的相位差谱,arctan[]表示反正切计算,然后对Φ1(ω)进行相位展开处理得到绝对或相对相位差谱Φ(ω);
(7)根据相位差谱计算得到Lamb波所选模式的群速度,具体过程如下:
利用Φ(ω)将Lamb波所选模式的群速度计算为:
其中cg(ω)为Lamb波所选模式在角频率ω下的群速度,L为该模式信号从激励经散射体传播至传感器整个散射路径的距离。
2.根据权利要求1所述的一种复杂结构Lamb波群速度频域测量方法,其特征在于,步骤(1)所述结构为铝板结构,尺寸为800mm×800mm×4mm。
3.根据权利要求1所述的一种复杂结构Lamb波群速度频域测量方法,其特征在于,步骤(2)的具体过程如下:在结构中引入单个散射体,并利用传感器采集结构当前状态下的Lamb波传感信号。
CN201710951005.6A 2017-10-13 2017-10-13 一种复杂结构Lamb波群速度频域测量方法 Active CN107884058B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710951005.6A CN107884058B (zh) 2017-10-13 2017-10-13 一种复杂结构Lamb波群速度频域测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710951005.6A CN107884058B (zh) 2017-10-13 2017-10-13 一种复杂结构Lamb波群速度频域测量方法

Publications (2)

Publication Number Publication Date
CN107884058A CN107884058A (zh) 2018-04-06
CN107884058B true CN107884058B (zh) 2019-06-25

Family

ID=61781626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710951005.6A Active CN107884058B (zh) 2017-10-13 2017-10-13 一种复杂结构Lamb波群速度频域测量方法

Country Status (1)

Country Link
CN (1) CN107884058B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208387B (zh) * 2019-05-23 2021-11-23 江苏大学 一种基于焊缝特征导波柔性传感器的信号聚焦方法
CN111175379B (zh) * 2020-01-21 2022-08-02 南京邮电大学 一种基于同步压缩小波变换的Lamb波板结构健康监测方法
CN111780856B (zh) * 2020-06-01 2022-03-29 哈尔滨工业大学 基于瑞利散射谱的相位谱分析的光纤分布式振动测量方法
CN111812579B (zh) * 2020-09-02 2023-03-14 南昌大学 超精密渡越时间测量方法及系统
CN113504314B (zh) * 2021-05-17 2022-10-11 南京航空航天大学 环境变化场合下的高可靠性超声导波损伤因子计算方法
CN113295933A (zh) * 2021-05-25 2021-08-24 云南电网有限责任公司电力科学研究院 一种雷电电磁波传播的群速度提取方法
CN113588414B (zh) * 2021-08-03 2023-04-25 重庆大学 基于超声频谱能量衰减的螺栓轴向应力检测方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6330924A (ja) * 1986-07-25 1988-02-09 Hitachi Ltd 弾性波応用座標入力装置
US6360609B1 (en) * 2000-02-23 2002-03-26 Massachusetts Institute Of Technology Method and system for interpreting and utilizing multimode dispersive acoustic guided waves
CN101451977A (zh) * 2008-12-30 2009-06-10 南京航空航天大学 一种基于双元压电传感阵列和时间窗函数的无基准Lamb波损伤监测方法
CN102243214A (zh) * 2011-04-25 2011-11-16 中国人民解放军后勤工程学院 一种基于群速度失配的超声兰姆波二次谐波时域信号的测量方法
CN102565201A (zh) * 2011-10-25 2012-07-11 中国人民解放军理工大学 一种基于测量波数曲线的Lamb波频散补偿方法
CN103344709A (zh) * 2013-06-06 2013-10-09 南京航空航天大学 基于线性映射频散补偿的压电超声相控阵多损伤成像方法
CN103808810A (zh) * 2014-03-14 2014-05-21 哈尔滨工业大学 一种基于小波时频图和群速度曲线的钢轨裂纹声发射特征提取方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO327139B1 (no) * 2006-05-30 2009-05-04 Clampon As Fremgangsmate og system for bestemmelse av tap i materialtykkelse i en fast struktur

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6330924A (ja) * 1986-07-25 1988-02-09 Hitachi Ltd 弾性波応用座標入力装置
US6360609B1 (en) * 2000-02-23 2002-03-26 Massachusetts Institute Of Technology Method and system for interpreting and utilizing multimode dispersive acoustic guided waves
CN101451977A (zh) * 2008-12-30 2009-06-10 南京航空航天大学 一种基于双元压电传感阵列和时间窗函数的无基准Lamb波损伤监测方法
CN102243214A (zh) * 2011-04-25 2011-11-16 中国人民解放军后勤工程学院 一种基于群速度失配的超声兰姆波二次谐波时域信号的测量方法
CN102565201A (zh) * 2011-10-25 2012-07-11 中国人民解放军理工大学 一种基于测量波数曲线的Lamb波频散补偿方法
CN103344709A (zh) * 2013-06-06 2013-10-09 南京航空航天大学 基于线性映射频散补偿的压电超声相控阵多损伤成像方法
CN103808810A (zh) * 2014-03-14 2014-05-21 哈尔滨工业大学 一种基于小波时频图和群速度曲线的钢轨裂纹声发射特征提取方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于Lamb波群速度修正的复合材料结构损伤定位";高飞 等;《航空制造技术》;20171001(第19期);第60-66页

Also Published As

Publication number Publication date
CN107884058A (zh) 2018-04-06

Similar Documents

Publication Publication Date Title
CN107884058B (zh) 一种复杂结构Lamb波群速度频域测量方法
CN103344709A (zh) 基于线性映射频散补偿的压电超声相控阵多损伤成像方法
Ding et al. A new method for waveform analysis for estimating AE wave arrival times using wavelet decomposition
CN103698748A (zh) 复合材料中基于信号波速与衰减补偿的2d-music冲击定位方法
CN106768109A (zh) 基于回波上升段峰值拟合和基于回波能量点定位的气体超声波流量计信号处理方法
CN102565201B (zh) 一种基于测量波数曲线的Lamb波频散补偿方法
CN206905879U (zh) 一种分离式叶片固有频率检测装置
CN102435844A (zh) 一种频率无关的正弦信号相量计算方法
CN101226070B (zh) 一种基于外差式随机共振的涡街频率检测方法
CN105300508A (zh) 一种超声波信号传播时间的双阈值检测电路及检测方法
CN102589627A (zh) 一种用于超声波流量计的绝对传播时间测量方法
CN101874744B (zh) 用于长骨分析的超声导波参数测量方法
CN110068613B (zh) 一种结构导波响应群速度频散测试方法
CN102866261B (zh) 检测超声波在测流速中飞行时间的方法
Zhu et al. A two-step impact localization method for composite structures with a parameterized laminate model
CN114459649A (zh) 一种基于压电换能器阵列的无基线数据平面应力场在线监测方法、系统、设备和介质
CN105717198A (zh) 一种面向结构冲击定位的单频-重估计music方法
CN104748704A (zh) 薄壁结构超声共振测厚频谱分析内插校正方法
Mu et al. Echo signal envelope fitting based signal processing methods for ultrasonic gas flow-meter
CN102590598B (zh) 基于多阈值比较的超声波信号过零点预测方法
CN103075981A (zh) 一种超声波测厚方法
CN105841794B (zh) 基于优化的耦合型光纤传感器平板超声源定位方法及系统
CN102565541A (zh) 一种用于电学层析成像系统的递推解调方法
CN104458913A (zh) 材料性能退化的非线性导波评估方法与装置
CN103743818B (zh) 基于波的能流图的损伤诊断方法和实施该方法的损伤诊断系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant