仲曉敏，韓慶邦，蔣 謇，黃建偉，吳 寧

(河海大學(xué)物聯(lián)網(wǎng)工程學(xué)院，江蘇 常州 213022)

0 引 言

近年來，基于超聲的非線性效應(yīng)，如諧波的產(chǎn)生、波束混疊以及諧振頻率的改變而發(fā)展的非線性超聲無損檢測[1]和評價方法，對于克服傳統(tǒng)超聲的不足具有積極的意義，是對傳統(tǒng)線性超聲檢測[2]手段的有效補充。國內(nèi)外很多學(xué)者在這方面做了大量的研究，XIANG Y X等[3]對非線性超聲導(dǎo)波檢測技術(shù)進行了研究，王成會等[4]研究了超聲波作用下氣泡的非線性振動，SHUI Y等[5]研究了固體板中的波的非線性傳播。當(dāng)波與接觸界面相互作用時，產(chǎn)生二階及高階諧波，這就構(gòu)成了一個相對較新的研究領(lǐng)域，即接觸聲非線性[6-8]。與傳統(tǒng)材料的非線性[9]相比，接觸聲非線性會產(chǎn)生高次諧波，能提供更有價值的接觸狀態(tài)信息，可用于檢測部分閉合缺陷的材料[10]。

本文基于BIWA S[11]研究的一維基頻簡諧波在固-固界面的傳播特性以及RICHARDSON[12]理論，引入包含基波和諧波的復(fù)合入射波，從界面剛度特性和其它聲學(xué)參數(shù)出發(fā)，建立一個一維的非線性接觸剛度模型，推導(dǎo)出反射波、透射波的表達式，分析反射系數(shù)、透射系數(shù)和諧波的幅值比，詳細討論了諧波傳播時的接觸剛度參數(shù)與接觸應(yīng)力之間的關(guān)系。

1 理論建模

1.1 理論模型

建立如圖1所示的波傳播模型，一維彈性縱波沿x軸傳播，接觸面的參考平面位置分別設(shè)置為x=X-,x=X+，參考平面之間的間隙距離為X+-X-。假設(shè)在區(qū)域x＜X-,x＞X+內(nèi)兩個線性彈性固體具有相同的屬性(密度ρ和縱向剛度E)，固-固界面有一定的粗糙度。在區(qū)域x＜X-，入射一維縱波finc，得到反射波fref，在區(qū)域x＞X+，得到透射波ftra。

圖1 一維波通過接觸面?zhèn)鞑サ氖疽鈭DFig.1 One-dimensional wave propagation through contact interface

在沒有波入射的平衡狀態(tài)下，設(shè)兩個固體的接觸靜應(yīng)力為σ0，間隙距離為h0。一維縱波沿x軸正向傳播，運動方程以及應(yīng)力位移關(guān)系為[11]

其中，u(x,t)表示聲波沿x方向傳播的位移，σ(x,t)表示應(yīng)力變化，t表示時間。方程(1)、(2)的解為

其中，函數(shù)finc(x-ct)、fref(x+ct)及ftra(x-ct)分別表示入射波、反射波和透射波，聲波波速c=(λ/ρ)1/2，λ為彈性常數(shù)。

當(dāng)波與接觸界面相互作用時，間隙距離隨時間而變化，即：

應(yīng)力滿足：

即應(yīng)力在界面上是連續(xù)的，并且服從于間隙距離h與接觸應(yīng)力σ(h)的非線性函數(shù)關(guān)系。

1.2 反射波與透射波

為求解上述方程，引入以下變量Z(t)、Y(t)[12]：

其中，Z代表界面兩端位移和的一半，Y表示界面兩端位移差，對式(8)、(9)求導(dǎo)并推導(dǎo)[12]得：

經(jīng)一系列推導(dǎo)，透射波和反射波的表達式為

2 復(fù)合諧波入射特性分析

