分層半空間表面非線性瑞利波的激發(fā)?

賈 璐 閻守國 張碧星 黃 娟

(1 中國科學(xué)院聲學(xué)研究所 聲場聲信息國家重點(diǎn)實(shí)驗(yàn)室 北京 100190)

(2 中國科學(xué)院大學(xué) 北京 100049)

0 引言

分層結(jié)構(gòu)由于其具有高強(qiáng)度、高剛度等優(yōu)點(diǎn)被廣泛應(yīng)用于各類實(shí)際的工程領(lǐng)域。但是分層結(jié)構(gòu)的近表面容易率先出現(xiàn)材料性能退化和微裂紋等損傷，這可能導(dǎo)致構(gòu)件嚴(yán)重的損壞，甚至造成災(zāi)難性事故。因此，研究一種無損評估(Non-destructive examination, NDE)方法以盡早發(fā)現(xiàn)缺陷對于確保分層結(jié)構(gòu)的安全性和完整性非常重要。相比于線性超聲波檢測方法，固體結(jié)構(gòu)中高階彈性常數(shù)對材料性能的變化更加敏感，高階彈性常數(shù)的改變將導(dǎo)致聲波中高次諧波發(fā)生變化，因此可以通過對高次諧波的測量反映材料的性能。Landau 等[1]建立的固體中非線性聲波理論為非線性超聲檢測技術(shù)的研究奠定理論基礎(chǔ)。

非線性聲學(xué)根據(jù)高次諧波的來源可分為經(jīng)典非線性聲學(xué)和接觸非線性聲學(xué)，其中接觸非線性主要指結(jié)構(gòu)中接觸類微損傷與聲波相互作用產(chǎn)生的非線性特征，國內(nèi)的劉曉宙等[2?5]對其做了大量研究。另一類經(jīng)典非線性問題則與材料晶格的非簡諧效應(yīng)有關(guān)，本文討論的問題是基于經(jīng)典非線性聲學(xué)理論。Cantrel 等[6]由離散晶格模型建立各向同性介質(zhì)中的非線性彈性波方程，通過攝動法對一維非線性縱波進(jìn)行求解，得到諧波幅值與非線性系數(shù)的關(guān)系，并展開相關(guān)實(shí)驗(yàn)研究。張世功等[7]利用有限元仿真和實(shí)驗(yàn)研究分析一維非線性聲波傳播過程，提出具有一定物理意義的二次諧波隨傳播距離變化的數(shù)學(xué)關(guān)系。稅國雙等[8]運(yùn)用非線性縱波對列車外圓彈簧的疲勞損傷進(jìn)行實(shí)驗(yàn)檢測，初步實(shí)現(xiàn)非線性超聲對特殊結(jié)構(gòu)的檢測。錢祖文[9]推導(dǎo)了彈性介質(zhì)中二階勢函數(shù)的波動方程，并研究各向同性介質(zhì)中縱波和橫波非線性作用。在非線性瑞利波的研究方面，Zaboloskaya 等[10?11]基于哈密頓理論公式對各向同性固體介質(zhì)表面?zhèn)鞑サ姆蔷€性瑞利波進(jìn)行理論研究，并通過數(shù)值方法驗(yàn)證。Herrmann等[12]基于聲表面波的縱波分量推導(dǎo)出瑞利波非線性系數(shù)，并且提出一種測量金屬模型中瑞利波二次諧波的方法。然而這些體波和非頻域表面波理論部分大多近似為一維非線性縱波，并且局限于均勻單一結(jié)構(gòu)中。近年來，頻散導(dǎo)波由于其能量大、衰減小等優(yōu)勢，引起了學(xué)者廣泛研究。Lima 等[13]采用微擾法和互易定理探究板中Lamb波二次諧波產(chǎn)生條件和傳播規(guī)律；鄧明晰等[14?16]分別通過界面非線性聲波發(fā)射法和導(dǎo)波模式展開法對平面固體結(jié)構(gòu)中導(dǎo)波的非線性問題進(jìn)行研究；Liu 等[17]提出了產(chǎn)生具有較強(qiáng)累積效應(yīng)二次諧波的基波的選擇標(biāo)準(zhǔn)；Xiang 等[18]對Lamb 波在復(fù)合金屬材料中傳播時產(chǎn)生的二次諧波問題進(jìn)行了相關(guān)的理論和試驗(yàn)研究。

