地震下岸坡坝段侧向抗滑稳定研究范文

摘要:岸坡坝段的抗滑稳定性关乎大坝结构运行安全。为研究地震作用下岸坡坝段的侧向抗滑稳定性，结合某工程采用ANSYS有限元软件建立岸坡坝段的三维有限元模型，分别采用拟静力法和反应谱法进行对比分析，结果表明:二者计算所得岸坡坝段地基面最大应力均小于地基允许承载能力4．0MPa;稳定计算均满足抗滑稳定要求，反应谱法计算结果较拟静力法偏大。该对比结果可为实际工程设计和安全评估提供重要参考。

关键词:地震;岸坡坝段;抗滑稳定;拟静力法;反应谱法

前言

近年来中国水电事业高速发展，在大坝建设方面取得了显著的成效，在大坝建设中遇到的许多问题都得到合理的解决。就重力坝而言，岸坡坝段稳定对整个坝体的稳定影响较大，对于岸坡坝段稳定方面的问题研究的比较多，其中李永明［1］分别用刚体极限平衡法和有限元法2种方法，并且按分项系数法的原理对某重力坝岸坡坝段侧向稳定进行了分析研究;丁旭［2］等人对蟒河口水库岸坡坝段进行分区，并利用刚体极限平衡法，按可能失稳情况进行稳定应力计算，最后对结果进行了分析;宫照光［3］等人分别采用有限元和材料力学2种方法，对某重力坝岸坡坝段上游面裂缝对岸坡坝段稳定的影响进行了分析。地震因素是不可控制因素，对岸坡坝段的影响较大，容易造成岸坡坝段失稳。在地震作用下岸坡坝段动力稳定性的问题研究的比较少，王超［4］等人采用概率思想，并结合确定性矢量和安全系数方法，对设计基准期内随机地震条件下的重力坝岸坡坝段动力抗滑稳定性问题进行了初步探讨;邓良军［5］等人采用数值计算法和模型实验法对金安桥水电站岸坡坝段在地震作用下进行了动力分析。而对于地震作用下岸坡坝段侧向抗滑稳定问题的研究就更少了，为了对这一问题做进一步研究和分析，本文将结合某实际工程，选取其岸坡坝段，采用AN-SYS有限元法，建立岸坡坝段三维有限元模型，在地震作用下进行应力计算，根据地基面上的作用效应和抗力计算结果，分析岸坡坝段在地震作用下的侧向抗滑稳定性。

1计算方法和模型

1．1抗滑稳定计算方法根据SL319－2005《混凝土重力坝设计规范》中规定，坝体抗滑稳定计算主要核算坝基面滑动条件，采用抗剪断强度公式计算坝基面的抗滑稳定安全系数［6］。式中:K'为按抗剪断公式计算的抗滑稳定安全系数;f'为坝体混凝土与坝基接触面的抗剪断摩擦系数;∑W为作用于坝体上全部荷载(包括扬压力)对滑动平面的法向分值，kN;c'为坝体混凝土与坝基接触面的抗剪断凝聚力，kPa;A为坝基接触面截面积，m2;∑P为作用于坝体上全部荷载对滑动平面的切向分值，kN。

1．2地震分析方法1．2．1拟静力法根据SL203－97《水工建筑物抗震设计规范》中规定，采用拟静力法计算重力坝地震作用效应时，只需考虑水平地震惯性力，计算公式为:Fi=αhξGEiαi/g(2)式中:Fi为作用在质点i的水平向地震惯性力代表值;αh为水平向设计地震加速度代表值;ξ为地震作用的效应折减系数，取ξ=0．25;GEi为集中在质点i的重力作用标准值;αi为质点i的动态分布系数;g为重力加速度［6－7］。动态分布系数应按照下式进行计算:式中:n为坝体计算质点总数;H为坝高，溢流坝的H应算至闸墩顶;hi、hj分别为质点i、j的高度;GE为产生地震惯性力的建筑物总重力作用的标准值。采用拟静力法计算重力坝地震作用效应时，水深h处的地震动水压力代表值应按下式进行计算:式中:pw(h)为作用在直立迎水坝面水深h处的地震动水压力代表值;ψ(h)为水深h处的地震动水压力分布系数;ρw为水体质量密度标准值;H0为水深。1．2．2反应谱法地震反应谱多指单自由度体系在给定的地震作用下最大加速度反应与体系自振周期的关系曲线。反应谱法的实质则为将计算结构的动力问题转化为静力问题进行求解，在计算地震响应时，将地面运动特征同结构动力特性一同纳入考虑的抗震计算方法［8－9］。单自由度体系动力方程为:式中:m为质量;x为位移;k为体系刚度;xg为地面运动水平位移;c=2mωζ，ζ为阻尼比，ω=(k/m)1/2。将各变量代入式(5)可得:综上所述，通过反应谱法对结构进行动力响应分析，即将各震型作用效应采用平方和开放的组合法求解。

