曾勝洋 賈璐 張書增? 李雄兵 王猛

1)(中南大學(xué)交通運(yùn)輸工程學(xué)院,長沙 410075)

2)(清華大學(xué)航天航空學(xué)院,應(yīng)用力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)

3)(清華大學(xué),柔性電子技術(shù)研究中心,北京 100084)

為解決非線性聲表面波的求解難題,本文從二階非線性各向同性介質(zhì)的超彈性本構(gòu)方程出發(fā),采用位移勢函數(shù)法,建立二維表面波的非線性勢函數(shù)控制方程;通過微擾法推導(dǎo)非線性表面波的準(zhǔn)線性解和絕對非線性系數(shù),討論表面波二次諧波解的主要組成部分;并建立模擬非線性表面波傳播的有限元模型,位移幅值的仿真結(jié)果與理論符合良好,驗(yàn)證了本文非線性表面波理論的準(zhǔn)確性.根據(jù)微擾解的數(shù)值結(jié)果,探討了非線性表面波的傳播以及非線性系數(shù)的特性,結(jié)果表明: 表面波二次諧波由累積項(xiàng)及非累積項(xiàng)組成,前者與表面波縱波分量自相互作用相關(guān),但當(dāng)初始條件和傳播距離相同時(shí),該部分諧波幅值比純縱波的二次諧波幅值大;此外,縱波和表面波的非線性系數(shù)存在正比關(guān)系,該比例關(guān)系由材料的二階彈性系數(shù)確定.本文探究的非線性表面波的傳播特性及其絕對非線性系數(shù)的定義表達(dá)式,對指導(dǎo)非線性表面波的實(shí)際應(yīng)用具有一定意義.

1 引言

近年來,非線性超聲理論和技術(shù)得到了廣泛的研究,并在工業(yè)無損檢測領(lǐng)域開展了相關(guān)應(yīng)用[1].材料的非線性系數(shù)對高階彈性常數(shù)變化和聲-微結(jié)構(gòu)作用十分敏感,因此通過測量非線性超聲的高次諧波,可有效檢測和評價(jià)疲勞損傷、熱損傷、粘接分層和微裂紋等缺陷信息[2,3].當(dāng)前非線性超聲檢測主要采用縱波[4]、表面波[5]、蘭姆波[6]以及導(dǎo)波[7]等波形,其中非線性表面波具有傳播距離遠(yuǎn)、能量集中和單邊測量的特性,更適用于實(shí)際檢測應(yīng)用,因此相關(guān)理論和實(shí)驗(yàn)研究備受關(guān)注.

在非線性表面波的應(yīng)用方面,國內(nèi)外學(xué)者開展了一系列工作.例如,通過非線性表面波檢測玻璃的表面刮擦損傷,體現(xiàn)了其快速評價(jià)表面損傷的潛力[8];相關(guān)實(shí)驗(yàn)表明非線性表面波系數(shù)與鈦合金的塑性形變程度呈正相關(guān)[9];此外,利用空耦探頭測量混凝土中的非線性表面波,可有效區(qū)分因堿硅反應(yīng)導(dǎo)致的不同程度的損傷[10].然而,當(dāng)前非線性表面波的實(shí)驗(yàn)研究普遍測量損傷試塊的相對非線性系數(shù),并需要相同實(shí)驗(yàn)條件下的初始試樣測量結(jié)果作基準(zhǔn),通過觀測該系數(shù)的變化做出相對性評價(jià)[11].相關(guān)研究表明,通過測量材料的絕對非線性系數(shù),不僅可以對材料彈性常數(shù)等特性進(jìn)行評價(jià),還能夠?qū)Σ牧系膿p傷狀態(tài)等進(jìn)行定量或定性分析[3].因此利用絕對表面波非線性系數(shù)進(jìn)行材料特性或損傷評價(jià)具有更大的優(yōu)勢,而研究非線性表面波的相關(guān)理論并提出對應(yīng)測量方法是保證其應(yīng)用的基礎(chǔ).

