吳榮興,王曉明,鄭 東
(1.寧波職業(yè)技術學院 建筑工程學院,浙江 寧波 315800; 2.寧波大學 機械與力學學院,浙江 寧波 315211)
各類聲表面波氣體、質量和偏場傳感器已被廣泛應用于各種復雜環(huán)境下的檢測工作[1~3]。主要工作原理是通過聲學和電學參數的測量來推定所處環(huán)境因素或物質特性[3~7]。過去認為聲表面波在液體檢測中由于存在嚴重能量衰減,無法檢測其有效信號,目前證明可以將聲表面波技術應用于微流體的特性檢測[8]。在慣性傳感器領域,理論研究表明哥氏力對聲表面波的傳播有抑制作用,并且存在著對應關系,可以用來檢測旋轉場的旋轉頻率[9,10]。大部分聲表面波陀螺效應傳感器所建的模型都是基于半無限大各向同性彈性體結構,而真實的壓電聲波器件往往采用各向異性的壓電材料組成的層狀板型結構[11,12]。雖然實驗研究結果顯示陀螺效應對聲表面波傳播有著顯著影響,但是實驗樣機的分辨率和精度都無法滿足傳感器要求[13]。因此,完整建立更為真實的壓電聲波陀螺效應傳感器模型和精確分析成為一項迫切的工作。
這里分析了聲表面波在考慮旋轉場作用下的半無限石英晶體和無限大石英晶體板內的傳播特性,獲得了波速、板厚、旋轉軸角度和旋轉頻率之間的對應關系,為實際聲表面波陀螺效應傳感器的研制提供了理論指導。
聲表面波在考慮旋轉場作用下的石英晶體板內傳播如圖1所示。旋轉場的旋轉頻率為Ω,旋轉軸與x1軸的夾角為θ。首先考慮底部是半無限大石英晶體的情況。
圖1 帶旋轉場的石英晶體板
當聲表面波在半無限大石英晶體內沿著x1方向傳播時[14,15],位移和電勢可以假設為
u1=A1ekβx2eik(x1-ct),u2=A2ekβx2eik(x1-ct),
u3=A3ekβx2eik(x1-ct),φ=A4ekβx2eik(x1-ct)
(1)
式中u1(u2,u3),φ,A1(A2,A3,A4),k,β,x1(x2),c和t分別為聲表面波的位移、電勢、振幅、波數、衰減系數、坐標、波速和時間變量。
基于式(1)的位移假設,獲得應變和電場分量為
S1=iku1,S2=kβu2,S3=0,
S4=kβu3,S5=iku3,S6=kβu1+iku2,
E1=-ikφ,E2=-kβφ,E3=0
(2)
式中Sp(p=1,2,…,6)和E1(E2,E3)為各方向的應變和電場分量。
基于應變和電場分量式(2),可以得到應力和電位移分量為
Tp=k[(cp 6β+icp1)u1+(cp 2β+icp 6)u2+
(cp 4β+icp 5)u3+(e2pβ+ie1p)φ],
Dj=k[(ej 6β+iej 1)u1+(ej 2β+iej 6)u2+
(ej 4β+iej 5)u3+(εj 2β+iεj 1)φ],
p=1,2,3,4,5,6,j=1,2,3
(3)
式中Tp,Dj,cpq(q=1,2,…,6),ejp和εij(1,2,3)分別為各方向應力、電位移、石英晶體材料的彈性常數、壓電常數和介電常數[8]。
考慮旋轉場作用下半無限大石英晶體的運動方程為[12]
(4)
將應力和電位移表達式(3)代入運動方程式(4),經整理后得到關于A1(A2,A3,A4)的方程組。如這些振幅存在非零解,必須要求振幅A1(A2,A3,A4)的系數行列式的值為零。這樣可以獲得關于衰減系數β的4對解,取其沿著深度方向呈衰減的衰減系數[14,15],重寫位移和電勢如下
(5)
式中A4n(n=1,2,3,4)和βn為新定義的振幅和獲得的衰減系數。
根據新位移和電勢式(5),基于式(3),重寫應力和電位移表達式。當聲表面波在如圖1所示的半無限大石英晶體內沿著x1方向傳播時,必須滿足力學邊界條件[14,15]。同時對于石英晶體這一壓電彈性體,除了滿足力學邊界條件外,還應該滿足電學邊界條件[14,15]。全部邊界條件如下
T2=T4=T6=D2=0,x2=0
(6)
將新的應力和電位移表達式代入邊界條件表達式(6),可以得到振幅A4n的四個方程組。若聲表面波存在非零解,必然要求其系數行列式的值為零,這樣就建立了聲表面波的波速方程??梢苑謩e繪制出繞x1軸和x2軸的旋轉頻率和聲表面波波速的關系曲線如圖2所示。這里選用的石英晶體切型為傳感器領域常用的ST切型,其材料參數見文獻[8]。
圖2 半無限大石英晶體的波速—旋轉頻率曲線
在實際工作環(huán)境中,并不能確定旋轉場的旋轉軸與聲表面波傳播方向的角度??梢赃M一步繪制在一定旋轉頻率下,即Ω/c0=500時,聲表面波波速和旋轉軸角度之間的關系如圖3所示。
圖3 半無限大石英晶體的波速—旋轉軸角度的關系曲線
從圖3可以發(fā)現,隨著旋轉軸角度的變化,聲表面波波速也發(fā)生變化。在實際應用中,不僅需要檢測出旋轉場的旋轉頻率,也需要檢測出旋轉軸角度。結合圖2和圖3,可以進行聯立求解,確定旋轉場的旋轉頻率和旋轉軸角度[14,15]。
本研究初步納入了2015年1月至2017年1月期間就診于上海市浦東醫(yī)院眼科的患者82例(100眼),于2015年1月獲得上海市浦東醫(yī)院倫理委員會批準(批號:LW2015-001),遵循赫爾辛基宣言。所有患者均告知檢查目的及檢查內容,手術前均簽署手術知情同意書。
當聲表面波在圖1所示的無限大石英晶體板內傳播時,其位移將會呈現出板波的形式[16]。因此,運動方程式(4)獲得的衰減系數將全部應用到新的位移和電勢假設。按照式(5)和式(3)可以寫出新的位移、電勢、應力和電位移表達式。
對于如圖1所示的半無限大石英晶體板的上下表面,其應力和電位移邊界條件如下
T2=T4=T6=D2=0,x2=0,
T2=T4=T6=D2=0,x2=-h
(7)
式中h為石英晶體板的板厚。為了計算簡單,對石英晶體板厚度進行如下歸一化處理
H=h/ξ,kξ=2π,kh=2πH
(8)
式中H和ξ分別為石英晶體板的歸一化厚度和聲表面波波長,這里無限大板的材料仍然取ST切石英晶體作為研究對象。
從圖4~圖5可以發(fā)現,聲表面波在無限大石英晶體板內傳播時,出現了兩個聲表面波波速,并且隨著石英晶體板厚度的增加,波速的上下分支逐漸接近。這一計算結果符合過去的研究和實驗結果[14,15]。同樣隨著旋轉頻率的增加,石英晶體板內各聲表面波波速都在減小。與半無限大石英晶體的結果類似,繞軸旋轉時的波速減小趨勢更快,并且相應的線性關系可以應用于旋轉場的傳感[10~12]。這里建立的模型更加符合實際傳感器的層狀板型結構。
圖4 半無限大石英晶體板的波速—旋轉頻率曲線(三個波長)
圖5 半無限大石英晶體板的波速—旋轉頻率曲線(四個波長)
通過對考慮陀螺效應的聲表面波傳播方程的求解,獲得了聲表面波波速、旋轉頻率和旋轉軸角度之間的關系曲線圖。研究結果表明:旋轉場對于聲表面波的傳播都有著顯著影響,特別是繞聲表面波傳播方向旋轉時影響效果更為明顯。當聲表面波在無限大石英晶體板傳播時,隨著板厚的增加,聲表面波波速逐漸趨向于無限大的情況。這里所建的模型為考慮壓電效應的層狀結構,更符合實際的聲表面波壓電器件分析模型。