相控聚焦超聲與圓孔交互作用的散射特性研究

周進節(jié)，鄭陽，張宗健

?

相控聚焦超聲與圓孔交互作用的散射特性研究

周進節(jié)1，鄭陽2，張宗健2

(1. 中北大學(xué)機械與動力工程學(xué)院，山西太原030051；2. 中國特種設(shè)備檢測研究院，國家質(zhì)量監(jiān)督檢驗檢疫總局無損檢測與評價重點實驗室，北京100029)

相較于傳統(tǒng)超聲檢測，相控陣技術(shù)的聲束聚焦控制能力使得相控超聲檢測的分辨率、信噪比和靈敏度等性能得到有效提升。但針對相控超聲與缺陷的交互作用認識不足，導(dǎo)致目前檢測技術(shù)無法對缺陷進行準確定性和精確定量。為此，建立了一維線陣相控超聲檢測圓孔缺陷的有限元模型，重點研究相控聚焦波束與圓孔缺陷交互作用的聲場散射特性。結(jié)果表明，圓孔散射L模態(tài)的能量主要集中于正反射區(qū)和正透射區(qū)，且在透射區(qū)的能量大于反射區(qū)的能量，而在與入射聲束垂直的方向上能量極其微弱，幾乎沒有散射能量；散射L模態(tài)沿入射方向呈軸對稱分布；對于不同直徑的圓孔缺陷，直徑越大，在相同方向上散射波的能量越大；對于相同直徑的圓孔，聚焦波與平面波入射時，散射波的分布規(guī)律相同，但聚焦波散射的能量明顯高于平面波入射時散射的能量。

相控；聚焦超聲；圓孔；散射特性；散射矩陣

0 引言

相控陣超聲以其便于聲束控制、檢測范圍大、結(jié)果直觀等特點成為超聲無損檢測領(lǐng)域極具應(yīng)用前景的技術(shù)之一[1-3]。隨著電子技術(shù)的迅猛發(fā)展，商用超聲相控陣的硬件電路已能支持256陣元同時激勵和接收[4]。近年來，許多學(xué)者致力于研究超聲相控陣成像算法，旨在利用相控陣檢測信號提高成像質(zhì)量，代表性的算法有全聚焦算法[5]、逆波場外推算法[6]和波數(shù)空間算法[7]等。但這些成像算法將各種缺陷均簡化為完美散射體，認為所有缺陷散射聲場在各方向分布規(guī)律相同，從而導(dǎo)致目前算法成像結(jié)果中圓孔缺陷顯示為一片扇形區(qū)間，無法成像出圓孔缺陷的實際形狀，更不能對圓孔進行精確的定量分析。而實際上，缺陷并非完美散射體，缺陷的散射信息是隨入射波種類、方向，以及缺陷形貌、尺寸變化的[8-10]。

超聲相控陣是通過控制各通道延遲參數(shù)來實現(xiàn)聚焦檢測的，為實現(xiàn)超聲相控陣高分辨率成像以對缺陷進行準確定性、精確定量，研究了相控陣合成聚焦聲束與缺陷交互作用規(guī)律。聚焦聲束與缺陷的交互作用主要體現(xiàn)在交互作用后的散射場分布上，從而可通過研究缺陷散射場分布來研究聚焦聲束與缺陷的交互作用。

1 仿真模型的建立

建立一維線性相控陣超聲與孔缺陷交互作用的二維有限元模型，如圖1所示，在待檢試件表面設(shè)置多個等間距分布的壓強載荷來模擬相控陣傳感器激勵產(chǎn)生超聲波，線性相控陣陣元數(shù)為32個，單個陣元寬度為0.5 mm，陣元中心間距為0.6 mm。待檢試件為鋼材，材料參數(shù)為密度7.85×103 kg/m3、楊氏模量2.1×105 MPa、泊松比0.26，幾何參數(shù)為長50 mm、寬40 mm。在研究圓孔缺陷的散射特性時，當(dāng)試件待檢區(qū)域較大時，會增加計算時間，但當(dāng)待檢區(qū)域較小時，邊界回波可能與缺陷散射信號混疊，增大散射信息提取的難度。綜合考慮上述因素，本文擬采用較小區(qū)域模型進行計算，并將模型外邊沿設(shè)置為吸收邊界，以極大程度地減小甚至消除邊界反射回波[11]。吸收邊界由25層0.2 mm寬的阻尼層構(gòu)成，各阻尼層是在待檢材料上增設(shè)瑞利阻尼系數(shù)項而形成的，且從內(nèi)邊界至外邊界各層阻尼系數(shù)值逐層增加。為了獲取到缺陷處的散射信號，以缺陷中心為圓心設(shè)置一個半徑為10 mm的檢測圓，并在其上等間隔設(shè)置32個監(jiān)測點，通過提取監(jiān)測點的信號獲得散射信息。

