基于ETM的水电站厂房抗震分析

　　摘要:耐震时程法是一种动力弹塑性抗震分析方法，该方法的关键在于耐震加速度时程曲线的合成。基于我国水工抗震设计规范的标准设计反应谱人工合成4条耐震加速度时程曲线，并以某实际工程的水电站厂房为例，分别对其进行7条天然地震动下的增量动力分析和4条耐震时程曲线输入下的非线性分析，对比研究了水电站厂房在不同地震强度下的最大层间位移角、上下游墙顶点最大相对位移，并对水电站厂房在设计地震作用下的典型节点的相对地面峰值位移与峰值加速度进行了对比。结果表明:耐震时程方法可有效地预测水电站厂房在不同等级地震作用下的地震响应，其结果处于增量动力分析结果的包络值之内，且在设计地震作用下，耐震时程方法在一定的容许范围之内可有效地对厂房结构作出地震响应评价，较增量动力分析，耐震时程方法避免了多次的调幅计算，大大提高了计算效率。

　　关键词:耐震时程分析;水电站厂房结构;增量动力分析;抗震分析

　　1研究背景

　　我国水资源丰富，但分布不均，主要集中于我国西南地区。近年来在我国水资源丰富的西南地区修建了一大批巨型水电站，如溪洛渡、锦屏、大岗山、白鹤滩水电站等。该地区属于高地震烈度区，水电工程选址难以避让，厂房抗震安全受到了严重威胁。因此开展对水电站厂房的抗震安全研究及评价具有重要意义。传统的抗震分析方法有:反应谱法、时程分析法、pushover分析方法和增量动力分析(IncrementalDynamicAnalysis，IDA)方法等。

　　这些方法优缺点明确，pushover方法是一种静力分析方法，该方法具有预测结构地震响应的直观性、简便性以及对结构输入地震动的无依赖性[1-4]，但该方法不能有效地考虑高阶振型及地震动特性的影响，对结构地震响应的结果会造成一定的误差。IDA方法以大量的时程分析结果作为结构预测地震的基础，该方法有效地考虑地震动作用的随机性并且能够较为准确地预测结构在不同强度等级地震作用下的地震响应[5]，但由于每条地震波需进行多次地震峰值加速度调幅，造成大型复杂结构的非线性计算，效率过低。Estekanchi等[6]、Riahi等[7]、Nozari等[8]将pushover方法与IDA方法相结合，提出了耐震时程法(EnduranceTimeMethod，ETM)，该方法考虑时程分析中的动力效应并且结合了pushover方法中的简便性，在有效预测结构地震响应的同时，大大提高了计算效率。

　　目前国内外对于ETM的应用主要集中于混凝土框架结构、混凝土重力坝以及钢结构领域[9-11]，对于水工领域内的水电站厂房结构还未涉及。由于水电站厂房规模巨大，结构复杂，地震作用对其影响较大，发生震损将严重影响水电站枢纽的运行并威胁站内工作人员的安全。因此有必要开展ETM对于水电站厂房结构的抗震分析的适用性研究。

　　本文结合ETM理论，基于我国水工抗震设计规范的标准设计反应谱，生成了4条用于水电厂房抗震的耐震时程加速度曲线(Endurancetimeaccelerograms，ETA)，综合考虑结构的材料，荷载等参数，建立了水电站厂房三维有限元模型，并对水电站厂房分别进行了IDA和ETM分析，对比研究了结构在不同地震动强度下的最大层间位移角、上下游墙顶点最大相对位移，并对比结构在设计地震(7度)作用下典型节点处相对地面的峰值位移和峰值加速度的地震响应，对比讨论了IDA与ETM分析下的结构加速度分布系数，以验证ETM在水电站抗震的适用性。对比结果表明，ETM在节约了时间成本，提高了计算效率的同时，可以有效地对水电站厂房结构地震响应作出评价。

　　2ETA合成

　　2.1ETA合成的基本原理

