CN105510444B - 基于超声导波的双向时间反演损伤成像方法 - Google Patents
基于超声导波的双向时间反演损伤成像方法 Download PDFInfo
- Publication number
- CN105510444B CN105510444B CN201510852380.6A CN201510852380A CN105510444B CN 105510444 B CN105510444 B CN 105510444B CN 201510852380 A CN201510852380 A CN 201510852380A CN 105510444 B CN105510444 B CN 105510444B
- Authority
- CN
- China
- Prior art keywords
- signal
- wave
- damage
- time reversal
- guided wave
- 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
Abstract
本发明公开了一种基于超声导波的双向时间反演损伤成像方法,包括:采用调制的激励信号对损伤实物进行主动应力波探测,在被测区域的投射端与反射端两端阵列上同步采集信号;对采集到的信号进行频散补偿分离,分离出纵向模态导波和弯曲模态导波;提取分离后的信号中的纵向模态导波,做时间反演处理,放大反演后信号;将反演后信号导入有限元模型中进行仿真反演,提取模型中所有点云所有时刻的位移值,进行时间、空间聚焦处理,得到振动云图;定义空间曲线将振动云图转为三维彩色云图,以三维彩色点云形式进行损伤成像。该方法从整体上分析损伤的特征,创新性地提取不同模态导波进行透射、反射两端同时做时间反演分析,具有重要的学术价值和实际意义。
Description
技术领域
本发明涉及无损检测技术领域,尤其涉及一种基于超声导波的细长状结构的双向时间反演损伤成像方法。
背景技术
在交通运输、建筑桥梁、输电线缆等各领域中,往往使用如绞线、管道、钢丝绳等各类复杂长状类结构,这些复杂结构在工作中会长期受到外部环境的影响而在产生各种损伤,如磨损、锈蚀、断股等,这些损伤造成极大的隐患,威胁到人员及设备财产的安全。因此对这类复杂长状结构的检测意义重大,也出现许多检测方法。
目前按是否采用应力波进行损伤检测,分为非应力波损伤检测方法和应力波损伤检测方法两大类。非应力波检测方法主要有光学法、电流法、电磁法、X射线法和电涡流法等,能够很好对被测物裸露部分进行缺陷检测,然而此类方法往往需采用与被测物间相对运动方式来实现,细长状结构的结点结构覆盖区域成为检测盲区,应用受到很大限制。应力波检测方法主要有声发射法、超声波反射法、穿透式导波法。声发射法可采用感应器测量承力索绞线结构变化时产生弹性波从而可分析损伤情况,但属于被动检测方式,受限于只有在结构产生损伤时才能起到作用,是瞬间响应,然而为后期超声波反射法、穿透式导波法研究奠定了基础。超声波反射法采用换能器在材料表面激发询问脉冲应力波,该应力波在内部与材料的微结构相互作用后到达同边接收换能器,然后对接收到的波形信号进行分析,提取出其中能反映结构损伤因子。该法不但能对承力索绞线进行损伤检测和特性分析,而且补充了声发射被动检测方式,可以用于动态承力索绞线缺陷的主动监测,然而由于容易受到传播反射过程的介质损耗等影响,在大范围结点结构存在下的监测对象应用中受到了一定限制。穿透式导波法主要采用压电换能器在材料表面激发应力波,该应力波在内部与材料的微结构(材料基体和损伤区域)相互作用后继续传播,到达置于大型结构另一端的一个或多个接收透射波信号换能器。该法由于激励信号可设可控、传播距离远、接收信号覆盖面广等,在结构健康监测领域越来越受重视,特别适用于杆、管等细长状类特殊结构缺陷检测。因此,面对复杂结构损伤检测中出现的问题以及超声导波的特殊传输特性,将透射法与反射法结合起来,取两者优点是值得探索的一条新途径,并具有重要的学术价值和实际意义。
时间反演是指对采集到的信号,在时域上将其反转的一种逆序处理方法。在频域上,这种处理等效于相位共轭。时间反演具有空间聚焦和时间聚焦两方面特性。时间反演的时间聚焦,指各路时间反演信号的最大能量在复杂媒质中经过多径传输后,会在同一时刻到达目标接收点。时间反演的空间聚焦,指在没有任何关于目标接收点先验知识的情况下,时间反演信号会自适应地聚焦到目标接收点所在位置。时间反演技术利用时间反演信号的聚焦特性有效降低了频散效应和非均匀介质的影响。目前人们已经利用时间反演聚焦信号的峰值、能量、时空聚焦本特征,应用时间反演技术处理各种损伤测方法检测到的损伤信号后对多种结构进行了损伤检测。
然而,传统时间反演方法,要求信号从损伤到接收阵列单元之间的路径较短且不尽相同,且所需的接收阵列单元数量尽可能多以弥补传播损耗,损伤目标点才可能成为空间最大信号所在的点被检测出来。而对于细长状实心结构,导波多以平面波的传播,同侧径向单元间的路径差异小,传统时间反演方法难以实现。
本发明将传统的超声导波透射法、反射法与时间反演法相结合,有效补偿非均匀复杂环境、媒质引起的信号衰减,及细长状结构传感阵列聚焦困难等问题,且具备一定的抗干扰和抗噪声性能等优点。故此,在初步识别的基础上,实现细长型结构的高精度和高准确率的损伤检测,也是各类复杂长状类结构损伤检测方法今后非常重要的研究方向。
发明内容
为解决上述技术问题,本发明的目的是提供一种基于超声导波的双向时间反演损伤成像方法。
本发明的目的通过以下的技术方案来实现:
A采用调制的激励信号对损伤实物进行主动应力波探测,并在被测区域的投射端与反射端两端阵列上同步采集信号;
B对采集到的信号进行频散补偿分离,分离出纵向模态导波和弯曲模态导波;
C提取分离后的信号中的纵向模态导波,做时间反演处理,并放大反演后信号;
D将反演后信号导入有限元模型中进行仿真反演,并提取模型中所有点云所有时刻的位移值,进行时间、空间聚焦处理,得到振动云图;
E定义空间曲线将振动云图转为三维彩色云图,以三维彩色点云形式进行损伤成像。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:
1)采用基于编码激励的主动超声导波探测,能够有效地增强信号幅度,提高信噪比,实现对细长状结构的长距离检测。
2)通过结合投射法、反射法与时间反演法,实现了对结构中损伤目标的聚焦检测,提高了损伤目标检测的分辨率与准确性。
3)将检测结果进行三维成像处理,更加直观清晰。
附图说明
图1是本发明提供的方法流程图;
图2是本发明的超声导波双向时间反演检测原理图;
图3是三维仿真实施案例模型图;
图4是同步采集传感器阵列两端的信号图;
图5是信号分离后提取的纵向模态导波图;
图6是案例模型Z方向双向时间反演聚焦振动分布图;
图7是案例模型中定义的RGB三维空间;
图8是案例模型中损伤三维点云图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合实施例及附图对本发明作进一步详细的描述。
如图1所示,展示了基于超声导波的双向时间反演损伤检测方法,该方法通过选择一个细长状的检测模型案例进行说明(如图2中( a) 和( b) 所示):
步骤10采用调制的激励信号对损伤实物进行主动应力波探测,并在被测区域的投射端与反射端两端阵列上同步采集信号;
在激励换能器阵列上加载初始激励信号:所述初始激励信号为汉宁窗调制的5峰值正弦窄带信号,在实施案例中,选取75Khz,幅值3Vpp,经功率放大器放大50x,所述换能器整理采用基于压电效应工作的长度型压电片,采用耦合剂粘贴与钢棒表面,实现振动与电信号的转换。
对损伤两端——透射端与反射端的信号进行同步采集:波形发生器的 sync同步接口与数据采集卡的外触发接口相连,使得信号激励时采集卡开始采集,保证每次采样的时间起点均一致。数据采集卡为多通道采集卡,各个通道之间无延迟。透射端与反射端阵列分别接入数据采集卡的不同通路,保证在每一次激励时,两端信号均是同步采集具体实施案例见图3、4。
参见图3,为模型案例中的纵向模态导波在损伤处的模态转换情况。低频纵向模态在长杆状钢棒中的波速最快,遇到损伤后,反射端、透射端均会接收到转换后的纵向模态及弯曲模态导波。
参见图4,为模型案例中反射端、透射端提取到的节点位移值,以此表征长杆状钢棒表面振动情况。
步骤20对采集到的信号进行频散补偿分离,分离出纵向模态导波和弯曲模态导波;
对采集到的信号做频散补偿以分离模态:频散补偿可以理解为频散过程的逆过程,已知理论频散传递函数即可对导波频散信号加以补偿。设激励信号F(ω)激发出纵向模态L和弯曲模态F,在距离d上产生的传递函数为HL(ω) 和HF(ω),则接收到的混合模态导波信号为:GLF(ω)=F(ω)HL(ω)+F(ω)HF(ω),先用对GLF(ω)进行补偿,则得对上式进行傅里叶逆变换得到时域信号,将F(ω)从补偿后的信号中去除,再将结果变为频域,再用HL(ω)进行反补偿,即则可得到分离后的弯曲模态F的频谱,再对其进行傅里叶逆变换,可得对应的时域波形。同理,可分离出纵向L模态导波。
步骤30提取分离后的信号中的纵向模态导波,做时间反演处理,并放大反演后信号;
14对分离后的信号做时间反演处理:时间反演是指对采集到的信号,在时域上将其反转的一种逆序处理方法。在频域上,这种处理等效于相位共轭。在实施案例中,将分离后的信号导入数字处理软件中,按照信号到达阵列的顺序进行前后反转,生成新的反演信号。具体实施例如图5所示。
参见图5,为模型案例中透射端、反射端经模态分离后的所提取的纵向模态导波,并进行时间反演后的结果。
步骤40将反演后信号导入有限元模型中进行仿真反演,并提取模型中所有点云所有时刻的位移值,进行时间、空间聚焦处理,得到振动云图;
将反演信号作为新的激励信号,导入有限元仿真模型中进行激励反演:根据被测物相关信息,在有限元软件中建立与实物想对应的模型,定义其几何参数、材料参数、仿真参数,将步骤4中的反演信号分别加载到仿真模型中与实物接收阵列相对应的位置,设置所需提取的特征参量,实施案例中选择提取节点的位移作为特征量。提交作业进行运算得到结果文件。
将结果文件中的节点位移导出,导入数据处理软件做时间、空间聚焦处理:结果文件中,包含了每个时刻每个节点的位移值,将同一时刻同一节点坐标下,在不同端激励时产生的位移值相乘,再将所有时刻下的计算结果叠加,将叠加的结果作为此节点坐标的位移值。具体实施案例如图6,为案例模型Z方向双向时间反演聚焦的振动分布图。
步骤50定义空间曲线将振动云图转为三维彩色云图,以三维彩色点云形式进行损伤成像。
是自定义一个空间曲线将时间、空间聚焦后的位移值映射到RGB三维坐标中。处理后过后的位移值,是从一个最小值到最大值的线性递增集合,不利于成像显示。定义一个空间坐标,参见图7,xyz三个坐标轴分别代表颜色 RGB值,从中定义一条空间曲线,将位移值与此曲线映射,可使位移值中的分布于明显的色阶中,则损伤位置及其区域可明显被区分出来。将案例中图 6中信号进行处理得到图8的案例模型中损伤三维点云图。
上述实施例具有较高的成像准确率及抗干扰能力,解决了传统方法存在的准确率低、可视化程度低等不足。在复杂环境下,实现了对绞线、管道、钢杆等各类复杂长状结构的损伤三维成像过程。并从整体上分析损伤的特征,创新性地提取不同模态导波进行透射、反射两端同时做时间反演分析,具有重要的学术价值和实际意义。
虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (6)
1.基于超声导波的双向时间反演损伤成像方法,其特征在于,所述方法包括:
A采用调制的激励信号对损伤实物进行主动应力波探测,并在被测区域的投射端与反射端两端阵列上同步采集信号;所述损伤实物的形状为细长状实心结构;
B对采集到的信号进行频散补偿分离,分离出纵向模态导波和弯曲模态导波;
C提取分离后的信号中的纵向模态导波,做时间反演处理,并放大反演后信号;
D将反演后信号导入有限元模型中进行仿真反演,并提取模型中所有点云所有时刻的位移值,进行时间、空间聚焦处理,得到振动云图;
E定义空间曲线将振动云图转为三维彩色云图,以三维彩色点云形式进行损伤成像。
2.如权利要求1所述的基于超声导波的双向时间反演损伤成像方法,其特征在于,所述步骤A具体包括:采用主动超声导波探测技术进行初始激励时,加载到阵列上的初始激励信号选用调制窄带信号,中心频率选择被测物理论频散曲线中较低频处的平缓位置,并在被检测区域两端布置接收阵列以同步采集反射端、透射端阵列上的信号。
3.如权利要求1所述的基于超声导波的双向时间反演损伤成像方法,其特征在于,所述步骤B具体包括:对于采集信号,采用频散补偿的方法进行补偿、剔除、反补偿,逐一分离出纵向模态导波、弯曲模态导波。
4.如权利要求1所述的基于超声导波的双向时间反演损伤成像方法,其特征在于,所述步骤C具体包括:选取反射端、透射端分离后的信号中的纵向模态,进行滤波后做时间反演处理,建立与损伤实物对应的无损伤有限元仿真模型,将反演信号导入有限元模型中对应损伤实物的接收位置,再次激励仿真反演。
5.如权利要求1所述的基于超声导波的双向时间反演损伤成像方法,其特征在于,所述步骤D中,提取模型所有点云所有时刻的位移值作为表面振动特征量,进行时间、空间聚焦处理,得到振动云图。
6.如权利要求1所述的基于超声导波的双向时间反演损伤成像方法,其特征在于,所述步骤E具体包括:定义空间曲线,将二维振动值映射到RGB三维空间坐标中,对损伤区域有明显色阶过渡,将振动云图转为三维色阶云图,以三维彩色点云的形式进行损伤成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510852380.6A CN105510444B (zh) | 2015-11-27 | 2015-11-27 | 基于超声导波的双向时间反演损伤成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510852380.6A CN105510444B (zh) | 2015-11-27 | 2015-11-27 | 基于超声导波的双向时间反演损伤成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105510444A CN105510444A (zh) | 2016-04-20 |
CN105510444B true CN105510444B (zh) | 2018-04-13 |
Family
ID=55718575
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510852380.6A Active CN105510444B (zh) | 2015-11-27 | 2015-11-27 | 基于超声导波的双向时间反演损伤成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510444B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106018551A (zh) * | 2016-05-03 | 2016-10-12 | 中国计量大学 | 一种基于多通道时间反转法的铝管缺陷检测定位方法 |
CN106596737B (zh) * | 2016-11-29 | 2019-12-27 | 大连理工大学 | 一种基于Lamb波波结构的多模态信号分离方法 |
CN107300787B (zh) * | 2017-06-12 | 2019-07-16 | 中国科学院上海光学精密机械研究所 | 在双层生物散射组织间三维光场的调控方法和装置 |
CN108226293B (zh) * | 2017-12-26 | 2021-05-07 | 中国航空工业集团公司基础技术研究院 | 一种用于夹芯结构超声反射透射并行检测的信号采集装置 |
CN108362358B (zh) * | 2017-12-29 | 2020-01-07 | 广州供电局有限公司 | 瓷套管液位检测方法、装置、计算机设备和存储介质 |
CN108680649A (zh) * | 2018-05-30 | 2018-10-19 | 广州特种承压设备检测研究院 | 一种内衬防腐蚀管道的缺陷成像方法 |
CN109781849B (zh) * | 2019-01-25 | 2021-02-12 | 江苏大学 | 一种基于频域分模态叠加损伤成像方法 |
CN110441398B (zh) * | 2019-07-17 | 2022-03-18 | 复旦大学 | 一种基于多层介质声速模型的合成孔径超声成像方法 |
CN111812207A (zh) * | 2020-07-21 | 2020-10-23 | 大连理工大学 | 基于超声导波转换模态提取的无基准损伤诊断成像方法 |
CN113567560B (zh) * | 2021-08-06 | 2022-12-13 | 东南大学 | 一种基于超声导波的含附属结构管道的损伤检测方法 |
CN114384152B (zh) * | 2022-01-13 | 2023-09-01 | 山东大学 | 基于搜索点匹配的超声导波损伤定位方法及系统 |
CN114878041B (zh) * | 2022-05-06 | 2023-09-01 | 中国石油大学(华东) | 一种利用双向超声探头测量在役油气管道应力的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101982943A (zh) * | 2010-10-11 | 2011-03-02 | 复旦大学 | 基于时域重排的超声导波频散补偿与多模式分离算法 |
CN103105434A (zh) * | 2013-01-17 | 2013-05-15 | 华南理工大学 | 基于主动声发射法的长形金属结构损伤检测方法及系统 |
CN104034801A (zh) * | 2014-05-13 | 2014-09-10 | 华南理工大学 | 基于合成时反的结构损伤迭代聚焦成像监测方法 |
CN104049036A (zh) * | 2014-07-01 | 2014-09-17 | 华南理工大学 | 基于时反聚焦峰值的结构损伤多峰值回归评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5893538B2 (ja) * | 2012-09-20 | 2016-03-23 | 日立Geニュークリア・エナジー株式会社 | ガイド波を用いた非破壊検査方法及び装置 |
-
2015
- 2015-11-27 CN CN201510852380.6A patent/CN105510444B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101982943A (zh) * | 2010-10-11 | 2011-03-02 | 复旦大学 | 基于时域重排的超声导波频散补偿与多模式分离算法 |
CN103105434A (zh) * | 2013-01-17 | 2013-05-15 | 华南理工大学 | 基于主动声发射法的长形金属结构损伤检测方法及系统 |
CN104034801A (zh) * | 2014-05-13 | 2014-09-10 | 华南理工大学 | 基于合成时反的结构损伤迭代聚焦成像监测方法 |
CN104049036A (zh) * | 2014-07-01 | 2014-09-17 | 华南理工大学 | 基于时反聚焦峰值的结构损伤多峰值回归评估方法 |
Non-Patent Citations (2)
Title |
---|
A new multichannel time reversal focusing method for circumferential Lamb waves and its applications for defect detection in thick-walled pipe with large-diameter;Zenghua Liu et al.;《Ultrasonics》;20140515;第54卷;第1967-1976页 * |
长骨中超声导波的频散反演方法研究;许凯亮等;《声学技术》;20131231;第32卷(第6期);第34-37页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105510444A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105510444B (zh) | 基于超声导波的双向时间反演损伤成像方法 | |
Kudela et al. | Identification of cracks in thin-walled structures by means of wavenumber filtering | |
CN102053254B (zh) | 一种激光超声波检测系统及其检测方法 | |
CN103492855B (zh) | 使用非聚焦超声的超声测振 | |
CN104020495B (zh) | 一种基于探地雷达的地下管线参数自识别方法 | |
CN109143208B (zh) | 一种船舶声纳平台自噪声试验测试系统及水平测试方法 | |
CN103698748B (zh) | 复合材料中基于信号波速与衰减补偿的2d‑music冲击定位方法 | |
CN106404911B (zh) | 用于板状结构检测的真时延单模态Lamb波相控阵系统 | |
CN108037410B (zh) | 一种基于可控响应功率的局部放电超声波定位方法与装置 | |
CN105136905B (zh) | 基于时空四维宽频阵列的高速钢轨超声探测成像方法与装置 | |
CN106198742B (zh) | 相控阵电磁超声方法和装置 | |
CN108645498A (zh) | 基于相敏光反射和卷积神经网络深度学习的冲击定位方法 | |
CN104655714A (zh) | 基于宽频磁波反射通路参数辨识的检测与成像方法及装置 | |
CN107271095A (zh) | 一种基于微波干涉的索力遥感测量方法 | |
CN104407340A (zh) | 拖曳线列阵阵形标定装置及方法 | |
CN105223271A (zh) | 一种有限空间条件下水声无源材料隔声量时反聚焦测量方法 | |
CN104034801B (zh) | 基于合成时反的结构损伤迭代聚焦成像监测方法 | |
Druet et al. | Passive guided waves measurements using fiber Bragg gratings sensors | |
CN109781864A (zh) | 基于导波信号稀疏分解方法的信号降噪与缺陷检测方法及缺陷信号重构方法 | |
AU2017376406A1 (en) | Seismic surveys with increased shot point intervals for far offsets | |
Chen et al. | Pipeline two-dimensional impact location determination using time of arrival with instant phase (TOAIP) with piezoceramic transducer array | |
CN106769734B (zh) | 一种超声波聚焦式河流泥沙浓度在线测量方法 | |
CN102262242B (zh) | 地震勘探资料采集中精确定位干扰源位置的方法 | |
CN103376443A (zh) | 探地雷达探测地面干扰快速消除方法 | |
CN105022067B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200731 Address after: 510640 five mountain road, Guangzhou, Guangdong, Tianhe District, South China University of Technology Co-patentee after: GUANGDONG GOWORLD Co.,Ltd. Patentee after: SOUTH CHINA University OF TECHNOLOGY Address before: 510640 five mountain road, Guangzhou, Guangdong, Tianhe District, South China University of Technology Patentee before: SOUTH CHINA University OF TECHNOLOGY |
|
TR01 | Transfer of patent right |