CN110794033B - 一种精准控制声场幅值与偏振方向的多波聚焦方法 - Google Patents

一种精准控制声场幅值与偏振方向的多波聚焦方法 Download PDF

Info

Publication number
CN110794033B
CN110794033B CN201911110028.XA CN201911110028A CN110794033B CN 110794033 B CN110794033 B CN 110794033B CN 201911110028 A CN201911110028 A CN 201911110028A CN 110794033 B CN110794033 B CN 110794033B
Authority
CN
China
Prior art keywords
wave
waves
amplitude
array element
excitation
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
CN201911110028.XA
Other languages
English (en)
Other versions
CN110794033A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201911110028.XA priority Critical patent/CN110794033B/zh
Publication of CN110794033A publication Critical patent/CN110794033A/zh
Application granted granted Critical
Publication of CN110794033B publication Critical patent/CN110794033B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/04Analysing solids
    • 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/22Details, e.g. general constructional or apparatus details
    • G01N29/221Arrangements for directing or focusing the acoustical waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids

Abstract

一种精准控制声场幅值与偏振方向的多波聚焦方法,属于超声无损检测领域,该方法根据介质已知横波声速cs与纵波声速cl建立声场指向性函数来描述介质内任意一点的声场,对介质表面一阵元施加两个激励脉冲,两脉冲的时间差由目标焦点的声程与介质声速确定,两脉冲的激励幅值由声场指向性函数与目标焦点处的聚焦波幅值和偏振方向角γ确定。通过改变脉冲的激励时间差t0、幅值,可实现目标焦点处横波与纵波的同时聚焦,对介质内任意点声场幅值与偏振方向进行精准控制。本方法无需预模拟,单阵元即可多波聚焦,克服了利用单一波聚焦的缺点,实现较大扫描区域聚焦同时可进行多偏振方向扫描,有效提高信噪比和分辨率,对实际工程运用具有意义。

Description

