地震动强度指标一直是众多学者研究的热点之一,刘恢先[1]认为选择地震动强度指标的方向主要有两个:一个是地震动地面运动指标,如PGA、PGV等;另一个是采取地震动的反应谱值,如结构第一自振周期对应的谱加速度。众多学者通过研究提出的地震动强度指标目前已有30种之多[2]。增量动力分析[3](Incremental Dynamic Analysis, IDA)是近年来受到众多学者关注的方法,该方法相比于静力弹塑性方法和动力弹塑性方法具有明显的优势。但因为其计算量大和研究积累的成果较少,所以,还未能成为抗震分析的主流方法,目前主要被广泛应用于结构基于性能的抗震评估中[4-6]。但随着计算机科学技术的飞速发展,IDA将是未来建筑结构弹塑性分析的重要发展方向。
目前众多学者在进行增量动力分析时,都选用PGA或者Sa(T1, ξ)作为地震动强度指标进行调幅,因为这两种指标概念清晰,且计算简单,满足了比例调整鲁棒性、充分性和地震危险性可计算性的要求,但在有效性方面,两种指标的适用性一直存在较大争议。对于以第一振型为主导的低层、多层结构,用Sa(T1, ξ)进行调幅计算得到的结果离散性较小[7],这一结论已获得了较一致的认同。但对于高层或者超高层建筑则尚未有一致的结论。Guan等[8]以框架核心筒为研究对象,对地震动强度指标进行研究,结果表明,PGV作为长周期地震动强度指标时效果较好,而PGA不适用于长周期结构。卢啸等[9]考虑到超高层结构高阶振型对结构有较大影响的特点,基于弯-剪耦合简化模型,提出适用于超高层结构的地震动强度指标Sa,并通过两栋高层结构的倒塌实例分析验证了该指标的相关性和离散程度。周颖等[10]对某型钢混凝土框架-钢筋混凝土核心筒超高层建筑进行了IDA分析,对比不同地震动参数IDA曲线的条件对数标准差平均值,得到对于高层建筑,考虑前3阶谱值参数的S123离散型最小、PGV次之、PGA离散型较大的结论。候喆辰[11]对高层剪力墙结构实际工程进行增量动力分析,对比分析结果的标准差,得到选用PGA作为地震动强度指标时计算结果的离散性比选用PGV要小的结论。郝永蕾[12]对某30层型钢混凝土框架-混凝土核心筒混合结构进行IDA分析,选择分位数曲线法和条件对数标准差对3个参数所得结果进行离散性对比分析,得到选用PGA作为地震动参数时离散性最小,PGV离散性稍大,Sa(T1,ξ)离散性较大的结论。由此可见,目前对于高层和超高层建筑的地震动强度指标并没有较一致的看法,笔者提出适用于超高层建筑结构增量动力分析的地震动强度指标Sa%,并通过实际工程证明其指标的有效性。
由于超高层建筑结构在设计时大多采用基于性能的抗震设计方法,具有较高的安全储备,在大震作用下,结构的主要抗侧力构件基本保持弹性或进入塑性的程度不高[13-14],非弹性反应谱是地震动特性和结构材料非线性相互耦合作用的结果,其中混合了结构的非线性影响的成分[10],故笔者暂不考虑采用非弹性谱作为地震动强度指标。基于目前地震动强度指标的研究现状和地震动强度指标在进行增量动力分析时需要满足的特性,纳入几个目前常用的地震动强度指标,并纳入概念简单、使用方便、且能够考虑高阶振型影响从而适用于超高层结构的地震动强度指标。选用已有的地震动强度指标如下:
1) 地震动的峰值加速度PGA。
2) 地震动的峰值速度PGV。
3) 阻尼比为4%时,结构第一自振周期对应的谱加速度Sa(T1, 4%)。
4) 阻尼比为4%时,地震动加速度谱峰值PSA。
5) 阻尼比为4%时,地震动速度谱峰值PSV。
6) 文献[3]中提出的考虑高阶振型参与的幂函数乘积形式的地震动强度指标。它以结构主要周期对应的弹性加速度反应谱值为底数,以各周期的振型参与质量系数比为指数的各振型结果的乘积。
当考虑前2阶振型时,其表达式为
式中: $\alpha=m_{1} /\left(m_{1}+m_{2}\right) ; \beta=m_{2} /\left(m_{1}+m_{2}\right) $。
当考虑前3阶振型时,其表达式为
式中:α=m1/(m1+m2+m3);β=m2/m1+m2+m3);γ=m3/(m1+m2+m3)。其中,m1、m2和m3分别为对应于前3阶自振周期T1、T2和T3的振型参与质量系数。
7) 文献[9]中提出的适用于超高层建筑的改进地震动强度Sa。该指标能够反映结构中高阶振型的影响,而且形式简便,便于应用。其表达式为
式中:Sa(Ti)为第i周期对应的谱加速度;n为所考虑的结构平动振型数,与结构的基本周期相关;Sa的物理含义是结构前n阶平动自振周期对应的谱加速度的几何平均值,反映了结构地震响应中参与显著的结构振型数量。
在第6)、7)地震动强度指标的基础上,提出适用于超高层结构的地震动强度指标Sa%。该指标物理含义为考虑结构前a%振型质量参与系数的自振周期对应谱加速度的几何平均值。
提出此地震动强度指标的理由如下:
1) S12和S123考虑高阶振型参与影响的能力有限。对于许多复杂高层结构,若考虑更高阶振型参与的指标形式S1234和S12345,会因为更高阶振型的质量参与系数占比不高使该指标数值变化很小,这会在一定程度上弱化高阶振型对结构的影响。但其依据质量参与系数来确定参与振型数量的方法取得了一定的效果,所以,Sa%也以平动方向的质量参与系数来确定振型的参与数量。
2) Sa考虑的结构前n阶振型主要由经验公式所确定,该经验公式由弯曲-剪切耦合的二维连续化模型进行时程分析得到,并非三维的超高层模型。且用两个超高层模型进行验证,模型A的验证效果较好,模型B中该指标的效果不如PGV。综上来看,该指标确定振型数量的经验公式存在局限性,因为该经验公式只与结构的第一自振周期有关。而即使是第一自振周期相近的超高层结构,高阶振型的参与程度也因受到诸多因素的影响而有较大差别。但该指标计算各自振周期对应谱加速度的几何平均值的计算方法相比于S12和S123能更有效地反映出结构高阶振型的参与。
在分析上述两种指标优点和不足的基础上,提出以平动方向前a%振型质量参与系数为依据,确定参与的自振周期数,并计算对应谱加速度的几何平均值。
笔者研究的内容是使用不同地震动强度指标进行IDA调幅计算,将计算结果得到的数据进行离散性分析,此项研究主要涉及两个问题:一是如何将由PGA调幅计算得到的数据转化成为各项地震动指标调幅计算下得到的结果,二是采用何种方法来评定数据的离散性。
由于IDA计算需要大量的非线性时程分析,算例结构复杂,每次时程分析耗时较长,且考虑的地震动强度指标较多,但地震动强度指标都具有比例调整鲁棒性,可以将PGA进行调幅后计算得到的结果转化成等效于用其他指标调幅后计算所得到的结果,PGA调幅值与各地震动指标下相对应的调幅值的具体转换步骤如下:1)计算各条地震动作用下结构各个地震动指标的强度值;2)计算出各条地震动作用下其他地震动强度指标与PGA的比值;3)由于各条地震动PGA的调幅值是按照50、100、200、300、400、600、800、1 000 cm/s2进行,所以,将与PGA对应的各个调幅值(50、100、200、300、400、600、800、1 000 cm/s2)相乘,最终得到该条地震动作用下PGA调幅值相对应的各个地震动指标的强度值。由上述内容可知,只要计算出每条地震动作用下各地震动指标与PGA的比值,将该比值与PGA相应调幅值相乘,就可以把PGA调幅后计算得到的响应结果相应地分配给其他地震动强度指标。
为了评定数据的离散型,采用分位数曲线法判断IDA计算数据的离散程度,然后再计算IDA计算数据对于地震动强度指标的条件对数标准差平均值来直观判断数据的离散程度。分位数曲线法假定每条IDA曲线均服从对数正态分布,然后通过对数据百分位数的判定来衡量数据所在的相对位置。计算方法为计算不同DM值相对应的不同IM值的均值μM和不同IM值的对数标准差δM,得到(DM,μM)、(DM,μM×exp(+δM))、(DM,μM×exp(-δM))3条曲线,分别为50%、84%、16%比例曲线,可知不同地震动强度指标的离散程度。为了更加直观地看出数据的离散程度,计算每个θmax和Vmax的分位数曲线离散程度的误差
式中:Δ84%DM、Δ16%DM分别为84%分位数曲线和16%分位数曲线相对于50%分位数曲线的误差;IM84%、IM50%、IM16%分别为每个DM值对应下84%、50%、16%分位数曲线上对应的IM值。计算出某地震动指标相对于θmax和Vmax的分位数曲线离散程度误差后,取平均值即可得到离散程度平均误差。
同时也计算出每个地震动指标下IM对于DM的条件对数标准差,并计算出条件对数标准差的平均值,平均值越小代表离散程度越小,计算数据的有效性越好。
深圳汇德大厦项目位于深圳北站附近,包含2#、4#裙房和1#、3#主楼。其中,1#主楼结构高度为248.75 m,地上总共58层,地下3层。项目所在地抗震设防烈度为7度,场地类别是Ⅱ类,设计地震分组为第一组,场地特征周期是0.35 s,设计基本地震加速度值为0.10g。1#主楼是钢框架-钢筋混凝土核心筒混合结构体系,在12~16层、24~28层、35~39层、46~50层设置斜柱实现外框架相对于核心筒斜切角式的缩进,使建筑表现出旋转向上的效果。
汇德大厦1#主楼第一转换区以下的外框架平面尺寸为45.1 m×45.1 m的正方形,核心筒尺寸为23.1 m×23.7 m的矩形,外框架与核心筒的边界并不平行,通过4个转换区外框斜切角的转换后,外框平面尺寸变为38.4 m×38.4 m,且外框架边与核心筒边平行,部分楼层结构平面布置如图 1所示。核心筒的厚度和平面布置变化较多,布置方式如图 2所示,混凝土强度等级与厚度变化如表 1所示。外框架采用钢管混凝土柱,主要钢材牌号为Q345B,斜柱转换层附近的部分柱牌号采用了Q390GJC。梁构件由于截面尺寸众多,不再列举。
综合考虑软件的计算精度和效率,最终选择PERFORM-3D对结构进行增量动力分析。对于梁柱单元的模拟,PERFORM-3D中可采用弦转角模型、塑性铰模型和塑性区模型来模拟梁柱单元。弦转角模型虽然计算简单,但由于其反弯点在跨中,所以只适用于对称的结构。塑性铰模型由弹性段和塑性铰组成,塑性铰集中于一点,无法模拟塑性区的长度,研究人员需对塑性铰的特点有较深入的了解。塑性区模型将构件分为弹性段和塑性段,可以模拟塑性区段的刚度退化,更加符合实际结构的破坏状态,PERFORM-3D中使用纤维单元对塑性区段进行模拟,模型将混凝土梁划分为10条混凝土纤维和2条钢筋纤维;将钢梁划分为12条纤维,腹板10条纤维,两翼缘各1条纤维;钢管混凝土柱划分为48条纤维,混凝土划分为36条纤维,钢管划分为12条纤维,参考文献[15]的建议,取截面高度的0.5倍作为塑性区的长度。对于剪力墙,剪力墙的宏观模型为等代杆系模型和纤维截面模型,对于高层建筑这类具有大量剪力墙构件的结构,使用此类宏观模型可以兼顾计算精度与效率,剪力墙构件采用shear wall宏观单元模拟,沿截面高度把混凝土和钢筋各分配8条纤维,沿墙肢高度每层划分一个单元。钢材选取非屈曲型钢材的材料本构,骨架曲线拟采用三折线本构来模拟钢筋与钢材材料。混凝土可以分为非约束混凝土、箍筋约束混凝土和钢管混凝土。一般混凝土构件中非箍筋加密区部分的混凝土为非约束混凝土,本构关系按照《混凝土结构设计规范》(GB 50010—2010)附录C中的混凝土本构关系确定;箍筋约束混凝土采用Kent-Park模型;由于目前没有非常权威的钢管混凝土本构关系,也采用Kent-Park模型来定义,相较于箍筋约束混凝土乘以1.1的强度提高系数。
为了保证PERFORM-3D软件建立结构弹塑性分析模型的准确性,根据项目的实际施工图,先使用YJK软件进行建模并施加相应的荷载,通过计算得到结构的配筋和模态,并根据施工图修改剪力墙等主要构件的配筋,为PERFORM-3D模型的建立与校核做好准备。算例中结构的YJK模型和PERFORM-3D模型如图 3所示。
PERFORM-3D模型与YJK模型的荷载信息、配筋信息、几何信息以及质量信息都应该相接近,PERFORM-3D模型在墙单元中设置了内嵌梁来模拟连梁和墙的刚接状态,故需要将两个模型的模态结果进行对比来验证PERFORM-3D模型是否与YJK模型存在较大偏差。两种软件所建模型的前3阶周期对比结果见表 2,从表 2可以看出,两个模型的自振周期十分接近,说明PERFORM-3D模型相比于YJK模型没有存在较大的偏差,此外,还对两个模型的能量误差以及多遇与罕遇地震作用下基底剪力进行校核,校核结果表明结构整体的各项信息在合理范围内,PERFORM-3D模型具有较高的准确性,可以用于结构的增量动力分析。
通过分析大量的强震观测记录,《抗规》规定设计地震动用弹性加速度反应谱表示,设计特征周期定义为“抗震设计用的地震影响系数曲线中,反映地震震级、震中距和场地类别等因素的下降段起始点对应的周期值”[16]。当选择多条地震波对结构进行抗震分析时,根据场地特征周期与地震动的特征周期接近这一原则进行选波,可以从平均意义上得到比较符合工程所在场地特性的地震波。
综合上述分析,所选地震波的特征周期应与场地特征周期相接近,同时,选择的地震波还要满足《高规》中有效时长的要求。为保证结构进行增量动力分析时的随机性,不选用来自不同测站但是同一次地震的记录。选取12条天然地震波并生成3条人工地震波进行结构的增量动力分析,所选地震波信息如表 3所示,各条地震波与规范的归一化加速度反应谱,即地震动力系数曲线如图 4所示。
增量动力分析(IDA)方法首先应确定地震动的强度指标,不同的地震动强度指标进行IDA计算的结果具有不同的离散性。选定地震动强度指标后就要进行结构性能参数的选取,选取的参数有结构最大基底剪力、结构最大层间位移角。在确定地震动的强度指标和结构性能参数后,就应确定地震动的调幅方法。
目前,PGA和Sa(T1,ξ) 是最常用的地震动强度指标,而Sa(T1,ξ)对于长周期超高层结构并不适用。由于PGA概念简单,获取地震动时程记录后就可以直接获得,中国规范目前选定PGA作为地震动强度指标,先采用PGA进行调幅计算,然后再对常用的地震动强度指标和新提出的Sa%进行调幅计算,Sa%是结构前a%振型质量参与系数的自振周期所对应谱加速度的几何平均值,其首先要计算各阶振型所对应的质量参与系数,从而得到Sa%所对应的振型,Sa%最终的计算结果如表 4所示。对于PGA的调幅,由于中国《抗规》中9度烈度区罕遇地震加速度时程的最大值为620 cm/s2,且有发生“特大地震”的可能,所以,IDA分析中PGA上限值取1 000 cm/s2,为了保证结构获得较精确的计算结果并具有高效的计算效率,PGA调幅按照50、100、200、300、400、600、800、1 000 cm/s2进行。PGA在400 cm/s2之前,考虑到结构由弹性逐渐进入弹塑性阶段,构件的刚度和强度均有不同程度的退化,如果步长选取过长,会导致这一阶段计算结果的准确度不高,所以,这一阶段调幅步长选取相对较小。PGA在400 cm/s2之后已属于罕遇地震,在结构设计基准期内发生的概率非常小,并且当发生这类大地震时,结构进入塑性的程度较高,每次弹塑性时程分析都需要大量的时间,综合考虑这两点因素,故适当增加了步长。对于其他地震动指标相对于PGA的调幅值,可按前述方法进行计算,在地震动Y1作用下,PGA与各个地震动强度指标相对应的调幅值如表 5所示。为使用分位数曲线法,保证每条地震动最终都能将θmax计算至0.01以上,Vmax都能计算至100 000 kN以上。
在PERFORM-3D模型中,瑞利阻尼与模态阻尼是叠加的,结构的阻尼采用统一的模态阻尼加较小的瑞利阻尼来模拟,模态阻尼比统一取0.04。每条地震动按规范建议采取双向输入,X向和Y向的地震波峰值取为1∶0.85。
主要研究结构X向的地震动响应,以下结构的响应(基底剪力、层间位移角)都是X方向。依据IDA计算结果,可以得到地震波在各地震动强度指标下的基底剪力和最大层间位移角的IDA曲线,限于篇幅,仅给出地震波在PGA下的基底剪力和最大层间位移角的IDA曲线,如图 5和图 6所示。根据IDA计算的结果和前述方法,可以得到各地震动强度指标对于θmax和Vmax的分位数曲线,限于篇幅,只给出PGA和PGV对于θmax的分位数曲线,如图 7和图 8所示。基于θmax计算得到的分位数曲线离散程度平均误差如表 6所示,基于Vmax计算得到的分位数曲线离散程度平均误差如表 7所示。计算各地震动强度指标相对于θmax和Vmax的条件对数标准差平均值如表 8所示。
从表 6、表 7和表 8中可以得到以下结论:
1) 在最常用的3个地震动强度指标PGA、PGV、Sa(T1, 4%)中,PGV得到的分位数曲线离散程度最小,PGA次之,Sa(T1, 4%)最大,这也与众多学者认为Sa(T1, ξ)只适用于第一振型反应为主导的低层、规则结构的结论保持一致。对于Vmax这一结构性能指标,PGV基于Vmax得到的分位数曲线离散程度相较于θmax有所增大,与PGA接近。PSA和PSV的概念虽然简单,但离散程度大于PGA和PGV,实用价值不高。
2) S12、S123、S12345得到的分位数曲线离散程度也较大,相较于Sa(T1, 4%)未能取得良好的效果,这是因为算例结构较复杂,高阶振型的质量参与系数很小,所以,从一定程度上弱化了高阶振型对结构的影响。S3(S70%)得到的分位数曲线离散程度小于S123(S70%)和S12345(S80%)也验证了这一点。
3) 随着所考虑质量参与系数的百分比不断增加,Sa%离散程度也不断减小。S90%得到的分位数曲线离散程度相比于S80%要小,这也说明,若以S90%作为地震动强度指标调幅进行IDA计算,得到的θmax和Vmax结果相对文中的其他地震动指标的有效性较好。
通过以上研究可以看出,提出的Sa%在考虑结构X方向平动前90%振型质量参与系数的自振周期对应的谱加速度的几何平均值S90%作为地震动强度指标调幅进行IDA计算时,可以有效减少θmax和Vmax数据的离散性,提高计算数据的有效性,从而适当减少时程分析的次数。低层和多层等较规则结构由于第一平动振型所占的质量参与系数较大,往往达到70%~80%甚至更多,所以,Sa(T1, ξ)等效于S70%或S80%,从而具有优于PGA的有效性,但对于超高层结构则不适用。考虑高阶振型参与的幂函数乘积形式的地震动强度指标S12、S123和S12345的效果也不理想。PGV提高数据有效性的效果仅次于S90%,其相比于PGA具有一定优势,也值得考虑使用。
在适用于超高层建筑结构的地震动强度指标S123和Sa的基础上,提出考虑高阶振型影响并适用于超高层结构的地震动强度指标Sa%;根据某实际超高层建筑结构的IDA计算结果,使用分位数曲线法,并以分位数离散程度平均误差和条件对数标准差均值最小为依据,验证提出的地震动强度指标S90%相对于其他地震动指标的有效性最好。得到的主要结论如下:
1) 对于最常用的几个地震动强度指标,PGV得到的分位数曲线离散程度最小,PGA次之,Sa(T1, 4%)、PSA和PSV的离散程度较大,不适用于高阶振型影响较大的复杂超高层结构。
2) 相比于现有的若干个常用的和适用于超高层结构的地震动强度指标,在已有研究的基础上提出的S90%在用于超高层混合结构的增量动力分析时,能更有效减少θmax和Vmax数据的离散性,从而提高计算数据的有效性。