一維聲子晶體帶隙研究概述

邱學(xué)云

(文山學(xué)院數(shù)理系,云南文山663000)

一維聲子晶體帶隙研究概述

邱學(xué)云

(文山學(xué)院數(shù)理系,云南文山663000)

文章對聲子晶體的概念、基本特征和分類進行簡要概述.總結(jié)分析一維聲子晶體的三種簡化模型和一維聲子晶體的計算方法及其帶隙研究成果.展望一維聲子晶體在低頻振動、噪聲控制、抗振防震方面的應(yīng)用可能,為深入研究一維聲子晶體提供依據(jù).

一維聲子晶體;能帶結(jié)構(gòu);帶隙

半導(dǎo)體中的電子與周期分布的原子勢場相互作用,使得半導(dǎo)體形成電子帶隙,即電子能帶間的頻率范圍,也叫電子禁帶,人為設(shè)計能夠控制電子的流動.近年來,能帶理論突破了以固有材料為研究對象的限制,進入了通過能帶設(shè)計來模擬實際晶格以獲得新型功能材料和器件的新階段.在這些材料中存在能夠禁止某種經(jīng)典波傳播的頻率范圍,這些頻率范圍稱為帶隙.具有經(jīng)典波帶隙的周期性復(fù)合材料或結(jié)構(gòu)統(tǒng)稱波晶體.通常把存在電磁波帶隙,介電常數(shù)周期分布的材料或結(jié)構(gòu)稱光子晶體,把存在彈性波帶隙,彈性常數(shù)及密度周期分布的材料或結(jié)構(gòu)稱聲子晶體.

文章對聲子晶體的概念、基本特征、和分類進行簡要概述,總結(jié)分析一維聲子晶體的計算方法及其帶隙研究成果,展望一維聲子晶體的應(yīng)用可能,為深入一維聲子晶體的應(yīng)用研究提供依據(jù).

1 聲子晶體概述

1.1 聲子晶體概念和基本特征

1987 人們年發(fā)現(xiàn)了一種新型光子材料--光子晶體[1-2],在光子晶體研究的基礎(chǔ)上,提出了聲子晶體的概念.M.S.Kushwaha等1993年第1次提出聲子晶體概念時,就對電子晶體、光子晶體和聲子晶體的有關(guān)特性進行了比較[3],通過比較研究,得出它們非常相似性[4].聲子晶體是一種具有彈性波帶隙的周期性結(jié)構(gòu)材料,即在帶隙頻率范圍內(nèi)的彈性波不能夠傳播.基于這一特性,聲子晶體在隔聲、降噪、減振方面有廣闊的應(yīng)用可能.

聲子晶體內(nèi)部材料的彈性常數(shù)、密度、楊氏模量等材料參數(shù)周期性變化和材料結(jié)構(gòu)參數(shù)晶格尺寸、組份比的不同,聲子晶體的彈性波帶隙也就不同.帶隙分完全帶隙和不完全帶隙,在特定的頻率范圍內(nèi),波在波矢的所有方向上都不能傳播的稱為完全帶隙.頻率范圍內(nèi)只允許某些方向上的波通過,其它方向禁止通過的帶隙具有方向性,稱為不完全帶隙.一般來說,聲子晶體中各組份比越大,入射波將被散射得越強烈,就越容易產(chǎn)生帶隙[5].

1.2 聲子晶體的分類

聲子晶體一般由兩相或以上的彈性介質(zhì)復(fù)合組成.根據(jù)聲子晶體結(jié)構(gòu)在笛卡爾坐標(biāo)系中周期排列形式的不同,分為一維、二維及三維結(jié)構(gòu).如圖1所示,從左到右分別為一維(層狀)、二維(柱狀)、三維(正方體)聲子晶體結(jié)構(gòu)示意圖.聲子晶體結(jié)構(gòu)常用的組份形態(tài)有層狀、柱狀、長方體、正方體、球形、橢球形等,材料可以是實心或中空的,可以是不同物相不同組份的復(fù)合.

圖1 聲子晶體結(jié)構(gòu)示意圖

2 一維聲子晶體簡化模型

一維聲子晶體一般有層狀結(jié)構(gòu)和桿狀結(jié)構(gòu)兩種.對一維聲子晶體帶隙進計算分析時,通過調(diào)控組元材料的結(jié)構(gòu)參數(shù)或材料參數(shù)來調(diào)控聲子晶體的帶隙情況.對于一維聲子晶體,當(dāng)相鄰兩個離散單元為同種材料時沿x方向的拉壓剛度為2ES/(dj+1+dj),不同種材料時拉壓剛度為2EAEBS/(EAdj+1+EBdj),E為彈性模量.對于原胞中包含兩種以上材料疊合的情形也可以按上述方法進行簡化處理.[6]

