層狀周期結構動力衰減域特性研究

程志寶，石志飛，向宏軍

(北京交通大學 土木建筑工程學院，北京 100044)

作為復合材料/結構(簡稱周期性材料或周期性結構)的聲子晶體材料/結構因其具有周期排列特性而廣受關注［1-3］。其存在的頻率禁帶(衰減域)及通帶(非衰減域)特性濾波性質，可為各種濾波器研制、工程結構隔震/減振降噪材料設計與開發提供新途徑［4］。

在土木工程領域，各種形式動力作用直接影響建筑物的安全性與舒適性。尤其地震給人類生命財產的安全帶來極大威脅［5］。為減弱地震對建筑結構造成的損害，出現了各種形式的結構減、隔震技術。其中基礎隔震為實踐證明且較有效的抗震措施［6］。石志飛等［7-9］將固體物理學中有關周期性結構特性引入土木工程領域，提出的研制周期性隔震基礎設想及有益探索表明，利用周期性基礎進行基礎隔震可行。

本文通過推導層狀周期結構色散方程，建立衰減域邊界控制方程并給出物理解釋;分析層狀周期結構帶隙的產生條件，研究各種因素對衰減域三個控制參數(衰減域起始頻率、衰減域寬度、衰減因子)影響。用有限元軟件，數值研究有限層狀周期結構的諧響應及瞬態響應，并制作隔震試件完成振動臺實驗，進一步證實層狀周期結構的動力衰減特性。該研究可為將層狀周期結構應用于工程結構隔震、減振提供基礎。

1 基本理論－色散帶隙及衰減因子

兩種均勻材料組成的一維層狀周期結構見圖1(a)。設結構沿x-y平面內無限大，沿z方向無限周期性交替排布，兩種材料界面完全粘結。取x-z平面內沿z方向傳播的剪切波Sx進行研究，分析周期性結構中波的色散特性。利用周期性理論，取其中一個單元見圖1(b)。

圖1 層狀周期結構及計算單元示意圖Fig.1 Configuration of layer periodic structure and the Wigner-Seitz cell

設u為質點沿x的振動位移，ρ為材料密度;λ和μ為材料的拉梅常數。剪切波Sx的控制方程為:

應用傅里葉變換，兩種材料在局部坐標系下分別求解式(1)，可得其位移與剪應力穩態解:

其中:k=α+iβ為復數波矢，α為相位因子，β為衰減因子［12］。式(3)的系數矩陣行列式等于零，由此可得色散曲線解析表達式為:

采用工程中常用的混凝土與橡膠材料構造層狀周期結構，材料常數見表1。混凝土層與橡膠層厚度均取0.2 m。周期結構的色散曲線見圖2。在0～35 Hz范圍內存在兩個完整帶隙6.565～15.005 Hz與17.81～30.10 Hz。由于工程結構多為0～50 Hz的低頻振動，因此上述帶隙可用于工程結構的隔震、減振。

本文算例的衰減因子隨頻率變化關系見圖3。在禁帶范圍內，衰減因子關于帶隙中心頻率大致對稱，第一濾波帶中心衰減因子為3.54，帶隙邊界對應的衰減因子為零。

表1 材料參數Tab.1 Material coefficients

2 層狀周期結構衰減域參數控制

工程結構振動具有低頻寬帶特性，周期結構衰減域(禁帶)應起始頻率低、衰減域寬。尤其第一衰減域，起始頻率越低且衰減域越寬，對結構隔震、減振越有利。衰減域應徹底濾掉某些頻率范圍的波，故希望衰減因子盡可能大。第一衰減域起始頻率(FLBF)、衰減域寬度(FWAZ)及衰減域衰減因子(β)構成衰減域三個控制參數。以下分別討論周期結構的物理參數及幾何參數對第一衰減域控制參數的影響。

圖2 色散曲線Fig.2 Dispersion curves

圖3 衰減因子Fig.3 Attenuation parameters

圖4 混凝土層厚度對第一衰減域的影響Fig.4 Effects of the thickness of the concrete layer on the first AZ

2.1 第一衰減域起始頻率及寬度

欲獲得較低的第一衰減域起始頻率，即使色散曲線初始斜率盡可能小。將原點處切線方程ω=Ck代入式(4)，并化簡得:

式(5)對h1或h2求導得:

可見，當材料的物理屬性不變時，隨著單元組分材料厚度增大，原點處相速度減小，即第一衰減域的起始頻率降低。混凝土層厚度變化時，FLBF與FWAZ的變化趨勢見圖4。FLBF隨混凝土層厚度的增加而減小，但FWAZ隨著混凝土層厚度的增加會增加。

對物理因素ζ與γ，兩者耦合于材料中。ζ為兩種材料波阻抗差異，ζ=μγ為某一材料的波阻抗。將ζ與γ視為兩獨立參數進行分析:

例如，必修2有關艾弗里證明DNA是遺傳物質的實驗中，教材所展示的是分別將分離的蛋白質、多糖、DNA等加入R型菌的方法（做加法），而新版教材將改變為更接近艾弗里真實實驗的逐步排除的方法（做減法）。教師應當補充這些內容給學生。

隨ζ的增大，原點處相速度減小。FLBF與FWAZ隨ζ的變化見圖5。FLBF隨ζ的增加而減小;FWAZ隨ζ的增加而增加。隨剪切波速的減小，γ增大，原點處斜率減小，第一衰減域起始頻率降低。FLBF與寬度FWAZ隨橡膠層材料的剪切波速單獨變化關系見圖6。可看出，隨橡膠材料剪切波速的增加，FWAZ增加急劇;而FLBF則增加緩慢。

圖5 第一衰減域隨材料波阻抗差系數的變化Fig.5 Effects of the wave impedance coefficient on the first AZ

圖6 第一衰減域隨橡膠材料剪切波波速的變化Fig.6 Effects of the shear wave velocities of the rubber layer on the first AZ

圖7 幾何參數對最大衰減因子的影響Fig.7 Effects of geometrical parameters on the max attenuation parameter

2.2 第一衰減域衰減因子

衰減域的衰減特性由衰減因子控制。利用周期性結構阻隔波動傳播時，通常希望某些頻段波通過衰減域時能被充分過濾，即希望有大的衰減因子。為此，分析各種物理參數及幾何參數對衰減因子影響。材料厚度及波阻抗比對最大衰減因子影響見圖7。由圖7看出，其它物理參數不變時，隨混凝土層厚度增大，衰減因子總在減小;隨橡膠層厚度增大，衰減因子先增大后減小，變化幅度不大。衰減因子隨材料波阻抗差系數的增大而增大。

3 數值模擬與實驗驗證

將無限周期結構頻率帶隙特性應用于工程結構隔震、減振，需研究有限周期結構動力特性。應用Ansys有限元分析軟件，對有限層狀周期結構進行分析。采用Soild-45實體單元建模，平面尺寸為8 m×8 m。同時通過實際地震動實驗，證明層狀周期結構對地震動具有的良好阻隔效果。

3.1 諧響應分析

有限層狀周期結構，周期數取9，混凝土層與橡膠層厚度取0.2 m。約束模型底部y，z方向位移，沿x方向施加幅值為1、頻率0～35 Hz以內的簡諧位移荷載，頻率間隔0.1 Hz，計算不同頻率下有限周期結構動力穩態響應。

經2、4個周期單元濾波后的位移響應函數見圖8。由圖看出，有限周期結構存在明顯的衰減域，兩個衰減域分別與第一帶隙及第二帶隙對應。ω=10.785 Hz時層狀周期結構的瞬態解與穩態解比較見圖9。其中，瞬態解為周期結構底部輸入單一頻率正弦激勵且t∈［0，5 s］時，周期結構各層響應幅值。由圖9看出，有限周期結構中振動幅值沿z軸變化的數值解與無限周期結構中振動幅值沿z軸變化的理論解完全吻合。在不計入瞬態效應情況下，衰減域中的振動主要集中在前1～2個周期范圍內，經2個周期后振動幅值減小到5%左右;計入瞬態效應后衰減幅度略有減小，在頂層自由邊界處存在鞭梢放大效應，但經3個周期后振動仍可衰減到20%左右，鞭梢效應影響并不大。

3.2 組合頻率時程分析

由于衰減域內振動傳播主要集中在前2～3個周期內，瞬態效應及自由邊界影響有限。為此分析含4個周期數的層狀結構在組合頻率下第一衰減域內結構反應。設歸一化組合頻率位移時程為:

其中:ui為(0，1)之間單頻率幅值隨機數，φi為(0.2π)之間相位隨機數，umax為組合位移時程最大值。約束底部y、z方向位移而沿x方向輸入組合位移時程，時間步長Δt=0.005 s。0～5 s內頂層位移響應見圖10。可以看出，外部激勵經四個周期傳播后，減小80%以上。

圖8 位移響應函數Fig.8 Displacement response function

圖9 瞬態與穩態衰減對比Fig.9 Attenuation comparison

圖10 組合頻率位移時程對比Fig.10 Comparison of the displacement response

圖11 振動臺實驗Fig.11 Shaking table test

3.3 實驗驗證及數值分析

為進一步探明層狀周期結構基礎減振、隔震特性，進行兩方面實驗:① 頻域測試，了解結構在不同激勵頻率下的響應特性;② 輸入真實地震動，測試其中包含的多種頻率成份，觀察地震動主頻段落在衰減域內時周期基礎的隔震濾波效果。本試驗取后者。

