基于有限元數(shù)值方法的金屬材料內(nèi)壁缺陷超聲頻譜特性研究

孫繼華，趙揚，南鋼洋，馬健，劉帥，巨陽

(山東省科學(xué)院激光研究所，山東省無損檢測工程技術(shù)研究中心，山東 濟南 250103)

?

【新材料】

基于有限元數(shù)值方法的金屬材料內(nèi)壁缺陷超聲頻譜特性研究

孫繼華，趙揚，南鋼洋，馬健，劉帥，巨陽

(山東省科學(xué)院激光研究所，山東省無損檢測工程技術(shù)研究中心，山東 濟南 250103)

根據(jù)金屬材料內(nèi)壁缺陷的特點，建立了吸收邊界和自由邊界相結(jié)合的模擬模型，對含有不同寬度內(nèi)壁凹陷形平底槽缺陷的模型進(jìn)行了超聲頻譜特性模擬研究。通過對波場快照、時域信號和頻譜分析得到了不同尺寸平底槽對聲波傳播頻譜特性的影響以及中心頻率、帶寬、半波高寬度和峰值等呈規(guī)律性變化的頻域特征量。研究結(jié)果對超聲檢測頻譜分析技術(shù)的理論和實驗提供了參考。

金屬材料；超聲檢測；模擬；內(nèi)壁缺陷；頻域特征量

隨著現(xiàn)代超聲檢測技術(shù)的發(fā)展，模擬研究越來越多地應(yīng)用到相關(guān)領(lǐng)域中。板材、管道和壓力容器等金屬材料的內(nèi)壁由于腐蝕、外力以及摩擦等作用，易形成凹陷形的內(nèi)壁缺陷，這種缺陷對金屬材料的性能影響極大，不及時發(fā)現(xiàn)就會帶來安全隱患。超聲無損檢測具有無損、快捷、成本低及對人體無害等特點[1]，是對該類型缺陷較為可行、有效的檢測方法，但傳統(tǒng)超聲檢測方法通常在時域中對檢測結(jié)果進(jìn)行定性分析，而超聲檢測信號中通常攜帶了大量的結(jié)構(gòu)信息，不同結(jié)構(gòu)的時頻特性不同[2]，所以目前對檢測信號的分析正在從時域法向頻域法發(fā)展，頻域分析法為超聲無損檢測技術(shù)注入了活力，一定程度上起到了變革作用[3-4]。通過對超聲檢測信號進(jìn)行頻譜分析可以獲得缺陷在頻域中的各種特征，進(jìn)而對缺陷進(jìn)行進(jìn)一步的識別。

數(shù)值模擬方法可以應(yīng)用于超聲無損檢測技術(shù)的研究，是理論和實驗重要的驗證手段，其中有限元模擬方法理論基礎(chǔ)比較成熟，具有通用性和靈活性，在模擬聲場方面具有獨特的優(yōu)點，可以正確地模擬金屬材料中超聲傳播情況，是一種有效的數(shù)值模擬工具[5-6]。國內(nèi)外學(xué)者近十幾年來應(yīng)用有限元數(shù)值模擬方法對超聲檢測進(jìn)行了一系列的研究[7-11], 但針對金屬材料內(nèi)壁缺陷超聲頻譜特性方面的研究還比較少。為了深入地認(rèn)識超聲波在有凹陷形的內(nèi)壁缺陷金屬中的傳播情況，本文采用有限元數(shù)值模擬方法研究超聲波在金屬材料中的傳播，通過對所得信號進(jìn)行頻譜分析，得到了不同尺寸內(nèi)壁平底槽對超聲波頻譜特性的影響，模擬結(jié)果可以為超聲檢測頻譜分析技術(shù)的理論和實驗研究提供參考。

1 建立有限元數(shù)值模型

1.1 金屬材料中超聲檢測有限元數(shù)值模擬原理[12]

在各向同性彈性固體中，聲波的Navier控制方程的矢量形式為：

(1)

式中，λ和μ為Lamé常數(shù)，ρ為介質(zhì)密度，u為位移。

按照有限元數(shù)值方法計算步驟，首先選取合適網(wǎng)格對計算區(qū)域離散化，并在單元內(nèi)選擇合適的插值函數(shù)得到單元結(jié)點位移、速度和加速度；利用伽遼金(Galerkin)法，由式(1)可導(dǎo)出整體有限元運動方程：

(2)

1.2 建立模型

為了模擬金屬材料中的超聲波傳播，選用在民用、軍工等行業(yè)比較常用的鋁板為母材模型，鋁板的縱波速度vP為6 200 m/s。本文從有限元數(shù)值模擬方法基本理論出發(fā)，采用邊長為0.1 mm的三角形單元對一個長22 mm，厚度為16 mm的鋁板材料進(jìn)行有限元剖分[13]，建立的有限元模擬模型如圖1所示。圖1中鋁板的上表面和下表面采用自由邊界用于模擬鋁-空氣界面，兩側(cè)采用吸收邊界用于模擬鋁板無限延伸[11]。在模型底部設(shè)置了3種深度均為0.8 mm，寬度則分別為4、5、6 mm的平底槽缺陷，用于模擬材料的內(nèi)壁缺陷。

圖1 有限元模型Fig.1 Finite element model

為了使模擬結(jié)果與實際超聲檢測時Φ10 mm超聲探頭產(chǎn)生的平面縱波聲場一致，本文在模擬模型上表面中心10 mm區(qū)域內(nèi)加載高斯分布的應(yīng)力邊界條件以激發(fā)出脈沖超聲縱波[14]，以模擬超聲波激勵源，該激勵源同時可以作為接收源，用于模擬接收超聲波信號并以A-掃描形式顯示出來。

2 模擬結(jié)果與分析

2.1 波場快照圖分析

模擬過程中超聲波由上表面激發(fā)出向下傳播，對超聲波分別傳播到0.51、2.58、3.39、4.38 μs時刻的全波場y方向位移波場快照進(jìn)行分析，不同時刻波場快照如圖2所示，從圖上能明顯地看出模擬激發(fā)的超聲波在鋁板中的傳播情況。

圖2 不同時刻波場快照圖Fig.2 Wave field snapshots at different moments

本文模擬的是垂直入射的超聲縱波，由圖2可以清晰地看到在0.51 μs時，模擬激發(fā)的超聲縱波P已經(jīng)開始向下傳播，在2.58 μs時超聲波已經(jīng)傳播到底面的平底槽位置，在3.39 μs和4.38 μs時可以看到超聲波碰到平底槽表面和底面后發(fā)生的反射和散射現(xiàn)象。從圖上可以清晰地觀察到很多類型的超聲波，其中信號比較明顯的波形有反射縱波R-P和模式轉(zhuǎn)換的橫波M-S等，此后時間里這些不同類型的超聲波會按自己的速度和方向繼續(xù)傳播。從圖2上也可以看出，本文建立的有限元模型中加了吸收邊界條件的兩個側(cè)面沒有明顯的超聲反射波，說明模型加的吸收邊界條件效果很好，達(dá)到了模擬鋁板的無限延伸的目的。所以通過以上波場快照圖分析，直觀地顯示了超聲波在帶有底面凹槽鋁中傳播的情況，模擬的結(jié)果與實際相符。

2.2 時域信號分析

為了獲得時域A-掃描信號，我們在上表面激發(fā)出超聲波的區(qū)域內(nèi)對y方向位移進(jìn)行掃描，可以得到圖3所示的不同缺陷尺寸下A-掃描曲線圖，其中圖3a是底面沒有缺陷時的A-掃描曲線圖，圖3b、c、d是底面有不同尺寸的平底槽時的A-掃描曲線圖。

從鋁板上表面激發(fā)出的超聲波向下傳播過程中遇到鋁板底面和平底槽時就會發(fā)生反射和散射，從圖3a可以看出上表面激發(fā)的始波(P)和底面回波(R-P)信號，這和實際的超聲檢測儀顯示的A-掃描形式一致。為了確定獲得信號的正確性，我們把圖3a中無缺陷底面回波R-P的接收時間t=5.115 μs，代入公式d=tvP/2，可以計算出鋁板厚度為15.9 mm，與實際模型厚度16 mm誤差僅為0.6%，采用同樣方法計算出的各種平底槽深度也與實際模型數(shù)值相符，說明模擬的超聲波傳播過程是和實際相符的。

2.3 頻譜分析

為了獲得超聲檢測信號頻域特征，對底面回波信號進(jìn)行頻譜分析。如圖3b中所示，我們選取合適長度的時間閘門分別截取圖3a、b、c、d中A-掃描曲線上的底面回波，經(jīng)頻譜分析得到幅度譜圖如圖4所示，同時計算得到頻域特征量如表1所示。

由圖4可以看出，當(dāng)鋁板底面有不同尺寸的平底槽時接收到的回波幅度比沒有平底槽時的幅度大，而且不同尺寸平底槽的回波幅度也不同。由表1可以看出，當(dāng)鋁板底面有寬度為4 mm和5 mm的平底槽時，頻域特征量中的中心頻率和峰值頻率比沒有平底槽時要大；而當(dāng)有寬度為6 mm的平底槽時，其中心頻率和峰值頻率比沒有平底槽時要小。由表1還可以看到當(dāng)有不同寬度的平底槽時，頻域特征量中的半高寬度、帶寬等特征量比沒有平底槽時都要小，這也反映了平底槽對不同頻域特征量的影響。

圖3 A-掃描曲線Fig.3 The A-scanning curve

圖4 回波幅度譜Fig.4 Amplitude spectrum of echo

表1 頻域特征量Table 1 Frequency domain characteristics

為了得到鋁板中不同尺寸的平底槽對超聲頻域特征量的影響，我們把平底槽寬度d以1 mm為間隔從1 mm增加到12 mm分別進(jìn)行模擬，通過對結(jié)果進(jìn)行計算分析，可以得到峰值、中心頻率、半高寬度和帶寬等的變化趨勢如圖5所示。

圖5 頻域特征量變化圖Fig.5 Variation chart of frequency domain characteristics

從圖5可以看出，峰值變化的曲線是隨著平底槽寬度的增加首先呈上升趨勢，當(dāng)平底槽寬度在6 mm左右時出現(xiàn)極大值，隨后就開始變小。頻域特征量中心頻率在平底槽寬度為2 mm左右時有一個極大值，然后開始變??；在8 mm左右出現(xiàn)最小值，此后又開始變大。半高寬度變化曲線在平底槽寬度為4 mm左右有極小值，另外在11 mm左右有極大值，隨后開始變小。帶寬和半高寬度的變化趨勢基本一致，只是實際數(shù)值比半高寬度的要小。

通過對頻域特征量變化規(guī)律的分析可知，隨著鋁板中平底槽尺寸的變化，會引起各種頻域特征量規(guī)律性的變化，這可以為缺陷識別提供一定的參考。

3 結(jié)論

本文通過有限元數(shù)值方法模擬計算了超聲波在有平底槽缺陷鋁板中的傳播情況，通過波場快照、時域A-掃描圖和頻譜等分析了鋁板中不同尺寸內(nèi)壁平底槽缺陷對超聲傳播的影響，主要得出以下結(jié)論：

(1)通過模擬過程中的波場快照圖，能清晰地觀察到超聲波在鋁板中的傳播以及遇到底面和平底槽時的反射、散射和模式轉(zhuǎn)換等各種現(xiàn)象，模擬結(jié)果與實際情況相符。

(2)由A-掃描圖上無缺陷底面回波時間計算出的模擬區(qū)域厚度與模型厚度一致，說明采用有限元數(shù)值模擬方法可以正確地模擬超聲在鋁板中的傳播情況。

(3)通過對鋁板不同尺寸平底槽回波進(jìn)行頻譜分析，得到各種頻域特征量峰值、中心頻率、半高寬和帶寬等呈規(guī)律性的變化，可以為缺陷識別提供一定的參考。

通過以上分析可知，采用有限元數(shù)值方法模擬超聲波在金屬材料中的傳播是可行的，通過頻譜分析獲得的結(jié)果可以為實際中的內(nèi)壁凹陷形缺陷識別提供依據(jù)，并為超聲檢測頻譜分析技術(shù)的理論和實驗研究提供參考。但是本文只對金屬內(nèi)壁較規(guī)則的平底槽進(jìn)行了模擬計算，下一步將開展針對金屬內(nèi)壁不同形狀缺陷的模擬計算研究，以期獲得不同類型缺陷對超聲頻譜特性的影響。

[1]齊子誠，鄭子軍，唐盛明，等.滾動式鋼板超聲檢測設(shè)備研制[J].兵器材料科學(xué)與工程,2016,39(2):76-79.

[2]艾春安,蔡笑風(fēng),劉繼方,等.短時傅里葉變換聲-超聲檢測信號缺陷識別[J].中國測試,2015,41(4):29-31.

[3] 劉立群,趙振南.頻譜分析技術(shù)在超聲檢測中的應(yīng)用[J].一重技術(shù),2013(5):57-61.

[4]蔣志峰,吳作倫,吳瑞明,等.材料缺陷超聲檢測中譜分析技術(shù)的研究[J].兵器材料科學(xué)與工程,2009,32(4):117-120.

[5]FELICE M V， VELICHKO A， WILCOX P D,e tal. Obtaining geometries of real cracks and using an efficient finite element method to simulate their ultrasonic array response[J].Insight: Non-Destructive Testing and Condition Monitoring, 2014, 56(9): 492-498 .

[6] 趙軍輝，魏勤，袁媛,等. Lamb波檢測板中裂紋的有限元模擬[J].壓電與聲光,2013,35(3):320-324.

[7]XUE T,LORD W,UDPA S. Numerical analysis of the radiated fields of ultrasonic transducers [J]. Journal of Nondestructive Evaluation, 1995,14(3): 137-146.

[8]DATTA D,KISHORE N N. Features of ultrasonic wave propagation to identify defects in composite materials modelled by finite element method [J]. NDT & E International,1996,29(4): 213-223.

[9]KISHORE N N,SRIDHAR I,IYENGR N G R.Finite element modeling of the scattering ultrasonic waves by isolated flaws[J]. NDT & E International,2000, 33(5):297-305.

[10]BLANLOEUIL P, MEZIANE A, BACON C. 2D finite element modeling of the non-collinear mixing method for detection and characterization of closed cracks [J]. NDT & E International,2015, 76: 43-51.

[11]孫朝明.超聲橫波檢測的有限元模擬計算與分析[J].壓電與聲光,2014,36(6):1018-1021.

[12] 李樹榜，李書光，劉學(xué)鋒.裂紋超聲散射的有限元模擬[J].無損檢測,2007,29(1):27-31.

[13]LIU T,LIU K,ZHANG J. Triangular grid method for stress-wave propagation in 2-D orthotropic materials[J]. Archive of Applied Mechanics, 2005, 74(7):477-488.

[14] 韓秀甫,滕永平,吳迪,等.超聲換能器暫態(tài)聲場的理論分析與計算[J].安慶師范學(xué)院學(xué)報(自然科學(xué)版),2004,10(3):25-28.

Finite element based research on ultrasonic spectral characteristic of inner wall defect of metal materials

SUN Ji-hua, ZHAO Yang, NAN Gang-yang, MA Jian, LIU Shuai，JU Yang

(Shandong Provincial Research Center of Non-destructive Testing Engineering Technology,Laser Institute, Shandong Academy of Sciences, Jinan 250103, China)

∶We construct an absorption boundary and free boundary combined model for inner wall defect of metal materials. We further simulate ultrasonic spectral characteristic of concave groove of different width inner wall. We also acquire the impact of different sized groove on spectral characteristic of ultrasonic propagation by wave field snapshots, time-domain signal and spectral analysis, and such regularly variant frequency domain characteristics as central frequency, bandwidth, half width and peaks. The results provide a reference for the theory and experiment of ultrasonic testing spectral analysis.

∶ metal materials; ultrasonic testing; simulation; inner wall defect; frequency domain characteristic quantity

10.3976/j.issn.1002-4026.2016.05.012

2016-05-05

山東省科學(xué)院杰出青年計劃( 201505)；濟南市高校院所自主創(chuàng)新基金( 201401220，201401209)

孫繼華(1979—)，男，助理研究員，研究方向為材料的無損檢測與評價研究。

TP391.9；TB52+6

A