縱波非線性特性已得到深入的研究和廣泛認(rèn)知[12,13],但非線性表面波非線性的研究仍存在不同的認(rèn)識,受當(dāng)前相關(guān)理論差異的影響,難以實(shí)現(xiàn)表面波絕對非線性系數(shù)的計(jì)算和進(jìn)一步的實(shí)驗(yàn)測量.針對非線性聲波問題主要有兩類研究方法,一是數(shù)值法,以有限差分法或有限元法為代表[14],目前多用于復(fù)雜結(jié)構(gòu)的聲信號模擬與定性分析[15];另一類是解析法,通常可以給出理想問題的準(zhǔn)確解或近似解.前者先將聲波方程離散化,輸入相關(guān)材料參數(shù),再求解出瞬態(tài)聲場;雖可獲得特定復(fù)雜問題的數(shù)值解,但難以建立各物理參數(shù)之間的具體聯(lián)系,即探究相關(guān)問題的物理本質(zhì).相比之下,求非線性表面波問題的解析解,對于非線性表面波特性、計(jì)算乃至測量表面波絕對非線性系數(shù)具有更重要的意義.

針對利用解析法進(jìn)行非線性表面波求解的難題,有學(xué)者通過能流密度[16]等方法進(jìn)行推導(dǎo),但重點(diǎn)研究的是表面波的諧波產(chǎn)生和波形演化等特性,研究中包含了多組表面波非線性系數(shù),且并未給出具體的檢測方法.Herrmann等[17]基于聲表面波的縱波分量推導(dǎo)出表面波非線性系數(shù),并給出該系數(shù)的測量表達(dá)式,在工程實(shí)踐中得到一定應(yīng)用.但Masurkar和Tse[18]指出僅利用表面波縱波分量推導(dǎo)并不能得到真實(shí)的非線性系數(shù),他們認(rèn)為橫波分量對表面波非線性存在貢獻(xiàn)且不容忽視,并提出一種基于測量諧波幅值的表面波非線性系數(shù),然而該方法與數(shù)值法類似,無法解釋諧波隨距離增加的累積效應(yīng),且不能描述表面波的諧波傳播特性.

針對當(dāng)前研究中的不足,本文在準(zhǔn)線性理論成立條件下開展非線性表面波傳播特性研究,并解決表面波非線性系數(shù)定義等問題.為了能夠應(yīng)用解析法對非線性表面波方程進(jìn)行求解,引入位移勢函數(shù)描述非線性表面波的波動控制方程,結(jié)合表面波的邊界條件和微擾法求得表面波基波和二次諧波解,進(jìn)而研究表面波非線性系數(shù)的定義和測量表達(dá)式.為驗(yàn)證理論求解的正確性,引入有限元模型仿真對非線性表面波的傳播進(jìn)行模擬.通過數(shù)值計(jì)算分析表面波、縱波二次諧波的傳播特性并探究幾種非線性系數(shù)的聯(lián)系,為后續(xù)非線性表面波技術(shù)的發(fā)展和應(yīng)用提供理論依據(jù).

2 理論

2.1 二階非線性表面波傳播的控制方程

本文主要研究在半無限大、各向同性的固體表面?zhèn)鞑サ淖畹湫偷穆暠砻娌?即瑞利波.如圖1 的笛卡爾坐標(biāo)系所示,非線性表面波沿x1方向傳播,并在x3方向迅速衰減,忽略聲波傳播時(shí)衍射的影響,即各物理量不隨坐標(biāo)x2變化.當(dāng)該表面波傳播時(shí)滿足平面應(yīng)變條件,質(zhì)點(diǎn)的運(yùn)動方程表示如下:

圖1 表面波傳播坐標(biāo)系Fig.1.The coordinate for the propagation of the surface wave.