圖1 相控陣超聲與孔缺陷交互作用的二維有限元模型

2 相控超聲與圓孔缺陷的交互作用

相控陣超聲通過相位延遲可合成不同類型的聲束，按合成超聲波前可分為平面波束和聚焦波束。為此，在相控傳感器的正下方設(shè)置圓孔缺陷，直徑分別為0.5、1、1.5和2，對比研究相控陣超聲平面波束和聚焦波束檢測圓孔缺陷時的散射縱波分布規(guī)律，其中為超聲縱波波長。

2.1 平面波束在圓孔缺陷處的散射特性

在有限元模型中，以相同延遲時間同時激勵32個陣元，可合成平行于激勵表面的平面縱波，簡稱平面波。研究平面波束在不同直徑圓孔處的散射特性時，保持平面波從90°方向入射且激勵能量不變。圖2是相控超聲平面波束與圓孔缺陷交互作用的仿真結(jié)果圖。

(a) 圓孔直徑為0.5(b) 圓孔直徑為

(c)圓孔直徑為1.5(d)圓孔直徑為2

圖2 相控陣超聲平面波與圓孔缺陷交互作用的仿真云圖

Fig.2 Simulation images of phased array focused plane wave interacted with the hole defects of different diameters

如圖2(a)所示，相控陣超聲平面波在經(jīng)過圓孔缺陷時，平面波與圓孔交互作用產(chǎn)生了兩種不同模態(tài)的散射波，即散射縱波(L模態(tài))和散射橫波(S模態(tài))。隨著傳播時間的推移，散射縱波和散射橫波將會在時域上分開。從圖2中可知，相控陣超聲平面波束入射到圓孔缺陷時，在缺陷的反射區(qū)(正反射方向左右90°范圍區(qū)域)存在散射縱波，并且隨著圓孔直徑增加，散射L模態(tài)的能量增加，在透射區(qū)(正透射方向左右90°范圍區(qū)域)不能明顯看出是否存在散射縱波。

2.1.1 周向散射系數(shù)

從信號監(jiān)測圓上提取不同直徑圓孔缺陷的周向散射信息，為了便于對比分析，對不同直徑圓孔缺陷的周向散射系數(shù)進行歸一化處理，并繪制圓孔缺陷散射L模態(tài)的周向散射系數(shù)圖，如圖3所示。

圖3 平面波入射到圓孔缺陷的周向散射系數(shù)

從圖3可知：(1)在反射區(qū)和透射區(qū)均存在散射縱波；(2) 散射L模態(tài)的能量隨著圓孔直徑增加而增加，且透射區(qū)的能量增幅大于反射區(qū)的能量增幅；(3)圓孔散射L模態(tài)的散射能量主要集中于90°及270°方向左右，且透射區(qū)的能量大于反射區(qū)的能量，而在與入射方向垂直的0°及180°方向上能量極其微弱，幾乎沒有散射能量；(4) 散射L模態(tài)關(guān)于入射方向呈軸對稱分布。

2.1.2 全散射系數(shù)矩陣

全散射系數(shù)矩陣用于反映所有入射角下的散射場分布，對于圓孔缺陷可通過對稱變換獲取其他入射角下的周向散射系數(shù)，以此構(gòu)成全散射系數(shù)矩陣。將各圓孔對應(yīng)的全散射系數(shù)矩陣進行歸一化處理，并分別繪制為散射矩陣圖，如圖4所示。

散射矩陣圖中，橫坐標為平面波的入射角度，縱坐標為散射角度，色彩值表示散射系數(shù)值。為了便于描述，引入以下特征詞：“脊線”和“脊帶”代表幅值明顯高于相鄰區(qū)域的線狀區(qū)和帶狀區(qū)；“阱線”和“阱帶”代表幅值明顯低于周圍區(qū)域的線狀區(qū)和帶狀區(qū)。從圖4可知：(1) 隨著圓孔直徑的增加，關(guān)于矩陣副對角線對稱的脊帶寬度逐漸變窄且能量變高；(2) 在副對角線上的脊帶隨著圓孔直徑的增加在逐漸變寬，幅值能量逐漸增大，圓孔直徑增加到2時出現(xiàn)脊線。

