浅析河床式水电站厂房结构的地震响应分析方法研究

'第12卷第5期南水北调与水利科技Vo1．12No．52014年1O月South-to-NorthWaterTransfersandWaterScience8LTechnology0ct．2014DOI：10．13476／j．cnki．nsbdqk．2014．05．006河床式水电站厂房结构的地震响应分析方法于倩倩。，王海军，练继建(1．水电水利规划设计总院，北京100120；2．天津大学水利工程仿真与安全国家重点实验室，天津300072)摘要：以某河床式水电站工程为例，在分析厂房整体结构特点的基础上，建立了三维有限元模型，并应用《水工建筑物抗震设计规范》(DL5037-2000)和《欧洲结构抗震设计规范》中三种不同的反应谱，进行了河床式水电站厂房结构的地震响应分析。将设计反应谱转换成20S的人造地震波，通过时程分析，对反应谱分析得到的结构应力和位移状况进行验证对比，讨论在规范规定的不同反应谱分析方法下河床式水电站厂房结构地震响应的异同。期待研究成果可以为我国抗震规范的修订提供部分借鉴。关键词：河床式水电站厂房；地震响应；动力特性；反应谱法；时程分析中图分类号：TV31文献标志码：A文章编号：1672—1683(2014)05—0022—05Researchonseismicresponseanalysismethodoftherun-of-riverhydropowerstationYUQian-qian～．WANGHai-jun．LIANJi-jian(1．ChinaRenewableEnergyEngineeringInstitute，Beijing100120，China；2．StateKeyLaboratoryofHydraulicEngineeringSimulationandSafety，＆inUniversity，nnin300072，China)Abstract：Basedontheanalysisofstructuralcharacteristicsofarun-of-riverhydropowerstation，thefiniteelementanalysissoft—wareAnsyswasusedtObuildathree-dimensiona1structuralmodeltoanalyzetheseismicresponseoftherun-of-riverhydro—powerstation．ThreedifferentresponsespectrumsfromboththeEuropeancode，Eurocode8，andtheChineseCode，Specifica—tionsforSeismicDesignofHydraulicStructures(DI5037—2000)wereapplied．Thedynamicbehaviorsofthestructureunderearthquakeloadswereinvestigatedusingthetimehistoryanalysisduringthe20一secondsartificialseismicwavewhichwasex—tractedfromtheelasticresponsespectrum．Afteranalyzingandcomparingthestressesanddisplacementsofthepowerhousestructurewhichwereachievedfromdifferentcalculationtheoriesundertheseismicloadingcondition，themaincharacteristicsoftherun-of-riverhydropowerstationsduringearthquakeactivitieswereillustratedandsomesuitablesuggestionsforseismica—nalysiswereprovided．Keywords：run-of-riverhydropowerstation；seismicresponse；dynamiccharacteristic；spectrumanalysis；time-historyanalysis我国水能资源丰富，但地区分布很不均衡，其中67．8化存在突变。此外，厂房的上下游面均与水体直接接触，形位于西南地区。作为水利工程建设主要区域，西南地区处在成临空面，使得厂房结构在地震动力载荷作用下的响应问题欧亚地震带，受三大板块挤压，断裂带十分发育，地震灾害严变得更加复杂l3]。本文以某河床式水电站厂房结构为例，重。近年来，在汶川、玉树等破坏性大地震中，作为能源供应通过三维有限元方法，应用水工建筑物抗震设计规范(DI中心的水电站厂房结构遭受严重破坏，无法提供必要的电力5O37—2000)和欧洲结构抗震设计规范(Eurocode8)中三种保证，严重影响了救灾工作的顺利进行，甚至威胁到周边人不同类型的反应谱对结构的地震响应进行分析，并应用时程民的生命与财产安全[】]。作为一种常见的厂房结构型式，法进行结果对比验证和评价。河床式水电站厂房在大流量、低水头的河段能有效提高水资l河床式水电站厂房结构的自振特性研究源的利用率，增大发电量。但是这种类型的水电站厂房本身既是发电建筑物又是挡水建筑物，流道的空间所占的比重本文选取某水电站厂房一号机组坝段进行三维有限元大，呈现厚薄不匀的空间曲面形态，结构物的质量与刚度变建模，地基计算范围取厂房结构向上、下游各延伸60m，底收稿日期：2013—12—10修回日期：2014—07—07网络出版时间：2014—0827网络出版地址：http：／／www．cnki．net／kcms／doi／10．13476／j．cnki．nsbdqk．2014．05．006．html基金项目：国家自然科学基金创新研究群体科学基金项目(51021004)；国家自然科学基金青年科学基金(50909072)；教育部新教师基金(20090032120082)作者简介：于倩倩(1988一)，女，河北石家庄人，硕士研究生，主要从事水电站厂房结构工程方面的研究。E-mail：675215036@qq．com．22．哉H硷磺宄 于倩倩等·河床式水电站厂房结构的地震响应分析方法部延度不小于60in。坐标原点选在机组中心线处，规定y震规范以国内外各类地基上的多个地震加速度记录为基础，轴与机组中心线重合，向上为正；X轴为水平方向，顺水流方取多次地震的平均值作为代表值，绘制反应谱，经修正得到向为正；Z轴为水平垂直水流方向；三轴符合右手螺旋定则。一条具有代表性的、偏于安全的外包曲线，作为水电站厂房对于边界约束条件，取地基的下表面全约束，四周表面法向抗震设计依据。由于本次研究的河床式水电站厂房结构的约束，机组段之间因考虑分缝定为自由面。厂房结构三维模地震设计烈度为7度，水平向设计地震加速度代表值一型有限元模型见图l、图2。0．101g，设计反应谱见图3。设计反应谱最大值代表值为一2．25，最小值代表值Pm。一O．45，特征周期为—O．20S。根据工程的地震设计烈度和工程等级，此次计算只考虑顺水流方向的水平地震作用l6]。225q100图1厂房结构三维模型045Fig．13一Dmodelofthehydropowerstation010230r／s图3《水工建筑物抗震设计规范》设计反应谱曲线Fig．3ThedesignresponsespectrumintheDesignCodeforHydropowerHouse2．1．2《欧洲结构抗震设计规范))(Eurocode8)中的线弹性反应谱法图2厂房结构有限元模型与我国反应谱曲线不同，由于考虑了地震烈度对地基特Fig．2Finiteelementmodelofthehydropowerstation性和拐点频率有直接影响，《欧洲结构抗震设计规范》(Euro—通过有限元计算，该水电站厂房结构模型的前3O阶频code8)中的加速度响应谱由四段曲线组成。针对地震面波率E783及部分振型见表1。从表1中可以看出，结构的自振频震级Ms对结构物地震响应影响的差异，该规范共有两种推率值比较密集，相邻阶之间的频率值相差较小，而从第一阶荐反应谱，其中，Type1适用于震级小于5．5的情况，Type2到第三十阶结构的频率变化范围为1．42～11．83Hz。由振是适用于震级大于5．5的情况。本文根据工程设计烈度和型分析可知结构第一阶振型为上部排架结构的顺河向振动；地基基本状况，选择Type2的反应谱中GroundTypeA，谱第七阶出现厂房结构整体横河向振动，主要是进水口的振动；第十阶出现厂房整体顺河向振动；第十二阶出现厂房结中相应的特征参数为：土壤系数s一1；特征周期一0．05构整体扭振；第二十八阶出现尾水平台的横河向振动。可以S，Tc—O．25s；一l_2S，选择阻尼比为5A0，则r／=1。只考看出，河床式水电站厂房具有二元并联体系特征，即具有两虑顺水流方向水平地震作用，反应谱见图4。组自振特性：一组是主要反应厂房排架结构自振特性，一组主要反应大块混凝土结构的自振特性。计算结果符合一般r一25。河床式水电站厂房自振特性规律E33。，1＼0．625表1结构自振频率1+30T＼0．7510一Tab．1Naturalfrequencyofthestructureduringfreevibration阶数频率值／Hz阶数频率值／Hz阶数频率值／Hz卜1L42115．O9217．13005025123丁21．46125．8O227．2131．90135．82237．22图4《欧洲建筑抗震设计规范》(Eurocode8)41．92145．83247．87线弹性反应谱曲线52．1O156．23258．64Fig．4TypicalelasticresponsespectruminEurocode862．13166．3l268．942．1．3《欧洲建筑抗震设计规范》(Eurocode8)中的72．65176．48279．0883．79186．88281O．51设计反应谱法93．81196．9O291l_53反应谱理论中结构被假定为弹性状态，但通过对单自由104．122O7．123O11．83度体系的研究发现，结构的延性可以通过折减结构线弹性响应谱中的地震力来实现。对于自振周期较长的结构采用与2河床式水电站厂房结构的地震响应研究线弹性响应有相等的位移来实现，而自振周期短的结构则可2．1河床式厂房结构地震的反应谱分析采用与线弹性响应释放相同的能量来实现。根据以上原则，2．1．1《水工建筑物抗震设计规范》中的线弹性反利用系数可将结构线弹性反应谱折减为考虑延性等弹塑性应谱法因素的设计反应谱。但是由于我国的具体国情，对河床式水针对水电站厂房结构的特点，我国现行的水工建筑物抗电站厂房结构的地震响应情况研究较为缺乏，在谱分析过程g8试斌|交；·23· 第12卷总第74期·南水北调与水利科技·2014年第5期中并没有考虑结构本身的塑性特点，现行的《水工建筑物抗表2不同类型反应谱计算结果震设计规范》以及设计手册中也没有相关内容。综上，考虑Tab．2Seismicresponsesfromspectrumanalysis到水电站厂房结构的经济合理性和钢筋混凝土材料在地震荷载下的延性以及地震能量扩散特点，本文借鉴《欧洲结构抗震设计规范)(Eurocode8)中的规定，对河床式水电站厂房结构进行通过折减系数修正后的设计反应谱分析l9。根据规范，结构对地震中能量的散逸能力后，结构整体刚度减小，对地震力响应力小于线弹性分析过程中的响应力，而相应的位移会大幅度增加，所以通过性能系数(即折减系数)q对弹性分析进行折减。结构的性能系数q是结构在5的黏滞阻尼下对地震荷载的线弹性响应与设计地震力的比值，与结构材料、系统特点以及设计过程等方面有关，并且受到结构的动力稳定和结构尤其是结构连接处循环荷载作用下的疲劳损伤的限制。结合规范，由于混凝土结构能量散溢能力和塑性较差，属于低塑性类，q不大于1．5，且上部结构刚度和质量都远小于下部大体积混凝土结构，横河向结构注：1．表中对于欧洲标准设计反应谱法的位移计算结果是按照规范规定进行刚度远小于顺河向，所以不满足规范中的结构刚度质量分布修正后的结果；2．该表中应力结果为利用Ansys中应力均化后的计算结果。均匀、对称的规定，采用折减0．8后的性能系数q一1．2Ela~s]，证结构整体的安全，但是对于河床式水电站厂房上部结构，则设计反应谱见图5。根据线弹性反应谱所设计的结构依据地震工况下的计算结E乙果显得较为笨重。但是地震工况可能不为该结构计算的控2083‘Il制工况，结果需结合排架柱的纵向钢筋配置研究。0．521从位移结果分析，结构的最大位移矢量与顺河向位移均IO00＼0．625出现于上游排架柱顶端，但由于在计算过程中没有考虑屋架0．2结构，使整体刚度降低，位移值较大。三种反应谱法得到的iI1005025l20I763【J【l位移值接近，这是由于结构的非线性即塑性、延性在7度的地震荷载中没有完全体现，混凝土结构的绝大部分依然处在图5《欧洲建筑抗震设计规范))Eurocode8设计反应谱曲线Fig．5Typicaldesignresponsespectrum线弹性变化范围。所以对高地震烈度地区应用系数折减理论，需要考虑材料的非线性条件，如根据混凝土、钢筋混凝土(ductilitymodifiedspectrum)inEurocode8在受到动力荷载时的动态响应和损伤演化进行分析。2．1．4反应谱分析计算结果应用系数折减后的设计反应谱计算的位移结果更符合应用三种不同类型的反应谱对河床式水电站厂房结构结构的在地震工况下的实际响应特性，能更加快捷的对结构进行分析，通过静动叠加后，动应力所占的比例减小至原应的减振抗震、控制变形保证厂房结构在设计地震下正常运行力值的0．35倍，结构中的应力和位移主要是静荷载作用的做出判断，但该设计反应谱对结构延塑性和能量散溢能力的结果。对河床式水电站厂房结构而言，自振特性主要受低阶计算结果可靠性、精确性则需要通过时程分析结果进行对比模态影响，即刚度较小的排架结构处振型变化较为明显，使分析。得在地震过程中动应力影响以及绝大多数的结构破坏都产总体来说，通过三种反应谱的计算分析，该河床式厂房生在上部结构中。结构体系的应力、位移分布规律在三种反在地震荷载的作用下依然能满足应力和位移的基本要求，并应谱分析下类似，其中位移最大值在上游排架柱顶端，应力保证结构在设计地震荷载下处于安全状态，即该河床式水电最大值为排架结构刚度变化处，最大值出现的原因是地震波站厂房结构的设计合理。与结构的自振频率产生共振。三种反应谱的计算结果最大值见表2。2．2河床式厂房结构地震的时程分析从应力角度分析，三种反应谱计算的应力最大值出现在为对比河床式水电站厂房结构在地震工况下的反应谱上游排架柱中非常小的区域内(应力集中区)，原因是此处排分析结果，本文采用时程分析法，对结构进行地震工况分析。架柱的横截面积发生变化，造成刚度突变；另外由于没有考研究区缺乏地震实测记录，因此本文通过厂房设计反应谱曲虑上部屋架结构，使厂房结构的刚度减小，排架结构受弯后线拟合得到了人工地震波即采用加速度时程曲线。按照结出现了较为明显的二次效应，弯矩值随变形量增加而进一步构处基本地基参数条件，水平顺河向的人工拟合地震波见图加大，刚度变化处应力集中更加明显。将应力计算结果在6，假设该次地震共持续20S。Ansys中进行均化处理后，发现应用三种谱分析，结构的材从图7可以看出，地震波与标准反应谱基本一致，所以料未充分利用，即大多数的混凝土结构仍为线弹性状态。对可以用于对结构进行抗震分析。计算得到结构瑞利阻尼系于承担挡水并能起到雍高水位作用的河床式水电站厂房结数a=0．4516，口一0．0053。根据时程分析结果，结构的最大构，下部混凝土结构处在线弹性阶段能够有效防止开裂，保位移矢量和x向位移也在上游排架柱顶端达到最大值，分．24．试验碍}交 干倩倩等·河床式水电站厂房结构的地震响应分析方法于对结构的抗震情况做出快速的判断，但是该反应谱没有考虑到不同地震烈度对谱曲线的影响，所以对计算的精度有一定影响，尤其不能满足高地震烈度带的精细化分析需要。i盛其次，结构的有限元计算结果表明，三种反应谱计算得到的应力水平相当，其中，欧洲结构抗震设计规范中线弹性。iI飘反应谱法过于保守，水工建筑物抗震规范得到的应力结果与●f欧洲结构抗震设计规范中的设计反应谱结果相当。由于地震力折减，应用欧洲结构抗震设计规范设计反应谱分析得到图6顺河向人造地震波的位移更加真实，为更加快捷地对结构的裂缝和变形的提前Fig．6Artificialearthquakewaveinflowdirection控制提供了保证，适用于工程实际。此外，该方法中折减系数考虑结构的空间几何形式，有利于推广规范在不同形式结构体系中的使用。国内水电站厂房设计规范中尚未提到通m过系数折减修正反应谱法的概念，所以建议将此方法引入规呈范中，有助于对结构的非线性和能量散溢能力以及地震工况捌下的变形进行初步评估，使结构设计更为合理。但是，欧洲结构抗震设计规范主要针对的是工业与民用建筑，在水工建筑结构中应用时需要进行更细致的分析。综合以上所选用的地震响应分析方法可以看出，反应谱图7顺河向人造地震波与反应谱拟合Fig．7Verificationoftheartificialearthquakewave分析可同时考虑结构各频段振动的振幅最大值和频谱两个主要要素，计算方便，计算量小，可有效计算结构响应的最大别为2．891cm和2．657cm。相比反应谱结果，该位移结果值。而系数修正后的设计反应谱方法继承了线弹性反应谱与利用欧洲结构抗震设计规范中设计反应谱分析所得的计法的优点，较好地弥补了结构在线弹性反应谱方法中未考虑算值较为接近。的塑性状态以及结构中内力重新分布现象。通过对结构整体应力计算结果的分析可知，的最大值但是，反应谱实质依然是线性分析，计算精度取决于参出现在上游排架柱刚度变化处，最大值为3．51MPa；．值最与合并振型的数目，且其忽略了基础与土层之间的相互作大应力为5．82MPa，出现在排架结构上刚度变化处，结构配用、地震作用的随机性、以及时间持续性。而大量的震害调筋计算时需要重点注意，另外在排架结构根部区域和尾水管查表明，有些按反应谱理论设计的结构，在未超过设防烈度出口处根部区域也出现较多的，较大应力区；应力值最的地震中，也遭受到了严重的破坏，充分说明了持续时间要大值为1．83MPa，出现在排架柱与横向联结系的连接部位，素在设计中应该被考虑。时程分析方法能够记录结构响应另外在排架结构的顶端、副厂房和主厂房楼板也出现了部分的整个过程，即模拟结构在整个地震持续时间内各时刻的地接近于0．50MPa的应力区域。震响应，能成功的处理非线性相关问题，但计算量大，地震响根据计算结果可以发现，时程分析法计算的响应位移值应计算值较大依赖于地震波时程曲线的选取，且这种方法只与设计反应谱法计算的响应位移峰值在数值上很接近，相对反映结构在一条特定地震波作用下的性能，往往不具有普遍误差很小；时程方法计算得到的位移值大于按照反应谱法计性。所以，可将反应谱理论作为第一阶段抗震初步设计，然算得到的值，说明在持续性的地震动力作用下，结构的刚度后利用动态时程分析理论作为第二阶段精细抗震设计，结合会减小，位移值增大，所以对于水电站厂房结构的位移控制两者的优势来更为有效的为高震害地区的重要性较高的河需要引起关注。而时程分析法计算所得的应力值较小，进一床式厂房地震响应进行研究。步说明了如果按照反应谱理论进行结构的地震响应动力计算，采用的是对于结构任意周期的最大加速度，对结构的实3结论际应力是放大的。本文以某河床式水电站厂房结构为例，结合三维有限元2．3河床式厂房结构地震响应分析方法讨论分析软件Ansys对厂房结构的自振特性和地震作用下的动对比水工建筑物抗震设计规范(DL5037—2000)和欧洲力响应进行分析，认为河床式水电站厂房具有二元并联体系结构抗震设计规范(Eurocode8)中结构地震响应的反应谱分特征，且在三种不同反应谱计算中，该水电站厂房在地震荷析方法，可以发现以下情况。载的作用下均能较好的保证结构安全，即该河床式水电站厂首先，《欧洲结构抗震设计规范》在工程中的应用较复房结构的设计在地震状态下是合理的。通过时程分析发现，杂，需要根据区域地震烈度和场地地基特点确定三个特征周通过系数修正后的欧洲结构抗震设计规范设计反应谱法计期，保证结构在地震过程中谱反应的曲线段变化(即加速度、算得到的位移和应力值与时程分析法计算的响应峰值在数速度和位移恒定状态)；此外，将结构的阻尼作为变量(非恒值上很接近，相对误差很小，二者仿真结果具有一致性，即在定的5的阻尼)进行考虑，增加了计算量。而《水工建筑物本工程中，欧洲结构抗震设计规范设计反应谱法能较好地考抗震规范》只需要与场地地基类别和水工建筑物本身的类型虑到结构的延性和能量散溢能力，可有效节省计算资源。因有关的特征周期和最大值的代表值两个参数，使用简便，易此建议将此方法引入规范中，便于对结构地震响应进行有效诚l蠢i宽．． 第12卷总第74期·南水北调与水利科技·2014年第5期分析。[6]DI5037—2000，水工建筑物抗震设计规范[s]．(DI5037—2000，SpecificationsforSeismicDesignofHydraulicStructures参考文献(References)：rS]．(inChinese))[7]欧阳金惠，陈厚群，张超然，等．基于接触单元的三峡电站厂房Eli王海军，练继建，王日宣．水电站厂房结构地震响应非线性分析[J]．水电能源科学，2008(3)：88—91．(WANGHai-jun，LIANJi—振动分析[J]，水力发电学报，2008(5)：41—46．(OUYANGJin—hui，CHENHou—qun，ZHANGChao-ran。eta1．AnalysisonthejJan，WANGRi-xuan．Nonlinearanalysisofthepowerhouseinvibrationofthreegorgeshydropowerhouseonthebasisofseismicactivities，WaterResourcesandPower，2008(3)：8891．(inChinese))double!nodecontactelements[J]．JournalofHydroelectricEn—gineering，2008(5)：41—46．(inChinese))E22林皋，陈健云．混凝土大坝的抗震安全评价[J]．水利学报，2001(2)：8-15．(IINGao，CHENGJian-yun．Seismicsafetyevalua—E8]SL266—2001，水电站厂房设计规范[s]．(si266—2001，De—tionoflargeconcretedams[J]．JournalofHydraulicEngineer—signCodeforHydropowerHouseES]．(inChinese))ing，2001(2)：8-15．(inChinese))[9]EN1998—1：2004，Eurocode8：Designofstructuresforearth—quakeresistance-Part1：Generalrules，seismicactionsand[3]王海军，练继建，杨敏，等．河床式水电站厂房结构模型试验研究[J]．水利水电技术，2005(10)：100—102．(WANGHai-jun，rulesforbuildings[S]．LIANJi-jian，YANGMing，eta1．StructuralmodelexperimentL1O]ENI998—2：2004，Eurocode8：Designofstructuresforearth—researchofthepowerhouseinriverchannel[j]．WaterRe—quakeresistance-Part2：Bridges[S]．sourcesandHydropowerEngineering，2005(10)：100—102．(inLl1]EN1998—3：2005，Eurocode8：Designofstructuresforearth—quakeresistance—Part3：Assessmentandretrofittingofbuild—Chinese))ings[S]．[4]练继建，秦亮，何成连．基于原型观测的水电站厂房结构振动分析[J]．天津大学学报：自然科学版，2006(2)：176～180．(LIANL12]EN1998—5：2004，Eurocode8：Designofstructuresforearth—Ji-jian，QINLiang，HELian—cheng．Structuralvibrationofhy—quakeresistanc~Part5：Foundations，retainingstructuresanddropowerhousebasedonprototypeobservation[J]．Journalofgeotechnicalaspects[S]．TiainUniversity：Natural，2006(2)：176—180．(inChinese))Ll3]EN1990：2002+A1：2005，Eurocode0：Eurocode：Basisof[5]彭新民，郑伟，秦亮．水电站厂房结构振动规律探讨[J]．水利水structuraldesign[s]．电技术，2005(11)：62—64．(PENGXin—min，ZHENGWei，QINL14]EN1992—1—1：2004，Eurocode2：Eurocode2：DesignofconcretestructuresPart1-1：Generalrulesandrolesforbuildings[s]．Liang．Researchonthestructuralvibrationofpowerhouse~J]．WaterResourcesandHydropowerEngineering，2005(11)：62—L15]ENl992—3：2006，Eurocode2：Designofconcretestructures-64．(inChinese))Part3：Liquidretainingandcontainmentstructures[S~．版权转让声明本刊已加入万方数据数字化期刊群(www．wanfangdata．com．cn)、中国知网(www．cnki．net)、维普资讯网(dx1．cqvip．corn)和龙源期刊网(www．qikan．com)等网站，并被中国核心期-Dj(遴选)数据库、中国期刊全文数据库、美国《化{学文摘》(cA)等数据库收录。凡本刊录用的稿件将通过因特网进行网络出版或提供信息服务，稿件一经录用，将一次{{性支付作者著作权使用报酬(即包括印刷版、光盘版和网络版各种使用方式的报酬)，作者将该论文的复制权、发行权、j{信息网络传播权、汇编权等在全世界范围内转让给本刊。若有异议，请在投稿时作文字说明，编辑部将酌情处理。{；特此声明!i{《南水北调与水利科技》编辑部t；iI．k．．+一—_卜”+”++++．．+++．．+．．++．+．+++．．+++．．+++．+．++．+．．++．+．++．+．+．．+．．+．+．++++．．+．++．·26·试研尧'