陳圣兵，韓小云，王 剛，溫激鴻

(1.國防科技大學(xué) 機(jī)電工程研究所，長沙 410073;2.國家光子/聲子晶體教育部重點(diǎn)實(shí)驗(yàn)室，長沙 410073)

聲子晶體概念的提出為減振降噪理論研究和工程設(shè)計(jì)提供了新的思路。聲子晶體是具有彈性波帶隙的周期性結(jié)構(gòu)功能材料，頻率落在帶隙范圍內(nèi)的振動(dòng)會(huì)被禁止［1］。根據(jù)帶隙頻率所對(duì)應(yīng)的波長與晶格常數(shù)的比例關(guān)系，聲子晶體可分為Bragg散射型和局域共振型。局域共振型聲子晶體中，單個(gè)散射體的共振特性起主導(dǎo)作用。通常是在基體材料中引入局域諧振子，通過振子的機(jī)械共振與基體中的長波行波相互耦合作用形成帶隙。

隨著智能材料的發(fā)展，將智能結(jié)構(gòu)引入聲子晶體，通過主動(dòng)控制，使聲子晶體的帶隙可調(diào)，已經(jīng)逐漸被人們重視。Ruzzene和Baz等人［2］首先將形狀記憶合金引入到周期結(jié)構(gòu)中，通過改變形狀記憶合金的彈性模量來調(diào)節(jié)帶隙。隨后他們又將壓電陶瓷材料及簡單控制電路周期性地附加在一維桿［3，4］、管路［5－7］和二維周期彈簧質(zhì)量系統(tǒng)［8］中，采用簡單的控制電路，對(duì)周期結(jié)構(gòu)的振動(dòng)帶隙進(jìn)行調(diào)節(jié)。李鳳明等［9，10］采用傳遞矩陣法研究了周期嵌有壓電材料的桿和梁結(jié)構(gòu)中的波動(dòng)傳播特性及其振動(dòng)主動(dòng)控制問題。張文群等［11］研究了負(fù)電容壓電分流阻尼系統(tǒng)的能量耗散特性，分析了不同連接方式電路的特點(diǎn)。孫浩和楊智春等［12］也對(duì)壓電分流阻尼系統(tǒng)中壓電片的拓?fù)鋬?yōu)化做了較為深入的研究。

本文提出在基體桿上周期性粘貼連有分流電路的壓電片，用分流電路中LC回路的電磁振蕩代替局域共振聲子晶體中傳統(tǒng)機(jī)械局域諧振子的振動(dòng)。利用壓電片的機(jī)電轉(zhuǎn)換能力和LC電路的諧振特性設(shè)計(jì)了局域共振聲子晶體桿結(jié)構(gòu)，運(yùn)用傳遞矩陣法計(jì)算了桿的帶隙結(jié)構(gòu)，并分析了桿的帶隙形成機(jī)理以及在LC電路中引入耗散電阻后對(duì)帶隙的影響。最后，采用商用有限元分析軟件ANSYS中的壓電模塊和電路模塊對(duì)桿的振動(dòng)特性進(jìn)行了仿真計(jì)算。

1 聲子晶體桿振動(dòng)的波動(dòng)模型

聲子晶體桿的結(jié)構(gòu)如圖1所示。在基體桿上周期性對(duì)貼壓電陶瓷片，形成周期結(jié)構(gòu)。因?yàn)閴弘娞沾善旧泶嬖陔娙?，所以采用并?lián)電容分流電路可以使振蕩電路中總電容等于外接電容和壓電片電容之和，與串聯(lián)電容分流電路相比，可以有效增大電路中總電容的調(diào)節(jié)范圍，而且可以使電路固有振蕩頻率與電容的調(diào)節(jié)關(guān)系更直觀。所以將壓電陶瓷片與電容C、電感L和電阻R并聯(lián)，形成并聯(lián)分流電路。

圖1 周期桿示意圖Fig.1 Schematic diagram of periodic rod

圖2 壓電分流電路Fig.2 Piezoelectric shunting circuit