2.2 聚焦波束在圓孔缺陷處的散射特性

保持各圓孔缺陷不變，利用相控陣聚焦延時法則，合成聚焦于圓孔缺陷處的聲束，簡稱聚焦波束。為比較聚焦聲束與缺陷的交互作用規(guī)律，采用相同的延遲參數(shù)和激勵信號，同時保持入射方向、聚焦點。圖5展示了相控超聲聚焦波束與圓孔缺陷交互作用的仿真結(jié)果圖。

(a) 圓孔直徑為0.5(b) 圓孔直徑為1

(c) 圓孔直徑為1.5(d) 圓孔直徑為2

圖4 相控超聲平面波與圓孔交互作用的散射系數(shù)矩陣圖

Fig.4 The scattering coefficient matrix chart of the plane wave interacted with the hole defects of different diameters

(a) 圓孔直徑為0.5(b) 圓孔直徑為

(c) 圓孔直徑為1.5(d) 圓孔直徑為2

圖5 90°入射的相控超聲聚焦波束與圓孔缺陷交互作用的仿真云圖

Fig.5 Simulation images of phased array focused ultrasonic beam interacted with the hole defects of different diameters for 90° incidence

如圖5所示，相控陣超聲聚焦波在經(jīng)過圓孔缺陷時，聚焦波束與圓孔缺陷交互作用產(chǎn)生了散射縱波和散射橫波，散射縱波主要出現(xiàn)在反射區(qū)，在透射區(qū)的分布不明顯，且隨著圓孔直徑增加，散射能量增加。此外，如圖5(a)所示，由于在仿真相控陣陣元產(chǎn)生超聲波的過程中，會同時激勵出兩種形態(tài)的超聲波，即縱波和橫波，因此在研究聚焦縱波與圓孔交互作用時，可以觀察到入射的橫波。

2.2.1 周向散射系數(shù)

提取聚焦波束入射到不同直徑圓孔缺陷的周向散射系數(shù)，并繪制圓孔缺陷散射L模態(tài)的周向散射系數(shù)圖，如圖6所示。分析圖6可發(fā)現(xiàn)，聚焦波束在圓孔缺陷處的散射L模態(tài)分布規(guī)律與平面波束入射時的散射場分布相似。

圖6 聚焦波入射到圓孔缺陷的周向散射系數(shù)圖

2.2.2 全散射系數(shù)矩陣

利用對稱變換，獲取聚焦波束入射到不同直徑的圓孔缺陷的全散射系數(shù)矩陣，并分別繪制為全散射系數(shù)矩陣圖，如圖7所示。

從圖7可知：(1) 散射矩陣圖上的脊線、脊帶、阱線和阱帶均平行于散射矩陣的副對角線；(2) 隨著圓孔直徑的增加，關(guān)于散射矩陣副對角線對稱的脊帶變窄，能量變大；(3) 位于副對角線上的脊帶隨著圓孔直徑增加逐漸變寬，能量逐漸增大。

(a) 圓孔直徑為0.5(b) 圓孔直徑為1

(c) 圓孔直徑為1.5(d) 圓孔直徑為2

圖7 相控超聲聚焦波束與圓孔交互作用的全散射系數(shù)矩陣圖

Fig.7 The scattering coefficient matrix chart of the interaction of the focusing beam with the hole defects of different diameters

2.3 平面波與聚焦波在圓孔處散射特性的比較

對兩種波束在同一缺陷處的周向散射系數(shù)進行同一尺度歸一化處理，并繪制在同一個周向散射系數(shù)分布的直角坐標圖中，分析平面波束和聚焦波束在同一圓孔缺陷處的散射場能量，如圖8所示。

(a) 圓孔直徑為0.5(b) 圓孔直徑為1

(c) 圓孔直徑為1.5(d) 圓孔直徑為2

圖8 不同入射波在圓孔處的散射系數(shù)統(tǒng)一歸一化分布圖