1．3计算模型利用有限元ANSYS软件建立岸坡坝段三维有限元模型。水流方向为X轴方向，向下游为正，铅直方向为Y轴方向，竖直向上为正，垂直水流方向(坝轴线方向)为Z轴方向，向右岸为正，坐标原点在该坝段计算模型坝踵高程最低点处。计算基础深度取2倍坝高，上游自坝踵位置处向上游取2倍坝高，下游自坝趾位置处向下游取2倍坝高，左右岸方向上自坝肩或坝体横缝处分别向左右岸延伸2倍坝高。计算域内单元划分基本采用八节点六面体实体单元，部分通过四面体实体单元进行过渡，坝体和地基均采用单元SOLID45进行模拟，地基与坝体相接部位设有5m和3m的过渡区域。坝体和5m过渡区域采用八节点六面体实体单元，3m过渡区域采用四面体实体单元进行过渡。坝体和5m过渡区网格尺寸控制在1．0m以内，3m过渡区域以外地基单元边长为5～10m，如图1、2所示。

2计算分析

2．1拟静力法计算成果2．1．1岸坡坝段基底应力计算成果该坝段在计算工况(正常运行+地震)情况下用材料力学法的拟静力法计算时，其基底X向最大拉应力出现在坝体上游边缘及坝肩右侧边缘处，为0．6MPa;Y向最大拉应力出现在右侧坝肩处，为0．2MPa;Z向最大拉应力出现在右侧坝肩上游边缘处，为1．0MPa;切向最大剪应力出现在坝体上游边缘处，为0．6MPa。基底最大拉应力值、最大压应力值及其出现部位见表1。2．1．2岸坡坝段侧向抗滑稳定计算成果在计算岸坡坝段的侧向抗滑稳定中，将该坝 段的坝基面整体进行分块计算，该坝段的抗滑稳定计算结果见表2。

2．2反应谱法计算成果2．2．1岸坡坝段基底应力计算成果该坝段在计算工况(正常运行+地震)情况下用反应谱法计算时，其基底X向最大拉应力出现在坝体上游边缘及坝肩右侧边缘处，为0．4MPa;Y向最大拉应力出现在右侧坝肩处，为0．3MPa;Z向最大拉应力出现在右侧坝肩上游边缘处，为1．5MPa;切向最大剪应力出现在坝体上游边缘处，为0．6MPa。最大拉应力值、最大压应力值及其出现部位见表3。2．2．2岸坡坝段侧向抗滑稳定计算成果在计算岸坡坝段的侧向抗滑稳定中，将该坝段的坝基面整体进行分块计算，该坝段的抗滑稳定计算结果见表4。通过表1和表3对比分析可知，二者计算所得基底应力X向最大拉应力差幅0．2MPa基本接近，均出现在坝体上游边缘及坝肩右侧边缘;Y向最大拉应力基本一致，均出现在右侧坝肩处;但Z向最大拉应力，拟静力法计算结果为1．0MPa，反应谱法计算结果为1．5MPa，二者计算结果差幅0．5MPa偏大，均出现在右侧坝肩上游边缘处;切向最大剪应力差幅0．2MPa基本接近，均出现在坝体上游边缘处。2种方法计算所得应力虽有些偏差，但出现的部位一致，且均小于Ⅲ类地基岩体允许承载力4．0MPa，满足应力控制要求。对比分析表2和表4可知，根据抗剪断试验资料，其中抗剪断摩擦系数f'取1．0，抗剪断凝聚力c'取800kPa，在同样的分块情况及荷载组合下，侧向抗滑稳定计算结果二者均大于2．3，满足抗滑稳定要求，二者切向力计算结果基本一致，但反应谱法计算所得法向力结果较拟静力法偏大，从而反应谱法计算得到的抗滑稳定系数偏大。

3结语

2种方法计算结果对比分析得到以下结论:(1)二者计算得到的基底应力均小于Ⅲ类地基岩体允许承载力4．0MPa，满足应力控制要求。(2)二者计算得到的侧向抗滑稳定系数均大于2．3，满足抗滑稳定要求。(3)反应谱法计算所得应力及抗滑稳定系数较拟静力法偏大。

参考文献:

［1］李永明．高重力坝岸坡坝段侧向稳定分析研究［J］．云南水力发电，2015，31(01):46－49．

［2］丁旭，罗畅，王秦皓，王陆．蟒河口碾压混凝土重力坝岸坡坝段稳定及边坡开挖问题的研究［J］．水利水电工程设计，2009，28(04):46－47．

［3］宫照光，高鹏，周秋景，李海枫．某重力坝岸坡坝段上游面裂缝对工作性态影响分析［J］．中国水利水电科学研究院学报，2014，12(03):332－336．

［4］王超，张社荣，于茂．峡谷区重力坝岸坡坝段动力稳定性评估方法［J］．岩土力学，2014，35(05):1383－1390．

［5］邓良军，李双宝，吴余生．金安桥水电站大坝抗震分析和抗震措施设计［J］．水力发电，2011，37(01):26－28．

［6］邓东平，李亮，赵炼恒．地震作用下边坡稳定性分析的拟静力法研究［J］．中南大学学报(自然科学版)，2014，45(10):3578－3588．

［7］刘杰，李建林，张玉灯，屈建军，李剑武，李映霞．基于拟静力法的大岗山坝肩边坡地震工况稳定性分析［J］．岩石力学与工程学报，2009，28(08):1562－1570．

［8］杨璐，陈虹，岳永志，李明飞．反应谱法与时程分析法抗震分析对比［J］．沈阳工业大学学报，2016，38(03):331－336．

［9］苏亮，邵寿磊．多点地震输入的反应谱法研究［J］．工程力学，2014，31(12):83－88．

作者：邵明磊 杜慧超 单位：中国电建集团西北勘测设计研究院有限公司