2.1 嚴(yán)格周期性聲子晶體模型

如下圖2所示是一維嚴(yán)格周期性聲子晶體簡化模型.圖(a)是一維二組元結(jié)構(gòu),圖(b)是一維三組元結(jié)構(gòu).不同彈性常數(shù)和密度的材料A、B或材料A、B、C沿x方向交替排列形成一維二組元(三組元)聲子晶體的簡化模型.用集中質(zhì)量法可以將其分解為有限個集中質(zhì)量,各個集中質(zhì)量間之間的連接簡化為無質(zhì)量的彈簧連接,聲子晶體原胞由材料A和材料B組成,理論上是無限多自由度的連續(xù)介質(zhì)系統(tǒng).因此,其原胞可簡化為有限個自由度的彈簧振子結(jié)構(gòu),而一維聲子晶體則簡化為周期彈簧振子結(jié)構(gòu).[7]如下圖3所示是一維二組元聲子晶體的離散示意圖,我們只要找到離散后的集中質(zhì)量和彈簧剛度與連續(xù)介質(zhì)材料參數(shù)之間的關(guān)系,就可以計算連續(xù)介質(zhì)的一維二組元聲子晶體的帶隙結(jié)構(gòu).

2.2 層厚遞變式一維準(zhǔn)周期結(jié)構(gòu)聲子晶體模型

曹永軍等[8]提出一種層厚遞變式一維準(zhǔn)周期結(jié)構(gòu)聲子晶體模型.如圖4所示,系統(tǒng)共有N個周期2N層,分3種情況:第一種情況是第1個周期中兩種介質(zhì)的厚度為dA和dB,其后每個周期的厚度同時依次遞增或遞減一個Δd;第二種情況是第1個周期中兩種介質(zhì)的厚度為dA和dB,其后每個周期中介質(zhì)A的厚度不變,而介質(zhì)B的厚度依次遞增或遞減一個Δd;第三種情況是第1個周期中兩種介質(zhì)的厚度為dA和dB,其后每個周期中介質(zhì)B的厚度不變,而介質(zhì)A的厚度依次遞增或遞減一個Δd.通過計算彈性波通過該一維準(zhǔn)周期結(jié)構(gòu)聲子晶體的透射系數(shù),并與周期結(jié)構(gòu)聲子晶體的透射系數(shù)進行比較.研究發(fā)現(xiàn)利用準(zhǔn)周期排列結(jié)構(gòu)可以有效地調(diào)節(jié)聲子晶體的帶隙寬度和所在的頻率范圍.

圖4 層厚遞變式準(zhǔn)周期結(jié)構(gòu)聲子晶體模型

2.3 一維功能梯度材料聲子晶體模型

宿星亮等[9]提出了一維功能梯度材料聲子晶體模型.如圖5所示,該模型研究的一維聲子晶體是由材料常數(shù)按指數(shù)形式分布的功能梯度材料沿x軸方向周期排列而成.該功能梯度材料單元兩端表面材料常數(shù)相同,在原胞中心a處達到材料常數(shù)峰值,周期排列后構(gòu)成材料常數(shù)宏觀上連續(xù)變化的一維結(jié)構(gòu).

圖5 一維功能梯度材料聲子晶體模型

上述三種一維聲子晶體模型:第一種是材料參數(shù)和結(jié)構(gòu)參數(shù)嚴(yán)格周期性變化的;第二種是材料參數(shù)不變,結(jié)構(gòu)參數(shù)層厚遞變式的準(zhǔn)周期性結(jié)構(gòu);第三種是材料參數(shù)按指數(shù)形式分布變化,結(jié)構(gòu)參數(shù)不變的準(zhǔn)周期性結(jié)構(gòu).三種模型各有優(yōu)點,都具有各自的實際應(yīng)用價值.

3 一維聲子晶體帶隙計算及其研究

國內(nèi)外有關(guān)聲子晶體帶隙研究的文獻主要是對聲子晶體的布拉格散射機理[10-12]、局域共振帶隙機理[13-18]和帶隙特性進行研究.局域共振帶隙機理由我國劉正猷教授等人2000年提出.局域共振帶隙機理認為,在特定的彈性波激勵下,聲子晶體結(jié)構(gòu)基體中的散射單元產(chǎn)生共振,并同彈性波相互作用,從而抑制其傳播.這為低頻振動與噪聲控制的應(yīng)用研究開辟了一條新思路[19].另外,王剛等人主要研究了不同周期結(jié)構(gòu)的帶隙計算方法和帶隙特性,研究也取得了很多的成果[20].下面是常用的一維聲子晶體帶隙的計算方法及其部分研究成果.

