馬涌泉　邱洪興

應(yīng)分析新方法.通過在基礎(chǔ)隔震結(jié)構(gòu)的頂部布置調(diào)諧質(zhì)量阻尼器來(lái)構(gòu)建減震結(jié)構(gòu)，分別采用BoucWen模型及其剛度退化模型模擬隔震層及各樓層的滯回特性.通過在精細(xì)積分法中引入復(fù)化Cotes積分，并結(jié)合虛擬激勵(lì)法，提出了求解減震結(jié)構(gòu)非平穩(wěn)隨機(jī)反應(yīng)的CCIM法.依據(jù)首次超越破壞準(zhǔn)則，建立了以結(jié)構(gòu)層間位移角為評(píng)價(jià)指標(biāo)的動(dòng)力可靠度極限狀態(tài)方程.通過分別采用CCIM、蒙特卡羅法和時(shí)域顯式蒙特卡羅法對(duì)減震結(jié)構(gòu)進(jìn)行隨機(jī)反應(yīng)分析，驗(yàn)證了CCIM具有高效率和高精度的特點(diǎn).以一座30層鋼框架結(jié)構(gòu)為算例，分別計(jì)算了減震結(jié)構(gòu)、基礎(chǔ)隔震結(jié)構(gòu)和未控制結(jié)構(gòu)在8度和9度罕遇地震作用下的隨機(jī)反應(yīng).結(jié)果表明：本文提出的“頂吸基隔”減震結(jié)構(gòu)的整體可靠度比基礎(chǔ)隔震結(jié)構(gòu)和未隔震結(jié)構(gòu)的都要高，該減震結(jié)構(gòu)具有極大的工程推廣價(jià)值.

關(guān)鍵詞：建筑物；抗震設(shè)計(jì)；地震反應(yīng)；調(diào)諧質(zhì)量阻尼器；精細(xì)積分法；動(dòng)力可靠度

中圖分類號(hào)：O328； TU352.1 文獻(xiàn)標(biāo)識(shí)碼：A

隨著中國(guó)城鎮(zhèn)化建設(shè)步伐的加快，各地興建了為數(shù)不少的高層建筑，如何提升它們的抗震性能，一直是學(xué)者們研究的重點(diǎn).基礎(chǔ)隔震作為一項(xiàng)有效的減震技術(shù)，在高層結(jié)構(gòu)中正得到普遍應(yīng)用，但其隔震層在震后會(huì)產(chǎn)生較大變形[1].雖然可以通過在隔震層中增設(shè)阻尼器或限位器來(lái)減小或限制其位移，但會(huì)引起主結(jié)構(gòu)層間位移的增大[2].在基礎(chǔ)隔震結(jié)構(gòu)中引入調(diào)諧質(zhì)量阻尼器（TMD）組成“TMD基礎(chǔ)隔震”結(jié)構(gòu)可以解決隔震層及主結(jié)構(gòu)位移均較大的難題，本文將這種結(jié)構(gòu)命名為“頂吸基隔”減震結(jié)構(gòu)（簡(jiǎn)稱減震結(jié)構(gòu)）.學(xué)者們對(duì)該結(jié)構(gòu)的減震效果進(jìn)行過一些研究[3-4]，但這些研究均是輸入確定的地震動(dòng)來(lái)求解結(jié)構(gòu)反應(yīng)，由于地震動(dòng)具有隨機(jī)性，因此這些研究成果缺乏參考價(jià)值；目前學(xué)者們僅開展了基礎(chǔ)隔震結(jié)構(gòu)隨機(jī)反應(yīng)的求解工作[5-6].然而他們所采用的精細(xì)積分法在計(jì)算過程中會(huì)存在矩陣求逆的問題，不僅計(jì)算量大，而且穩(wěn)定性也較差，甚至?xí)霈F(xiàn)逆矩陣不存在的情況.因此探尋減震結(jié)構(gòu)非平穩(wěn)隨機(jī)反應(yīng)的高效率和高精度的求解方法就顯得尤為重要.

本文采用BoucWen模型及其剛度退化模型描述隔震層及樓層的滯回特性，利用精細(xì)積分法和虛擬激勵(lì)法推導(dǎo)出復(fù)化Cotes精細(xì)積分法（CCIM），并驗(yàn)證了其精度和效率.使用CCIM分別對(duì)一座30層減震、基礎(chǔ)隔震和未隔震結(jié)構(gòu)進(jìn)行非平穩(wěn)隨機(jī)反應(yīng)與動(dòng)力可靠度分析，得出的結(jié)論可供類似結(jié)構(gòu)減震設(shè)計(jì)時(shí)參考.

經(jīng)過精細(xì)計(jì)算以后，可以使得其值非常接近精確值，但是其積分項(xiàng)解析形式的精度仍不能保證.文獻(xiàn)[10]將精細(xì)積分法引入隨機(jī)地震反應(yīng)分析，提出了時(shí)域顯式蒙特卡羅法（PTIM），它假定地震激勵(lì)在（tl，tl+1）內(nèi)線性變化，但當(dāng)?shù)卣鸺?lì)具有強(qiáng)非線性時(shí)這種算法會(huì)帶來(lái)較大誤差，如果通過減小時(shí)間步長(zhǎng)Δt來(lái)提高計(jì)算精度，勢(shì)必會(huì)增加計(jì)算量.文獻(xiàn)[11]在結(jié)構(gòu)隨機(jī)反應(yīng)分析中引入了Simpson積分，雖然提高了結(jié)構(gòu)反應(yīng)的精度，但由于該積分方法的固有缺陷，只能達(dá)到三次代數(shù)精度.為了解決上述問題，本文引入復(fù)化Cotes積分理論來(lái)求解結(jié)構(gòu)隨機(jī)反應(yīng)一般解的積分項(xiàng)，提出了一套基于復(fù)化Cotes理論的精細(xì)積分法，并將其命名為CCIM.該方法既能良好地處理非線性地震激勵(lì)，又能獲得比Simpson積分更高的精度，并且對(duì)Δt也不敏感CCIM的建立過程如下：將復(fù)化Cotes積分引入式（30），積分區(qū)間分點(diǎn)為

上式即為求解結(jié)構(gòu)反應(yīng)的顯式表達(dá)式.結(jié)合一階矩和二階矩的運(yùn)算特點(diǎn)，可得Z（tl）的期望及方差分別為：

式（40）中的協(xié)方差矩陣可由隨機(jī)地震激勵(lì)的相關(guān)函數(shù)構(gòu)成，相關(guān)系數(shù)可由式（18）求得.

相鄰時(shí)刻的系數(shù)矩陣存在以下遞推關(guān)系：

綜上所述，按照式（38）構(gòu)造出該時(shí)刻的ρl和Yl便可求得tl時(shí)刻的結(jié)構(gòu)反應(yīng)Z（tl）；如果求解（tl，tl+1）時(shí)間段內(nèi)的Z（tl），只需結(jié)合式（41）便可較快地獲得各時(shí)刻點(diǎn)的系數(shù)矩陣.

上述CCIM法的推導(dǎo)雖然是以剪切型結(jié)構(gòu)為例，但是此算法是以精細(xì)積分法和復(fù)化Cotes積分理論為基本框架，因此它對(duì)結(jié)構(gòu)的動(dòng)力特性矩陣沒有任何特殊要求，故此算法完全適用于求解彎剪型或彎曲型結(jié)構(gòu)的非平穩(wěn)隨機(jī)地震反應(yīng).

2.3結(jié)構(gòu)層間位移角可靠度分析

由CCIM得到隔震層和各樓層的位移方差σ2x和速度方差σ2，便可求得其層間位移角方差σ2θ，從而依據(jù)首次超越破壞準(zhǔn)則，求得其層間位移角的動(dòng)力可靠度.對(duì)于給定的層間位移角限值θb，隔震層和各樓層的層間位移角θ（t）在時(shí)間Λ內(nèi)超越限值θb的次數(shù)的期望可表示為：

3CCIM精度與效率的驗(yàn)證

