CN105717201B - 基于声场波数空间谱的轴对称体缺陷检测重构方法 - Google Patents

基于声场波数空间谱的轴对称体缺陷检测重构方法 Download PDF

Info

Publication number
CN105717201B
CN105717201B CN201610051850.3A CN201610051850A CN105717201B CN 105717201 B CN105717201 B CN 105717201B CN 201610051850 A CN201610051850 A CN 201610051850A CN 105717201 B CN105717201 B CN 105717201B
Authority
CN
China
Prior art keywords
mrow
msub
delta
msup
mfrac
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
Application number
CN201610051850.3A
Other languages
English (en)
Other versions
CN105717201A (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.)
North University of China
Original Assignee
North University of China
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 North University of China filed Critical North University of China
Priority to CN201610051850.3A priority Critical patent/CN105717201B/zh
Publication of CN105717201A publication Critical patent/CN105717201A/zh
Application granted granted Critical
Publication of CN105717201B publication Critical patent/CN105717201B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/449—Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
    • 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/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • 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/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/0289—Internal structure, e.g. defects, grain size, texture

Abstract

本发明涉及基于声场波数空间谱的轴对称体缺陷检测重构方法，采用轴对称体超声检测系统，实现轴对称体的全面检测，利用超声应力波的传播规律和声场波数空间谱，通过有限差分迭代法计算得出反演声场和一个检测截面二维重构结果，最后利用体绘制法重构出缺陷的三维形状。本发明将所检测截面采用有限元法进行单元分割，将超声回波信号的反褶信号加载在超声阵元处，通过建立平面坐标x、y变量与声场波数kx，ky之间的二维傅里叶变换关系，根据应力波的传递原理，计算传播过程的声压，最终将所有阵元的声波传递声场进行叠加，重构出检测截面的缺陷形状；考虑了多阵元信号的相互作用，大大提高了重构精度，利用声场波数空间谱和快速傅里叶变换提高了运算速度。

Description

