结构思享【10】 | 基于EPA的地震动记录调整方法及应用实例

发表时间:2022-07-19 15:28

图片

图片

结构思享汇,将不定期分享RBS在技术领域的研究及设计成果,借此机会与同行们互相学习切磋,共同提高,以期促进建筑结构设计行业进一步发展。



基于EPA的地震动记录调整方法及应用实例


图片


此文发表于《建筑结构》杂志

本文作者:常磊 廖耘 王亚勇





摘要

分析了基于PGA与基于EPA两种地震动记录调整方法的物理含义、调整过程及存在的缺陷。提出一种与地震动记录特性相适应的分段EPA地震动调整方法,并应用于某588m高地标建筑的时程分析中,取得良好的效果。




图片


第一作者

常磊 CHANG LEI

RBS主任工程师

高级工程师 / 工学博士




1 前言


在对结构进行弹性和弹塑性时程分析时,如何调整地震动记录的能量使其与建筑物所在场地的设防烈度相符,是一个重要的问题。过往通常采取的调整方法,是将地震动的峰值加速度(即PGA)与规范中该场地地震动峰值加速度相比得到一个调整系数,再将整条地震动记录按此系数进行缩放(以下简称PGA调整方法)。

但需要注意的是,现行建筑抗震设计规范[1](以下简称抗规)中建议的时程分析用地震动峰值加速度,是基于设计反应谱(地震影响系数曲线α(T),下称规范谱)中的平台段值αmax,除以地震动加速度反应谱的平均放大系数2.25得到的(如7度多遇地震下为0.08x1000/2.25=35Gal)。这一参数表征的是统计意义上的多条地震动记录平均有效峰值加速度,是由加速度谱曲线换算而来(即EPA),与单条地震动记录加速度时程曲线上的峰值加速度,在概念上是有本质区别的。

因此,抗规建议的地震动记录调整方法,其实是基于EPA的调整方法,而过往通常采用PGA调整方法并不能保证将地震动能量调整到与设防烈度相符。以下将对这两种方法的调整结果差异进行对比,并对EPA调整方法进行改进,使其能够应用于实际工程。




2 基于PGA和基于EPA的地震动记录调整


2.1 PGAEPA的概念定义


对地震动时程a(t),其加速度反应谱记为Sa(T),则其峰值加速度PGA为地震动时程内的最大绝对加速度值,定义按式(1)。


图片

根据上式,地震动加速度时程的PGA有可能因少数高频脉冲的峰值很高而导致GPA很高,而这些加速度脉冲对长周期结构的地震响应的影响较小,因此PGA与地震动作用下结构的损伤水平并没有太高的相关性,所以多位研究者提出了地震动参数EPA,主要有基于第二加速度峰值定义的[2]、基于弹性体系结构反应幅值进行定义的[3],以及基于加速度时程能量进行定义的[4]。但目前应用较为广泛的EPA定义是5%阻尼比的单自由度体系在周期范围0.1~0.5s的谱加速度反应的平均值除以2.5[5](我国抗规中取2.25),即按式(2),式中”avg”表示取均值。

图片


2.2 基于PGA和EPA的地震动记录调整过程


以人工地震动记录L750-1和天然地震动记录L952为例,按“抗规”7度(0.1g)大震下峰值加速度220Gal,基于PGA的调整非常简单,直接将地震动等比放大到峰值为220Gal即可,PGA调整后的地震动加速度时程曲线如图1、图2所示。


图片

图片


