Seismic Fragility Analysis of RC Frame Structure Based on Energy Balance
-
摘要: 基于“强柱弱梁”的屈服机制,依据能量平衡方法设计了某6层RC框架结构,采用震级-震中距条带地震动记录选取方法,选取12条随机地震动,利用Perform-3D有限元分析软件对结构进行增量动力(IDA)分析,得到了结构的地震易损性曲线、破坏状态概率曲线以及结构破坏概率矩阵。分析结果表明:该方法设计的结构能够形成预设的“强柱弱梁”屈服机制,可以保证结构中梁充分参与耗能,同时结构具有较强的抗倒塌能力,可以满足“小震不坏,中震可修,大震不倒”的性能要求。Abstract: 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.
-
引言
关于场地地震反应的分析已有大量研究成果,研究表明土壤在地震作用下会表现出材料非线性效应ADDIN EN.CITE.DATA(Joyner等,1975;Huang等,2001;Arslan等,2006;Hosseini等,2012)。等效线性化方法ADDIN EN.CITE.DATA(Schnabel等,1972;Idriss等,1992;Bardet等,2000;王笃国等,2016)是一种频域方法,通过在不同土体应变条件下选择等效阻尼比和剪切模量,将非线性问题转化为线性问题。当采用材料非线性本构模型描述土体非线性时,需采用时间积分算法求解非线性动力有限元方程。时间积分算法可分为隐式方法和显式方法。隐式算法每时刻需求解线性代数方程组,计算效率相对较低,如Wilson-θ法和Newmark法等。显式算法无需求解线性代数方程组,适合于强非线性和自由度数目较大的问题。研究者已提出多种显式时间积分算法ADDIN EN.CITE.DATA(Chung等,1994;王进廷等,2002;Belytschko等,2014)。作者近期提出一种二阶精度的单步显式算法,该算法适合变时步问题,在线弹性范围内稳定性较好。本文将该算法推广至求解非线性动力有限元方程中,并将其应用于地震波垂直入射时非线性地震反应分析。
1. 非线性动力有限元方程的显式时间积分算法
设已知非线性体系第${t_i}$时步的受力状态,求解第${t_{i + 1}}$时步的非线性结构动力学方程:
$${\boldsymbol{M}}{{\boldsymbol{\ddot u}}_{i + 1}}{\boldsymbol{ + C}}{{\boldsymbol{\dot u}}_{i + 1}} + {\boldsymbol{f}}_{i + 1}^S{\boldsymbol{ = }}{{\boldsymbol{f}}_{i + 1}}$$ (1) 式中M、C、${{\boldsymbol{f}}^S}$和${\boldsymbol{f}}$分别表示非线性体系的质量矩阵、阻尼矩阵、内力向量和外荷载向量;u表示位移,点号对时间t求导,i+1表示第${t_{i + 1}}$时刻。第i+1时刻时间步长为:
$${\boldsymbol{\Delta }}{t_i} = {t_{i + 1}} - {t_i}$$ (2) 文献显式方法求解非线性方程(1)的过程如下,第i+1时刻位移${{\boldsymbol{u}}_{i + 1}}$为:
$${{\boldsymbol{u}}_{i + 1}} = {{\boldsymbol{u}}_i} + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}{{\boldsymbol{\dot u}}_i} + \frac{{\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}^2}}{2}{{\boldsymbol{\ddot u}}_i}$$ (3) 第i+1时刻位移增量$\mathit{\Delta }{{\boldsymbol{u}}_i}$、内力增量$\mathit{\Delta }{\boldsymbol{f}}_i^S$和内力全量${\boldsymbol{f}}_{i + 1}^S$分别为:
$$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{{\boldsymbol{u}}_i} = {{\boldsymbol{u}}_{i + 1}} - {{\boldsymbol{u}}_i}$$ (4) $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{f}}_i^S = {\boldsymbol{f}}(\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{{\boldsymbol{u}}_i})$$ (5) $${\boldsymbol{f}}_{i + 1}^S = {\boldsymbol{f}}_i^S + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{f}}_i^S$$ (6) 第i+1时刻预估速度${{\boldsymbol{\dot {\tilde u}}}_{i + 1}}$、预估加速度${{\boldsymbol{\ddot {\tilde u}}}_{i + 1}}$、速度${{\boldsymbol{\dot u}}_{i + 1}}$和加速度${{\boldsymbol{\ddot u}}_{i + 1}}$分别为
$${{\boldsymbol{\dot {\tilde u}}}_{i + 1}} = {{\boldsymbol{\dot u}}_i} + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}{{\boldsymbol{\ddot u}}_i}$$ (7) $${{\boldsymbol{\ddot {\tilde u}}}_{i + 1}} = {{\boldsymbol{M}}^{ - 1}}({{\boldsymbol{f}}_{i + 1}} - {\boldsymbol{C\dot {\tilde u}}}_{i + 1}^{} - {\boldsymbol{f}}_{i + 1}^S)$$ (8) $${{\boldsymbol{\dot u}}_{i + 1}} = {{\boldsymbol{\dot u}}_i} + \frac{{\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{t_i}}}{2}({{\boldsymbol{\ddot u}}_i} + {{\boldsymbol{\ddot {\tilde u}}}_{i + 1}})$$ (9) $${{\boldsymbol{\ddot u}}_{i + 1}} = {{\boldsymbol{M}}^{ - 1}}({{\boldsymbol{f}}_{i + 1}} - {\boldsymbol{C\dot u}}_{i + 1}^{} - {\boldsymbol{f}}_{i + 1}^S)$$ (10) 式(3)—式(10)为求解式(1)的显式算法。算法中需由位移增量计算内力增量,目前常用的应力计算方法包括向前欧拉法、向后欧拉法和完全隐式计算法等ADDIN EN.CITE.DATA(Sloan等,1992;2001;Ahadi等,2003)。下面给出式(5)由位移增量计算内力增量的过程,即一种带误差控制的修正欧拉算法。
对于每个有限单元,由位移增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{u}}_i^e$计算应变增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e$的表达式为:
$$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e = {{\boldsymbol{B}}^e}\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{u}}_i^e$$ (11) 式中Be为应变矩阵。将ti时刻单元应变增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e$赋值给子步应变增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$,ti时刻单元应力${\boldsymbol{ \pmb{\mathit{ σ}} }}_i^e$赋值给${\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e$,初始化子步应变增量和应力状态分别为:
$$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e \leftarrow \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e$$ (12) $${\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e \leftarrow {\boldsymbol{ \pmb{\mathit{ σ}} }}_i^e$$ (13) 每个子步中应力增量计算思路见图 1,具体计算公式如下:
$${\boldsymbol{D}}_1^e = {\boldsymbol{D}}({\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e)$$ (14) $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e = {\boldsymbol{D}}_1^e\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (15) $${\boldsymbol{D}}_2^e = {\boldsymbol{D}}({\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e)$$ (16) $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_2^e = {\boldsymbol{D}}_2^e\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (17) $$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_s^e = \frac{{\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_2^e}}{2}$$ (18) 式中${{\boldsymbol{D}}^e}$为单元应力-应变关系矩阵。判断每个子步中应力增量$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{{\boldsymbol{ \pmb{\mathit{ σ}} }}_s}$是否符合精度要求的误差判断式为:
$${e_r} = \frac{{\left\| {\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_1^e - \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_2^e} \right\|}}{{\left\| {{\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_s^e} \right\|}}$$ (19) 判断误差er是否小于预先给定的判断值st,条件不满足时,缩小子步应变增量为:
$$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e \leftarrow A\sqrt {{{{s_t}} \mathord{\left/ {\vphantom {{{s_t}} {{e_r}}}} \right. } {{e_r}}}} \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (20) 式中A为误差峰值系数。采用缩小的子步应变增量重新进行式(14)—式(19)的计算与判断,循环直至满足精度要求,更新剩余应变增量和应力状态分别为:
$$\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e \leftarrow \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_i^e - \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ ε}} }}_s^e$$ (21) $${\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e \leftarrow {\boldsymbol{ \pmb{\mathit{ σ}} }}_{i + 1}^e + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }}{\boldsymbol{ \pmb{\mathit{ σ}} }}_s^e$$ (22) 利用更新剩余应变增量和应力状态循环执行式(14)—式(20),直至剩余应变增量小于等于零结束。
利用求得的第i+1时刻单元应力可得到单元应力增量和内力增量分别为:
$$ \Delta \boldsymbol{\sigma }_i^e = \boldsymbol{\sigma }_{i + 1}^e - \boldsymbol{\sigma }_i^e $$ (23) $$ \Delta {\boldsymbol{f}}_i^S{\rm{ = }}\sum\limits_e {\int {{{\boldsymbol{B}}^{e{\rm{T}}}}\boldsymbol{\Delta }{\boldsymbol{\sigma }}_i^e{\bf{d}}A} } $$ (24) 2. 地震波垂直入射时场地非线性地震反应分析
本节将上述非线性有限元方程的显式时间积分算法应用于地震波垂直入射时场地非线性地震反应分析中。假定基岩为线弹性半空间,考虑基岩上覆土层的材料非线性,不考虑土体阻尼。在土层下部设置黏性边界条件模拟半空间基岩的辐射阻尼,并在该处以等效结点力的方式实现地震动输入。
计算模型见图 2,选取A点作为观测点。土体非线性材料本构模型选取邓肯-张模型,土体线弹性参数见表 1,未给出配套的非线性参数,故算例中的非线性参数参考实际情况选取,后续研究中将使用更真实表现土体非线性行为的本构模型及真实工程场地参数。算例中的大气压参数取100kPa,内摩擦角增量取0°。入射地震动分别选取狄拉克脉冲和实测地震动(Gilroy Array #3,Coyote Lake, 1979)。入射狄拉克脉冲见图 3,观测点结果见图 4,实测地震动见图 5,观测点结果见图 6。图 4、图 6中给出采用中心差分法的计算结果作为参考解,由图 4、图 6可知,本文算法与中心差分法计算结果吻合较好,说明本文算法的有效性。
表 1 土层参数Table 1. Parameters of soils土质 深度/
m$\rho $/
(g/cm3)cs /
(m/s)v
-EN
-Rf
-c/
(MPa)θ/(°) D
-F
-人工填土 0—1.0 1.9 140 0.33 0.33 0.758 0.084 26.9 1.06 0.021 全新世砂土 1.0—5.1 1.9 140 0.32 0.33 0.758 0.084 26.9 1.06 0.021 全新世砂土 5.1—8.3 1.9 170 0.32 0.36 0.768 0.120 31.0 1.11 0.015 更新世粘土 8.3—11.4 1.9 190 0.40 0.44 0.822 0.188 28.4 1.01 0.012 更新世粘土 11.4—17.2 1.9 240 0.30 0.44 0.822 0.188 28.4 1.01 0.012 更新世砂土 17.2—22.2 2.0 330 0.26 0.51 0.840 0.300 30.0 1.02 0.011 基岩 >22.2 2.0 330 0.26 - - - - - - 表 1中ρ、cs、v、EN、Rf、c、θ为模型参数,分别表示密度、剪切波速、泊松比、无量纲幂次、破坏比、土的内聚力、土的摩擦角。D、F为试验常数。
3. 结论
本文发展一种求解材料非线性结构动力学方程的显式时间积分算法,并应用于地震波竖直入射时非线性地震反应分析中,通过算例验证了该方法的有效性。该显式算法具有无需对角阻尼矩阵、单步、稳定性良好等优点。本文考虑了邓肯-张非线性弹性本构模型,下步研究可考虑将该显式算法扩展到弹塑性本构模型及更能反映土层真实变形的本构模型中。
-
表 1 12条地震动记录
Table 1. Records of 12 seismic waves
地震动 台站 震级MW 震中距R/km 峰值加速度aPGA/g Northridge Beverly Hills-14145 Mulhol 6.7 13.39 0.416 Cape Mendocino Rio Dell Overpass-FF 7.1 22.64 0.549 Superstition Hills Poe Road(temp) 6.5 11.2 0.45 Loma Prieta Capitola 6.9 15.2 0.53 Duzce, Turkey Bolu 7.1 41.27 0.728 Friuli,Italy Tolmezzo 6.5 15.8 0.35 Imperial Valley Delta 6.5 22 0.35 Northridge Canyon Country-WLC 6.7 12.4 0.48 San Fernando LA-Hollywood 6.5 22.8 0.21 Superstition Hills El Centro Imp.Co. 6.5 18.2 0.36 Hector Mine Hector 7.1 11.7 0.34 Kocaeli Turkey Duzce 7.1 15.4 0.36 表 2 不同性能状态下层间位移角限值
Table 2. Definition of drift angle between layers in different performance states
正常使用 基本可使用 修复后使用 生命安全 倒塌 [θ]≤1/550 1/550<[θ]≤1/220 1/220<[θ]≤1/100 1/100<[θ]≤1/50 [θ]>1/50 表 3 楼(屋)面荷载、结构构件参数
Table 3. Roof load and sizes of main structural members
楼层 柱/mm×mm 梁/mm×mm 恒载/kN·m-2 活载/kN·m-2 1-3 600×600 300×600 5.5 3.5 4-5 500×500 250×500 5.5 3.5 6 500×500 250×500 6.5 0.5 表 4 结构破坏概率矩阵
Table 4. Damage probability matrix of the structure
地震水准${S_a}({T_1}, \xi)$/g 性能水平/% 正常使用 基本可使用 修复后使用 生命安全 倒塌 能量方法设计 规范设计 能量方法设计 规范设计 能量方法设计 规范设计 能量方法设计 规范设计 能量方法设计 规范设计 多遇地震 87.39 86.6 12.61 13.4 0 0 0 0 0 0 设防地震 12.19 9.18 61.15 54.73 19.92 26.79 6.31 8.61 0.43 0.69 罕遇地震 1.2 0.76 20.82 13.59 18.89 20.8 55.97 60.62 3.12 4.23 -
白久林, 欧进萍, 2012.基于能量平衡的钢筋混凝土框架结构抗震塑性设计方法.建筑结构学报, 33(10):22-31. http://d.old.wanfangdata.com.cn/Periodical/jzjgxb201210004 白久林, 金双双, 欧进萍, 2017.防屈曲支撑-钢筋混凝土框架结构基于能量平衡的抗震塑性设计.建筑结构学报, 38(1):125-134. http://d.old.wanfangdata.com.cn/Periodical/jzjgxb201701014 陈昉健, 易伟建, 2015.近场地震作用下锈蚀钢筋混凝土桥墩的IDA分析.湖南大学学报(自然科学版), 42(3):1-8. http://d.old.wanfangdata.com.cn/Periodical/hndxxb201503001 龚思礼, 2003.建筑抗震设计手册.北京:中国建筑工业出版社. 黄悠越, 2012.基于构件性能的RC框架结构层间位移角性能指标限值研究.广州:华南理工大学. 李刚, 程耿东, 2004.基于性能的结构抗震设计——理论、方法与应用.北京:科学出版社. 缪志伟, 马千里, 叶列平, 2013.钢筋混凝土框架结构基于能量抗震设计方法研究.建筑结构学报, 34(12):1-10. http://d.old.wanfangdata.com.cn/Periodical/jzjgxb201312001 缪志伟, 叶列平, 2014.钢筋混凝土框架-剪力墙结构基于能量抗震设计方法的耗能需求计算与设计流程.建筑结构学报, 35(1):10-18. http://d.old.wanfangdata.com.cn/Periodical/jzjgxb201401002 施炜, 叶列平, 陆新征等, 2011.不同抗震设防RC框架结构抗倒塌能力的研究.工程力学, 28(3):41-48, 68. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201100215417 叶列平, 缪志伟, 程光煜等, 2014.建筑结构基于能量抗震设计方法研究.工程力学, 31(6):1-12, 20. http://cdmd.cnki.com.cn/Article/CDMD-10703-1015037045.htm 叶列平, 马千里, 缪志伟, 2009.结构抗震分析用地震动强度指标的研究.地震工程与工程振动, 29(4):10-22. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd200904002 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 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. doi: 10.1016/j.engstruct.2004.02.012 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. doi: 10.1061/(ASCE)0733-9445(2002)128:8(1046) Liao W. C., 2010. Performance-based plastic design of earthquake resistant reinforced concrete moment frames. Michigan:University of Michigan, 39-47. doi: 10.1007%2Fs11803-013-0181-1 -