3.1 傳遞矩陣法

傳遞矩陣法[21-22]是從連續(xù)狀態(tài)參數(shù)應(yīng)力、質(zhì)點位移等的基本方程入手,結(jié)合界面連續(xù)條件,得出單個周期的傳遞矩陣.通過引入周期邊界條件得到相應(yīng)的色散關(guān)系即能帶結(jié)構(gòu),同時通過有限個傳遞矩陣相乘可得到有限周期傳輸特性.該方法計算量較小,可計算一維聲子晶體帶隙,不能直接處理二維和三維聲子晶體的帶隙.

郁殿龍等[23]利用該法研究了表面局域態(tài)對一維聲子晶體中水平剪切波傳輸特性的影響情況,研究表明共振峰的極值與入射角度和聲子晶體層數(shù)有關(guān),合適的入射角度和層數(shù)可以使聲波完全透射.當(dāng)入射角在一定范圍內(nèi)連續(xù)變化時,在較寬頻率范圍內(nèi)均出現(xiàn)較大透過率.聲子晶體的這一特性可以應(yīng)用于高性能的阻抗匹配材料和聲波濾波器中.

蔣澤等[24]應(yīng)用廣義傳輸矩陣法,建立了聲波傳播特性的理論分析模型,得到了其聲波場的平面波解,給出了數(shù)值實現(xiàn)方案.其研究表明,該方法可精確地模擬彈性波通過一維有限厚的嚴(yán)格周期結(jié)構(gòu)、準(zhǔn)周期結(jié)構(gòu)以及完全無序結(jié)構(gòu)的傳播特性.

3.2 平面波展開法

對于平面波展開法[25-26]是因為聲子晶體具有周期性,可將相關(guān)參數(shù)按傅里葉級數(shù)展開,結(jié)合Bloch定理,把聲子晶體波動方程放到倒格矢空間以平面波疊加的形式展開,將聲子晶體波動方程的求解轉(zhuǎn)化成本征值問題,從而得到頻率與波矢之間的色散關(guān)系即聲子晶體的帶隙結(jié)構(gòu).該法可用于計算一維、二維、三維聲子晶體中固體/固體、液體(氣體)/液體(氣體)的復(fù)合結(jié)構(gòu),但在計算、液體(氣體)/固體結(jié)構(gòu)時存在困難.應(yīng)用該法時當(dāng)組元材料參數(shù)差異較大時計算量大,收斂慢,但是隨著計算機的更新?lián)Q代,這個問題已經(jīng)得到改善.宿星亮,等[9]應(yīng)用此法研究了由功能梯度材料周期復(fù)合而成的一維聲子晶體中存在的彈性波帶隙特征,結(jié)果表明功能梯度材料聲子晶體較常規(guī)材料聲子晶體在相同范圍內(nèi)能夠出現(xiàn)更多階帶隙結(jié)構(gòu).

肖偉等[27]提出用波傳播法來研究一維聲子晶體的帶隙特性,將該方法與平面波展開法進行比較,發(fā)現(xiàn)平面波展開法隨著波數(shù)的增加而逐漸收斂于波傳播法的結(jié)果.如將波傳播法應(yīng)用于一維二組元和一維四組元聲子晶體禁帶特性的計算中,在相同計算精度下波傳播法的計算時間大約為平面波展開法的1/50和1/100.當(dāng)考慮到粘彈性材料的頻變特性時,波傳播法能直接得到聲子晶體的禁帶特性,在相同的計算精度下波傳播法的計算量大約要比經(jīng)過迭代改進的平面波展開法的計算量小兩個數(shù)量級.

閆志忠等[28]發(fā)展了一種基于小波的一維聲子晶體彈性波帶隙計算方法,將彈性波場在小波基上展開,得到一個關(guān)于自適應(yīng)計算小波積分的一般矩陣特征值問題.將該方法應(yīng)用到二元體系的聲子晶體,與傳統(tǒng)平面波展開法相比,該方法的計算結(jié)果與之相符合,而且可在得到同樣計算精度的條件下,顯著降低計算量,提高計算速度.

王剛等[29]采用迭代法改進了一維聲子晶體帶隙特性計算的平面波展開算法,以使其適用于組成材料粘彈性所導(dǎo)致的彈性常數(shù)隨頻率非線性變化的特性.在將該算法應(yīng)用于丁腈橡膠和鋼組成的一維周期結(jié)構(gòu)聲子晶體振動帶隙的研究中,理論計算和振動測試結(jié)果吻合理想.

