Lamb波損傷散射及損傷成像的模擬

李鴻源，徐 鴻，田振華

（華北電力大學(xué) 能源動力與機械工程學(xué)院，北京 102206）

在航空航天、石油化工、交通運輸、能源電力領(lǐng)域內(nèi)對部件的完整性進行監(jiān)測和評估十分重要。超聲導(dǎo)波檢測具有傳播距離長的特點，特別適合作為各類部件結(jié)構(gòu)健康監(jiān)測（SHM）的方法，如金屬或復(fù)合材料板、管線、鐵軌、鋼絞線等［1－7］。

Lamb波是存在于板狀材料中的導(dǎo)波，具有傳播距離遠且對損傷敏感的特點，適合對板狀構(gòu)件進行大范圍快速檢測及在線監(jiān)測。但是Lamb波同時具有多模態(tài)及頻散特點，傳播過程中，波場較復(fù)雜。當Lamb波經(jīng)過損傷時，在損傷處發(fā)生反射、透射、散射及模態(tài)轉(zhuǎn)換等，使傳播結(jié)構(gòu)發(fā)生變化，導(dǎo)波波場復(fù)雜性進一步加劇，這將造成損傷特征信號提取困難。因此需要研究Lamb波場與損傷相互作用，以便于損傷識別、特征提取。文章將基于有限元方法，對Lamb波與損傷間的交互作用進行模擬研究。

對于損傷定位，目前已有多種方法，如時間反轉(zhuǎn)法［8］、相控陣列方法［9］及概率的成像方法［10］。文章在分析了Lamb波及損傷波的傳播過程的基礎(chǔ)上，通過傳感器陣列實現(xiàn)對損傷的成像檢測。

1 平板中的Lamb波

1.1 Lamb波傳播模型

Lamb波是在具有上下界面的波導(dǎo)結(jié)構(gòu)中傳播的一種應(yīng)力波，它是由橫波和縱波在上下界面反射、疊加形成的。均質(zhì)平板中Lamb波特征方程為［11］：

式（1）中p和q分別為：

式中h，k，ω分別為平板半壁厚、波數(shù)和角頻率；cl和ct為材料中縱波及橫波的波速；α＝0對應(yīng)于Lamb波對稱模態(tài)，α＝π／2對應(yīng)于非對稱模態(tài)。

通過對特征方程進行數(shù)值求解，可以獲得Lamb波頻散曲線，如圖1所示。通過頻散曲線可以看出，Lamb傳播具有多模態(tài)特點，模態(tài)數(shù)隨頻厚積增大而增多；即使在低頻區(qū)域（A1模態(tài)截止頻厚積之前），也存在兩種Lamb波模態(tài)，即A0模態(tài)和S0模態(tài)。另外，Lamb波傳播具有頻散特點，即傳播速度隨頻厚積變化而變化，若波包內(nèi)不同頻率成分傳播速度不同，將導(dǎo)致波包時域上持續(xù)時間隨傳播距離增加而增長，幅值隨傳播距離增加而顯著降低。

在損傷檢測過程中，Lamb的多模態(tài)和頻散會造成檢測信號復(fù)雜性增加，難于分析出損傷特征信號。因此在實際應(yīng)用中，通常需激發(fā)單一Lamb波模態(tài)作為檢測信號，并通過頻率選擇盡量降低信號的頻散。

1.2 Lamb波模態(tài)選擇

為了控制模態(tài)數(shù)，超聲導(dǎo)波激發(fā)頻率通常選擇在A1模態(tài)截止頻率以前。在此頻率范圍內(nèi)僅存在A0和S0模態(tài)。

圖2為300kHz下，1.5mm厚鋁板中A0和S0模態(tài)位移分布。從圖2可以看出，A0模態(tài)波長相對較短，整個波場以離面位移為主，呈反對稱分布；S0模態(tài)波長相對較長，整個位移場以面內(nèi)位移為主，呈對稱結(jié)構(gòu)。

在低頻區(qū)域，S0模態(tài)頻散較小，傳播速度較快，同時其波場以面內(nèi)位移為主，在外部為液體環(huán)境下能量衰減較小。因此選擇低頻區(qū)的Lamb波對稱模態(tài)（S0模態(tài)）作為損傷檢測信號。

2 有限元模擬

此處將基于有限元方法研究Lamb波的傳播過程，以及損傷對Lamb波傳播的影響。

2.1 有限元模型

基于Ansys有限元分析軟件，建立了350mm×350mm×1.5mm鋁板模型，如圖3所示。以鋁板中心為坐標原點，在坐標為（50mm，50mm）的D點建立直徑為6.35mm的圓形通孔模擬損傷，信號激發(fā)點位于坐標為（－25mm，－25mm）的T點，P點為信號提取點。

