基于遗传算法的岩土边坡稳定性分析方法.pdf

基于遗传算法的岩土边坡稳定性分析方法 毛忠义 张 军 （辽宁蒲石河抽水蓄能有限公司，丹东市振安区太平湾 118216） [摘 要] 在岩土边坡稳定性分析中，最常用的搜索临界滑移面的方法是基于传统的非线性优化方法。非 线性优化方法主要缺陷在于当边坡土层是非均质复杂条件下，无法保证搜索到安全系数的全局最小解，而 只能搜索到局部极小解。基于遗传算法，建立了搜索岩土边坡稳定性分析最小安全系数和滑移面坐标与半 径的数值方法。问题的可行解在变量搜索区间内搜索，包括滑弧的圆心坐标和半径。数值计算结果表明， 与其它的基于梯度搜索方法相比，遗传算法能够搜索到更好的解。 [关键词] 边坡稳定 遗传算法 滑移面 安全系数 1 引言 边坡稳定性评估是岩土工程的重要课题，在土木工程及水电工程建设中,经常遇到岩土边坡工程 的稳定性问题。为了确定边坡潜在的最危险滑动面及其对应的边坡稳定性安全系数,已有许多方法被 提出,其中包括极限平衡法、极限分析法、变分法、有限元法等。多年来,各类工程技术人员在研究 边坡的滑动机理、发展规律及其稳定性分析计算方法等方面,均取得了较大的进展。边坡稳定问题涉 及到水利、港工、采矿、铁路及各类工业与民用建设工程，因此研究边坡的稳定性问题具有重要的 意义。边坡稳定性是大坝和露天矿以及高速公路、铁路、机场、高层建筑深基坑开挖等土木工程建 设中十分重要的问题。现在岩土边坡的开挖深度已高达数百米，边坡稳定性分析结果的正确与否直 接关系到边坡工程的成败，甚至还关系人们的生命财产安全。目前用来进行边坡稳定性分析的方法 有很多，但在实际应用中多数的边坡稳定性分析还是使用传统的极限平衡处理方法。所有平衡方法 的一个共同的难点是它们都要以失效的土条能被条分这种假设为基础。所有这些方法的大体过程可 以概括为：1）给定一个滑面；2）对滑面的剪压静力分析；3）计算安全系数k，其中k被定义为整个 滑面的抗剪强度与实际产生的剪应力之比；4）通过多次搜索确定k 及其所对应的临界滑面。刘珊 min 珊 （2006）采用STAB边坡稳定计算软件分析了蕉东闸桥工程施工中的边坡稳定问题[1]。张鲁渝 （2003） 采用有限元强度折减系数法计算土坡稳定安全系数的精度，通过对106 个算例的比较分析, 表明折 减系数法所得稳定安全系数比简化Bishop 法平均高出约5.7%, 且离散度极小, 这不仅验证了文中 所提措施的有效性, 也说明了将折减系数法用于分析土质边坡稳定问题是可行的[2]。杜运兴（2006） 提出了基于遗传算法的边坡稳定有限元分析方法，该方法采用Marc 有限元软件计算边坡的应力,考 虑了几何和材料非线性. 采用一条折线模拟边坡滑动面. 根据遗传算法原理编制程序,搜索对应最 小安全系数折线上的节点,从而确定最危险的滑动面[3]。殷建华(2004)采用刚体有限元上限分析方法 和非线性规划方法研究了求解最小安全系数的求解问题[4]，曹烘(2004)基于突变理论模型研究高边坡 的局部和整体稳定性问题[5] [6] ，周翠英(2004) 基于GIS技术研究区域性滑坡发生概率方法 ，李守巨 [7] （2003）采用模拟退火算法研究搜索岩土边坡最小安全系数的搜索方法 。 计算表明，利用传统的非线性最优化技术解决边坡稳定性分析，会得多个数值解。这些技术的 主要缺点就是其非鲁棒性，即在复杂多样的土质条件下，安全系数的局部最小解不是全局最小解。 本文的目的就是要说明对于复杂的土工边坡，存在多个局部最小解，但是安全系数的全局最小解及 其所对应的临界滑面只有一个。传统的最优化方法无法收敛于安全系数的全局最小解[10] 。而遗传算 法是一种基于遗传学原理的优化方法，它可以通过一连串的进化找到全局最优解。遗传算法拥有强 大的全局收敛性和并行搜索能力，可以避免局部最优解问题[11] 。 2 稳定性分析计算的条分法 条分法先假定若干可能的剪切面——滑裂面。然后将滑裂面以上土体分成若干垂直土条，对作 用于各土