CN106645401B - 一种基于频率波数估计的损伤定位与重构方法及系统 - Google Patents
一种基于频率波数估计的损伤定位与重构方法及系统 Download PDFInfo
- Publication number
- CN106645401B CN106645401B CN201610911831.3A CN201610911831A CN106645401B CN 106645401 B CN106645401 B CN 106645401B CN 201610911831 A CN201610911831 A CN 201610911831A CN 106645401 B CN106645401 B CN 106645401B
- Authority
- CN
- China
- Prior art keywords
- wave
- wave number
- distribution map
- measured
- frequency
- 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
- 230000006378 damage Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000009826 distribution Methods 0.000 claims abstract description 98
- 235000019687 Lamb Nutrition 0.000 claims abstract description 65
- 230000002159 abnormal effect Effects 0.000 claims abstract description 32
- 239000006185 dispersion Substances 0.000 claims description 14
- 238000005457 optimization Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 241000406668 Loxodonta cyclotis Species 0.000 claims 1
- 230000003902 lesion Effects 0.000 description 10
- 238000001514 detection method Methods 0.000 description 8
- 230000002068 genetic effect Effects 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 3
- 230000035772 mutation Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 239000004411 aluminium Substances 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 210000000349 chromosome Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010429 evolutionary process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
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/041—Analysing solids on the surface of the material, e.g. using Lamb, Rayleigh or shear waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/646—Specific applications or type of materials flaws, defects
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提供一种基于频率波数估计的损伤定位与重构的方法及系统。该方法包括:首先获取待测对象的Lamb波波场图,并依据所述Lamb波波场图,获取待测对象的Lamb波波数分布图,对所述波数分布图进行分析,获取异常波数分布图和正常波数分布图,进而确定待测对象的损伤位置,同时根据所述正常波数分布图,确定正常波数频厚积的关系,又根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。与现有技术相比,该方法及系统不仅可以检测实际正在使用的结构损伤的位置,同时还检测出发生损伤的部分的厚度信息。
Description
技术领域
本发明涉及损伤定位与重构技术领域,特别是涉及一种基于频率波数估计的损伤定位与重构方法及系统。
背景技术
随着工程结构对安全和可靠性的要求日益增高,对材料结构进行定期、有效的无损检测与评估成为结构制造和维护过程中必不可少的环节。在传统的超声检测技术中,如:水浸式超声检测系统,通常利用压电材料制成的超声换能器激励超声波,并通过声耦合介质如水、油等将波传入结构当中,该方法能够用于结构制造、焊接等过程中的损伤检测,但是,鉴于扫描过程中对声耦合介质的需求以及探头与结构间距离的要求,该方法不能用于实际正在使用的结构。
为了解决上述技术问题,发明了具有非接触、对被测构件的尺寸和形状没有特殊要求,能够实现对复杂表面的结构进行检测等特点的激光超声检测技术,其技术通常利用激光脉冲对构件表面进行激励,在构件中激发出超声导波,同时通过传感器对超声导波进行接收,提取结构中的特征参数,定性判断构件中是否存在缺陷,从而实现对构件进行无损检测。
但是,现有的激光超声检测方法虽然能够检测出损伤位置,却无法检测出损伤在厚度方向上的变化参数,从而无法对损伤程度进行评估,如在金属腐蚀的检测时,只能定位腐蚀位置,而无法判断腐蚀程度。
发明内容
本发明的目的是提供一种基于频率波数估计的损伤定位与重构方法及系统,它可以克服现有技术上述的不足,不仅实现损伤定位,同时还检测出发生损伤的部分的厚度信息。
为实现上述目的,本发明提供了如下方案:
本发明提供一种基于频率波数估计的损伤定位与重构方法,该方法包括:
获取待测对象的Lamb波波场图;
根据所述Lamb波波场图,获取待测对象的Lamb波波数分布图;
根据所述波数分布图,获取异常波数分布图和正常波数分布图;
根据所述异常波数分布图,确定待测对象损伤位置;
根据所述正常波数分布图,确定正常波数频厚积的关系;
根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分;
根据所述重构的待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。
可选的,所述获取待测对象的Lamb波波场图,具体包括:
根据待测对象的尺寸,确定扫描范围、扫描点间隔以及采集频率;
根据所述确定的扫描范围、扫描点间隔以及采集频率,扫描待测对象;
采集所述待测对象的各扫描点处的波信号;
根据所述采集的波信号,获取Lamb波波场图。
可选的,所述获取待测对象的Lamb波波数分布图,具体包括:
选取任一频率,对所述获取的Lamb波波场图进行三维短空间傅里叶变换,获取Lamb波波数分布图。
可选的,所述确定正常波数频厚积的关系,具体包括:
提取所述正常波数分布图下不同频率波数信息;
根据所述获取的不同频率下的波数信息,利用最小二乘法拟合出Lamb波频散曲线;
根据所述Lamb波频散曲线,确定正常波数频厚积的关系。
可选的,所述重构待测对象的损伤位置部分,具体包括:
根据所述异常波数分布图和所述正常波数频厚积的关系,利用最优化算法,重构待测对象的损伤位置部分。
本发明还提供一种基于频率波数估计的损伤定位与重构系统,该方系统包括:
波场图获取模块,用于获取待测对象的Lamb波波场图;
分布图获取模块,用于根据所述Lamb波波场图,获取待测对象的Lamb波波数分布图;
异常波数分布图和正常波数分布图获取模块,用于根据所述波数分布图,获取异常波数分布图和正常波数分布图;
位置确定模块,用于根据所述异常波数分布图,确定待测对象损伤位置;
关系确定模块,用于根据所述正常波数分布图,确定正常波数频厚积的关系;
重构模块,用于根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分;
厚度确定模块,用于根据所述重构的待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。
可选的,所述波场图获取模块,具体包括:
确定子模块,用于根据待测对象的尺寸,确定扫描范围、扫描点间隔以及采集频率;
扫描子模块,用于根据所述确定的扫描范围、扫描点间隔以及采集频率,扫描待测对象;
采集子模块,用于采集所述待测对象的各扫描点处的波信号;
波场图获取子模块,用于根据所述采集的波信号,获取Lamb波波场图。
可选的,所述分布图获取模块,具体包括:
分布图获取子模块,用于选取任一频率,对所述获取的Lamb波波场图进行三维短空间傅里叶变换,获取Lamb波波数分布图。
可选的,所述关系确定模块,具体包括:
提取子模块,用于提取所述正常波数分布图下不同频率波数信息;
拟合子模块,用于根据所述获取的不同频率下的波数信息,利用最小二乘法拟合出Lamb波频散曲线;
关系确定子模块,根据所述Lamb波频散曲线,确定正常波数频厚积的关系。
可选的,其特征于,所述重构模块,具体包括:
重构子模块,重构子模块,用于根据所述异常波数分布图和所述正常波数频厚积的关系,利用最优化算法,重构待测对象的损伤位置部分。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供一种基于频率波数估计的损伤定位与重构方法及系统,首先获取待测对象的Lamb波波场图,并依据所述Lamb波波场图,获取待测对象的Lamb波波数分布图,对所述波数分布图进行分析,获取异常波数分布图和正常波数分布图,进而确定待测对象损伤位置;同时根据所述正常波数分布图,确定正常波数频厚积的关系,从而解决了理论求解波数需要已知材料参数的缺点,提高了波数计算的效率,又根据所述异常波数分布图和所述正常波数频厚积的关系,重构损伤位置部分的待测对象,进而确定所述待测对象损伤位置的厚度信息。与现有技术相比,该方法及系统不仅实现了实际正在使用的结构损伤的位置,同时还检测出发生损伤的部分的厚度信息,且该方法及系统实验稳定性好,无需预知材料属性,为航空航天复杂结构的损伤检测与评估提供了新方法及系统。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中的基于频率波数估计损伤定位与重构方法流程图;
图2为本发明实施例中的基于频率波数估计损伤定位与重构系统结构图;
图3为本发明实施例中的基于频率波数估计损伤厚度检测流程图;
图4为本发明实施例中的激光超声成像系统示意图;
图5为本发明实施例中的损伤厚度检测遗传算法运算流程图;
图6为本发明实施例中的损伤厚度检测结果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
参见图1所示,图1为本发明实施例中的基于频率波数估计损伤定位与重构方法流程图,包括:
步骤101:获取待测对象的Lamb波波场图,所述采集的波场信号是宽频信号。
所述获取待测对象的Lamb波波场图,具体包括:
根据待测对象的尺寸,确定扫描范围、扫描点间隔以及采集频率;
根据所述确定的扫描范围、扫描点间隔以及采集频率,扫描待测对象;
采集所述待测对象的各扫描点处的波信号;
根据所述采集的波信号,获取Lamb波波场图。
步骤102:根据所述Lamb波波场图,获取待测对象的Lamb波波数分布图。
所述获取待测对象的Lamb波波数分布图,具体包括:
选取任一频率,对所述获取的Lamb波波场图进行三维短空间傅里叶变换,获取Lamb波波数分布图。
步骤103:根据所述波数分布图,获取异常波数分布图和正常波数分布图,具体包括:依据损伤部位会导致厚度减小,进而使得损伤部位的波数变大,从而确定异常波数分布图和正常波数分布图。
步骤104:根据所述异常波数分布图,确定待测对象损伤位置。
步骤105:根据所述正常波数分布图,确定正常波数频厚积的关系;所述频厚积为波数与频率厚度乘积。
所述确定正常波数频厚积的关系,具体包括:
提取所述正常波数分布图下不同频率波数信息;
根据所述获取的不同频率下的波数信息,利用最小二乘法拟合出Lamb波频散曲线;
根据所述Lamb波频散曲线,确定正常波数频厚积的关系。
步骤106:根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分。
所述重构待测对象的损伤位置部分,具体包括:
根据所述异常波数分布图和所述正常波数频厚积的关系,利用最优化算法,重构待测对象的损伤位置部分。
步骤107:根据所述重构的待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。
本发明提供的上述实施例方法,不仅实现了实际正在使用的结构损伤的位置,同时还检测出发生损伤的部分的厚度信息。
实施例二
参见图2所示,图2为本发明实施例中的基于频率波数估计损伤定位与重构系统结构图,包括:波场图获取模块201、分布图获取模块202、异常波数分布图和正常波数分布图获取模块203、位置确定模块204、关系确定模块205、重构模块206、厚度确定模块207。
所述波场图获取模块201,用于获取待测对象的Lamb波波场图,具体包括:
确定子模块,用于根据待测对象的尺寸,确定扫描范围、扫描点间隔以及采集频率;
扫描子模块,用于根据所述确定的扫描范围、扫描点间隔以及采集频率,扫描待测对象;
采集子模块,用于采集所述待测对象的各扫描点处的波信号;
波场图获取子模块,用于根据所述采集的波信号,获取Lamb波波场图。
所述分布图获取模块202,用于根据所述Lamb波波场图,获取待测对象的Lamb波波数分布图,具体包括:分布图获取子模块,用于选取任一频率,对所述获取的Lamb波波场图进行三维短空间傅里叶变换,获取Lamb波波数分布图。
异常波数分布图和正常波数分布图获取模块203,用于根据所述波数分布图,获取异常波数分布图和正常波数分布图。
所述位置确定模块204,用于根据所述异常波数分布图,确定待测对象损伤位置。
所述关系确定模块205,用于根据所述正常波数分布图,确定正常波数频厚积的关系,具体包括:
提取子模块,用于提取所述正常波数分布图下不同频率波数信息;
拟合子模块,用于根据所述获取的不同频率下的波数信息,利用最小二乘法拟合出Lamb波频散曲线;
关系确定子模块,根据所述Lamb波频散曲线,确定正常波数频厚积的关系。
所述重构模块206,用于根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分,具体包括:重构子模块,用于根据所述异常波数分布图和所述正常波数频厚积的关系,利用最优化算法,重构待测对象的损伤位置部分。
所述厚度确定模块207,用于根据所述重构的待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。
本发明提供的上述实施例系统,不仅实现了实际正在使用的结构损伤的位置,同时还检测出发生损伤的部分的厚度信息。
实施例三
参见图3所示,图3为本发明实施例中的基于频率波数估计损伤厚度检测流程图,包括:
步骤301:波场测量,具体包括:
根据试验件的大小和检测的精度要求,计算出激励所需的空间分辨率以及各激励点在试验件上的位置坐标;
激励单元调整好激光的激励光路后,激励单元发射激光脉冲;
激光脉冲发出的同时,激励单元产生一个触发信号,触发数据采集单元进行数据采集工作,以确保激励与采集的时间同步;
传感器接收到的响应信号经过数据采集单元滤波,放大等处理后被采集卡采集,并保存在控制单元中;
控制单元对信号进行数据进行可视化处理,获得Lamb波传播波场图。
其中,波场测量检测系统示意图包括的硬件由激励单元、数据采集单元和控制单元三部分构成,所述控制单元包括控制系统401以及扫描仪控制器402,所述数据采集单元包括采集设备403、所述激励单元包括激光控制器404以及激光器405。所述激励单元为该系统的核心部分,由Q开关Nd:YAG激光器和振镜式激光反射扫描仪组成。所述数据采集单元由传感器、信号调理设备、数据采集卡组成。控制单元为基于PXI平台的设备。
步骤302:波数测量与损伤位置成像,具体包括:
根据步骤301测量的波场,可以获得了Lamb波的波场W(x,y,t),并在测得的波场上加上二维矩形窗,可以将波场转化为一系列的短空间波场。其中以原点为中心,直径为D的矩形窗函数可以表示式(1)为:
则以(x0,y0)为中心的短空间波场可表示式(2)为:
ws(x,y,t)=w(x,y,t)Ω(x-x0,y-y0,t) (2)
利用三维傅里叶变换,将短空间波场变换到频率波数域,表示式(3)为:
W(kx,ky,ω)=∫∫∫ws(x,y,t)e-j(kxx+kyy-ωt)dxdydt (3)
选取合适的频率ω0,可以获得该频率下的波数分布W(kx,ky,ω0),找到W(kx,ky,ω0)值最大时所对应的波数kx、ky,则位置(x0,y0)处A0模态的波数为重复上述步骤,可以获得任一位置(x,y)处的波数,从而得到波数在检测区域的分布k(x,y),并绘制出波数分布的强度图。由于损伤处的波数会发生变化,因此强度图中波数异常区域即为损伤区域。
步骤303:频散曲线拟合,具体包括:
Lamb波的波数与板的频厚积的关系可由函数k=g(ω·d)来表示。在第二步的波数测量中,我们可以得到健康区域处不同频率ω所对应的波数k,选取20组ω和k的值,利用最小二乘法对其进行多项式拟合,可以获得由频厚积计算波数的公式。
步骤304:损伤厚度重构,具体包括:
损伤深度的重构可以采用多种智能算法,如人工神经网络、遗传算法等,这里使用遗传算法进行重构。首先对需要重构的参数d用二进制整数进行编码,作为遗传算法运算过程中的个体,将步骤302获得Lamb波波数分布图与步骤303中依据频散曲线计算的波数间的误差最小化,得到相应的优化结果,即为重构试验件损伤厚度。
参见4所示,图4为损伤厚度检测的遗传算法运算过程,具体包括:
步骤401:初始化,设置进化代数计数器t=0,设置最大进化代数T。随机生成M个个体作为初始群体P(0)。每个个体即代表一个可能的损伤厚度d。
步骤402:个体评价,计算群体P(t)中各个个体的适应度。在本实施例中,采用理论值与测量值之间的误差的平方和,即式(4):
作为目标函数来进行最小化。
步骤403:选择运算,将选择算子作用于群体。在本实施例中,采用与适应度成正比的概率来确定各个个体复制到下一代群体中的数量。
步骤404:交叉运算,将交叉算子作用于群体。在本实施例中,先对群体进行随机配对,其次随机设置交叉点位置,再相互交换配对染色体之间的部分基因。
步骤405:变异运算,将变异算子作用于群体。在本实施例中,首先随机确定各个个体的基因变异位置,然后依照某一概率将变异点的原有基因值取反。
步骤406:终止条件判断,t是否大于T,若大于T,则以进化过程中所得到的具有最大适应度的个体作为最优解输出,终止计算,执行步骤408;若小于等于T,则执行步骤407。
步骤407:命令t=t+1。
步骤408:厚度输出:根据上述步骤,得到群体中每个具有最大适应度的个体,即得到试验件的损伤厚度。
本发明采用上述公开的实施例,对厚度为2mm的铝板进行扫描,以1mm为间隔,扫描一条200mm的直线,扫描点共计201点,在100-150mm处有一深为1.4mm的凹坑。结果如图6所示,实线为板实际厚度,虚线为采用上述处理方法获得的铝板厚度,可以看出检测结果与实际厚度相符。该实施例采用的方法不仅实现了实际正在使用的结构损伤的位置,同时还检测出发生损伤的部分的厚度信息,且该方法实验稳定性好,无需预知材料属性。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (6)
1.一种基于频率波数估计的损伤定位与重构方法,其特征在于,所述方法包括:
获取待测对象的Lamb波波场图;
根据所述Lamb波波场图,获取待测对象的Lamb波波数分布图;具体包括:选取任一频率,对所述获取的Lamb波波场图进行三维短空间傅里叶变换,获取Lamb波波数分布图;其中,所述Lamb波波数分布图为Lamb波波数分布的强度图;
根据所述波数分布图,获取异常波数分布图和正常波数分布图;
根据所述异常波数分布图,确定待测对象损伤位置;
根据所述正常波数分布图,确定正常波数频厚积的关系;具体包括:提取所述正常波数分布图下不同频率波数信息;根据所述获取的不同频率下的波数信息,利用最小二乘法拟合出Lamb波频散曲线;根据所述Lamb波频散曲线,确定正常波数频厚积的关系;
根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分;
根据所述重构的待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。
2.根据权利要求1所述的方法,其特征在于,所述获取待测对象的Lamb波波场图,具体包括:
根据待测对象的尺寸,确定扫描范围、扫描点间隔以及采集频率;
根据所述确定的扫描范围、扫描点间隔以及采集频率,扫描待测对象;
采集所述待测对象的各扫描点处的波信号;
根据所述采集的波信号,获取Lamb波波场图。
3.根据权利要求1所述的方法,其特征于,所述重构待测对象的损伤位置部分,具体包括:
根据所述异常波数分布图和所述正常波数频厚积的关系,利用最优化算法,重构待测对象的损伤位置部分。
4.一种基于频率波数估计的损伤定位与重构系统,其特征在于,所述系统包括:
波场图获取模块,用于获取待测对象的Lamb波波场图;
分布图获取模块,用于根据所述Lamb波波场图,获取待测对象的Lamb波波数分布图;具体包括:分布图获取子模块,用于选取任一频率,对所述获取的Lamb波波场图进行三维短空间傅里叶变换,获取Lamb波波数分布图;
异常波数分布图和正常波数分布图获取模块,用于根据所述波数分布图,获取异常波数分布图和正常波数分布图;
位置确定模块,用于根据所述异常波数分布图,确定待测对象损伤位置;
关系确定模块,用于根据所述正常波数分布图,确定正常波数频厚积的关系;具体包括:提取子模块,用于提取所述正常波数分布图下不同频率波数信息;拟合子模块,用于根据所述获取的不同频率下的波数信息,利用最小二乘法拟合出Lamb波频散曲线;关系确定子模块,根据所述Lamb波频散曲线,确定正常波数频厚积的关系;
重构模块,用于根据所述异常波数分布图和所述正常波数频厚积的关系,重构待测对象的损伤位置部分;
厚度确定模块,用于根据所述重构的待测对象的损伤位置部分,确定所述待测对象损伤位置的厚度信息。
5.根据权利要求4所述的系统,其特征在于,所述波场图获取模块,具体包括:
确定子模块,用于根据待测对象的尺寸,确定扫描范围、扫描点间隔以及采集频率;
扫描子模块,用于根据所述确定的扫描范围、扫描点间隔以及采集频率,扫描待测对象;
采集子模块,用于采集所述待测对象的各扫描点处的波信号;
波场图获取子模块,用于根据所述采集的波信号,获取Lamb波波场图。
6.根据权利要求4所述的系统,其特征于,所述重构模块,具体包括:
重构子模块,用于根据所述异常波数分布图和所述正常波数频厚积的关系,利用最优化算法,重构待测对象的损伤位置部分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610911831.3A CN106645401B (zh) | 2016-10-19 | 2016-10-19 | 一种基于频率波数估计的损伤定位与重构方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610911831.3A CN106645401B (zh) | 2016-10-19 | 2016-10-19 | 一种基于频率波数估计的损伤定位与重构方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106645401A CN106645401A (zh) | 2017-05-10 |
CN106645401B true CN106645401B (zh) | 2019-07-16 |
Family
ID=58856357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610911831.3A Active CN106645401B (zh) | 2016-10-19 | 2016-10-19 | 一种基于频率波数估计的损伤定位与重构方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106645401B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108801823B (zh) * | 2018-06-25 | 2019-08-23 | 南京航空航天大学 | 一种多尺度的复合材料结构局部疲劳评估方法及系统 |
CN110018235B (zh) * | 2019-03-12 | 2021-07-13 | 上海大学 | 一种基于时域补偿的欠采样Lamb波波数搜索方法 |
CN110376282B (zh) * | 2019-07-15 | 2021-06-08 | 北京航空航天大学 | 一种基于椭圆概率与贝叶斯估计的Lamb波损伤定位法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4274288A (en) * | 1979-07-23 | 1981-06-23 | Rockwell International Corporation | Method for measuring the depth of surface flaws |
CN102043016A (zh) * | 2010-11-05 | 2011-05-04 | 上海交通大学 | 基于兰姆波的自主式损伤识别成像方法 |
CN102128881A (zh) * | 2010-12-22 | 2011-07-20 | 南京邮电大学 | 利用信号分解的Lamb波工程结构损伤监测方法 |
CN102445496A (zh) * | 2011-10-20 | 2012-05-09 | 南京航空航天大学 | 一种基于Lamb波的板状结构无基准快速损伤检测方法 |
CN103134857A (zh) * | 2013-02-21 | 2013-06-05 | 南京邮电大学 | 利用Lamb波反射场的工程结构裂纹损伤监测评估方法 |
CN105372326A (zh) * | 2015-10-20 | 2016-03-02 | 南京航空航天大学 | 一种基于Lamb波波数扫描的空间-波数滤波器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10724994B2 (en) * | 2015-12-15 | 2020-07-28 | University Of South Carolina | Structural health monitoring method and system |
-
2016
- 2016-10-19 CN CN201610911831.3A patent/CN106645401B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4274288A (en) * | 1979-07-23 | 1981-06-23 | Rockwell International Corporation | Method for measuring the depth of surface flaws |
CN102043016A (zh) * | 2010-11-05 | 2011-05-04 | 上海交通大学 | 基于兰姆波的自主式损伤识别成像方法 |
CN102128881A (zh) * | 2010-12-22 | 2011-07-20 | 南京邮电大学 | 利用信号分解的Lamb波工程结构损伤监测方法 |
CN102445496A (zh) * | 2011-10-20 | 2012-05-09 | 南京航空航天大学 | 一种基于Lamb波的板状结构无基准快速损伤检测方法 |
CN103134857A (zh) * | 2013-02-21 | 2013-06-05 | 南京邮电大学 | 利用Lamb波反射场的工程结构裂纹损伤监测评估方法 |
CN105372326A (zh) * | 2015-10-20 | 2016-03-02 | 南京航空航天大学 | 一种基于Lamb波波数扫描的空间-波数滤波器 |
Non-Patent Citations (3)
Title |
---|
双层层压材料中激光超声检测方法的应用;朱倩等;《激光与光电子学进展》;20160303;第53卷;第031402-2至031402-3页"3基于波数分析的损伤识别方法" |
基于Lamb波的复合材料结构损伤成像研究;严刚、周丽;《仪器仪表学报》;20070430;第28卷(第4期);第583-589页 |
结构裂纹损伤的Lamb波层析成像监测与评估研究;王强等;《机械工程学报》;20160331;第52卷(第6期);第30-36页 |
Also Published As
Publication number | Publication date |
---|---|
CN106645401A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110389170B (zh) | 基于Lamb波成像的列车部件裂纹损伤检测方法和系统 | |
Kudela et al. | Identification of cracks in thin-walled structures by means of wavenumber filtering | |
Merazi Meksen et al. | Automatic crack detection and characterization during ultrasonic inspection | |
US8525831B2 (en) | Method and apparatus for three-dimensional visualization and analysis for automatic non-destructive examination of a solid rotor using ultrasonic phased array | |
EP2567190B1 (en) | A method for inspecting an object by means of ultrasound | |
CN106645401B (zh) | 一种基于频率波数估计的损伤定位与重构方法及系统 | |
CN107014668A (zh) | 一种基于压电和智能涂层传感器的疲劳裂纹综合监测方法 | |
Zhang et al. | Monte Carlo inversion of ultrasonic array data to map anisotropic weld properties | |
CN109900804A (zh) | 一种基于超声导波的金属材料裂纹定量监测方法 | |
KR101774514B1 (ko) | 불규칙들을 측정할 때 saft 분석을 개선하기 위한 방법 및 디바이스 | |
CN104422732A (zh) | 超声波探伤传感器以及超声波探伤方法 | |
Tola et al. | Bolt looseness detection based on ultrasonic wavefield energy analysis using an Nd: YAG pulsed laser scanning system | |
CN106501285B (zh) | 无损检测预应力管道的压浆密实度的设备及检测方法 | |
CN109900789A (zh) | 一种树木纵截面内部缺陷成像方法 | |
CN104040329B (zh) | 用于检测在检查对象内部的缺陷的方法和装置 | |
CN103837580B (zh) | 一种基于超声和电磁超声相结合的双模无损检测方法 | |
Anderson et al. | Assessment of crack detection in heavy-walled cast stainless steel piping welds using advanced low-frequency ultrasonic methods | |
US11402356B2 (en) | Object identification system and method | |
CA3040331A1 (en) | Method and device for detecting and characterizing a reflecting element in an object | |
Merazi Meksen et al. | A method to improve and automate flat defect detection during ultrasonic inspection | |
Liu et al. | Sequential dynamic aperture focusing strategy for transmissive ultrasonic phase array tomography | |
Dziewierz et al. | An application-specific design approach for 2D ultrasonic arrays | |
CN106198745B (zh) | 基于反射能量比参量的锚固缺陷长度识别方法 | |
JP2017500553A (ja) | 断片の表面を再構築する方法 | |
Kim et al. | Improvement of wood CT images by consideration of the skewing of ultrasound caused by growth ring angle |
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 |