黃佳　尹進(jìn)　張昭　張洪武

摘要：基于聲子晶體特性提出一種層疊式方柱型聲子晶體單胞結(jié)構(gòu)，通過沖擊響應(yīng)譜分析和隨機(jī)振動分析考察其周期性結(jié)構(gòu)減振特性，并與相應(yīng)非周期結(jié)構(gòu)對比.計算結(jié)果表明該聲子晶體周期性結(jié)構(gòu)對于沖擊載荷和隨機(jī)振動載荷均有較好的減振特性.在三組分材料參數(shù)中影響單胞帶隙特性的主要因素是外層材料密度和中間層材料彈性模量.聲子晶體結(jié)構(gòu)的減振效果隨周期數(shù)的增加而愈加明顯.組分材料力學(xué)性能參數(shù)和周期型結(jié)構(gòu)周期數(shù)是聲子晶體結(jié)構(gòu)減振設(shè)計的主要因素.

關(guān)鍵詞：聲子晶體； 單胞結(jié)構(gòu)； 減振設(shè)計； 隔振材料； 低頻帶隙

中圖分類號： TH113； O328

文獻(xiàn)標(biāo)志碼：B

Abstract：A phononic crystal unit cell structure with stacked square columns is designed on the basis of the phononic crystal characteristics. The vibration reduction characteristics of the periodic structure are investigated through the shock response spectrum analysis and the random vibration analysis and compared with that of the corresponding non-periodic structure. The results show that the periodic structure of phononic crystal demonstrates better vibration reduction characteristics in the shock load and random vibration load problems. In the parameters of 3-components material， the density of the material in the outer layer and the elastic modulus of the material in the middle layer are the main factors that affect the band gap of unit cell. The vibration reduction of phononic crystal structure is obviously increased with the increase of the periodic numbers. The mechanical performance parameters of each component material and the periodic numbers of periodic structure are the main factors for the vibration reduction design of phononic crystal structure.

Key words：phononic crystal； unit cell structure； vibration reduction design； vibration isolation material； low frequency band gap

0引言

在航空航天、船舶等領(lǐng)域，沖擊和隨機(jī)振動是結(jié)構(gòu)設(shè)計中需要考慮的重要力學(xué)環(huán)境，如火箭發(fā)動機(jī)點(diǎn)火產(chǎn)生的瞬態(tài)沖擊[1-2]、水下武器對船舶攻擊時的非接觸爆炸[3]等，因此考慮沖擊載荷和隨機(jī)振動載荷作用的減振設(shè)計，對于結(jié)構(gòu)服役安全具有重要意義.

聲子晶體是一種空間周期性排布的多組分復(fù)合材料，具有阻擋某些頻率范圍內(nèi)彈性波和聲波傳播的特性，對應(yīng)的頻率范圍稱為帶隙或禁帶.[4]該周期性復(fù)合材料的最小單元稱為單胞.由聲子晶體單胞構(gòu)成的有限周期結(jié)構(gòu)稱為聲子晶體結(jié)構(gòu).作為由天然晶體中電子帶隙類推而來的周期結(jié)構(gòu)，聲子晶體可阻擋彈性波傳播的特性，使得彈性波特別是聲波在聲子晶體中的傳播特性研究受到越來越廣泛的關(guān)注.[5-6]聲子晶體的禁帶特性使得其在減振隔振[7]、降噪等領(lǐng)域都有潛在、廣闊的應(yīng)用前景.值得注意的是，美國國防部高級計劃局大力資助聲子晶體在減振降噪方面的研究，使聲子晶體的應(yīng)用范圍擴(kuò)大到飛行器發(fā)動機(jī)底座、航空電子設(shè)備底盤等.[4]對于聲子晶體的研究大多以帶隙特性為重點(diǎn)，以梁板類結(jié)構(gòu)[8-10]、理論分析和一維彎曲、扭轉(zhuǎn)及二維面內(nèi)面外振動帶隙[11-15]為主.