圖3 有限元模型示意圖

Lamb波傳播過程仿真屬于瞬態(tài)動力學(xué)分析，模型單元尺寸和時間步長均會影響計算精度和時間。通常網(wǎng)格尺寸和積分步長需要滿足如下關(guān)系：

式中Lmax為單個網(wǎng)格上兩節(jié)點間的最大距離；λmin為最小波長；cmin是Lamb波的最小群速度；nmin是每個波長內(nèi)的最少網(wǎng)格數(shù)，通常取8～10［12］。Lmin，cl分別為最小網(wǎng)格尺寸和縱波波速。根據(jù)式（3）和（4），將模型劃分為邊長1.25mm的solid185單元，時間步長選為0.12μs。

為了在板中激發(fā)單一S0模態(tài)，在激發(fā)點T處上下表面及板內(nèi)各節(jié)點均勻施加面內(nèi)位移載荷。激發(fā)信號采用經(jīng)Hanning窗調(diào)制的5周期正弦波，中心頻率為300kHz，其表達式如下：

式中A為信號幅值；fc為中心頻率。

研究主要關(guān)注損傷與S0模態(tài)間的相互作用。為了去除邊界造成的復(fù)雜影響，首先基于無損傷模型對S0模態(tài)在鋁板中的傳播過程進行仿真。隨后在無損傷模型基礎(chǔ)上，將以D為中心的6.35mm范圍內(nèi)的單元去除，建立損傷模型，基于損傷模型再次對S0模態(tài)的傳播行為進行仿真。兩次仿真結(jié)果之差即為單純由損傷引起的波場變化。

2.3 模擬結(jié)果分析

基于上述方法，獲得的S0模態(tài)傳播波場如圖4所示。

圖4（a）是S0模態(tài)在無損傷板中的傳播過程，分別顯示了37.6，45.7，51.8和58.0μs時整個波場的位移云圖。從中可以看出，S0模態(tài)逐漸向四周傳播，并且位移場強度逐漸減弱。當S0模態(tài)到達邊界時，將在邊界處發(fā)生邊界反射，形成邊界反射波S邊界。

圖4（b）顯示了S0模態(tài)在有損傷板中的傳播過程。當S0模態(tài)經(jīng)過損傷時，將發(fā)生反射及透射，形成以損傷為中心、向四周傳播的位移場，記為S損傷。

為了便于分析損傷對S0模態(tài)傳播的影響，圖4（c）顯示了單純由損傷引起的波場S損傷。從中可以看出，損傷波場以損傷為中心，逐漸向四周擴散。同時，可以觀察到S0模態(tài)在損傷處反射時，發(fā)生模態(tài)轉(zhuǎn)換，兩種模態(tài)由于存在速度差，在傳播過程中逐漸分離。為了確定這兩個波包所對應(yīng)的模態(tài)，在距離損傷點D140mm處分別提取有損和無損兩種情況下的時域信號，兩者之差為損傷反射的時域信號，如圖5所示。第一個波包包絡(luò)最大點對應(yīng)的時間為51.6μs，第二個波包對應(yīng)的傳播時間為68.7μs。

圖5 損傷反射信號及包絡(luò)

根據(jù)頻散特征方程可知，300kHz下S0和SH0模態(tài)波長分別為18.06和10.48mm，傳播速度分別為5.39和3.14mm／μs。因此可判斷外層波包對應(yīng)為反射的S0模態(tài)，內(nèi)層波包為S0模態(tài)，在損傷處發(fā)生模態(tài)轉(zhuǎn)換產(chǎn)生了SH0模態(tài)。

3 損傷成像

通過上述有限元模型分析可知，S0模態(tài)在損傷處會發(fā)生透射、反射，形成損傷波場。以下研究基于不同位置接收信號對損傷定位的方法。

3.1 損傷成像算法

圖6 Lamb波傳播路徑概圖

假設(shè)板內(nèi)C點處存在損傷，如圖6所示。Lamb波從激發(fā)端T到c（x，y），再到接收端R的傳播時間為：

式中cg為導(dǎo)波群速度，dTD和dDR分別是激發(fā)點與損傷間和損傷與接收點間距離，分別為：

因此c點反射強度對應(yīng)在T點激發(fā)，在R點接收的損傷信號在tm處能量。時域信號能量通常用其包絡(luò)表示，信號包絡(luò)定義為：

式中 H（t）是信號u（t）的 Hilbert變換；e（t）是信號u（t）的包絡(luò)。

根據(jù)式（6）和（9）可建立平面內(nèi)任意位置（x，y）處的損傷反射場強值Em（x，y）為：

E（x，y）＝e（t）｜t＝tm（x，y） （10）

損傷圖像中各像素點值I（x，y）可使用多對激發(fā)－接收組合獲得，其計算式為：