一种精准控制声场幅值与偏振方向的多波聚焦方法
技术领域
本发明涉及一种精准控制声场幅值与偏振方向的多波聚焦方法,其属于超声无损检测领域。
背景技术
在超声检测中常使用聚焦声束进行检测,使用聚焦声束会使声能更为集中,中心轴线上的声压增强,同时可改善声束指向性,对提高探伤灵敏度、分辨力和信噪比均为有利。在常规的超声聚焦检测中,人们普遍采用单一类型波包的速度计算时间延迟来进行聚焦发射和聚焦接收处理。这种聚焦方式只能使一种类型的波包聚焦,导致其检测方向只限于该类型波的主传播方向。超声多波聚焦是一种利用时间反转自适应聚焦特性使多种成分的波包同时聚焦的超声检测方法。声波在固体介质中传播时,存在纵波、横波等不同类型的波,同时利用介质中的横波与纵波进行多波聚焦时,不仅能得到更加丰富的介质和缺陷信息,而且可实现大范围聚焦。另一方面在多波聚焦的基础上,通过改变纵、横波的激发强度,可以在一定角度范围内,任意控制目标焦点处的声偏振方向和幅值,形成多个偏振方向的扫描,实现对缺陷形状和方向的多角度扫描和检测。
该方法属于一个比较新颖的思想,目前,在多波聚焦方面张碧星等人基于数值模拟,在预设聚焦点处施加预设方向振动的激励位移,得到应加载在激励阵元上的激励信号,将该信号加载到激励阵元上激发,在预设聚焦点处实现多波聚焦。该方法的偏振方向的控制需要通过预模拟来得到激励信号,每改变一次偏振方向都需要进行一次预模拟。如果实验中的介质与数值模拟中的介质参数存在偏差,则用该方法进行实际实验时会存在误差。有人提出了一种质点偏振方向控制与扫描检测的方法,该方法的特点是使用了两个换能器,直接通过矢量合成的方法计算预设聚焦点处的声偏振方向,然后控制两个换能器的时间延迟及激励脉冲进行声偏振方向的控制,与前一个方法对比不需要预模拟。但是该方法对换能器的位置十分敏感,当换能器位置与预设位置出现极小偏差时,声偏振方向将偏离预设方向或不再是线偏振,便可导致聚焦效果变差。
发明内容
本方法的目的是提供一种精准控制声场幅值与偏振方向的多波聚焦方法。利用横波和纵波进行多波聚焦,实现了对聚焦波的声偏振方向和幅值的精准控制。与以上提到的两种方法相比,本方法从原理角度,通过公式计算声偏振方向并加以控制,无需预模拟;并且利用单阵元探头进行检测,避免了阵元位置对聚焦效果的影响。本方法通过改变脉冲的激励时间差t0与激励幅值,可实现目标焦点处横波与纵波的同时聚焦,对介质内任意点声场幅值与偏振方向进行精准控制。可通过控制声偏振实现对多个偏振方向的聚焦扫描,并在较大的扫描区域内都能获得很好的聚焦效果,有效提高信噪比和分辨率,对实际工程运用具有意义。
一种精准控制声场幅值与偏振方向的多波聚焦方法包括以下步骤:
(1)建立声场指向性理论公式,激励阵元宽度为2a,θ为目标焦点相对原点的方位角,即tanθ=x/y,离阵元中心距离为r一点(x,y)处纵波和横波的位移势为:
Figure GDA0002943097750000021
Figure GDA0002943097750000022
其中
Figure GDA0002943097750000031
kp和ks分别为纵波和横波的波数,Φp和Φs分别为纵波和横波的位移势,P(θ)和S(θ)为纵波和横波的指向性因子,其表达式如(3)与(4)所示,a为换能器阵元在x方向上尺寸的一半;
Figure GDA0002943097750000032
Figure GDA0002943097750000033
Figure GDA0002943097750000034
其中ξ在公式(3)与(4)中分别为kp sinθ与ks sinθ;
(2)已知介质其横波与纵波声速cs与cl,用于阵元发射的两个脉冲的激励时间差t0=r/cs-r/cl,阵元发射的激励信号中第一个脉冲表示为S波,第二个脉冲表示为P波;
则S波的位移函数为:
Figure GDA0002943097750000035
P波的位移函数为:
Figure GDA0002943097750000036
其中,D为一个脉冲的持续时间,qs为S波的放大倍数,qp为P波的放大倍数;
(3)在固体介质中,这两个波会分别激发出对应的横波和纵波,将S波激发的横波命名为S-S波,将P波激发的纵波命名为P-P波,则根据时间反转法的原理,在设置的激励时间差t0条件下,S-S波与P-P波会在目标焦点(x,y)处聚焦合成,根据步骤(1)中的公式推导得:
聚焦点处S-S波在X方向的位移函数为:
Figure GDA0002943097750000037
在Y方向的位移函数为:
Figure GDA0002943097750000041
P-P波在X方向的位移函数为:
Figure GDA0002943097750000042
在Y方向的位移函数为:
Figure GDA0002943097750000043
其中,ts=r/cs、tp=r/cl,在聚焦点处,有t0+tp=ts的关系,此时两个脉冲的位移函数相位相同;
则聚焦波在X方向的位移函数为:
Figure GDA0002943097750000044
在Y方向的位移函数为:
Figure GDA0002943097750000045
则目标焦点处聚焦波的偏振角为γ,单位为度,偏振角的正切值即为目标焦点处X方向位移值与Y方向的位移值之比;目标焦点处的聚焦波幅值为X方向位移值与Y方向位移值的矢量叠加之和;
在步骤(2)与(3)的基础上,根据公式(12)与(13),通过不断改变qs和qp,可对介质内任意点实现声偏振方向和幅值的精准控制;针对目标焦点处的偏振方向和幅值,根据公式推导出对应的qs和qp,可得到阵元应激励脉冲的波形;发射该激励脉冲,则在目标焦点处S-S波与P-P波相位一致,实现聚焦。
本发明的效果和益处是:这种精准控制声场幅值与偏振方向的多波聚焦方法,可任意控制聚焦波的幅值和声偏振方向,实现介质内任意位置的横波与纵波聚焦。该方法具有聚焦检测中声能集中、指向性好、灵敏度和分辨率高的优点,可通过控制声偏振实现对多个偏振方向的聚焦扫描,并在较大的扫描区域内都能获得很好的聚焦效果,在超声检测中的缺陷取向识别方面具有实用意义。
附图说明
图1为具体实施方式中数值模拟的模型图。
图2为应施加在阵元上的激励脉冲图。
图3为探针接收到的位移结果图。
图4为根据结果作出的声偏振极图。
图5为改变条件下的声偏振极图。
具体实施方式
(1)在介质表面布置一单阵元换能器,阵元宽度为1mm,以阵元中心作为(0,0)点,目标焦点的位置为(-30mm,-30mm),模型图如图1所示,计算出从换能器中心到目标焦点的声程r=42.426mm,材料介质为钢,cs=3140m/s、cl=5770m/s,利用t0=r/cs-r/cl,计算出用于阵元发射的两个脉冲的激励时间差to=6.159μs,目标焦点处预设幅值为0.00890μm,偏振方向为110°。
(2)根据声场指向性函数计算对应的P(θ)和S(θ),然后根据发明内容中的公式(12)与(13),由给定的偏振角和入射角推导出qs/qp=0.886,则在设定的幅值条件下,计算得qs=0.828,qp=0.935。由此得到应施加在阵元上的激励脉冲的波形,如图2所示。
(3)使阵元激励该脉冲,则在目标焦点处理论上接收波包的顺序是:S-P波、S-S波与P-P波的聚焦波、P-S波,在目标焦点处设置两个探针,分别用于接收X方向和Y方向的位移,探针所接收到的整体波形如图3所示,与理论结果相符,对聚焦波进行偏振分析,作出的偏振极图如图4所示,偏振角为109.68°,聚焦波幅值为0.00885μm。与理论值相比,偏振角度误差为0.3%,偏振幅值误差为0.6%。
(4)修改偏振角度为120°,重新计算出qs=0.883,qp=0.573。重新模拟的偏振极图如图5所示,偏振角为119.78°,聚焦波幅值为0.00886μm,误差均在允许范围内。
(5)根据以上步骤,通过不断改变qs与qp的值,可对聚焦点处进行多偏振方向的扫描,将计算好的S波和P波激发条件输入阵元,可以在目标焦点处接收到S-S波和P-P波,此时两个波实现聚焦,聚焦波的质点偏振方向和幅值均与理论值相同,实现了声偏振方向和幅值的双重控制,对超声聚焦检测的实际运用具有一定的价值。