3.3 集中質(zhì)量法

集中質(zhì)量法[30-31]是基于振動力學(xué)中連續(xù)系統(tǒng)的離散化思想,在聲子晶體中將各組元連續(xù)介質(zhì)中的質(zhì)量集中到有限個節(jié)點或截面上,把有限個節(jié)點或截面視為有限多個自由度的彈簧振子結(jié)構(gòu),即將聲子晶體彈性波帶隙的計算簡化為計算周期彈簧振子結(jié)構(gòu)的彈性波帶隙.其本質(zhì)是將無限自由度系統(tǒng)轉(zhuǎn)化成有限自由度系統(tǒng)近似求解.這種方法特別適合計算大彈性常數(shù)差組份復(fù)合而成的一維聲子晶體,并且這種方法可以更加直觀地描述聲子晶體內(nèi)部作用機理,這對聲子晶體帶隙的產(chǎn)生機理揭示將起到重要作用.

劉鐵權(quán)等[32]將一維聲子晶體的原胞簡化為有限多個自由度的彈簧振子結(jié)構(gòu)后,在辛對偶變量體系下探討晶格振動,引入辛數(shù)學(xué)方法確定波矢與本證值的色散關(guān)系.通過本證值計數(shù)法計算特征頻率,從而得到禁帶區(qū)間.研究認為與傳統(tǒng)集中質(zhì)量法相比,該算法的計算結(jié)果與之吻合很好,且提高了計算精度和計算效率,在低頻處收斂性更好,可以借鑒參考.

結(jié)束語

深入研究一維聲子晶體的帶隙特性,將會發(fā)現(xiàn)許多新的物理現(xiàn)象,從理論計算尋找到具有帶隙起始頻率低且有一定帶寬的周期結(jié)構(gòu),分析其在噪聲控制、低頻防震方面的應(yīng)用可能,理論設(shè)計出具有隔聲、降噪、減振性能的一維聲子晶體模型和具有低頻防震性能的工程模型,可為工程應(yīng)用做貢獻.

作為一種新材料,各種結(jié)構(gòu)聲子晶體的應(yīng)用都還處于展望階段,但由于聲子晶體所具有的特殊性質(zhì)使得其在航空航天、電子器件、人造臟器、汽車發(fā)動機、聲功能器件、等諸多方面都有廣泛的應(yīng)用前景.

[1]Yablonovitch E.Inhibited spontaneous emission in solidstate physics and electronics.Phys[J].Rev.Lett.,1987,58(20):2059

[2]John S.Strong localization of photons in certain disordered dielectric[J].Phys.Rev.Lett.,1987,58(23):2486.

[3]Kushwaha M S,Halevi P,DobrzynsiL.Acoustic band structure of periodic elastic composites[J].Physical Review Letters,1993,71(13):2022-2025.

[4]溫熙森,等著.光子/聲子晶體理論與技術(shù)[M].北京:科學(xué)出版社,2006,3-4.

[5]Vasseur J O,Djafarirouhani B,DobrZynski L,et al.Acoustic band gaps in fibre composite materials of boron nitride structure[J].J.Phys.Condens.Mat.,1997,9(35):7327.

[6]溫激鴻,王剛,郁殿龍,等.聲子晶體振動帶隙及減振特性研究[J].中國科學(xué)E輯,2007,37(9):1126-1139.

[7]溫熙森,溫激鴻,郁殿龍,等.聲子晶體[M].北京:國防工業(yè)出版社,2009;101.

[8]曹永軍,胡曉穎,周培勤.一維準(zhǔn)周期結(jié)構(gòu)聲子晶體帶隙特性的研究[J].內(nèi)蒙古師范大學(xué)學(xué)報(自然科學(xué)漢文版),2008,37(1):53-57.

[9]宿星亮,高原文.一維功能梯度材料聲子晶體彈性波帶隙研究[J].功能材料,2010,41(增刊):368-370.

[10]LIU Z,CHAN C T,SHENGP.Three-component elastic wave band-gap material[J].Phys.Rev.B,2002,65(16):165116.

[11]GOFFAUX C,DEHESA S J.Two-dimensi onal phononic crystals studied using a variati onalmethod:app licati on to lattices of l ocally res onantmaterials[J].Phys.Rev.B,2003,67(14):144301.

[12]H IRSEKORN M.Small-size sonic crystalswith strong attenuation bands in the audible frequency range[J].Appl.Phys.Lett,2004,84(17):3364-3366.

[13]GOFFAUX C,DEHESA S J,YEYAT IA L,et al.Evidence of fano-like interference phenomena in l ocallyres onantmaterials[J].Phys.Rev.Lett,2002,88(22):225502.

