CN104897780B - 一种利用声发射信号能量对声发射源进行定位的方法 - Google Patents

一种利用声发射信号能量对声发射源进行定位的方法 Download PDF

Info

Publication number
CN104897780B
CN104897780B CN201510272307.1A CN201510272307A CN104897780B CN 104897780 B CN104897780 B CN 104897780B CN 201510272307 A CN201510272307 A CN 201510272307A CN 104897780 B CN104897780 B CN 104897780B
Authority
CN
China
Prior art keywords
mrow
msup
msub
mtr
mtd
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.)
Expired - Fee Related
Application number
CN201510272307.1A
Other languages
English (en)
Other versions
CN104897780A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201510272307.1A priority Critical patent/CN104897780B/zh
Publication of CN104897780A publication Critical patent/CN104897780A/zh
Application granted granted Critical
Publication of CN104897780B publication Critical patent/CN104897780B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明涉及一种利用声发射信号能量对声发射源进行定位的方法,属于力学传感器定位技术领域。其具体操作步骤为:①在被检测物上布置多个声发射传感器;②声发射传感器实时采集声发射信号能量;③利用声发射信号能量确定声发射源的位置坐标。本发明方法利用到弹性波衰减与传播距离的关系,直接对声发射源位置进行定位,整个计算过程中不需要测量波速,因而避免了弹性波波速测量的偏差对声发射源定位造成的影响。

Description

