李鴻源，徐 鴻，田振華

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

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

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

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

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對(duì)應(yīng)于Lamb波對(duì)稱(chēng)模態(tài)，α＝π／2對(duì)應(yīng)于非對(duì)稱(chēng)模態(tài)。

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

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

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

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

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

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

2 有限元模擬

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

2.1 有限元模型

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

圖3 有限元模型示意圖

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

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

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

式中A為信號(hào)幅值；fc為中心頻率。

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

2.3 模擬結(jié)果分析

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

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

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

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

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

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

3 損傷成像

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

3.1 損傷成像算法

圖6 Lamb波傳播路徑概圖

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

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

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

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

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

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

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

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

3.2 損傷成像

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

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

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

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

4 結(jié)論

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

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

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

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

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

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

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

［5］徐鴻，李鴻源，王冰.管道中具有累積效應(yīng)的二階諧縱向?qū)Рㄉ牲c(diǎn)的數(shù)值驗(yàn)證［J］.固體力學(xué)學(xué)報(bào)，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.