• ISSN 1673-5722
  • CN 11-5429/P

基于中远场小震结构台阵观测数据的地震响应分析

尹建华 冀昆 崔建文 温瑞智 任叶飞

孙柏涛, 张博, 闫培雷. 基于IDA方法的高层剪力墙结构抗震能力影响因素分析[J]. 震灾防御技术, 2018, 13(4): 851-859. doi: 10.11899/zzfy20180412
引用本文: 尹建华, 冀昆, 崔建文, 温瑞智, 任叶飞. 基于中远场小震结构台阵观测数据的地震响应分析[J]. 震灾防御技术, 2018, 13(4): 911-920. doi: 10.11899/zzfy20180418
Sun Baitao, Zhang Bo, Yan Peilei. Analysis of Influencing Factors of Seismic Capability of High-rise Shear Wall Structure by Using Incremental Dynamic Method[J]. Technology for Earthquake Disaster Prevention, 2018, 13(4): 851-859. doi: 10.11899/zzfy20180412
Citation: Yin Jianhua, Ji Kun, Cui Jianwen, Wen Ruizhi, Ren Yefei. Analysis of Structural Response by Using Strong Ground Motion Recordings from Middle and Far-field Small Earthquake Events[J]. Technology for Earthquake Disaster Prevention, 2018, 13(4): 911-920. doi: 10.11899/zzfy20180418

基于中远场小震结构台阵观测数据的地震响应分析

doi: 10.11899/zzfy20180418
基金项目: 

国家重点研发计划项目 2017YFC1500801

国家自然科学基金 51778589

国家自然科学基金 51878632

黑龙江省自然科学基金 E2017065

详细信息
    作者简介:

    尹建华, 男, 生于1984年。博士研究生。研究方向:工程地震与强震动观测研究。E-mail:yinjhhust@163.com

    通讯作者:

    温瑞智, 男, 生于1968年。研究员。研究方向:工程地震与强震动观测研究。E-mail:ruizhi@iem.ac.cn

