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

场地对地震动放大作用的方向相关性研究

卢滔 郭林松 佟鑫 卢玉林 霍敬妍

孙柏涛, 张博, 闫培雷. 基于IDA方法的高层剪力墙结构抗震能力影响因素分析[J]. 震灾防御技术, 2018, 13(4): 851-859. doi: 10.11899/zzfy20180412
引用本文: 卢滔,郭林松,佟鑫,卢玉林,霍敬妍,2022. 场地对地震动放大作用的方向相关性研究. 震灾防御技术,17(3):481−489. doi:10.11899/zzfy20220308. doi: 10.11899/zzfy20220308
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: Lu Tao, Guo Linsong, Tong Xin, Lu Yulin, Huo Jingyan. Study on Azimuth-related Site Amplification Effect of Seismic Ground Motion[J]. Technology for Earthquake Disaster Prevention, 2022, 17(3): 481-489. doi: 10.11899/zzfy20220308

场地对地震动放大作用的方向相关性研究

doi: 10.11899/zzfy20220308
基金项目: 中央高校基本科研业务费专项资金创新团队资助计划项目(ZY20160110);河北省自然科学基金(E2021512002)
详细信息
    作者简介:

    卢滔,男,生于1979年。博士,副教授,硕士生导师。主要从事工程地震理论和应用研究。E-mail:lutao@cidp.edu.cn