一种利用声发射信号能量对声发射源进行定位的方法
技术领域
本发明涉及一种利用声发射信号能量对声发射源进行定位的方法,属于力学传感器定位技术领域。
背景技术
工程材料在应用过程或由于载荷的复杂性,材料内部会出现微损伤,如裂纹或空洞。在外载下,这些微缺陷会进一步扩展而导致材料或结构发生失效破坏。如何检测或识别这些微缺陷并对结构的损伤程度及损伤发展趋势进行评估是工程领域重要的问题。
声发射检测技术与超声或其它无损检测方法的本质区别是,声发射传感器接收的信号是被检测对象自身发出的,材料内部的缺陷主动参与检测过程,具有其它方法不可替代的优越性。
通过声发射信号对微缺陷源进行定位,研究人员们提出了许多定位方法,如时差定位法、基于小波分析的定位方法等。其中,应用广泛的时差定位算法在计算过程中需要事先测量弹性波在被测物中的传播速度,而由于弹性波的传播会受到介质细观结构的不均匀性、被测物体的尺寸、几何边缘等因素的影响,从而进一步影响定位精度。
本发明方法利用弹性波衰减与传播距离的关系,直接对声发射源位置进行定位,整个计算过程中不需要测量波速,因而避免了弹性波波速测量的偏差对声发射源定位造成的影响。
发明内容
本发明的目的是提出一种利用声发射信号能量对声发射源进行定位的方法,用于确定被检测物上微缺陷源的位置。
本发明的目的是通过以下技术方案实现的。
本发明提出的一种利用声发射信号能量对声发射源进行定位的方法,其特征在于:其具体操作步骤为:
步骤一、在被检测物上布置n个声发射传感器。
当n个声发射传感器被布置在二维平面上时,n≥4;当n个声发射传感器被布置在三维空间中,n≥5。
步骤二、声发射传感器实时采集声发射信号能量。
在步骤一操作的基础上,n个声发射传感器实时采集声发射信号能量。
步骤三、确定声发射源的位置坐标,用符号(x,y,z)表示。
在步骤二操作的基础上,使用每个声发射传感器采集到的声发射信号能量,建立由n个关系式组成的声发射传感器和声发射源位置的方程组,如公式(1)所示。
其中,k为与测量电路及检测到的声发射信号相关的参数;α为衰减系数;Ei为第i个声发射传感器采集到的信号能量,i∈[1,n];(x,y,z)为声发射源的位置坐标;(xi,yi,zi)为第i个声发射传感器的位置坐标。
通过求解如公式(1)所示的方程组,即可得到未知量声发射源的位置坐标(x,y,z),从而确定被检测物上微缺陷源的位置。
利用声发射传感器采集到的声发射信号能量,建立如公式(1)所示的由n个关系式组成的声发射传感器和声发射源位置的方程组的推导过程为:
第1步:被检测物的材料内发生微缺陷时产生弹性波,声发射传感器检测到的正弦阻尼波电压峰值(用符号V′表示)与微缺陷扩展量的关系如公式(2)所示。
其中,Y是与微缺陷形状相关的参数;c是微缺陷大小;Δl是微缺陷扩展量;E为材料的弹性模量;v为材料的泊松比;P为应力;随着弹性波的传播距离变大,将p的衰减形式表示为指数形式,如公式(3)所示。
其中,P0为声源声压;x0为弹性波的传播距离。
第2步:对于突发型声发射信号,将声发射传感器输出的信号视为正弦衰减信号,如公式(4)所示。
V(t)=V′·e-βtsinω0t (4)
其中,V(t)是声发射传感器检测到的正弦阻尼波电压;β是检测得到的声发射信号衰减系数,对于确定的传感器,其为常数,并且β>0;t是时间;ω0是正弦衰减信号的圆频率。
第3步:被检测物的材料或结构受到载荷作用产生许多声发射信号,声发射信号的总能量(用符号E′表示)用单次声发射信号的电压表示,如公式(5)所示。
其中,R是电压测量线路的输入阻抗;T是正弦阻尼波的周期。
第4步:将公式(4)代入公式(5),化简后可得声发射信号的总能量E′与正弦阻尼波电压峰值V′的二次方成正比,如公式(6)所示。
其中,A为与正弦阻尼波形状有关的参数。
第5步:由于正弦阻尼波电压峰值V′与声源声压P0的成正比,比例系数用符号q表示,因此公式(6)可表示为公式(7)。
其中,x0为弹性波的传播距离。
第6步:令则声发射信号的总能量E′随弹性波传播距离x0的衰减可表示为公式(8)。
公式(8)可进一步表示为公式(9)。
(k-ln E′)=2αx0 (9)
根据n个声发射传感器采集到的信号能量E1,E2...En可以得到如公式(1)所示的定位方程组。
有益效果
本发明提出的一种利用声发射信号能量对声发射源进行定位的方法与已有技术相比较,其优点为:本发明方法利用到弹性波衰减与传播距离的关系,直接对声发射源位置进行定位,整个计算过程中不需要测量波速,因而避免了弹性波波速测量的偏差对声发射源定位造成的影响。
附图说明
图1为本发明具体实施方式中使用利用声发射信号能量对声发射源进行定位的方法检测混凝土长方体上的声发射源进行位置测量的操作流程示意图;
图2为本发明具体实施方式中在混凝土长方体上布置8个声发射传感器的位置示意图。
具体实施方式
下面结合附图和具体实施例对本发明技术方案作进一步说明。
本实施例中的被检测物为混凝土长方体,人为设置声发射源的位置为(-7.158,39.70,-16.21)。使用本发明提出的利用声发射信号能量对声发射源进行定位的方法检测混凝土长方体上(-7.158,39.70,-16.21)位置的声发射源进行位置测量,其操作流程如图1所示,具体操作步骤为:
步骤一、在混凝土长方体上布置8个声发射传感器,其坐标分别为(a,0,c)、(0,-b,c)、(-a,0,c)、(0,b,c)、(a,0,-c)、(0,-b,-c)、(-a,0,-c)、(0,b,-c),其中,a=50mm,b=50mm,c=75mm,具体如图2所示。图2中的数字1至数字8分别表示8个声发射传感器的序号。
步骤二、步骤一中所述8个声发射传感器实时采集声发射信号能量,如表1所示。
表1声发射传感器采集得到的能量
传感器编号 传感器能量(mv.us)
1 2881
2 2577
3 2993
4 3163
5 3059
6 2761
7 3234
8 3491
步骤三、使用每个声发射传感器采集到的信号能量,建立由8个关系式组成的声发射传感器和声发射源位置的方程组,如公式(10)所示。
由公式(10)可得公式(11)。
对公式(10)中的8个方程分别加减,可得公式(12)。
将公式(12)改写成公式(13)。
其中,用符号m代替公式(12)中第1个方程右侧的分子部分;用符号n代替公式(12)中第2个方程右侧的分子部分;用符号p代替公式(12)中第3个方程右侧的分子部分。
通过公式(13)进一步可得公式(14)。
用公式(10)中的第1个和第7个方程相除,得到公式(15)。
并将公式(14)带入公式(15),得到公式(16)。
其中,x(1)为用公式(10)中的第1个和第7个方程计算得到的声发射源位置的横坐标值;A1、B1、Δ1可通过公式(17)得到。
其中,C1=a2+c2
同理,利用公式(10)中的第2个和第8个方程相除,用公式(10)中的第3个和第5个方程相除,用公式(10)中的第4个和第6个方程相除,可分别求得用公式(10)中的第2个和第8个方程计算得到的声发射源位置的横坐标值x(2)、用公式(10)中的第3个和第5个方程计算得到的声发射源位置的横坐标值x(3)、用公式(10)中的第4个和第6个方程计算得到的声发射源位置的横坐标值x(4),并对x(1)至x(4)求平均值,即可求得声发射源位置的横坐标x=-7.3672mm。同样方法可求得声发射源位置的横坐标y=42.1944mm、z=-17.2974mm。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (1)

