通过大型结构有限元分析软件ansys对一边坡进行稳定性分析，并通过杀死单元的方法来模拟边坡开挖。采用强度折减法来进行边坡稳定性分析，并且得到了边坡在开挖前后的x方向的位移变化以及塑性应变的变化。在边坡工程当中，在边坡上部进行基坑开挖是能够经常遇到的，这对于边坡来说既是一种卸载的过程，同时也是破坏了其原有结构的稳定。那么对于边坡的开挖进行稳定性分析是十分有必要的。本文立足于探求在边坡上进行基坑开挖前后的稳定性分析，采用有限元分析软件ansys对参考文献中国内某矿的边坡进行模型建立。
1有限元模型的建立
本例所采用的边坡的单元是材料库中的PLANE82单元，边坡的实例选取国内某矿，边坡模型建立选用弹性和弹塑性两种材料，定义的两种边坡围岩截面的材料属性参数如表1
表1边坡模型围岩材料属性参数
材料类别弹性模量/pa泊松比v容重kN/m3粘聚力/pa内摩擦角φ（。）
围岩13.2e+110.242700--
围岩23.0e+110.2525009e+542
模型建立以后进行单元离散和网格划分对模型左右边界限制其水平位移，下边界限制其水平及竖向位移，施加自重作用。设置边坡稳定收敛条件，求解到不收敛为止。在对正常边坡求解完成以后，通过杀死单元的方法来模拟边坡上部的基坑开挖，使其在同等条件进行求解，一直到边坡不收敛为止。
2有限元强度折减原理
在ansys对于平面应变下的强度问题，可以用Drucker—Prager法则来代替摩尔—库伦准则，而对于D—P法则模型需要三种材料的力学参数：粘聚力C、内摩擦角φ、膨胀角ψ，其中膨胀角ψ在本例中均设为0。而材料的粘聚力C和内摩擦角φ均使用强度折减法进行计算。
计算公式如下式所示：
式中，C和φ为边坡土体的原始的粘聚力和内摩擦角，C′和φ′是折减后的粘聚力和内摩擦角，F是强度折减系数。将折减后的参数输入模型求解计算，若完成求解，则此时边坡是稳定的；就增大折减系数F，直到计算至不收敛为止，此时的折减系数就是稳定系数。
3计算结果分析
运用ansys对边坡土体在不同的折减系数下的稳定性进行分析，求解完成以后，进入后处理查看器中，可以查看边坡的变形图、X方向位移云图以及塑性应变云图。随着折减系数的增大，边坡的X方向位移也在增大，塑性应变发展迅速，当塑性区发展形成一个贯通区域时，求解计算不收敛，此时认为边坡发生了破坏。
图1基坑开挖前F=2.4的塑性应变云图图2基坑开挖后F=2.4的塑性应变云图
通过对于基坑开挖前后边坡求解的结果，发现两者都是在F=2.6的时候计算不收敛，F=2.4时候的塑性应变云图如图1、2所示。
根据基坑开挖前后的塑性应变图可以看出，两者的塑性应变云图一致的，此时的塑性区即将贯通至坡顶。在总体变化上，基坑开挖以后的塑形图围绕着基坑边缘更快的向上贯穿，但是在基坑底部的土层仍未发生塑性区的开展。
将基坑开挖以前的边坡在不同折减系数下求解得到的X方向位移的最大值以及最小值的结果列于表2：
表2开挖基坑前后的边坡X方向位移值
折减系数开挖前X方向位移开挖后X方向位移
Min（mm）Max（mm）Min（mm）Max（mm）
F=1.0-32.65446.855-26.32650.357
F=1.2-32.66946.847-26.33750.362
F=1.4-32.87746.977-26.35850.371
F=1.6-46.20853.659-36.2157.078
F=1.8-61.77152.122-49.05455.025
F=2.0-76.54346.368-61.22348.027
F=2.2-90.98635.385-73.47534.80
F=2.4-103.79823.144-89.80619.212
4结语
根据基坑开挖前后的边坡的X方向的位移云图以及塑性应变云图可以得到以下结论：
4.1无论是基坑开挖前还是基坑开挖后，在X方向上的位移值在不同的折减系数的情况下有很大的改变。其中对于-X方向的位移值，随着F的增大，其位移的绝对值都是一直在增大。而且基坑开挖以后的在-X方向的位移值还要小于未开挖基坑的值，这是由于基坑开挖对于边坡相当于是一个卸载的过程。
4.2对于X方向的位移值，在刚开始的时候，位移值都是增大，当达到一定的强度系数以后，X方向的位移值开始下降。由于基坑的存在，开挖基坑以后的边坡在最开始出现塑性区开展时（F=1.4），是出现在基坑边角处；对于基坑开挖前的边坡，虽然也是在F=1.4的时候出现塑性区开展，但是是在坡脚处进行开展。
4.3对于基坑开挖后，F=2.6时，计算结果不收敛，塑性区贯穿整个边坡。当F=2.4时候，计算结果仍是收敛的，这与未开挖基坑时的结果一致。由于此时的强度折减系数就是稳定系数，说明基坑开挖后的边坡仍是稳定的。