边坡角度优化设计的二分法.docxVIP

边坡角度优化设计的二分法 1 基于数值计算方法的岩质边坡稳定分析 露天矿边坡设计对露天矿的生产安全和经济效益有重大影响。过小的边坡角,将增加剥岩量,影响矿山的经济效益,而过大的边坡角,将导致岩石塌落和滑坡事故的发生,严重的影响矿山的正常生产,因此边坡角的优化设计应达到最优状态,既能保证边坡的稳定性,又能提高矿山的经济效益。 目前边坡角的设计方法主要是极限平衡方法,可用于分析具有少量节理切割的少量块体,对于强节理化岩体显得无能为力。因此许多学者将更多的数值分析方法(如刚性元法、等效连续模型、离散单元法(DEM)、块体理论、DDA)引入用于分析边坡的稳定性。 在此基础上,许多学者将更丰富的理论基础引入到边坡角度的优化设计中,和数值计算方法共同使用,取得了一定的成效。例如采用离散元数值计算方法与神经网络预测方法进行了大量的工作,并取得了一定的理论成果;对复杂岩质边坡,利用有限元强度折减法由程序自动求得岩质边坡的滑动面以及相应的安全系数,为岩质边坡稳定分析开辟了新的途径;基于弹塑性有限元分析的计算结果,提出二维情况下边坡潜在滑移线满足一个常微分方程组初值问题,为自动搜索边坡的潜在滑移线提供了有效的方法。 基于前人的研究发现,本文拟采用数学中的二分法对边坡角度进行优化设计,并采用有限差分数值计算方法对边坡稳定性进行校验,模拟计算过程中采用的力学参数不变,进而保证了计算结果的可比性。 2 bb满足[a,bx]的夹逼区间[a,b]的[a,b]的夹逼区间[a,b]的夹逼区间] 二分法属于数学问题,如图1,其原理为:设单变量函数B (x)在区间[a,b]上为单调连续函数,若B (a)·B (b)=0,则B x)=0时的根ε必存在于区间[a,b]中,此时[a,b]就称为B (x)的夹逼区间。取夹逼区间[a,b]中点时,有以下三种情况: 边坡角度的优化设计根据二分法原理,具体形势如下:设边坡的稳定性函数f (x)在边坡角度区间[0°,90°]上为单调连续函数,x1为临界边坡角,则有: 边坡角度具体优化流程见图2。 3 实例分析 3.1 基于flac2c有限差分程序的边坡稳定性分析 首钢水厂铁矿边坡为高陡深凹露天边坡,本文取其中具有代表性的最高剖面进行边坡角优化设计及稳定性分析,所取剖面边坡包含F36、F11和F8三条断层,图3给出了剖面边坡的岩性和断层的分布图。 根据上述边坡角度的优化流程进行优化设计,边坡角度误差为0.5°。根据生产技术要求,边坡需在+50m水平设置运输平台,平台宽40m,由于边坡较高,运输平台上部边坡服务年限较长,下部边坡服务年限相对较短,因此边坡以运输平台为界,上下边坡分别采用不同的边坡优化角度,即上部边坡较缓,下部边坡较陡。 本文采用FLAC2D有限差分程序验算优化后的边坡稳定性,边坡计算模型高度从水平—600m起,一直模拟到地表,高度925m (图3)。边坡稳定性主要根据边坡的破坏场、位移量以及位移速度的变化趋势综合判断,边坡岩体采用莫尔一库仑屈服准则,计算采用的岩体力学参数见表1。 3.2 加固后边坡局部变化的监测 下面给出了边坡角度(见图4～图9)优化过程中有代表性的几组数据: 方案1 (45°～45°):上部边坡角45°,下部边坡角45°;上、下部分界位于+50m水平左右,以下同; 方案2 (47°～47°):上部边坡角47°,下部边坡角47°; 方案3 (43°～43°):上部边坡角43°,下部边坡角43°; 方案4 (42°～45°):上部边坡角42°,下部边坡角45°; 方案5(43°～45°):上部边坡角43°,下部边坡角45°; 方案6(44°～45°):上部边坡角44°,下部边坡角45° 从有限差分数值模拟计算出的破坏场可以得出,边坡角为(45°～45°)时,边坡的破坏范围较小;边坡角度提高到(47°～47°)时,坡脚出现剪切滑移带,易造成边坡局部垮落,影响边坡的局部稳定性;边坡角为(43°～43°)时,坡脚暴露的岩体力学性质较差,坡脚出现剪切滑移带。当边坡角为(42°～45°)、(43°～45°)和(44°～45°)时,边坡的破坏范围与坡角为(45°～45°)时基本相同,坡脚不出现剪切滑移带。因此,需要通过边坡的位移量以及位移速度变化趋势对边坡的角度进行进一步确定。 表2给出了坡面上典型部位的水平位移和水平移动速度的监测结果,数据表明边坡角为(43°～43°)和(47°～47°)时,边坡位移量大幅度增加是边坡失隐的隐患,与破坏场的判断结果一致。 坡角为(45°～45°)、(42°～45°)、(43°～45°)和(44°～45°)时,边坡产生的最大水平位移量基本相同,典型位置处的位移量基本保持不变,但边坡角为(45°～45°)、(42°～45°)和(44°～45°)时,典型位置处有三点以上的移动速度不趋