基于声场波数空间谱的轴对称体缺陷检测重构方法
技术领域
本发明属于轴对称体缺陷检测方法技术领域，具体涉及一种基于声场波数空间谱的轴对称体缺陷检测重构方法。
背景技术
轴对称结构是材料加工、产品生产的主要结构，轴对称体广泛应用于国防军工、航空航天、交通运输等领域，一般由金属或合金材料经温挤压或铸造工艺生产而成，在生产过程中由于温度的瞬间变化、介质与模具的摩擦、交变载荷等原因会造成工件内部和表面产生不同程度缺陷，这些缺陷会影响产品正常、安全的使用。因此需要在生产过程中和产品使用前对其进行无损检测和评价。
超声波探伤方法具有指向性好、高灵敏度、性能稳定的特点，是金属和合金材料缺陷检测的常用方法，能够检测材料表面和内部的缺陷，对如何从众多的超声回波信号中提取缺陷信息，并对缺陷进行重构是该领域的难点。现有的方法一般是根据单一回波提取缺陷特征信号，并根据特征信号的时间延时和超声波的传播速度，计算得出缺陷的边界位置，最终把所有的缺陷表面连接起来就得到缺陷的最终轮廓。这种方法能够初略地重构出缺陷形状，但很难准确地进行重构，主要原因是超声波传播存在扩散现象，偏离传播中心线的缺陷也可能会有反射回波，但按照上述方法进行缺陷重构时却把缺陷的边界确定到传播中心线上，会带来很大的误差。
发明内容
本发明的目的是解决现有轴对称体缺陷检测重构方法存在缺陷重构结果误差大的技术问题，提供一种基于声场波数空间谱的轴对称体缺陷检测重构方法。
为解决上述技术问题，本发明采用的技术方案为：
基于声场波数空间谱的轴对称体缺陷检测重构方法，其采用以下的步骤进行轴对称体的缺陷检测重构：
1)利用轴对称体超声检测系统，采集轴对称体每个位置的超声回波信号，截取有效区间[tmin,tmax]内的信号进行反褶运算得到反褶信号sz,θ(t)，其中z表示轴向检测位置，θ表示圆周方向旋转角度,tmin和tmax分别为超声回波信号中的近表面回波开始位置和远表面回波结束位置，超声回波信号的采样间隔为Δt；
2)按顺序对轴对称体其中一检测截面z＝z0，采用有限元法进行单元分割，x、y方向的步长分别为Δx、Δy；初始化各有限元单元节点处的声压pz(x,y)＝0，对于超声探头所在位置的声压，进行回波信号加载pz(x,y)＝pz(x,y)+sz,θ(0)；
3)对各有限元单元节点处的声压pz(x,y)进行二维傅里叶变换得到声压波数空间谱Pz(kx,ky)＝F[pz(x,y)]，初始化应力波传播的初始密度并计算初始应变其中kx、ky分别为x、y方向的波数，c为纵波声速，p和P分别为pz(x,y)和Pz(kx,ky)的简写，j2＝-1；
4)通过有限差分迭代法计算下一时刻t+Δt时刻的应变
λ为边界匹配参数；计算t+Δt时刻的应变梯度
5)计算t+Δt时刻的密度并得出声压对于超声探头所在位置的声压，进行回波信号加载p＝p+sz,θ(t+Δt)，进行二维傅里叶变换得到声压波数空间谱P＝F[p]；
6)循环执行步骤4)-5)，直到sz,θ(t)加载完毕，并对得到的最后声压图进行阈值分割，得到z＝z0的二维重构结果；
7)循环执行步骤2)-6)，得到所有截面的二维重构结果，采用体绘制法进行三维重构即得轴对称体缺陷的三维重构结果。
本发明采用以上技术方案，与背景技术相比，本发明具有以下优点：
1)根据声波的传播原理和波动方程，考虑了超声传播过程的扩散现象，避免在重构原理上就带入了重构误差，采用该方法大大提高了重构精度；
2)采用基于声场波数空间谱的方法将超声波传播平面直角坐标系，转换为声场波数空间坐标系，并利用快速傅里叶变换及其相关性质，提高了运算速度。
为表明本发明具有以上优点，分别采用传统时间延时方法和本发明的方法对同一轴对称体的单一截面缺陷进行重构，其中，图3是采用本发明的方法得到的单一截面缺陷二维重构结果，图4是采用传统时间延时方法得到的单一截面缺陷二维重构结果，通过比较可知，采用本发明所述方法大大提高了重构精度。
附图说明
图1为本发明所使用的轴对称体超声检测系统的结构示意图；
图2为本发明基于声场波数空间谱的轴对称体缺陷重构软件的流程图；
图3为采用本发明得到的单一截面缺陷二维重构结果；
图4为传统时间延时方法得到的单一截面缺陷二维重构结果；
图5为本发明通过体绘制得到的最终缺陷三维重构结果。
具体实施方式
如图1所示，本实施例所使用的轴对称体超声检测系统，包括轴对称体旋转卡盘2、轴对称体旋转电机3、“┚”型检测台支撑座4、超声探头移动电机5、超声探头移动导轨6、超声探头7、电机控制卡8、工控机9、超声发射接收卡10及缺陷检测与处理软件，所述轴对称体旋转电机3设在“┚”型检测台支撑座4竖档内侧的上部，轴对称体旋转卡盘2安装在轴对称体旋转电机3的转轴上，轴对称体旋转卡盘2通过卡口卡住待检轴对称体1，超声探头移动电机5设在“┚”型检测台支撑座4的底座上，超声探头移动导轨6安装在超声探头移动电机5的转轴上，超声探头7安装在超声探头移动导轨6上，超声发射接收卡10和电机控制卡8插在工控机9上，并通过PCI总线相连；超声发射接收卡10通过数据线与超声探头7相连，电机控制卡8通过控制线与超声探头移动电机5和轴对称体旋转电机3相连；缺陷检测与处理软件安装在工控机9上。
本实施例中的一种基于声场波数空间谱的轴对称体缺陷检测重构方法，其采用以下的步骤进行轴对称体的缺陷检测重构：
1)利用轴对称体超声检测系统，采集轴对称体每个位置的超声回波信号，工控机9通过电机控制卡8发出指令控制超声探头移动电机5和超声探头移动导轨6，移动超声探头7位于检测位置处，工控机9通过电机控制卡8发出指令控制轴对称体旋转电机3和轴对称体旋转卡盘2，带动待检轴对称体1旋转，工控机9控制超声发射接收卡10和超声探头7采集超声回波信号，截取有效区间[tmin,tmax]内的信号进行反褶运算得到反褶信号sz,θ(t)，其中z表示轴向检测位置，θ表示圆周方向旋转角度,tmin和分别为超声回波信号中的近表面回波开始位置和远表面回波结束位置，超声回波信号的采样间隔为Δt。本实施例中，超声探头采用中心频率为5MHz的直探头，轴对称体为直径110mm，长度100mm的铝圆柱体，内部带有直径为5mm，长度为5mm的孔状缺陷，检测时轴向步进Δz＝1mm，圆周方向步进角度Δθ＝1°，超声回波信号的采样间隔为Δt＝0.01μs，tmin＝10μs，tmax＝50μs，所检测的圆周截面共100个，一个圆周截面共采集360个回波信号，每个反褶信号的序列长度为4001点；
2)按顺序对轴对称体其中一检测截面z＝z0(本案例共100个检测截面)，采用有限元法进行单元分割，x、y方向的步长分别为Δx、Δy，本实施例中，Δx＝0.1mm，Δy＝0.1mm，初始化各有限元单元节点处的声压pz(x,y)＝0，对于超声探头所在位置的声压，进行回波信号加载pz(x,y)＝pz(x,y)+sz,θ(0)；
3)对各有限元单元节点处的声压pz(x,y)进行二维傅里叶变换得到声压波数空间谱Pz(kx,ky)＝F[pz(x,y)]，初始化应力波传播的初始密度并计算初始应变其中kx、ky分别为x、y方向的波数，c为纵波声速取6320m/s，p和P分别为pz(x,y)和Pz(kx,ky)的简写，j2＝-1；
4)通过有限差分迭代法计算下一时刻t+Δt时刻的应变
λ为边界匹配参数，取0.5；计算t+Δt时刻的应变梯度
5)计算t+Δt时刻的密度并得出声压对于超声探头所在位置的声压，进行回波信号加载p＝p+sz,θ(t+Δt)，进行二维傅里叶变换得到声压波数空间谱P＝F[p]；
6)循环执行步骤4)-5)，直到sz,θ(t)加载完毕，并对得到的最后声压图进行阈值分割，得到z＝z0的二维重构结果,见图3；
7)循环执行步骤2)-6)，得到所有截面的二维重构结果，采用体绘制法进行三维重构即得轴对称体缺陷的三维重构结果，见图5。
如图2所示，本发明基于声场波数空间谱的轴对称体缺陷重构软件的流程为：采用轴对称超声检测系统采集每个位置的超声回波信号，截取有效区间内的反褶信号；取出第1个检测截面的360个反褶信号，采用有限元的方法对第1个检测截面进行单元分割，并初始化各有限元单元节点处的声压，反褶信号第0点的数据加载对应超声探头处；对声压进行二维傅里叶变换得到声压波数空间谱，并初始化应力波传播的初始密度和初始应变；根据初始的数据通过差分迭代法计算第1次迭代的应变、应变梯度、密度、声压分布，并将反褶信号第1点的数据加载对应超声探头处，同时计算声压波数空间谱；根据第1次迭代、加载得到的数据通过差分迭代法计算第2次迭代的应变、应变梯度、密度、声压分布，并将反褶信号第2点的数据加载对应超声探头处，同时计算声压波数空间谱；依次类推进行第2-3999次迭代；根据第3999次迭代的数据计算得到第4000次迭代的应变、应变梯度、密度、声压分布，并对此时得到的声压图进行阈值分割，得到二维重构结果；对所有100个检测截面按先后顺序进行上述计算，得到每个检测截面的二维重构结果，采用体绘制方法进行三维重构，得到缺陷三维形状及所在的位置。