式中N為激發(fā)－接收組合數(shù)。

3.2 損傷成像

為了獲得損傷圖像，以激發(fā)點T為中心，30mm為半徑，30°為間隔，布置12個接收傳感器，如圖6所示。

根據(jù)上述損傷成像算法，利用傳感器R－12接收信號可獲得單路徑損傷分布圖像，如圖7所示。圖中圓圈代表真實損傷位置。從中可以看出，基于單一傳感器可以判定多個損傷可能位置，無法對損傷定位。

圖8為12個傳感器獲得的損傷圖像疊加，從中可以看出，多路徑損傷圖像可以指示出損傷位置。在此基礎(chǔ)上，對圖8引入75%閾值限定，可更加清晰地顯示損傷位置，如圖9所示。

通過上述損傷成像過程可知，僅基于一對激發(fā)－接收傳感器組合時，損傷圖像呈單圓環(huán)，無法指示損傷位置。當激發(fā)－接收傳感器組合增加時，不同的圓環(huán)相互疊加，在實際損傷位置，圖像強度疊加增強。通過限定最低閾值的方法可以清晰地顯示出損傷位置。

4 結(jié)論

Lamb波對稱模態(tài)S0具有傳播速度快、低頻散的特點，并且以面內(nèi)位移為主，是損傷檢測常用的導(dǎo)波模態(tài)?；谟邢拊椒?，對S0模態(tài)在鋁板中的傳播過程進行了仿真研究，并分析了圓孔損傷對其傳播的影響。結(jié)果表明，在圓孔處S0模態(tài)會發(fā)生反射、透射，形成損傷波場，并同時伴隨有模態(tài)轉(zhuǎn)換現(xiàn)象。損傷波場以損傷為中心，向四周傳播，同時出現(xiàn)SH0模態(tài)。兩導(dǎo)波模態(tài)傳播速度不同，在傳播過程中逐步分離。

為了定位損傷，在激發(fā)傳感器四周建立了環(huán)形接收陣列，并基于損傷成像算法獲得損傷圖像。研究表明，僅基于單一傳感器路徑，損傷圖像呈圓環(huán)形，無法定位損傷?；趥鞲衅麝嚵?，損傷圖像可以準確指示損傷位置。

然而，當板狀結(jié)構(gòu)中存在裂紋且不斷擴展時，Lamb波場的變化將更為復(fù)雜。裂紋對Lamb波的傳播的影響，裂紋擴展與波場變化之間的關(guān)系，有待下一步的研究。

［1］武新軍，徐江，沈功田.非接觸式磁致伸縮導(dǎo)波管道無損檢測系統(tǒng)的研制［J］.無損檢測，2010，32（3）：166－70.

［2］徐鴻，王冰，姜秀娟.超聲導(dǎo)波用于管道缺陷檢測的數(shù)值研究［J］.動力工程，2009，29（11）：1018.

［3］王強，袁慎芳.復(fù)合材料板脫層損傷的時間反轉(zhuǎn)成像監(jiān)測［J］.復(fù)合材料學(xué)報，2009，26（3）：99－104.

［4］劉增華，張易農(nóng)，張慧昕，等.基于磁致伸縮效應(yīng)在鋼絞線中激勵接收縱向?qū)РB(tài)的試驗研究［J］.機械工程學(xué)報，2010，46（6）：71.

［5］徐鴻，李鴻源，王冰.管道中具有累積效應(yīng)的二階諧縱向?qū)Рㄉ牲c的數(shù)值驗證［J］.固體力學(xué)學(xué)報，2011，32（1）：82－7.

［6］Cawley P，Lowe M J S，Alleyne D N，et al.Practical long range guided wave testing：applications to pipes and rail［J］.Materials Evaluation，2003，61（1）：66－74.

［7］Rose J L，Avioli M J，Mudge P，et al.Guided wave inspection potential of defects in rail［J］.NDT ＆E International，2004，37（2）：153－61.

［8］Wang C H，Rose J T，Chang F K.A synthetic timereversal imaging method for structural health monitoring［J］.Smart Mater Struct，2004，13（2）：415.

［9］Yu L，Giurgiutiu V.In situ 2－D piezoelectric wafer active sensors arrays for guided wave damage detection［J］.Ultrasonics，2008，48（2）：117.

［10］Su Z Q，Wang X M，Cheng L，et al.On selection of data fusion schemes for structural damage evaluation［J］.Struct Health Monit，2009，8（3）：223.

［11］Rose J，Nagy P.Ultrasonic waves in solid media［J］.The Journal of the Acoustical Society of America，2000（107）：1807.

［12］Yang C H，Ye L，Su Z Q，et al.Some aspects of numerical simulation for Lamb wave propagation in composite laminates［J］.Compos Struct，2006，75（1）：267.