實驗在臺灣地震研究中心完成。實驗中制作的兩縮尺三層鋼框架結構完全相同，其一直接固定在振動臺上，另一個固定在周期性基礎上，見圖11。鋼框架柱0.06 m ×0.06 m;梁截面 0.03 m ×0.06 m;梁、柱截面鋼板厚0.0032 m。樓板厚0.003 m。框架平面0.6 m×0.6 m，層高0.3 m。每層配重20 kg。周期性基礎含三層混凝土及兩層橡膠，層厚均0.2 m，平面1.0 m×1.0 m。

由于模型在平面內對稱，實驗時振動臺僅在水平向一個方向輸入1975 Oroville地震加速度記錄［13］。該地震加速度記錄的傅里葉譜見圖12，從圖中看出地震動主頻段處于衰減域范圍內。有無隔震基礎情況下上部結構頂層水平加速度響應實驗結果見圖13。從圖中可見有層狀周期基礎情況，上部結構加速度響應減小近50%。本文實驗模型，應用ANSYS進行數值模擬，結構阻尼比取5%。有周期基礎情況的上部結構頂層加速度響應數值模擬結果與實驗測試結果見圖14，由圖可見兩者符合較好。因此，工程中利用周期性結構具有的衰減域特性進行隔震、減振可行。

圖12 地震加速度記錄的傅里葉譜Fig.12 Fourier spectrum of the earthquake acceleration record

圖13 頂層加速度響應Fig.13 Acceleration response at top floor

圖14 結果對比Fig.14 Results comparisons

4 結論

本文給出層狀周期結構衰減域邊界控制方程，并分析影響衰減域的三個特性參數，對層狀周期結構動力衰減域特性進行數值模擬與實驗驗證。結論如下:

(1)材料物理屬性不變時，第一衰減域起始頻率(FLBF)與寬度(FWAZ)隨混凝土層厚度的增加而減小;FLBF隨橡膠層厚度的增加而減小，但FWAZ隨橡膠層厚度的增加而增加，有利于工程應用。

(2)剪切波速一定情況下，隨材料波阻抗差別的增大FLBF會減小，而FWAZ與衰減因子會增大。隨材料剪切波速的增加，FWAZ增加;FLBF隨橡膠層材料剪切波速的增加而緩慢增加。工程中應據具體衰減域要求確定層狀結構的各種參數。

(3)頻率在衰減域范圍內的外部激勵，在層狀周期性結構中傳播主要集中在前1～2個周期。經2個周期后振動可衰減到5%左右，計入瞬態效應后，響應衰減有所降低，但仍達80%以上。

(4)彈性波或振動在層狀周期結構中傳播類似于半波反射，層狀周期基礎可有效降低地震動作用下上部結構響應，周期結構在工程結構的隔震、減振中具有廣闊的應用前景。

［1］ Kushwaha M S，Halevi P，Dobrzynski L，et al.Acoustic band structure of periodic elastic composites［J］.Physical Review Letters，1993，71(13):2022-2025.

［2］Liu Z Y，Zhang X X，Mao Y，et al.Locally resonant sonic materials［J］.Science，2000，289:1734-1736.

［3］溫熙森.聲子晶體［M］.北京:國防工業出版社，2009.

［4］ Xiang H J，Shi Z F.Analysis of flexural vibration band gaps in periodic beams using differential quadrature method［J］.Computers and Structures，2009，87(23-24):1559-1566.

［5］ Ye L P，Lu X Z，Qu Z，et al.Analysis on building seismic damage in the Wenchuan earthquake［C］.The 14thWorld Conference on Earthquake Engineering， Beijing，China，2008.

［6］李宏男，霍林生.結構多維減震控制［M］.北京:科學出版社，2008.

［7］Jia G F，Shi Z F.A new seismic isolation method and its feasibility study［J］. Earthquake Engineering and Engineering Vibration，2010，9(1):75-82.

［8］ Xiang H J，Shi Z F，Bao J.Seismic isolation of buildings with a new type of periodic foundations［C］.ASCE Earth and Space 12th Conference，Honolulu，Hawaii，2010.

［9］ Shi Z F，Cheng Z B，Xiong C.A new seismic isolation method by using a periodic foundation［C］.ASCE Earth and Space 12th Conference，Honolulu，Hawaii，2010.

［10］溫激鴻.聲子晶體振動帶隙及減振特性研究［D］.長沙:國防科技大學，2005.

［11］ Delph T J，Herrmann G，Kaul R K.Harmonic wave propagation in a periodically layered，infinite elastic body:antiplane strain［J］.Journal of Applied Mechanics，1978，45(2):343.

［12］ Mead D J.Wave propagation in continuous periodic structures:Research contributions from Southampton［J］.Journal of Sound and Vibration，1996，190(3):495-524.

［13］ Peer.Peergroud motion database ［C］. http://peer.berkeley.edu/peer_ground_motion_database，2011.