Claims (1)

1.基于声场波数空间谱的轴对称体缺陷检测重构方法，其特征在于：采用以下的步骤进行轴对称体的缺陷检测重构：
1)利用轴对称体超声检测系统，采集轴对称体每个位置的超声回波信号，截取有效区间[tmin,tmax]内的信号进行反褶运算得到反褶信号sz,θ(t)，其中z表示轴向检测位置，θ表示圆周方向旋转角度,tmin和tmax分别为超声回波信号中的近表面回波开始位置和远表面回波结束位置，超声回波信号的采样间隔为Δt；
2)按顺序对轴对称体其中一检测截面z＝z0，采用有限元法进行单元分割，x、y方向的步长分别为Δx、Δy；初始化各有限元单元节点处的声压pz(x,y)＝0，对于超声探头所在位置的声压，进行回波信号加载pz(x,y)＝pz(x,y)+sz,θ(0)；
3)对各有限元单元节点处的声压pz(x,y)进行二维傅里叶变换得到声压波数空间谱Pz (kx,ky)＝F[pz(x,y)]，初始化应力波传播的初始密度并计算初始应变 <mrow> <msub> <mi>u</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> </mfrac> <mo>&amp;CenterDot;</mo> <msup> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <msub> <mi>jk</mi> <mi>x</mi> </msub> <msup> <mi>e</mi> <mrow> <msub> <mi>jk</mi> <mi>x</mi> </msub> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> </msup> <mi>P</mi> <mo>&amp;rsqb;</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <msub> <mi>u</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> </mfrac> <mo>&amp;CenterDot;</mo> <msup> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <msub> <mi>jk</mi> <mi>y</mi> </msub> <msup> <mi>e</mi> <mrow> <msub> <mi>jk</mi> <mi>y</mi> </msub> <mi>&amp;Delta;</mi> <mi>y</mi> </mrow> </msup> <mi>P</mi> <mo>&amp;rsqb;</mo> <mo>;</mo> </mrow> 其中kx、ky分别为x、y方向的波数，c 为纵波声速，p和P分别为pz(x,y)和Pz(kx,ky)的简写，j2＝-1；
4)通过有限差分迭代法计算下一时刻t+Δt时刻的应变 <mrow> <msub> <mi>u</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;lambda;</mi> <mo>&amp;CenterDot;</mo> <mo>{</mo> <msub> <mi>&amp;lambda;u</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> </mfrac> <mo>&amp;CenterDot;</mo> <msup> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <msub> <mi>jk</mi> <mi>x</mi> </msub> <msup> <mi>e</mi> <mrow> <msub> <mi>jk</mi> <mi>x</mi> </msub> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> </msup> <mi>P</mi> <mo>&amp;rsqb;</mo> <mo>}</mo> <mo>,</mo> </mrow> <mrow> <msub> <mi>u</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;lambda;</mi> <mo>&amp;CenterDot;</mo> <mo>{</mo> <msub> <mi>&amp;lambda;u</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> </mfrac> <mo>&amp;CenterDot;</mo> <msup> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <msub> <mi>jk</mi> <mi>y</mi> </msub> <msup> <mi>e</mi> <mrow> <msub> <mi>jk</mi> <mi>y</mi> </msub> <mi>&amp;Delta;</mi> <mi>y</mi> </mrow> </msup> <mi>P</mi> <mo>&amp;rsqb;</mo> <mo>}</mo> <mo>,</mo> </mrow> λ为边界匹配参数；计算t+Δt时刻的应变梯度 <mrow> <msub> <mi>&amp;Delta;u</mi> <mi>y</mi> </msub> <mo>=</mo> <msup> <mi>F</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <msub> <mi>jk</mi> <mi>y</mi> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jk</mi> <mi>y</mi> </msub> <mi>&amp;Delta;</mi> <mi>y</mi> </mrow> </msup> <mi>F</mi> <mo>&amp;lsqb;</mo> <msub> <mi>u</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>&amp;rsqb;</mo> <mo>;</mo> </mrow>
5)计算t+Δt时刻的密度 <mrow> <msub> <mi>&amp;rho;</mi> <mi>x</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>-</mo> <msub> <mi>&amp;Delta;u</mi> <mi>x</mi> </msub> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>&amp;Delta;u</mi> <mi>x</mi> </msub> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> <mo>,</mo> </mrow> <mrow> <msub> <mi>&amp;rho;</mi> <mi>y</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;lambda;</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>-</mo> <msub> <mi>&amp;Delta;u</mi> <mi>y</mi> </msub> <msub> <mi>&amp;rho;</mi> <mn>0</mn> </msub> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <msub> <mi>&amp;Delta;u</mi> <mi>y</mi> </msub> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> <mo>,</mo> </mrow> 并得出声压对于超声探头所在位置的声压，进行回波信号加载p＝p+sz,θ (t+Δt)，进行二维傅里叶变换得到声压波数空间谱P＝F[p]；
6)循环执行步骤4)-5)，直到sz,θ(t)加载完毕，并对得到的最后声压图进行阈值分割，得到z＝z0的二维重构结果；
7)循环执行步骤2)-6)，得到所有截面的二维重构结果，采用体绘制法进行三维重构即得轴对称体缺陷的三维重构结果。
CN201610051850.3A 2016-01-26 2016-01-26 基于声场波数空间谱的轴对称体缺陷检测重构方法 Active CN105717201B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610051850.3A CN105717201B (zh) 2016-01-26 2016-01-26 基于声场波数空间谱的轴对称体缺陷检测重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610051850.3A CN105717201B (zh) 2016-01-26 2016-01-26 基于声场波数空间谱的轴对称体缺陷检测重构方法