目前，對聲子晶體低頻減振特性的三維振動特性研究相對較少.研究三維聲子晶體的結(jié)構(gòu)特性對其進(jìn)一步工程化應(yīng)用具有現(xiàn)實意義.聲子晶體的尺寸和禁帶的頻率范圍是需要重點(diǎn)考慮的因素，此外，聲子晶體周期數(shù)對減振效果的影響以及聲子晶體單胞材料參數(shù)變化對帶隙特性的影響也不容忽視.鑒于此，建立一種具有隔振效果的層疊式方柱型聲子晶體三維模型，考察其周期結(jié)構(gòu)對沖擊載荷和隨機(jī)振動載荷的隔振效果，并與非周期結(jié)構(gòu)進(jìn)行對比.采用多重多級子結(jié)構(gòu)方法[16-18]考察其三維模型的振動帶隙特性和單胞材料參數(shù)對帶隙特性的影響，對比分析周期結(jié)構(gòu)和非周期結(jié)構(gòu)傳輸特性，考察周期數(shù)對結(jié)構(gòu)傳輸特性的影響.

1基于有限元的帶隙計算方法

聲子晶體研究中沿用晶格理論描述其空間結(jié)構(gòu)的周期性.周期結(jié)構(gòu)的最小重復(fù)單元稱為原胞，理想周期結(jié)構(gòu)是由原胞沿著3個不共面的基本平移矢量周期分布構(gòu)成，這3個基本平移矢量稱為基矢[19]，記為

2.2沖擊響應(yīng)譜振動分析

考察2種結(jié)構(gòu)的沖擊響應(yīng)譜振動特性.分別將聲子晶體結(jié)構(gòu)和非周期結(jié)構(gòu)的4根立柱底面設(shè)置為載荷輸入面，沿z方向輸入某固體火箭發(fā)動機(jī)點(diǎn)火時加速度沖擊響應(yīng)譜[2]，頻率范圍為10～370 Hz，輸入功率譜最大值約100g（g為重力加速度），可計算得到整體結(jié)構(gòu)響應(yīng).因上部結(jié)構(gòu)鋼板的響應(yīng)為考察對象，故單獨(dú)取鋼板模型及其結(jié)果進(jìn)行查看，見圖4.鋼板中部z方向加速度響應(yīng)很小，約為0.23g，整個鋼板最大加速度響應(yīng)也很小，約為0.33g.

取非周期結(jié)構(gòu)上部的鋼板響應(yīng)為考察對象，單獨(dú)取出鋼板模型及其結(jié)果進(jìn)行查看，見圖5.由此可知：鋼板中部z方向加速度響應(yīng)較大，約為32.4g；整個鋼板z方向最大加速度響應(yīng)位于鋼板中心部位.可見，聲子晶體結(jié)構(gòu)上部鋼板中心位置z方向加速度響應(yīng)僅為一般結(jié)構(gòu)的0.71%，鋼板最大加速度僅為一般結(jié)構(gòu)的1.02%.與非周期結(jié)構(gòu)相比，聲子晶體結(jié)構(gòu)響應(yīng)值更小，即彈性波穿過聲子晶體結(jié)構(gòu)后具有更小的能量.因此，本文提出的聲子晶體結(jié)構(gòu)具有更好的抗沖擊能力，可在一定程度上減弱沖擊載荷對其周圍結(jié)構(gòu)的不利影響.

2.3隨機(jī)振動功率譜分析

