CN104237391B - 一种相控阵超声探伤***聚焦延迟法则计算方法 - Google Patents
一种相控阵超声探伤***聚焦延迟法则计算方法 Download PDFInfo
- Publication number
- CN104237391B CN104237391B CN201410536042.7A CN201410536042A CN104237391B CN 104237391 B CN104237391 B CN 104237391B CN 201410536042 A CN201410536042 A CN 201410536042A CN 104237391 B CN104237391 B CN 104237391B
- Authority
- CN
- China
- Prior art keywords
- array element
- array
- phase
- delay
- virtual
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 238000004364 calculation method Methods 0.000 title abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims description 12
- 239000000523 sample Substances 0.000 claims description 6
- 238000000819 phase cycle Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 241000764238 Isis Species 0.000 claims 1
- 229920000535 Tan II Polymers 0.000 claims 1
- 238000007689 inspection Methods 0.000 claims 1
- 238000004611 spectroscopical analysis Methods 0.000 claims 1
- 230000001052 transient effect Effects 0.000 abstract description 4
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 238000012876 topography Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种相控阵超声探伤***聚焦延迟法则计算方法,属于超声波探伤技术领域。该方法包括以下步骤:创建虚拟阵元、计算聚焦点处相位序列、计算校正相位、计算延迟相位、计算聚焦延迟时间法则。该方法采用先计算延迟相位再换算得到聚焦延迟时间法则的方法,不同于传统的采用射线跟踪法计算延迟时间法则的方法。利用本发明的方法,可以实现对多种阵元,多种分界面形式条件下的聚焦延迟法则进行计算,相比瞬态模型具有计算效率高的优势,与射线跟踪法相比又有计算精度高的优势。
Description
技术领域
本发明属于超声波探伤技术领域,特别涉及一种基于简谐波模型的相控阵超声探伤***聚焦延迟法则计算方法。
背景技术
超声检测在无损检测技术领域具有广泛的应用,随着电子技术的高速发展及超声相控阵理论及成像算法的不断改进,超声相控阵检测技术得到了广泛的应用,超声相控阵探伤技术方法,包括:计算延迟法则、选择性激发阵元发射超声波、接收反射回波和采用适当的算法生成检测结果图像。超声相控阵探伤***根据计算得到的延迟法则,设置激发阵元的延迟时间,以达到控制声束在指定点实现能量聚集的效果,在聚焦点处声压能量最大,而在其他区域声压能量小于聚焦点处的声压能量,通过采用不同的延迟法则,实现在检测区域内多个地方进行聚焦,达到聚焦扫查的目的,这样就可以减少阵元的机械移动次数,提高超声无损检测效率。
已有的超声相控阵的计算延迟法则方法,均采用射线跟踪法进行计算,或采用瞬态声压模型的瑞利积分法进行计算。射线跟踪法是建立在假设声束是从阵元的中心点发射出来的基础之上的,并假设远处声场为球面波分布,但是对于几何形状不规则的阵元或者阵元尺寸大于波长的阵元,这种方法并不适用,会产生一定的误差。采用瞬态声压模型的瑞利积分法进行计算虽然可以获得较高精度的延迟时间,但是会大大增加计算时间,并不实用。
发明内容
本发明的目的是为克服已有技术的不足之处,针对超声相控阵探伤***,实现更高精度但又高效率地计算其聚焦延迟法则,提出一种基于简谐波模型的相控阵超声探伤***聚焦延迟法则计算方法。
本发明的技术方案如下:
本发明所述的相控阵超声探伤***聚焦延迟法则计算方法,它包括以下步骤:
1)创建虚拟阵元
设相控阵超声探头阵列共有M个阵元,M为大于等于2的整数;其中第m个阵元的参数化坐标变量为其中m=1,2,…,M,阵元m的几何中心为阵元m的局部坐标系(x’y’z’)m与全局坐标系xyz的旋转关系矩阵为阵元m的表面尺寸形状描述函数为设q=1,2,…,M-1,重复执行如下步骤:
在第q个阵元和第q+1个阵元之间创建K(q)个虚拟阵元,K(q)为大于等于2的整数,所创建的K(q)个虚拟阵元中,第k个虚拟阵元的编号为Wq(k)=k,其中k=1,2,…,K(q),第Wq(k)个虚拟阵元的参数化变量为几何中心为局部坐标系为所述的局部坐标系与全局坐标系xyz的旋转关系矩阵为第个虚拟阵元的表面尺寸形状描述函数为所述的参数化变量和为直角坐标、角度坐标或球坐标系中的一种或几种的组合,该参数化变量具有线性运算性质;所述的表面尺寸形状描述函数为矩形、圆环形、梯形、平行四边形或圆柱面形,该表面尺寸形状描述函数具有能够确定所述表面尺寸形状的大小和位置的性质;
2)计算聚焦点处相位序列
聚焦点的坐标为设q=1,2,…,M-1,重复执行如下步骤:
对第q个阵元和第q+1个阵元之间所创建的虚拟阵元,采用简谐波模型计算第Wq(k)个虚拟阵元(20)在聚焦点处声压值的复数表达式其中,为实数部分,为虚数部分,相位角
3)计算校正相位αq′
设循环累加计数器为Dq=0,对于第q个阵元和第q+1个阵元,其中q=1,2,…,M-1,执行如下步骤:
对步骤2)中计算得到的序列设Wq(k)=1,2,…,K(q)-1,循环执行如下步骤:
若且则Dq=Dq+1,若且则Dq=Dq-1,所述的γ1为大于0且小于π的实数,所述的γ2为小于0且大于-π的实数,所述的γ3为小于0且大于-π的实数,所述的γ4为大于0且小于π的实数,其中N为大于等于零的整数;
对第q个阵元,其中q=2,3,…,M,计算其校正相位αq′=αq+2π·Dq,其中Dq为相位偏差量,α1′=α1;
4)计算延迟相位Δαq
对步骤(3)得到的校正相位αq′,其中q=1,2,…,M,计算其相位延迟Δαq=max{αq′}-αq′,其中max{αq′}表示校正相位αq′的最大值,其中q=1,2,…,M;
5)计算延迟时间Δtq
对步骤4)得到的相位延迟Δαq,计算第q个阵元的时间延迟法则为:Δtq=Δαq/ω,其中ω为所述简谐波模型的角频率。
上述技术方案中,所述的相控阵超声探头阵列为一维线性阵元、一维凸型阵元、一维凹形阵元、一维柱形阵列阵元、一维圆台形阵列阵元、圆环形阵列阵元或二维平面阵列阵元。
本发明所述的简谐波模型为瑞利积分法、瑞利-索末菲尔德积分法、角谱法、多元高斯叠加法、非近轴近似多元高斯叠加法、离散点源法、弹性动力学有限积分法和有限元法中的一种或多种的组合。
本发明与现有技术相比,具有以下优点和突出性效果:本发明采用简谐波模型计算校正相位的方式,实现了超声相控阵延迟法则的计算,利用本发明的方法,可以实现对多种阵元,多种分界面形式条件下的聚焦延迟法则进行计算,采用本方法相比瞬态模型具有计算效率高的优势,本方法计算过程中考虑了阵元的几何形状,与射线跟踪法相比又有计算精度高的优势。
附图说明
图1是本发明所述的相控阵超声探伤***聚焦延迟法则计算方法的流程图。
图2是本发明的第一种实施例的一维线性阵列和圆柱面形界面示意图。
图3是图2所示实施例的阵元表面形貌示意图。
图4是图2所示实施例的虚拟阵元分布示意图。
图5是图2所示实施例的相位计算结果。
图6是图2所示实施例的校正相位计算结果。
图7是图2所示实施例的延迟时间法则计算结果。
图8是本发明的第二种实施例的一种凸型阵列示意图。
图9是图8所示实施例的阵元表面形貌示意图。
图10是图8所示实施例的虚拟阵元分布示意图。
图11是本发明的第三种实施例的一种环形阵列示意图。
图12是图11所示实施例的阵元表面形貌示意图。
图13是图12所示实施例的虚拟阵元分布示意图。
在图1至图13中:10-阵元;20-虚拟阵元。
具体实施方式
下面结合附图及实施例进一步详细说明本发明的具体工作流程。
本发明所述的相控阵超声探伤***聚焦延迟法则计算方法的一种实施例,如图1、图2、图3、图4、图5和图6所示;
1)创建虚拟阵元
设相控阵超声探头阵列共有M个阵元10,M为大于等于2的整数;其中第m个阵元的参数化坐标变量为其中m=1,2,…,M,阵元m的几何中心为阵元m的局部坐标系(x’y’z’)m与全局坐标系xyz的旋转关系矩阵为阵元m的表面尺寸形状描述函数为设q=1,2,…,M-1,重复执行如下步骤:
在第q个阵元和第q+1个阵元之间创建K(q)个虚拟阵元20,K(q)为大于等于2的整数,所创建的K(q)个虚拟阵元20中,第k个虚拟阵元20的编号为Wq(k)=k,其中k=1,2,…,K(q),第Wq(k)个虚拟阵元20的参数化变量为几何中心为局部坐标系为所述的局部坐标系与全局坐标系xyz的旋转关系矩阵为第个虚拟阵元20的表面尺寸形状描述函数为所述的参数化变量和为直角坐标、角度坐标或球坐标系中的一种或几种的组合,该参数化变量具有线性运算性质;所述的表面尺寸形状描述函数为矩形、圆环形、梯形、平行四边形或圆柱面形,该表面尺寸形状描述函数具有能够确定所述表面尺寸形状的大小和位置的性质;本实施例中,设图1中所示的一维线性相控阵超声探头阵列共有M=32个阵元10,阵列的倾斜角度为10°,探头布置在半径R=30[mm],圆心为(0,0,-30)[mm]的圆柱曲面前,介质1的密度为ρ1=1.18g/cm3,纵波声速为c1=2700m/s,介质2的密度为ρ2=7.85g/cm3,纵波声速为c2=5900m/s;相邻两个阵元之间的中心距离为0.7[mm],阵元的宽度为0.6[mm];其中第m个阵元的参数化坐标变量为(m=1,2,…,M),阵元几何中心为其中阵元局部坐标系x’y’z’m与全局坐标系xyz的旋转关系矩阵为
本实施例中,如图3所示,阵元表面尺寸形状为矩形,阵元表面尺寸描述函数为其中w=0.6[mm]为矩形的宽度,沿x’方向,h=12[mm]为矩形的长度,沿y’方向;如图3所示,本实施例中,对于序列q=1,…,31,在第q个阵元和第q+1个阵元之间创建K(q)=50个虚拟阵元20,编号为q(k)=1,…,K(q);第q(k)个虚拟阵元20的参数化变量为第q个阵元和第Wq(1)个虚拟阵元20重合,第q+1个阵元和第K(q)个虚拟阵元20重合,第Wq(k)个虚拟阵元20的几何中心为局部坐标系x’y’z’q(k)与全局坐标系xyz的旋转关系矩阵为虚拟阵元20的表面尺寸形状描述函数为
2)计算聚焦点处相位序列
聚焦点的坐标为设q=1,2,…,M-1,重复执行如下步骤:
对第q个阵元和第q+1个阵元之间所创建的虚拟阵元,采用简谐波模型计算第Wq(k)个虚拟阵元(20)在聚焦点处声压值的复数表达式其中,为实数部分,为虚数部分,相位角
本实施例中,设需要聚焦的点坐标为计算得到的相位序列如图5所示;
3)计算校正相位αq′
设循环累加计数器为Dq=0,对于第q个阵元和第q+1个阵元,其中q=1,2,…,M-1,执行如下步骤:
对步骤(2)中计算得到的序列(q(k),αq(k)),从q(k)=1到q(k)=K(q)-1循环执行如下步骤:
若找到其中满足αq(k)>γ1且αq(k)+1<γ2的元素则Dq=Dq+1,若找到满足αq(k)<γ3且αq(k)+1>γ4的元素则Dq=Dq-1,所述的γ1为大于0且小于π的实数,所述的γ2为小于0且大于-π的实数,所述的γ3为小于0且大于-π的实数,所述的γ4为大于0且小于π的实数,其中N为大于等于零的整数;
对第q个阵元,其中q=2,3,…,M,计算其校正相位αq′=αq+2π·Dq,其中Dq为相位偏差量,α1′=α1;
本实施例中,γ1=4π/5,γ2=-4π/5,γ3=-4π/5,γ4=4π/5;本实施例的校正相位计算结果如图6所示;
4)计算延迟相位Δαq
对步骤(3)得到的校正相位αq′,其中q=1,2,…,M,计算其相位延迟Δαq=max{αq′}-αq′,其中max{αq′}表示校正相位αq′的最大值,其中q=1,2,…,M;
5)计算延迟时间Δtq
对步骤(4)得到的相位延迟Δαq,计算第q个阵元的时间延迟法则为:Δtq=Δαq/ω。本实施例的延迟时间法则计算结果如图7所示。
本发明所述的另一种实施例的创建虚拟阵元20流程,如图8和图9所示,如图8所示的一种凸型超声相控阵阵列共有M=16个阵元10,阵列的圆弧中心坐标为C=(0,0,0)[mm],圆弧半径为R=30[mm],相邻两个阵元之间的夹角为5.625°,阵元的宽度为1[mm];其中第m个阵元的参数化坐标变量为(m=1,…,M),阵元几何中心为阵元局部坐标系x’y’z’m与全局坐标系xyz的旋转关系矩阵为
本实施例中,如图3所示,与前一个实施例相同,阵元表面尺寸形状为矩形,阵元表面尺寸描述函数为其中w=1.5[mm]为矩形的宽度,沿x’方向,h=12[mm]为矩形的长度,沿y’方向;如图9所示,本实施例中,对于序列q=1,…,16,在第q个阵元和第q+1个阵元之间创建K(q)=50个虚拟阵元20,编号为Wq(k)=1,…,K(q);第Wq(k)个虚拟阵元20的参数化变量为第q个阵元和第Wq(1)个虚拟阵元20重合,第q+1个阵元和第K(q)个虚拟阵元20重合,第Wq(k)个虚拟阵元20的几何中心为局部坐标系x’y’z’q(k)与全局坐标系xyz的旋转关系矩阵为虚拟阵元20的表面尺寸形状描述函数为
本发明所示的第三种实施例创建虚拟阵元20的流程,如图10和图11所示,如图8所示的一种圆环形超声相控阵阵列共有M=10个阵元,阵列的中心坐标为C=(0,0,0)[mm],其中第m个阵元的参数化坐标变量为sm=(m+2)[mm],(m=1,…,M),该参数化坐标变量表示第m个阵元的中心半径,
阵元局部坐标系(x’y’z’)m与全局坐标系xyz的旋转关系矩阵为
本实施例中,如图11所示,阵元表面尺寸形状为圆环形,阵元表面尺寸描述函数为
其中Rm1=sm+0.45[mm]为圆环的外径函数,Rm2=sm-0.45[mm]为圆环的内径函数;如图12所示,本实施例中,对于序列q=1,…,16,在第q个阵元和第q+1个阵元之间创建K(q)=11个虚拟阵元20,编号为Wq(k)=1,…,K(q);第q(k)个虚拟阵元20的参数化变量为第q个阵元和第Wq(1)个虚拟阵元20重合,第q+1个阵元和第K(q)个虚拟阵元20重合,第Wq(k)个虚拟阵元20的几何中心为局部坐标系(x’y’z’)q(k)与全局坐标系xyz的旋转关系矩阵为虚拟阵元20的表面尺寸形状描述函数为
本发明的关键点在于首先在两个相邻阵元之间创建虚拟探头,然后采用简谐波模型计算相位的计算方法,获得两个相邻阵元之间的相位变化过程,从而对实际阵元计算得到的相位进行校正,得到校正相位,利用本发明的方法,可以实现对多种阵元,多种分界面形式条件下的聚焦延迟法则进行计算,采用本方法相比瞬态模型具有计算效率高的优势,本方法计算过程中考虑了阵元的几何形状,与射线跟踪法相比又有计算精度高的优势。
Claims (3)
1.一种相控阵超声探伤***聚焦延迟法则计算方法,其特征在于:它包括以下步骤:
1)创建虚拟阵元
设相控阵超声探头阵列共有M个阵元(10),M为大于等于2的整数;其中第m个阵元的参数化坐标变量为其中m=1,2,…,M,阵元m的几何中心为阵元m的局部坐标系(x’y’z’)m与全局坐标系xyz的旋转关系矩阵为阵元m的表面尺寸形状描述函数为设q=1,2,…,M-1,重复执行如下步骤:
在第q个阵元和第q+1个阵元之间创建K(q)个虚拟阵元(20),K(q)为大于等于2的整数,所创建的K(q)个虚拟阵元(20)中,第k个虚拟阵元(20)的编号为Wq(k)=k,其中k=1,2,…,K(q),第Wq(k)个虚拟阵元(20)的参数化变量为几何中心为局部坐标系为所述的局部坐标系与全局坐标系xyz的旋转关系矩阵为第个虚拟阵元(20)的表面尺寸形状描述函数为所述的参数化变量和为直角坐标、角度坐标或球坐标系中的一种或几种的组合,该参数化变量具有线性运算性质;所述的表面尺寸形状描述函数为矩形、圆环形、梯形、平行四边形或圆柱面形,该表面尺寸形状描述函数具有能够确定所述表面尺寸形状的大小和位置的性质;
2)计算聚焦点处相位序列
设聚焦点的坐标为设q=1,2,…,M-1,重复执行如下步骤:
对第q个阵元和第q+1个阵元之间所创建的虚拟阵元(20),采用简谐波模型计算第Wq(k)个虚拟阵元(20)在聚焦点处声压值的复数表达式其中,为实数部分,为虚数部分,相位角其中式为:
3)计算校正相位αq′
设循环累加计数器为Dq=0,对于第q个阵元和第q+1个阵元,其中q=1,2,…,M-1,执行如下步骤:
对步骤2)中计算得到的序列设Wq(k)=1,2,…,K(q)-1,循环执行如下步骤:
若且则Dq=Dq+1,若且的则Dq=Dq-1,所述的γ1为大于0且小于π的实数,所述的γ2为小于0且大于-π的实数,所述的γ3为小于0且大于-π的实数,所述的γ4为大于0且小于π的实数,其中N为大于等于零的整数;
对第q个阵元,其中q=2,3,…,M,计算其校正相位αq′=αq+2π·DΣ,其中Dq为相位偏差量,α1′=α1;
4)计算延迟相位Δαq
对步骤3)得到的校正相位αq′,其中q=1,2,…,M,计算其相位延迟Δαq=max{αq′}-αq′,其中max{αq′}表示校正相位αq′的最大值,其中q=1,2,…,M;
5)计算延迟时间Δtq
对步骤4)得到的相位延迟Δαq,计算第q个阵元的时间延迟法则为:Δtq=Δαq/ω,其中ω为所述简谐波模型的角频率。
2.根据权利要求1所述的相控阵超声探伤***聚焦延迟法则计算方法,其特征在于,相控阵超声探头阵列为一维线性阵元、一维凸型阵元、一维凹形阵元、一维柱形阵列阵元、一维圆台形阵列阵元、圆环形阵列阵元或二维平面阵列阵元。
3.根据权利要求1所述的相控阵超声探伤***聚焦延迟法则计算方法,其特征在于,所述的简谐波模型为瑞利积分法、瑞利-索末菲尔德积分法、角谱法、多元高斯叠加法、非近轴近似多元高斯叠加法、离散点源法、弹性动力学有限积分法和有限元法中的一种或多种的组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410536042.7A CN104237391B (zh) | 2014-10-11 | 2014-10-11 | 一种相控阵超声探伤***聚焦延迟法则计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410536042.7A CN104237391B (zh) | 2014-10-11 | 2014-10-11 | 一种相控阵超声探伤***聚焦延迟法则计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104237391A CN104237391A (zh) | 2014-12-24 |
CN104237391B true CN104237391B (zh) | 2017-05-17 |
Family
ID=52225886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410536042.7A Active CN104237391B (zh) | 2014-10-11 | 2014-10-11 | 一种相控阵超声探伤***聚焦延迟法则计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104237391B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107362466B (zh) * | 2017-06-29 | 2019-12-10 | 哈尔滨医科大学 | 一种用于声动力治疗的超声聚焦方法 |
CN107727742B (zh) * | 2017-10-11 | 2020-07-31 | 哈尔滨工业大学 | 一种电磁超声相控阵*** |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0006197B1 (de) * | 1978-06-16 | 1983-01-12 | F. HOFFMANN-LA ROCHE & CO. Aktiengesellschaft | Ultraschall-Bildgerät und dessen Verwendung für die medizinische Ultraschall-Diagnostik |
FR2930642B1 (fr) * | 2008-04-23 | 2010-06-11 | Commissariat Energie Atomique | Traducteur ultrasonore de contact, a elements multiples, sabot flexible et profilometre |
CN101403729A (zh) * | 2008-04-30 | 2009-04-08 | 硕德(北京)科技有限公司 | 相控阵超声检测中基于集总相关法的相位畸变校正方法 |
CN101403730A (zh) * | 2008-04-30 | 2009-04-08 | 硕德(北京)科技有限公司 | 相控阵超声检测中基于联合相关法的相位畸变校正方法 |
CN102967659B (zh) * | 2012-10-30 | 2014-11-26 | 广东电网公司电力科学研究院 | 相控阵超声探头在多层介质中探伤时声场分布的计算方法 |
CN103412045B (zh) * | 2013-07-31 | 2015-06-10 | 北京时代之峰科技有限公司 | 一种具有楔块的超声波检测中延迟量计算方法及装置 |
-
2014
- 2014-10-11 CN CN201410536042.7A patent/CN104237391B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104237391A (zh) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8955384B2 (en) | Test probe as well as family of test probes for the non-destructive testing of a workpiece by means of ultrasonic sound and testing device | |
CN102967659B (zh) | 相控阵超声探头在多层介质中探伤时声场分布的计算方法 | |
Zernov et al. | A refinement of the Kirchhoff approximation to the scattered elastic fields | |
CN104237391B (zh) | 一种相控阵超声探伤***聚焦延迟法则计算方法 | |
US9117038B2 (en) | System and method for optimizing the design of an ultrasonic transducer | |
Franco et al. | Acoustic beam modeling of ultrasonic transducers and arrays using the impulse response and the discrete representation methods | |
US11480550B2 (en) | Method and device for mapping components for detecting elongation direction | |
Höhne et al. | SAFT imaging of transverse cracks in austenitic and dissimilar welds | |
Cheng et al. | A new algorithm for spatial impulse response of rectangular planar transducers | |
Rabe et al. | Application of the total focusing method for quantitative nondestructive testing of anisotropic welds with ultrasound | |
Wise et al. | Staircase-free acoustic sources for grid-based models of wave propagation | |
Zhong et al. | A modified convolution model for calculating the far field directivity of a parametric array loudspeaker | |
Beissner | On the lateral resolution of focused ultrasonic fields from spherically curved transducers | |
Guo et al. | Simulation of acoustic fields emitted by ultrasonic phased array in austenitic steel weld | |
Katch et al. | Focal depth localization for highly focused transducers in isotropic materials | |
Jeong et al. | Transmission Phase Control of Annular Array Transducers for Efficient Second Harmonic Generation in the Presence of a Stress-Free Boundary | |
Danilov et al. | Simulation of the Operation of a Normal Transducer with a Phased Array in the Transmission Mode | |
Chekroun et al. | Measurement of dispersion curves of circumferential guided waves radiating from curved shells: Theory and numerical validation | |
Wright et al. | Prediction of acoustic radiation from axisymmetric surfaces with arbitrary boundary conditions using the boundary element method on a distributed computing system | |
Liaptsis et al. | Focal law calculations for annular phased array transducers | |
Li et al. | Parameter optimization of linear phased array transducer for defect detection | |
Roberts | Model development for ultrasound inspection of anisotropic materials | |
Lupien et al. | Improved titanium billet inspection sensitivity through optimized phased array design, Part I: Design technique, modeling and simulation | |
US20240094169A1 (en) | Determination of geometric dimensions | |
Shi et al. | Simulation of ultrasonic phased array field in weld bead |
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 |