寧榮輝 朱石堅(jiān) 翁雪濤 張振海

(海軍工程大學(xué)動(dòng)力工程學(xué)院 武漢 430033)

0 引 言

結(jié)構(gòu)中的振動(dòng)通常是以彈性波的形式傳播的，彈性波在周期性介質(zhì)中傳播時(shí)，會(huì)產(chǎn)生彈性波帶隙[1]，由于彈性波帶隙的存在，使得周期結(jié)構(gòu)在減振降噪、聲濾波器和新型換能器等方面具有廣泛的應(yīng)用前景[2]，因此受到了人們的廣泛關(guān)注.

彈性波在一維周期結(jié)構(gòu)中的傳播特性已經(jīng)研究較多.溫激鴻等[3]對(duì)周期細(xì)直梁的彎曲振動(dòng)特性進(jìn)行了研究；祁鵬山等[4]研究了一維雙周期結(jié)構(gòu)的帶隙特性；李雁飛等[5]對(duì)周期管路的傳播特性和低頻特性進(jìn)行了研究.在研究周期結(jié)構(gòu)的帶隙特性時(shí)，都是針對(duì)無限周期結(jié)構(gòu)進(jìn)行的.而實(shí)際運(yùn)用的周期結(jié)構(gòu)都是有限的.因此對(duì)有限周期結(jié)構(gòu)進(jìn)行相關(guān)研究，確定衰減頻帶的位置和帶寬是十分有必要的.

目前周期結(jié)構(gòu)帶隙特性的計(jì)算方法主要有傳遞矩陣法，平面波展開法，時(shí)域有限差分法[6]，多重散射法[7]，有限元法[8]以及譜單元法[9]，其中平面波展開法、時(shí)域差分法、多重散射法主要用于無限長和半無限長結(jié)構(gòu)的帶隙特性計(jì)算，而有限元法主要用于有限長周期結(jié)構(gòu)的帶隙特性分析，其計(jì)算精度受網(wǎng)格密度和網(wǎng)格質(zhì)量影響較大；而譜單元法結(jié)合了有限元方法的離散和整合的特點(diǎn)，將結(jié)構(gòu)的運(yùn)動(dòng)方程通過Fourier變換轉(zhuǎn)到頻域上，在頻域上對(duì)結(jié)構(gòu)的各個(gè)頻率下運(yùn)動(dòng)進(jìn)行求解，然后將所有不同頻率的波進(jìn)行疊加即得到整體的運(yùn)動(dòng)解.對(duì)于每一個(gè)頻率來說，其解得的結(jié)果是精確的，因此譜單元法適合用于對(duì)有限周期結(jié)構(gòu)的計(jì)算.

本文利用譜單元法結(jié)合Bloch理論對(duì)一維周期梁結(jié)構(gòu)的單個(gè)元胞進(jìn)行了分析，得到無限周期結(jié)構(gòu)的復(fù)能帶結(jié)構(gòu)，然后結(jié)合有限長周期結(jié)構(gòu)的邊界條件，計(jì)算了周期梁結(jié)構(gòu)彎曲振動(dòng)傳遞率，利用彎曲振動(dòng)傳遞率分析了有限周期結(jié)構(gòu)梁的周期數(shù)、晶格長度、梁的截面對(duì)結(jié)構(gòu)的衰減頻帶的影響.

1 周期梁結(jié)構(gòu)振動(dòng)特性計(jì)算方法

一維周期結(jié)構(gòu)梁單個(gè)元胞的示意圖見圖1，單個(gè)元胞一般是由兩種截面形狀不同或材料參數(shù)不同的梁組成，假設(shè)梁的參數(shù)滿足長度遠(yuǎn)大于截面尺寸條件，則各梁滿足歐拉-伯努利梁條件.兩段梁分別長為L1，L2，單個(gè)元胞的長度為L=L1+L2，截面形狀為矩形，截面尺寸分別為b1，h1和b2，h2，元胞左右兩端的位移分別為向量uL，uR，受力也分別為向量fL、fR，中間連接處的位移和受力分別為向量uI，fI.

圖1 周期梁結(jié)構(gòu)單個(gè)元胞示意圖

基于譜元法的建模思想，單個(gè)元胞中的兩段梁可以分別處理為譜單元，考慮梁的彎曲問題時(shí)，因此上述的位移和力向量都有兩個(gè)分量，第一段梁的動(dòng)力學(xué)方程為

(1)

其中：剛度矩陣的各個(gè)子矩陣分別為

(2)

式中：

α=

β=

(3)

同理可得第二段梁的動(dòng)力學(xué)方程，然后由位移連續(xù)條件有：uI1=uI2=uI，fI1=-fI2.于是元胞的動(dòng)力學(xué)方程為

