基于COMSOL的聲表面波器件仿真*

阮 鵬,陳智軍,2,付大豐,張亦居,王萌陽(yáng)

(1.南京航空航天大學(xué)自動(dòng)化學(xué)院,江蘇南京 210016; 2.南京大學(xué)近代聲學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210093)

聲表面波(Surface Acoustic Wave,SAW)是一種沿固體表面?zhèn)鞑サ膹椥圆?能量集中在固體表面1～ 2個(gè)波長(zhǎng)的范圍內(nèi)[1],可利用沉積在壓電基片表面的叉指換能器(Inter-digital Transducer,IDT)來(lái)有效激發(fā)[2].SAW器件在通信、傳感、射頻識(shí)別等領(lǐng)域都有著極為廣泛的應(yīng)用,最常見(jiàn)的結(jié)構(gòu)為“半無(wú)限壓電基片-叉指換能器”形式[3-4].隨著研究的不斷深入,出現(xiàn)了各種改進(jìn)的SAW器件結(jié)構(gòu),譬如在叉指換能器表面濺射一層薄膜的樂(lè)甫波(Love Wave)器件[5].由于適合液體檢測(cè),樂(lè)甫波器件成為聲表面波器件的研究熱點(diǎn)[6-7].

在對(duì)聲表面波器件的研究過(guò)程中,器件仿真工作是器件設(shè)計(jì)和制作的基礎(chǔ),通過(guò)仿真可以提高設(shè)計(jì)效率、減少制造過(guò)程中的材料浪費(fèi)等.聲表面波器件的仿真可基于脈沖響應(yīng)模型、等效電路模型、耦合模型等方法,通過(guò)編寫(xiě)代碼進(jìn)行數(shù)值分析[8].隨著ANSYS、ANSOFT、COMSOL等大型有限元分析軟件的出現(xiàn),可利用上述軟件對(duì)SAW器件進(jìn)行仿真分析,從而省去人工編寫(xiě)代碼的繁瑣流程[9].

與ANSYS、ANSOFT等有限元仿真軟件相比,COMSOL具有以下顯著特點(diǎn):①完全開(kāi)放的架構(gòu),用戶可以自由定義所需專業(yè)的偏微分方程;②專業(yè)的計(jì)算模型庫(kù),內(nèi)置各種常用的物理模型,用戶可以選擇參考并進(jìn)行必要的修改;③豐富的后處理功能,可根據(jù)用戶需要進(jìn)行各種數(shù)據(jù)、曲線、圖片及動(dòng)畫(huà)的輸出與分析[10].

本文基于COMSOL建立了SAW器件的有限元分析模型,得到了直觀的SAW傳播形式;仿真分析了SAW器件的頻率特性,探討了IDT在不同敷金比、高度比時(shí)的SAW器件頻率變化.在此基礎(chǔ)上,研究了在IDT表面濺射一層薄膜時(shí)的樂(lè)甫波器件頻率特性;當(dāng)薄膜上負(fù)載液體時(shí),通過(guò)器件頻率對(duì)液體密度的靈敏程度來(lái)對(duì)薄膜厚度進(jìn)行了優(yōu)化.其研究結(jié)果可以為聲表面波器件的精確設(shè)計(jì)提供仿真依據(jù).

1 模型的建立

聲表面波器件是一種典型的諧振式器件,特征頻率是器件的主要參數(shù)[11-12].當(dāng)聲表面波器件用做傳感器時(shí),主要通過(guò)特征頻率隨待測(cè)量的變化來(lái)實(shí)現(xiàn)檢測(cè)功能.在有限元分析中,模態(tài)是機(jī)械結(jié)構(gòu)的固有振動(dòng)特性,每一個(gè)模態(tài)具有特定的特征頻率和模態(tài)振型.因此,通過(guò)有限元軟件COMSOL的模態(tài)分析模塊可仿真研究聲表面波器件的頻率特性.

圖1 聲表面波器件結(jié)構(gòu)Fig.1 The structure of surface acoustic w ave device