式中Pij,j表示?Pij/?xj,ui,tt表示?2ui/?t2(該 簡寫標(biāo)記方式將貫穿本文);Pij為第一類皮奧拉-基爾霍夫應(yīng)力(Piola-Kirchhoff stress,PK 應(yīng)力),ui為xi方向的位移,t為時(shí)間,ρ0為材料密度;式中i,j滿足愛因斯坦求和約定.針對常見金屬材料的超彈性本構(gòu)特性,采用Landau和Liftshitz[19]提出的應(yīng)變能表達(dá)式來描述:

式中W是彈性體每單位質(zhì)量比應(yīng)變能;λ和μ是材料的拉梅常數(shù);A,B和C是材料的三階彈性常數(shù);e是格林應(yīng)變張量.方程(1)中的PK 應(yīng)力與應(yīng)變能之間的關(guān)系為

式中F是變形梯度,Fijδij+ui,j.將方程(2)和方程(3)代入到方程(1),僅考慮到二階非線性項(xiàng),可得非線性表面波質(zhì)點(diǎn)位移的運(yùn)動方程:

式中α=1,β=3 或α=3,β=1.方程(4)等號右側(cè)表示為非線性項(xiàng),其中包含A,B和C的項(xiàng)體現(xiàn)了材料本構(gòu)的非線性;包含λ和μ的項(xiàng)體現(xiàn)了應(yīng)力應(yīng)變關(guān)系的幾何非線性.當(dāng)此項(xiàng)為0 時(shí),方程(4)表示為線性表面波的質(zhì)點(diǎn)位移運(yùn)動方程.表面波在自由邊界傳播時(shí),即在x30 處,其應(yīng)力滿足:

式中P13和P33是應(yīng)力分量,二者均包括線性項(xiàng)和非線性項(xiàng),在方程中分別用上標(biāo)L和NL 表示,具體表達(dá)式為

方程(4)—方程(10)構(gòu)成了二維非線性表面波質(zhì)點(diǎn)運(yùn)動位移場問題的控制方程.

2.2 微擾法求解

引入位移勢函數(shù)和微擾法推導(dǎo)非線性表面波位移場的準(zhǔn)線性解,即二次諧波幅值遠(yuǎn)小于基波幅值,且諧波幅值處于累積范圍中.引入縱波勢φ和橫波勢?,他與位移有如下關(guān)系:

根據(jù)二階微擾近似法,將方程的解表示為線性解和二階非線性解部分,即??(1)+?(2),φφ(1)+φ(2),且?(2)??(1),φ(2)?φ(1),結(jié)合方程(4)—方程(11),可將描述質(zhì)點(diǎn)位移的運(yùn)動方程表示為

方程在x30 滿足的邊界條件為

同理,將方程(11)代入方程(8)和方程(10)可得應(yīng)力分量非線性項(xiàng)的具體表達(dá)式:

方程(12)、方程(13)和方程(16)描述了線性表面波的位移運(yùn)動方程和邊界條件問題,方程的解析解在相關(guān)文獻(xiàn)[20]中描述如下:

二次勢函數(shù)諧波分量是基波解的修正,即方程(14)、方程(15)和方程(17)的解.求解之前,首先需獲得僅含勢函數(shù)φ(2)和?(2)的解耦方程,可將線性表面波解方程(22)代入方程(14)和方程(15),接著方程(14)對x1求導(dǎo),方程(15)對x3求導(dǎo),二者相加得到方程(25);方程(14)對x3求導(dǎo),方程(15)對x1求導(dǎo),二者相減得到方程(26):

其中R1-R4表示如下:

其中,κcl/ct為縱橫波聲速比,βl-3-(2A+6B+2C)/(λ+2μ)為縱波非線性系數(shù),βt-(λ+2μ+A/2+B)/μ為橫波非線性系數(shù)[21].方程(25)和方程(26)為非齊次偏微分方程,可使用疊加原理[22]得到特解如下:

至此,方程(22)、方程(31)和方程(32)組成了非線性表面波傳播的準(zhǔn)線性解,可見僅縱波位移勢φ(1)自諧振產(chǎn)生的二次分量(即平方項(xiàng))具有累積效果,這一結(jié)論與非線性表面波超聲中可檢測的二次諧波來自于縱波分量的自相互作用[23]的認(rèn)知相符合.

需要注意到前面求解線性問題時(shí)曾利用線性邊界條件方程(16)確定表面波聲速;而表面波在各向同性材料是非頻散波[16],且微擾法暗含瑞利波波速不變的假設(shè),故不再對非線性邊界條件方程(17)進(jìn)行研究.

2.3 基于微擾解的表面波非線性系數(shù)

忽略方程(4)含下標(biāo)β的項(xiàng)并令α1 即可得到一維非線性縱波的質(zhì)點(diǎn)位移運(yùn)動方程如下:

因此,縱波非線性系數(shù)測量表達(dá)式為

通過對比縱波非線性系數(shù)的定義式及測量表達(dá)式,開展對表面波非線性系數(shù)的研究.在非線性表面波實(shí)驗(yàn)中,測量信號一般是x30 表面的面外位移分量.結(jié)合方程(18)、方程(19)和方程(22)—(27),位移u3可采用縱波勢函數(shù)幅值A(chǔ)φ表示:

式中,右邊第一部分A1代表基波,第二部分A2代表累積效應(yīng)的二次諧波,第三部分代表無累積效應(yīng)的二次諧波.其中Dconst表示如下:

在實(shí)際應(yīng)用中,換能器的接收信號中二次諧波相對較弱,且二次以上的高次諧波更難獲得,所以求得的攝動解仍具有參考意義.此外,需要說明的是,本文僅考慮在一定的傳播距離內(nèi)進(jìn)行測量,即可檢測隨距離增加具有累積效應(yīng)的二次諧波成分.因此可以參考方程(35)給出非線性表面波的測量表達(dá)式:

結(jié)合方程(27)、方程(36)和方程(38),可以得到非線性表面波實(shí)驗(yàn)中測量的非線性系數(shù)(定義式)為

實(shí)際上,消去表達(dá)式中的波數(shù)kr,以及η和ξ中包含的角頻率ω和密度ρ0,可以看到方程(39)與縱波非線性系數(shù)βl相似,僅由材料的二、三階彈性常數(shù)確定.

3 仿真分析

本節(jié)基于COMSOL MultiphysicsR有限元軟件驗(yàn)證第2 節(jié)中得到的非線性表面波解析解,主要應(yīng)用該軟件中“Murnaghan material”材料模型探究材料非線性與聲學(xué)非線性問題.該模型需要指定材料二、三階彈性常數(shù)數(shù)值以確定材料應(yīng)力應(yīng)變關(guān)系,與前述理論推導(dǎo)得到非線性表面波僅由材料的二、三階彈性常數(shù)確定的結(jié)論相符,且已有研究證實(shí)了利用該軟件驗(yàn)證非線性超聲理論問題的可行性[24].

3.1 表面波有限元模型

忽略表面波傳播過程中的衰減和衍射影響,即建立非衰減材料中的二維有限元模型,如圖2 所示,選用純鋁作為仿真對象,具體材料參數(shù)如表1 所列.為了增強(qiáng)非線性聲波的仿真效果,提高信噪比,仿真中將材料三階彈性常數(shù)設(shè)置為表1 中材料三階彈性常數(shù)的10倍[25](即βl173.699,βr292.197).根據(jù)方程(24)計(jì)算得到線性表面波聲速為2892.22 m/s,采用位移幅值為10 nm,頻率為2 MHz 的8 個(gè)周期正弦波信號作為點(diǎn)聲源.模型寬度為40 mm,點(diǎn)源設(shè)置在左端點(diǎn)10 mm 處,避免將邊界網(wǎng)格奇點(diǎn)作為激發(fā)源,同時(shí)降低反向傳播的聲波的干擾,且能保證足夠的正向傳播距離;模型厚度大于表面波傳播距離的2 倍,以消除反射的縱/橫波的影響.求解區(qū)域材料屬性設(shè)置為超彈性,采用直接求解器中速度最快的Pardiso 求解器求解.最后,在表面點(diǎn)源后以1 mm 間隔放置探針記錄x3方向位移,采用FFT 算法將記錄的時(shí)域位移信號轉(zhuǎn)換為頻域信號進(jìn)行分析.

