考慮主余震序列影響的低延性鋼筋混凝土框架易損性分析

韓建平，李 軍

(1.蘭州理工大學(xué)甘肅省土木工程防災(zāi)減災(zāi)重點實驗室，甘肅，蘭州 730050；2.蘭州理工大學(xué)防震減災(zāi)研究所，甘肅，蘭州 730050)

大量震害研究表明，主震發(fā)生后伴隨著大量的余震作用，這些相繼發(fā)生的余震會對主震受損結(jié)構(gòu)造成進一步的破壞，甚至導(dǎo)致主震受損結(jié)構(gòu)的倒塌，特別是低延性結(jié)構(gòu)[1]。我國目前有大量老舊的低延性鋼筋混凝土結(jié)構(gòu)，這類結(jié)構(gòu)沒有經(jīng)過抗震設(shè)計，抗震設(shè)防水平不足或抗震設(shè)計不當(dāng)、抗震構(gòu)造措施欠缺，加之這類建筑已經(jīng)服役多年，使得這類結(jié)構(gòu)成為抗震設(shè)防的薄弱環(huán)節(jié)。研究表明[2]，此類既有低延性結(jié)構(gòu)配箍率低、結(jié)構(gòu)耗能能力弱，極易發(fā)生脆性破壞，因此該類結(jié)構(gòu)是抗震性能研究的重點。較為精確地模擬這類低延性結(jié)構(gòu)的破壞機制，評估其在主余震序列作用下的抗震性能有利于為該類結(jié)構(gòu)的加固與改造提供理論依據(jù)。

目前已有學(xué)者對主余震作用下結(jié)構(gòu)的抗震性能進行了分析研究，針對低延性鋼筋混凝土框架結(jié)構(gòu)或構(gòu)件進行了試驗和數(shù)值模擬分析。針對箍筋不足的低延性鋼筋混凝土柱，Elwood[3]提出極限狀態(tài)材料模型，通過定義破壞面來模擬低延性框架柱剪切及軸壓破壞；Shoraka等[4]對Elwood提出的極限狀態(tài)材料模型通過軸-剪-彎耦合作用(axial-shearflexure interaction, ASFI)分析對柱的剪切破壞點進行了理論性的界定，對低延性柱的彎剪耦合破壞取得了較好的模擬效果，Ghannoum和Moehle[5]則以框架柱塑性鉸區(qū)轉(zhuǎn)角為指標(biāo)，定義了鋼筋混凝土柱彎剪破壞的破壞準(zhǔn)則，對框架柱倒塌階段顯著的強度和剛度退化行為進行了較為準(zhǔn)確的模擬。在梁柱節(jié)點的研究中，Lowes等[6]提出了考慮節(jié)點核心區(qū)剪切破壞和節(jié)點梁筋黏結(jié)滑移的節(jié)點宏觀單元，并應(yīng)用修正壓力場理論(modified compression field theory，MCFT)定義節(jié)點區(qū)剪切域的行為，Mitra等[7]分析了該模型，論證了低延性框架節(jié)點區(qū)剪切主要由斜壓桿機制傳遞。在整體框架結(jié)構(gòu)的研究中，Yavari等[8]對沒有進行抗震設(shè)計結(jié)構(gòu)的振動臺試驗研究表明，節(jié)點核心區(qū)和柱的剪切破壞對結(jié)構(gòu)整體變形和失效模式產(chǎn)生顯著影響，周洲等[9]對主余震作用下鋼筋混凝土框架進行了易損性分析和安全評估，于曉輝等[10―11]采用多條真實主余震和人工構(gòu)造主余震序列作為輸入，通過對比單自由度體系在主余震序列和主震單獨作用下的結(jié)構(gòu)損傷，建立了主余震序列地震動的增量損傷譜，分析了不同主余震序列對結(jié)構(gòu)的增量損傷，張沛洲等[12]提出一種分析結(jié)構(gòu)在主余震序列作用下抗震性能的方法，并用該方法分析了余震對鋼筋混凝土框架結(jié)構(gòu)抗震性能的影響。綜上可知，梁柱節(jié)點區(qū)和柱的剪切破壞狀態(tài)和柱的強度和剛度退化直接決定著低延性結(jié)構(gòu)的抗震性能，而傳統(tǒng)的主余震作用下低延性結(jié)構(gòu)抗震性能的分析并沒有考慮梁柱節(jié)點和柱的剪切破壞以及柱強度和剛度的退化行為，使得分析結(jié)果容易低估主余震作用對低延性結(jié)構(gòu)安全潛在的威脅。

針對這一問題，本文以某未進行抗震設(shè)計的低延性鋼筋混凝土框架結(jié)構(gòu)為研究對象，考慮梁柱節(jié)點和框架柱剪切破壞，分別采用真實主余震序列和基于重復(fù)法構(gòu)造的主余震序列，評估該結(jié)構(gòu)在主余震序列作用下的抗震性能，以期為低延性結(jié)構(gòu)抗震性能評估及加固改造提供理論依據(jù)。

1 低延性框架結(jié)構(gòu)梁柱節(jié)點和柱的破壞機制

低延性框架梁柱節(jié)點典型的破壞模式有：① 框架梁或柱屈服前框架梁柱節(jié)點區(qū)發(fā)生剪切破壞，結(jié)構(gòu)表現(xiàn)出強度和剛度的快速退化；② 框架梁屈服后框架梁柱節(jié)點區(qū)發(fā)生剪切破壞，結(jié)構(gòu)表現(xiàn)出強度和剛度的逐漸退化，結(jié)構(gòu)具有中等的延性行為；③ 梁底鋼筋在梁柱節(jié)點區(qū)伸入長度小，導(dǎo)致梁屈服之前發(fā)生錨固破壞。對于低延性框架柱通常發(fā)生彎剪破壞，導(dǎo)致側(cè)向承載能力降低。典型的低延性框架柱的破壞模式有：① 框架柱彎曲屈服之前一個或多個臨界斜裂縫的突然擴展以及柱混凝土被破碎或斜劈裂破壞；② 框架柱彎曲屈服之后，導(dǎo)致柱塑性鉸區(qū)抗剪承載能力降低而發(fā)生柱的彎剪破壞。

依據(jù)低延性鋼筋混凝土框架梁柱節(jié)點和柱的破壞機制，節(jié)點區(qū)的剪切破壞用一個轉(zhuǎn)動彈簧來模擬，低延性柱彎剪破壞以及強度和剛度退化用Hysteretic材料模型。

低延性梁柱節(jié)點區(qū)應(yīng)力-應(yīng)變關(guān)系采用文獻[13]簡化的壓桿-拉桿模型(strut and tie model，STM)計算。根據(jù)文獻[14]，低延性柱的初始剪切強度Vn和剪切剛度退化斜率可分別通過式(1)和式(2)計算，并假設(shè)柱殘余剪切強度為初始剪切強度Vn的20%。

式中：k為強度退化系數(shù)；λ＝0.75或1，分別對應(yīng)輕型骨料混凝土或一般混凝土；fc為混凝土抗壓強度；為彎矩M和剪力V與柱截面有效高度d的乘積之比為作用在柱上的軸向荷載；Ag為柱截面面積；Av為箍筋面積；fyv為箍筋屈服強度；s為箍筋間距。

式中：dc為截面核心尺寸(箍筋截面中心到柱截面中心的距離)；L為柱高。

2 低延性鋼筋混凝土框架結(jié)構(gòu)的設(shè)計與建模

2.1 低延性結(jié)構(gòu)的設(shè)計信息

選取西北地區(qū)未經(jīng)過抗震設(shè)計的6層低延性鋼筋混凝土框架結(jié)構(gòu)為研究對象。不考慮風(fēng)荷載和抗震設(shè)防，鋼筋混凝土現(xiàn)澆樓板，樓板厚100 mm，屋面和樓面荷載標(biāo)準(zhǔn)值分別為 6.0 kN/m2、4.5 kN/m2；屋面和樓面活荷載標(biāo)準(zhǔn)值分別為0.5 kN/m2、2.0 kN/m2；考慮填充墻的作用，標(biāo)準(zhǔn)層框架梁上施加7.5 kN/m的線荷載，考慮女兒墻的作用，屋面梁施加3 kN/m的線荷載；梁柱縱向受力鋼筋均采用熱軋鋼筋 HRB335，其余鋼筋采用HPB235，混凝土采用C30。結(jié)構(gòu)的立面及梁柱配筋如圖1所示。

2.2 OpenSees有限元模型

根據(jù)上述方法精確模擬低延性框架梁柱節(jié)點區(qū)剪切破壞并結(jié)合纖維模型建立 OpenSees[15]有限元分析模型，部分細節(jié)如圖2所示。該模型中，梁、柱采用非線性梁柱單元(nonlinear beam column element)模擬；節(jié)點采用轉(zhuǎn)動彈簧以考慮節(jié)點區(qū)剪切變形的影響，節(jié)點區(qū)轉(zhuǎn)動彈簧采用可以考慮強度退化、剛度退化和捏攏效應(yīng)的三線性 Pinching4材料模型。其 16個骨架曲線參數(shù)設(shè)置參考宋孟超[13]簡化的STM，其余參數(shù)參照Mitra等[7]和Paspuleti[16]的推薦取值。柱端通過零長度單元考慮柱的彎剪破壞以及強度和剛度退化并應(yīng)用Hysteretic材料模擬。梁與節(jié)點連接處縱向鋼筋的粘結(jié)滑移通過梁端設(shè)置零長單元(zero-length element)，采用文獻[17]建立鋼筋粘結(jié)滑移關(guān)系，用Bond_SP01材料來模擬。關(guān)于纖維截面材料模型，鋼筋采用考慮包辛格效應(yīng)的隨動強化模型(Steel02)模擬；混凝土采用能夠考慮剛度退化和箍筋約束作用的 Kent-Scott-Park模型(Concrete01)。

圖1 算例低延性RC框架結(jié)構(gòu)布置Fig.1 Structural configuration of the investigated low-ductile RC frame

圖2 算例低延性RC框架OpenSees數(shù)值模型細節(jié)Fig.2 Modelling details of OpenSees numerical model for the investigated low-ductile RC frame

2.3 低延性梁柱節(jié)點模型驗證

為驗證上述模型對低延性鋼筋混凝土框架節(jié)點模擬的可靠性，選取 Li等[18―19]完成的試件 A1和試件AL1擬靜力試驗結(jié)果，進行滯回性能模擬分析，分析結(jié)果如圖3所示。

圖3 低延性RC框架節(jié)點滯回性能試驗與數(shù)值模擬對比Fig.3 Comparisons of the simulated and tested hysteretic curves for low-ductile RC frame joints

為了定量分析圖2模型在低延性結(jié)構(gòu)數(shù)值分析中的準(zhǔn)確性，根據(jù)圖3分別計算選取的試件A1和試件 AL1每次循環(huán)往復(fù)荷載作用下試驗結(jié)果和數(shù)值分析結(jié)果的累計耗能，計算可得試件A1滯回曲線累積耗能的誤差在6%~23.2%，試件AL1累積耗能的誤差在 1%~19.4%。該數(shù)值分析結(jié)果的初始剛度和極限承載力與試驗結(jié)果吻合較好，且循環(huán)往復(fù)荷載作用下數(shù)值分析結(jié)果的累積耗能與試驗結(jié)果的誤差也較小，因此圖2模型能較準(zhǔn)確地模擬低延性鋼筋混凝土框架節(jié)點由剪切作用引起的強度和剛度退化現(xiàn)象。

3 主余震序列的選取與構(gòu)造

為了全面分析不同主余震序列對低延性結(jié)構(gòu)抗震性能評估的影響，采用真實主余震序列和基于重復(fù)法構(gòu)造的主余震序列作為地震動輸入。基于文獻[20]選取的 65條真實主余震序列信息如表1所示。采用真實主余震序列進行分析時，考慮余震次數(shù)對低延性結(jié)構(gòu)抗震性能的影響。在主余震序列構(gòu)造時，主余震序列的主震記錄和余震記錄之間設(shè)置90秒的時間間隔，其目的是使結(jié)構(gòu)在主震后經(jīng)自由振動達到平衡位置，如圖4所示。采用重復(fù)法構(gòu)造主余震序列[21]進行分析時，考慮余震方向性的影響，如圖5所示。

為了量化余震對低延性結(jié)構(gòu)抗震性能的影響，本文亦對主震單獨作用下結(jié)構(gòu)的抗震性能進行分析。為了方便論述，將所考慮的5種地震作用工況縮寫為：主震(MS)、基于重復(fù)法構(gòu)造的主余震序列(余震為正向)(MA-REP-P)、基于重復(fù)法構(gòu)造的主余震序列(余震為反向)(MA-REP-N)、真實記錄的主震+1條余震記錄(MA-REC-1)、真實的主震+2條余震記錄(MA-REC-2)。

圖4 選取的真實主余震序列示例Fig.4 Example of selected real mainshock-aftershock sequence

圖5 基于重復(fù)法構(gòu)造的主余震序列及余震方向的考慮Fig.5 The mainshock-aftershock sequence constructed based on the repeated method and the consideration of aftershock direction

表1 選取的真實主余震序列信息Table 1 Summary of the selected real mainshock-aftershock sequences

續(xù)表

4 主余震作用下低延性鋼筋混凝土框架易損性分析

4.1 主余震作用下結(jié)構(gòu)易損性分析方法

采用譜加速度Sa作為地震動強度指標(biāo)，根據(jù)文獻[22]，地震易損性函數(shù)可以表示為：

式中：P[D＞C｜IM]表示強度為IM的地震作用下地震需求D大于抗震能力C的超越概率；SP表示地震需求D的中位值；SC表示抗震能力的中位值；βC和βD|IM分別表示抗震能力C和地震需求D的對數(shù)標(biāo)準(zhǔn)差。

與傳統(tǒng)的易損性分析不同的是本文采用主余震序列的譜加速度Sa,MA為主余震序列強度指標(biāo)。建立Sa,MA和主余震作用下地震需求DMA之間的統(tǒng)計關(guān)系：

式中：a和b為回歸系數(shù)；為結(jié)構(gòu)在強度為Sa,MA的主余震序列作用下的地震需求DMA的中位值；為結(jié)構(gòu)在強度為Sa,MA的主余震序列作用下的地震需求DMA的對數(shù)標(biāo)準(zhǔn)差；DMA,i為結(jié)構(gòu)在第i條主余震序列作用下的需求。

主余震序列作用下的結(jié)構(gòu)地震易損性函數(shù)：

4.2 主余震作用下低延性結(jié)構(gòu)概率地震需求參數(shù)計算

采用前述主余震序列作為地震動輸入，對低延性結(jié)構(gòu)算例進行非線性時程分析。結(jié)構(gòu)最大層間位移角θmax為工程需求參數(shù)，根據(jù)式(4)和式(5)，求低延性結(jié)構(gòu)在主余震作用下的概率地震需求參數(shù)a、b和如圖6所示。

圖6結(jié)果表明，主余震序列作用下低延性結(jié)構(gòu)需求要大于單獨主震作用下的結(jié)構(gòu)需求。基于重復(fù)法構(gòu)造的主余震序列作用下，結(jié)構(gòu)需求最大，改變余震加速度時程的方向時，結(jié)構(gòu)需求不同。由該低延性結(jié)構(gòu)需求中位值與主余震序列強度指標(biāo)的擬合結(jié)果可知，基于重復(fù)法構(gòu)造的主余震序列作用下的對數(shù)線性擬合結(jié)果斜率最大，這說明基于重復(fù)法構(gòu)造的主余震序列較真實主余震序列有更大的破壞力。兩次余震的真實主余震序列作用下，對應(yīng)結(jié)構(gòu)需求略大于一次余震作用下的結(jié)構(gòu)需求中位值的擬合結(jié)果，這表明多次余震作用對該低延性結(jié)構(gòu)的附加損傷較小。

圖6 不同主余震序列作用下算例低延性RC框架ln(θm a x )_ l n(Sa)擬合結(jié)果Fig.6 Fitting results of l n(θm a x )_ l n(Sa)of investigated low-ductile RC frame under different mainshock-aftershock sequences

4.3 極限狀態(tài)的劃分及易損性曲線

考慮低延性RC框架梁柱節(jié)點的剪切破壞及柱的強度、剛度退化行為，根據(jù)文獻將結(jié)構(gòu)的損傷狀態(tài)劃分為：輕微破壞、中等破壞、嚴(yán)重破壞和完全破壞四種[23]。定義這四種損傷狀態(tài)對應(yīng)的抗震能力中位值和對數(shù)標(biāo)準(zhǔn)差，如表3所示。

表3 損傷狀態(tài)的劃分及其能力參數(shù)Table 3 Damage states and the corresponding capacity parameters

通過式(6)，結(jié)合圖6中的概率地震需求參數(shù)擬合結(jié)果及表3中的抗震能力參數(shù)，計算結(jié)構(gòu)在僅有主震作用和不同主余震序列作用下的易損性曲線，如圖7所示。圖7結(jié)果表明，在同一地震作用工況下，隨著低延性結(jié)構(gòu)損傷狀態(tài)的提高，對應(yīng)損傷狀態(tài)的超越概率減小。這一結(jié)果與一般鋼筋混凝土結(jié)構(gòu)在地震作用下的易損性分析結(jié)果一致。

圖7 不同主余震序列作用下算例低延性RC框架易損性曲線Fig.7 Fragility curves of the investigated low-ductile RC frame under different mainshock-aftershock sequences

為了對比不同主余震序列作用對低延性結(jié)構(gòu)易損性的影響，分別給出不同工況下低延性結(jié)構(gòu)輕微破壞、中等破壞、嚴(yán)重破壞和完全破壞四種損傷狀態(tài)的易損性曲線，如圖8所示，其橫坐標(biāo)Sa(T1, 5%)表示阻尼比為5%的結(jié)構(gòu)基本周期對應(yīng)譜加速度，縱坐標(biāo)Pf表示不同損傷狀態(tài)的超越概率。

圖8結(jié)果表明，對同一損傷狀態(tài)，主余震序列作用下該低延性結(jié)構(gòu)的易損性曲線均在僅主震作用下的易損性曲線之上。這表明低延性結(jié)構(gòu)在經(jīng)歷主震之后，再次經(jīng)歷余震作用將增加低延性結(jié)構(gòu)的損傷，也表明低延性結(jié)構(gòu)在主震作用下有一定的損傷，將會降低結(jié)構(gòu)的抗震能力。

對輕微破壞狀態(tài)，僅有主震和不同主余震序列作用下，結(jié)構(gòu)的易損性曲線接近。隨著損傷程度的提高，對中等破壞、嚴(yán)重破壞和完全破壞狀態(tài)，結(jié)構(gòu)在主余震序列作用下的易損性曲線明顯高于僅主震作用下的易損性曲線，相應(yīng)損傷狀態(tài)的超越概率增加明顯。

對比真實主余震記錄和基于重復(fù)法構(gòu)造的主余震序列作用下低延性結(jié)構(gòu)的易損性曲線可發(fā)現(xiàn)，除了輕微損傷狀態(tài)，真實主余震記錄作用下低延性結(jié)構(gòu)的易損性曲線明顯低于重復(fù)法構(gòu)造的主余震序列作用下的易損性曲線。這在一定程度上說明基于重復(fù)法構(gòu)造的主余震序列不能很好地反應(yīng)余震的特性及其對結(jié)構(gòu)的潛在破壞力。對嚴(yán)重破壞和完全破壞狀態(tài)，基于重復(fù)法構(gòu)造的不同余震方向的主余震序列作用下超越概率也有明顯的不同。當(dāng)基本周期對應(yīng)的譜加速度取0.25g時，在完全破壞損傷狀態(tài)下，基于重復(fù)法構(gòu)造的主余震序列其余震為正向和余震為反向的超越概率較僅在主震作用下的超越概率的增量百分比分別為 156.67%和117.83%，真實記錄的主震和一條余震記錄和真實的主震和兩條余震記錄作用下的超越概率增量百分比分別為43.86%和74.86%。對比主震后兩次余震和主震后一次余震作用下的易損性曲線可知，對輕微破壞狀態(tài)，余震作用的次數(shù)對低延性結(jié)構(gòu)的易損性影響不大，隨著結(jié)構(gòu)損傷程度的提高，相對于一次余震作用，兩次余震作用下相應(yīng)損傷狀態(tài)的超越概率明顯增大。綜上，余震對低延性結(jié)構(gòu)的潛在破壞力在結(jié)構(gòu)地震易損性分析中不可忽視。隨著結(jié)構(gòu)損傷程度的提高，余震對結(jié)構(gòu)的破壞力也將提高?；谥貜?fù)法構(gòu)造的主余震序列會高估主余震序列對結(jié)構(gòu)的破壞力。

圖8 不同主余震序列作用下算例低延性RC框架易損性曲線對比Fig.8 Comparison of fragility curves of the investigated low-ductile RC frame under different mainshock-aftershock sequences

5 結(jié)論

基于 OpenSees有限元分析軟件，以一未進行抗震設(shè)計的低延性RC框架結(jié)構(gòu)為研究對象，根據(jù)低延性結(jié)構(gòu)的特點，較為精確地模擬了低延性結(jié)構(gòu)梁柱節(jié)點剪切破壞、柱的強度和剛度退化行為以及梁端縱向鋼筋的粘結(jié)滑移。分別采用主震和不同主余震序列為地震動輸入，對該結(jié)構(gòu)進行動力分析，獲得不同損傷狀態(tài)對應(yīng)的地震易損性曲線，進而分析不同主余震序列作用對其抗震性能的影響。分析結(jié)果表明：

(1)本文提出的考慮低延性結(jié)構(gòu)梁柱節(jié)點的剪切破壞、柱的強度和剛度退化行為以及梁端縱向鋼筋的粘結(jié)滑移的模擬方法能較好地模擬低延性結(jié)構(gòu)梁柱節(jié)點剪切破壞和由剪切作用引起的剛度和強度退化現(xiàn)象。

(2)余震對低延性結(jié)構(gòu)會造成附加損傷，主余震序列作用下低延性 RC框架結(jié)構(gòu)對應(yīng)損傷狀態(tài)的超越概率高于僅主震作用下結(jié)構(gòu)對應(yīng)損傷狀態(tài)的超越概率，且隨著結(jié)構(gòu)損傷程度的提高，提高越明顯。

(3)基于重復(fù)法構(gòu)造的主余震序列作用下對應(yīng)損傷狀態(tài)的超越概率高于真實主余震序列作用下的超越概率，說明基于重復(fù)法構(gòu)造的主余震序列會高估其對結(jié)構(gòu)的破壞力。