Analysis of Structural Response by Using Strong Ground Motion Recordings from Middle and Far-field Small Earthquake Events

  • 摘要: 近年来我国布设的强震动观测结构台阵越来越多,并在数次地震中捕获了一批宝贵的观测记录,这些记录既是对应用日益广泛的结构数值模拟结果的检验和参数标定的依据,也是结构健康监测的重要数据。本文选取昆明防灾减灾中心大楼观测台阵捕获的3次中远场小震记录作为研究对象,对比ANSYS结构数值模拟和观测值之间在时程以及谱特性等方面的响应差别,证明了在弹性范围下二者结果具有较好的一致性。最后,基于观测记录,采用高效简单的经验谱比法进行了结构模态参数快速识别,3次地震下的自振频率计算结果相对稳定,与数值建模计算结果的差值可控制在0.17Hz以内。
  • 近年来,随着中国城镇化进程的加快,高层建筑数量逐年增加,其中剪力墙结构凭借其整体性能好、刚度大、在水平作用下侧向变形小的优势,逐渐成为中国城镇的主要结构类型。中国80%的城镇位于Ⅶ度地震设防区,一旦发生破坏性地震,将造成巨大的经济损失和人员伤亡。为了评估现有高层钢筋混凝土结构的抗震能力,世界地震工程界进行了大量的研究,其中增量动力分析法是1种有效的分析方法。

    增量动力分析方法(Incremental Dynamic Analysis,IDA)考虑了结构抗震能力和需求的不确定性和随机性因素,能较合理地评估不同性能水准下结构抗震能力。并且,通过IDA法可以分析结构从弹性到弹塑性、逐渐损伤直至破坏倒塌的全过程,进而可进一步研究控制结构破坏程度的条件和寻找防止结构倒塌的有效措施(吕西林等,2012)。

    2000年,美国开创性地在FEMA-350(Federal Emergency Management Agency,2000)中引入了IDA的概念,以此作为未来新建建筑结构抗震设计的标准。Vamvatsikos等(2002)在前人的研究的基础上,总结并提出了增量动力分析方法,创造性地提出了hunt & fill点跟踪算法,为实际工程中增量动力分析方法的应用创造了先例。21世纪初以来,该方法已经被广泛用于结构抗震性能评估。

    对于一般框架结构而言,影响其抗震能力的主要因素有柱轴压比、梁柱配筋率、梁剪跨比等,而对于高层剪力墙结构来说,由于主要承重构件的变化导致影响其抗震能力的主要因素为抗震设防等级、剪力墙轴压比、剪力墙厚度、结构高度、剪力墙配筋率等,其中在结构高度以及设防等级对剪力墙结构抗震能力的影响方面研究甚少。因此,本文在增量动力分析的基础上,讨论了抗震设防等级和高度对剪力墙结构抗震能力的影响,以期为研究工程结构的地震破坏和损失提供一定依据。

    增量动力分析方法将1个或多个地震动记录输入到结构模型中,然后对地震动参数(如PGV、PGA等)按照一定的规则进行调幅。通过这种方式,获得了1组不同强度的地震动记录,继而在这组记录的激励下,可以得到该结构的弹塑性分析结果。根据分析结果,可得到地震动强度指标IM(Intensity Measures)和结构损伤指标DM(Damage Measures)之间的关系曲线,即IDA曲线。通过该曲线,可以研究不同强度地震作用下结构损伤的全过程及相关失效概率。

    对结构进行IDA分析的主要流程和步骤为(卜一等,2009汪梦甫等,2010):

    (1)严格按照抗震设计规范建立弹塑性模型;

    (2)按照一定的规则选取有代表性的地震动记录(至少10条);

    (3)将选取的地震动记录按照既定的调幅方法进行调幅并输入到结构中进行计算分析;

    (4)将得到的IMDM点连接成线,即为结构的IDA曲线。

    考虑到断层距、结构周期以及场地类型等因素的影响,谢礼立等(2003)基于地震动潜在破坏势提出了最不利地震动的概念,并且综合考虑中国规范中有关场地类型的划分以及地震动强度等级的相关规定,给出了确定最不利地震动的方法和对应于不同类型的场地和结构的最不利地震动记录,从中选取了10条地震动作为地震动输入记录,如表 1所示。

    表 1  最不利设计地震动详细信息
    Table 1.  The detailed information of the most unfavorable design ground vibration
    编号 震级 年份 地点 断层距/km PGA/g
    F2 6.69 1994 北岭 54.45 0.154
    F3 6.61 1971 圣费尔南多 22.63 0.32
    F4 5.01 1979 帝王谷#10 12.96 0.054
    F5 7.36 1952 克恩县 38.89 0.18
    F7 6.95 1940 埃尔森特罗 6.09 0.281
    F10 6.36 1983 科林加 29.48 0.274
    F11 6.53 1979 帝王谷#6 1.35 0.449
    F12 6.19 1966 帕克菲尔德 12.9 0.272
    N1 6.7 1988 澜沧 17.6 0.207
    N4 7.8 1976 宁河 67 0.15
    下载: 导出CSV 
    | 显示表格

    表 1的地震动记录中,Ⅰ类场地的有F2、N1;Ⅱ类场地的有F3、F4、F5;Ⅲ类场地的有F7、F12;Ⅳ类场地的有F10、F11、N4。

    远场地震动记录的断层距大于20km,根据表 1中的信息可知属于远场地震动记录的有F2、F3、F5、F10,其余的均为近场地震动记录,由此可见选取的地震动记录充分考虑了断层距、场地类型等因素的影响。

    为了使结构分析结果具有普遍意义,有必要选择典型的结构,因此本文选择了剪力墙结构比较普遍的住宅楼,如图 1所示。

    图 1  结构平面布置
    Figure 1.  The plan layout of structure

    选取的工程为上海金桥地铁保障房,钢筋混凝土剪力墙结构体系,该建筑结构地上31层,总高度86.8m,层高2.8m,高宽比为5.36,安全等级为二级,丙类建筑,Ⅶ度设防,基本地震加速度值为0.10g,设计地震分组为第1组,Ⅳ类场地,特征周期为0.9s。

    为达到研究目的,本文以图 1为模板,严格按照中国规范规定分别建立了Ⅵ、Ⅶ、Ⅷ度3种不同抗震设防,15层、20层、25层、30层4种不同房屋高度的12个模型(中华人民共和国住房和城乡建设部,2010;中华人民共和国建设部,2010中华人民共和国住房和城乡建设部,中华人民共和国国家质量监督检验检疫总局,2010中华人民共和国住房和城乡建设部,2012),以此来研究不同抗震设防等级以及不同结构高度对建筑结构抗震能力的影响。

    本文选择Etabs进行结构的弹塑性分析。梁、柱采用三维框架单元模拟,剪力墙和楼板采用三维壳体单元模拟,为了反映混凝土材料在地震荷载作用下非线性阶段刚度退化的影响,其恢复力模型采用Takeda模型,而基于钢筋材料的特性,恢复力模型则采用考虑包辛格效应的Kinematic模型,强屈比取值为1.2,极限塑性应变取0.025。在进行弹塑性分析时采用显式积分法。

    逐一对所选择的地震动进行调幅,调幅后的PGA分别为0.005g、0.205g、0.455g、0.755g、1.105g和1.505g,地震动强度按XY向1:0.85输入,并按照文献(韩建平等,2007)的判断准则确定继续调幅或是终止分析。进而,可得到结构性能参数(最大层间位移角θmax)与地震动强度(地面峰值加速度PGA)之间的关系曲线,即IDA曲线。图 2绘出了12个不同模型的IDA曲线,图中3条IDA曲线分别为50%、16%、84%分位的IDA曲线。从图中可以看出,16%、84%分位的IDA曲线与结构的均值曲线,即50%分位的IDA曲线偏差并不大。

    图 2  IDA曲线
    Figure 2.  IDA curves of the structure

    工程结构的破坏可分为5个等级:基本完好,轻微破坏,中等破坏,严重破坏,倒塌。本文使用最大层间位移角作为损伤量化指标,结合邓明科等(2008)韩淼等(2011)关于高层钢筋混凝土剪力墙结构破坏状态的划分原则,给出本文不同破坏状态水平下的量化指标限值,如表 2所示。

    表 2  量化指标限值
    Table 2.  Limits of the quantitative indicators
    破坏等级 基本完好 轻微破坏 中等破坏 严重破坏 倒塌
    最大层间位移角θmax <1/950 1/950—1/500 1/500—1/175 1/175—1/90 >1/90
    下载: 导出CSV 
    | 显示表格

    根据上述破坏状态定义原则,对各模型的IDA曲线进行分析,为了得到一般性结论,选用50%分位曲线为研究对象,通过线性插值确定各个破坏状态下结构所能经受的地震强度指标值,具体数据见表 3

    表 3  结构各破坏状态对应的能力值(50%分位曲线)
    Table 3.  The corresponding ability value of each failure state of the structure
    模型 轻微破坏 中等破坏 严重破坏 倒塌
    θmax PGA/g θmax PGA/g θmax PGA/g θmax PGA/g
    Ⅵ度15层 0.00105 0.1222 0.002 0.2402 0.00571 0.7437 0.01111 1.242
    Ⅵ度20层 0.1195 0.2331 0.7289 1.218
    Ⅵ度25层 0.1174 0.2274 0.718 1.1753
    Ⅵ度30层 0.114 0.2191 0.7029 1.1328
    Ⅶ度15层 0.1483 0.2981 0.8138 1.3803
    Ⅶ度20层 0.1417 0.2857 0.8006 1.3413
    Ⅶ度25层 0.1348 0.2713 0.7858 1.3014
    Ⅶ度30层 0.1297 0.2592 0.7725 1.2784
    Ⅷ度15层 0.173 0.3647 0.9619 1.6712
    Ⅷ度20层 0.1654 0.3507 0.9449 1.6077
    Ⅷ度25层 0.1524 0.3284 0.9342 1.5273
    Ⅷ度30层 0.1505 0.3206 0.9223 1.4911
    下载: 导出CSV 
    | 显示表格

    根据表 3的数据,可以定量地分析出结构随抗震设防等级的提高在各个性能点处抗震能力提高的程度,为了更直观的表现,将以上数据绘制在图 3中。

    图 3  各破坏性能点随抗震设防等级变化的规律
    Figure 3.  The regular graph of the failed performance point varying with the earthquake resistance level

    图 3可知,无论结构高度如何,随着抗震设防等级的增加,各破坏状态点处的地震指标(PGA)都相应的增加,虽然轻微破坏和中等破坏所对应的PGA值有所增加,但是增加的幅度较小,而严重破坏和倒塌所对应的PGA值却大幅度增加。由此可见,相同地震动强度下,不同抗震设防等级的结构发生轻微破坏、中等破坏的概率变化并不大,但是抗震设防等级越低,发生严重破坏和倒塌的概率越大。

    根据表 3的数据,还可定量地分析出结构随高度的增加在各个性能点处抗震能力降低的程度,为了更直观地表现,将以上数据绘制在图 4中。

    图 4  各破坏性能点随结构高度变化的规律
    Figure 4.  The regular graph of the failed performance point varying with the structure height

    将结构层数信息用高度来表示,见表 4

    表 4  结构层数所对应高度
    Table 4.  The height value corresponding to the number of structural layers
    结构层数 15 20 25 30
    结构高度/m 41.2 55.2 69.2 83.2
    下载: 导出CSV 
    | 显示表格

    表 4图 4可以看出,Ⅵ、Ⅶ和Ⅷ度设防的结构,随着结构高度的增加,各破坏状态点处的地震指标PGA值都相应的降低,虽然轻微破坏、中等破坏和严重破坏所对应的PGA值有所降低,但是降低的幅度较小,而倒塌性能点所对应的PGA值却大幅度降低。由此可见,相同地震动强度下,不同结构高度的结构发生轻微破坏、中等破坏和严重破坏的概率变化并不大,但是结构高度的增加,会导致发生倒塌的概率显著增加。

    本文以Ⅵ、Ⅶ和Ⅷ度3种不同抗震设防等级和15层、20层、25层、30层4种不同结构高度共计12个剪力墙结构模型为基本算例,通过IDA对比分析,对不同抗震设防等级和不同结构高度的剪力墙结构的抗震性能进行了评估,总结出以下结论:

    (1)在对比分析不同抗震设防等级下剪力墙结构模型时发现,随着结构抗震设防等级的增加,结构的抗震能力随之增强,虽然在较小最大层间位移角处抗震设防等级的增加,并未导致PGA显著提高,但随着最大层间位移角的增大,抗震设防等级的提高所带来的抗震能力会有显著的提高。即抗震设防等级越高,结构在刚度、强度及延性等方面的表现越好,在遭受同等强度的地震动作用下结构也就越难发生轻微破坏、中等破坏和严重破坏。抗震设防等级的增大对增强结构抗震能力的效果很明显。

    (2)在对比分析不同结构高度下剪力墙结构模型时发现,随着结构高度的增加,结构的抗震能力随之降低,结构高度对抗震能力的影响并没有抗震设防等级对结构抗震能力的影响明显。但随着最大层间位移角的增加,较高的房屋其地震危险性还是比较大的。

  • 图  1  平面布置图(a)和侧视图(b)

    Figure  1.  Set-up in floor plan (a) and side view (b)

    图  2  结构三维效果图(a)和模型图(b)

    Figure  2.  3D model image (a) and sketch figure (b)

    图  3  观测台阵室内环境

    Figure  3.  Equipment layouts in observation room

    图  4  姚安地震各层加速度值

    Figure  4.  Simulated and observed acceleration time-histories at different floors in the Yao'an earthquake

    图  5  姚安余震地震各层加速度值

    Figure  5.  Simulated and observed acceleration time-histories at different floors in the Yao'an aftershock

    图  6  楚雄地震各层加速度值

    Figure  6.  Simulated and observed acceleration time-histories at different floors in the Chuxiong earthquake

    图  7  第9层加速度记录时频分析图

    Figure  7.  Time-frequency analysis plot of acceleration at floor#9

    图  8  姚安地震加速度反应谱

    Figure  8.  Acceleration response spectrum of the Yao'an earthquake

    图  9  姚安余震加速度反应谱

    Figure  9.  Acceleration response spectrum of the Yao'an aftershock

    图  10  楚雄地震加速度反应谱

    Figure  10.  Acceleration response spectrum of the Chuxiong earthquake

    图  11  3次地震下H/V谱比曲线

    Figure  11.  The H/V spectral ratio curves of three earthquake events

    表  1  梁、柱截面尺寸

    Table  1.   Section size for column and beam

    层数 梁/mm 外柱/mm 内柱/mm
    1—5 350×650 600×600 700×700
    6—9 350×650 500×500 600×600
    下载: 导出CSV

    表  2  各地震事件相关参数

    Table  2.   Earthquake parameters of three events

    地震事件 震级 震源深度/km 震中距/km 发生时间 B1层EW向PGA/cm·s-2 B1层NS向PGA/cm·s-2
    姚安 6.3 6 175 2009年7月9日19时19分14秒 4.8 6.6
    姚安余震 4.7 13 175 2009年7月9日20时57分31秒 1.2 1.1
    楚雄 5.2 20 93 2010年2月25日12时56分51秒 1.9 2.2
    下载: 导出CSV

    表  3  基于经验公式结构自振频率

    Table  3.   Structure natural frequency with empirical formula

    自振频率/Hz 模拟值 张相庭经验公式 公茂盛修正公式
    1阶 1.180 0.946 1.206
    2阶 2.237 2.056
    3阶 2.425 3.051
    下载: 导出CSV

    表  4  基于谱比法下的结构自振频率

    Table  4.   Structural natural frequency identification result by using spectra ratio method

    姚安 姚安余震 雄楚 模拟值
    EW/V NS/V H/V EW/V NS/V H/V EW/V NS/V H/V
    自振频率/Hz 1.221 1.685 1.685 1.270 1.343 1.306 1.343 1.086 1.086 1.180
    相差值/Hz 0.041 0.505 0.505 0.09 0.163 0.126 0.163 -0.094 -0.094
    下载: 导出CSV
  • 程文瀼, 王铁成, 颜德姮, 2005.混凝土结构-上册-混凝土结构设计原理.3版.北京:中国建筑工业出版社.
    公茂盛, 谢礼立, 欧进萍, 2010.基于结构强震记录框剪结构自振周期经验公式.土木工程学报, 43(12):35-40. http://d.old.wanfangdata.com.cn/Conference/7378874
    金星, 韦永祥, 张红才等, 2009.基于强震观测的隔震结构地震反应分析.地震工程与工程振动, 29(2):19-28. http://d.old.wanfangdata.com.cn/Periodical/dzgcygczd200902004
    李国强, 李杰, 苏小卒, 2008.北京:建筑结构抗震设计.2版.北京:中国建筑工业出版社.
    李鸿晶, Celebi M.2001.一幢钢筋混凝土建筑结构的强地震反应观测.地震工程与工程振动, 21(4):35-39. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK200101560007
    梁远森, 许红, 王云昌, 2005.高层建筑结构的自振周期的计算与实测.河南科学, 23(5):699-703. doi: 10.3969/j.issn.1004-3918.2005.05.023
    罗桂纯, 刘澜波, 齐诚等, 2011.基于地脉动和地铁振动的钢筋混凝土建筑结构响应分析.地球物理学报, 54(10):2808-2715. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201110028
    罗桂纯, 李小军, 王玉石等, 2014.关于中村(Nakamura)方法分析结构响应有效性的讨论.地震学报, 36(3):491-499. doi: 10.3969/j.issn.0253-3782.2014.03.014
    王飞, 刘英华, 任志林等, 2015.北京市防震减灾中心结构地震反应观测及振动特性识别.震灾防御技术, 10(3):539-546. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20150307&journal_id=zzfyjs
    王卫争, 2007.高层建筑结构地震反应观测台阵优化布设方法研究.大连: 大连理工大学. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y1092397
    张相庭, 1997.高层建筑抗风抗震设计计算.上海:同济大学出版社.
    周雍年, 2001.强震观测的发展趋势和任务.世界地震工程, 17(4):19-26. doi: 10.3969/j.issn.1007-6069.2001.04.004
    Fäcke A., Parolai S., Richwalski S. M., et al., 2006. Assessing the vibrational frequencies of the Cathedral of Cologne (Germany) by means of ambient seismic noise analysis. Natural Hazard, 38 (1-2):229-236. doi: 10.1007/s11069-005-8616-2
    Gallipoli M. R., Mucciarelli M., Castro R. R., et al., 2004. Structure, soil-structure response and effects of damage based on observations of horizontal-to-vertical spectral ratios of microtremors. Soil Dynamics and Earthquake Engineering, 24 (6):487-495. doi: 10.1016/j.soildyn.2003.11.009
    Nagata S., Takabatake D., Sato H., et al., 2015. Study on response characteristics and analytical models in substation steel structures based on seismic observation records. IEEJ Transactions on Power and Energy, 135 (12):743-752. doi: 10.1541/ieejpes.135.743
  • 加载中
图(11) / 表(4)
计量
  • 文章访问数:  132
  • HTML全文浏览量:  20
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-03-20
  • 刊出日期:  2018-12-01

目录

/

返回文章
返回