圖2 模態(tài)分析二維模型Fig.2 2D model of analysis

聲表面波器件的結(jié)構(gòu)如圖1所示.IDT沉積在壓電基片表面,由周期性排列并與匯流條交替連接的多對(duì)電極構(gòu)成.利用周期性邊界條件,可把IDT簡(jiǎn)化為由一對(duì)電極組成的周期結(jié)構(gòu),從而建立簡(jiǎn)化的二維模型以減小計(jì)算量.聲表面波器件的模態(tài)分析二維模型如圖2所示.橫、縱坐標(biāo)的單位:μm,基片寬度與聲表面波波長(zhǎng) λ一致;由于聲表面波的能量主要集中在基片表面1～2個(gè)波長(zhǎng)深度內(nèi),隨著基片深度的增加其振幅呈指數(shù)衰減,因此基片高度可簡(jiǎn)化為2～5倍波長(zhǎng);電極中心間距p=λ/2.

模型邊界條件的設(shè)定如下:

1)上表面邊界條件 ΓT為自由邊界條件;

2)下表面邊界條件 ΓB為固定邊界條件;

3)左右邊界條件 ΓT、ΓR為一對(duì)周期性邊界條件,即左右邊界上對(duì)應(yīng)節(jié)點(diǎn)的3個(gè)自由度(x方向位移u,y方向位移v,電勢(shì) Φ)相等.

二維模型各參數(shù)如表1所示.

表1 模態(tài)分析二維模型參數(shù)Tab.1 Param eter of 2D analysismodel

基片材料為壓電陶瓷PZT-5H,其密度ρ為7 500 kg/m3,壓電、彈性、介電常數(shù)分別為

IDT采用金屬鋁電極,其密度 ρIDT=2 700,kg/m3,楊氏模量E=70e9, Pa,泊松比v=0.33.

建立完備的幾何模型還需要進(jìn)行網(wǎng)格劃分.網(wǎng)格劃分越密,計(jì)算結(jié)果越精確,但是計(jì)算量會(huì)大大增加,因此合適的網(wǎng)格劃分是求解的關(guān)鍵.COMSOL軟件的網(wǎng)格劃分可以分為初始化網(wǎng)格、細(xì)化網(wǎng)格和自由網(wǎng)格三種方法.初始化網(wǎng)格由軟件根據(jù)幾何模型粗略進(jìn)行自動(dòng)網(wǎng)格劃分,根據(jù)模型各處結(jié)構(gòu)的復(fù)雜程度決定網(wǎng)格尺寸;細(xì)化網(wǎng)格在初始化網(wǎng)格的基礎(chǔ)上由軟件進(jìn)行更密集的劃分;自由網(wǎng)格對(duì)用戶開(kāi)放,用戶可以設(shè)定最大單元尺寸、單元增長(zhǎng)率等參數(shù),并可以決定局部細(xì)化的區(qū)域.

網(wǎng)格劃分后的二維模型如圖 3所示.首先對(duì)整體模型結(jié)構(gòu)進(jìn)行初始化網(wǎng)格,但由于聲表面波集中在基片表面1～2個(gè)波長(zhǎng)內(nèi)傳播,因此再對(duì)基片表面進(jìn)行局部細(xì)化.

圖3 模型的網(wǎng)格劃分Fig.3 Meshing of model

2 無(wú)電極壓電基片結(jié)構(gòu)

在實(shí)際應(yīng)用中,通常通過(guò)IDT來(lái)有效激發(fā)出聲表面波.如果不考慮聲表面波激發(fā)原理,僅考慮其傳播時(shí)需要滿足的邊界條件,同樣可以通過(guò)COMSOL的模態(tài)分析模塊來(lái)研究無(wú)電極壓電基片結(jié)構(gòu)的聲表面波特征頻率,此時(shí)只需在上述模型中去掉叉指電極即可.

