《表3 均质边坡强度参数：基于GGSO和有限元的土坡稳定性计算研究》

  提示：宽带有限、当前游客访问压缩模式

本系列图表出处文件名：随高清版一同展现

《基于GGSO和有限元的土坡稳定性计算研究》

1. 获取 高清版本忘记账户？点击这里登录

1. 下载图表忘记账户？点击这里登录

为验证本文方法的适用性，首先采用均质土坡算例来验算。该均质土坡几何结构如图4所示。土坡的土体的物理力学参数取值如表3所列，其中γ为土体容重、c为土体黏聚力、为土体内摩擦角。使用本文方法时，首先在COMSOL中划分计算网格（本例采用三角形单元，有限元网格划分见图5），然后通过MATLAB软件编制GGSO算法，产生若干潜在滑面，并通过MATLAB调用COMSOL计算这些潜在滑面上的应力值，尔后按照式（7）计算每个潜在滑面的稳定性系数，并从中选择最小的稳定性系数值及其对应的滑面作为边坡的最终稳定性系数和最可能滑面[2-6]。本例中，GGSO的种群数目设置为20，迭代次数为100。按照本文方法，确定式（8）中的圆心o为（44.93m，26.35 m），半径r为20.31 m，进而获得的边坡稳定性系数为1.301，作者同样采用Geo-studio软件计算该土坡的稳定性系数，得到结果为1.322。这里有必要说明，在土坡稳定性分析中，往往获取其真实滑面是比较困难的，而根据优化理论，在搜索所得滑面无明显异常的情况下（即滑面形状通常不会呈凸形），只需搜索得到的稳定性系数低于另一方法多得的稳定性系数，则可认为该方法更优。图4同时还给出了两个稳定性系数所对应的滑动面，可见搜索所得滑面是合理的。根据上文分析，可见GGSO搜索到的稳定性系数要更接近于土坡的真实稳定性系数。