考察聲子晶體結(jié)構(gòu)對隨機(jī)振動載荷環(huán)境的減振效果.分別將聲子晶體結(jié)構(gòu)和非周期結(jié)構(gòu)4根立柱底面設(shè)置為載荷輸入面.隨機(jī)振動功率譜輸入選取某運(yùn)載火箭貯箱增壓系統(tǒng)隨機(jī)振動工作條件[20]，載荷沿z方向.選取上部鋼板結(jié)構(gòu)中心點(diǎn)z方向加速度響應(yīng)結(jié)果進(jìn)行分析.聲子晶體結(jié)構(gòu)和非周期結(jié)構(gòu)鋼板中心位置z方向功率譜密度（Power Spectral Density， PSD）見圖6.PSD描述隨機(jī)振動條件下結(jié)構(gòu)振動能量密度在頻域上的分布特性.由圖6可知，在給定的隨機(jī)振動條件下，聲子晶體結(jié)構(gòu)和非周期結(jié)構(gòu)上部鋼板振動響應(yīng)能量密度主要集中在20 Hz以上，其中非周期結(jié)構(gòu)PSD最大值為1 641.33 Hz-1，最小值為1.66×10-6 Hz-1；聲子晶體結(jié)構(gòu)PSD最大值為7.03 Hz-1，最小值為7.64×10-20 Hz-1.在該頻率范圍內(nèi)，聲子晶體結(jié)構(gòu)對應(yīng)的PSD普遍顯著小于非周期結(jié)構(gòu)，即隨機(jī)振動使聲子晶體結(jié)構(gòu)產(chǎn)生的振動在相應(yīng)頻率上分布的能量密度小于非周期結(jié)構(gòu).值得注意的是，在60 Hz以上的頻率范圍中，聲子晶體結(jié)構(gòu)產(chǎn)生隨機(jī)振動的PSD甚至可以比非周期結(jié)構(gòu)低20個數(shù)量級.

為使分析結(jié)果更直觀，計算2種結(jié)構(gòu)的PSD均方根值的積分（Cumulative Root-Mean-Square， CRMS）.PSD的CRMS表征能量概念.計算得到聲子晶體結(jié)構(gòu)和非周期結(jié)構(gòu)鋼板中心位置z方向PSD的CRMS結(jié)果見圖7.作為對隨機(jī)振動響應(yīng)量能量大小的衡量標(biāo)準(zhǔn)，PSD的CRMS的大小可以反映結(jié)構(gòu)對隨機(jī)振動的響應(yīng)情況.在給定的隨機(jī)振動輸入條件下，聲子晶體結(jié)構(gòu)上部鋼板的響應(yīng)能量遠(yuǎn)小于非周期結(jié)構(gòu)的，其中在50 Hz以上的頻率范圍中，聲子晶體結(jié)構(gòu)的PSD的CRMS值一直小于非周期結(jié)構(gòu)，且二者差值越來越大，這與PSD及其CRMS之間的積分關(guān)系一致.由2種結(jié)構(gòu)PSD及其CRMS的結(jié)果對比可知，聲子晶體結(jié)構(gòu)能夠顯著減小隨機(jī)振動對上部結(jié)構(gòu)的影響，說明所提出的聲子晶體結(jié)構(gòu)具有更好的減振特性.

3單胞帶隙特性計算與分析

3.1單胞帶隙特性計算

為進(jìn)一步對本文提出的聲子晶體隔振特性進(jìn)行研究，計算該聲子晶體單胞帶隙特性，并分析其影響因素.對前文中的單胞形式進(jìn)行帶隙特性計算，得到帶隙特性，見圖8.由此可知：在80～310 Hz出現(xiàn)寬度很大的禁帶，說明層疊型方柱聲子晶體單胞具有良好的帶隙特性，可以使頻率范圍較寬且頻率較低的彈性波大幅衰減.這可以解釋該聲子晶體的周期型結(jié)構(gòu)能夠?qū)ο鄳?yīng)的沖擊和隨機(jī)振動問題產(chǎn)生明顯減振效果的現(xiàn)象.

3.2材料參數(shù)對單胞帶隙特性的影響

為探討該型聲子晶體單胞帶隙的可調(diào)控性，研究單胞材料參數(shù)對單胞帶隙特性的影響.當(dāng)中心質(zhì)量塊與中間層材料不變時，考察外層金屬材料參數(shù)對該聲子晶體單胞帶隙特性的影響.外層金屬主要材料參數(shù)有彈性模量、密度和泊松比.分別計算單胞帶隙上下邊界和帶隙中心頻率隨3個材料參數(shù)的變化規(guī)律，計算結(jié)果見圖9a～9c.由此可以看出：外層金屬的彈性模量和泊松比對帶隙幾乎沒有影響，外層金屬的密度只影響帶隙的下邊界；隨著外層金屬密度變大為1 600～31 600 kg/m3，帶隙下邊界先降低后趨于不變，為66～146 Hz，從而帶隙中心頻率也呈現(xiàn)先降低而后趨于不變的趨勢，范圍為230～190 Hz.