1.一种利用声发射信号能量对声发射源进行定位的方法,其特征在于:其具体操作步骤为:
步骤一、在被检测物上布置8个声发射传感器,其坐标分别为(a,0,c)、(0,-b,c)、(-a,0,c)、(0,b,c)、(a,0,-c)、(0,-b,-c)、(-a,0,-c)、(0,b,-c);
步骤二、声发射传感器实时采集声发射信号能量;
在步骤一操作的基础上,声发射传感器实时采集声发射信号能量;
步骤三、确定声发射源的位置坐标,用符号(x,y,z)表示;
在步骤二操作的基础上,使用每个声发射传感器采集到的声发射信号能量,建立由8个关系式组成的声发射传感器和声发射源位置的方程组,如公式(10)所示;
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>b</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>b</mi> <mo>+</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>8</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mn>4</mn> <msup> <mi>&amp;alpha;</mi> <mn>2</mn> </msup> <mo>&amp;lsqb;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>b</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中,k为与测量电路及检测到的声发射信号相关的参数;α为衰减系数;Ei为第i个声发射传感器采集到的信号能量,i∈[1,8];(x,y,z)为声发射源的位置坐标;
由公式(10)可得公式(11);
<mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>8</mn> </msub> </mrow> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>8</mn> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
对公式(10)中的8个方程分别加减,可得公式(12);
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>4</mn> <mi>&amp;alpha;</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mn>1</mn> <mrow> <mn>8</mn> <mi>a</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>8</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>2</mn> <mi>k</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>8</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>4</mn> <mi>&amp;alpha;</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mn>1</mn> <mrow> <mn>8</mn> <mi>b</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>8</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>2</mn> <mi>k</mi> <mrow> <mo>(</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>8</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>z</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>4</mn> <mi>&amp;alpha;</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mn>1</mn> <mrow> <mn>16</mn> <mi>c</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>-</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <msup> <mi>ln</mi> <mn>2</mn> </msup> <msub> <mi>E</mi> <mn>8</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>2</mn> <mi>k</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>3</mn> </msub> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>4</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>5</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>6</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>7</mn> </msub> <mo>+</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>8</mn> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
将公式(12)改写成公式(13);
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mfrac> <mi>m</mi> <mrow> <mn>32</mn> <mi>&amp;alpha;</mi> <mi>a</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mi>n</mi> <mrow> <mn>32</mn> <mi>&amp;alpha;</mi> <mi>b</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>z</mi> <mo>=</mo> <mfrac> <mi>p</mi> <mrow> <mn>64</mn> <mi>&amp;alpha;</mi> <mi>c</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
其中,用符号m代替公式(12)中第1个方程右侧的分子部分;用符号n代替公式(12)中第2个方程右侧的分子部分;用符号p代替公式(12)中第3个方程右侧的分子部分;
通过公式(13)进一步可得公式(14);
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mi>n</mi> </mrow> <mrow> <mi>b</mi> <mi>m</mi> </mrow> </mfrac> <mi>x</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>z</mi> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mi>p</mi> </mrow> <mrow> <mn>2</mn> <mi>c</mi> <mi>m</mi> </mrow> </mfrac> <mi>x</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
用公式(10)中的第1个和第7个方程相除,得到公式(15);
<mrow> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>k</mi> <mo>-</mo> <mi>ln</mi> <mi> </mi> <msub> <mi>E</mi> <mn>7</mn> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>=</mo> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>-</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <msup> <mrow> <mo>(</mo> <mi>a</mi> <mo>+</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>c</mi> <mo>+</mo> <mi>z</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
并将公式(14)带入公式(15),得到公式(16);
<mrow> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mn>2</mn> <msub> <mi>B</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mn>1</mn> </msub> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;PlusMinus;</mo> <msqrt> <msub> <mi>&amp;Delta;</mi> <mn>1</mn> </msub> </msqrt> </mrow> <mrow> <mn>2</mn> <msub> <mi>A</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
其中,x(1)为用公式(10)中的第1个和第7个方程计算得到的声发射源位置的横坐标值;A1、B1、Δ1可通过公式(17)得到;
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <msup> <mi>b</mi> <mn>2</mn> </msup> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>m</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>a</mi> <mn>2</mn> </msup> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>n</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>a</mi> <mn>2</mn> </msup> <msup> <mi>b</mi> <mn>2</mn> </msup> <msup> <mi>p</mi> <mn>2</mn> </msup> </mrow> <mrow> <msup> <mi>b</mi> <mn>2</mn> </msup> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>m</mi> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>B</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>+</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> <mi>m</mi> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;Delta;</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>4</mn> <mrow> <mo>(</mo> <msup> <msub> <mi>B</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>-</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <msub> <mi>C</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>16</mn> <msup> <msub> <mi>B</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <msub> <mi>S</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
其中,C1=a2+c2
同理,利用公式(10)中的第2个和第8个方程相除,用公式(10)中的第3个和第5个方程相除,用公式(10)中的第4个和第6个方程相除,可分别求得用公式(10)中的第2个和第8个方程计算得到的声发射源位置的横坐标值x(2)、用公式(10)中的第3个和第5个方程计算得到的声发射源位置的横坐标值x(3)、用公式(10)中的第4个和第6个方程计算得到的声发射源位置的横坐标值x(4),并对x(1)至x(4)求平均值,即可求得声发射源位置的横坐标x;同样方法可求得声发射源位置的坐标值y和z,从而确定被检测物上微缺陷源的位置。
CN201510272307.1A 2015-05-25 2015-05-25 一种利用声发射信号能量对声发射源进行定位的方法 Expired - Fee Related CN104897780B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510272307.1A CN104897780B (zh) 2015-05-25 2015-05-25 一种利用声发射信号能量对声发射源进行定位的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510272307.1A CN104897780B (zh) 2015-05-25 2015-05-25 一种利用声发射信号能量对声发射源进行定位的方法

Publications (2)

Publication Number Publication Date
CN104897780A CN104897780A (zh) 2015-09-09
CN104897780B true CN104897780B (zh) 2018-04-03

Family

ID=54030577

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510272307.1A Expired - Fee Related CN104897780B (zh) 2015-05-25 2015-05-25 一种利用声发射信号能量对声发射源进行定位的方法

Country Status (1)

Country Link
CN (1) CN104897780B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842345B (zh) * 2016-04-15 2018-09-11 重庆邮电大学 一种光纤布拉格光栅传感网络结构多裂纹定位方法
CN106248800A (zh) * 2016-07-12 2016-12-21 清华大学 一种基于圆环试验的评价混凝土温度断裂的方法
CN106198753B (zh) * 2016-08-29 2019-03-22 中国科学院武汉岩土力学研究所 一种提高声发射定位时空演化过程精度的方法
CN106500925B (zh) * 2016-10-25 2018-09-11 云南电网有限责任公司电力科学研究院 一种gis设备气体泄漏的检测系统、装置及方法
CN106706760A (zh) * 2016-12-20 2017-05-24 北京工业大学 全向性双圆形阵列的复合材料板声发射源定位方法
CN107918109B (zh) * 2017-11-08 2021-02-09 重庆大学 一种声发射源的平面定位装置及控制方法
CN108519436B (zh) * 2018-04-03 2021-01-01 北京理工大学 一种基于声发射信号统计分析的混凝土损伤评估方法
CN110542721B (zh) * 2019-09-29 2020-06-05 中南大学 一种基于随机传感器网络的声发射源定位方法和系统
CN111398433B (zh) * 2020-04-17 2020-12-25 中南大学 一种基于线性加权最小二乘法的声发射源定位方法及系统
GB2597756B (en) * 2020-08-03 2022-11-23 Crane John Uk Ltd Determining remaining lifetime of a seal based on accumulation of an acoustic emission energy
CN114152679A (zh) * 2021-10-28 2022-03-08 航天材料及工艺研究所 一种超低温液体环境下的钛合金气瓶声发射二维平面定位方法
CN115616090B (zh) * 2022-12-20 2023-03-31 中国铁路设计集团有限公司 一种基于声信号的列车部件质检方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4592034A (en) * 1982-11-15 1986-05-27 Cornell Research Foundation, Inc. Acoustic emission source location on plate-like structures using a small array of transducers
CN101477194A (zh) * 2009-02-17 2009-07-08 东南大学 一种转子碰摩声发射源定位方法
CN102721749A (zh) * 2012-07-09 2012-10-10 中国人民解放军后勤工程学院 受限空间中储罐底部的声发射检测装置和方法
CN102890265A (zh) * 2012-08-09 2013-01-23 西北工业大学 一种基于水下声传感器网络的目标被动定位方法
CN102928817A (zh) * 2012-10-18 2013-02-13 东南大学 一种利用时延估计进行转子碰摩声发射源定位的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4592034A (en) * 1982-11-15 1986-05-27 Cornell Research Foundation, Inc. Acoustic emission source location on plate-like structures using a small array of transducers
CN101477194A (zh) * 2009-02-17 2009-07-08 东南大学 一种转子碰摩声发射源定位方法
CN102721749A (zh) * 2012-07-09 2012-10-10 中国人民解放军后勤工程学院 受限空间中储罐底部的声发射检测装置和方法
CN102890265A (zh) * 2012-08-09 2013-01-23 西北工业大学 一种基于水下声传感器网络的目标被动定位方法
CN102928817A (zh) * 2012-10-18 2013-02-13 东南大学 一种利用时延估计进行转子碰摩声发射源定位的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种声发射源的新型平面定位方法研究;龚斌等;《声学技术》;20060430;第107-109页 *
声发射检测中利用能量进行定位的新方法;朱祥军;《中国测试》;20110131;2 波的衰减特性和能量分析,3 能量定位法 *

Also Published As

Publication number Publication date
CN104897780A (zh) 2015-09-09

Similar Documents

Publication Publication Date Title
CN104897780B (zh) 一种利用声发射信号能量对声发射源进行定位的方法
CN108444653B (zh) 一种面阵型超声波气体泄漏点定位检测装置及方法
CN202916242U (zh) 气液两相流流型在线监测装置
CN107013813B (zh) 一种供水管道泄漏量估算系统及方法
CN103968883B (zh) 淤泥探测方法及淤泥探测仪
CN102865952A (zh) 一种混凝土工作应力的无损测试方法
CN103926313B (zh) 一种基于超声检测的复合材料孔隙率数值评估方法
CN107356677B (zh) 基于旅行时层析成像和逆时偏移成像的超声无损检测方法
CN104374532A (zh) 航天器在轨泄漏定向方法
CN110274956A (zh) 一种基于Kolmogorov熵的管道超声导波检测方法
CN114235962B (zh) 一种面向各向异性结构的超声导波成像方法及系统
CN109163769B (zh) 一种管道流量电磁阵列传感器的检测方法
CN105403878A (zh) 一种基于时延敏感核的海洋声层析方法
CN104101648A (zh) 基于李雅普诺夫指数的超声导波定位缺陷的方法
CN109187769B (zh) 一种基于超声导波定量化检测管道缺陷的方法
CN108036200B (zh) 基于重心频率衰减特性的单声发射传感器泄漏源定位方法
CN204758542U (zh) 一种金属结构的裂纹的检测装置
CN107063526A (zh) 基于临界折射纵波的钢构件绝对应力分布的检测方法
CN205352653U (zh) 一种低速风洞内的超声精确测距系统
CN106769734B (zh) 一种超声波聚焦式河流泥沙浓度在线测量方法
CN107504936A (zh) 淤泥检测系统及淤泥检测方法
CN203053483U (zh) 淤泥探测仪
CN109946023A (zh) 一种管道气体泄漏判别装置及判识方法
CN103604570A (zh) 一种超声波气密性检测方法及装置
CN103616436B (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180403

Termination date: 20190525