當(dāng)分層結(jié)構(gòu)的總厚度遠(yuǎn)大于檢測聲波波長時，認(rèn)為在結(jié)構(gòu)表面?zhèn)鞑サ氖侨鹄?，這種分層結(jié)構(gòu)中的瑞利波區(qū)別于經(jīng)典瑞利波，具有頻散和多模的特性。通過對瑞利波的頻散特性進(jìn)行反演可以獲得介質(zhì)信息，瑞利波的這一特性已經(jīng)在地質(zhì)工程和無損檢測中得到廣泛應(yīng)用，但對于非線性的頻散瑞利波的研究，由其激發(fā)和傳播機(jī)制的復(fù)雜性，在以往的研究中很少涉及。

本文針對均勻分層半空間結(jié)構(gòu)中瑞利波二次諧波的激發(fā)和傳播規(guī)律進(jìn)行研究。采用二階微擾近似法和模態(tài)分解得到結(jié)構(gòu)中二次諧波的位移解析式，根據(jù)瑞利波頻散曲線得到基波與二倍頻波相速度相等的匹配點(diǎn)，針對不同分層結(jié)構(gòu)中的匹配模式的二次諧波傳播特性進(jìn)行理論分析和數(shù)值計(jì)算，為實(shí)際檢測中選擇合適的聲源頻率和瑞利波模式提供理論依據(jù)。

1 理論基礎(chǔ)

1.1 非線性波動方程

建立二維直角坐標(biāo)系Oxz，各向同性N層半空間結(jié)構(gòu)分布在z≥0范圍內(nèi)，z <0為真空，如圖1所示，z= 0 表示第一層介質(zhì)的上表面，層與層之間滿足應(yīng)力和位移連續(xù)的邊界條件。自由表面條件下的非線性波動方程和邊界條件分別為[1]

圖1 分層半空間示意圖Fig.1 The draft of a layered half-space

其中，u表示位移矢量，nz是z方向上的單位向量，ρ為無形變時固體的密度，λ和μ表示拉梅系數(shù)，A、B和C代表landau 形式的三階彈性系數(shù)，PL為第一階Piola-Kirchhoff應(yīng)力張量的線性項(xiàng)，PNL為第一階Piola-Kirchhoff應(yīng)力張量的非線性項(xiàng)，F(xiàn)是第一階Piola-Kirchhoff應(yīng)力張量非線性項(xiàng)的散度，它們分別表示為

其中，

腳標(biāo)i、j、k、l、m均取1、2、3 且滿足愛因斯坦求和約定。

1.2 非線性方程的求解

在弱非線性條件下，利用微擾法，將質(zhì)點(diǎn)位移矢量u(x,z,t)近似為

其中，u(0)(x,z,t)和u(1)(x,z,t)分別對應(yīng)聲場中的基波和二次諧波，并且二者滿足如下關(guān)系：

將(4)式代入式(1)后可得以下線性方程：

其中，F(xiàn)(u(0))表示在固體非線性條件下由基波位移產(chǎn)生的驅(qū)動力。

由式(6a)結(jié)構(gòu)中頻率為ω的第l階導(dǎo)波模式位移可寫為如下形式：

其中k=ω/cl為波數(shù)，cl表示第l階導(dǎo)波模式的傳播速度，后文中u(l0)=u(l0)(z)。

根據(jù)Auld 理論[19]，將式(7a)中二次諧波聲場表示為角頻率為2ω時導(dǎo)波模式的線性疊加：

其中，(z)表示線性條件中頻率2ω處的第m階導(dǎo)波模式位移在z方向的分布，其定義與式(8)中的(z)相同，體現(xiàn)的是某一特定結(jié)構(gòu)下線性導(dǎo)波模式的固有特性，Am(x)是微擾條件下由基波振動產(chǎn)生的二次諧波模式的振幅，并且其值與傳播距離x有關(guān)。根據(jù)互易定理，第m階導(dǎo)波模式對應(yīng)的Am可表示為如下形式[19]：

其中，

其中，是頻率為2ω第m階導(dǎo)波模式的質(zhì)點(diǎn)速度，zj+和zj?表示第j層介質(zhì)的上表面和下表面，nx為x方向上的單位向量，?表示取共軛，Pmm是導(dǎo)波模式沿x正方向的平均功率流(沿y方向取單位長度)，和分別為表面牽引力和體積力所提供的面源和體源。值得注意的是，在層狀半間結(jié)構(gòu)中，最后一層介質(zhì)為無限大的半空間結(jié)構(gòu)，此時式(11a)～(11c)的積分下限為負(fù)無窮大，與板中導(dǎo)波的展開形式不同。