(4)

化簡得到只含有左右兩端的位移和受力的表達(dá)式為

(5)

根據(jù)周期結(jié)構(gòu)的Bloch定理，元胞的邊界位移和力向量滿足以下關(guān)系：

uR=e-iqLuL，fR=-e-iqLfL

(6)

將式(6)代入式(5)得到只含有位移的特征方程：

(7)

求解式(7)的特征值，可以得到兩對(duì)大小相同符合相反的Bloch波矢解±q1，±q2以及e-iqL相對(duì)應(yīng)的特征向量.取e-iq1L，e-iq2L，eiq1L，eiq2L對(duì)應(yīng)的特征向量分別記為φ1，φ2，φ3，φ4.于是有

uL=a1φ1+a2φ2+a3φ3+a4φ4

(8)

利用式(5)和Bloch定理，可得

a1Φ1+a2Φ2+a3Φ3+a4Φ4

(9)

對(duì)于有限周期的梁結(jié)構(gòu)，其周期數(shù)記為N，由Bloch定理可得

uNR=a1e-iNq1Lφ1+a2e-iNq2Lφ2+

a3eiNq1Lφ3+a4eiNq2Lφ4

(10)

式(8)和式(10)寫成矩陣形式為

(11)

同理可得

(12)

綜合式(11)和式(12)可得

(13)

式中:

(14)

2 有限長周期梁結(jié)構(gòu)帶隙特性及振動(dòng)傳遞率的計(jì)算

根據(jù)Bloch波矢解可以得到周期結(jié)構(gòu)的帶隙，由此得出的帶隙是由該元胞組成的無限周期結(jié)構(gòu)的帶隙，但是實(shí)際上結(jié)構(gòu)都是有限的，因此僅通過Bloch波矢并不能有效的反應(yīng)有限周期結(jié)構(gòu)對(duì)彈性波傳播的作用.對(duì)此本文通過計(jì)算有限周期梁結(jié)構(gòu)的振動(dòng)傳遞率來反應(yīng)有限周期梁結(jié)構(gòu)對(duì)彈性波的作用，結(jié)合有限周期梁結(jié)構(gòu)的邊界條件和式(13)可以計(jì)算得到有限長周期結(jié)構(gòu)梁的振動(dòng)傳遞率，其計(jì)算過程為

假設(shè)有限長周期梁結(jié)構(gòu)左端受振幅為wL0的簡諧橫向位移激勵(lì)，此外不受其他任何激勵(lì)和約束，則有邊界條件：

uL={wL0,θL0}T，fL={VL,0}T，fNR=0

(15)

將式(15)代入式(13)可得

(16)

式中取二階方陣AN、BN使得

(17)

則可得有限結(jié)構(gòu)的位移傳遞率為

(18)

本文以鋁制的變截面周期結(jié)構(gòu)梁為例，依據(jù)元胞模型和鋁的結(jié)構(gòu)參數(shù)，利用式(11)計(jì)算得到周期結(jié)構(gòu)梁在無限周期條件下的彎曲振動(dòng)帶隙.然后結(jié)合周期結(jié)構(gòu)的周期數(shù)和邊界條件，利用式(18)計(jì)算得到有限周期下梁的彎曲振動(dòng)傳遞率.結(jié)構(gòu)模型的具體參數(shù)如下：元胞周期長L=0.1 m，元胞中兩段梁長度比值L1:L2=1:1，截面尺寸b1=0.01 m，h1=0.01 m，b2=0.01 m，h2=0.005 m，鋁的材料參數(shù)為密度ρ=2799 kg/m3，彈性模量E=72 GPa.

計(jì)算得到的無限周期結(jié)構(gòu)梁的復(fù)能帶結(jié)構(gòu)見圖2，當(dāng)歸一化波矢qL的實(shí)部為0或±π，虛部不為0時(shí)，描述的本征波是衰減波，即該頻率是處于衰減范圍，從圖中可以看到，無限周期結(jié)構(gòu)梁在10 000 Hz以下存在兩個(gè)帶隙，分別為1 160～1 640 Hz和5 200～7 440 Hz.

圖2 無限周期結(jié)構(gòu)梁復(fù)能帶結(jié)構(gòu)圖