Claims (1)

1.一种精准控制声场幅值与偏振方向的多波聚焦方法,其特征在于,包括以下步骤:
(1)建立声场指向性理论公式,激励阵元宽度为2a,θ为目标焦点相对原点的方位角,即tanθ=x/y,离阵元中心距离为r一点(x,y)处纵波和横波的位移势为:
Figure FDA0002950800880000011
Figure FDA0002950800880000012
其中
Figure FDA0002950800880000013
kp和ks分别为纵波和横波的波数,Φp和Φs分别为纵波和横波的位移势,P(θ)和S(θ)为纵波和横波的指向性因子,其表达式如(3)与(4)所示,a为换能器阵元在x方向上尺寸的一半;
Figure FDA0002950800880000014
Figure FDA0002950800880000015
Figure FDA0002950800880000016
其中ξ在公式(3)与(4)中分别为kpsinθ与kssinθ;
(2)已知介质其横波与纵波声速cs与cl,用于阵元发射的两个脉冲的激励时间差t0=r/cs-r/cl,阵元发射的激励信号中第一个脉冲表示为S波,第二个脉冲表示为P波;
则S波的位移函数为:
Figure FDA0002950800880000017
P波的位移函数为:
Figure FDA0002950800880000021
其中,D为一个脉冲的持续时间,qs为S波的放大倍数,qp为P波的放大倍数;
(3)在固体介质中,这两个波会分别激发出对应的横波和纵波,将S波激发的横波命名为S-S波,将P波激发的纵波命名为P-P波,则根据时间反转法的原理,在设置的激励时间差t0条件下,S-S波与P-P波会在目标焦点(x,y)处聚焦合成,根据步骤(1)中的公式推导得:
聚焦点处S-S波在X方向的位移函数为:
Figure FDA0002950800880000022
在Y方向的位移函数为:
Figure FDA0002950800880000023
P-P波在X方向的位移函数为:
Figure FDA0002950800880000024
在Y方向的位移函数为:
Figure FDA0002950800880000025
其中,ts=r/cs、tp=r/cl,在聚焦点处,有t0+tp=ts的关系,此时两个脉冲的位移函数相位相同;
则聚焦波在X方向的位移函数为:
Figure FDA0002950800880000026
在Y方向的位移函数为:
Figure FDA0002950800880000031
则目标焦点处聚焦波的偏振角为γ,单位为度,偏振角的正切值即为目标焦点处X方向位移值与Y方向的位移值之比;目标焦点处的聚焦波幅值为X方向位移值与Y方向位移值的矢量叠加之和;
(4)在步骤(2)与(3)的基础上,根据公式(12)与(13),通过不断改变qs和qp,对介质内任意点实现声偏振方向和幅值的精准控制;针对目标焦点处的偏振方向和幅值,根据公式推导出对应的qs和qp,得到阵元应激励脉冲的波形;发射该激励脉冲,则在目标焦点处S-S波与P-P波相位一致,实现聚焦。
CN201911110028.XA 2019-11-13 2019-11-13 一种精准控制声场幅值与偏振方向的多波聚焦方法 Active CN110794033B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911110028.XA CN110794033B (zh) 2019-11-13 2019-11-13 一种精准控制声场幅值与偏振方向的多波聚焦方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911110028.XA CN110794033B (zh) 2019-11-13 2019-11-13 一种精准控制声场幅值与偏振方向的多波聚焦方法

