CN104535655B - 一种射线追踪式超声Lamb波缺陷层析成像方法 - Google Patents
一种射线追踪式超声Lamb波缺陷层析成像方法 Download PDFInfo
- Publication number
- CN104535655B CN104535655B CN201410683016.7A CN201410683016A CN104535655B CN 104535655 B CN104535655 B CN 104535655B CN 201410683016 A CN201410683016 A CN 201410683016A CN 104535655 B CN104535655 B CN 104535655B
- Authority
- CN
- China
- Prior art keywords
- ray
- emat
- extrapolation
- defect
- imaging
- 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
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/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
-
- 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/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/0672—Imaging by acoustic tomography
-
- 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
- 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/22—Details, e.g. general constructional or apparatus details
- G01N29/24—Probes
- G01N29/2412—Probes using the magnetostrictive properties of the material to be examined, e.g. electromagnetic acoustic transducers [EMAT]
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Electromagnetism (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明属无损检测领域,涉及一种射线追踪式超声Lamb波缺陷层析成像方法;先将待测材料成像区域划分为多个小网格,并在此区域两侧设置相同数量的发射EMAT和接收EMAT。然后令所有发射端激发A0模态Lamb波,对侧所有接收端同时接收。对全部检测波形的分析结果提取并记录相应走时,再使用联合迭代重建算法确定每个网格的慢度及缺陷分布,将得到的缺陷分布作为已知条件,建立RT外推公式,利用外推法逐步对射线路径进行追踪并加以修正,如此反复多次,直至通过计算获得较高的成像精度。本发明避免了基于直射线理论的传统TOF跨孔层析成像方法成像精度低的问题,计算准确、高效、快速,为提高缺陷成像精度提供重要依据,具有广阔的应用前景。
Description
技术领域
本发明涉及一种射线追踪式超声Lamb波缺陷层析成像方法,属于无损检测技术领域。
背景技术
利用超声Lamb波缺陷层析成像技术可以快速、有效地获得缺陷的轮廓和尺寸等具体信息。该技术继承了传统超声Lamb波检测的诸多优点。Lamb波走时(Time-of-flight,TOF,指Lamb波在收发换能器对之间的传播时间)跨孔层析成像技术,利用线性换能器阵列进行扇束投影和成像,是一种基于迭代法的高效、快速的层析成像方法。传统的Lamb波TOF跨孔层析成像均是基于直射线理论,即假设Lamb波射线在传播过程中遇到缺陷后仍然沿直线传播。但是,当平板中存在强散射缺陷时,Lamb波射线传播方向会发生偏转,即Lamb波射线实际将沿弯曲路径传播。此时,已有的直射理论将不再适用,并将导致层析成像的精度下降。上述问题是制约Lamb波走时跨孔层析成像技术发展的一个瓶颈问题。中国专利文献公开了一种“刚制储罐底板的导波在线检测方法”该技术涉及一种钢质储罐底板的导波在线检测方法,采用超声探头阵列收发导波,通过滤波反投影算法重建层析图像,可快速实现对储罐底板的全面检测和缺陷分析,但该技术局限于矩阵变换理论,只能判断出缺陷的存在位置及严重程度,而对其具体尺寸和轮廓并不能达到十分精确的估计。
发明内容
本发明的目的是提供一种射线追踪式超声Lamb波缺陷层析成像方法,该方法按如下步骤进行:
步骤1:在待测材料上选定成像区域,将所述成像区域划分成N1×N2个小网格;在成像区域的一侧设置M个发射EMAT(电磁声换能器),对侧相同位置处设置M个接收EMAT;其中,N1、N2、M为正整数;
步骤2:激励第1~M个发射EMAT全向激发A0模态Lamb波,M个接收EMAT依次接收Lamb波;
步骤3:利用伪Wigner-Ville分布(PWVD)对M×M个EMAT检测波形进行时频分析和模态识别;
步骤4:提取M×M个检测波形的时频分析结果,记录M×M个A0模态Lamb波的走时T1~TM×M;
步骤4:提取M×M个检测波形的时频分析结果,记录M×M个A0模态Lamb波的走时T1~TM×M;
步骤5:使用联合迭代重建算法确定每个小网格的慢度及缺陷分布;
步骤6:在成像区域中建立三维直角坐标系,然后建立初始发射角度,利用外推法逐步对经过缺陷区域的全部射线路径进行追踪,对于射线i前进方向上的一点P(x,y),得到Lamb波的相速度值cp和和其中,x、y为P点沿x方向、y方向的坐标值
步骤7:将得到的cp、和代入射线追踪外推公式,找到新的外推点;
步骤8:判断新找到的外推点坐标值是否到达或超出子网格区域边界的坐标值,若否,则返回步骤6;若是,则结束对射线i的射线追踪,记录此时外推点的坐标作为终点坐标;
步骤9:判断终点坐标是否到达或接近接收EMAT,若是,则得到射线i的预期最短走时路径;若否,则改变初始发射角的值,返回步骤6,直至找到连接收发EMAT的射线i的最短走时路径;
步骤10:重复步骤6~9,直至找到经过缺陷区域的全部射线的最短走时路径。
步骤11:将得到的全部射线的最短走时路径作为已知量代入步骤5联合迭代重建算法,再次得到网格慢度,以及慢度曲线在缺陷区域的突变情况,判断缺陷大小及分布。
步骤12:判断缺陷分布的成像精度是否满足要求,若是,则得到预期结果,若否,则重复步骤2~11。
收发EMAT直径为20~80mm,相邻EMAT中心间距为20~100mm。
射频功率放大器激励发射EMAT,射频功率放大器的激发频率为50~500kHz。
激发的Lamb波为单一的A0模态Lamb波。
联合迭代重建算法确定每个网格的慢度及缺陷分布的方法为:
其中,Sj为待求的第j个网格的慢度;Lij为第i条投影射线在第j个网格中的长度;Ti为第i条投影射线的实测走时;n为正整数,且n=N1×N2;m为正整数,且m=M×M。
初始发射角度的范围是0~180度。
利用外推法逐步对射线路径进行追踪的方法为:以外推点为中心获取16个像素;得到相应的16个cp值;将cp值代入二维三次样条插值公式,对Lamb波cp的分布情况进行二维曲面拟合,求解每一步外推所需的和求解公式为:
其中:cp是相速度,x,y分别表示16个像素中心点沿x轴方向、y轴方向的坐标值;Ckl表示与16个像素相对应的待求解系数;k,l=0,1,2,3。
根据射线追踪外推公式,寻找新的外推点的方法为:
其中,yk表示经过k步外推获得的射线i上的离散点沿y轴方向的坐标值;yk+1、yk-1分别表示上述离散点前、后两点沿y轴方向的坐标值;Δx表示射线外推的固定步长,且
改变初始发射角的值的方法为:每次增加或减少2~6度。
本发明的有益效果是:本发明将RT迭代算法模型与TOF跨孔层析成像技术相结合,综合两者优点,通过正反问题的交替反复求解,对正问题中的射线路径加以修正,以还原缺陷的真实分布。本发明计算准确、高效、快速,解决了基于直射线理论的传统层析成像方法成像精度低的问题,具有广阔的应用前景;同时避免了基于直射线理论的传统TOF跨孔层析成像方法成像精度低的问题,计算准确、高效、快速,为提高缺陷成像精度提供重要依据。
附图说明
图1为射线追踪式超声Lamb波层析成像流程图。
图2为成像区域网格划分示意图。
具体实施方式
本发明提出的方法,其实质是建立一种全新的射线追踪(Ray Tracing,RT)迭代算法模型,并将其与TOF跨孔层析成像技术相结合,综合两者优点,通过正反问题的交替反复求解,对正问题中的射线路径加以修正,以获得更为真实的缺陷分布情况。下面结合实施例对本发明做进一步说明:
步骤1:取厚度为1~5mm的待测铝板,在上表面选择一块面积为A×B的矩形区域作为成像区域进行实验,其中,A、B为正实数,且均大于30mm。在所述区域内加工一个直径大于等于30mm、深为1~5mm的标准人造腐蚀坑缺陷。然后将所述成像区域划分成N1×N2个小网格(即N1×N2个像素)。在成像区域的一侧等间距设置M个发射EMAT(电磁声换能器),对侧相同位置处设置相同数量的接收EMAT。所述的EMAT直径为20~80mm,同侧相邻EMAT中心间距为20~100mm。
步骤2:选择一台射频功率放大器用于激励EMAT,令发射端中所有M个发射EMAT均被全向激发A0模态Lamb波,所有M个接收EMAT进行接收,激发频率为50~500kHz。所述的Lamb波为单一的A0模态Lamb波,不含有A0以外其他模态Lamb波。
步骤3:利用伪Wigner-Ville分布(PWVD)对M×M个EMAT检测波形进行时频分析和模态识别;
步骤4:提取所有M×M个检测波形的时频分析结果,并记录M×M个A0模态Lamb波的走时,记为T1~TM×M。
步骤5:使用联合迭代重建算法(SIRT)公式确定每个小网格的慢度(速度的倒数)及缺陷的分布:
其中,Sj为待求的第j个小网格的慢度;Lij为第i条投影射线在第j个小网格中的长度;Ti为第i条投影射线的实测走时,n=N1×N2,m=M×M;
Lij通过几何关系计算得出:已知成像区域面积,小网格边长,收发EMAT数量、位置和间距;按照直射线理论计算得出射线路径,可知射线从发射端到接收端经过的每一个小网格中的长度。
步骤6:首先在成像区域建立一个三维直角坐标系,其中,x、y代表射线上的一点在二维平面上的位置,z代表该点的相速度cp,然后进行子网格区域划分,每个子网格区域包含4×4=16个像素,接着建立一个初始发射角度,范围在0~180度之间;
利用外推法逐步对经过缺陷区域的全部射线路径进行追踪;
对于第i条射线,利用外推法逐步对其射线路径进行追踪。根据Lamb波的传播特性可知,当遇到腐蚀缺陷时,其相速度cp会随板厚的变化而变化。因此,在外推点区域,以外推点为中心,获取4×4区域的16个像素,得到对应16个A0模态Lamb波的相速度cp值,代入二维三次样条插值(2D-CSI)公式:
构成16个线性方程,对Lamb波cp的分布情况进行二维曲面拟合,进而精确求解每一步外推所需的和
其中,cp是相速度,x,y表示沿追踪方向x、y的坐标值;Ckl(k,l=0,1,2,3)表示16个待求解系数,16个像素中每一个像素都对应一个x、y值。
步骤7:将步骤6得到的cp、和代入射线追踪(Ray Tracing,RT)外推公式,完成一步外推:
其中,yk表示经过k步外推获得的射线i上的离散点沿y轴方向的坐标值;
yk+1、yk-1分别表示上述离散点前、后两点沿y轴方向的坐标值;
Δx表示射线外推的固定步长,具体有:
步骤8:判断新找到的外推点坐标值是否到达或超出子网格区域边界坐标值,若不是,则返回步骤5;若是,则结束对射线i的射线追踪,记录下终点坐标(xf,yf)。
步骤9:判断终点坐标(xf,yf)是否到达或接近目标接受位置,若是,则得到射线i的预期最短走时路径;若不是,则将初始发射角每次增加或减少2~6度,再重复步骤6~8,直到找到连接收发EMAT的射线i的最短走时路径。
步骤10:重复步骤6~9,直至找到经过缺陷区域的全部射线的最短走时路径。
步骤11:将得到的全部射线的最短走时路径作为已知量代入步骤5所述的联合迭代重建算法(SIRT)公式,再次对网格慢度进行求解,得到所有射线在各个小网格内的慢度,然后根据慢度曲线在缺陷区域的突变情况即可判断缺陷大小及分布。如此,如果成像精度不高,可重复步骤2~11,获得更高的成像精度。
可变例1
步骤1:待测铝板厚度为3mm,成像区域的面积为边长为640mm的正方形区域;标准人造腐蚀坑缺陷的直径为30mm、深为2mm;缺陷中心点位于(320mm,460mm)。成像区域被划分为128×128个小网格(即128×128个像素)。在成像区域的一侧等间距的设置14个发射EMAT位置,对侧相同位置处设置相同数量的接收EMAT位置,所述的EMAT直径为30mm,相邻EMAT中心间距为45.7mm。
步骤2:射频功率放大器选择AG1024,用于激励EMAT,令发射端中第一个发射EMAT全向激发A0模态Lamb波,所有14个接收EMAT依次进行接收,激发频率为290kHz。重复操作,令发射端其他13个发射EMAT均激发过一次A0模态Lamb波。所述的Lamb波为单一的A0模态Lamb波,不含A0以外其他模态Lamb波。
步骤3:利用伪Wigner-Ville分布(PWVD)对M×M个EMAT检测波形进行时频分析和模态识别;
步骤4:提取所有196个检测波形的分析结果并记录A0模态Lamb波的走时,记为T1~T196。
步骤5:使用联合迭代重建算法(SIRT)确定每个网格的慢度(速度的倒数)及缺陷分布:
其中,Sj为待求的第j个网格的慢度;
Lij为第i条投影射线在第j个网格中的长度;
Ti为第i条投影射线的实测走时;
可知,n=128×128=16384,m=14×14=196。
计算出的缺陷直径值约为40mm。
步骤6:在成像区域建立一个三维直角坐标系,然后进行子网格区域划分,每个子网格区域包含4×4=16个像素;建立一个30度的初始发射角,利用外推法逐步对射线路径进行追踪。在外推点区域均匀获取16个像素,得到对应16个cp值,代入二维三次样条插值(2D-CSI)公式,构成16个线性方程,对Lamb波cp的分布情况进行二维曲面拟合,精确求解每一步外推所需的和
其中,cp是相速度,x,y表示追踪方向;Ckl(k,l=0,1,2,3)表示16个待求解系数,16个像素中每一个像素都对应一个x、y值。
步骤7:将步骤6得到的cp、和代入RT外推公式,完成一步外推:
其中,yk表示经过k步外推获得的射线i上的离散点沿y轴方向的坐标值;
yk+1、yk-1分别表示上述离散点前、后两点沿y轴方向的坐标值;
Δx表示射线外推的固定步长,具体有:
步骤8:判断新找到的外推点坐标值是否到达或超出子网格区域边界坐标值,若不是,则返回步骤5;若是,则结束对射线i的RT,记录下终点坐标(xf,yf)。
步骤9:判断终点坐标(xf,yf)是否到达或接近目标接受位置(即第i对EMAT中的接收EMAT),若是,则得到射线i的预期最短走时路径;若不是,则每次增加或减少3度,从而改变初始发射角的值,重复步骤6~8,直到找到连接收发EMAT的射线i的最短走时路径。
步骤10:重复步骤6~9,直至找到经过缺陷区域的全部射线的最短走时路径。
步骤11:将得到的全部射线的最短走时路径作为已知量代入步骤5所述的联合迭代重建算法(SIRT)公式,再次对网格慢度进行求解,得到所有射线在各个小网格内的慢度,然后根据慢度曲线在缺陷区域的突变情况即可判断缺陷大小及分布。
如此,重复步骤2~11,以获得更高的成像精度;通常的,缺陷长,宽,深的量化误差小于10%时认为成像精度满足要求。本实施例最后得到的缺陷直径值约为32mm,与实际值吻合良好,量化误差小于10%。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (10)
1.一种射线追踪式超声Lamb波缺陷层析成像方法,包括
步骤1:在待测材料上选定成像区域,将所述成像区域划分成N1×N2个小网格;在成像区域的一侧设置M个发射EMAT,对侧相同位置处设置M个接收EMAT;其中,N1、N2、M为正整数;
步骤2:激励第1~M个发射EMAT全向激发A0模态Lamb波,M个接收EMAT依次接收Lamb波;
其特征在于,所述方法还包括
步骤3:利用伪Wigner-Ville分布对M×M个EMAT检测波形进行时频分析和模态识别;
步骤4:提取M×M个检测波形的时频分析结果,记录M×M个A0模态Lamb波的走时T1~TM×M;
步骤5:使用联合迭代重建算法确定每个小网格的慢度及缺陷分布;
步骤6:在成像区域中建立三维直角坐标系,建立初始发射角度,利用外推法逐步对经过缺陷区域的全部射线路径进行追踪,对于射线i前进方向上的一点P(x,y),得到Lamb波的相速度值cp和和其中,x、y为P点沿x方向、y方向的坐标值;
步骤7:将得到的cp、和代入射线追踪外推公式,找到新的外推点;
步骤8:判断新找到的外推点坐标值是否到达或超出子网格区域边界坐标值,若否,则返回步骤6;若是,则结束对射线i的射线追踪,记录此时外推点的坐标作为终点坐标;
步骤9:判断终点坐标是否到达或接近接收EMAT,若是,则得到射线i的预期最短走时路径;若否,则改变初始发射角的值,返回步骤6,直至找到连接收发EMAT的射线i的最短走时路径;
步骤10:重复步骤6~9,直至找到经过缺陷区域的全部射线的最短走时路径;
步骤11:将得到的全部射线的最短走时路径作为已知量代入步骤5联合迭代重建算法,再次得到网格慢度,以及慢度曲线在缺陷区域的突变情况,判断缺陷大小及分布。
2.根据权利要求1所述的方法,其特征在于:收发EMAT直径为20~80mm,相邻EMAT中心间距为20~100mm。
3.根据权利要求1所述的方法,其特征在于:射频功率放大器激励发射EMAT,射频功率放大器的激发频率为50~500kHz。
4.根据权利要求1所述的方法,其特征在于:激发的Lamb波为单一的A0模态Lamb波。
5.根据权利要求1所述的方法,其特征在于:所述联合迭代重建算法确定每个网格的慢度及缺陷分布的方法为:
其中,
Sj为待求的第j个网格的慢度;
Lij为第i条投影射线在第j个网格中的长度;
Ti为第i条投影射线的实测走时;
n为正整数,且n=N1×N2;
m为正整数,且m=M×M。
6.根据权利要求1所述的方法,其特征在于:所述初始发射角度的范围是0~180度。
7.根据权利要求1所述的方法,其特征在于:利用外推法逐步对射线路径进行追踪的方法为:
以外推点为中心获取16个像素;
得到相应的16个cp值;
将cp值代入二维三次样条插值公式,对Lamb波cp的分布情况进行二维曲面拟合,求解每一步外推所需的和
求解公式为:
其中:cp是相速度,
x,y分别表示16个像素中心点沿x轴方向、y轴方向的坐标值;
Ckl表示与16个像素相对应的待求解系数;
k,l=0,1,2,3。
8.根据权利要求1所述的方法,其特征在于:根据射线追踪外推公式,寻找新的外推点的方法为:
其中,
yk表示经过k步外推获得的射线i上的离散点沿y轴方向的坐标值;
yk+1、yk-1分别表示上述离散点前、后两点沿y轴方向的坐标值;
Δx表示射线外推的固定步长,且
9.根据权利要求1所述的方法,其特征在于:改变初始发射角的值的方法为:增加或减少2~6度。
10.根据权利要求1所述的方法,其特征在于:所述方法还包括,
步骤12:判断缺陷分布的成像精度是否满足要求,若是,则得到预期结果,若否,则重复步骤2~11。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410683016.7A CN104535655B (zh) | 2014-11-24 | 2014-11-24 | 一种射线追踪式超声Lamb波缺陷层析成像方法 |
US14/948,398 US10197534B2 (en) | 2014-11-24 | 2015-11-23 | Method and device for testing defect based on ultrasonic lamb wave tomography |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410683016.7A CN104535655B (zh) | 2014-11-24 | 2014-11-24 | 一种射线追踪式超声Lamb波缺陷层析成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104535655A CN104535655A (zh) | 2015-04-22 |
CN104535655B true CN104535655B (zh) | 2017-06-30 |
Family
ID=52851225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410683016.7A Active CN104535655B (zh) | 2014-11-24 | 2014-11-24 | 一种射线追踪式超声Lamb波缺陷层析成像方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10197534B2 (zh) |
CN (1) | CN104535655B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044212A (zh) * | 2015-06-30 | 2015-11-11 | 清华大学 | 一种多模态超声Lamb波复杂缺陷层析成像方法 |
CN105136906A (zh) * | 2015-08-28 | 2015-12-09 | 电子科技大学 | 复合材料单向板超声层析成像减少伪影的方法 |
US20170191966A1 (en) * | 2016-01-04 | 2017-07-06 | General Electric Company | Distributed circle method for guided wave based corrosion detection in plate-like structures |
CN105486759B (zh) * | 2016-01-12 | 2018-04-10 | 清华大学 | 一种Lamb波时频能量密度析出走时提取方法 |
CN106198725B (zh) * | 2016-06-22 | 2019-04-30 | 江苏大学 | 一种基于特征导波的对接焊缝缺陷检测系统及检测方法 |
EP3510564B1 (en) * | 2016-09-06 | 2022-11-02 | Virtamed AG | Ray-tracing methods for realistic interactive ultrasound simulation |
JP6882729B2 (ja) * | 2016-10-05 | 2021-06-02 | 株式会社Core技術研究所 | 内部品質評価方法、及び内部品質解析システム |
CN106525975B (zh) * | 2016-11-02 | 2019-03-01 | 清华大学 | 金属板实际复杂缺陷磁声阵列导波散射成像方法 |
CN107045003A (zh) * | 2016-12-20 | 2017-08-15 | 北京航空航天大学 | 一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法 |
CN108318584B (zh) * | 2018-01-25 | 2023-09-15 | 上海市民防地基勘察院有限公司 | 单孔声波联合跨孔ct检测桩身质量和倾斜度的方法 |
CN109900789B (zh) * | 2019-03-22 | 2020-05-08 | 江南大学 | 一种树木纵截面内部缺陷成像方法 |
CN111624251B (zh) * | 2020-03-18 | 2023-01-31 | 北京工业大学 | 一种应用于大型板壳结构的智能兰姆波缺陷定位方法 |
CN111398427A (zh) * | 2020-04-03 | 2020-07-10 | 中瓴埃斯科(重庆)环保产业有限公司 | 大型储罐底板成像方法 |
CN114371217A (zh) * | 2021-11-27 | 2022-04-19 | 北京工业大学 | 一种基于密集阵列的智能兰姆波缺陷定位方法 |
CN114324591B (zh) * | 2021-12-31 | 2023-09-01 | 上海理工大学 | 基于Radon变换的Lamb波波包飞行时间检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353718A (zh) * | 2011-07-11 | 2012-02-15 | 南京航空航天大学 | 用于复合材料板结构损伤监测的Lamb波损伤概率成像方法 |
WO2012026054A1 (ja) * | 2010-08-27 | 2012-03-01 | 株式会社日立製作所 | 内部欠陥検査方法及びその装置 |
CN102507739A (zh) * | 2011-11-15 | 2012-06-20 | 北京理工大学 | 碳纤维复合材料冲击损伤的超声导波缺陷成像方法 |
CN104132998A (zh) * | 2014-08-06 | 2014-11-05 | 北京科技大学 | 一种基于超声扫描显微镜的内部微观缺陷检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6128092A (en) * | 1999-07-13 | 2000-10-03 | National Research Council Of Canada | Method and system for high resolution ultrasonic imaging of small defects or anomalies. |
KR100784072B1 (ko) * | 2003-09-22 | 2007-12-10 | 김형윤 | 구조물의 건전성 감시용 센서 및 시스템 |
US7783433B2 (en) * | 2007-05-22 | 2010-08-24 | Honeywell International Inc. | Automated defect detection of corrosion or cracks using SAFT processed Lamb wave images |
JP5253424B2 (ja) * | 2010-01-05 | 2013-07-31 | 日立Geニュークリア・エナジー株式会社 | 超音波探傷方法及び超音波探傷装置 |
US10025272B2 (en) * | 2013-01-25 | 2018-07-17 | General Electric Company | Ultrasonic holography imaging system and method |
US20170191966A1 (en) * | 2016-01-04 | 2017-07-06 | General Electric Company | Distributed circle method for guided wave based corrosion detection in plate-like structures |
-
2014
- 2014-11-24 CN CN201410683016.7A patent/CN104535655B/zh active Active
-
2015
- 2015-11-23 US US14/948,398 patent/US10197534B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012026054A1 (ja) * | 2010-08-27 | 2012-03-01 | 株式会社日立製作所 | 内部欠陥検査方法及びその装置 |
CN102353718A (zh) * | 2011-07-11 | 2012-02-15 | 南京航空航天大学 | 用于复合材料板结构损伤监测的Lamb波损伤概率成像方法 |
CN102507739A (zh) * | 2011-11-15 | 2012-06-20 | 北京理工大学 | 碳纤维复合材料冲击损伤的超声导波缺陷成像方法 |
CN104132998A (zh) * | 2014-08-06 | 2014-11-05 | 北京科技大学 | 一种基于超声扫描显微镜的内部微观缺陷检测方法 |
Non-Patent Citations (2)
Title |
---|
Ultrasonic Lamb wave tomography using scanned EMATs and wavelet processing;K.S.HO et al.;《Nondestructive Testing and Evaluation》;20070331;第22卷(第1期);第19-34页 * |
各向同性薄板中横穿孔缺陷的超声兰姆波层析成像;张海燕等;《声学学报》;20070131;第32卷(第1期);第83-90页 * |
Also Published As
Publication number | Publication date |
---|---|
US10197534B2 (en) | 2019-02-05 |
CN104535655A (zh) | 2015-04-22 |
US20160146762A1 (en) | 2016-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104535655B (zh) | 一种射线追踪式超声Lamb波缺陷层析成像方法 | |
Malyarenko et al. | Ultrasonic Lamb wave diffraction tomography | |
Flynn et al. | Maximum-likelihood estimation of damage location in guided-wave structural health monitoring | |
CN106772641B (zh) | 估测平均土壤含水量的方法 | |
CN106124632A (zh) | 一种基于超声波的混凝土密实性评估方法 | |
Zhang et al. | Monte Carlo inversion of ultrasonic array data to map anisotropic weld properties | |
MXPA05006835A (es) | Metodos para determinar los parametros de yacimiento y pozo de sondeo utilizando tomografia de volumen de fresnel. | |
CN109696480B (zh) | 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法 | |
CN105488795A (zh) | 一种复合材料损伤识别方法 | |
CN110261485A (zh) | 一种超声波测量材料内部各处弹性模量及泊松比的方法 | |
CN100464185C (zh) | 混凝土超声层析成像算法 | |
CN107389793A (zh) | 一种铝合金多孔结构孔边腐蚀损伤监测方法 | |
CN105403878A (zh) | 一种基于时延敏感核的海洋声层析方法 | |
CN105212961A (zh) | 一种声辐射剪切波波速检测方法及系统 | |
He et al. | A fast acoustic emission beamforming localization method based on Hilbert curve | |
CN106645420B (zh) | 基于费马原理的棒材超声线型阵列成像检测方法 | |
CN103018333A (zh) | 分层物体的合成孔径聚焦超声成像方法 | |
CN105044212A (zh) | 一种多模态超声Lamb波复杂缺陷层析成像方法 | |
WO2020191896A1 (zh) | 一种树木纵截面内部缺陷成像方法 | |
CN113468804B (zh) | 一种基于矩阵束和深度神经网络的地下管道识别方法 | |
CN1967177A (zh) | 一种在隧道中测取掌子面前方岩体波速的方法 | |
CN102253122B (zh) | 一种基于柔性超声阵列换能器的多声束自动扫描成像方法 | |
CN115452670A (zh) | 一种流动水体中泥沙含量的测定方法及装置 | |
CN115754007A (zh) | 一种基于声发射技术和层析成像技术的损伤检测方法 | |
CN113433214B (zh) | 一种复合板Lamb波时间反转多路径成像损伤定位方法 |
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 |