當(dāng)中心質(zhì)量塊與外層金屬材料不變時，考察中間層材料參數(shù)對該聲子晶體單胞帶隙特性的影響.

中間層主要材料參數(shù)包括彈性模量、密度和泊松比，分別計算單胞帶隙上下邊界和帶隙中心頻率隨3個材料參數(shù)的變化規(guī)律，計算結(jié)果見圖9d～9f.由此可知，中間層的彈性模量對帶隙的上下邊界均有影響，中間層彈性模量（0.5～3 MPa）越大，帶隙上下邊界越高，上邊界為223～540 Hz，下邊界為56～137 Hz，且上邊界增高的程度隨彈性模量增大而更顯著.因此，中間層彈性模量越大，帶隙寬度Δf越大，為167～403 Hz，帶隙中心頻率fc越高，為139.5～338.5 Hz，帶隙寬度變化程度更為顯著，即相對帶隙rbg=Δf/fc越大；中間層泊松比變大，帶隙中心頻率稍有提高，但帶隙寬度變化不大；中間層的材料密度對帶隙幾乎沒有影響.可見，中間層材料對于該單胞帶隙特性影響較為顯著的參數(shù)為彈性模量.

4聲子晶體結(jié)構(gòu)傳輸特性計算分析

4.1傳輸特性計算

為驗證單胞帶隙特性并更直觀地展示其衰減彈性波的作用，考察聲子晶體周期結(jié)構(gòu)和非周期結(jié)構(gòu)的傳輸特性.分別將聲子晶體結(jié)構(gòu)和非周期結(jié)構(gòu)4根立柱底面設(shè)置為載荷輸入面，輸入載荷為z方向單位加速度，頻率范圍為0～400 Hz.考察鋼板中心點(diǎn)位置z方向加速度響應(yīng)a/a0.計算得到2種結(jié)構(gòu)傳輸特性曲線見圖10.

參考結(jié)構(gòu)（非周期結(jié)構(gòu)）相對加速度對數(shù)值在0～400 Hz的大部分頻率范圍內(nèi)都大于0，即參考結(jié)構(gòu)在對應(yīng)頻率范圍內(nèi)會將載荷輸入放大.聲子晶體結(jié)構(gòu)在本文考慮的頻率范圍內(nèi)的絕大頻域上產(chǎn)生的相對加速度對數(shù)值都小于0，即聲子晶體結(jié)構(gòu)在對應(yīng)頻率范圍內(nèi)會將載荷輸入衰減.其中，在30～400 Hz頻率范圍內(nèi)存在寬頻率段衰減現(xiàn)象，衰減范圍與上文計算得到的聲子晶體單胞帶隙特性一致，且有限周期聲子晶體結(jié)構(gòu)在比帶隙范圍更大的頻域范圍內(nèi)也存在使載荷輸入衰減的效果，因此聲子晶體具有良好的減振作用.

4.2周期數(shù)對振動特性的影響

為考察聲子晶體周期數(shù)量對周期結(jié)構(gòu)隔振效果的影響，分別建立1周期和2周期聲子晶體模型并計算其傳輸特性，有限元模型見圖11.

分別對1，2和3周期時聲子晶體結(jié)構(gòu)傳輸特性進(jìn)行計算，得到其傳輸特性曲線，見圖12.3種周期數(shù)的聲子晶體結(jié)構(gòu)z方向加速度均在50 Hz處開始衰減，在50～310 Hz范圍內(nèi)有很顯著的衰減效果，這

與上述計算得到的單胞帶隙特性一致.3種聲子晶體結(jié)構(gòu)相對加速度最小值均出現(xiàn)在200～250 Hz范

