基于强度折减技术的糯扎渡心墙堆石坝稳定性分析.docxVIP

基于强度折减技术的糯扎渡心墙堆石坝稳定性分析 边坡稳定分析是土力学理论的重要组成部分，也是土木工程技术的宝贵基础。其现有的计算方法主要包括限制平衡法、限制元法和限制分析法。近年来随着计算机技术的发展, 有限元法在岩土工程中得到了广泛应用, 是分析边坡应力、变形和稳定性态的比较理想的方法。 目前强度折减技术特别适合于用有限元 (FEM) 数值计算方法来实现, 随着计算机技术的发展和数值计算技术的提高, 强度折减分析方法正成为边坡稳定分析研究的新趋势。笔者利用强度折减技术对糯扎渡心墙堆石坝进行安全稳定计算, 讨论各种破坏标准的优缺点, 为强度折减理论的进一步推广提供一定的科学依据。 1 临界强度参数的确定及弹塑性分析 Duncan (1996) 指出边坡的安全系数可以定义为边坡达到临界破坏状态时, 对土体的剪切强度进行折减的程度, 即定义安全系数是土体的实际剪切强度与临界破坏时折减后的剪切强度的比值。基于这种强度储备概念的安全系数Fs可定义为:当土体材料的抗剪强度参数c和φ分别用其临界抗剪强度参数cc和φc代替后, 结构处于临界破坏状态, 其中cc=c/Fs,φc=arctan (tanφ/Fs) 。在求解Fs时, 常应用c′=c/R,φ′=arctan (tanφ/R) , 假设外荷载不变, 来折减土体的强度指标c、φ, 然后对边坡进行弹塑性有限元分析, 通过不断地增加折减系数, 反复进行应力应变分析, 直至边坡达到临界破坏, 此时的折减系数就是安全系数Fs。 弹塑性非线性有限元分析需要考虑材料参数剪胀性对计算结果的影响。在弹塑性分析中如果采用同屈服函数形式相同的塑性势函数, 那么剪胀角同内摩擦角的关系将影响流动法则。为了保持在强度折减计算过程中流动法则不变, 对材料参数剪胀角也进行折减, 且折减公式采用同内摩擦角折减公式相同的形式, 如ψ′=arctan (tanψ/R) 。 2 抗剪强度与坡顶水平位移关系 (1) 采用解的不收敛性 (Zienkiewicz和Tayor, 1989) 作为破坏标准。在指定的收敛准则下算法不能收敛, 表示应力分布不能满足摩尔—库仑破坏准则和总体平衡要求, 意味着出现破坏, Ugai (1989) 指定迭代次数为500次, 超过这个限值认为土体破坏。这个迭代次数限值有很大的人为因素。 (2) Matsui基于剪应变屈服区域同边坡最终破坏面是一致的现象, 认为强度折减过程中边坡的破坏和剪应变的发展直接相关, 定义边坡破坏为剪应变破坏区域从坡角发展到坡顶, 而事实上塑性区从坡脚到坡顶贯通并不一定意味着破坏, 还要看是否产生很大的且无限发展的塑性变形和位移。根据应变分析的结果, 以规定数值的剪应变贯通整个边坡为破坏标准。该判别标准受土体剪应变计算精度以及剪应变发展不规律的限制。 (3) Lane认为边坡破坏和算法的不收敛是同时发生的, 而且伴随着结点位移 (或位移增量) 的迅速增加。依据这种现象, 将折减系数R与Eδmax/ (γH2) (无量纲的位移) 以图形形式显示, 其中δmax是收敛点的最大节点位移,H是边坡的高度, 依靠位移的网格图以及向量图来确定安全系数和显示边坡破坏机制的特性。 为了确定描述边坡破坏状态的客观标准, 绘出抗剪强度的折减系数同边坡顶点水平位移的关系曲线。发现当抗剪强度折减到一定程度后, 土体位移会持续增大, 折减系数与坡顶水平位移关系曲线趋于平缓。根据这一现象, 可以将边坡的“破坏程度”定义为强度折减系数与坡顶水平位移关系曲线斜率的函数, 斜率愈小说明边坡愈接近破坏。直到该斜率达到某一固定值时, 即认为边坡达到“破坏”状态。 以上各种标准均存在一定的局限性, 其中利用水平位移的变化趋势来判断坝坡是否破坏同样具有片面性, 所以笔者提出新的破坏标准——综合位移矢量标准。根据水平位移和竖向位移得出综合位移的数值和方向, 综合位移的方向采用位移矢量与纵轴的夹角来表示, 分别绘制出数值与角度随折减系数的变化趋势, 定义Sm=ΔR/ΔU,Sn=ΔR/Δα,R为折减系数,U为综合位移的数值,α为综合位移的角度。当Sm≤10-2且Sn≤10-2时, 认为坝体达到破坏状态。 下面以糯扎渡工程的计算结果为例, 进一步解释这种破坏标准的表达方式。 3 坝体边坡稳定性分析 糯扎渡水电站位于云南省澜沧江中下游河段, 是澜沧江中下游河段梯级规划二库八级电站的第五级, 也是梯级中的两大水库之一。控制流域面积14.47万km2, 多年平均年径流量546亿 m3。工程采用心墙堆石坝, 最大坝高261.50 m, 总库容237.03亿m3。 针对糯扎渡直心墙堆石坝的典型剖面, 用强度折减法和有限元分析软件进行坝体边坡的稳定性分析。坝体受重力、地震力以及水的作用力, 计算过程中重力分62级逐步加载, 模拟整个土