Study on Azimuth-related Site Amplification Effect of Seismic Ground Motion

  • 摘要: 现有场地对地震动的放大作用分析中,常采用观测得到正交方向的HVSR和SBSR的矢量合成结果表征场地放大作用,而不考虑放大作用的方向相关性。为研究不同地震作用下场地放大作用的方向相关性,采用多向HVSR和多向SBSR及谱比峰值方向相关分布概率,以11次地震事件获得的地表、井下基岩记录作为数据,开展响嘡场地效应台阵场地对地震动放大作用的方向相关性分析。研究结果表明,常用的矢量合成法分析结果与多向谱比上包络曲线法分析结果存在差异,前者存在忽略部分峰值的可能;多向谱比的极坐标等高线图可较直观地体现场地对地震动放大作用的方向相关性特征,但不同多向谱比表现存在差异;多向谱比峰值方向分布概率可较好地定量分析放大作用峰值在不同方向上的差异,但不同地震事件存在差异,不同的多向谱比存在差异。可采用多向谱比上包络曲线法分析场地放大作用,除可避免矢量合成法引起的误估外,也可更好地分析实际场地频率-方向相关的放大作用。
  • 近年来,随着中国城镇化进程的加快,高层建筑数量逐年增加,其中剪力墙结构凭借其整体性能好、刚度大、在水平作用下侧向变形小的优势,逐渐成为中国城镇的主要结构类型。中国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  多向HVSR极坐标等高线图和上包络曲线

    Figure  1.  Azimuth-related HVSR polar contour map and upper envelope

    图  2  多向HVSR峰值点分布散点图和方位角相关统计分布概率图

    Figure  2.  Peak value point distribution scatter of azimuth-related HVSR in 0-5 Hz range and the histogram of azimuth-related distribution probability

    图  3  根据弱震记录计算得到的HVSR和SBSR均值曲线

    Figure  3.  Mean HVSR and SBSR curve from different algorithms and different events

    图  4  HVSR均值、SBSR均值和理论传递函数及理论SBSR的对比

    Figure  4.  Comparison of the mean of HVSR, SBSR, theoretical transfer function, and SBSR

    图  5  全部地震记录多向HVSR和多向SBSR峰值点角度相关分布概率

    注:图中红色为HVSR数据,蓝色为SBSR,紫色为二者相同重叠部分

    Figure  5.  Peak value point distribution probability of azimuth-related HVSR and SBSR based on 11 earthquake events

    图  6  单次地震记录多向HVSR峰值点角度相关分布概率和全部记录分布概率对比

    注:彩色柱状图为单次地震记录结果,灰色为全部地震记录结果

    Figure  6.  Peak value point distribution probability of azimuth-related HVSR from single event and all events

    图  7  单次地震记录多向SBSR峰值点角度相关分布概率和全部记录分布概率对比

    注:彩色柱状图为单次地震记录结果,灰色为全部地震记录结果

    Figure  7.  Peak value point distribution probability of azimuth-related SBSR from single event and all events

    图  8  2种计算方法得到的典型地震记录HVSR和SBSR差异比较

    Figure  8.  Comparison of HVSR and SBSR by 2 different algorithms from the records in typical earthquake event

    图  9  典型地震记录HVSR和SBSR峰值点角度相关分布概率差异比较

    Figure  9.  Comparison of peak value point distribution probability of azimuth-related HVSR and SBSR from the records in typical earthquake event

    表  1  本文选用记录对应的地震事件相关参数

    Table  1.   Parameter of the earthquake events selected in this work

    发震时间/(年-月-日)震中位置震级ML震源深度/km震源关系PGA/Gal
    纬度/°经度/°地点震中距离/km方位角/°EW向NS向
    2001-07-11北纬39.75东经118.28陡河3.5536.4278.96.228.30
    2001-10-22北纬39.78东经118.70滦县3.18.9018.4111.58
    2001-12-28北纬39.67东经118.65滦县4.2125.4232.297.96161.71
    2002-01-30北纬39.68东经118.25唐山3.338.7266.93.104.75
    2002-02-09北纬39.77东经118.37陡河4.029.3285.513.6118.51
    2002-05-06北纬39.80东经118.77滦县3.812.628.448.8060.00
    2002-05-19北纬39.37东经117.98丰南4.771.9239.65.846.83
    2002-06-06北纬39.70东经118.75滦县2.764.390.05.806.90
    2002-06-20北纬39.88东经118.78卢龙2.8921.118.98.957.32
    2002-07-20北纬39.85东经118.75沙河驿3.31017.214.46.3012.40
    2002-11-13北纬39.77东经118.63滦县3.0159.8322.48.657.55
    下载: 导出CSV
  • 李红光, 冷崴, 2019. H/V谱比方法在强震记录应用中几个影响因素的初步研究. 地震工程学报, 41(1): 169—176 doi: 10.3969/j.issn.1000-0844.2019.01.169

    Li H. G. , Leng W. , 2019. Preliminary study on the influencing factors of the horizontal-to-vertical spectral ratio applied in strong motion. China Earthquake Engineering Journal, 41(1): 169—176. (in Chinese) doi: 10.3969/j.issn.1000-0844.2019.01.169
    卢滔, 周正华, 周雍年等, 2006. 关于Nakamura方法有效性的讨论. 地震工程与工程振动, 26(1): 43—48 doi: 10.3969/j.issn.1000-1301.2006.01.007

    Lu T. , Zhou Z. H. , Zhou Y. N. , et al. , 2006. Discussion on validation of Nakamura's technique. Earthquake Engineering and Engineering Vibration, 26(1): 43—48. (in Chinese) doi: 10.3969/j.issn.1000-1301.2006.01.007
    荣棉水, 李小军, 王振明等, 2016. HVSR方法用于地震作用下场地效应分析的适用性研究. 地球物理学报, 59(8): 2878—2891 doi: 10.6038/cjg20160814

    Rong M. S. , Li X. J. , Wang Z. M. , et al. , 2016. Applicability of HVSR in analysis of site-effects caused by earthquakes. Chinese Journal of Geophysics, 59(8): 2878—2891. (in Chinese) doi: 10.6038/cjg20160814
    姚鑫鑫, 任叶飞, 温瑞智等, 2019. 强震动记录H/V谱比法计算处理的若干关键环节. 震灾防御技术, 14(4): 719—730 doi: 10.11899/zzfy20190403

    Yao X. X. , Ren Y. F. , Wen R. Z. , et al. , 2019. Some technical notes on the data processing of the spectral ratio based on the strong-motion records. Technology for Earthquake Disaster Prevention, 14(4): 719—730. (in Chinese) doi: 10.11899/zzfy20190403
    张照鹏, 荣棉水, 卢滔等, 2019. HVSR谱比法应用于强震数据时基本假定合理性讨论. 地震研究, 42(4): 538—545 doi: 10.3969/j.issn.1000-0666.2019.04.012

    Zhang Z. P. , Rong M. S. , Lu T. , et al. , 2019. Discussion on rationality of basic assumptions when applying HVSR method to strong earthquake data. Journal of Seismological Research, 42(4): 538—545. (in Chinese) doi: 10.3969/j.issn.1000-0666.2019.04.012
    D'Amico V. , Picozzi M. , Baliva F. , et al. , 2008. Ambient noise measurements for preliminary site-effects characterization in the urban area of Florence, Italy. Bulletin of the Seismological Society of America, 98(3): 1373—1388. doi: 10.1785/0120070231
    Del Gaudioa V. , Zhao B. , Luo Y. H. , et al. , 2019. Seismic response of steep slopes inferred from ambient noise and accelerometer recordings: the case of Dadu River valley, China. Engineering Geology, 259: 105197. doi: 10.1016/j.enggeo.2019.105197
    Konno K. , Ohmachi T. , 1998. Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor. Bulletin of the Seismological Society of America, 88(1): 228—241. doi: 10.1785/BSSA0880010228
    Lee C. T. , Cheng C. T. , Liao C. W. , et al. , 2001. Site classification of Taiwan free-field strong-motion stations. Bulletin of the Seismological Society of America, 91(5): 1283—1297.
    Li X. J. , Li N. , Rong M. S. , et al. , 2022. Novel evaluation method for site effect on earthquake ground motion based on modified horizontal to vertical spectral ratio. Frontiers in Earth Science, 10: 938514. doi: 10.3389/feart.2022.938514
    Nakamura Y. , 1989. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Quarterly Reports of RTRI, 30(1): 25—33.
    Rong M. S. , Wang Z. M. , Woolery E. W. , et al. , 2016. Nonlinear site response from the strong ground-motion recordings in western China. Soil Dynamics and Earthquake Engineering, 82: 99—110. doi: 10.1016/j.soildyn.2015.12.001
    Wang Z. M. , Carpenter N. S. , Woolery E. W. , 2019. Horizontal-to-vertical spectral ratio of S-waves and SH-wave transfer functions at the vertical seismic and strong-motion arrays in the Central United States. Journal of Applied Geophysics, 162: 64—71. doi: 10.1016/j.jappgeo.2018.10.017
    Wen R. Z. , Ren Y. F. , Shi D. C. , 2011. Improved HVSR site classification method for free-field strong motion stations validated with Wenchuan aftershock recordings. Earthquake Engineering and Engineering Vibration, 10(3): 325—337. doi: 10.1007/s11803-011-0069-x
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  322
  • HTML全文浏览量:  150
  • PDF下载量:  49
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-06-23
  • 刊出日期:  2022-09-30

目录

/

返回文章
返回