圍內(nèi)，1周期結(jié)構(gòu)加速度響應(yīng)最小值較輸入小6個數(shù)量級，2周期結(jié)構(gòu)小10個數(shù)量級，3周期結(jié)構(gòu)小

12個數(shù)量級.

由此可見，周期數(shù)分別為1，2和3時，聲子晶體結(jié)構(gòu)在帶隙范圍內(nèi)均有減振效果，周期數(shù)越大，減振效果越顯著.在對于減振要求不高的情況下，單胞形式的1周期聲子晶體結(jié)構(gòu)即可實現(xiàn)一定的減振功能.

5結(jié)論

基于聲子晶體概念建立層疊式方柱型隔振材料單胞，具有能夠衰減較寬的低頻彈性波的特性.通過對該型單胞的帶隙特性及其周期結(jié)構(gòu)振動特性分析，得到以下結(jié)論.

1）層疊式方柱型聲子晶體單胞的有限周期結(jié)構(gòu)具有很好的隔振效果，可以大幅衰減振動頻率在其單胞帶隙范圍內(nèi)的彈性波，較一般非周期結(jié)構(gòu)有明顯的隔振效果.

2）層疊式方柱型聲子晶體單胞具有良好的帶隙特性，可以使相應(yīng)頻率范圍內(nèi)的彈性波大幅衰減.

3）聲子晶體結(jié)構(gòu)中影響帶隙特性的主要設(shè)計參數(shù)為外層材料密度和中間層材料彈性模量.

4）周期聲子晶體結(jié)構(gòu)具有良好的減振特性，并且隨著周期數(shù)的增加，減振效果愈加明顯.

參考文獻(xiàn)：

[1]易磊， 寇軍強(qiáng)， 李衛(wèi)鵬， 等. 戰(zhàn)術(shù)固體火箭發(fā)動機(jī)點(diǎn)火沖擊試驗研究[J]. 彈箭與制導(dǎo)學(xué)報， 2010， 30（6）： 150-151.

YI Lei， KOU Junqiang， LI Weipeng， et al. Experiment research on ignition impulsion of solid rocket motor[J]. J Projectiles， Rockets， Missiles & Guidance， 2010， 30（6）： 150-151.

[2]王翠榮， 施廣富， 郭軍. 固體火箭發(fā)動機(jī)沖擊信號響應(yīng)譜分析[J]. 固體火箭技術(shù)， 2003， 26（2）： 57-60.

WANG Cuirong， SHI Guangfu， GUO Jun. Analysis of impact signal response spectrum of solid rocket motors[J]. J Solid Rocket Technol， 2003， 26（2）： 57-60.

[3]石菲. 船舶發(fā)動機(jī)排氣消聲器抗沖擊性能研究[D]. 哈爾濱： 哈爾濱工程大學(xué)， 2006.

[4]郁殿龍. 基于聲子晶體理論的梁板類周期結(jié)構(gòu)振動帶隙特性研究[D]. 長沙： 國防科學(xué)技術(shù)大學(xué)， 2006.

[5]LIU Z， ZHANG X， MAO Y， et al. Locally resonant sonic materials[J]. Scince， 2000（289）： 1734-6.

[6]MALDOVAN M. Sound and heat revolutions in phononics[J]. Scince， 2013（503）： 209-217.

[7]溫激鴻， 王剛， 郁殿龍， 等. 聲子晶體振動帶隙及減振特性研究[J]. 中國科學(xué)： E輯： 技術(shù)科學(xué)， 2007， 37（9）： 1126-1139.

WEN Jihong， WANG Gang， YU Dianlong， et al. Vibration attenuation and band gap characteristics of phononic crystals[J]. Scientia Sinica： Ser E： Technol Sci， 2007， 37（9）： 1126-1139.

[8]宋卓斐， 王自東， 王艷林， 等. 一維桿狀聲子晶體的帶隙特性[J]. 振動與沖擊， 2010， 29（2）： 145-148.

SONG Zhuofei， WANG Zidong， WANG Yanlin， et al. Band gap property of a one dimension rod phononic crystal[J]. J Vibration & Shock， 2010， 29（2）： 145-148.