表1 純鋁和三種鋁合金的密度和彈性常數(shù)Table 1.Densities and elastic coefficients of pure aluminum and three aluminum alloys.

圖2 非線性表面波二維有限元模型示意圖Fig.2.Schematic diagram of two-dimensional (2D)finite element model for nonlinear surface wave.

為保證計(jì)算精度和可靠性,將近表面黃色區(qū)域(厚度大于4 倍波長)最大網(wǎng)格尺寸dmax設(shè)置為最大解析頻率4 MHz 時(shí)表面波波長的1/20,即3.615×10—2mm,最小網(wǎng)格尺寸dmin設(shè)置為dmax/2,其余區(qū)域劃分為預(yù)定義極細(xì)自由三角形網(wǎng)格.時(shí)間步長 Δt由關(guān)聯(lián)最小網(wǎng)格尺寸和時(shí)間步長的CFL參數(shù)確定:

CFL 一般取值0.2,則 Δt取值為1.25×10—9s.

3.2 有限元模型檢驗(yàn)

為了驗(yàn)證本文激勵(lì)信號波形、材料參數(shù)和網(wǎng)格時(shí)間步長設(shè)置的有效性,開展非線性縱波傳播數(shù)值模擬的預(yù)驗(yàn)證實(shí)驗(yàn).采用線源和周期性邊界條件,使縱波能以平面波形式傳播.相關(guān)傳播模型如圖3 所示,試塊厚度設(shè)為1 個(gè)縱波波長以節(jié)省計(jì)算量.與3.1 節(jié)相比,仿真材料參數(shù)數(shù)值設(shè)置不變.時(shí)間步長不變,網(wǎng)格參數(shù)按照縱波波長做相應(yīng)調(diào)整,探針布置在試塊中間記錄x1方向位移.

圖3 非線性縱波二維有限元模型示意圖Fig.3.Schematic diagram of 2D finite element model for the nonlinear longitudinal wave.

探針波形信號以及縱波中基波和二次諧波隨傳播距離的變化情況如圖4 所示.基波幅值在傳播過程保持不變,擬合幅值是10 nm,與所模擬縱波為平面波的設(shè)置相符;二次諧波幅值隨距離線性增加,需要說明的是,非0 截距來自數(shù)值離散誤差.

圖4 縱波仿真結(jié)果圖 (a)傳播距離15 mm 處探針的信號圖;(b)基波和二次諧波幅值隨距離變化圖,藍(lán)色點(diǎn)是仿真結(jié)果,黃色實(shí)線是線性擬合結(jié)果Fig.4.Simulation results for the nonlinear longitudinal wave: (a)Typical signal at 15 mm propagation distance;(b)plots of fundamental wave and second harmonic amplitude versus propagation distance,where blue points denote the simulation results,and the yellow solid line denotes the fitting line.

根據(jù)方程(34),二次波幅值-傳播距離的斜率理論系數(shù)應(yīng)為利用仿真結(jié)果擬合值以及βl和kl的理論值可反算出基波幅值,再將其與仿真基波擬合幅值進(jìn)行對比驗(yàn)證.預(yù)測基波幅值為9.98 nm,而仿真結(jié)果為10 nm.縱波超聲非線性具有較為成熟的理論,仿真結(jié)果與理論一致,由此驗(yàn)證了該仿真模型在非線性超聲波傳播特性模擬中的有效性.

3.3 非線性表面波仿真

