CN107807173A - 基于幂函数的平板结构突发型声发射源定位方法 - Google Patents

基于幂函数的平板结构突发型声发射源定位方法 Download PDF

Info

Publication number
CN107807173A
CN107807173A CN201711319012.0A CN201711319012A CN107807173A CN 107807173 A CN107807173 A CN 107807173A CN 201711319012 A CN201711319012 A CN 201711319012A CN 107807173 A CN107807173 A CN 107807173A
Authority
CN
China
Prior art keywords
mrow
acoustic emission
msup
power function
peak
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
Application number
CN201711319012.0A
Other languages
English (en)
Other versions
CN107807173B (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.)
Shijiazhuang Tiedao University
Original Assignee
Shijiazhuang Tiedao University
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 Shijiazhuang Tiedao University filed Critical Shijiazhuang Tiedao University
Priority to CN201711319012.0A priority Critical patent/CN107807173B/zh
Publication of CN107807173A publication Critical patent/CN107807173A/zh
Application granted granted Critical
Publication of CN107807173B publication Critical patent/CN107807173B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/14Investigating 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 using acoustic emission techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/48Processing the detected response signal, e.g. electronic circuits specially adapted therefor by amplitude comparison
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/01Indexing codes associated with the measuring variable
    • G01N2291/015Attenuation, scattering