　　ETM以动力推覆进行结构抗震性能评估，该方法对结构输入一条地震动强度随持时增大的加速度时程曲线，以实现结构从线性到非线性，直至结构破坏失稳的全过程分析。ETM的关键在于合成ETA，要求该加速度曲线在某一时程内，其目标加速度反应谱和目标位移反应谱与时间呈线性增大的关系。Sa(T，t)=ttTargetSac(T)(1)Su(T，t)=ttTargetSac(T)T24π2(2)式中:tTarget为目标时间点;Sac(T)为预先指定的反应谱;Sa(T，t)和Su(T，t)分别为目标加速度反应谱与目标位移反应谱。在一定精度条件下耐震加速度时程曲线上每个点不能同时满足公式(1)和公式(2)，利用Matlab编程对该问题进行优化求解如公式(3):minF(ag)=∫0Tmax∫0tmax[Sa(T，t)-SaT(T，t)2+α[Su(T，t)-SuT(T，t)]{2}dtdT(3)

　　2.2基于我国水工抗震设计规范合成

　　ETA在合成ETA时，需要预先指定的反应谱，一般情况下取规范反应谱或地震动记录形成的反应谱[4]。本文采用我国水工抗震设计规范中的标准设计反应谱作为预先指定的反应谱Sac(T)，具体参数详见水工抗震设计规范[12]。水平设计峰值加速度ah的取值大小将影响ETA的强度大小，故ah的取值大小不受规范限制，由需要合成ETA的大小进行调整。

　　基于水工抗震规范的标准设计反应谱，取目标时间点tTarget=10s，水平峰值加速度代表值ah=0.3g，场地的特征周期Tg=0.20s，阻尼比ζ=0.05，合成了4条持时为30s的ETA。同时给出了每条ETA中0～10、0～20、0～30s的反应谱与目标规范谱的拟合情况，可以看出每条ETA均与规范谱具有良好的拟合。需要注意的是初始地震动的形状对ETA合成有着重要影响，本文采用SIMQKE软件人工合成初始地震动，使其只有上升段，没有下降段与平台段与ETA形状保持一致。

　　从ETA可以看出，每条ETA在频谱特性上具有一定的差异，并且ETA的强度是随时间线性增大而一般地震加速度时程的强度并非随时间线性增大[13]用一条ETA进行结构分析，难免会造成分析结果的离散性，所以本文采用4条ETA对水电站厂房结构进行计算分析。给出了4条ETA在0～10s的反应谱与规范谱的的相对误差，可以看出4条ETA的0～10s反应谱的中位值与规范谱的相对误差控制在20%以内，故采用4条ETA的中位值对结构进行计算分析时，其反应谱误差较一条ETA要低。

　　3有限元模型

　　本文以某实际工程的坝后式水电站厂房为研究对象，并对中间标准机组段建立基础-厂房三维有限元模型。坐标原点位于地基上游底部角点处，并规定顺水流向为X轴向下游为正;垂直水流向为Y轴向左岸为正;竖直方向为Z轴向上为正;沿厂房建基面分别向Z轴负向、上下游以及左右岸两侧延深1倍厂房高度从而建立结构的地基。

　　地基、厂房下部大体积混凝土结构以及上下游墙均采用C3D8单元离散，楼板、风罩、尾水管、钢蜗壳采用采用S4单元离散，梁、柱采用B33单元模拟，钢网架以及厂房上下游墙的配筋采用T3D2单元模拟。水电站厂房上下游墙顶部与顶部网架通过铰接相连，水轮发电机组、吊车、屋面荷载、流道内的动水压力等荷载在相应的位置以附加质量模拟。厂房两侧设为自由边界，忽略相邻机组段的相互影响。

　　基础底面全约束，四周法向约束，逐步积分法采用H.H.T法求解。厂房主体混凝土的力学参数如表1所示。混凝土材料采用ABAQUS有限元软件中的混凝土损伤塑性模型(CDP模型)，本文利用《混凝土结构设计规范》(GB50010-2010)规定的C25等级混凝土应力-应变关系，根据Najar损伤理论计算混凝土的拉伸压缩损伤因子，建立出C25等级混凝土的损伤塑性模型。考虑到厂房的上部结构是抗震关注的重点，所以仅对有限元模型的上下游墙进行了配筋，上下游墙内外两侧分别布置2层钢筋，竖向钢筋布置为Φ36@200，横向钢筋布置为Φ20@20。钢筋单元(T3D2)采用*EMBEDED方式嵌入到混凝土单元中，采用双线性随动强化本构模拟其力学行为。地基岩体及钢构件假定为线弹性模型。

　　4水电站厂房抗震分析

　　4.1地震波的选取

　　为对比分析ETM与IDA结构地震响应预测的差异，本文以水工抗震设计规范中的标准设计反应谱为目标谱，从地震数据库中心选取6条天然地震波，并基于水工抗震设计规范的标准设计反应谱人工合成一条地震波。给出了6条天然地震波调幅后与水工设计规范标准设计反应谱的对比情况，地震动反应谱的中位值与规范谱能够较好地吻合。将上述地震波记录进行调幅，每条地震波调整PGA为(0.1～1.0)g之间的地震波记录，每条地震波记录调幅10次，调幅增量为0.1g。

　　5ETM分析结果

　　5.1典型节点处的地震响应

　　本文中实际工程的水电站厂房结构是按照7度(0.1g)进行设计，对IDA结果取0.1g峰值加速度对应的结果，对于ETM分析结果同样取0.1g峰值加速度对应的结果。为更加全面地反映厂房结构地震响应结果，在厂房结构有限元模型中选取5个典型节点作为评价对比的依据，分别为节点8076(下游墙顶部中部)、节点6626(下游墙与发电机楼板连接处中部)、节点8235(发电机层楼板跨中某一点)、节点8068(风罩靠下游侧某一点)和节点3001(座环靠下游侧某一点)，分别读取其相对地面峰值加速度值和相对峰值位移值作为厂房结构的地震响应，以应证ETM在水电站厂房结构抗震中的应用。分别给出了厂房结构5个关键节点处的时程分析与ETM分析的相对地面峰值位移与峰值加速度响应，分别对时程分析结果与ETA分析结果求其响应的平均值，并对时程分析与ETA响应的平均值进行差值分析。

　　对于相对地面峰值位移与峰值加速度响应的相对差值，节点3001(座环靠下游侧某一点)分别为9.57%、17.81%，节点6626(发电机层楼板同下游墙连接处中部)分别为2.284%、17.42%，节点8076(下游墙顶中部)分别为6.202%、19.93%，节点8068(风罩靠下游侧某一点)分别为1.377%、18.61%，节点8235(发电机层楼板跨中某一点)分别为2.549%、17.45%。结果表明，在一定的误差容许范围内，ETM可以有效的预测水电站厂房结构地震响应。

　　5.2加速度分布系数

　　在水电站厂房的抗震设计中，水平地震加速度沿厂房高度的分布是设计过程中关注的重点[17]，为对比研究ETM在水电站厂房结构抗震分析的适用性，将时程分析与ETM分析的厂房结构加速度分布对比研究，给出了厂房结构在设计地震(7度)作用下时程分析与ETM分析下的水平地震加速度沿高度的变化，时程分析的中位值与ETM分析的中位值近似的接近。当水电站厂房的高度大于发电机楼板(33.80m)时，由于鞭梢效应的存在，使得加速度分布系数开始显著的增大。对于水平地震加速度沿厂房高度的分布，时程分析与ETM分析的变化趋近似相同，这说明ETM可以有效地运用在水电站厂房的抗震设计中。

　　6结论

　　本文基于我国水工抗震设计规范的标准设计反应谱，人工合成了4条持时为30s的ETA，并建立了某实际工程水电站厂房的三维有限元分析模型，分别对水电站厂房进行6条天然地震波加1条人工波的IDA以及4条ETA的ETM分析，得到以下结论:

　　(1)在基于水工设计规范生成ETA的过程当中，由于ETA的地震强度等级持续增大，所以设计水平地震动峰值ah可不按照规范选取，而是根据需要合成ETA的强度自行调整。如本文中以水平地震动峰值ah=0.3g，人工合成目标时间点tTarget=10s、持时30s的4条ETA曲线。(2)人工合成ETA曲线时，初始地震动的形状将会影响ETA的合成，严重时有可能造成ETA曲线的失真。本文在合成ETA曲线的过程中，采用SIMQKE人工合成初始地震动，通过控制参数使初始地震动的波形只有上升段，没有平台段和下降段，与ETA曲线的形状相似。

　　(3)对比IDA与ETM分析结果的最大层间位移角、上下游墙顶点最大相对位移，尽管由于输入地震动的随机性导致其结果出现了一定的离散性，但IDA与ETM分析结果非常接近，而ETM分析只进行了4次运算就可以得出与IDA70次运算相近似的结果，大大提高了计算的效率。(4)研究对比了在设计地震下(7度)作用下，时程分析与ETM分析的相对地面峰值位移与峰值加速度响应以及加速度分布系数，结果表明在一定的容许误差范围内ETM分析可以有效地对水电站厂房进行抗震响应评价。

　　(5)当水电厂房遭遇大震(9度)作用时，对比研究了时程分析与ETM分析的损伤状况，结果表明尽管由于输入地震动的持时和频谱特性对损伤结果造成了一定的误差，但作为一种简化的抗震分析方法是可以接受的。

　　综上所述，对于水电站厂房结构，ETM在预测不同强度等级地震作用下的地震响应时能与IDA的结果较为吻合，可以大大提高计算效率。对于研究水电站厂房的地震易损性[18-20]具有极大的潜力，但需要更近一步的研究。

　　参考文献:

　　[1]ATC-40，Seismicevaluationandretrofitofconcretebuildings[S].RedwoodCity(CA):AppliedTechnologyCouncil，1996.

　　[2]马千里，叶列平，陆新征，等.采用逐步增量弹塑性时程分析方法对RC框架结构推覆分析侧力模式的研究[J].建筑结构学报，2008，29(2):132-140

　　[3]吴京，梁仁杰，王春林，等.基于非线性静力分析的地震强度直接计算方法[J].建筑结构学报，2011，32(9):44-49.

　　[4]白久林，欧进萍.基于耐震时程法的钢筋混凝土框架结构抗震性能评估[J].工程力学，2016，33(10):86-96.

　　[5]VAMVATSIKOSD，CORNELLCA.Incrementaldynamicanalysis[J].EarthquakeEngineeringandStructuralDynamics，2002，29(10):2551-2562

　　[6]ESTEKANCHIHE，VALAMANESHV，VAFAIA.Applicationofendurancetimemethodinlinearseismicanalysis[J].EngineeringStructure，2007，29(10):2551-2562.

　　作者：朱少坤1，宋志强1，张大伟2

　　建筑论文投稿刊物：建筑结构学报宗旨在报道和交流建筑结构领域中代表我国学术水平的最新研究成果，反映本学科发展最新动态和趋势，推动国内外的学术交流，主要刊登建筑结构、抗震防振、地基基础等学科的基础理论研究、应用研究和科学实验技术的学术论文，研究报告及最新进展动态，为我国建筑科学技术研究的发展服务。