鄒立華，郭 潤，黃 凱，許志旭

(福州大學(xué) 土木工程學(xué)院，福州 350108)

當(dāng)不同建筑由于建筑用地或使用功能要求而相鄰布置,或同一建筑由于變形及構(gòu)造要求，用變形縫將其分成若干單體時，形成相鄰結(jié)構(gòu)。其中，當(dāng)相鄰結(jié)構(gòu)部分或全部單體為隔震結(jié)構(gòu)時，形成相鄰隔震結(jié)構(gòu)。相鄰結(jié)構(gòu)各單體間的間距通常很小，在強(qiáng)震作用下，相鄰結(jié)構(gòu)單體，特別是相鄰隔震結(jié)構(gòu)，常常由于水平位移過大而發(fā)生碰撞。相鄰結(jié)構(gòu)的碰撞可能導(dǎo)致建筑物破壞甚至倒塌，如1985年的墨西哥地震中，在遭受嚴(yán)重破壞或倒塌的300多棟建筑中，40%以上有碰撞現(xiàn)象，其中15%因?yàn)榕鲎捕顾鶾1-2]。

墨西哥地震以后，強(qiáng)震作用下相鄰結(jié)構(gòu)，特別是相鄰隔震結(jié)構(gòu)的碰撞破壞引起了學(xué)者的關(guān)注并開始進(jìn)行系統(tǒng)的研究。Nagarajaiah等[3]在分析Northridge地震中相鄰結(jié)構(gòu)的碰撞后發(fā)現(xiàn)，具有隔震單體的相鄰結(jié)構(gòu)更容易發(fā)生碰撞。為此，Matsagar等[4]研究了不同類型的隔震結(jié)構(gòu)在地震作用下與相鄰結(jié)構(gòu)的碰撞反應(yīng),結(jié)果表明：隔震結(jié)構(gòu)的隔震效果由于其與相鄰結(jié)構(gòu)碰撞而受到較大影響。Panayiotis等[5]對隔震結(jié)構(gòu)在不同地震波作用下與其具有固接基礎(chǔ)的相鄰建筑的碰撞進(jìn)行研究后認(rèn)為：隔震結(jié)構(gòu)與相鄰結(jié)構(gòu)碰撞響應(yīng)主要取決于結(jié)構(gòu)的動力特性及地震波的特性。

由于通常的橡膠隔震支座水平限位能力不足，結(jié)構(gòu)在強(qiáng)震作用下產(chǎn)生過大位移而更容易與相鄰結(jié)構(gòu)產(chǎn)生碰撞。為此，鄒立華等[6-7]在普通疊層橡膠隔震支座的基礎(chǔ)上，在支座上部豎向荷載施加前，用柔性拉索預(yù)先對支座施加豎向預(yù)壓應(yīng)力，形成預(yù)應(yīng)力隔震支座。預(yù)應(yīng)力隔震支座為變剛度支座，除具有可靠抗拉及減小結(jié)構(gòu)豎向變形作用外，還具有良好的水平限位作用。

本文應(yīng)用非線性Hertz-damp碰撞單元，研究預(yù)應(yīng)力隔震支座對相鄰建筑碰撞的影響及碰撞響應(yīng)的影響因素。

1 分析模型

圖1 相鄰隔震結(jié)構(gòu)碰撞分析模型

如圖1示為由A、B兩等高單體組成的相鄰結(jié)構(gòu)，單體A為一般多層(nA層)建筑，單體B為帶預(yù)應(yīng)力隔震支座的多層(nB層)隔震結(jié)構(gòu)。在強(qiáng)震作用下其碰撞分析模型由結(jié)構(gòu)模型、接觸點(diǎn)的碰撞模型及預(yù)應(yīng)力隔震支座模型組成。

1.1 結(jié)構(gòu)模型

將單體A、B均簡化為多自由度剪切模型，其任意層質(zhì)量、阻尼和剛度分別為mAi、cAi、kAi和mBi、cBi、kBi。為研究方便，假定單體A與地面固接，單體B為隔震結(jié)構(gòu)。兩單體間距為gp，假定碰撞發(fā)生在各樓層質(zhì)點(diǎn)處且為對心碰撞。

1.2 碰撞模型

常用的碰撞模型有線彈簧模型、Kelvin模型[8]、Hertz模型[9]和Hertz-damp模型[10-12]。

圖2 Hertz-damp模型

其中，Hertz-damp模型[13]被認(rèn)為是目前最合理的碰撞模型。它由一個非線性彈簧和一個非線性阻尼并聯(lián)而成(如圖2所示)，能較好地模擬碰撞時碰撞力的非線性變化及能量耗散。質(zhì)點(diǎn)碰撞時，其可表示為

(1)

(2)

ch=ξ(ui-uj-gp)3/2

(3)

式中:阻尼常數(shù)ξ可表達(dá)為

(4)

1.3 預(yù)應(yīng)力隔震支座模型

圖3為預(yù)應(yīng)力橡膠隔震支座示意圖。與普通橡膠隔震支座不同，預(yù)應(yīng)力橡膠隔震支座為水平向變剛度支座。一般情況下，阻尼和豎向剛度變化不大。因此，可以在水平向?qū)⑵浜喕癁榉蔷€性彈簧和線性阻尼系統(tǒng)，且不考慮豎向隔震影響。則在任意時刻其水平剛度可以由下式求出[7-8]：

圖3 預(yù)應(yīng)力橡膠隔震支座

(5)

式中:kb為預(yù)應(yīng)力橡膠隔震支座初始水平剛度，kv為預(yù)應(yīng)力橡膠隔震支座豎向剛度，H為施加預(yù)應(yīng)力前預(yù)應(yīng)力疊層橡膠隔震支座總高度，Lp為預(yù)應(yīng)力筋工作時長度，up為支座水平位移，Δv為豎向壓縮長度。a為修正系數(shù)，可按下式計算：

α=-0.27S1+1.85

(6)

式中，S1為第一形狀系數(shù)。

若初始預(yù)應(yīng)力與上部結(jié)構(gòu)豎向荷載相等時，同參數(shù)的PRB與RB水平剛度隨水平位移變化關(guān)系如圖4示。PRB為水平變剛度支座，其水平剛度值隨水平位移增大而增大，此外，還與預(yù)應(yīng)力值大小、上部荷載大小有關(guān)，詳細(xì)分析見文獻(xiàn)[8]。

圖4 支座剛度隨水平位移變化圖

2 運(yùn)動方程及求解

根據(jù)以上模型，可以得到強(qiáng)震作用下相鄰結(jié)構(gòu)碰撞時各單體的運(yùn)動方程。

單體A：

(7)

單體B：

(8)

聯(lián)立式(7)和式(8)可得：

(9)

式中：

方程組(9)變形后得出如下等效形式：

(10)

式中:

由于式(10)中，C為非經(jīng)典阻尼，不能用振型疊加法求解，故本文采用Newmark-β法求解。

3 算例及分析

某辦公樓A為6層框架結(jié)構(gòu)，與6層的隔震結(jié)構(gòu)B相鄰。結(jié)構(gòu)A層高為4 m，各樓層質(zhì)量為3×104kg，各樓層剛度均為8.0×107N/m，結(jié)構(gòu)阻尼比為0.05；結(jié)構(gòu)B層高4 m，各樓層質(zhì)量為3×104kg，各樓層剛度均為8.0×107N/m，結(jié)構(gòu)阻尼比為0.05。無特別說明時，防震縫寬度為0.03 m。

隔震結(jié)構(gòu)的隔震層分別采用普通橡膠支座和預(yù)應(yīng)力厚層橡膠隔震支座，其中預(yù)應(yīng)力橡膠支座的豎向剛度僅在計算其水平剛度時考慮，隔震層參數(shù)如表1所示。

表1 隔震結(jié)構(gòu)隔震層參數(shù)

為便于比較，PRB的初始水平剛度kb取與RB水平剛度相同的值，支座其它參數(shù)見文獻(xiàn)[8]PRB1。假設(shè)混凝土材料碰撞剛度kh=8.68×103kN/m，支座第一形狀系數(shù)S1=2.35，結(jié)構(gòu)體系模型如圖2所示。

選取EI-Centro地震波，加速度峰值為0.3 g，利用Newmark-β逐步積分法，對結(jié)構(gòu)體系進(jìn)行計算分析。

3.1 碰撞影響分析

為研究碰撞對相鄰隔震結(jié)構(gòu)單體地震響應(yīng)的影響，假設(shè)結(jié)構(gòu)B采用預(yù)應(yīng)力支座，則在地震作用下，分別計算出各單體在碰撞和無碰撞時的位移和加速度響應(yīng)如圖5～圖8所示。由圖5和圖6可以看出，結(jié)構(gòu)A在發(fā)生碰撞時，其位移峰值為0.049 m，是不碰撞時的0.74倍；而加速度峰值為9.94 m/s2，是不碰撞時的1.12倍。由圖7和圖8可以看出，結(jié)構(gòu)B在發(fā)生碰撞時，其位移峰值為0.145 m，比無碰撞時稍??；加速度峰值為9.57 m/s2，是無碰撞時的1.08倍。

這是因?yàn)榻Y(jié)構(gòu)碰撞時，各單體碰撞接觸后相互扶持，從而減小了位移響應(yīng)；但碰撞時的碰撞脈沖使單體的加速度響應(yīng)增大。

圖5 結(jié)構(gòu)A頂層位移時程

圖6 結(jié)構(gòu)A頂層加速度時程

圖7 結(jié)構(gòu)B頂層位移時程

圖8 結(jié)構(gòu)B頂層加速度時程

為研究PRB和RB對相鄰隔震結(jié)構(gòu)碰撞的影響，假設(shè)結(jié)構(gòu)B支座分別采用PRB和RB，分別計算在地震作用下碰撞時各單體的地震響應(yīng)如圖9～圖10所示。

圖9和圖10分別為結(jié)構(gòu)B隔震層位移和加速度碰撞時程。由圖可以看出，PRB隔震層的位移峰值為0.11m，僅為RB隔震層的61%；加速度峰值為9.28m/s2，為RB隔震層隔震結(jié)構(gòu)的0.89倍。

圖9 結(jié)構(gòu)B隔震層碰撞位移時程

圖10 結(jié)構(gòu)B隔震層碰撞加速度時程

圖11 結(jié)構(gòu)A頂層碰撞位移時程

圖12 結(jié)構(gòu)A頂層碰撞加速度時程

圖13 結(jié)構(gòu)B頂層碰撞位移時程

圖14 結(jié)構(gòu)B頂層碰撞加速度時程

圖11和圖12分別為發(fā)生碰撞時結(jié)構(gòu)A頂層位移和加速度時程。由圖可以看出，采用PRB時兩側(cè)的位移峰值分別為0.048 m和0.042 m，分別為采用RB時兩側(cè)位移的1.55倍和0.6倍；采用PRB時的加速度峰值響應(yīng)為9.94m/s2，是采用RB時響應(yīng)峰值的0.88倍。

圖13和圖14分別為發(fā)生碰撞時結(jié)構(gòu)B頂層位移和加速度時程。由圖可以看出，采用PRB時的位移峰值為0.14 m，是采用RB時位移響應(yīng)峰值的0.64倍；加速度峰值為9.57 m/s2，是采用RB時響應(yīng)峰值的的0.87倍。

這主要是因?yàn)?，與RB相比，PRB對結(jié)構(gòu)水平位移有明顯的限位作用，因而在同樣條件下，可以減小單體間碰撞的頻度和脈沖強(qiáng)度，從而減小單體之間的碰撞響應(yīng)。但對于與之相鄰的非隔震結(jié)構(gòu)A，由于采用PRB時結(jié)構(gòu)B的水平位移減小，結(jié)構(gòu)A的容許運(yùn)動空間增大，因而，在碰撞的一側(cè)，結(jié)構(gòu)A的位移響應(yīng)反而增大。

3.2 碰撞影響參數(shù)分析