Publications (2)

Publication Number Publication Date
CN105717201A CN105717201A (zh) 2016-06-29
CN105717201B true CN105717201B (zh) 2018-04-27

Family

ID=56154833

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610051850.3A Active CN105717201B (zh) 2016-01-26 2016-01-26 基于声场波数空间谱的轴对称体缺陷检测重构方法

Country Status (1)

Country Link
CN (1) CN105717201B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106525966B (zh) * 2016-09-22 2019-01-25 中北大学 一种轴对称工件类椭圆缺陷重构方法
KR102648455B1 (ko) * 2019-03-13 2024-03-15 제이에프이 스틸 가부시키가이샤 초음파 탐상 방법, 초음파 탐상 장치, 강재의 제조 설비, 강재의 제조 방법, 및 강재의 품질 관리 방법

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2328119A1 (de) * 1973-06-01 1974-12-05 Ronald Dr Grisell Apparat und methode fuer die abbildung der inneren struktur eines dreidimensionalen festen und/oder fluessigen objekts
US5535627A (en) * 1994-01-14 1996-07-16 The Board Of Regents Of Oklahoma State University Roll structure acoustic gage and method
US5932793A (en) * 1996-08-01 1999-08-03 Gas Research Institute Apparatus and method for determining thermophysical properties using an isochoric approach
CN1285510A (zh) * 1999-08-18 2001-02-28 中国科学院物理研究所 一种利用声波探测物体内部结构的方法
CN1307675A (zh) * 1998-06-26 2001-08-08 塞德拉公司 利用声压测量管道中的流体参数
CN101023376A (zh) * 2004-07-23 2007-08-22 比约恩·A·J·安杰尔森 超声成像
CN101726537A (zh) * 2002-04-24 2010-06-09 西德拉企业服务公司 测量管道中微粒/流体混合物的参数的装置和方法
CN102141543A (zh) * 2010-12-28 2011-08-03 天津大学 基于传声器阵列的激光焊接质量检测的方法和装置
CN102636568A (zh) * 2012-03-26 2012-08-15 湖南致力工程检测技术有限公司 一种检测混凝土内部缺陷的有限元超声成像方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015031395A1 (en) * 2013-08-26 2015-03-05 The Regents Of The University Of Colorado Imaging through scattering media with high signal to noise ratio and resolution

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2328119A1 (de) * 1973-06-01 1974-12-05 Ronald Dr Grisell Apparat und methode fuer die abbildung der inneren struktur eines dreidimensionalen festen und/oder fluessigen objekts
US5535627A (en) * 1994-01-14 1996-07-16 The Board Of Regents Of Oklahoma State University Roll structure acoustic gage and method
US5932793A (en) * 1996-08-01 1999-08-03 Gas Research Institute Apparatus and method for determining thermophysical properties using an isochoric approach
US6065328A (en) * 1996-08-01 2000-05-23 Gas Research Institute Apparatus and method for determining thermophysical properties using an isobaric approach
CN1307675A (zh) * 1998-06-26 2001-08-08 塞德拉公司 利用声压测量管道中的流体参数
CN1285510A (zh) * 1999-08-18 2001-02-28 中国科学院物理研究所 一种利用声波探测物体内部结构的方法
CN101726537A (zh) * 2002-04-24 2010-06-09 西德拉企业服务公司 测量管道中微粒/流体混合物的参数的装置和方法
CN101023376A (zh) * 2004-07-23 2007-08-22 比约恩·A·J·安杰尔森 超声成像
CN102141543A (zh) * 2010-12-28 2011-08-03 天津大学 基于传声器阵列的激光焊接质量检测的方法和装置
CN102636568A (zh) * 2012-03-26 2012-08-15 湖南致力工程检测技术有限公司 一种检测混凝土内部缺陷的有限元超声成像方法

