Study on Azimuth-related Site Amplification Effect of Seismic Ground Motion
-
摘要: 现有场地对地震动的放大作用分析中,常采用观测得到正交方向的HVSR和SBSR的矢量合成结果表征场地放大作用,而不考虑放大作用的方向相关性。为研究不同地震作用下场地放大作用的方向相关性,采用多向HVSR和多向SBSR及谱比峰值方向相关分布概率,以11次地震事件获得的地表、井下基岩记录作为数据,开展响嘡场地效应台阵场地对地震动放大作用的方向相关性分析。研究结果表明,常用的矢量合成法分析结果与多向谱比上包络曲线法分析结果存在差异,前者存在忽略部分峰值的可能;多向谱比的极坐标等高线图可较直观地体现场地对地震动放大作用的方向相关性特征,但不同多向谱比表现存在差异;多向谱比峰值方向分布概率可较好地定量分析放大作用峰值在不同方向上的差异,但不同地震事件存在差异,不同的多向谱比存在差异。可采用多向谱比上包络曲线法分析场地放大作用,除可避免矢量合成法引起的误估外,也可更好地分析实际场地频率-方向相关的放大作用。
-
关键词:
- 多向HVSR /
- 多向SBSR /
- 谱比包络 /
- 场地效应 /
- 峰值方向相关分布概率
Abstract: In recent analysis on site amplification effect of seismic ground motion, the square root of HVSR and SBSR obtained from seismic observation in orthogonal direction was often used to characterize the site amplification effect, however, the azimuth-related site amplification effect of seismic ground motion is not considered. To discuss this effect, the azimuth-related HVSRs, SBSRs, and distribution probability of peak values on the spectral ratios were selected, and the surface and borehole seismic records on Xiangtang Array site in 11 earthquake events were collected, finally, the azimuth-related site amplification effect was analyzed. The results show: (1) the results from the square root algorithm have the risk of ignoring some frequencies which are different with the proposed upper envelope algorithm in this work, (2) the contour of azimuth-related spectral ratios in polar coordinates could directly represent the directionality characteristics of site amplification effect, but the results from different multi-directional spectral ratios are different, (3) the azimuth-related of peak values could characterize in quantity the difference in azimuth of site amplification effect, but the probabilities were various with the different multi-directional spectral ratio and the data from different event. According to the results and discussion in this work, the azimuth-related spectral ratios could be used in the analysis of site amplification effect, for avoiding the under-estimate of amplification effect by the square root algorithm, and the method could be used to understand the frequency-direction correlation of site amplification effect better. -
引言
近年来,随着中国城镇化进程的加快,高层建筑数量逐年增加,其中剪力墙结构凭借其整体性能好、刚度大、在水平作用下侧向变形小的优势,逐渐成为中国城镇的主要结构类型。中国80%的城镇位于Ⅶ度地震设防区,一旦发生破坏性地震,将造成巨大的经济损失和人员伤亡。为了评估现有高层钢筋混凝土结构的抗震能力,世界地震工程界进行了大量的研究,其中增量动力分析法是1种有效的分析方法。
增量动力分析方法(Incremental Dynamic Analysis,IDA)考虑了结构抗震能力和需求的不确定性和随机性因素,能较合理地评估不同性能水准下结构抗震能力。并且,通过IDA法可以分析结构从弹性到弹塑性、逐渐损伤直至破坏倒塌的全过程,进而可进一步研究控制结构破坏程度的条件和寻找防止结构倒塌的有效措施(吕西林等,2012)。
2000年,美国开创性地在FEMA-350(Federal Emergency Management Agency,2000)中引入了IDA的概念,以此作为未来新建建筑结构抗震设计的标准。Vamvatsikos等(2002)在前人的研究的基础上,总结并提出了增量动力分析方法,创造性地提出了hunt & fill点跟踪算法,为实际工程中增量动力分析方法的应用创造了先例。21世纪初以来,该方法已经被广泛用于结构抗震性能评估。
对于一般框架结构而言,影响其抗震能力的主要因素有柱轴压比、梁柱配筋率、梁剪跨比等,而对于高层剪力墙结构来说,由于主要承重构件的变化导致影响其抗震能力的主要因素为抗震设防等级、剪力墙轴压比、剪力墙厚度、结构高度、剪力墙配筋率等,其中在结构高度以及设防等级对剪力墙结构抗震能力的影响方面研究甚少。因此,本文在增量动力分析的基础上,讨论了抗震设防等级和高度对剪力墙结构抗震能力的影响,以期为研究工程结构的地震破坏和损失提供一定依据。
1. IDA分析方法的基本理论
1.1 IDA分析方法的定义
增量动力分析方法将1个或多个地震动记录输入到结构模型中,然后对地震动参数(如PGV、PGA等)按照一定的规则进行调幅。通过这种方式,获得了1组不同强度的地震动记录,继而在这组记录的激励下,可以得到该结构的弹塑性分析结果。根据分析结果,可得到地震动强度指标IM(Intensity Measures)和结构损伤指标DM(Damage Measures)之间的关系曲线,即IDA曲线。通过该曲线,可以研究不同强度地震作用下结构损伤的全过程及相关失效概率。
1.2 IDA分析方法的步骤
对结构进行IDA分析的主要流程和步骤为(卜一等,2009;汪梦甫等,2010):
(1)严格按照抗震设计规范建立弹塑性模型;
(2)按照一定的规则选取有代表性的地震动记录(至少10条);
(3)将选取的地震动记录按照既定的调幅方法进行调幅并输入到结构中进行计算分析;
(4)将得到的IM、DM点连接成线,即为结构的IDA曲线。
1.3 地震动记录的选取
考虑到断层距、结构周期以及场地类型等因素的影响,谢礼立等(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 在表 1的地震动记录中,Ⅰ类场地的有F2、N1;Ⅱ类场地的有F3、F4、F5;Ⅲ类场地的有F7、F12;Ⅳ类场地的有F10、F11、N4。
远场地震动记录的断层距大于20km,根据表 1中的信息可知属于远场地震动记录的有F2、F3、F5、F10,其余的均为近场地震动记录,由此可见选取的地震动记录充分考虑了断层距、场地类型等因素的影响。
2. 典型剪力墙结构IDA分析
2.1 典型剪力墙结构结构设计
为了使结构分析结果具有普遍意义,有必要选择典型的结构,因此本文选择了剪力墙结构比较普遍的住宅楼,如图 1所示。
选取的工程为上海金桥地铁保障房,钢筋混凝土剪力墙结构体系,该建筑结构地上31层,总高度86.8m,层高2.8m,高宽比为5.36,安全等级为二级,丙类建筑,Ⅶ度设防,基本地震加速度值为0.10g,设计地震分组为第1组,Ⅳ类场地,特征周期为0.9s。
为达到研究目的,本文以图 1为模板,严格按照中国规范规定分别建立了Ⅵ、Ⅶ、Ⅷ度3种不同抗震设防,15层、20层、25层、30层4种不同房屋高度的12个模型(中华人民共和国住房和城乡建设部,2010;中华人民共和国建设部,2010;中华人民共和国住房和城乡建设部,中华人民共和国国家质量监督检验检疫总局,2010;中华人民共和国住房和城乡建设部,2012),以此来研究不同抗震设防等级以及不同结构高度对建筑结构抗震能力的影响。
2.2 增量动力分析
本文选择Etabs进行结构的弹塑性分析。梁、柱采用三维框架单元模拟,剪力墙和楼板采用三维壳体单元模拟,为了反映混凝土材料在地震荷载作用下非线性阶段刚度退化的影响,其恢复力模型采用Takeda模型,而基于钢筋材料的特性,恢复力模型则采用考虑包辛格效应的Kinematic模型,强屈比取值为1.2,极限塑性应变取0.025。在进行弹塑性分析时采用显式积分法。
逐一对所选择的地震动进行调幅,调幅后的PGA分别为0.005g、0.205g、0.455g、0.755g、1.105g和1.505g,地震动强度按X、Y向1:0.85输入,并按照文献(韩建平等,2007)的判断准则确定继续调幅或是终止分析。进而,可得到结构性能参数(最大层间位移角θmax)与地震动强度(地面峰值加速度PGA)之间的关系曲线,即IDA曲线。图 2绘出了12个不同模型的IDA曲线,图中3条IDA曲线分别为50%、16%、84%分位的IDA曲线。从图中可以看出,16%、84%分位的IDA曲线与结构的均值曲线,即50%分位的IDA曲线偏差并不大。
2.3 结构破坏等级与性能指标
工程结构的破坏可分为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 2.4 各性态点的确定
根据上述破坏状态定义原则,对各模型的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 2.5 抗震设防等级对剪力墙结构抗震能力的影响
根据表 3的数据,可以定量地分析出结构随抗震设防等级的提高在各个性能点处抗震能力提高的程度,为了更直观的表现,将以上数据绘制在图 3中。
由图 3可知,无论结构高度如何,随着抗震设防等级的增加,各破坏状态点处的地震指标(PGA)都相应的增加,虽然轻微破坏和中等破坏所对应的PGA值有所增加,但是增加的幅度较小,而严重破坏和倒塌所对应的PGA值却大幅度增加。由此可见,相同地震动强度下,不同抗震设防等级的结构发生轻微破坏、中等破坏的概率变化并不大,但是抗震设防等级越低,发生严重破坏和倒塌的概率越大。
2.6 结构高度对剪力墙结构抗震能力的影响
根据表 3的数据,还可定量地分析出结构随高度的增加在各个性能点处抗震能力降低的程度,为了更直观地表现,将以上数据绘制在图 4中。
将结构层数信息用高度来表示,见表 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 由表 4和图 4可以看出,Ⅵ、Ⅶ和Ⅷ度设防的结构,随着结构高度的增加,各破坏状态点处的地震指标PGA值都相应的降低,虽然轻微破坏、中等破坏和严重破坏所对应的PGA值有所降低,但是降低的幅度较小,而倒塌性能点所对应的PGA值却大幅度降低。由此可见,相同地震动强度下,不同结构高度的结构发生轻微破坏、中等破坏和严重破坏的概率变化并不大,但是结构高度的增加,会导致发生倒塌的概率显著增加。
3. 结语
本文以Ⅵ、Ⅶ和Ⅷ度3种不同抗震设防等级和15层、20层、25层、30层4种不同结构高度共计12个剪力墙结构模型为基本算例,通过IDA对比分析,对不同抗震设防等级和不同结构高度的剪力墙结构的抗震性能进行了评估,总结出以下结论:
(1)在对比分析不同抗震设防等级下剪力墙结构模型时发现,随着结构抗震设防等级的增加,结构的抗震能力随之增强,虽然在较小最大层间位移角处抗震设防等级的增加,并未导致PGA显著提高,但随着最大层间位移角的增大,抗震设防等级的提高所带来的抗震能力会有显著的提高。即抗震设防等级越高,结构在刚度、强度及延性等方面的表现越好,在遭受同等强度的地震动作用下结构也就越难发生轻微破坏、中等破坏和严重破坏。抗震设防等级的增大对增强结构抗震能力的效果很明显。
(2)在对比分析不同结构高度下剪力墙结构模型时发现,随着结构高度的增加,结构的抗震能力随之降低,结构高度对抗震能力的影响并没有抗震设防等级对结构抗震能力的影响明显。但随着最大层间位移角的增加,较高的房屋其地震危险性还是比较大的。
-
表 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.5 5 36.4 278.9 6.22 8.30 2001-10-22 北纬39.78 东经118.70 滦县 3.1 — 8.9 0 18.41 11.58 2001-12-28 北纬39.67 东经118.65 滦县 4.2 12 5.4 232.2 97.96 161.71 2002-01-30 北纬39.68 东经118.25 唐山 3.3 — 38.7 266.9 3.10 4.75 2002-02-09 北纬39.77 东经118.37 陡河 4.0 — 29.3 285.5 13.61 18.51 2002-05-06 北纬39.80 东经118.77 滦县 3.8 — 12.6 28.4 48.80 60.00 2002-05-19 北纬39.37 东经117.98 丰南 4.7 — 71.9 239.6 5.84 6.83 2002-06-06 北纬39.70 东经118.75 滦县 2.7 6 4.3 90.0 5.80 6.90 2002-06-20 北纬39.88 东经118.78 卢龙 2.8 9 21.1 18.9 8.95 7.32 2002-07-20 北纬39.85 东经118.75 沙河驿 3.3 10 17.2 14.4 6.30 12.40 2002-11-13 北纬39.77 东经118.63 滦县 3.0 15 9.8 322.4 8.65 7.55 -
李红光, 冷崴, 2019. H/V谱比方法在强震记录应用中几个影响因素的初步研究. 地震工程学报, 41(1): 169—176 doi: 10.3969/j.issn.1000-0844.2019.01.169Li 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.007Lu 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/cjg20160814Rong 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/zzfy20190403Yao 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.012Zhang 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 -