2.1 攝動分析

這里引入包含基波和二次諧波的復(fù)合波入射，表達式為

其中，A和B表示為基波及諧波的振幅，ω為角頻率。考慮到波的位移較小時，間隙距離改變也很小，利用泰勒展開：

其中：K1表示線性剛度，表示接觸面非線性剛度，(二階剛度)。

將式(14)、 (15)代入式(11)得關(guān)于Y的方程：

考慮方程的近似解[13]Y=Y1+Y2，Y2是二階微量，且Y1、Y2滿足以下方程：

解得：

其中，δ3=arctan[a/(4ω)]，?1=arctan[a/(2ω)]，δ4=θ1-δ1-δ2，δ5=θ2+δ1-δ2，θ1=arctan[a/(3ω)]，θ2=arctan(a/ω)，系數(shù)D、E見附錄I。

透射和反射波的表達式分別表示為

其中，M、N分別為透射波中一倍頻分量、二倍頻分量的振幅，Ι、L分別為反射波中一倍頻分量、二倍頻分量的振幅，表達式見附錄I，相位Δ1,Δ2,Δ3,Δ4見附錄I。

從式(20)和式(21)可以清楚地看到，透射波和反射波都出現(xiàn)了角頻率3ω、4ω，這是復(fù)合波在界面處由于非線性相互作用產(chǎn)生的高次諧波。

2.2 反射、透射系數(shù)與非線性特性

定義透射、反射波基頻分量的絕對振幅與入射波振幅的比值，分別為

其中，反射系數(shù)R1、R2分別為一倍頻和二倍頻下，反射波中振幅與入射波對應(yīng)振幅之比，透射系數(shù)T1、T2定義同上，R、T沒有非線性效應(yīng)。

為研究非線性效應(yīng)，定義非線性系數(shù)β1、γ1是透射波及反射波中三次諧波分量振幅與基波分量振幅的比值，β2、γ2是透射波及反射波中四次諧波與二倍頻分量振幅的比值。經(jīng)推導(dǎo)可得：

DRINKWATER等[14]對各種粗糙表面進行建模，采用超聲波入射實驗獲得剛度K1與接觸應(yīng)力σ0平方根成比例的結(jié)論，將粗糙表面擬合成平滑表面，即線性剛度與應(yīng)力關(guān)系在界面上是連續(xù)的，假定一個簡單的應(yīng)力與剛度的冪律關(guān)系，其表達式為

其中，C和m是正常數(shù)。將和與式(25)相結(jié)合，給出二階剛度與接觸壓力的函數(shù)，即：

3 數(shù)值結(jié)果及分析

針對有限范圍內(nèi)的粗糙接觸面，傳播介質(zhì)為金屬鋁，ρ為2 700 kg·m-3，速度c為6 420 m·s-1，對冪律函數(shù)進行計算分析，選取常數(shù)C為6×1010Pa-1/2m-1，常數(shù)m為0.5。分析相關(guān)系數(shù)時，頻率選取5、10、15 MHz三種情況，入射波振幅選取5 nm和20 nm進行對比，以上參數(shù)選自文獻[11]。

圖2為不同頻率下反射系數(shù)與壓力的關(guān)系圖。從圖2(a)、2(b)中可以看出，反射系數(shù)R1、R2隨著頻率的增大而增大，隨接觸壓力的增大而減小。

圖3為不同頻率下透射系數(shù)與壓力的關(guān)系圖。從圖3(a)、3(b)可以看出，透射系數(shù)T1、T2隨著頻率的增大而減小，隨接觸壓力的增大而增大。

圖2 振幅為20 nm時，不同頻率下，反射系數(shù)與壓力的關(guān)系圖Fig.2 The relationship between reflection coefficient and contact pressure at different frequencies (the amplitude of incident wave is 20 nm)