Also Published As

Publication number Publication date
CN105717201A (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
Malyarenko et al. Ultrasonic Lamb wave diffraction tomography
CN106125047B (zh) 一种气流环境下的声源定位方法
CN106770664A (zh) 一种基于全聚焦成像算法改进边缘缺陷检测的方法
CN102608205B (zh) 基于变波速相位迁移的多层物体无损检测超声成像方法
CN104688224B (zh) 一种应用于声学非均匀媒介磁声耦合成像重建方法
CN203432946U (zh) 用于颗粒两相体系超声层析成像的超声换能器旋转装置
Han et al. Time–frequency beamforming for nondestructive evaluations of plate using ultrasonic Lamb wave
CN107356677B (zh) 基于旅行时层析成像和逆时偏移成像的超声无损检测方法
CN110261485A (zh) 一种超声波测量材料内部各处弹性模量及泊松比的方法
CN107478325A (zh) 一种非稳态平面声源的自由场还原方法
CN105717201B (zh) 基于声场波数空间谱的轴对称体缺陷检测重构方法
Guan et al. Material damage diagnosis and characterization for turbine rotors using three-dimensional adaptive ultrasonic NDE data reconstruction techniques
CN106646375A (zh) 一种基于all‑fft相位差的岩石声发射源定位方法
CN106199064A (zh) 用于管道内轴向流场成像的迭代超声层析成像方法
CN116577415A (zh) 一种构件缺陷的超声成像检测系统及方法
Rahani et al. Modeling of transient ultrasonic wave propagation using the distributed point
Dolmatov et al. Applying the algorithm of calculation in the frequency domain to ultrasonic tomography of layered inhomogeneous media using matrix antenna arrays
Jin et al. Frequency domain synthetic aperture focusing technique for variable-diameter cylindrical components
CN101819182B (zh) 重构非均匀介质中缺陷形状的方法
Dolmatov et al. Frequency-domain imaging algorithm for ultrasonic testing by application of matrix phased arrays
CN117233266A (zh) 一种基于循环神经网络的全波形反演导波层析成像方法
Dolmatov et al. Fourier-domain post-processing technique for Digital Focus Array imaging with Matrix phased array for ultrasonic testing of ITER components
CN106525966B (zh) 一种轴对称工件类椭圆缺陷重构方法
Dolmatov et al. Digital coherent signal processing with calculations in frequency domain for solving ultrasound tomography problems using matrix antenna arrays with nonequidistant arrangement of elements
Li et al. Imaging of lamination defect in metallic plate based on local wavenumber domain analysis

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Chen Youxing

Inventor after: Wang Jian

Inventor after: Zhao Xia

Inventor after: Wu Qizhou

Inventor after: Wang Zhaoba

Inventor after: Jin Yong

Inventor after: Li Haiyang

Inventor before: Chen Youxing

Inventor before: Zhao Xia

Inventor before: Wu Qizhou

Inventor before: Wang Zhaoba

Inventor before: Jin Yong

Inventor before: Li Haiyang

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant