高 捷，陳友興，王召巴

(中北大學(xué)信息與通信工程學(xué)院，山西 太原 030051)

眾所周知，在國(guó)防軍工、航空航天、交通運(yùn)輸?shù)阮I(lǐng)域，金屬柱狀物體的應(yīng)用一直很廣泛，例如兵器中的彈體、汽車(chē)的零部件等。然而，在這些柱狀物體的生產(chǎn)過(guò)程中，往往由于溫度的瞬時(shí)變化、介質(zhì)與模具的摩擦等原因，造成物體內(nèi)部或表面產(chǎn)生不同程度的缺陷，這些缺陷常常會(huì)降低產(chǎn)品的使用效果，或在使用過(guò)程中受到各種外在因素的影響而產(chǎn)生結(jié)構(gòu)缺陷損傷程度積累問(wèn)題，引起重大的安全隱患。因而，在生產(chǎn)過(guò)程以及使用過(guò)程中，對(duì)產(chǎn)品進(jìn)行無(wú)損檢測(cè)具有十分重大的意義。

超聲檢測(cè)是目前應(yīng)用最廣的一種缺陷檢測(cè)技術(shù)，由于超聲波的穿透能力強(qiáng)、靈敏度高，且對(duì)人體無(wú)害，已被廣泛應(yīng)用于工業(yè)及高技術(shù)產(chǎn)業(yè)中。超聲波的特點(diǎn)是傳播能量大、方向性好，而且在介質(zhì)內(nèi)傳播過(guò)程中遇到缺陷時(shí)會(huì)產(chǎn)生界面反射或者引起聲速和能量衰減的變化，利用這一特性就可以達(dá)到檢測(cè)缺陷的目的［1］。

超聲波傳播理論復(fù)雜，傳播過(guò)程抽象，況且隨著工業(yè)技術(shù)的不斷發(fā)展，定性檢測(cè)已不能滿(mǎn)足實(shí)踐的要求，需要定量地確定缺陷的幾何特征(形狀、尺寸和取向)和材料性能，這使得超聲波檢測(cè)人員對(duì)于儀器顯示的波形數(shù)據(jù)正確解釋需要豐富的經(jīng)驗(yàn)，導(dǎo)致檢測(cè)結(jié)果會(huì)有一定的主觀(guān)性。在平時(shí)的檢測(cè)中，就需要制作大量的對(duì)比試件進(jìn)行模擬檢測(cè)，分析波形，用實(shí)際檢測(cè)波形與之對(duì)比，由此來(lái)提供檢測(cè)結(jié)果評(píng)判的準(zhǔn)確率。但是試件制作費(fèi)時(shí)費(fèi)力，而且試件的存儲(chǔ)也比較麻煩。為了研究超聲波傳播理論，更好地利用超聲波進(jìn)行檢測(cè)，并且避免以上問(wèn)題，通過(guò)引入有限元分析軟件直觀(guān)的反應(yīng)波的傳播過(guò)程和聲場(chǎng)特性，為缺陷評(píng)判提供理論數(shù)據(jù)和模型參考，提高缺陷檢測(cè)的準(zhǔn)確性。

本文同時(shí)采用有限元和實(shí)驗(yàn)兩種手段對(duì)均勻介質(zhì)中的缺陷檢測(cè)進(jìn)行研究。

1 超聲檢測(cè)方法

主要研究的對(duì)象為圓柱體，如圖1 所示檢測(cè)示意圖。當(dāng)檢測(cè)無(wú)缺陷圓柱體時(shí)，如圖1(a)所示，從探頭發(fā)射出的超聲脈沖波可以直接到達(dá)試件的底部，不會(huì)發(fā)生反射、折射現(xiàn)象。當(dāng)檢測(cè)有缺陷圓柱體時(shí)，超聲波當(dāng)沒(méi)有遇到缺陷時(shí)，是以常速度沿著固定的路線(xiàn)不受干擾地向前傳播的，但當(dāng)遇到缺陷時(shí)，波與它的原來(lái)路徑發(fā)生了偏離，產(chǎn)生折射、反射現(xiàn)象，因此，可以從中提取出缺陷特征回波信號(hào)，如圖1(b)所示。

圖1 圓柱體水浸檢測(cè)示意圖

2 研究對(duì)象以及實(shí)驗(yàn)對(duì)象

實(shí)驗(yàn)選用兩個(gè)高為200 mm、底面半徑為50 mm 的鋁合金圓柱體試件作為研究對(duì)象，其中一個(gè)在構(gòu)件底面的中心位置含有半徑為5 mm 的圓柱形空穴缺陷，如圖2(a)所示;另一個(gè)在底面的中心有一條光滑的長(zhǎng)為100 mm，寬5 mm 的凹槽，如圖2(b)所示。

圖2 實(shí)驗(yàn)試件

使用圖1 所示的超聲檢測(cè)方法，利用圖2(c)所示的實(shí)驗(yàn)系統(tǒng)對(duì)圖2(a)(b)的兩個(gè)試件分別進(jìn)行了超聲信號(hào)采集，得到了兩個(gè)試件的時(shí)域波形，從時(shí)域波形中讀出始波與缺陷回波的時(shí)間差，依據(jù)公式2l=t ×v(l 為超聲波傳播的距離，t為探頭接收到反射回波的時(shí)間，v 為超聲波傳播的速度)計(jì)算得到構(gòu)件下表面與缺陷之間的距離，就可以得出缺陷的位置。對(duì)兩個(gè)試件缺陷回波信號(hào)進(jìn)行傅里葉變換，得到其頻域變化情況。其中，圓形缺陷、槽形缺陷的時(shí)域波形和其頻域變化情況如圖3 所示。

圖3 實(shí)驗(yàn)試件超聲回波時(shí)域圖、頻譜圖

3 有限元模擬仿真結(jié)果

被檢工件為鋁，密度為2 700 kg/m3，彈性模量為70 GPa，泊松比為0.33，截面半徑為50 mm，圓孔缺陷半徑為5 mm，缺陷的截面模型如圖4(a)所示。為獲取比較細(xì)膩的波的傳播圖像，需對(duì)模型劃分很細(xì)的網(wǎng)格。在波的傳播計(jì)算中，單元在每個(gè)步長(zhǎng)計(jì)算時(shí)沿波傳播方向的長(zhǎng)度應(yīng)小于波長(zhǎng)。根據(jù)參考文獻(xiàn)［8］，選取單元尺寸小于激勵(lì)頻率下最短波長(zhǎng)的十分之一。如果單元尺寸大于波長(zhǎng)的十分之一，誤差將迅速增大。但如果單元尺寸選得太小，對(duì)于同樣面積的模型，單元數(shù)量將以尺寸縮小程度2 次冪的速度增長(zhǎng)，且時(shí)間步長(zhǎng)更小，因此計(jì)算的成本將迅速增大。同時(shí)，當(dāng)單元尺寸小到一定程度時(shí)，舍入誤差將隨之增加。因而，在選取單元尺寸時(shí)，并非越小越好。所以，本文設(shè)置的單元尺寸是0.1 mm。

圖4 有限元模型及激勵(lì)

利用有限元模擬建立的模型參數(shù)與試驗(yàn)試件一致，也采用縱波垂直入射，其激勵(lì)源為漢寧窗調(diào)制的信號(hào)，激勵(lì)頻率為2.5 MHz，與實(shí)驗(yàn)室探頭發(fā)生的激勵(lì)信號(hào)一致。激勵(lì)波形如圖4(b)所示。

在有限元計(jì)算中，從理論上講，時(shí)間步長(zhǎng)越小，計(jì)算結(jié)果越精確，但計(jì)算時(shí)間也相應(yīng)增長(zhǎng)。在本文中，根據(jù)所施加激勵(lì)的中心頻率，選用的時(shí)間步長(zhǎng)為0.1 us。采用有限元的計(jì)算，2.5 MHz 的縱波在介質(zhì)鋁中傳播過(guò)程的示意圖如圖5所示。

圖5 縱波在介質(zhì)中傳播過(guò)程示意圖

從圖5 中可以看出，波在均勻的介質(zhì)鋁中傳播，當(dāng)沒(méi)有遇到圓形缺陷時(shí)，是以常速度沿著固定的路線(xiàn)不受干擾地向前傳播的，擔(dān)當(dāng)遇到圓形缺孔時(shí)，波的傳播路線(xiàn)發(fā)生了變化，即在其他的、未受干擾的入射波的作用下，圓形的孔缺陷起到了一個(gè)次生波源的作用，形成了波的散射，而且波與它的原來(lái)路徑發(fā)生了偏離，產(chǎn)生了衍射現(xiàn)象。

圖6 有限元仿真時(shí)域圖、頻譜圖