基于EPA的地震动记录调整,则需用杜哈梅积分(Duhamel's integral)先求出上述地震动的加速度谱,如图3所示。


图片



如前所述,要计算EPA值,首先需要将加速度谱进行平滑处理,采用低通滤波器过滤高频成分,常用的有Gauss低通滤波器以及Butterworth低通滤波器等(可以考虑高阶),以下采用一种较为简单的跟随平均处理方法,取定数据点数n0,对于需要平滑处理的第i点反应谱数据,取i点前后各n0/2个数据值(均为平滑处理前数据值,数据点不足n0/2时按端点数据值补足),求和计算平均值。L750-1、L952分按该方法处理后的平滑反应谱与初始反应谱对比分别如图4、图5所示,EPA值分别为198Gal和187Gal。

由于EPA值低于220Gal,L750-1和L952地震动需分别再放大1.11倍和1.18倍才能满足要求,此时这两条地震动记录的PGA值将达到244Gal和260Gal。可以看到,上述方法虽然操作简单,但其对地震动能量的判断却存在明显的误差。


图片

图片





3 EPA调整方法的局限和改进


上述两条地震动记录,采用EPA调整方法效果较为理想,是由于这两条记录的加速度谱峰值段正好位于0.1s~Tg范围内,与式(2)所定义的取值范围吻合。

但不同强震加速度记录的特性差异很大,如下图的天然记录L0055和L2572(初设地震动峰值均为220Gal),加速度时程分别如图6、图7所示,其加速度谱峰值段就不在0.1s~Tg区间。如果还是按照0.1s~Tg区间的加速度谱平均值去确定EPA调整系数,则两条地震动需分别放大1.83倍和1.94倍之多,PGA达到403Gal和427Gal,已经接近8度罕遇地震的水平,这显然是不合理的,如图8所示。


图片

图片

图片


研究表明,不同的强震加速度记录所包含的频率成分不同,有的呈现高频显著,有的呈现中、低频显著,从而导致由加速度反应谱推算出来的拟速度谱和位移谱差别很大。根据这一现象,可将加速度反应谱区分为“加速度型”、“速度型”和“位移型”三类。出现这一问题的原因,在于式(2)的EPA调整方法,只针对加速度谱峰值在0.1s~Tg区间的加速度型地震动,而不适用于速度型和位移型地震动。图9给出了各条地震动记录的速度谱和位移谱曲线与规范谱的对比,可以看到,L0055在Tg~5Tg区间,速度谱曲线明显超出规范谱,是典型的速度型地震动记录;而L2572在5Tg之后的长周期段速度谱与位移谱均超出规范谱,是位移型地震动记录。因此,上述EPA的调整方法,只适用于加速度谱峰值在0.1s~Tg区间的加速度型地震动,而不适用于速度型和位移型地震动。


图片


参考ASCE的做法[6],可以对基于EPA的调整方法进行改进,根据地震动记录的类型分段计算EPA并调整。对于加速度型、速度型和位移型地震动,其EPA分别取0.1s~Tg、Tg~5Tg、5Tg~6s区间的加速度谱平均值或最大值,相应的EPA目标值也根据规范谱在相同区间的最大值α(Tg)来计算。对人工地震动、加速度型天然地震动、速度型天然地震动以及位移型天然地震动,其EPA值和目标值分别按式(6)~式(9)计算。其中,考虑天然地震动记录加速度反应谱的波动性比人工模拟地震动更大,其EPA计算时取平台段区间平滑谱的最大值除以1.414,对应的计算式中分母为1.414*2.25 = 3.18。


图片


对速度型地震动,分析表明,上述分段EPA调整方法有时仍然会导致调整系数过高或过低。考虑到规范谱在Tg~5Tg区间变化较大,参照ASCE[6]和FEMA[7]的方法,当按上述方法调整后结构弹性时程基底剪力与CQC值差异仍超过35%时,其EPA计算值还需予以修正,修正系数取为Tg~5Tg周期段地震动平滑谱均值与规范谱均值的比值,按式(10)计算。



图片



图片


速度型地震动L0055和位移型地震动L2572,按上式(8)和式(9)调整之后,应分别乘以调整系数1.14和0.65,调整后PGA值分别为250.8Gal和143Gal。调整后的平滑加速度谱对比如图12所示,已没有此前过分放大地震动的问题,且在速度段和位移段也仍然保留了地震动记录的原有特性。




4 应用实例


某588m高地标超高层建筑,结构形式为带巨柱和加强层的框架核心筒结构,设置7道环桁架层,其中4处设有伸臂,结构外观以及核心筒和加强层立面如图13所示。计算的基本周期为T=9.4s,为超长周期结构,对具有长周期脉冲的位移型地震动记录的响应必然较大。

该项目按规范要求选择了2组人工地震动记录和5组天然地震动记录,分别是:L750-1、L750-4;L0055、L952、L2572、LMEX002及LMEX026,其中人工地震动记录和天然地震动记录L952为加速度型,L0055、LMEX002为速度型,L2572、LMEX026为位移型,初始PGA均调整至220Gal。

对各组地震动的主分量,按PGA=220Gal直接计算,和基于上述分段EPA方法调整后,计算得到的结构基底剪力如表1所示。


图片

图片

可以看到,采用传统的PGA调整方法,由于受到个别脉冲加速度的影响,加速度型和速度型地震动计算的基底剪力往往偏小,而位移型地震动的基底剪力则偏大很多,7组中有3组都不能采用。而采用分段EPA方法调整后,各组基底剪力均与CQC结果更为接近,且两组长周期的位移型地震动都能得到保留,这对于考察588m高建筑的地震破坏情况是非常重要的。


5 结语


从地震工程学的角度来看,峰值加速度PGA对地震动记录的能量反映是非常片面的,因此过往常用的PGA调整方法虽然简单,但调整后的地震动能量与设防烈度存在明显偏差;

基于加速度谱得到的有效峰值加速度EPA,能更好地反映地震动能量,但传统的EPA调整方法,只适用于短周期的加速度型地震动,不能满足实际工程的需要。

本文提出的按地震动记录的类型,区分加速度型、速度型和位移型分段进行EPA调整的方法,适用于各种不同周期的建筑,可以较准确地将地震动能量调整到与规范谱接近,使得时程分析中能保留更多有特色的地震动记录,值得在更多的案例中研究推广。

但上述分段EPA调整方法暂时还仅适用于水平地震动,对于竖向地震动,考虑其目前研究尚还不足[8],仍然建议按PGA方法调整。


参   考   文   献


[1] GB50011-2010 建筑抗震设计规范[S]. 北京:中国建筑工业出版社,2010.

[2] Gupta I. D. Defining effective peak acceleration via order statistics of acceleration peaks[J]. European Earthquake Engineering. 1994, 2: 3-11.

[3] Sabetta F., Bommer J. J. Modification of the spectral shapes and subsoil conditions in Eurocode 8[C]. Proceedings of the Twelfth European Conference on Earthquake Engineering. London: 2002. Paper No. 518.

[4] Sarma S. K., Yang K. S. An evaluation of strong motion records and a new parameter A95[J]. Earthquake Engineering & Structural Dynamics. 1987, 15(1): 119-132.

[5] Council Applied Technology. Tentative provisions for the development of seismic regulations for buildings[R]. Redwood City, California, 1978.

[6] ASCE STANDARD ASCE/SEI 7-10 Minimum Design Loads for Buildings and Other Structures[S]. Virginia:American Society of Civil Engineers,2010.

[7] NEHRP Recommended Seismic Provisions for New Buildings and Other Structures   (FEMAP-1050-1) [S]. Washington D.C.:BUILDING SEISMIC SAFETY COUNCIL,2015.

[8] 王亚勇. 有效峰值加速度EPA和地面峰值加速度PGA的异化及后果——兼述设计地震动参数αmax, β,PGA, Tg的定义与关系[C]// 《建筑结构》杂志社. 第二届复杂建筑结构弹塑性分析技术交流会报告专家PPT汇编. 苏州: 《建筑结构》杂志社, 2016: 1–4.







END