根據(jù)式(10)可得當(dāng)導(dǎo)波模式滿足如下條件時，將產(chǎn)生隨傳播距離具有積累效應(yīng)的二次諧波：(1)基頻波和某一模式二倍頻波相速度匹配，即km= 2kl；(2)基波和二倍頻波之間存在非零的能量流傳遞，即+≠0。在數(shù)值計(jì)算中和的幅值均大于等于零，因此僅討論考慮和不同時為零的情況。已有研究表明當(dāng)導(dǎo)波模式在深度方向上具有對稱性時，存在== 0 的情況[20]，而瑞利波在分層結(jié)構(gòu)中傳播時其幅值在深度方向上逐漸衰減、不具有對稱性[21]，即體驅(qū)動力和面驅(qū)動力均不為零。因此在理論上滿足相速度匹配的點(diǎn)，均能產(chǎn)生具有累積效應(yīng)的二次諧波，而當(dāng)km≠ 2kl時二次諧波不具備積累效應(yīng)，幅度將隨傳播距離周期震蕩。需要說明的是，本文采用攝動近似法對非線性瑞利波的傳播特性進(jìn)行理論分析時僅考慮了基波和二次諧波，忽略了二次諧波以上的高次諧波，得到以上結(jié)論，并未考慮實(shí)際傳播過程中二次諧波能量向更高次諧波傳遞以及聲波在介質(zhì)中傳播時的衰減，后續(xù)研究中將結(jié)合實(shí)驗(yàn)進(jìn)行詳細(xì)討論。

由于不具備距離積累效應(yīng)的二次諧波很難被觀察到，因此通常只關(guān)心相速度匹配條件下瑞利波二次諧波的特性，定義參數(shù)βm表示匹配條件(km=2kl)下，第l階基波模式產(chǎn)生的瑞利波二次諧波強(qiáng)度：

2 數(shù)值計(jì)算

本節(jié)以實(shí)際應(yīng)用中使用較為廣泛的雙層半空間和含有低速夾層三層半空間結(jié)構(gòu)為例，數(shù)值計(jì)算并分析二次諧波傳播規(guī)律，選擇出模型中適用于實(shí)際檢測的聲源頻率和基波模式。分層結(jié)構(gòu)中的材料相關(guān)參數(shù)[10,22]如表1所示。

表1 材料參數(shù)Table 1 Parameters of materials

2.1 雙層半空間結(jié)構(gòu)

建立由環(huán)氧樹脂和鋼復(fù)合的雙層半空間結(jié)構(gòu)模型，第一層介質(zhì)為環(huán)氧樹脂其厚度h= 0.01 m。根據(jù)瑞利波頻散方程[21]繪制該模型的頻散曲線如圖2 所示，S 系列和s系列分別表示的基頻瑞利波和二倍頻瑞利波的頻散曲線，它們的交點(diǎn)即為相速度匹配點(diǎn)。圖中計(jì)算βm時，假設(shè)基階頻散曲線上的每一點(diǎn)都存在匹配模式(匹配模式的速度與基階模式相等，頻率為基頻的2 倍)，根據(jù)公式(10)，此時匹配模式只能為其自身，即m=l，最終得到與基頻瑞利波頻散曲線相對應(yīng)的βm，對其取對數(shù)后再除以lgβm中的最大值進(jìn)行歸一化處理，并以顏色的變化在頻散曲線上表示。圖2 中相速度匹配點(diǎn)處的βm則反映實(shí)際匹配條件下對應(yīng)基波模式能產(chǎn)生具有線性累積效應(yīng)二次諧波分量的強(qiáng)度。

根據(jù)(9)式可知，伴隨基頻瑞利波傳播時所發(fā)生的二次諧波聲場可視為一系列二倍頻瑞利波疊加而成，如點(diǎn)P對應(yīng)頻率處，豎直虛線L與二倍頻瑞利波頻散曲線存在一系列交點(diǎn)，這些交點(diǎn)對應(yīng)的二倍頻瑞利波模式的疊加即表示伴隨S2 模式傳播過程中所發(fā)生的二次諧波聲場。根據(jù)理論部分的分析可知，相速度匹配的模式均能激發(fā)出具有累積效應(yīng)的二次諧波，如圖2 所示該模型中P點(diǎn)和Q點(diǎn)對應(yīng)的參數(shù)βm較大，因此選擇這兩點(diǎn)對應(yīng)的頻率和基波模式進(jìn)行理論計(jì)算，同時為了對比結(jié)果，同時計(jì)算βm值較小的R點(diǎn)對應(yīng)模式產(chǎn)生的二次諧波分量，具體參數(shù)如表2 所示。需要說明的是，由于實(shí)際檢測中通常只獲得法向位移信號，因此本文僅對法向位移分量進(jìn)行研究。

圖2 雙層半空間中瑞利波頻散曲線及參數(shù)βm 分布圖Fig.2 The dispersion curve of Rayleigh waves and distribution of βm in a two-layer half-space

表2 匹配點(diǎn)對應(yīng)基波模式的具體參數(shù)Table 2 Parameters of fundamental modes at matching points