圖8中，紅色曲線代表聚焦波，藍色曲線代表平面波。對比兩種不同的相控陣超聲波束分別從90°方向入射到不同直徑圓孔的散射系數(shù)，可以看出聚焦波在圓孔缺陷處L-L模態(tài)的散射能量整體比平面波在圓孔處的散射能量大，特別是在90°與270°方向上聚焦波的散射能量明顯大于平面波的散射能量。

在同一圓孔缺陷處，將兩種不同的入射波的周向散射系數(shù)分別進行獨立的歸一化處理，對比分析兩種不同波束在同一圓孔缺陷處的散射場能量分布情況，如圖9所示。從圖中可以看出：圓孔直徑為0.5和1時，兩種入射波在圓孔處的散射場分布曲線重合，即散射場的分布趨勢相同；在圓孔直徑為1和2時，兩入射波在90°方向附近的散射能量略有差異，但散射場分布趨勢相同，而在其他角度處的散射場分布相同。整體而言，聚焦波和平面波在圓孔缺陷處的散射L模態(tài)的分布趨勢相同。

(a) 圓孔直徑為0.5(b) 圓孔直徑為1

(c) 圓孔直徑為1.5(d) 圓孔直徑為2

圖9 不同入射波在圓孔缺陷處的散射系數(shù)單獨歸一化分布

綜上分析，相控陣超聲聚焦波和平面波與圓孔缺陷交互作用散射L模態(tài)的能量有明顯差異，聚焦波在同一圓孔處的散射場能量明顯高于平面波在圓孔處的散射場能量，但兩者在同一圓孔缺陷處的散射L模態(tài)的分布趨勢是相同的。

3 結(jié)論

本文對相控超聲與圓孔缺陷交互作用的散射場分布特性進行了研究，得出的主要結(jié)論有：

(1) 圓孔缺陷處散射L模態(tài)的能量主要集中在正發(fā)射區(qū)和正透射區(qū)，且透射區(qū)的能量大于反射區(qū)能量，在與入射方向垂直的方向上，散射L模態(tài)的能量極其微弱；

(2) 同一相控陣超聲聲束入射下，圓孔缺陷隨著直徑的增加，散射場的分布趨勢基本不變，但散射能量在增強；

(3) 相控陣超聲聚焦波和平面波入射到同一圓孔缺陷時，聚焦波的散射場能量明顯大于平面波的散射場能量，但兩者在同一缺陷處的散射場分布趨勢相同。

從以上研究結(jié)論中可發(fā)現(xiàn)，散射聲場中包含有圓孔缺陷的特征信息，通過合成不同的相控聲束能夠改善檢測分辨率，這為利用散射聲場的分布特性進行圓孔缺陷反演重構(gòu)，實現(xiàn)缺陷高分辨率成像及準確的定性定量分析提供了理論基礎(chǔ)，對相控陣超聲檢測工藝的制定也具有重要的指導(dǎo)意義。

[1] 施克仁, 郭寓岷, 霍健, 等. 相控陣超聲成像檢測[M]. 北京: 高等教育出版社, 2010: 31-47. SHI Keren, GUO Yumin, HUO Jian, et al. Phased Array Utrasonic Imaging and Testing[M]. Beijing: Higher Education Press, 2010: 31-47.

[2] 周正干, 李洋, 周文彬. 相控陣超聲后處理成像技術(shù)研究、應(yīng)用和發(fā)展[J]. 機械工程學(xué)報, 2016, 52(6):1-11. ZHOU Zhenggan, LI Yang, ZHOU Wenbin. Ultrasonic Phased Array Post- processing Imaging Techniques: A Review[J]. Journal of Mechanical Engineering, 2016, 52(6):1-11.

[3] Drinkwater B W, Wilcox P D. Ultrasonic arrays for non-destructive evaluation: a review[J]. NDT&E International, 2006, 39(7): 525-541.

[4] Zetec, Inc. OVERIEW: Ultra-high performance phased array UT instruments[EB/OL]. [2016-05-18]. http://www.zetec.com/ products/ ultrasound/ instrumentation/ dynaray- dynaray- lite/ dynaray – key - features/.

[5] Holmes C, Drinkwater B W, Wilcox P D. Post-processing of the full matrix of ultrasonic transmit-receive array data for non- destructive evaluation[J]. NDT & E International, 2005, 38(8): 701-711.

[6] P?rtzgen N, Gisolf D, Blacquière G. Inverse wave field extrapolation: a different NDI approach to imaging defects[J]. Ultrasonics, Ferroelectrics, and Frequency Control, IEEE Transactions on, 2007, 54(1): 118-127.