1.1 壓電方程和桿的結(jié)構(gòu)動(dòng)力學(xué)方程

壓電分流電路如圖2所示，假設(shè)壓電片除去沿長度的兩端外其余表面均自由，且壓電片只受沿厚度方向的電場(chǎng)作用，則d－型壓電方程簡化為［13］:

其中T1、S1為壓電片在長度方向的應(yīng)力和應(yīng)變，D3為壓電片厚度方向表面電位移，E3為壓電片在厚度方向的電場(chǎng)強(qiáng)度，為彈性柔順常數(shù)，為厚度方向恒應(yīng)力介電常數(shù)，d31為壓電常數(shù)。

將電場(chǎng)和電位移公式代入方程(1)解得［4］:

其中:s為拉普拉斯算子，As為壓電片電極面積，hp為壓電片厚度，Cp為壓電片電容，

將s=iω代入式(2)得簡諧振動(dòng)下壓電片應(yīng)力與應(yīng)變的關(guān)系為:

其中i為虛數(shù)單位，ω為振蕩頻率。

即在分流電路作用下，壓電片的等效彈性模量變?yōu)?

其中Cto為電路中總電容，Cto=C+Cp。

取桿的一個(gè)基本周期單元，如圖3所示?；w桿的截面為矩形，厚度為hb，寬度為lb，截面面積為Ab，彈性模量為Eb，密度為ρb;壓電片的長度為b，寬度為lp，截面面積為Ap，密度為ρp。M表示未貼壓電片的桿，N表示貼有壓電片的桿。設(shè)y(x，t)為桿在x處的振動(dòng)位移，則桿的振動(dòng)方程為［14］:

式中:ρ、E為桿x處的密度和彈性模量，y為桿x處的振動(dòng)位移，y=X(x)eiωt。

圖3 單個(gè)周期單元Fig.3 Sketch map of the cell

由方程(5)解得，桿的第n個(gè)周期單元中M部分的振動(dòng)幅值在局部坐標(biāo)系下為:

其中:

桿的第n個(gè)周期單元中N部分振動(dòng)幅值在局部坐標(biāo)系下為:

其中:

1.2 桿單元的邊界條件

在周期單元n中N部分的左邊界處，位移X(x)、力EAX'(x)的連續(xù)性可得:

整理為矩陣形式為:

式中 ψn=［An1An2］T，γn=［Bn1Bn2］T。

在周期單元 n中 N的右邊界處，位移X(x)、力EAX'(x)的連續(xù)性可得:

整理為矩陣形式為:

式中:ψn=［A(n+1)1A(n+1)2］T，γn=［Bn1Bn2］T。

將式(8)代入式(9)得:

1.3 桿中波傳播的動(dòng)力學(xué)特性

桿是由相同的周期單元組成的，根據(jù) Bloch定理［15］:

其中μ為波傳播常數(shù)。

將式(10)代入式(11)得特征值方程:

其中:I為單位矩陣。

由方程可得，對(duì)于給定頻率ω，可以求出波傳播常數(shù)μ。當(dāng)頻率位于某些范圍時(shí)，波傳播常數(shù)實(shí)部為零，這表明周期桿中這些頻率范圍內(nèi)的彈性波將無衰減傳播。而波傳播常數(shù)實(shí)部不為零時(shí)，周期桿中彈性波僅存在衰減的傳播模式，這些頻率范圍形成帶隙。

1.4 桿中局域共振帶隙形成機(jī)理研究

壓電片與分流電路相互作用，其等效剛度變?yōu)槭?4)所示，整理為:

則桿上貼有壓電片的N部分的等效彈性模量為:

聯(lián)立式(13)、式(14)得，當(dāng) ω =ωzr時(shí)，Re(Ezr)=0。

式中:ω0—— 分流電路固有頻率。

此時(shí)桿中N部分的等效彈性模量實(shí)部為0，那么彈性波顯然不能在桿中傳播，必然形成帶隙。因此桿的局域共振帶隙在ωzr處，比電路振蕩固有頻率高。

