2. 兰州理工大学 西部土木工程防灾减灾教育部工程研究中心,甘肃 兰州 730050;
3. 山东建筑大学 建筑结构鉴定加固与改造省重点实验室,山东 济南 250101
2. Western Eng. Research Center of Disaster Mitigation in Civil Eng. of Ministry of Education, Lanzhou Univ. of Technol., Lanzhou 730050, China;
3. Shandong Provincial Key Lab. of Appraisal and Retrofitting in Building Structures, Shandong Jianzhu Univ., Jinan 250101, China
地震易损性分析是一种基于概率的结构抗震性能评估方法,它从概率意义上描述地震动强度与结构损伤状态之间的关联性,从而定量评估结构的抗震性能[1]。目前,易损性分析理论被广泛应用于各类结构的抗震性能评估中[2–5],特别是钢筋混凝土框架结构的地震易损性研究更是得到国内外学者的广泛关注[6–7]。
如何合理地考虑地震作用和结构中存在的强不确定性是近年来工程结构地震易损性研究的重点,而地震动输入的不确定性则是影响结构响应分析结果的主要因素[8]。在非线性动力分析和增量动力分析中,结构地震响应分析结果的不确定性和地震动记录的选择与调整及用于地震动调整的强度指标(intensity measure,IM)有关。不同的地震动强度指标与结构最大响应的相关性差异很大,按照不同的强度指标选择和调整地震动输入也会使得所关心的工程需求参数(engineering demand parameters,EDP)分析结果的离散性很大[9]。Luco[10]提出,地震动强度指标的选择既要考虑其充分性,也要考虑其有效性。Baker等[11]的研究表明,多个IM组成的向量型IM比单一IM包含更多的地震动信息,能够更加充分地描述地面运动对结构的影响。另外一种选择和调整地震动记录的方法是基于谱匹配的方法,通过对地震动的谱形进行非线性修改,以期达到减小谱的变异性、改变谱的偏差,从而使其与目标谱匹配[12]。与远场地震动相比,目前对近场地震动特性及其对结构的影响研究还不够充分[13]。
基于此,作者首先选取实际近场和远场地震动记录各20条,运用谱匹配方法对其进行目标设计谱匹配,得到谱匹配后的近场和远场地震动记录各20条。选择地震动最大峰值速度(peak ground velocity,PGV)和结构基本周期对应加速度反应谱
引入偏差,即输入调整或匹配的新的一组地震动获得的结构工程需求参数(EDP,engineering demand parameter)的几何平均值与输入选定的原始地震动获得的结构工程需求参数的几何平均值的比值,其可以有效地表示地震动选择与调整方法的合理性和可靠性[12]。几何平均值表示为:
$ED{P_{{\rm{geo}}}} = \exp \left[ {\frac{1}{n}\sum\limits_{i = 1}^n {\ln ED{P_i}} } \right]$ | (1) |
因此,偏差表示为:
$bias = \frac{{ED{P_{{\rm{geo}}}}}}{{ED{P_{{\rm{geo,ref}}}}}}$ | (2) |
式中,
如图1所示,选取一按现行规范设计的5层钢筋混凝土框架结构宿舍楼为研究对象,工程所在地抗震设防烈度为8度,设计基本地震加速度为0.2g,场地类别为Ⅱ类,设计地震分组为第3组,梁柱混凝土强度等级为C35,纵向受力钢筋为HRB400,箍筋为HPB300。
采用OpenSees有限元软件进行建模分析。建模时,梁柱单元均采用基于纤维截面的DispBeamColumn单元,考虑
![]() |
图1 算例结构 Fig. 1 The structure for case study |
2.2 地震动记录选择与调整
综合考虑中国抗震设计规范的相关要求,从PEER地震动数据库中选取了一组相当于Ⅱ类场地的20条断层距在10 km以内,震级
![]() |
图2 原始和谱匹配后地震动记录的加速度反应谱 Fig. 2 Acceleration response spectra of original and spectrum-matched ground motion record |
3 单一地震动强度下的结构地震响应分析
以原始及匹配后的近场、远场地震动记录分别作为输入,对算例结构模型进行单一强度下的非线性动力分析,提取最大层间位移角和最大残余层间位移角,分析其离散性及偏差,并分析其与
计算不同输入下结构最大层间位移角和最大残余层间位移角等地震响应分析结果的均值、离散系数和谱匹配后的偏差如表1所示。
表1结果表明,与输入原始地震动记录的结构地震响应结果相比,输入匹配后的近场地震动和远场地震动记录所得结构地震响应结果的离散性均有明显降低,且最大层间位移角的离散性明显小于最大残余层间位移角的离散性。结构最大层间位移角的偏差明显优于最大残余层间位移角的偏差。
表1 结构地震响应分析结果的均值、离散系数及偏差 Tab. 1 Mean values, dispersion coefficients and biases of analyzed seismic response results of the structure |
![]() |
3.2 相关性分析
输入原始近场、远场地震动下结构最大层间位移角和最大残余层间位移角与
![]() |
图3 输入原始近场地震动,结构地震响应结果与PGV和Sa(T1)的相关性 Fig. 3 Correlation between seismic responses of the structure under original near-field ground motion and PGV, Sa(T1) |
![]() |
图4 输入原始远场地震动,结构地震响应结果与PGV和Sa(T1)的相关性 Fig. 4 Correlation between seismic responses of the structure under original far-field ground motion and PGV, Sa(T1) |
输入谱匹配后的近场、远场地震动下结构最大层间位移角和最大残余层间位移角与
![]() |
图5 输入谱匹配的近场地震动,结构地震响应结果与PGV和Sa(T1)的相关性 Fig. 5 Correlation between seismic responses of the structure under spectrum-matched near-field ground motion and PGV, Sa(T1) |
![]() |
图6 输入谱匹配的远场地震动,结构地震响应结果与PGV和Sa(T1)的相关性 Fig. 6 Correlation between seismic responses of the structure under spectrum-matched far-field ground motion and PGV, Sa(T1) |
4 IDA分析
以
基于IDA结果,分别计算4组地震动记录作用下的16%、50%、84%分位数的IDA曲线。限于篇幅,本文仅给出以
图7结果表明:
1)在
2)弹塑性阶段,以
3)对比近场和远场地震动记录谱匹配前后的50%分位数IDA曲线,以
![]() |
图7 50%分位数IDA曲线 Fig. 7 IDA curves with 50% fractiles |
5 结构易损性分析
参考中国现行规范给出的4种损伤状态的最大层间位移角控制指标:基本完好(1/550,弹性限值)、轻微破坏(1/250)、中等破坏(1/120)和不严重破坏(1/60)[14],以及结构倒塌状态最大层间位移角建议控制指标(1/20)[15],利用IDA结果,采用正态分布模型进行参数估计,得到以最大层间位移角表示的五种损伤状态的易损性曲线。
表2给出了不同损伤状态50%超越概率相对应的
1)就近、远场地震动而言,以
表2 不同类型地震动下五种损伤状态50%超越概率相应的
|
![]() |
![]() |
图8 基于最大层间位移角的不同损伤状态相应地震易损性曲线 Fig. 8 Seismic fragility curves of different damage states based on maximum inter-story drift ratios |
2)随着结构损伤程度的增大,相应损伤状态50%超越概率对应的
3)当结构处于弹性阶段或进入弹塑性程度较弱时,地震动未进行谱匹配,对应
4)地震动谱匹配前后相比,以
5)无论近场还是远场地震动,地震动谱匹配前后,各损伤状态下50%超越概率对应的
以5种损伤状态下结构的最大层间位移角为参考控制指标,通过IDA分析得到相应5种损伤状态,不同地震动记录作用下结构的最大残余层间位移角,并进行统计分析,分析结果如表3所示。同样,图9仅给出中等破坏、不严重破坏和倒塌3种状态下结构的最大残余层间位移角分布,横轴中“1”代表原始近场地震动记录,“2”代表谱匹配后的近场地震动记录,“3”代表原始远场地震动记录,“4”代表谱匹配后远场地震动记录。
表3 不同地震动下相应5种损伤状态的最大残余层间位移角统计特性 Tab. 3 Statistical values of maximum residual inter-story drift ratios corresponding to five different damage states under different types of ground motion |
![]() |
![]() |
图9 不同损伤状态相应的最大残余层间位移角的分布 Fig. 9 Maximum residual inter-story drift ratios corresponding to different damage states |
1)随着结构损伤程度的增大,最大残余层间位移角逐渐增大。当结构进入强非线性阶段后,最大残余层间位移角的离散性随结构损伤程度的增大而逐渐降低。
2)结构倒塌前,近场地震动作用下各损伤状态对应的最大残余层间位移角平均值较远场地震动作用下的平均值大;达到倒塌状态时,近场地震动作用下结构的最大残余层间位移角平均值小于远场地震动作用下的平均值。
3)无论是近场还是远场地震动,地震动反应谱匹配前后,结构各损伤状态对应的最大残余层间位移角无明显偏差,同时谱匹配后最大残余层间位移角的离散性无明显减小。
6 结 论以实际近场、远场地震动记录及谱匹配后的地震动记录为输入,以
1)实际近、远场地震动作用下,结构最大层间位移角和最大残余层间位移角与
2)谱匹配后的近、远场地震动作用下,结构最大层间位移角和最大残余层间位移角分析结果的离散性明显降低且最大层间位移角无明显偏差。
3)选取不同的地震动强度指标,可能会得到不同的易损性分析结论。如以
4)结构倒塌前,近场地震动作用下各损伤状态对应的最大残余层间位移角平均值较远场地震动作用下的平均值大;达到倒塌状态时,近场地震动作用下结构的最大残余层间位移角平均值小于远场地震动作用下的平均值。
[1] |
Li Ning,Li Zhongxian,Li Yang. Seismic fragile surface analysis of reinforced concrete bridges subjected to near-fault ground motions[J]. Earthquake Engineering and Engineering Dynamics, 2014, 34(S1): 307-312. [李宁,李忠献,李杨. 近断层地震动作用下桥梁结构易损性曲面分析[J]. 地震工程与工程振动, 2014, 34(S1): 307-312. DOI:10.13197/j.eeev.2014.S0.307.lin.047] |
[2] |
He Yibin,Li Yan,Shen Pusheng. Performance-based seismic fragility analysis of tall hybrid structures[J]. Engineering Mechanics, 2013, 30(8): 142-147. [何益斌,李艳,沈蒲生. 基于性能的高层混合结构地震易损性分析[J]. 工程力学, 2013, 30(8): 142-147. DOI:10.6052/j.issn.1000-4750.2012.04.0302] |
[3] |
Rossetto T,Elnashai A. Derivation of vulnerability functions for European-type RC structures based on observational data[J]. Engineering Structures, 2003, 25(10): 1241-1263. DOI:10.1016/S0141-0296(03)00060-9 |
[4] |
Schotanus M I J,Franchin P,Lupoi A P,et al. Seismic fragility analysis of 3D structures[J]. Structural Safety, 2004, 26(4): 421-441. DOI:10.1016/j.strusafe.2004.03.001 |
[5] |
Karim K R,Yamazaki F. Effect of earthquake ground motions on fragility curves of highway bridge piers based on numerical simulation[J]. Earthquake engineering and Structural Dynamics, 2001, 30(12): 1836-1856. DOI:10.1002/eqe.97 |
[6] |
Lü Dagang,Jin Guoqing,Yu Xiaohui. Modified seismic fragility analysis of reinforced concrete frame structures in consideration of collapse probability[J]. Journal of Building Structures, 2016, 37(9): 26-62. [吕大刚,金国庆,于晓辉. 考虑倒塌概率修正的钢筋混凝土框架结构地震易损性分析[J]. 建筑结构学报, 2016, 37(9): 26-62.] |
[7] |
Deierlein G G,Liel A B,Haselton C B,et al.ATC-63 methodology for evaluating seismic collapse safety of archetype buildings [C]// Structures 2008:Crossing Borders.Vancouver,BC,Canada:ASCE,2008.DOI:10.1061/41016 (314)48
|
[8] |
Katsanos E I,Sextos A G,Manolis G D. Selection of earthquake ground motion records:A state-of-the-art review from a structural engineering perspective[J]. Soil Dynamics and Earthquake Engineering, 2009, 30(4): 157-169. DOI:10.1016/j.soildyn.2009.10.005 |
[9] |
Han Jianping,Chen Jiqiang,Yan Qing,et al. State of the art of ground motion intensity measures accounting for the influence of spectral shape[J]. Engineering Mechanics, 2015, 32(10): 9-17. [韩建平,陈继强,闫青,等. 考虑谱形影响的地震动强度指标研究进展[J]. 工程力学, 2015, 32(10): 9-17. DOI:10.6052/j.issn.1000-4750.2014.07.ST06] |
[10] |
Luco N.Probabilistic seismic demand analysis,SMRF connection fractures,and near-source effects [D].Stanford:Stanford University,2001.
|
[11] |
Baker J W,Cornell C A. A vector-valued ground motion intensity measure consisting of spectral acceleration and epsilon[J]. Earthquake Engineering and Structural Dynamics, 2005, 34(10): 1193-1217. DOI:10.1002/eqe.474 |
[12] |
Seifried A E.Response spectrum compatibilization and its impact on structural response assessment [D].Stanford:Stanford University,2013.
|
[13] |
Tothong P,Cornell C A. Structural performance assessment under near-source pulse-like ground motions using advanced ground motion intensity measures[J]. Earthquake Engineering and Structural Dynamics, 2008, 37(7): 1013-1037. DOI:10.1002/eqe.792 |
[14] |
GB50011-2010建筑抗震设计规范[S].北京:中国建筑工业出版社,2010.
|
[15] |
Han Jianping,Huang Linjie,Sun Xiaoyun. Global collapse resistance capacity investigation of infilled RC frame considering in-plane and out-of-plane interaction of infill walls[J]. Engineering Mechanics, 2016, 33(9): 146-154. [韩建平,黄林杰,孙小云. 考虑填充墙平面内外相互作用评估RC框架-填充墙结构抗整体性倒塌能力[J]. 工程力学, 2016, 33(9): 146-154. DOI:10.6052/j.issn.1000-4750.2015.02.0098] |