CN105486759A - 一种Lamb波时频能量密度析出走时提取方法 - Google Patents
一种Lamb波时频能量密度析出走时提取方法 Download PDFInfo
- Publication number
- CN105486759A CN105486759A CN201610018399.5A CN201610018399A CN105486759A CN 105486759 A CN105486759 A CN 105486759A CN 201610018399 A CN201610018399 A CN 201610018399A CN 105486759 A CN105486759 A CN 105486759A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- energy density
- lamb wave
- discrete
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0234—Metals, e.g. steel
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0427—Flexural waves, plate waves, e.g. Lamb waves, tuning fork, cantilever
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/26—Scanned objects
- G01N2291/263—Surfaces
- G01N2291/2632—Surfaces flat
Landscapes
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种Lamb波时频能量密度析出走时提取方法,包括步骤:在待测材料检测区域两侧分别设置发射EMAT和接收EMAT,令发射端激发所需模态Lamb波,对侧接收端接收Lamb波;得到检测波形后,对检测波形数据进行离散短时傅里叶变换,得到检测波形的时频能量密度分布;根据时频能量密度分布,建立基于数据中心频率的时频能量密度析出方法,提取出数据中心频率下能量密度随时间变化的离散曲线;通过对离散曲线的拟合建立时域能量密度曲线,提取时域能量密度曲线各峰值所对应的时间,作为所需模态Lamb波的走时。本发明具有如下优点:对窄带Lamb波走时提取精度高、计算高效快速。
Description
技术领域
本发明涉及无损检测技术领域,具体涉及一种Lamb波时频能量密度析出走时提取方法。
背景技术
利用超声Lamb波缺陷层析成像技术可以快速、有效地获得缺陷的轮廓和尺寸等具体信息,该技术继承了传统超声Lamb波检测的诸多优点。Lamb层析成像技术在缺陷两侧并排布置换能器阵列,从检测的Lamb波波形中提取投影数据,根据投影数据和一定的重建算法,反演出缺陷的形状、尺寸等具体信息。Lamb波走时是常用的投影数据,是指Lamb波在收发换能器对之间的传播时间。Lamb走时作为层析成像重建算法的输入量,其提取的准确度直接决定了缺陷成像质量及缺陷尺寸量化精度。传统的Lamb波走时提取包括基于Lamb波时域波形峰值的提取方法、希尔伯特-黄变换(Hilbert-huangtransform,HHT)以及其他时频分析方法等,其走时提取的准确度不高,尤其是在走时提取局部区域的波形振动剧烈,无法确定真正意义上的走时,并带有很多人为因素,给走时提取带来很大误差,严重破坏了缺陷重建质量以及降低了缺陷尺寸量化精度。上述问题是制约Lamb波层析成像技术发展及缺陷成像质量的一个瓶颈问题。中国专利文献公开了“一种地质雷达波速层析探测的走时数据快速采集方法”,该技术涉及一种地质雷达波速层析探测的走时数据快速采集方法,对地质雷达波速层析探测的走时数据进行快速、连续的采集,但仍然采用时域波形峰值时间的走时提取方式,走时提取的准确度低,误差较大,降低了层析探测的精度。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的一个目的在于提出一种Lamb波时频能量密度析出走时提取方法。
为了实现上述目的,本发明的实施例公开了一种Lamb波时频能量密度析出走时提取方法,包括如下步骤:S1:在待测材料检测区域两侧分别设置发射EMAT和接收EMAT;S2:设定Lamb波检测的工作点,设定所述工作点的工作频率为fc,令发射端激发所需模态Lamb波,对侧接收端接收Lamb波,对接收到的Lamb波检测信号进行放大和窄带滤波,滤波中心频率为fc,然后再进行采集和存储,得到检测波形数据;S3:利用离散短时傅里叶变换对检测波形数据进行时频分析,得到检测数据的时频能量密度分布,时频能量密度的时间坐标为ti,频率坐标为fk,在某一组确定的时间点和频率点(ti,fk)处的能量密度为E(ti,fk);其中,i=1,2,…,M,k=1,2,…,N,M、N为正整数;S4:在时频能量密度分布中,对于某一离散的时间点ti,时间点ti对应的频率坐标为fk,其中k为自然数,判断是否存在频率点fk的值与中心频率fc的值相等,若是,提取出对应于中心频率fc处的能量密度E(ti,fc),并进行S7;若不是,进行S5;S5:对于离散时间点ti对应的频率坐标fk,在中心频率fc两侧分别提取与其最接近的频率点fl和fh,其中fl<fc<fh,并分别提取出在(ti,fl)和(ti,fh)处的能量密度E(ti,fl)和E(ti,fh);S6:根据E(ti,fl)和E(ti,fh),采用线性插值法,求取对应于离散时间点ti中心频率fc的能量密度E(ti,fc);S7:判断是否所有的离散时间点ti对应于中心频率fc处的能量密度E(ti,fc)均已获得,若是,则进行S8;若不是,则判断离散时间点变为ti+1,并返回S4;S8:利用所有的离散时间点t和中心频率处的能量密度E(t,fc),通过对离散时间点处能量密度的拟合,建立中心频率的时域能量密度曲线E(t,fc);S9:提取时域能量密度曲线E(t,fc)各峰值所对应的时间tp,作为所需模态Lamb波的走时,并结束走时提取,其中,p=1,2,…,M1,M1为正整数且M1<M。
根据本发明实施例的Lamb波时频能量密度析出走时提取方法,通过求取检测波形数据的时频能量密度分布,建立数据中心频率时频能量密度析出方法,并采用建立的时域能量密度曲线准确提取Lamb波的走时,计算准确、高效、快速,解决了传统Lamb波走时提取方法提取精度低、导致缺陷层析成像质量差的问题。
另外,根据本发明上述实施例的Lamb波时频能量密度析出走时提取方法,还可以具有如下附加的技术特征:
进一步地,EMAT采用圆形回折线圈设计,通过工作频率和导线间距的匹配关系,激发和接收所需模态的Lamb波。
进一步地,用高阶滤波电路对接收EMAT得到的Lamb波信号进行窄带滤波。
进一步地,利用离散短时傅里叶变换对检测波形数据进行时频分析,得到检测数据的时频能量密度分布,设Lamb波检测波形数据为x(m),其中m=1,2,…,M,则其离散短时傅里叶变换为:
其中,DSTFT(i,k)为Lamb波检测波形数据的离散短时傅里叶变换结果;g(i)为离散短时傅里叶变换的窗函数。检测数据的时频能量密度分布为:E(ti,fk)=|DSTFT(i,k)|2。
进一步地,提取中心频率两侧频率点的方法为:fl=argmin|fk-fc|,fk<fc,fh=argmin|fk-fc|,fk>fc。
进一步地,采用线性插值法求取离散的时间点ti对应中心频率处的能量密度E(ti,fc):
进一步地,通过最小二乘拟合建立基于中心频率的时域能量密度曲线E(t,fc),设拟合函数为时频能量密度为:
进一步地,提取时域能量密度曲线的峰值时间作为所需模态Lamb波的走时tp:tp=argmaxE(t,fc)。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明一个实施例的Lamb波时频能量密度析出走时提取流程图;
图2是本发明一个实施例的实验结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
参照下面的描述和附图,将清楚本发明的实施例的这些和其他方面。在这些描述和附图中,具体公开了本发明的实施例中的一些特定实施方式,来表示实施本发明的实施例的原理的一些方式,但是应当理解,本发明的实施例的范围不受此限制。相反,本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
以下结合附图描述根据本发明实施例的Lamb波时频能量密度析出走时提取方法。
图1是本发明一个实施例的Lamb波时频能量密度析出走时提取流程图。请参考图1,一种Lamb波时频能量密度析出走时提取方法,包括如下步骤:
S1:在待测材料检测区域两侧分别设置发射EMAT和接收EMAT。
在本发明的一个示例中,取一件厚度为1~10mm的待测钢板,在钢板两侧分别布置发射EMAT和接收EMAT,两EMAT直径为10~80mm,两EMAT中心间距为800mm。
S2:设定Lamb波检测的工作点,设定所述工作点的工作频率为fc,令发射端激发所需模态Lamb波,对侧接收端接收Lamb波,对接收到的Lamb波检测信号进行放大和窄带滤波,滤波中心频率为fc,然后再进行采集和存储,得到检测波形数据。
接着S1相应的示例,用射频功率放大器激励发射EMAT,用接收EMAT接收Lamb波,激发频率为256kHz;采用窄带滤波电路对接收到的Lamb波检测信号进行滤波,滤波中心频率为256kHz,再进行采集和存储,得到检测波形数据。
S3:利用离散短时傅里叶变换对检测波形数据进行时频分析,得到检测数据的时频能量密度分布,时频能量密度的时间坐标为ti,频率坐标为fk,在某一组确定的时间点和频率点(ti,fk)处的能量密度为E(ti,fk);其中,i=1,2,…,M,k=1,2,…,N,M、N为正整数。
接着S2相应的示例,利用离散短时傅里叶变换(Discreteshort-timeFouriertransform,DSTFT)对检测波形数据进行时频分析,时频能量密度的时间坐标为ti(i=1,2,…,M),频率坐标为fk(k=1,2,…,N),在某一组确定的时间点和频率点(ti,fk)处的能量密度为E(ti,fk);其中,M、N为正整数;设Lamb波检测波形数据为x(m),其中m=1,2,…,M,则其离散短时傅里叶变换为:
其中,DSTFT(i,k)为Lamb波检测波形数据的离散短时傅里叶变换结果;g(i)为离散短时傅里叶变换的窗函数。
求取检测数据的时频能量密度分布
E(ti,fk)=|DSTFT(i,k)|2。
S4:在时频能量密度分布中,对于某一离散的时间点ti,时间点ti对应的频率坐标为fk,其中k为自然数,判断是否存在频率点fk的值与中心频率fc的值相等,若是,提取出对应于中心频率fc处的能量密度E(ti,fc),并进行S7;若不是,进行S5。
S5:对于离散时间点ti对应的频率坐标fk,在中心频率fc两侧分别提取与其最接近的频率点fl和fh,其中fl<fc<fh,并分别提取出在(ti,fl)和(ti,fh)处的能量密度E(ti,fl)和E(ti,fh)。
接着S4,对于离散时间点ti对应的频率坐标fk(k=1,2,…,N),在中心频率fc两侧分别提取与其最接近的频率点fl和fh,其中fl<fc<fh
fl=argmin|fk-fc|,fk<fc
fh=argmin|fk-fc|,fk>fc
分别提取出在(ti,fl)和(ti,fh)处的能量密度E(ti,fl)和E(ti,fh)。
S6:根据E(ti,fl)和E(ti,fh),采用线性插值法,求取对应于离散时间点ti中心频率fc的能量密度E(ti,fc)。
接着S5的示例,根据E(ti,fl)和E(ti,fh),采用线性插值法,求取对应于离散时间点ti中心频率fc的能量密度E(ti,fc)
S7:判断是否所有的离散时间点ti对应于中心频率fc处的能量密度E(ti,fc)均已获得,若是,则进行S8;若不是,则判断离散时间点变为ti+1,并返回S4。
S8:利用所有的离散时间点t和中心频率处的能量密度E(t,fc),通过对离散时间点处能量密度的拟合,建立中心频率的时域能量密度曲线E(t,fc)。
接着S7,利用所有的离散时间点ti(i=1,2,…,M)和中心频率处的能量密度E(ti,fc),通过对离散时间点处能量密度的最小二乘拟合,建立中心频率的时域能量密度曲线E(t,fc),设拟合函数为时频能量密度为:
S9:提取时域能量密度曲线E(t,fc)各峰值所对应的时间tp,作为所需模态Lamb波的走时,并结束走时提取,其中,p=1,2,…,M1,M1为正整数且M1<M。
接着S8的示例,提取时域能量密度曲线E(t,fc)各峰值所对应的时间tp(p=1,2,…,M1)
tp=argmaxE(t,fc)
作为所需模态Lamb波的走时,并结束走时提取,其中,M1为正整数且M1<M。
本发明的Lamb波时频能量密度析出走时提取方法,还公开了另外一个实施例,具体如下:
步骤1:取一件厚度为4mm的待测铝板,在铝板两侧分别布置发射EMAT和接收EMAT,两EMAT直径均为40mm,两EMAT中心间距为800mm。
步骤2:用射频功率放大器AG124激励发射EMAT,产生A0模态Lamb波,用接收EMAT接收Lamb波,激发频率为128kHz;采用窄带滤波电路对接收到的Lamb波检测信号进行滤波,滤波中心频率为128kHz,再进行采集和存储,得到检测波形数据。
步骤3:利用离散短时傅里叶变换(Discreteshort-timeFouriertransform,DSTFT)对检测波形数据进行时频分析,时频能量密度的时间坐标为ti(i=1,2,…,M),频率坐标为fk(k=1,2,…,N),在某一组确定的时间点和频率点(ti,fk)处的能量密度为E(ti,fk);其中,M=636,N=458;设Lamb波检测波形数据为x(m),其中m=1,2,…,M,则其离散短时傅里叶变换为:
其中,DSTFT(i,k)为Lamb波检测波形数据的离散短时傅里叶变换结果;g(i)为离散短时傅里叶变换的窗函数,采用Hamming窗,窗函数时宽为159。
求取检测数据的时频能量密度分布
E(ti,fk)=|DSTFT(i,k)|2
步骤4:在时频能量密度分布中,对于某一离散的时间点ti,其对应的频率坐标为fk(k=1,2,…,N),判断是否存在频率点fk的值与中心频率fc=128kHz的值相等,若是,提取出对应于中心频率fc处的能量密度E(ti,fc),并进行步骤7;若不是,进行步骤5。
步骤5:对于离散时间点ti对应的频率坐标fk(k=1,2,…,N),在中心频率fc两侧分别提取与其最接近的频率点fl和fh,其中fl<fc<fh
fl=argmin|fk-fc|,fk<fc
fh=argmin|fk-fc|,fk>fc
分别提取出在(ti,fl)和(ti,fh)处的能量密度E(ti,fl)和E(ti,fh)。
步骤6:根据E(ti,fl)和E(ti,fh),采用线性插值法,求取对应于离散时间点ti中心频率fc的能量密度E(ti,fc)
步骤7:判断是否所有的离散时间点ti(i=1,2,…,M)对应于中心频率fc处的能量密度E(ti,fc)均已获得,若是,则进行步骤8;若不是,则考虑离散时间点ti+1,并返回步骤4。
步骤8:利用所有的离散时间点ti(i=1,2,…,M)和中心频率处的能量密度E(ti,fc),通过对离散时间点处能量密度的最小二乘拟合,建立中心频率的时域能量密度曲线E(t,fc),设拟合函数为时频能量密度为:
步骤9:提取时域能量密度曲线E(t,fc)各峰值所对应的时间tp(p=1,2,…,M1)
tp=argmaxE(t,fc)
作为所需模态Lamb波的走时,并结束走时提取,其中,M1为正整数且M1<M。
本实施例最后得到的Lamb波直达波(直接从发射EMAT到接收EMAT的Lamb波)的走时为281μs,A0模态Lamb波在4mm铝板、128kHz工作频率下的理论群速度为2944m/s,直达波传播距离为800mm,则其理论传播时间为272μs。走时提取得到的时间与理论传播时间的相对误差为3.3%,利用本发明得到的Lamb波走时精度较高。
另外,本发明实施例的Lamb波时频能量密度析出走时提取方法的其它构成以及作用对于本领域的技术人员而言都是已知的,为了减少冗余,不做赘述。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。
Claims (8)
1.一种Lamb波时频能量密度析出走时提取方法,其特征在于,包括如下步骤:
S1:在待测材料检测区域两侧分别设置发射EMAT和接收EMAT;
S2:设定Lamb波检测的工作点,设定所述工作点的工作频率为fc,令发射端激发所需模态Lamb波,对侧接收端接收Lamb波,对接收到的Lamb波检测信号进行放大和窄带滤波,滤波中心频率为fc,然后再进行采集和存储,得到检测波形数据;
S3:利用离散短时傅里叶变换对检测波形数据进行时频分析,得到检测数据的时频能量密度分布,时频能量密度的时间坐标为ti,频率坐标为fk,在某一组确定的时间点和频率点(ti,fk)处的能量密度为E(ti,fk);其中,i=1,2,…,M,k=1,2,…,N,M、N为正整数;
S4:在时频能量密度分布中,对于某一离散的时间点ti,时间点ti对应的频率坐标为fk,其中k为自然数,判断是否存在频率点fk的值与中心频率fc的值相等,若是,提取出对应于中心频率fc处的能量密度E(ti,fc),并进行S7;若不是,进行S5;
S5:对于离散时间点ti对应的频率坐标fk,在中心频率fc两侧分别提取与其最接近的频率点fl和fh,其中fl<fc<fh,并分别提取出在(ti,fl)和(ti,fh)处的能量密度E(ti,fl)和E(ti,fh);
S6:根据E(ti,fl)和E(ti,fh),采用线性插值法,求取对应于离散时间点ti中心频率fc的能量密度E(ti,fc);
S7:判断是否所有的离散时间点ti对应于中心频率fc处的能量密度E(ti,fc)均已获得,若是,则进行S8;若不是,则判断离散时间点变为ti+1,并返回S4;
S8:利用所有的离散时间点t和中心频率处的能量密度E(t,fc),通过对离散时间点处能量密度的拟合,建立中心频率的时域能量密度曲线E(t,fc);
S9:提取时域能量密度曲线E(t,fc)各峰值所对应的时间tp,作为所需模态Lamb波的走时,并结束走时提取,其中,p=1,2,…,M1,M1为正整数且M1<M。
2.根据权利要求1所述的方法,其特征在于:EMAT采用圆形回折线圈设计,通过工作频率和导线间距的匹配关系,激发和接收所需模态的Lamb波。
3.根据权利要求1所述的方法,其特征在于:用高阶滤波电路对接收EMAT得到的Lamb波信号进行窄带滤波。
4.根据权利要求1所述的方法,其特征在于:利用离散短时傅里叶变换对检测波形数据进行时频分析,得到检测数据的时频能量密度分布,设Lamb波检测波形数据为x(m),其中m=1,2,…,M,则其离散短时傅里叶变换为:
其中,
DSTFT(i,k)为Lamb波检测波形数据的离散短时傅里叶变换结果;
g(i)为离散短时傅里叶变换的窗函数。
检测数据的时频能量密度分布为:
E(ti,fk)=|DSTFT(i,k)|2。
5.根据权利要求1所述的方法,其特征在于:在时频能量密度分布中,提取中心频率两侧频率点的方法为:
fl=argmin|fk-fc|,fk<fc,
fh=argmin|fk-fc|,fk>fc。
6.根据权利要求1所述的方法,其特征在于:采用线性插值法求取离散的时间点ti对应中心频率处的能量密度E(ti,fc):
7.根据权利要求1所述的方法,其特征在于:通过最小二乘拟合建立基于中心频率的时域能量密度曲线E(t,fc),设拟合函数为时频能量密度为:
8.根据权利要求1所述的方法,其特征在于:提取时域能量密度曲线的峰值时间作为所需模态Lamb波的走时tp:
tp=argmaxE(t,fc)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610018399.5A CN105486759B (zh) | 2016-01-12 | 2016-01-12 | 一种Lamb波时频能量密度析出走时提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610018399.5A CN105486759B (zh) | 2016-01-12 | 2016-01-12 | 一种Lamb波时频能量密度析出走时提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105486759A true CN105486759A (zh) | 2016-04-13 |
CN105486759B CN105486759B (zh) | 2018-04-10 |
Family
ID=55673876
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610018399.5A Active CN105486759B (zh) | 2016-01-12 | 2016-01-12 | 一种Lamb波时频能量密度析出走时提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105486759B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106596724A (zh) * | 2016-12-08 | 2017-04-26 | 清华大学 | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 |
CN109212040A (zh) * | 2018-08-24 | 2019-01-15 | 清华大学 | Lamb波时频脊线分析走时提取方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101762635A (zh) * | 2008-12-25 | 2010-06-30 | 中国石油天然气股份有限公司 | 钢质储罐底板的导波在线检测方法 |
CN101995435B (zh) * | 2010-11-05 | 2012-07-04 | 上海交通大学 | 基于瞬时相位变化度的损伤检测方法 |
CN103239235B (zh) * | 2013-04-25 | 2014-09-10 | 清华大学深圳研究生院 | 一种多功能综合耳声发射检测仪 |
CN104535655B (zh) * | 2014-11-24 | 2017-06-30 | 清华大学 | 一种射线追踪式超声Lamb波缺陷层析成像方法 |
-
2016
- 2016-01-12 CN CN201610018399.5A patent/CN105486759B/zh active Active
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106596724A (zh) * | 2016-12-08 | 2017-04-26 | 清华大学 | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 |
CN106596724B (zh) * | 2016-12-08 | 2019-05-07 | 清华大学 | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 |
CN109212040A (zh) * | 2018-08-24 | 2019-01-15 | 清华大学 | Lamb波时频脊线分析走时提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105486759B (zh) | 2018-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100580445C (zh) | 一种工程结构损伤主动监测中Lamb波时间反转聚焦方法 | |
CN104807888B (zh) | 一种用于微裂纹长度测量的非共线混频超声检测方法 | |
CN104007176B (zh) | 一种复杂岩土工程介质的全波场检测系统及方法 | |
CN104101687B (zh) | 基于超声多普勒与电学多传感器的多相流可视化测试方法 | |
CN103969337B (zh) | 一种基于矢量全聚焦成像的超声阵列裂纹类缺陷方向识别方法 | |
CN102183226B (zh) | 基于多源信息融合的锚杆无损检测方法 | |
CN101571513A (zh) | 用于复合材料层合板质量评定的超声导波检测装置 | |
CN104090034B (zh) | 一种用于导波层析成像的电磁超声兰姆波换能器 | |
CN204154684U (zh) | 桥梁预应力管道注浆密实度检测仪 | |
CN103235039A (zh) | 基于线性调频技术的焊缝缺陷超声tofd检测方法 | |
CN107917957B (zh) | 一种板形结构的损伤检测方法 | |
CN103149274B (zh) | 一种混凝土缺陷检测方法 | |
CN106596724A (zh) | 一种窄带Lamb波检测数据稀疏压缩优化重建方法 | |
CN102589490A (zh) | 一种白车身减薄率超声波检测设备 | |
CN102507739A (zh) | 碳纤维复合材料冲击损伤的超声导波缺陷成像方法 | |
CN102788845A (zh) | 混凝土结构缺陷的巴克编码激励超声检测方法 | |
US20160034422A1 (en) | Technique for Non-Destructive Testing | |
CN106124625B (zh) | 空气耦合超声波高能检测方法及系统 | |
CN104142326A (zh) | 一种衰减系数检测方法 | |
CN103852521B (zh) | 一种用超声波检测多层介质耦合质量的方法 | |
CN105486759A (zh) | 一种Lamb波时频能量密度析出走时提取方法 | |
CN102249367A (zh) | 水处理中超声声强检测与控制方法 | |
Xu et al. | Quantification of guided mode propagation in fractured long bones | |
CN114235962A (zh) | 一种面向各向异性结构的超声导波成像方法及系统 | |
CN104698080A (zh) | 一种利用兰姆波对结构损伤进行状态监测的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |