基于能量方法设计的RC框架结构易损性分析1

王中阳 车佳玲 张尚荣 包 超



基于能量方法设计的RC框架结构易损性分析1

王中阳 车佳玲 张尚荣 包 超

（宁夏大学，土木与水利工程学院，银川 750021）

基于“强柱弱梁”的屈服机制，依据能量平衡方法设计了某6层RC框架结构，采用震级-震中距条带地震动记录选取方法，选取12条随机地震动，利用Perform-3D有限元分析软件对结构进行增量动力（IDA）分析，得到了结构的地震易损性曲线、破坏状态概率曲线以及结构破坏概率矩阵。分析结果表明：该方法设计的结构能够形成预设的“强柱弱梁”屈服机制，可以保证结构中梁充分参与耗能，同时结构具有较强的抗倒塌能力，可以满足“小震不坏，中震可修，大震不倒”的性能要求。

能量平衡方法 失效模式 易损性分析 地震需求分析 增量动力分析

引言

近年来，世界各地地震频发，给人们带来巨大的经济损失。通过总结历次震害的经验教训，人们对地震动特性及结构动力特性的理解逐渐加深，同时广泛重视对结构抗震设计方法的研究。我国现行的抗震设计规范提出了“三水准”设防目标与“两阶段”设计方法。然而，从汶川地震的震灾中很少看到规范中要求的“强柱弱梁”型破坏；同时，由于地震动具有极大的不确定性与随机性，极震区及其周边区域（如唐山地震、汶川地震）的实际地震烈度往往比设防烈度大得多（施炜等，2011），导致大量建筑发生倒塌，未能实现其性能要求；此外，基于弹性方法设计的结构无法反映出结构预定的屈服机制以及地震输入能量对结构的影响。在控制结构屈服机制的基础上采用能量方法的抗震设计能够有效解决上述问题（白久林等，2012）。能量方法抗震设计的概念自Housner等（1956）提出后，经过国内外诸多学者的研究，理论已基本成熟。Leelataviwat等（2002）和Liao等（2010）提出基于能量平衡的塑性设计方法对混凝土框架结构进行了系列研究；白久林等（2012，2017）采用能量平衡的抗震塑性设计方法分别对钢筋混凝土框架结构和防屈曲支撑-钢筋混凝土框架结构进行了研究；叶列平等（2014）和缪志伟等（2013，2014）指出合理的结构损伤耗能机制控制是实现基于能量抗震设计前提，并给出了出钢筋混凝土框架-剪力墙结构基于能量抗震设计的实施流程。目前，对用此方法设计的结构地震易损性的研究相对较少。

本文以“强柱弱梁”屈服机制作为RC框架结构失效模式，采用能量方法对某6层RC框架结构进行抗震设计，在考虑地震动随机性对地震易损性分析影响的基础上对其进行增量动力分析（简称IDA），建立结构在不同地震水平下结构的易损性方程，预测结构在不同地震水平下发生各级破坏的概率，同时验证此方法的合理性，为结构抗震设计提供参考。

1 基于能量方法抗震设计

图1 结构能量平衡图

图2 结构最优失效模式

强震作用下结构最优失效模式如图2所示。每层耗能梁段均形成塑性铰，达到最终破坏时仅底层柱端形成塑性铰。梁作为主要耗能构件消耗大部分地震能量，防止楼层柱被破坏。

基于能量方法的抗震设计流程主要分为2部分：确定地震作用下结构设计参数和进行构件设计，具体步骤如下：

（3）利用结构所得到的基底剪力与楼层剪力，假定各楼层梁端和底层柱底形成塑性铰，利用能量平衡原理求出各楼层梁端的内力大小，再对边柱和中柱取隔离体，确定各楼层柱的内力大小。

（4）按式（1）对梁柱截面进行抗震验算，截面配筋设计按照《混凝土设计规范》（GB 50010—2010）（中华人民共和国住房和城乡建设部，2011）与《建筑抗震设计规范》（GB 50011—2010）（中华人民共和国住房和城乡建设部等，2010）进行计算。

2 结构易损性分析

2.1 地震动的选择

地震易损性分析中主要考虑的不确定性包括地震动不确定性和结构不确定性。其中，地震动的不确定性对结构反应影响较大。10—20条地震动记录可精确评估中高层建筑的抗震性能（陈昉健等，2015）。本文选取12条随机地震动对结构进行IDA分析，所选地震动震级W、震中距均匀分布在一个较宽的范围：6.5＜W＜7.1，13km＜＜30km，其对应的峰值加速度也均匀分布在一个较宽的范围内，具体数据如表1所示。

表1 12条地震动记录

2.2 地震动强度指标及性能水平的确定

地震动强度参数对结构易损性分析及结果模拟的有效性起着重要的作用。多高层RC框架结构地震反应以基本振型反应为主，a（1，5%）作为输入地震动参数时，可以更合理地揭示地震动对这类结构的破坏作用（叶列平等，2009）。因此，本文选取a（1，5%）作为地震动强度指标，同时选取最大层间位移角max作为结构损伤指标。

取a（1，5%）的最大值为1g，分析时，首次取a（1，5%）为0.05g，以0.1g为增量在0.05—1g间取值。将结构最大层间位移角max达到10%或IDA曲线上切线斜率为弹性斜率20%所对应的地震强度较小点作为结构倒塌点。

FEMA356定义了立即使用（IO）、生命安全（LS）和防止倒塌（CP）3个性能水平下结构层间位移角的限值（FEMA356，2000）；《建筑抗震设计规范》（GB 50011—2010）将结构的破坏划分为5个等级，并给出了相应的层间位移角参考值（中华人民共和国住房和城乡建设部等，2010）。本文结合FEMA356、《建筑抗震设计规范》（GB 50011—2010）（简称抗规）相关规定及黄悠越（2012）的相关研究成果，确定了各状态的层间位移角限值［］（表2）。

表2 不同性能状态下层间位移角限值

3 算例分析

3.1 工程概况

某6层办公楼，首层层高4.5m，其余层高3.6m，结构平面布置图如图3所示。结构抗震设防烈度Ⅷ度（0.2g），设计地震分组为第2组，Ⅱ类场地；梁柱混凝土强度等级为C30，梁柱纵筋采用HRB400，箍筋采用HPB300；地面粗糙类别B类，基本风压0.65kN/m2，梁柱构件截面参数与楼（屋）面荷载见表3。

图3 结构平面布置图

3.2 有限元分析模型

本文进行增量动力分析（IDA分析）时采用Perform-3D有限元分析软件。该软件由伯克利的Powell教授开发，具有完备的材料、构件模型库，稳定可靠的算法，受到国内外学术界及工程界的广泛应用与认可。由于结构平面布置规则，仅取中间一榀横向框架进行分析。结构构件中梁、柱单元采用纤维模型，两端的塑性铰区定义为非弹性纤维截面段，塑性铰长度取梁、柱截面高度的0.5倍。核心区约束混凝土采用Kent-park本构模型，钢筋采用非屈曲弹塑性本构模型，钢筋屈服后的刚度为初始刚度的1%。各材料强度参数取其平均强度值。此外，框架柱考虑-效应并假设楼板平面内刚度无限大，结构阻尼采用Rayleigh阻尼，不考虑模态阻尼。

表3 楼（屋）面荷载、结构构件参数

3.3 概率地震需求分析

利用Perform-3D软件对结构进行弹塑性时程分析，得到不同地震加速度下的最大层间位移角，如图4所示。图中每列竖向数据点为相同地震加速度下结构响应。水平虚线为不同性能水平限值，从上到下依次为倒塌、生命安全、修复后使用、基本可使用、正常使用5个等级。

从图4可以看出，规范设计与能量方法设计的结构随着地震动加速度的增大，结构层间位移角随之增大，同时层间位移角的离散程度也变大。相同地震动水平下，能量方法设计结构的层间位移角总体上小于规范方法设计结构。

图5给出了2种设计方法设计的结构在设防地震、罕遇地震作用下结构的对数正态分布概率密度函数，图中竖向虚线对应不同的性能水平依次为正常使用、基本可使用、修复后可使用、生命安全、倒塌5个等级。由图可以看出，无论是能量方法设计的结构还是规范设计的结构，在设防地震下，结构地震反应主要分布在层间位移角为0.002—0.008的区间，2种结构基本处于基本可使用状态，二者概率相差不大；罕遇地震下，2种结构地震反应大致分布在层间位移角为0.0046—0.015的区间，结构处于修复后可使用状态的概率最高，处于生命安全状态的概率次之，处于基本可使用状态的概率最低，但二者的破坏概率有明显不同，能量方法设计结构的破坏概率小于按照规范设计结构的破坏概率。

图4 最大层间位移角与Sa（T1，5%）关系

图5 层间位移角概率密度分布图

3.4 地震需求概率模型

假设地震动强度指标与结构需求参数服从指数正态分布（龚思礼，2003），其表达式为：

式中：、为回归系数。

将图4中的数据取对数，再进行线性回归，如图6所示。

图6 结构概率地震需求模型

基于能量方法抗震设计结构的地震需求模型线性回归方程为：

基于规范抗震设计结构的地震需求模型线性回归方程为：

式（4）、（5）拟合相关系数分别为0.9758、0.9722，接近1，说明回归对散点拟合程度较好。

3.5 结构地震的易损性曲线

地震易损性曲线作为结构易损性分析的一种形式，可以获得结构地震反应参数达到指定破坏状态下结构抗震能力参数的超越概率，其表达式为：

由图7可以看出，随着地震动的增大，2种结构从基本完好到倒塌，结构的易损性曲线逐渐变得平缓，超越概率逐渐变小，其符合结构设计准则，同时表4给出了2种方法设计结构的各极限状态的超越概率。

由图7及表4可知，在多遇地震下，能量方法设计结构与规范设计结构超过正常使用状态的概率为12.61%和13.4%；设防地震下能量方法设计结构与规范方法设计结构达到“正常使用”状态、“基本可使用”状态、“修复后使用”状态、“生命安全”状态的概率分别为12.19%、61.15%、19.92%、6.31%和9.18%、54.73%、26.79%、8.61%；罕遇地震下能量方法设计结构与规范设计结构分别有55.97%和60.62%的概率达到生命安全状态。相同地震动作用下按规范设计的结构破坏状态比能量方法设计的结构破坏状态严重。

ATC-63（2008）报告建议：“设防大震下倒塌概率小于10%即认为达到大震性能的要求”，框架结构的限值为2.3，为反应结构抗倒塌能力指标，值越大结构抗倒塌能力越强，其表达式为：

图7 结构易损性曲线

由图7可知，能量方法设计结构50%为1.09g，本文大震取为0.4g，计算得到罕遇地震下能量方法设计结构的为2.7，接近倒塌状态的概率仅为3.12%，远远低于10%。

表4 结构破坏概率矩阵

3.6 弹塑性应变能Ep在构件中的分配

图8 梁、柱弹塑性应变能所占比例

4 结论

本文以“强柱弱梁”耗能机制为目标，采用基于能量方法的抗震设计对某6层RC框架结构进行设计，并对其进行IDA分析，得出如下结论：

（1）通过IDA分析，获得了结构在不同强度地震下的动力响应，建立了结构易损性曲线，进而得到结构在不同地震水准下的破坏概率，为结构震灾损伤评估提供参。

（2）基于能量方法设计结构在多遇地震下基本处于正常使用极限状态，设防地震下达到基本可使用状态的概率为61.15%，罕遇地震下达到接近倒塌极限状态的概率为3.12%，满足现行抗震设计规范规定的“小震不坏、中震可修、大震不倒”的设防目标。

（3）罕遇地震下，能量方法设计结构的抗倒塌能力储备系数高于ATC-63（2008）报告规定的值，具有较强的抗倒塌能力，同时结构中梁构件能充分发挥其在整体结构耗能的优势。

白久林，欧进萍，2012．基于能量平衡的钢筋混凝土框架结构抗震塑性设计方法．建筑结构学报，33（10）：22—31．

白久林，金双双，欧进萍，2017．防屈曲支撑-钢筋混凝土框架结构基于能量平衡的抗震塑性设计．建筑结构学报，38（1）：125—134．

陈昉健，易伟建，2015．近场地震作用下锈蚀钢筋混凝土桥墩的IDA分析．湖南大学学报（自然科学版），42（3）：1—8．

龚思礼，2003．建筑抗震设计手册．北京：中国建筑工业出版社．

黄悠越，2012．基于构件性能的RC框架结构层间位移角性能指标限值研究．广州：华南理工大学．

李刚，程耿东，2004．基于性能的结构抗震设计——理论、方法与应用．北京：科学出版社．

缪志伟，马千里，叶列平，2013．钢筋混凝土框架结构基于能量抗震设计方法研究．建筑结构学报，34（12）：1—10．

缪志伟，叶列平，2014．钢筋混凝土框架-剪力墙结构基于能量抗震设计方法的耗能需求计算与设计流程．建筑结构学报，35（1）：10—18．

施炜，叶列平，陆新征等，2011．不同抗震设防RC框架结构抗倒塌能力的研究．工程力学，28（3）：41—48，68．

叶列平，缪志伟，程光煜等，2014．建筑结构基于能量抗震设计方法研究．工程力学，31（6）：1—12，20．

叶列平，马千里，缪志伟，2009．结构抗震分析用地震动强度指标的研究．地震工程与工程振动，29（4）：10—22．

中华人民共和国住房和城乡建设部，中华人民共和国国家质量监督检验检疫总局，2010．GB 50011—2010建筑抗震设计规范．北京：中国建筑工业出版社．

中华人民共和国住房和城乡建设部，2011．GB 50010—2010混凝土结构设计规范．北京：中国建筑工业出版社．

ATC-63, 2008. Quantification of building seismic performance factors ATC-63 Project Report—90% Draft FEMA P695/April 2008. Redwood City , California: Applied Technology Council.

Erberik M. A., Elnashai A. S., 2004. Fragility analysis of flat-slab structures. Engineering Structures, 26(7): 937—948.

FEMA 356, 2000. Prestandard and commentary for the seismic rehabilitation of buildings. Washington, DC: Federal Emergency Management Agency.

Housner G. W., 1956. Limit design of structures to resist earthquakes. In: Proceedings of the 1st World Conference on Earthquake Engineering. Oakland, Calif: Earthquake Engineering Research Institute, 1—12.

Leelataviwat S., Goel S. C., Stojadinović B., 2002. Energy-based seismic design of structures using yield mechanism and target drift. Journal of Structural Engineering, 128(8): 1046—1054.

Liao W. C., 2010. Performance-based plastic design of earthquake resistant reinforced concrete moment frames. Michigan: University of Michigan, 39—47.

Seismic Fragility Analysis of RC Frame Structure Based on Energy Balance

Wang Zhongyang, Che Jialing, Zhang Shangrong and Bao Chao

（School of Civil Engineering and Hydraulic Engineering, Ningxia University, Yinchuan 750021, China）

The energy-based seismic design method of RC framed structure of the “strong column with weak beam” was proposed, and 6 layer RC frame structures were designed by selecting ground motion records using magnitude epicentral distance strips in this paper. We selected 12 ground motions and analyzed the RC frames structure through the IDA method by Perform-3D, and obtained the seismic fragility curves, failure probability curves and the structural failure probability matrices. The results showed that the structure designed by this method can form the preset “strong column with weak beam” yielding mechanism. It makes promise that the beam in the structure is fully involved in energy dissipation with strong anti-collapse capacity, not only can satisfy fully operational level under frequent earthquake and fulfill reparably operational level under rare earthquake, but also satisfy near collapse level under extremely rare earthquakes.

Energy balance method; Failure mode; Vulnerability analysis; Seismic demand analysis; Incremental dynamic analysis

王中阳，车佳玲，张尚荣，包超，2018．基于能量方法设计的RC框架结构易损性分析．震灾防御技术，13（3）：524—533．

10.11899/zzfy20180304

国家自然科学基金项目（51408328、51608283），宁夏青年科技人才托举工程

2017-11-20

王中阳，男，生于1991年。硕士研究生。主要从事钢筋混凝土结构抗震研究。E-mail：2550115272@qq.com

车佳玲，女，生于1985年。副教授。主要从事抗震相关研究。E-mail：che_jialing@126.com