Landscapes

  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Acoustics & Sound (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于幂函数的平板结构突发型声发射源定位方法,涉及声发射无损检测技术领域。所述方法包括以下步骤:在平板上建立坐标系,将四个声发射传感器布置于相应点位,进而搭建声发射波幅值衰减实验系统,对突发型声发射信号作小波降噪处理后,提取出四个信号的弯曲波首峰峰值,建立幅值衰减与传播距离的幂函数计算关系,最终得到声发射源的坐标值,该方法不仅能够实现突发型声发射源定位,而且与常用的时差定位不同,不需要时间同步和预知声速,定位更加快速。

Description

基于幂函数的平板结构突发型声发射源定位方法
技术领域
本发明涉及声发射无损检测领域,尤其涉及一种基于幂函数的平板结构突发型声发射源定位方法。
背景技术
声发射技术作为一种无损检测方法,已经广泛应用于航空航天、加工制造、交通运输等领域,其重要研究内容之一是由接收的信号反推到声发射源的问题,即源定位。根据声发射信号的产生率将其划分为突发型和连续型两种,其中突发型信号源定位采用的时差定位,是当前声发射源定位最为广泛的方法,但其需要假设材料声速为常数,且该关键参数易受多种因素影响而难以准确测量,致使时差定位方法在实际检测中存在一定的局限性。而对于连续型声发射信号应用的幅值衰减模型为指数函数,定位精度低,仅适合于区域定位,而且仅仅适应于连续型声发射信号,对于突发型声发射信号无能为力。
发明内容
本发明所要解决的技术问题是提供一种基于幂函数的平板结构突发型声发射源定位方法,不仅避免时差定位需要测量声速的桎梏,而且能够完成突发型信号的定位,并且具有定位准确、简单快捷的功能,从而提高声发射技术的更为广泛的应用性。
为解决上述技术问题,本发明提供了一种基于幂函数的平板结构突发型声发射源定位方法,主要包括下列步骤:
(1)在金属平板上建立直角坐标系,并在金属平板的相应点位分布四个声发射传感器;
(2)搭建声发射波幅值衰减定位系统,进行相关参数设定并调试;
(3)声发射波幅值衰减定位系统通过四个声发射传感器采集每组断铅事件对应的突发型信号;
(4)选取一个断铅事件所对应的四个声发射传感器信号,利用声发射波幅值衰减定位系统对发射传感器信号进行小波降噪处理,并提取出四个信号的扩展波首峰峰值;
(5)将断铅事件所对应的四个声发射传感器的距离采用信号扩展波的首峰峰值来表示,求得幂函数模型V=ax-b中的参数a、b值;
(6)计算出断铅位置的坐标值点P(x,y);
(7)重复(5)-(6)步骤,计算其它断铅位置的坐标值,完成声发射源定位。
所述步骤(1)包括在待检测的金属平板上建立直角坐标系并绘制网格,在以原点为中心的正方形四个顶点处分别布置一个型号为R15α的声发射传感器,进而建立断铅点P(x,y)与四个传感器的位置关系,
其中L1~L4为断铅点P(x,y)到四个传感器的距离,x、y断铅点的坐标值,m为相邻两个声发射传感器连线与坐标轴交点对应的坐标值。
所述步骤(5)包括提出声发射波在平板传播因纳入扩散衰减导致其幅值呈幂函数V=ax-b形式随声源距离衰减,将断铅点P(x,y)到四个传感器的距离L1~L4用接收信号的首峰峰值V1~V4来表示,分别表示为 进而可以求得幂函数V=ax-b模型中的参数a′、b′值。
所述步骤(6)包括代入步骤(5)参数a′、b′值,得到本次断铅信号的首峰峰值衰减模型V=a′x-b′,进而计算出本次断铅位置P的坐标值(x,y)。
采用上述技术方案所产生的有益效果在于:
(1)本发明采用了声发射波幅值衰减研究声源定位,克服了时差定位对声波速度的苛刻需求,从而取消对材料声速的测量或计算,有利于促进声发射技术的更为广泛应用;
(2)本发明提出了在平板传播的声发射波幅值衰减模型呈幂函数V=ax-b形式,并提供了相关参数的计算方法,能够定位出声发射波源的位置;
(3)本发明定义了声发射波首个弯曲波波峰幅值为首峰峰值,并采用它以评判声发射波的幅值衰减大小,提高了采用幅值衰减方法定位的精度和信度。
附图说明
通过阅读下面对本发明的示例性实施例的详细描述,同时结合附图考虑,可更好地理解本发明的上述实施例和其它实施例、目的、特性、优点、技术上和工业上的意义,各附图中:
图1是声发射波幅值衰减实验系统图;
图2是声发射幅值衰减源定位方法流程图;
图3是平板声发射传感器布设位置图;
图4是声发射系统采集的原始信号时域图(通道1);
图5是信号小波变换后的时间-频率-幅值三维图(通道1);
图6是信号小波变换后的时间-频率平面图(通道1);
图7是声发射波原始信号的弯曲波首波部分图(通道1);
图8是降噪后的声发射波信号弯曲波首波部分图(通道1)。
其中:1、铝合金平板;2、声发射传感器;3、前置放大器;4、计算机;5、信号分析仪。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
本发明提出了一种基于幅值幂函数衰减的平板结构突发型声发射源定位方法,所测平板材料为铝合金,采用断铅方法模拟声发射源,搭建声发射波衰减实验系统如图1所示,采用的定位方法流程如图2所示,具体验证步骤如下:
(1)如图3所示,在尺寸1000mm×1000mm×2mm的铝合金平板上,以平板中心为原点建立直角坐标系,分别在正方形的四个顶点,即坐标1(200,200)、2(-200,200)、3(-200,-200)、4(200,-200)处,各布置一个型号R15α的声发射传感器2;
(2)如图1所示,将布置好的四个声发射传感器2接入DS5系列全信息声发射信号分析仪5,设定通道门限固定为10dB、前置放大器2增益为40dB、峰值定义时间1000us、波击定义时间2000us、波击闭锁时间20000us,采样率为3MHz,并进行断铅检测各传感器安装情况和系统运行状态是否良好;
(3)在多个位置进行断铅,并采集相应不同位置断铅事件的信号数据,断铅位置选择在铝合金平板1的上表面,断铅作用力的方向与平板上表面垂直,断铅方向与铝合金平板表面成30°角;
(4)选取一次断铅事件,将此次断铅事件的四个声发射传感器2的信号数据,如图4所示,分别输入到具有小波变换功能的预设MATLAB程序(利用计算机4进行计算),得到声发射传感器信号小波变换后的时间-频率-幅值图,如图5、图6所示,可判断出该组信号需进行降噪处理;
(5)提取出四个声发射信号的扩展波波段,其信号幅值图如图7所示,经降噪处理后,如图8所示,定义该扩展波的首个波峰幅值为首峰峰值,并提取出这四个声发射传感器信号的首峰峰值,用作评判声发射波幅值的衰减大小的核心数据;
(6)将接收信号的首峰峰值V1~V4分别代入由幂函数模型V=ax-b推导出的公式中,分别计算出断铅点P(x,y)到四个声发射传感器2的距离L1~L4,利用编制的MATLAB程序,求得该次断铅事件的幅值幂函数模型中的参数a、b值,分别为a′、b′;
(7)将求得的参数a′、b′值代入到声发射幅值衰减模型中,得到本次断铅信号的首峰峰值衰减规律为V=a′x-b′,进而求得本次断铅位置P坐标x、y的数值解;
(8)若定位在其它位置的断铅点,可重复上述的步骤(4)-(7),均是首先提取出这四组信号的首峰峰值,然后求得对应于该次断铅试件的幅值衰减模型参数a、b值,最终得到相应断铅位置点的坐标值(x,y),即实现声源定位。
总之,本发明提出一种与时差定位不同的,基于幅值幂函数衰减的声发射源定位方法,从而避开测量声速的需求,提供出一种新的适合于突发型声发射源定位的方法。

Claims (4)

1.一种基于幂函数的平板结构突发型声发射源定位方法,其特征在于:主要包括下列步骤:
(1)在金属平板上建立直角坐标系,并在金属平板的相应点位分布四个声发射传感器;
(2)搭建声发射波幅值衰减定位系统,进行相关参数设定并调试;
(3)声发射波幅值衰减定位系统通过四个声发射传感器采集每组断铅事件对应的突发型信号;
(4)选取一个断铅事件所对应的四个声发射传感器信号,利用声发射波幅值衰减定位系统对发射传感器信号进行小波降噪处理,并提取出四个信号的扩展波首峰峰值;
(5)将断铅事件所对应的四个声发射传感器的距离采用信号扩展波的首峰峰值来表示,求得幂函数模型V=ax-b中的参数a、b值;
(6)计算出断铅位置的坐标值点P(x,y);
(7)重复(5)-(6)步骤,计算其它断铅位置的坐标值,完成声发射源定位。
2.根据权利要求1所述的基于幂函数的平板结构突发型声发射源定位方法,其特征在于:所述步骤(1)包括在待检测的金属平板上建立直角坐标系并绘制网格,在以原点为中心的正方形四个顶点处分别布置一个型号为R15α的声发射传感器,进而建立断铅点P(x,y)与四个传感器的位置关系,
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msup> <msub> <mi>L</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>-</mo> <mi>m</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mrow> <mi>y</mi> <mo>-</mo> <mi>m</mi> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>L</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>L</mi> <mn>3</mn> </msub> <mn>2</mn> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>+</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>L</mi> <mn>4</mn> </msub> <mn>2</mn> </msup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>+</mo> <mi>m</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中L1~L4为断铅点P(x,y)到四个传感器的距离,x、y断铅点的坐标值,m为相邻两个声发射传感器连线与坐标轴交点对应的坐标值。
3.根据权利要求2所述的基于幂函数的平板结构突发型声发射源定位方法,其特征在于:所述步骤(5)包括提出声发射波在平板传播因纳入扩散衰减导致其幅值呈幂函数V=ax-b形式随声源距离衰减,将断铅点P(x,y)到四个传感器的距离L1~L4用接收信号的首峰峰值V1~V4来表示,分别表示为(i=1、2、3、4),进而可以求得幂函数V=ax-b模型中的参数a′、b′值。
4.根据权利要求3所述的基于幂函数的平板结构突发型声发射源定位方法,其特征在于:所述步骤(6)包括代入步骤(5)参数a′、b′值,得到本次断铅信号的首峰峰值衰减模型V=a′x-b′,进而计算出本次断铅位置P的坐标值(x,y)。
CN201711319012.0A 2017-12-12 2017-12-12 基于幂函数的平板结构突发型声发射源定位方法 Active CN107807173B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711319012.0A CN107807173B (zh) 2017-12-12 2017-12-12 基于幂函数的平板结构突发型声发射源定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711319012.0A CN107807173B (zh) 2017-12-12 2017-12-12 基于幂函数的平板结构突发型声发射源定位方法

Publications (2)

Publication Number Publication Date
CN107807173A true CN107807173A (zh) 2018-03-16
CN107807173B CN107807173B (zh) 2019-12-10

Family

ID=61589110

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711319012.0A Active CN107807173B (zh) 2017-12-12 2017-12-12 基于幂函数的平板结构突发型声发射源定位方法

Country Status (1)

Country Link
CN (1) CN107807173B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110161126A (zh) * 2019-06-27 2019-08-23 中国工程物理研究院化工材料研究所 适用于宽温域的固体应力波幅值衰减系数测试装置及方法
CN112525999A (zh) * 2020-11-30 2021-03-19 北京强度环境研究所 一种基于频率成份的声发射信号相关性分析方法
CN115166057A (zh) * 2022-08-01 2022-10-11 江苏大学 一种基于声发射波速变化测量方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU873107A1 (ru) * 1979-04-11 1981-10-15 Anisimov Vladimir K Способ определени координат источников сигналов акустической эмиссии в издели х из листовых материалов
JP2001013118A (ja) * 1999-06-30 2001-01-19 Sumitomo Metal Ind Ltd 電磁超音波探触子
CN101477084A (zh) * 2008-12-22 2009-07-08 铁道第三勘察设计院集团有限公司 采用弹性波透射方式检测混凝土桩桩身完整性的方法
CN101769900A (zh) * 2010-01-25 2010-07-07 天津大学 声发射储罐罐底腐蚀检测中声源的定位方法
CN102944613A (zh) * 2012-11-16 2013-02-27 中国科学院半导体研究所 一种光纤声发射检测与定位系统
CN103472435A (zh) * 2013-09-23 2013-12-25 中国石油大学(华东) 海洋平台复杂结构系统声发射源两阶段定位方法
CN105425212A (zh) * 2015-11-18 2016-03-23 西北工业大学 一种声源定位方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU873107A1 (ru) * 1979-04-11 1981-10-15 Anisimov Vladimir K Способ определени координат источников сигналов акустической эмиссии в издели х из листовых материалов
JP2001013118A (ja) * 1999-06-30 2001-01-19 Sumitomo Metal Ind Ltd 電磁超音波探触子
CN101477084A (zh) * 2008-12-22 2009-07-08 铁道第三勘察设计院集团有限公司 采用弹性波透射方式检测混凝土桩桩身完整性的方法
CN101769900A (zh) * 2010-01-25 2010-07-07 天津大学 声发射储罐罐底腐蚀检测中声源的定位方法
CN102944613A (zh) * 2012-11-16 2013-02-27 中国科学院半导体研究所 一种光纤声发射检测与定位系统
CN103472435A (zh) * 2013-09-23 2013-12-25 中国石油大学(华东) 海洋平台复杂结构系统声发射源两阶段定位方法
CN105425212A (zh) * 2015-11-18 2016-03-23 西北工业大学 一种声源定位方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110161126A (zh) * 2019-06-27 2019-08-23 中国工程物理研究院化工材料研究所 适用于宽温域的固体应力波幅值衰减系数测试装置及方法
CN110161126B (zh) * 2019-06-27 2024-02-13 中国工程物理研究院化工材料研究所 适用于宽温域的固体应力波幅值衰减系数测试装置及方法
CN112525999A (zh) * 2020-11-30 2021-03-19 北京强度环境研究所 一种基于频率成份的声发射信号相关性分析方法
CN112525999B (zh) * 2020-11-30 2023-05-23 北京强度环境研究所 一种基于频率成份的声发射信号相关性分析方法
CN115166057A (zh) * 2022-08-01 2022-10-11 江苏大学 一种基于声发射波速变化测量方法

Also Published As

Publication number Publication date
CN107807173B (zh) 2019-12-10

Similar Documents

Publication Publication Date Title
US11900630B2 (en) Method for detecting leakage of water supply pipe based on ground-penetrating radar three-dimensional image attribute analysis
Aljets et al. Acoustic emission source location on large plate-like structures using a local triangular sensor array
CN106770664B (zh) 一种基于全聚焦成像算法改进边缘缺陷检测的方法
CN107807173A (zh) 基于幂函数的平板结构突发型声发射源定位方法
CN103675920B (zh) 道路基层隐含裂缝深度和水平位置的无损检测方法
CN109696480A (zh) 一种基于改进时间反转算法的玻璃纤维复合材料声发射源定位成像方法
CN107167521A (zh) 一种提高识别混凝土缺陷和损伤能力的检测方法
CN206090647U (zh) 一种桩基检测设备
CN106706760A (zh) 全向性双圆形阵列的复合材料板声发射源定位方法
CN109374748A (zh) 声发射源时反定位成像方法
CN105973752A (zh) 集成透射式和反射式超声过程层析成像方法
CN104101648A (zh) 基于李雅普诺夫指数的超声导波定位缺陷的方法
CN108519583A (zh) 适用于各向异性二维板的声发射源定位方法
CN105954358A (zh) 一种TR与Duffing系统相结合的超声导波小缺陷定位检测方法
CN103630605B (zh) 一种预应力锚索管道注浆检测方法
CN202837253U (zh) 一种火车车轮超声波探伤用试块组件
CN107642355B (zh) 基于超声波发射法的水力压裂裂缝监测系统及方法
CN205844271U (zh) 一种基于差频非线性超声检测金属薄板微裂纹的检测系统
CN103604869B (zh) 基于数值反演的无损检测模拟试块缺陷参数的识别方法
CN103376443B (zh) 探地雷达探测地面干扰快速消除方法
Freeseman et al. Characterization of concrete at various freeze-thaw damage conditions using SH-waves
CN108267509A (zh) 基于Lamb波的损伤图像化系统
CN103969334A (zh) 大尺寸混凝土构件缺陷应力波快速检测方法
US11946735B2 (en) Method for measuring deviation angle of fatigue microcrack based on nonlinear ultrasound
CN106770672A (zh) Lamb波主被动监测信息融合的工程结构损伤评估方法

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180316

Assignee: HEBEI XINDADI ELECTROMECHANICAL MANUFACTURING Co.,Ltd.

Assignor: SHIJIAZHUANG TIEDAO University

Contract record no.: X2024980008425

Denomination of invention: A Method for Localizing Sudden Acoustic Emission Sources in Flat Plate Structures Based on Power Functions

Granted publication date: 20191210

License type: Common License

Record date: 20240627

EE01 Entry into force of recordation of patent licensing contract