Publications (2)

Publication Number Publication Date
CN110794033A CN110794033A (zh) 2020-02-14
CN110794033B true CN110794033B (zh) 2021-05-04

Family

ID=69444508

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911110028.XA Active CN110794033B (zh) 2019-11-13 2019-11-13 一种精准控制声场幅值与偏振方向的多波聚焦方法

Country Status (1)

Country Link
CN (1) CN110794033B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049092B (zh) * 2021-03-26 2022-10-11 中北大学 约束条件下超声阵列辐射器辐射声场的计算方法及系统

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1694644A (zh) * 2002-06-28 2005-11-09 Oti眼技术股份有限公司 可调的深度分辨率及多功能光学映射装置
EP2422708A1 (en) * 2010-08-31 2012-02-29 General Electric Company Multi-focus ultrasound system and method
CN104729435A (zh) * 2015-03-13 2015-06-24 浙江大学 基于激光超声的板材厚度在线测量系统
CN105319271A (zh) * 2014-07-30 2016-02-10 中国科学院声学研究所 一种横、纵波联合超声相控阵检测方法
CN106990172A (zh) * 2017-03-31 2017-07-28 中国科学院声学研究所 一种超声检测与成像方法及其系统
CN107085041A (zh) * 2017-03-31 2017-08-22 中国科学院声学研究所 一种聚焦声场中声偏振方向的控制方法
KR20170095550A (ko) * 2016-02-15 2017-08-23 이선우 다기능 고강도 집속형 초음파 방출 카트리지 및 핸드피스
CN107205729A (zh) * 2015-02-18 2017-09-26 株式会社日立制作所 超声波成像装置以及超声波信号的处理方法
CN107860716A (zh) * 2017-10-30 2018-03-30 东北大学 一种基于激光超声的弹性常数的无损检测方法及设备
CN108645916A (zh) * 2018-04-23 2018-10-12 中国科学院声学研究所 一种基于多波聚焦的质点偏振方向控制新方法
CN109142537A (zh) * 2018-08-20 2019-01-04 中国科学院声学研究所 一种质点偏振方向控制与扫描检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013176045A1 (ja) * 2012-05-25 2013-11-28 富士フイルム株式会社 超音波信号処理装置および超音波信号処理方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1694644A (zh) * 2002-06-28 2005-11-09 Oti眼技术股份有限公司 可调的深度分辨率及多功能光学映射装置
EP2422708A1 (en) * 2010-08-31 2012-02-29 General Electric Company Multi-focus ultrasound system and method
CN105319271A (zh) * 2014-07-30 2016-02-10 中国科学院声学研究所 一种横、纵波联合超声相控阵检测方法
CN107205729A (zh) * 2015-02-18 2017-09-26 株式会社日立制作所 超声波成像装置以及超声波信号的处理方法
CN104729435A (zh) * 2015-03-13 2015-06-24 浙江大学 基于激光超声的板材厚度在线测量系统
KR20170095550A (ko) * 2016-02-15 2017-08-23 이선우 다기능 고강도 집속형 초음파 방출 카트리지 및 핸드피스
CN106990172A (zh) * 2017-03-31 2017-07-28 中国科学院声学研究所 一种超声检测与成像方法及其系统
CN107085041A (zh) * 2017-03-31 2017-08-22 中国科学院声学研究所 一种聚焦声场中声偏振方向的控制方法
CN107860716A (zh) * 2017-10-30 2018-03-30 东北大学 一种基于激光超声的弹性常数的无损检测方法及设备
CN108645916A (zh) * 2018-04-23 2018-10-12 中国科学院声学研究所 一种基于多波聚焦的质点偏振方向控制新方法
CN109142537A (zh) * 2018-08-20 2019-01-04 中国科学院声学研究所 一种质点偏振方向控制与扫描检测方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
Multi-mode ultrasonic waves focusing in a variable focus technique for simultaneous sound-velocity and thickness measurement;Xiaoyu Yang et al.;《Applied Acoustics》;20191022;第1-7页 *
Multi-Wave and Hybrid Imaging Techniques: A New Direction for Nondestructive Testing and Structural Health Monitoring;Yuhua Cheng et al.;《Sensors》;20131127;第16146-16189页 *
RESEARCH ON DETECTION OF WELD DEFECTS WITH MULTI-WAVE TOTAL FOCUSING METHOD;ZHANG Ping et al.;《2019 14th Symposium on Piezoelectrcity,Acoustic Waves and Device Applications》;20191104;第1-5页 *
Researches on the Ultrasonic Scattering Attenuation of Carbon Fibre Reinforced Plastics with 2D Real Morphology Void Model;DING et al.;《Acoustical Physics》;20171231;第63卷(第4期);第490-495页 *
Ultrasonic focusing and scanning with multiple waves;Zhang Bi-Xing et al.;《Chin. Phys. B》;20131231;第22卷(第1期);第014302-1-7页 *
光弹法测量超声换能器声场;金士杰等;《应用声学》;20140331;第33卷(第2期);第107-111页 *
凹面相控阵聚焦声场在液固界面上的反射和折射;张碧星,王文龙;《物理学报》;20080630;第57卷(第6期);第3613-3619页 *
激光超声横波声场指向性的研究;刘品潇等;《激光杂质》;20151231;第36卷(第9期);第40-43页 *
超声多波聚焦及声偏振方向控制方法;阎守国,张碧星;《应用声学》;20180131;第37卷(第1期);第63-69页 *
超声多波聚焦扫描方法研究;刘冬冬,张碧星;《声学技术》;20121031;第31卷(第5期);第31-32页 *