在模型建立完成后,選擇特征頻率處理器,確定求解范圍來(lái)進(jìn)行聲表面波特征頻率的分析.利用聲表面波能量集中在基片表面1～2個(gè)波長(zhǎng)深度內(nèi)的特點(diǎn),可以把聲表面波振型從很多各不相同的波動(dòng)振型中提取出來(lái),得到特征頻率.

聲波振型如圖4所示,圖中右邊的顏色刻度表示振動(dòng)位移大小,顏色從上到下表示振動(dòng)位移減弱.根據(jù)振動(dòng)位移隨基片深度的變化趨勢(shì),圖4(a)為體波振型,圖4(b)為聲表面波振型,其對(duì)應(yīng)的聲表面波特征頻率f為189.78MH z.通過(guò)多次仿真表明,每一個(gè)無(wú)電極壓電基片模型對(duì)應(yīng)唯一的聲表面波振型.

圖4 聲波振型Fig.4 Acoustic w ave vibration model

3 壓電基片-叉指換能器結(jié)構(gòu)

3.1 叉指電極敷金比分析

當(dāng)壓電基片上沉積叉指電極時(shí),以表 1所述的模型參數(shù)進(jìn)行仿真.與無(wú)電極壓電基片結(jié)構(gòu)不同的是,通過(guò)模態(tài)分析模塊可以提取到兩個(gè)符合聲表面波振型的特征頻率,定義為正特征頻率fsc+和反特征頻率fsc-分別如圖5(a)、(b)所示.

表2 聲表面波特征頻率隨MR的變化Tab.2 Resonant frequencies of Surface acoustic w ave's change with MR

正特征頻率和反特征頻率的出現(xiàn)是由于 IDT的電極效應(yīng)引起[13].定義叉指電極敷金比MR=a/P,當(dāng)=0時(shí),可以等效為基片上表面無(wú)電極情況,為自由化邊界條件,此時(shí)表面電荷密度 σ =0;當(dāng)MR=1時(shí),可以等效為基片上表面覆蓋一層金屬薄膜,為金屬化邊界條件,此時(shí)表面電勢(shì)Φ=0.通過(guò)仿真分析,得到正、反特征頻率隨的變化如表2所示.

當(dāng)0＜MR＜1,隨著MR的增大,正、反特征頻率逐漸減小.當(dāng)MR=0與MR=1時(shí),基片上表面分別處于自由化邊界條件和金屬化邊界條件,不存在IDT的電極效應(yīng),正特征頻率和反特征頻率退化為一個(gè)特征頻率.

3.2 叉指電極高度分析

在敷金比分析的基礎(chǔ)上,取=0.5時(shí)進(jìn)行IDT高度的分析,即在IDT寬度不變的前提下,研究IDT電極高度對(duì)拓征頻率的影響.定義叉指電極高度比

特征頻率隨叉指電極高度比Q的變化如表3所示.

圖5 正、反特征頻率振型Fig.5 The vibration m odel of resonan t and an ti-resonan t frequency

表3 特征頻率隨叉指電極高度比 Q的變化Tab.3 Resonant frequen cies'change w ith the heigh t of IDT

由上述分析可知,Q=0時(shí),正、反特征頻率相等,其原因在于Q=0與=0等效,相當(dāng)于無(wú)電極壓電基片結(jié)構(gòu).其它情況下叉指電極均激發(fā)出正、反特征頻率,并且當(dāng)Q值增大時(shí),正、返特征頻率會(huì)向低頻偏移.

4 壓電基片-叉指換能器-薄膜結(jié)構(gòu)

在“壓電基片-IDT”結(jié)構(gòu)的基礎(chǔ)上,在IDT表面濺射一層薄膜時(shí),該結(jié)構(gòu)稱為樂(lè)甫波器件.樂(lè)甫波器件是一種新型的聲表面波器件,因適于液體檢測(cè)而成為研究熱點(diǎn).

圖6 樂(lè)甫波三維振型Fig.6 3D vibration m odel of love w ave

圖7 樂(lè)甫波特征頻率 f隨薄膜厚度d的變化Fig.7 Resonant frequencies of love w ave's change w ith the thickness of th in film