影響相鄰隔震結(jié)構(gòu)碰撞響應(yīng)的因素較多，本文主要分析相鄰結(jié)構(gòu)的周期比、結(jié)構(gòu)的初始間隙(防震縫寬)、以及預(yù)應(yīng)力隔震支座預(yù)應(yīng)力筋松弛長度對結(jié)構(gòu)碰撞效應(yīng)的影響。

3.2.1 周期比對碰撞力的影響

假設(shè)結(jié)構(gòu)A、B的自振周期分別為TA、TB(隔震結(jié)構(gòu)為隔震后周期，預(yù)應(yīng)力支座按初始剛度)其周期比為TA/TB，其他參數(shù)不變，在采用PRB和RB時最大碰撞力隨周期比的變化如圖15所示。

由圖可以看出，當(dāng)周期比小于0.4時，隨著周期比的增大，碰撞力也增大；當(dāng)周期比大于0.4時，隨著周期比的增大，碰撞力減小。周期比為1時，碰撞力最小，但最小值并不等于0。這主要是因?yàn)榻Y(jié)構(gòu)的碰撞是由于單體間的非同步振動造成的，單體自振周期越接近，則單體越接近于同步振動。本文中單體A和B分別為隔震結(jié)構(gòu)和非隔震結(jié)構(gòu)，即使自振周期相同(周期比為1)，但振型并不相同，因而也沒有發(fā)生完全的同步振動，因而碰撞力并不為0。

由圖還可以看出，采用PRB支座時，其碰撞力比采用RB時要小得多，這說明，采用PRB支座能有效減小相鄰結(jié)構(gòu)的碰撞力。

3.2.2 初始間距對碰撞的影響

由于相鄰結(jié)構(gòu)的碰撞是由于單體間縫寬不能滿足其在地震作用時的自由振動引起的，因而其縫寬對碰撞反應(yīng)有重要影響。假設(shè)其它參數(shù)不變，最大碰撞力隨縫寬變化曲線如圖16所示。

由圖可以看出，無論是采用PRB還是RB支座，最大碰撞力都隨縫寬的增大而減小，與一般結(jié)構(gòu)不同，其中并不存在使碰撞力最小的最佳縫寬。這主要是因?yàn)楦粽鸾Y(jié)構(gòu)的變形主要是平動，因而縫寬越大，碰撞力越小。因此，為避免或減小碰撞現(xiàn)象的發(fā)生，相鄰隔震結(jié)構(gòu)之間必須留有足夠的縫寬。

3.2.3 預(yù)應(yīng)力索松弛長度對碰撞的影響

PRB的水平限位能力主要取決于支座的水平剛度，而水平剛度又與預(yù)應(yīng)力索的松弛長度(也即預(yù)應(yīng)力大小)密切相關(guān)，因而松弛長度可能對相鄰隔震結(jié)構(gòu)的碰撞產(chǎn)生影響。

為研究索松弛長度對碰撞的影響，引入無量綱響應(yīng)因子，規(guī)定以索松弛長度為36 mm時結(jié)構(gòu)的最大碰撞力、最大加速度以及最大位移值分別作為基值，其余松弛長度時計算所得的響應(yīng)值與基值的比值作為響應(yīng)因子，假設(shè)其它參數(shù)不變，響應(yīng)因子隨松弛長度的變化關(guān)系曲線如圖17所示。

圖15 碰撞力隨周期比變化曲線

圖16 碰撞力隨初始間距變化曲線

圖17 響應(yīng)因子隨松弛長度變化曲線

由圖可以看出，總體來說，隨著索松弛長度的增大，最大碰撞力、最大位移以及最大加速度也增大。但當(dāng)松弛長度在30～70 mm時增長較快，70 mm以后增速明顯減緩。這是因?yàn)?，?dāng)索的松弛長度達(dá)到一定值以后，PRB的性能接近于RB，其限位效果逐漸減小。

隨著松弛長度的增大，最大加速度的增長速度相對于最大碰撞力增長較慢，這是因?yàn)閱误w樓層的加速度取決于碰撞力和支座的隔震效果，當(dāng)索的松弛長度達(dá)到一定值以后，PRB的性能接近于RB，其隔震效果有所增加；而碰撞力的大小同時取決于單體的加速度和位移響應(yīng)，隨著松弛長度的增大，雖然加速度有所減小，但位移響應(yīng)則有所增大。

4 結(jié) 論

根據(jù)以上計算與分析,可以得到如下結(jié)論:

(1)與一般相鄰結(jié)構(gòu)相比，相鄰隔震結(jié)構(gòu)在地震作用下更容易發(fā)生碰撞。碰撞通常會使單體結(jié)構(gòu)加速度響應(yīng)增大，而位移響應(yīng)減小。

(2)由于PRB的水平限位作用，與RB相比，采用PRB的相鄰隔震結(jié)構(gòu)能有效減小單體之間的碰撞響應(yīng)。

(3)采用PRB的相鄰隔震結(jié)構(gòu)的碰撞與單體的周期比、單體間縫寬及PRB索的松弛長度等因素有關(guān)。一般來說，周期比越接近，碰撞響應(yīng)越小，但無論周期比多接近，都不會產(chǎn)生完全的同步振動；縫寬越大，碰撞響應(yīng)越小，且不存在最優(yōu)縫寬；預(yù)應(yīng)力索的松弛長度越大，碰撞響應(yīng)越大，但隨著松弛長度的增大，其影響力減小。

[1]Bertero V V.Observations on structural pounding[J].Proceedings International Conference on Mexico Earthquakes，ASCE，1987:264-278.

[2]Rosenblueth E，Meli R.The 1985 earthquake: Ca-uses and effects in Mexico City[J].Concrete International，1986，8: 23-34.

[3]Nagarajaiah S,Sun X H.Base-isolated FCC building: Impact response in Northridge earthquake[J].Journal of Structural Engineering，2001，127 (9):1063-1075.

[4]Matsagar V A，Jangid R S.Seismic response of base-isolated structures during impact with adjacent structures[J].Engineering Structures，2003，25: 1311-1323.

[5]Panayiotis C P，Petros K.On poundings of a seismically isolated building with adjacent structures during strong earthquakes[J].Earthquake Engineering Structural dynamics，2010，39:933-940.

[6]Zou L H，Huang K，Zhang W，et al.Prestress rubber isolator and its mechanical properties[J].Advanced Materials Research，2011，163-167: 4405-4414.

[7]鄒立華,饒宇,黃凱,等.預(yù)應(yīng)力厚層橡膠支座隔震性能研究[J].建筑結(jié)構(gòu)學(xué)報,2012，32(2)：76-82.

ZOU Li-hua，RAO Yu，HUANG Kai，et al.Research on isolating property of prestressed thick rubber bearings[J].Journal of Building Structures,2012，32(2)：76-82.

[8]Anagnostopoulos S A.Pounding of buildings in series during earthquakes[J].Earthquake Engineering and Structural Dynamics 1988，16(3): 443-456.

[9]Davis R O.Pounding of buildings modelled by an impact oscillator[J].Earthquake Engineering and Structural Dynamics，1992，21(3):253-274.

[10]Marhefka D W，Orin D E.A compliant contact model with nonlinear damping for simulation of robotic systems[J].IEEE Transactions on Systems，Man，and Cybernetics-Part A: Systems and Hu-Mans，1999，29: 566-572.

[11]Lankarani H M，Nikravesh P E.A contact force model with hysteresis damping for impact analysis of multibody systems[J].Journal of Mechanical Design (ASME)，1990，112: 369-376.

[12]Jankowski R.Non-linear viscoelastic modelling of earthquake-induced structural pounding[J].Earthquake Engineering and Structural Dynamics，2005，34: 595-611.

[13]Muthukumar S，DesRoches R.A Hertz contact model with non-linear damping for pounding simulation[J].Earthquake Engineering and Struc-tural dynamics，2006，35: 811-828.

[14]Panos T，Marc O E，John F S.Design of seismic restrainers for in-span hinges[J].Journal of Structural Engineering，1997，123(4): 469-478.