[7] Hunter A J, Drinkwater B W, Wilcox P D. The wavenumber algorithm for full-matrix imaging using an ultrasonic array[J]. Ultrasonics, Ferroelectrics, and Frequency Control, IEEE Transactions on, 2008, 55(11): 2450-2462.

[8] 鄭陽, 何存富, 周進節(jié), 等. 超聲蘭姆波斜入射裂紋時的散射特性[J]. 機械工程學(xué)報, 2013, 49(4): 6-12. ZHENG Yang, HE Cunfu, ZHOU Jìnjie, et al. Scattering Characteristic of Lamb Wave at Crack with Oblique Incident[J]. Journal of Mechanical Engineering, 2013, 49(4): 6-12.

[9] 鄭陽, 何存富, 周進節(jié), 等. 超聲 Lamb波在缺陷處的二維散射特性研究[J]. 工程力學(xué), 2013, 30(8): 236-243. ZHENG Yang, HE Cunfu, ZHOU Jìnjie, et al. 2D scattering characteristics of Lamb wave at defects[J]. Engineering Mechanics, 2013, 30(8): 236-243.

[10] 鄭陽, 周進節(jié). 蘭姆波在裂紋處的模態(tài)轉(zhuǎn)換及散射特性研究[J]. 工程力學(xué), 2014, 31(6): 21-29. ZHENG Yang, ZHOU Jìnjie. Mode Conversion And Scattering Properties of Lamb Waves Near Localized Cracking[J]. Engineering Mechanics, 2014, 31(6): 21-29.

[11] 吳斌, 張也弛, 鄭陽, 等. 超聲導(dǎo)波有限元仿真中吸收邊界設(shè)置及參數(shù)[J]. 北京工業(yè)大學(xué)學(xué)報, 2013, 39(12): 1777-1783.WU Bin, ZHANG Yechi, ZHENG Yang, et al. Modeling and Parameters of Absorbing Boundary for Ultrasonic-guided Wave in FE Simulation[J]. Journal of Beijing University of Technology, 2013, 39(12): 1777-1783.

Research on the scattering characteristics of phased array focused ultrasound interacted with hole

ZHOU Jin-jie1, ZHENG Yang2, ZHANG Zong-jian2

(1. North University of China, School of Mechanical and Power Engineering, Taiyuan 030051, Shanxi, China;2. China Special Equipment Inspection and Research Institute, Key Laboratory of nondestructive testing and evaluation of AQSIQ, Beijing 100029, China)

Compared with the traditional detection methods, the sound beam focusing control ability of phased array technique makes the resolution, signal-to-noise ratio, sensitivity and other performances of ultrasonic phased array detection effectively improved. But at present, lack of understanding the interaction of phased array focused ultrasound with defects causes the current detection techniques impossible to achieve accurately qualitative and precisely quantitative analysis of defects. Therefore, this paper establishes a finite element model of one-dimensional linear ultrasonic phased array for the detection of hole defectsand studies the scattering characteristics of the acoustic field caused by the interaction of the phased array focused beam with the hole defect. The results show that the energy of L mode scattering of the hole is mainly focused on the positive reflection area and the positive transmission area, and the energy in the transmission region is greater than that of the energy in the reflection region, and energy is extremely weak in the direction perpendicular to the direction of incidence, almost no scattering energy. The L mode scattering is axial symmetrical distribution along the direction of incidence. For hole defects of different diameters, the larger the diameter, the greater the energy of the scattering obtains in the same direction. For the same diameter of the hole, the distribution of the scattering energy is the same for focusing wave or plane wave incidence, but the scattering energy of focusing wave is significantly higher than that of plane wave.

phased array; focus ultrasound; hole; scattering characteristics; scattering coefficient matrix

TB52+5 TB553

A

1000-3630(2017)-02-0147-05

10.16300/j.cnki.1000-3630.2017.02.009

2016-07-07;

2016-10-07

質(zhì)檢公益性行業(yè)科研專項(201410026)、(201510068), 質(zhì)檢總局科技計劃項目(2014QK252)

周進節(jié)(1981－), 男, 湖北天門人, 博士, 講師, 研究方向為超聲無損檢測與評價、傳感器與檢測儀器開發(fā)、自動化檢測系統(tǒng)等。

鄭陽, E-mail: zhengyangchina@126.com。