器件模型在圖2的基礎(chǔ)上增加一層薄膜,薄膜材料選擇二氧化鈦(TiO2),薄膜上表面為自由邊界條件,左右邊界條件為一對(duì)周期性邊界條件.樂(lè)甫波三維振型如圖6所示.從圖6可以看出,振動(dòng)能量主要還是集中在壓電基片表面一個(gè)波長(zhǎng)深度左右,并且薄膜上的振動(dòng)也很明顯,這可以為液體的檢測(cè)提供很好的保證.

樂(lè)甫波特征頻率隨薄膜厚度的變化,如圖7所示.從圖7可知,與IDT高度的影響相似,由于薄膜的質(zhì)量加載效應(yīng),特征頻率同樣隨薄膜厚度的增大有向低頻偏移的趨勢(shì).

5 壓電基片-叉指換能器-薄膜-液體結(jié)構(gòu)

基于“壓電基片-IDT-薄膜”結(jié)構(gòu),在薄膜上負(fù)載液體時(shí),樂(lè)甫波三維振型如圖8所示.由圖8可知,振動(dòng)位移不會(huì)產(chǎn)生壓縮波進(jìn)入液體,從而不會(huì)向液體介質(zhì)輻射能量,因此樂(lè)甫波在傳播過(guò)程中不會(huì)產(chǎn)生太大衰減,表明了樂(lè)甫波用于液體檢測(cè)的優(yōu)勢(shì).

當(dāng)樂(lè)甫波器件用于液體密度檢測(cè)時(shí),定義靈敏度

圖8 負(fù)載液體的樂(lè)甫波三維振型Fig.8 Love wave 3D vibration m odel of load ing liquid

式中:ρ1和 ρ2為液體的密度;fρ1和fρ2分別為負(fù)載上述密度液體時(shí)對(duì)應(yīng)的樂(lè)甫波特征頻率.

靈敏度Sρ隨薄膜厚度的變化,如圖9所示.從圖9可知,靈敏度Sρ首先隨薄膜厚度d的增大而增大;當(dāng)d=5μm,即d/λ≈0.45時(shí),靈敏度Sρ達(dá)到最大;之后隨d的增大而減小,并最終趨于穩(wěn)定.

選擇薄膜d最優(yōu)厚度為5μm時(shí),樂(lè)甫波特征頻率f隨液體密度的變化,如圖10所示.從圖10可知,樂(lè)甫波器件用于液體密度檢測(cè)時(shí)具有較好的線性度.

圖9 靈敏度Sρ隨薄膜厚度d的變化Fig.9 Sensitivities'change w ith the thickness of thin film

圖10 樂(lè)甫波特征頻率 f隨液體密度的變化Fig.10 Resonant frequencies of love w ave's change w ith the density of liquid

6 結(jié) 論

利用有限元軟件COMSOL,按由淺入深的順序?qū)o(wú)電極壓電基片、壓電基片表面沉積叉指換能器、叉指換能器表面濺射薄膜、薄膜上負(fù)載液體的4種聲表面波器件結(jié)構(gòu)進(jìn)行了仿真分析,可得出以下結(jié)論:

1)通過(guò)有限元軟件COMSOL的模態(tài)分析模塊,利用聲表面波能量集中在基片表面1～2個(gè)波長(zhǎng)深度內(nèi)的特點(diǎn),可以把聲表面波振型從很多各不相同的波動(dòng)振型中提取出來(lái),并得到特征頻率.

2)對(duì)于無(wú)電極壓電基片結(jié)構(gòu),每一個(gè)模型對(duì)應(yīng)唯一的聲表面波振型.

3)叉指換能器的電極效應(yīng)會(huì)產(chǎn)生正、反特征頻率,而且正、反特征頻率隨著叉指換能器敷金比與高度比的增加而向低頻偏移.

4)對(duì)于“壓電基片-叉指換能器-薄膜”結(jié)構(gòu),樂(lè)甫波振動(dòng)能量主要還是集中在壓電基片表面大約一個(gè)波長(zhǎng)深度,并且薄膜上的振動(dòng)也很明顯.與此同時(shí),由于薄膜的質(zhì)量加載效應(yīng),特征頻率隨薄膜厚度的增大有向低頻偏移的趨勢(shì).