2 數(shù)值計(jì)算結(jié)果

基體桿的材料選用環(huán)氧樹脂，其彈性模量Eb=4.35 ×1010Pa，密度為 ρb=1180 kg/m3，厚度為 hb=0.01 m，寬度 lb=0.02 m，a=0.03 m，b=0.04 m。壓電陶瓷片選用 P ZT －5H，d31= －274×10－12C/m2，=3.01 × 10－8F/m，=16.5 ×10－12m2/N，厚度 hp=0.0005 m，lp=0.02 m。

2.1 桿的帶隙結(jié)構(gòu)

由以上參數(shù)仿真計(jì)算得到桿的帶隙結(jié)構(gòu)如圖4所示。圖4(a)為電感L=0，R=10 Ω時(shí)的帶隙結(jié)構(gòu)，此時(shí)分流電路不能產(chǎn)生電磁振蕩，其帶隙由壓電陶瓷片與基體桿間的Bragg散射引起，為Bragg帶隙。圖4(b)為電感 L=7 ×10－3H，C=0，R=10 Ω 即 Cto=Cp時(shí)的帶隙結(jié)構(gòu)，此時(shí)分流電路中能產(chǎn)生電磁振蕩，其帶隙包括局域共振帶隙和Bragg帶隙。

圖4 桿的帶隙結(jié)構(gòu)Fig.4 Band gaps of the rod

2.2 電感、電容和電阻對(duì)局域共振帶隙的影響

從式(15)可以看出，調(diào)節(jié)分流電路中電感L和電容C的值，可改變局域共振帶隙位置，如圖5所示。圖5(a)為帶隙隨電感L的變化，局域共振帶隙頻率隨著電感值的增大而降低。當(dāng)L=3×10－3H時(shí)，局域共振帶隙和Bragg帶隙位置重合，合并為一個(gè)帶隙。圖5(b)為帶隙隨電容的變化，局域共振帶隙頻率也是隨著電容值的增大而降低。

圖6為帶隙隨電阻的變化，從圖中可以看出隨著電阻的增大局域共振帶隙衰減幅值是逐漸減小的，因?yàn)殡娮璧脑龃髸?huì)使電路中電流減小，減弱了LC回路的電磁振蕩作用或者說振蕩回路的Q值下降。

為了驗(yàn)證理論計(jì)算的正確性，利用ANSYS有限元軟件對(duì)有限周期的桿進(jìn)行了有限元計(jì)算。取6個(gè)周期進(jìn)行有限元計(jì)算，其中基體桿用solid45單元進(jìn)行網(wǎng)格劃分，壓電陶瓷片采用solid5單元，電感和電容采用circuit94單元。在桿沒有粘貼壓電片的一端以恒定位移1×10－5m進(jìn)行激勵(lì)，利用ANSYS的諧響應(yīng)分析進(jìn)行有限元計(jì)算，得到振動(dòng)在桿中傳到另一端的位移響應(yīng)如圖7 所示，其中 R=10 Ω，C=0，L=7 ×10－3H。

從ANSYS有限元計(jì)算結(jié)果可以看出在8980 Hz～9520 Hz出現(xiàn)了局域共振帶隙，與圖4所示的理論計(jì)算結(jié)果基本一致。

圖5 帶隙隨電感和電容的變化Fig.5 Band gaps’variations with the inductance and capacitance

圖6 帶隙隨電阻的變化Fig.6 Band gaps’variations with the resistance

圖7 ANSYS計(jì)算傳輸特性Fig.7 Transmission properties calculated by ANSYS

3 結(jié)論

通過對(duì)含有壓電分流電路聲子晶體桿的理論計(jì)算和ANSYS分析可以得到以下結(jié)論:

(1)桿的帶隙是Bragg散射和局域共振共同產(chǎn)生。當(dāng)分流電路中電感等于零時(shí)，電路中不能產(chǎn)生電磁振蕩，帶隙由基體梁與壓電片之間的Bragg散射形成，為Bragg帶隙。當(dāng)分流電路中電感不等于零時(shí)，分流電路的電磁振蕩會(huì)形成局域共振，產(chǎn)生局域共振帶隙。