圖3 振幅為20 nm時，不同頻率下，透射系數(shù)與壓力的關(guān)系圖Fig.3 The relationship between transmission coefficient and contact pressure at different frequencies (the amplitude of incident wave is 20 nm)

下面選擇不同的入射波振幅和基波頻率分析非線性系數(shù)。選取入射波振幅分別為5 nm和20 nm，非線性系數(shù)與接觸壓力的關(guān)系如圖4所示。從圖4可以看出，在同一頻率下，當(dāng)入射波振幅選取20nm時，非線性系數(shù)更大，意味著在振幅較大的情況下，非線性效應(yīng)更顯著。

在入射波振幅為20 nm的情況下，入射頻率分別為5、10、15 MHz時，透射波的非線性系數(shù)與接觸壓力的關(guān)系如圖5所示。

從圖5(a)、5(b)可以看出，非線性系數(shù)β1，β2隨頻率的增大而增大，隨接觸壓力增大而減小。

圖4 頻率為5 MHz時，不同振幅下非線性系數(shù)與接觸壓力的關(guān)系圖Fig.4 The nonlinear coefficients (β1,β2 and γ1,γ2)under different amplitudes of incident wave versus the contact pressure at the frequency of 5 MHz

圖5 入射波振幅為20 nm時，不同頻率下透射波的非線性系數(shù)與接觸壓力的關(guān)系圖Fig.5 Relationship between nonlinear coefficients of transmission wave (β1and,β2)and contact pressure at different frequencies (the amplitude of incident wave is 20 nm)

在入射波振幅為20 nm時，入射頻率分別為 5、10、15 MHz時，反射波的非線性系數(shù)與接觸壓力的關(guān)系如圖6所示。從圖6可以看出，非線性系數(shù)γ1，γ2隨頻率增大而減小，隨接觸壓力增大而減小。

在入射波振幅為20 nm，入射頻率為5 MHz時，透射波與反射波中三次諧波與四次諧波非線性系數(shù)對比如圖7。從圖7可以看出，三次諧波比四次諧波的非線性系數(shù)更大，即三次諧波的非線性效應(yīng)更好。

圖6 入射波振幅為20 nm時，不同頻率下反射波的非線性系數(shù)與接觸壓力的關(guān)系圖Fig.6 Relationship between nonlinear coefficients of reflection wave (γ1and,γ2)and contact pressure at different frequencies (the amplitude of incident wave is 20 nm)

圖7 透射波和反射波中三次諧波與四次諧波非線性系數(shù)的對比圖Fig.7 Comparison between nonlinear coefficients of the third and fourth harmonics of transmission and reflection waves

4 結(jié) 論

本文基于固-固界面的線性剛度的冪律模型，計算一維入射復(fù)合波的反射系數(shù)、透射系數(shù)以及非線性系數(shù)，分析了接觸壓力和頻率對它們的影響，發(fā)現(xiàn)反射系數(shù)隨頻率的增大而增大，隨接觸壓力的增大而減小，而透射系數(shù)規(guī)律則相反；對比不同振幅下的非線性系數(shù)，發(fā)現(xiàn)選取較大振幅時，非線性現(xiàn)象更顯著。在透射波、反射波中，非線性系數(shù)隨接觸壓力的增大而減?。煌干洳ǚ蔷€性系數(shù)隨頻率的增大而增大，但變化幅度不明顯，在反射波中情況相反；在相同振幅、頻率下，三次諧波的非線性效應(yīng)優(yōu)于四次諧波。綜合考慮，反射波中三次諧波分量在檢測微小裂縫中的優(yōu)勢最大。

附錄I

Y2的表達式中的系數(shù)：

透射波表達式中的系數(shù)：

反射波表達式中的系數(shù)：

透射波表達式中的相位：

反射波表達式中的相位：

附錄II

式(11)的推導(dǎo)如下：

根據(jù)參考文獻[13]，可知：

化簡式(a)、(b)得：

由式(c)～(d)得：