圖6(a)(b)所示為圖4(a)有限元模型中A 點(diǎn)的回波信號(hào)圖，其中(a)圖為有限元仿真結(jié)果圖，(b)圖為缺陷處回波頻譜圖。對(duì)比圖3(a)(b)和圖6(a)(b)，可以發(fā)現(xiàn)，有限元仿真結(jié)果圖與實(shí)驗(yàn)采集的回波信號(hào)基本一致，由此可以驗(yàn)證仿真結(jié)果的準(zhǔn)確性。

運(yùn)用以上同樣的方法，用有限元ANSYS 仿真槽形缺陷試件，得到的時(shí)域回波信號(hào)以及其缺陷處回波信號(hào)的頻譜圖，如圖6(c)(d)所示。

對(duì)比圖6 和圖3，我們可以看出仿真結(jié)果和實(shí)驗(yàn)結(jié)果基本一致。對(duì)比圖3(b)和圖6(b)，圖3(d)和圖6(d)，我們可以看出，不同缺陷回波信號(hào)的頻譜圖是不一樣的，因此，我們可以根據(jù)缺陷處的回波信號(hào)及其頻譜圖來(lái)判斷試件中缺陷的位置、形狀及大小。

4 結(jié)論

本文通過(guò)有限元法和實(shí)驗(yàn)法兩種手段對(duì)鋁圓柱體內(nèi)缺陷進(jìn)行檢測(cè)，經(jīng)過(guò)仿真結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比，可以驗(yàn)證仿真結(jié)果的準(zhǔn)確性。對(duì)缺陷處回波信號(hào)進(jìn)行頻域處理，對(duì)比不同缺陷回波信號(hào)的頻譜圖，可以看出缺陷的形狀不同，其頻譜圖是可以明顯看出不同的，因此，我們可以根據(jù)缺陷回波信號(hào)的頻譜圖，基本確定缺陷的形狀。同時(shí)，我們也可以看出，ANSYS 能有效地仿真超聲波的傳播，構(gòu)建超聲場(chǎng)，表現(xiàn)超聲場(chǎng)特性。因此在實(shí)際檢測(cè)中可以用ANSYS 仿真超聲波傳播，對(duì)于開(kāi)發(fā)研究超聲波檢測(cè)新技術(shù)來(lái)說(shuō)，這點(diǎn)尤其重要，一方面通過(guò)仿真計(jì)算不斷地修正設(shè)計(jì)參數(shù)，達(dá)到預(yù)期的設(shè)計(jì)效果;另一方面可以避免大量制作模擬試塊進(jìn)行檢測(cè)驗(yàn)證;同時(shí)可以節(jié)約資金及時(shí)間，降低開(kāi)發(fā)成本，提高開(kāi)發(fā)效率，具有實(shí)用意義。

［1］張俊哲.無(wú)損檢測(cè)技術(shù)及其應(yīng)用［M］.北京:科學(xué)出版社，2010:1－2.

［2］Kosmatka B J，James M R.Damage Detection In Structures by Model Vibration Characterizati［J］.Journal of Structural Engineering，1999，12(12):1384－1392.

［3］顧雯琪，王召巴，陳友興.基于Born 近似法的鋁合金構(gòu)件缺陷的三維反演［J］.激光與光電子學(xué)進(jìn)展，2014，11:185－190.

［4］劉長(zhǎng)福，牛曉光，李中偉，等.基于ANSYS 的超聲縱_橫波傳播仿真計(jì)算［J］.無(wú)損檢測(cè)，2008，30(6):15－18.

［5］張偉志，鋼鐵，王軍.超聲波檢測(cè)計(jì)算機(jī)模擬和仿真的研究及應(yīng)用現(xiàn)狀［J］.應(yīng)用聲學(xué)，2003，3(1):41－46.

［6］郭文靜，陳友興，金永，等.基于多元高斯聲束模型的圓柱體超聲檢測(cè)聲場(chǎng)仿真［J］.應(yīng)用聲學(xué)，2013，32(5):354－360.

［7］Zheng Gangfeng，Wu Bin，He Cunfu Born.Inversion Technique for Ultrasonic Scattering Measurements［J］.Journal of Applied Sciences，2013，13(21):4399－4404.

［8］王勖成.有限單元法［M］.北京:清華大學(xué)出版社，2003:2－6.

［9］Hayashi T，Song W J，Rose J L.Guided Wave Dispersion Curves for a Bar Witan Arbitrary Cross-section a Rod and Rail Example［J］.Ultrasonics，2003，41:175－183.

［10］Yim H，Sohn Y.Numerical Simulation and Visualization of Elastic Waves Using Mass-Spring Lattice Model［J］.IEEE Trans UFFC，2000，47(3):549－558.