(2)局域共振帶隙的帶隙頻率為貼有壓電片部分等效彈性模量實(shí)部為零所對(duì)應(yīng)的頻率，比分流電路的固有頻率略高。通過改變電容和電感參數(shù)，可以調(diào)節(jié)局域共振帶隙的位置。

(3)分流電路中的耗散電阻會(huì)使電路中電流減小，電磁振蕩變?nèi)跫凑袷幓芈返腝值下降，從而減小了局域共振帶隙在帶隙頻率內(nèi)的振動(dòng)幅值、增加了帶寬。

總之，利用壓電材料的機(jī)電轉(zhuǎn)換作用，激勵(lì)分流電路，通過電磁振蕩形成局域共振，為局域共振型聲子晶體設(shè)計(jì)提供了新的思路。與傳統(tǒng)機(jī)械振子局域共振聲子晶體相比，由于分流電路中電感、電容和電阻可以方便的調(diào)節(jié)，因而可以根據(jù)需要主動(dòng)控制局域共振帶隙頻率的位置。而且分流電路不參與桿的機(jī)械振動(dòng)，減小了振動(dòng)結(jié)構(gòu)的質(zhì)量。

［1］溫熙森.光子/聲子晶體理論與技術(shù)［M］.北京:科學(xué)出版社，2006.

［2］Ruzzene M，Baz A.Attenuation and localization of wave propagation in periodic using shape memory inserts［C］.Proceedings of SPIE，2000，3991:389－407.

［3］Baz A.Active control of periodic structures［J］.J.Vibration and Acoustics，2001，123:472－479.

［4］ThorpO，Ruzzene M，Baz A.Attenuation and localization of wave propagation in rods with periodic shunted piezoelectric patches［C］.Proceedings of SPIE，2001，4331:218 －238.

［5］Ruzzene M.Active control of wave propagation in periodic fluid-loaded shells［J］.Smart Materials and Structures，2001，10:893－906.

［6］Ruzzene M，Gu Z，Baz A M.Wavelet analysis of wave propagation in shells with periodic stiffeners［C］.Smart Structures and Materials 2001:SmartStructures and Integrated Systems，Newport Beach，CA，USA，2001.

［7］Thorp O， RuzzeneM， BazA. Attenuation ofwave propagation in fluid-loaded shells with periodic shunted piezoelectric rings［J］.Smart Materials and Structures，2005，14:594－604.

［8］Kim Y，Baz A M.Active control of a two-dimensional periodic structure［C］.Smart Structures and Materials 2004:Damping and Isolation.Proceedings of the SPIE，2004，5386:329－339.

［9］李鳳明，汪越勝.壓電周期結(jié)構(gòu)振動(dòng)主動(dòng)控制研究［J］.振動(dòng)工程學(xué)報(bào)，2004，17:828－830.

［10］Chen A L，Li F M.Localization of flexural waves in disordered periodic piezoelectric beam［J］.Journal of Sound and Vibration，2007，304:863 －874.

［11］張文群，張 萌，吳新躍.負(fù)電容壓電分流阻尼系統(tǒng)的能量耗散特性［J］.振動(dòng)與沖擊，2008，27(10):70－73.

［12］Sun H，Yang Z C，Li K X，et al.Vibration suppression of a hard disk driver actuator arm using piezoelectric shunt damping with a topology-optimized PZT transducer［J］.Smart Materials and Structures，2009，18:1－13.

［13］欒桂冬，張金鐸，王仁乾.壓電換能器和換能器陣［M］.北京:北京大學(xué)出版社，2005.

［14］溫激鴻.聲子晶體振動(dòng)帶隙及減振特性研究［D］.長沙:國防科技大學(xué)，2005.

［15］溫激鴻，郁殿龍，等.周期結(jié)構(gòu)細(xì)直梁彎曲振動(dòng)中的振動(dòng)帶隙［J］.機(jī)械工程學(xué)報(bào)，2005，41(4):1－7.