圖3為5周期和15周期兩種情況下梁結(jié)構(gòu)的彎曲振動(dòng)傳遞率曲線，兩種情況下振動(dòng)傳遞率曲線都有兩個(gè)衰減頻帶，其中5周期結(jié)構(gòu)的衰減頻帶分別為900～2 240 Hz，5 020～8 280 Hz；15周期結(jié)構(gòu)的衰減頻帶分別為1 090～1 710 Hz，5 130～7 520 Hz，對(duì)比無限周期結(jié)構(gòu)的帶隙，可以發(fā)現(xiàn)，有限周期結(jié)構(gòu)的衰減頻帶要比無限周期的帶隙略寬，并且隨著結(jié)構(gòu)周期數(shù)目的增加，衰減頻帶逐漸向帶隙靠近.由圖3可知，隨著周期數(shù)的增加，有限結(jié)構(gòu)的衰減頻帶雖然減小了，但是衰減頻帶的振動(dòng)傳遞率更低了，因此在設(shè)計(jì)周期結(jié)構(gòu)時(shí)應(yīng)綜合考慮振動(dòng)的衰減程度和頻帶范圍來確定結(jié)構(gòu)的周期數(shù).

圖3 有限周期結(jié)構(gòu)梁振動(dòng)傳遞率

3 元胞參數(shù)對(duì)衰減頻帶影響

3.1 Bragg頻率計(jì)算

一維周期系統(tǒng)嚴(yán)格對(duì)應(yīng)的Bragg條件為

(19)

(20)

綜合式(19)～(20)可得

(21)

保持各段梁的截面為矩形并寬度保持不變(即b1、b2為常數(shù))，取兩段梁的長度比為L1:L2=ε，梁高的比為h1:h2=η，式(21)可以化簡為

(22)

由式(22)可知，當(dāng)ε和η不變時(shí)，由兩段梁組成的周期結(jié)構(gòu)的Bragg頻率與元胞長度的二次方成反比，與元胞的厚度成正比.對(duì)于有限周期結(jié)構(gòu)來說，其衰減頻帶的起始頻率和終止頻率應(yīng)該也保持這個(gè)規(guī)律.

3.2 元胞長度對(duì)衰減頻帶的影響

仍然考慮鋁制變截面梁結(jié)構(gòu)，首先保持ε=1和η=2不變，截面尺寸為b1=0.01 m，b2=0.01 m，h1=0.01 m，則h2=0.005 m，改變晶格長度L的值，計(jì)算周期數(shù)為5.

計(jì)算了元胞長度分別為0.08，0.1，0.12，0.15 m時(shí)的振動(dòng)傳遞率曲線見圖4.由圖4可知，隨著晶格長度的增加，衰減頻帶向低頻范圍移動(dòng)，且?guī)捊档?，但是在周期?shù)不變的情況下，衰減幅度改變并不明顯.

圖4 元胞長度對(duì)衰減頻帶的影響

表1為各個(gè)晶格長度情況下，第一衰減頻帶的起始頻率fs、截止頻率fe以及與元胞長度二次方的乘積.由表可知，衰減頻帶的起始頻率和截止頻率與元胞長度的二次方成反比.

表1 衰減頻帶與元胞長度的變化

3.3 元胞厚度對(duì)衰減頻帶的影響

依然保持ε=1和η=2不變，元胞長為0.1m，截面尺寸為b1=0.01m，b2=0.01 m，改變梁的厚度h1，h2則相應(yīng)變化，計(jì)算周期數(shù)為5.

計(jì)算了厚度h1分別為0.008，0.01，0.012，0.015 m時(shí)的振動(dòng)傳遞率見圖5，由圖5可知，隨著厚度的增加，衰減頻帶向高頻范圍移動(dòng)，且?guī)捲黾?，但是在周期?shù)不變的情況下，衰減幅度改變并不明顯.

圖5 元胞厚度對(duì)衰減頻帶的影響

表2為不同厚度情況下，第一衰減頻帶的起始頻率fs、終止頻率fe以及與厚度h1的比值.由表可知，在不同h1的情況下起始頻率與h1的比值幾乎保持不變，截止頻率與h1的比值也幾乎保持不變，即衰減頻帶的起始頻率和截止頻率與厚度h1成正比.

表2 衰減頻帶與厚度h1的變化

4 結(jié) 論

1) 有限周期梁結(jié)構(gòu)的衰減頻帶比無限周期的帶隙略寬，但是隨著周期數(shù)目的增加，衰減頻帶逐漸向帶隙范圍靠近.

2) 有限周期梁結(jié)構(gòu)衰減頻帶的衰減幅值與周期數(shù)有關(guān)，周期數(shù)越多衰減效果越好.

3) 當(dāng)有限周期梁結(jié)構(gòu)的材料、元胞內(nèi)各長度和厚度的比值確定時(shí)，不同元胞長度、不同厚度下的起止頻率的比值與元胞長度比值的二次方成反比與厚度的比值成正比.