龙源期刊网 http://www.qikan.com.cn
某重力坝复杂基础深层抗滑稳定分析
作者：王子健 陈志扬
来源：《科技创新导报》2017年第14期
DOI：10.16660/j.cnki.1674-098X.2017.14.045
摘 要：复杂基础的重力坝深层抗滑稳定计算是重力坝设计或复核的重难点。常规的刚体极限平衡法无法准确客观地寻找到最危险滑动面，该文以浙江省某水库重力坝为例，采用基于强度折减法的有限元分析，对正常蓄水位、设计洪水位和校核洪水位这三种工况分别进行了数值计算，计算稳定安全系数均能满足设计要求。通过本文实例，希望能为今后的类似工程设计提供一定的借鉴作用。
关键词：重力坝 抗滑稳定 有限元 强度折减法 复杂基础
中图分类号：TV22 文献标识码：A 文章编号：1674-098X（2017）05（b）-0045-05 Anti-sliding Stability Analysis of One Gravity Dam on Complex Foundation Wang Zijian Chen Zhiyang
（Zhejiang Design Institute of Water Conservancy & Hydro-electric Power， Hangzhou Zhejiang， 310002， China）
Abstract： Anti-sliding stability analysis of gravity dam on complex foundation is a key point in gravity dam design. Regular limited equilibrium method can’t help to find the most dangerous failure surface exactly and objectively. In this paper， anti-sliding stability analysis of a gravity dam in Zhejiang is simulated by FEM under three different conditions， based on strength reduction method. Results show that the calculating safety factors meet specification requirements. The study results provide a reference for similar projects.
Key Words： Gravity dam； Anti-sliding stability； FEM； Strength reduction method； Complex foundation
重力坝抗滑稳定计算一直是重力坝设计的重难点之一。规范中常规的抗滑计算是复核重力坝建基面的抗滑稳定安全系数，然而从实际情况来看，一般的大坝也不会出现在建基面上失稳的情况，更多地会出现深层滑动安全隐患。一般来说，大坝基础很少是完整的岩体，均或多或少地存在软弱结构层，深层抗滑稳定计算大多采用刚体极限平衡法，如何考虑多条软弱结构层的综合作用，寻找到最不利滑动面，具有较大的随机性和不确定性，同时也要求设计人员具有较丰富的工程经验。该文结合浙江省某水库重力坝的安全鉴定工作，对该重力坝进行了基于强
龙源期刊网 http://www.qikan.com.cn
度折减法的有限元数值计算，详细分析了该坝基础在三条不利软弱结构层共同作用时的深层抗滑稳定安全状况，为今后类似的工程提供借鉴作用。 1 工程概况
浙江省某水库始建于1964年12月。初建时水库库容850万m3，为小一型水库，主坝为浆砌块石重力坝，坝高29.3 m，坝顶高程66.96 m。1977—1990年进行水库扩建，大坝在原来的浆砌块石重力坝的基础上加高加宽，扩建为细骨料混凝土砌块石重力坝，扩建后坝高44.5 m，总库容2 838万m3，并在原溢洪道上修建副坝。2003年至2004年对主副坝进行除险加固，加固的主要内容：在主坝上游面高程47.16～73.16m间增设30 cm厚C25W6F50钢筋混凝土防渗面板，高程73.16～82.16 m采用HK-SN-918防水护面剂防渗，重打廊道内坝基排水孔。副坝上游面采用HK-SN-918防水护面剂防渗，增设上游侧防渗帷幕一道。新建放空预泄洞，进行电站改造及泄洪渠治理。2009年4月对主坝一、二期混凝土接触面进行化学灌浆处理。水库大坝现状见图1。
该重力坝建成时间较早，且后期又进行过加高加宽改造，而坝基又存在多条软弱结构带，因此，对该坝进行安全鉴定时，其抗滑稳定安全分析是整个工作的重点和难点。 2 抗滑稳定分析方法
根据《混凝土重力坝设计规范》（SL319-2005）规范要求，重力坝抗滑稳定计算主要核算坝基面滑动条件，一般采用抗剪断强度或者抗剪强度两种公式进行计算，计算公式如下： 抗剪断强度计算公式： （1）
抗剪强度计算公式： （2）
其中：K`为按抗剪断强度计算的抗滑稳定安全系数；K为按抗剪强度计算的抗滑稳定安全系数；f`为坝体混凝土与坝基接触面的抗剪断摩擦系数；f为坝体混凝土与坝基接触面的抗剪摩擦系数；c`为坝体混凝土与坝基接触面的抗剪断凝聚力；A为坝基接触面截面积；∑W为作用于坝体上全部荷载对滑动平面的法向分值∑P为作用于坝体上全部荷载对滑动平面的切向分值。
当坝基岩体内存在软弱结构面、缓倾角裂隙时，需核算深层抗滑稳定。根据滑动面的分布情况综合分析后，可分为单滑面、双滑面和多滑面计算模型，以刚体极限平衡法计算为主。然而，在实际存在复杂软弱结构面的坝基条件时，往往不同的滑动面选择计算所得到的结果相差很大，难以找到最不利的滑动面，具有较大的经验性和偶然性。在复杂基础抗滑稳定分析时，
龙源期刊网 http://www.qikan.com.cn
可采用有限元数值模拟的方式，对基础复杂软弱结构面进行数值模拟，分析出最不利的滑动条件，求得相应的滑动稳定安全系数。
该次计算采用有限元强度折减法，强度折减法最早是由Zienkiewics[1]提出，通过对材料原有强度的折减，找到结构物稳定的临界状态，临界状态所对应的折减系数即结构物的安全稳定系数。强度折减法应用于边坡稳定已有较多的工作[2-4]，对应于重力坝而言，是指将坝基岩体的真实抗剪强度除以一个折减系数F后再进行塑性计算，并以一定梯度逐渐增大F的大小，直到达到极限破坏状态为止，此时的折减系数即为重力坝的抗滑稳定系数。其计算公式为： c`=c/F （3）
Φ`=arctan（tanΦ/F） （4）
式中：c，Φ为岩土体真实的粘聚力和内摩擦角；c`，Φ`为折减后的粘聚力和内摩擦角。 此次计算采用摩尔库伦弹塑性本构模型，简称M-C模型。M-C模型基于理想弹塑性理论，由于加入了塑性理论，可以比较准确的描述岩土体的塑性破坏作用，相对于传统弹性模型而言有很大进步，考虑了大小主应力的影响，M-C模型屈服面为六棱锥形（如图2）。 3 建模计算与结果分析 3.1 重力坝模型及计算参数
选取主坝典型断面进行有限元分析，断面图（如图3）。 其中，主要的软弱结构带有：
#22软弱夹层，产状310°SW∠25°，宽度3～10 cm不等，呈扁片状花岗岩碎片，厚1～3 mm，其中部分为弱风化。该断层一直延伸至一期坝基下，从坝基断层缝内有渗水，流量估计10～15 mL/s，断层面凹凸不平，面上呈有黄绿色铁锰质染及绢云母化蚀变，性质差。 #26软弱夹层，产状310°SW∠16°，宽度2～5 cm，呈压扁片状，并有次生黄泥零星充填于空隙之间，厚1～2 mm，断层上下壁及破碎岩片均为黄褐色铁锰质染。
#27软弱夹层，产状280°～310°SW∠12°～25°，宽2～10 mm，局部为30～300 mm，性质自岩面至4 m深内为浅黄、浅灰、灰黑色粘性土，另为花岗岩压扁片状碎片，两侧为绿色蚀变矿物薄膜染。属于构造挤压错动，倾向下游，分布范围较广。
采用GeoStudio软件进行有限元建模计算，其中单元5 439个，节点5 528个，模型中考虑了#22、#26、#27三条软弱夹层，模型如图4。