5)在薄膜上負(fù)載液體時(shí),振動(dòng)位移不會(huì)產(chǎn)生壓縮波進(jìn)入液體,從而不會(huì)向液體介質(zhì)輻射能量,因此樂(lè)甫波在傳播過(guò)程中不會(huì)產(chǎn)生太大衰減,表明了樂(lè)甫波用于液體檢測(cè)的優(yōu)勢(shì).

6)當(dāng)樂(lè)甫波器件用于液體密度檢測(cè)時(shí),可通過(guò)器件特征頻率對(duì)液體密度的靈敏程度對(duì)薄膜厚度進(jìn)行優(yōu)化.仿真結(jié)果表明,樂(lè)甫波器件用于液體密度檢測(cè)時(shí)具有較好的線性度.

[1] Morgan D P.History of SAWdevices[C].In Proceedings of IEEE International frequency control symposium,Piscataw ay, NJ,USA,1998:439-460.

[2] 水永安.聲表面波與聲表面波器件講義[R].南京:南京大學(xué)電子科學(xué)與技術(shù)系,1998.

[3] Hashimoto K Y.Surface AcousticWave Devicesin Telecommunications:Modeling and Simulation[M].Berlin,Heidelberg: Springer-Verlag,2000.

[4] 陳明,范東遠(yuǎn),李歲勞.聲表面波傳感器[M].西安:西北工業(yè)大學(xué)出版社,1997.

[5] Tamarin O.Study of acoustic Love wave devices for real timebacteriophage detection[J].Sensors and Actuators,2003,B (91):275-284.

[6] 胡友旺,史建亮,賈宏光,等.Love波傳感器及其質(zhì)量靈敏度測(cè)試[J].傳感技術(shù)學(xué)報(bào),2006,19(5):2126-2129.

Hu Youwang,Shi Jianliang,Jia Hongguang,etal.Lovew ave sensor and themeasurementofmass sensitivity[J].Chinese Journal o f Sensors and A ctuators,2006,19(5):2126-2129.(in Chinese)

[7] 江城,倪世宏,李紅浪,等.Love波飛機(jī)結(jié)冰傳感器研究[J].航空學(xué)報(bào),2009,30(5):906-911.

Jiang Cheng,NiShihong,LiHonglang,etal.Study of lovew aveaircraft ice sensor[J].Acta Aeronautica et Astronautica Sinica,2009,30(5):906-911.(in Chinese)

[8] 吳佚卓.高頻無(wú)線無(wú)源聲表面波傳感器的仿真與測(cè)試研究[D].天津:天津理工大學(xué),2006.

[9] Peng Dasong,Yu Fengqi.A novelFEA simu lationmodel for RFIDSAWtag[J].IEEE Transactions on Ultrasonics,Ferroelectrics,and Frequency Control,2009,56(8):1753-1760.

[10] K rishnan R N,Nemade H B,Paily R.Simulation ofOne-PortSAWResonator using COMSOLMultiphysics[C].Boston, Proceedings of the COMSOL Users Conference,2006.

[11] 戴恩光,馮冠平.對(duì)SAW諧振器氣體傳感器的表面質(zhì)量負(fù)載與傳感器性能之間關(guān)系的研究[J].測(cè)試技術(shù)學(xué)報(bào), 1996,10(4):349-353.

Dai Enguang,Feng Guanping.The research in the relationships betw een the loading effects and the performance of sensor based upon SAWresonator[J].Journalof Test and Measurement Technology,1996,10(4):349-353.(in Chinese)

[12] Drafts B.A coustic wave technology sensors[J].IEEE Transactions on Microwave Theory and Techniques,2001,49(4): 795-802.

[13] 盧旭.基于COMSOL的聲表面波標(biāo)簽仿真[D].南京:南京航空航天大學(xué),2011.