[14]VASSEUR J O.Experimental and theoretical evidencefor the existence of abs olute acoustic band gap s in twodi mensi onal s olid phononic crystals[J].Phys.Rev.Lett,2000,86:3012.

[15]PABLOS GD.Theory and experi mental on elastic bandgaps[J].Phys.Rev.Lett,2000,84:4349.

[16]VASSEUR J O.Phononic crystal with l ow filling fracti on and abs olute acoustic band gap in the audible frequency range:a theoretical and experimental study[J].Phys.Rev.E,2002,65(5):6608-6613.

[17]H IRSEKORN M,DELSANTO P P,BATRA N K,et al.Modelling and si mulati on of acoustic wave p r opagati on in l ocally resonant s onicmaterial[J].Ultras onics,2004,42(129):231-235.

[18]HOK M,CHENG C K,YANG Z,et al.Br oadband locally resonant sonic shields[J].Appl.Phys.Lett,2003,83(26):5566-5568.

[19]L IU Z,ZHANG X,MAO Y,et al.Locally resonants onic materials[J].Science,2000,289:734-1736.

[20]王剛,溫激鴻,韓小云,等.二維聲子晶體帶隙計算中的時域有限差分方法[J].物理學(xué)報,2003,52(8)

[21]Camly R E,Djafari-Rouhani B,Dobrzynski L.Transverse elastic waves in periodically layered infinite and semi-infinite media[J].Phys.Rev.B.,1993,27(12):7318.

[22]Sigalas M M,Soukoulis C M.Elastic-wave propagation through disordered and/or absorptive layered systems[J].Phys.Rev.B.,1995,51(5):2780.

[23]郁殿龍,劉耀宗,邱靜.表面局域態(tài)對一維聲子晶體中水平剪切波傳輸特性的影響[J].人工晶體學(xué)報,2005,34(3):425-430.

[24]蔣澤,趙琳,周建超.一維聲子晶體中聲波傳播的理論分析[J].壓電與聲光,2007,29(6):638-640.

[25]Vasseur J O,DjafarirouhaniB,DobrZynskiL,et al.Acoustic band gaps in fibre composite materials of boron nitride structure[J].J.Phys.Condens.Mat.,1997,9(35):7327.

[26]Vasseur J O,Djafari-Rouhani B,Dobrzynski L,et al.Complete acoustic band gaps in periodic fibre reinforced composite materials:the carbon/epoxy composite and some metallic systems[J].J.Phys:Condens Mat,1994,(6):8759.

[27]肖偉,曾廣武.基于波傳播法的一維聲子晶體禁帶[J].機械工程學(xué)報,2007,43(1):125-128.

[28]閆志忠,汪越勝.一維聲子晶體彈性波帶隙計算的小波方法[J].中國科學(xué)G輯,2007,37(4):544-551.

[29]王剛,溫激鴻,劉耀宗,等.一維粘彈材料周期結(jié)構(gòu)的振動帶隙研究[J].機械工程學(xué)報,2004,40(7):47-50.

[30]溫激鴻,王剛,劉耀宗,等.基于集中質(zhì)量法的一維聲子晶體彈性波帶隙計算[J].物理學(xué)報,2004,53(10):3384.

[31]溫激鴻,王剛,劉耀宗,等.大彈性常數(shù)差二維聲子晶體帶隙計算中的集中質(zhì)量法[J].物理學(xué)報,2005,54(03):1247.

[32]劉鐵權(quán),鄧子辰,周加喜.基于辛數(shù)學(xué)方法的一維聲子晶體禁帶計算[J].振動與沖擊,2010,29(12):102-105.

Research Review on the Band Gap of One-dimensional Phononic Crystals

QIU Xue-yun

(Department of Math and Physics,Wenshan University,Wenshan Yunnan 663000,China)

The article briefly summarizes the concept of photonic crystals,its basic characteristics and classification.It analyzes and concludes the three simplified models of the one-dimensional phononic crystals,the calculation methods and the research result of the band gap.It offers a basis for further research on one-dimensional phononic crystals by forecasting the possible applications in low frequency vibration,noise control,anti-seismic and the earthquake prevention.

one-dimensional phononic crystals;band structure;band gap

O469

A

1008-9128(2011)02-0005-04

2011-03-12

國家自然科學(xué)基金項目(10664006);云南省教育廳科研基金項目(2010Y093)。

邱學(xué)云(1979-),男,云南宣威人,助理研究員,理學(xué)碩士。研究方向:聲子晶體、凝聚態(tài)物理。

[責(zé)任編輯 張燦邦]