本小節(jié)對非線性表面波進(jìn)行仿真研究.仿真中得到的非線性表面波的典型波形信號,以及不同位置處基波和二次波面外分量幅值隨傳播距離的變化如圖5 所示,二次波幅值隨傳播距離增加近似線性增長.

圖5 表面波仿真結(jié)果 (a)傳播距離15 mm 處探針的信號圖;(b)基波和二次諧波幅值隨距離變化圖,藍(lán)色點(diǎn)是仿真結(jié)果,黃色實(shí)線是擬合結(jié)果Fig.5.Simulation results for the nonlinear surface wave: (a)Typical signal at 15 mm propagation distance;(b)plots of fundamental wave and second harmonic amplitude versus propagation distance,where blue points denote the simulation results,and the yellow solid line denotes the fitting line.

首先,根據(jù)方程(41)表示的波形重疊法確定不同距離點(diǎn)信號的傳播時(shí)間,計(jì)算得到非線性表面波聲速[29]:

其中,f1和f2是距離為d的兩個(gè)時(shí)域信號序列.計(jì)算得到非線性聲表面波速度為2892 m/s,驗(yàn)證了2.2 節(jié)關(guān)于表面波傳播聲速不變的假設(shè).同時(shí),通過聲速可以確定圖5(a)所示波包為表面波.假設(shè)波包中的二次諧波僅源于縱波分量的純縱波非線性特性,根據(jù)3.2 節(jié)該條件下計(jì)算所得的二次諧波-距離斜率理論系數(shù)為0.121,與圖5(b)所示有明顯差異,表明該二次諧波應(yīng)來自非線性表面波傳播特性,且該過程與純縱波非線性效應(yīng)區(qū)別較大.最后,根據(jù)方程(38)和方程(39),二次波幅值-傳播距離斜率理論系數(shù)為據(jù)此反算出基波幅值為1.452 nm,與基波擬合幅值相比符合良好.綜上所述,有限元仿真結(jié)果驗(yàn)證了本文提出的非線性表面波理論的準(zhǔn)確性.

4 非線性表面波特性的數(shù)值計(jì)算與分析

4.1 諧波幅值的數(shù)值結(jié)果

首先利用微擾解的數(shù)值結(jié)果分析非線性表面波諧波傳播特性,并與非線性縱波諧波對比.根據(jù)方程(36)和方程(34),采用表1 中純鋁的材料參數(shù),分別計(jì)算不同初始條件和傳播距離的表面波面外分量,以及同等條件下縱波傳播時(shí)二次諧波幅值,結(jié)果如表2 所列.

表2 初始條件和傳播距離對面外分量和縱波二次諧波幅值的影響Table 2.The effects of initial conditions and propagation distance on the amplitudes of out-of-plane component and the second harmonic longitudinal wave.

由表2 可知,累積二次表面波幅值A(chǔ)2和縱波二次幅值均隨基波幅值、頻率或傳播距離的變化呈正相關(guān),具有相同的傳播特性.然而,在相同條件下,由于表面波波數(shù)較大,A2幅值會比高一個(gè)數(shù)量級;在實(shí)際中,低衰減和能量集中等特性將使二次表面波可測量性更強(qiáng).因此,非線性表面波技術(shù)相比縱波具有更高的信噪比.非積累二次諧波幅值幅值比A2小兩個(gè)數(shù)量級,且可能包含邊界耦合非線性[30]和微擾法誤差,因此實(shí)際測量中可以忽略.

4.2 非線性系數(shù)的數(shù)值結(jié)果

本小節(jié)將根據(jù)表1 中的材料參數(shù),計(jì)算縱波和表面波非線性系數(shù)數(shù)值結(jié)果并進(jìn)行對比分析.在此之前,首先列出當(dāng)前學(xué)界應(yīng)用較多的幾種表面波非線性系數(shù)定義,并納入比較范圍.Herrmann等[17]基于非線性表面波縱波分量的近似推導(dǎo)得出:

結(jié)合方程(38),其定義的表面波非線性系數(shù)可表示為