匹配點(diǎn)R、P和Q對應(yīng)二次諧波分量的法向位移隨傳播距離變化關(guān)系如圖3 所示。根據(jù)數(shù)值計(jì)算結(jié)果，在任意匹配點(diǎn)相速度匹配模式所對應(yīng)的二次諧波分量，隨傳播距離的增加而線性增長，而其余分量則沿傳播方向周期震蕩傳播。由圖3(c)可知，盡管匹配點(diǎn)Q對應(yīng)的參數(shù)βm最大，但同時該模式產(chǎn)生的周期震蕩二次諧波分量的幅值也較大，在一定傳播范圍內(nèi)對二次諧波的影響無法忽略，因此在匹配點(diǎn)Q對應(yīng)的基波第S4 階模式無法產(chǎn)生具有明顯線性增長的二次諧波，并且高頻信號將激發(fā)多個導(dǎo)波模式，使得聲波在結(jié)構(gòu)中傳播更為復(fù)雜，不利于聲波信號的分析。反之在匹配點(diǎn)R、P中，各自匹配模式所對應(yīng)的二次諧波分量不僅隨傳播距離x的增加而線性增長，而且其幅值遠(yuǎn)大于其他二次諧波分量的幅值，可忽略其他模式，此時伴隨基波第S1 階和第S2 階模式產(chǎn)生的二次諧波聲場隨傳播距離呈現(xiàn)明顯的線性增長關(guān)系。最后通過對比兩點(diǎn)的二次諧波幅值，可知在該雙層半空間模型中，頻率fP處基頻波第S2 階模式所產(chǎn)生的二次諧波適用于后續(xù)的實(shí)驗(yàn)中。

圖3 雙層半空間中匹配點(diǎn)對應(yīng)二次諧波分量的法向位移幅值隨傳播距離變化關(guān)系圖Fig.3 The influence of propagation distance on the normal displacement of second-harmonic components at matching points in a two-layer half-space

2.2 含有低速夾層的三層半空間結(jié)構(gòu)

含有低速夾層是指中間介質(zhì)的橫波速度低于上下層的橫波速度，數(shù)值計(jì)算中該結(jié)構(gòu)的材料從上到下依次為鋁、環(huán)氧樹脂和鋼，其中第一層介質(zhì)和第二次介質(zhì)厚度分別為0.03 m 和0.01 m，材料參數(shù)如表1所示。從頻散曲線圖4可知，該模型中同樣存在多個相速度匹配的頻率點(diǎn)，根據(jù)參數(shù)βm數(shù)值，選擇R點(diǎn)和P點(diǎn)對應(yīng)頻率和基波模式計(jì)算相應(yīng)二次諧波分量如圖5 所示。頻率fR、fP中與基頻波相速度匹配的二次諧波分量均隨傳播距離線性增長，但對比圖5(a)和圖5(b)可知在相同傳播距離內(nèi)匹配點(diǎn)P中線性增長的二次諧波分量幅值大，并且非匹配模型對二次諧波聲場的影響在較短傳播距離內(nèi)便可忽略不計(jì)，因此可以激勵聲源頻率為fR的第S2階基波模式對此結(jié)構(gòu)進(jìn)行非線性超聲檢測工作。

圖4 含低速夾層三層半空間中瑞利波頻散曲線及參數(shù)βm 分布圖Fig.4 The dispersion curve of Rayleigh waves and distribution of βm in a three-layer half-space with low-velocity layer

圖5 含低速夾層的三層半空間中匹配點(diǎn)對應(yīng)二次諧波分量的法向位移幅值隨傳播距離變化關(guān)系圖Fig.5 The influence of propagation distance on the normal displacement of second-harmonic components at matching points in a three-layer halfspace with low-velocity layer

3 結(jié)論

對于瑞利波傳播過程中的非線性二次諧波問題，首先由攝動理論得到零階和一階位移即基波和二次諧波分別滿足的運(yùn)動方程和邊界條件，再將一階位移表示為一系列二倍頻瑞利波模式的疊加，最后根據(jù)互易定理得到一階位移模式振幅表達(dá)式。理論分析和數(shù)值計(jì)算結(jié)果表明，與基波相速度匹配的二倍頻瑞利波模式所對應(yīng)的二次諧波分量在傳播方向上具有累積效應(yīng)，隨傳播距離線性增長；針對分層半空間結(jié)構(gòu)中基頻瑞利波存在多個相速度匹配的模式，甚至某一基頻波模式具有多個相速度匹配點(diǎn)問題，通過計(jì)算參數(shù)βm對激勵聲源頻率和瑞利波模式進(jìn)行選擇和優(yōu)化，為后續(xù)實(shí)際檢測工作提供理論指導(dǎo)。