為了驗(yàn)證CCIM的計(jì)算精度及效率，本文選用一座10層鋼框架結(jié)構(gòu)為算例.該算例位于Ⅱ類場(chǎng)地，抗震設(shè)防烈度為8度，設(shè)計(jì)基本地震加速度值為0.30g，處于設(shè)計(jì)地震第1組.首先對(duì)其進(jìn)行8度罕遇地震作用下的基礎(chǔ)隔震設(shè)計(jì)，獲得隔震層的各項(xiàng)最優(yōu)參數(shù)；然后對(duì)基礎(chǔ)隔震結(jié)構(gòu)進(jìn)行模態(tài)分析獲得其前10階自振頻率，以此獲得TMD系統(tǒng)的各項(xiàng)最優(yōu)參數(shù)（保證TMD系統(tǒng)的前10階自振頻率與基礎(chǔ)隔震結(jié)構(gòu)的分別相等或十分接近）；最后對(duì)減震結(jié)構(gòu)進(jìn)行模態(tài)分析，提取其前三階阻尼比，并取它們的平均值作為該結(jié)構(gòu)的阻尼比.分別運(yùn)用蒙特卡羅法[12]（MC）、CCIM和PTIM對(duì)減震結(jié)構(gòu)進(jìn)行8度罕遇地震作用下的非平穩(wěn)隨機(jī)地震反應(yīng)分析.

定義8和9度罕遇地震下的減震結(jié)構(gòu)分別為工況Ⅰ和Ⅱ，8和9度罕遇地震下的基礎(chǔ)隔震結(jié)構(gòu)分別為Ⅲ和Ⅳ，8和9度罕遇地震下的未隔震結(jié)構(gòu)分別為Ⅴ和Ⅵ.采用CCIM分別對(duì)6種工況進(jìn)行非平穩(wěn)隨機(jī)分析，得到頂層位移、第6層層間位移角及頂層加速度等方差時(shí)程如圖3所示.得到隔震層和各樓層的位移、層間位移角及加速度等方差峰值及層間位移角可靠度如圖4～圖5所示.

4.2結(jié)果分析

由圖3可知，在整個(gè)時(shí)程內(nèi)，減震結(jié)構(gòu)的位移方差、層間位移角方差和加速度方差比基礎(chǔ)隔震和未隔震結(jié)構(gòu)的方差都要小.6種工況的反應(yīng)方差均能在較短的時(shí)間內(nèi)收斂到一個(gè)穩(wěn)定的解，說明本文提出的CCIM具有良好的穩(wěn)定性，此方法可用于評(píng)估高層或多層結(jié)構(gòu)在非平穩(wěn)隨機(jī)地震作用下的抗震性能.

由圖4可知，減震結(jié)構(gòu)的隔震層和各樓層位移方差峰值、層間位移角方差峰值和加速度方差峰值均比基礎(chǔ)隔震和未隔震結(jié)構(gòu)（無(wú)隔震層）的要??；6種工況的上述三項(xiàng)值在第6層處均存在不同程度的轉(zhuǎn)折，其中減震結(jié)構(gòu)的轉(zhuǎn)折程度最不明顯.這是由于原結(jié)構(gòu)的第6層剛度有突變，形成了薄弱層.可見，通過在基礎(chǔ)隔震結(jié)構(gòu)頂部布置TMD而形成的減震結(jié)構(gòu)可以有效解決隔震層及主結(jié)構(gòu)在震后出現(xiàn)過大位移的難題.

由圖5可知，減震結(jié)構(gòu)無(wú)論是在8度還是9度罕遇地震下，其隔震層和各樓層的層間位移角可靠度均為100%；8度和9度罕遇地震下的基礎(chǔ)隔震結(jié)構(gòu)其隔震層和各樓層的層間位移角可靠度也均較高，根據(jù)式（46）算得其整體可靠度分別為95.53%和91.76%，但第6層處的可靠度有一定的降低，可見結(jié)構(gòu)剛度的突變對(duì)基礎(chǔ)隔震結(jié)構(gòu)的層間位移角可靠度有一定的影響；8度和9度罕遇地震下的未隔震結(jié)構(gòu)其各樓層的層間位移角可靠度很低，在第6層處達(dá)到了最低，根據(jù)式（45）算得其整體可靠度分別為58.82%和46.91%，因此在9度罕遇地震下，未隔震結(jié)構(gòu)的大多數(shù)層間位移角已超過了彈性層間位移角限值.可見在罕遇地震作用下，減震結(jié)構(gòu)的整體可靠度比基礎(chǔ)隔震和未隔震結(jié)構(gòu)的都要高得多.

5結(jié)論

1） 本文提出的CCIM其計(jì)算精度和計(jì)算效率要比PTIM的都要高，其計(jì)算結(jié)果逼近MC的結(jié)果，而CCIM所需的計(jì)算時(shí)間卻比MC的要少得多；CCIM能使結(jié)構(gòu)反應(yīng)方差在較短時(shí)間內(nèi)收斂到穩(wěn)定解.可見，本文提出的CCIM兼?zhèn)涓咝屎透呔鹊奶攸c(diǎn)，可望有極佳的工程應(yīng)用前景.

2） 減震結(jié)構(gòu)的樓層及隔震層位移、層間位移角和加速度方差比基礎(chǔ)隔震結(jié)構(gòu)的都要小，前者的整體可靠度比基礎(chǔ)隔震和未隔震結(jié)構(gòu)的都要高；薄弱層的存在對(duì)減震結(jié)構(gòu)的整體可靠度幾乎沒有影響.可見，本文提出的“頂吸基隔”減震結(jié)構(gòu)具有卓越的抗震性能，可望有良好的工程推廣價(jià)值.

參考文獻(xiàn)

[1]ARIGA T， KANNO Y， TAKEWAKI I. Resonant behaviour of baseisolated highrise buildings under longperiod ground motions[J]. The Structural Design of Tall and Special Buildings， 2006， 15（3）： 325-338.

[2]SPENCER B F， NAGARAJAIAH S. State of the art of structural control[J]. Journal of Structural Engineering， 2003， 129（7）： 845-856.

[3]TANIGUCHI T， KIUREGHIAN A D， MELKUMYAN M. Effect of tuned mass damper on displacement demand of baseisolated structures[J]. Engineering Structures， 2008， 30（12）： 3478-3488.

[4]DE I M， PETTI L， PALAZZO B. Combined control strategy base isolation and tuned mass damper： an effectiveness analysis of its application to nonlinear benchmark base isolated structure[C]//Proceedings of the 14th World Conference on Earthquake Engineering. Beijing： Science Press， 2008： 129-138.

[5]馬長(zhǎng)飛， 張亞輝， 譚平， 等. 三維非偏心基礎(chǔ)隔震結(jié)構(gòu)非平穩(wěn)隨機(jī)振動(dòng)分析[J]. 工程力學(xué)， 2013， 30（4）： 198-203.

MA Changfei， ZHANG Yahui， TAN Ping， et al. Nonstationary stochastic analysis for 3D noneccentric baseisolated buildings[J]. Engineering Mechanics， 2013， 30（4）： 198-203. （In Chinese）

[6]JANGID R S. Stochastic response of building frames isolated by leadrubber bearings[J]. Structural Control and Health Monitoring， 2010， 17（1）： 1-22.

[7]ERLICHER S， BURSI O S. BoucWenType models with stiffness degradation： thermodynamic analysis and applications[J]. Journal of Engineering Mechanics， 2008， 134（10）： 843-855.

[8]LIN Jiahao， ZHANG Wenshou， WILLIAMS F W. Pseudoexcitation algorithm for nonstationary random seismic responses[J]. Engineering Structures， 1994， 16（4）： 270-276.

[9]ZHONG Wanxie. On precise integration method[J]. Journal of Computational and Applied Mathematics， 2004， 163（1）： 59-78.

[10]ZHANG Yahui， LI Qiusheng， LIN Jiahao， et al. Random vibration analysis of longspan structures subjected to spatially varying ground motions[J]. Soil Dynamics and Earthquake Engineering， 2009， 29（4）： 620-629.

[11]POULIMENOS A G， FASSOIS S D. Parametric timedomain methods for nonstationary random vibration modelling and analysisa critical survey and comparison[J]. Mechanical Systems and Signal Processing， 2006， 20（4）： 763-816.

[12]PROPPE C， PRADLWARTER H J， SCHUELLER G I. Equivalent linearization and Monte Carlo simulation in stochastic dynamics[J]. Probabilistic Engineering Mechanics， 2003， 18（1）： 1-15.

[13]SEGHIR M. Generalized minimum variance gain scheduling controller for nonlinear structural systems under seismic ground motion[J]. Journal of Electrical Systems， 2012， 8（4）： 472-480.

[14]GB 50011-2010建筑抗震設(shè)計(jì)規(guī)范[S]. 北京： 中國(guó)建筑工業(yè)出版社， 2010： 172-176.

GB 50011-2010Code for seismic design of buildings[S]. Beijing： China Architecture and Building Press， 2010： 172-176. （In Chinese）

[15]薛素鐸， 王雪生， 曹資. 基于新抗震規(guī)范的地震動(dòng)隨機(jī)模型參數(shù)研究[J]. 土木工程學(xué)報(bào)， 2003， 36（5）： 5-10.