Shull等[31]基于Zabolotskaya[16]的理論體系,以質(zhì)點(diǎn)速度諧波幅值vn(n1,2,···)作為研究變量,參考流體情況,定義表面波非線性系數(shù)β11為

參數(shù)S11和ζ的詳細(xì)表達(dá)式見文獻(xiàn)[31],vn和有如下關(guān)系[16]:

Masurkar和Tse[18]對二階非線性應(yīng)力應(yīng)變關(guān)系做進(jìn)一步近似,提出一種基于諧波實(shí)際測量幅值的歸一化非線性系數(shù)δR.該系數(shù)實(shí)際是將信號按基波幅值歸一化,表示如下:

式中F(kr,η,ξ)是特征函數(shù),具體表達(dá)式見文獻(xiàn)[18].

將表1 中4 種金屬的材料參數(shù)代入方程(39)、方程(43)、方程(46)和方程(47)中計(jì)算上述相關(guān)的非線性系數(shù),結(jié)果如表3 所列.

表3 表1 中4 種材料的非線性系數(shù)Table 3.Nonlinear parameters of four materials listed in Table 1.

由表3 可知,βr數(shù)值比βl大,而的數(shù)值最小.這是由于推導(dǎo)過程主要是基于的質(zhì)點(diǎn)速度慢變幅值,且其推導(dǎo)中存在近似,不適用于實(shí)際測量對象是面外位移分量的情況.δR和βr數(shù)值相近,但對比方程(47)和方程(39),不考慮衰減和衍射效應(yīng)時(shí),δR值會隨初始條件和傳播距離的變化而變化;反之,當(dāng)衰減或衍射影響成為主導(dǎo)因素時(shí),其數(shù)值可能會出現(xiàn)單個(gè)或多個(gè)局部峰值,該現(xiàn)象還有待深入理論研究.此外,由方程(39)可知,表面波非線性系數(shù)與縱波非線性系數(shù)存在確定數(shù)值關(guān)系,與材料的二階彈性常數(shù)相關(guān),隨材料變化而改變.綜上,結(jié)果表明,非線性表面波傳播特性與縱波相似,但在初始條件相同時(shí)表面波諧波幅值比縱波諧波幅值高,因此在實(shí)際測量中信噪比更高.

5 總結(jié)

本文首先采用微擾法研究了準(zhǔn)線性條件下材料中的非線性表面波傳播特性,引入勢函數(shù)描述非線性表面波的波動控制方程,推導(dǎo)了基波和二次諧波解,并得到了表面波非線性系數(shù)定義式.結(jié)果顯示,由縱波位移勢自諧振產(chǎn)生的累積項(xiàng)幅值比縱波的二次諧波幅值大,而且具有能量集中、衰減低等特性,使非線性表面波檢測技術(shù)在實(shí)際應(yīng)用中比非線性縱波具有更高的信噪比;由橫波分量引起的諧波項(xiàng)幅值沒有累積效應(yīng),因此對表面波二次諧波的影響可以忽略.二維表面波有限元模型位移幅值仿真結(jié)果與理論值符合良好,證實(shí)了本文理論推導(dǎo)過程的有效性以及所得諧波微擾解的準(zhǔn)確性.最后,根據(jù)微擾解的數(shù)值結(jié)果探究了當(dāng)前幾種表面波相關(guān)的非線性系數(shù)之間的聯(lián)系,對比計(jì)算結(jié)果可知表面波非線性系數(shù)和縱波非線性系數(shù)存在正比關(guān)系,并且該系數(shù)只與材料二次彈性常數(shù)有關(guān).目前僅限于理論推導(dǎo)和數(shù)值模擬分析,后續(xù)將繼續(xù)研究實(shí)際測量相對于當(dāng)前理想情況的修正理論,從而研發(fā)一套可靠的測量設(shè)備及方法,實(shí)現(xiàn)利用表面波絕對非線性系數(shù)評價(jià)材料狀態(tài).