金佛山沥青混凝土心墙坝坝坡稳定性有限元剖析.docVIP

金佛山沥青混凝土心墙坝坝坡稳定性有限元剖析 　　摘要：针对重庆市金佛山水利工程沥青混凝土心墙堆石坝，采用有限元强度折减法，研究了上、下游坝坡在竣工期、蓄水期及水位骤降等不同工况下的稳定性。结果显示，在竣工期、蓄水期和水位骤降期的坝坡安全系数均符合规范要求。与极限平衡分析方法相比，强度折减法在滑弧面及塑性区的分析中弥补了极限平衡法中无法准确确定滑动面具体位置及滑体大小的缺陷，但是强度折减法计算结果偏安全，因此将两种方法结合进行稳定性分析，结果更为客观。 　　关键词：沥青混凝土心墙坝；坝坡稳定性；有限元强度折减法；判定准则 　　中图分类号：TV64文献标识码：A文章编号：1672-1683（2013）04-0206-04 　　1研究背景 　　目前，在岩土工程中广泛采用的边坡稳定性分析方法主要是传统的极限平衡法和有限元法，极限平衡法包括普通的条分法、毕肖普法、杨布条分法等，都需要做一定假设[1]。自有限元法引入边坡的稳定性分析以来，在边坡工程治理中得到了很大的发展[2-4]，实现了边坡在荷载作用下的应力、变形、滑动面等的计算[5]。其中，有限元强度折减法被部分学者用于土石坝的稳定分析：邹兵华[6]采用有限元强度折减法分析了均质中低坝型窑店煤矿拦渣坝的自重和堆渣两种工况的坝坡稳定性；彭凯[7]等提出采用坝体典型断面特征点的位移矢量与强度折减系数之比作为作为判别坝坡破坏的标准，当比值出现突变时即认为处于极限平衡状态。不过，尽管有限元强度折减法优点甚多，然而在土石坝中应用不广泛，原因之一是人们对有限元强度折减法准确性有疑虑。 　　本文将滑坡的变形破坏规律应用到土石坝的坝坡稳定性分析中，以金佛山沥青心墙堆石坝为例，采用毕肖普法、瑞典条分法分析有限元强度折减法的计算精度，探讨有限元强度折减法在填料复杂的高土石坝的稳定分析中的运用的可行性，并通过塑性区等值线云图对大坝坝坡安全性进行分析。 　　2工程概况 　　金佛山水利工程位于重庆市南川区、万盛区以及贵州省桐梓县，以灌溉、供水为主，兼顾发电等综合利用的大（Ⅱ）型水利工程。水库正常蓄水位836.00 m，相应库容9 614.6万m3，水库总库容1.03亿m3。水库枢纽工程由拦河大坝、溢洪道、泄洪放洪洞、坝后式电站、渠首工程等组成。按坝体结构要求，除沥青混凝土心墙外，大坝其他部位共分6个区，分别为上游堆石I区，下游堆石Ⅱ区和堆石Ⅲ区，下游砂岩大块石排水棱体Ⅳ区，及上下游过渡层Ⅰ和Ⅱ区，见图1。图1沥青混凝土心墙堆石坝坝体剖面图 　　Fig.1Profile of asphalt core concrete wall rock-fill dam 　　3计算方法 　　3.1.1有限元强度折减法 　　强度折减法原理是在理想弹塑性有限元计算中，通过不断折减岩土体的抗剪强度参数 c、φ，直至边坡达到临界破坏状态[2]，此时的强度折减系数就是边坡的稳定安全系数Fs，定义为土体实际的抗剪强度与折减后极限状况下抗剪强度的比值。土体的抗剪强度参数的折减公式为： 　　4计算结果分析 　　4.1失稳破坏判定准则 　　目前，判断边坡失稳破坏的标准通常包括：有限元数值计算的不收敛、塑性区的贯通、广义剪应变的贯通等[3]。吕庆等[11]认为在小应变假设中用数值计算不收敛作为判据，但是，计算不收敛的原因比较多，或者荷载过大，或者计算单元有奇异等，因此，以此为判据适用范围有一定的限制；栾茂田[12]建议采用塑性应变贯通作为判据，但此判据中主观因素占很大成分，并未排除弹性塑性应变的影响，破坏界限比较模糊。 　　分析坝坡失稳破坏的主要特征，不管其内部的变形机理是广义剪应变还是塑性应变，最终结果都是产生位移，位移是坝坡内部作用的外在表现。滑动主要是由剪应变和位移造成的。随着强度参数的不断折减，坝坡上的位移矢量和剪应变不断向坡脚处增大，因此，以坝顶特征点位移突变为失稳判据，意义明确、界限清晰。 　　4.2有限元计算结果分析 　　经过有限元强度折减法的坝坡稳定的分析，上下游坝坡均产生一定的位移和变形，与上游的蓄水期、水位骤降期及下游竣工期相比，竣工期有明显的滑动趋势。竣工期上游坝坡计算结果见图4。 　　图4显示了坝坡的位移塑性区总体从坝底向坝顶贯通，宏观上表现为坝坡滑弧，即在位移等值线云图中显示出来的滑弧。结合有限元强度折减系数与位移矢量关系图可以看出，折减系数为1.72时是位移矢量的拐点，之后位移矢量开始陡增；若以计算不收敛作为判定准则，安全系数为1.74。但是，从塑性区的贯通现象能够看出，塑性区贯通后坝坡自然位移快速增加，此时整个滑体开始滑动，即当Fs=1.72时，坝坡已经处于滑动过程中了。因此，将位移矢量与强度折减系数关系图中的“位移矢量陡增突变”作为判定准则，可以清楚地解释坝坡破坏启动的瞬间安