Also Published As

Publication number Publication date
CN110794033A (zh) 2020-02-14

Similar Documents

Publication Publication Date Title
Senesi et al. Experimental characterization of periodic frequency-steerable arrays for structural health monitoring
Drinkwater et al. Ultrasonic arrays for non-destructive evaluation: A review
Clay et al. Experimental study of phased array beam steering characteristics
Engholm et al. Adaptive beamforming for array imaging of plate structures using lamb waves
US20040066708A1 (en) Ultrasonic transmitting and receiving apparatus
Senesi et al. A frequency selective acoustic transducer for directional Lamb wave sensing
Trots et al. Synthetic aperture method in ultrasound imaging
CN108226304A (zh) 一种基于测量模型的超声相控阵线扫描灵敏度计算方法
CN110794033B (zh) 一种精准控制声场幅值与偏振方向的多波聚焦方法
CN110988850A (zh) 一种基于目标散射的换能器指向性测量方法及装置
Camacho et al. Auto-focused virtual source imaging with arbitrarily shaped interfaces
JP3635453B2 (ja) 超音波横波斜角探傷方法及び装置
Dixon et al. The wave-field from an array of periodic emitters driven simultaneously by a broadband pulse
Trots et al. Golay coded sequences in synthetic aperture imaging systems
JP3606132B2 (ja) 超音波探傷方法およびその装置
CN112684005A (zh) 基于二维矩阵换能器的全聚焦检测方法
Allevato et al. Spiral air-coupled ultrasonic phased array for high resolution 3D imaging
Rodrigues et al. Development of a 2-d array ultrasonic transducer for 3-d imaging of objects immersed in water
Li et al. Adaptive array processing for ultrasonic non-destructive evaluation
Stepinski et al. Designing 2D arrays for SHM of planar structures: a review
Li et al. A rapid approach to speckle noise reduction in ultrasonic non-destructive evaluation using matched filters
Berktay et al. Arrays of parametric receiving arrays
Kubrusly et al. Measurement of plate thickness based on the optimal unidirectional generation of ultrasonic waves using dual electromagnetic acoustic transducer
Engholm et al. Using 2-D arrays for sensing multimodal lamb waves
Xu et al. Ultrasonic array inspection for additive manufacturing components using full matrix capture

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