[9]盧天健， 高國欽， 馬守林. 二維四方排列半圓鋁管/空氣聲子晶體的禁帶特性[J]. 中國科學(xué)： E輯： 技術(shù)科學(xué)， 2009， 39（1）： 57-64.

LU Tianjian， GAO Guoqin， MA Shoulin. Properties of band gaps in 2D half Aluminum rods of hollows-air phononic crystal with square lattice[J]. Scientia Sinica： Ser E： Technol Sci， 2009， 39（1）： 57-64.

[10]丁紅星， 沈中華， 李加， 等. 復(fù)合蘭姆波聲子晶體中超寬部分禁帶[J]. 物理學(xué)報， 2012， 61（19）： 394-401 .

DING Hongxing， SHEN Zhonghua， LI Jia， et al. Large partial band-gaps for Lamb waves in multiple phononic crystals thin plates[J]. Acta Phisica Sinica， 2012， 61（19）： 394-401.

[11]SHENG P， ZHANG X， LIU Z， et al. Locally resonant sonic materials[J]. Physica B： Condensed Mat， 2003， 338（1-4）： 201-205.

[12]文岐華， 左曙光， 魏歡. 多振子梁彎曲振動中的局域共振帶隙[J]. 物理學(xué)報， 2012， 61（3）： 240-246.

WEN Qihua， ZUO Shuguang， WEI Huan. Locally resonant elastic wave band gaps in flexural vibration of multi-oscillators beam[J]. Acta Phisica Sinica， 2012， 61（3）： 240-246.

[13]舒海生， 劉少剛， 王威遠(yuǎn)， 等. 集中質(zhì)量邊界條件下聲子晶體桿的縱向振動傳遞特性研究[J]. 振動與沖擊， 2012， 31（19）： 113-117.

SHU Haisheng， LIU Shaogang， WANG Weiyuan， et al. Transmission characteristics of longitudinal vibration of a phononic crystal rod with concentrated mass boundary condition[J]. J Vibration & Shock， 2012， 31（19）： 113-117.

[14]YU D， WEN J， ZHAO H， et al. Vibration reduction by using the idea of phononic crystals in a pipe-conveying fluid[J]. J Sound & Vibration， 2008， 318（1-2）： 193-205.

[15]ZHANG S， WU J， HU Z. Low-frequency locally resonant band-gaps in phononic crystal plates with periodic spiral resonators[J]. J Appl Phys， 2013， 113（16）： 163511（1-8）.

[16]鐘萬勰. 一個多用途的結(jié)構(gòu)分析程序JIGFEX（一）[J]. 大連工學(xué)院學(xué)報， 1977， 17（3）： 19-42.

ZHONG Wanxie. A general-purpose program for structure analysis JIGFEX（Ⅰ）[J]. J Dalian Univ Technol， 1977， 17（3）： 19-42.

[17]鐘萬勰. 一個多用途的結(jié)構(gòu)分析程序JIGFEX（二）[J]. 大連工學(xué)院學(xué)報， 1977， 17（4）： 14-35.

ZHONG Wanxie. A general-purpose program for structure analysis JIGFEX（Ⅱ）[J]. J Dalian Univ Technol， 1977， 17（4）： 14-35.

[18]LU X， LIN J， ZHONG W. Subspace iteration method in multi-level substructure systems[J]. Comput & Structures， 1989， 33（2）： 459-462.

[19]溫熙森， 溫激鴻， 郁殿龍， 等. 聲子晶體[M]. 北京： 國防工業(yè)出版社， 2009： 29-41.

[20]王翀， 梁國柱. 隨機(jī)振動條件火箭貯箱增壓系統(tǒng)工作過程仿真[J]. 北京航空航天大學(xué)學(xué)報， 2013， 39 （10）： 1314-1318.

WANG Chong， LIANG Guozhu. Simulation of the working process of launch vehicle pressurization system in random vibration condition[J]. J Beijing Univ Aeronautics & Astronautics， 2013， 39 （10）： 1314-1318.

（編輯武曉英）