基于交替变量局部梯度法对土质边坡的三维稳定性分析

临界滑动面的获得是土质边坡稳定性分析中一个必要的过程,然而现有分析方法难以保证该临界滑面的准确性,而且部分分析方法在计算过程中会陷入局部最小值陷阱。基于单变量方法和最速梯度法,提出了交替变量局部梯度法,该方法能够突破最优化方法应用在边坡稳定性分析上的瓶颈。首先,通过斯宾塞极限平衡方法和网格搜索法（或者其他确定临界滑面的方法）,获得边坡初始椭球状临界滑面;其次,在该滑面上布置若干节点作为变量,目标函数为安全系数Fs关于空间节点坐标Zi的方程,为了使目标函数快速降低,沿负梯度方向循环优化每个节点,当满足一定精度要求时,即可获得三维土质边坡非圆弧状临界滑动面;最后,通过算例证明该方法可行,计算结果可靠。     

基于临界滑动场的岩质边坡临界滑动面搜索方法

在边坡稳定分析中临界滑动面搜索是关键问题。现行搜索算法包括数学规划法和智能搜索算法,均是以圆弧或圆弧加折线为假定滑动面形式在计算土质边坡基础上发展起来的,基本上不适用于岩质边坡计算。临界滑动场理论以余推力法和最优化理论为基础,采用新的结构面处理方法,可计算出与实际情况基本吻合的最大余推力方向场,能方便追踪出各可能出口点的滑动面,从而找出边坡整体临界滑动面。取一均质土质边坡算例对计算结果可信度进行验证,与强度折减法和Spencer极限平衡法所得计算结果相吻合。另取一露天矿边坡对该方法的实用性进行了验证。     

地下矿体采动下的露天采场边坡稳定性研究

对于进入中晚期的露天金属矿,为了延长矿山的寿命或露天转地采稳产过渡,通常选择露天与地下联合开采,那么矿体采动下的边坡稳定性问题则尤为突出。在边坡几何形态、上部荷载和材料性质等不变的情况下,极限平衡获得的边坡安全系数为定值,因此,不能用来反映地下矿体采动对边坡稳定性影响。为了克服常规极限平衡方法的上述缺陷,采用基于边坡应力场的矢量和法求解滑动面的安全系数,结合蚁群算法与遗传算法联合搜索边坡临界滑动面,以确定矿体开采前后边坡的安全系数和临界滑动面。最后,给出一个假定算例和两个工程算例说明上述安全系数求解和临界滑面搜索方法的有效性,并揭示了地下矿体开采对露天采场边坡的稳定性的影响。     

The large-scale landslide on the flank of caldera in South Sulawesi, Indonesia

An extraordinarily large-scale landslide with a volume of about 200 million m³, a width of about 1,600 m, and a height of about 750 m occurred on 2004 March 26, 13:45 local time, on a steep caldera wall on the northwest flank of Mt. Bawakaraeng (2,830 m) at the headwater of the Jeneberang River in South Sulawesi, Indonesia. The debris avalanche extended about 7 km from the headwater and buried the river valley, causing devastating damage. There are a great many calderas in the world, notably Japan. If a large-scale sector collapse were to occur in a heavily populated area, it would be a devastating disaster for the people living in the area. The aim of this paper is to outline such a landslide and explain its mechanism of occurrence. We evaluated the stability of the original slope before the landslide using the limit equilibrium method and the finite-element-based shear strength reduction method (SSRFEM) with the strength reduction factor. The limit equilibrium method showed that a rise in the groundwater level caused the landslide. Although the critical slip surface predicted by the SSRFEM was shallower than that of the actual slip surface, the end positions of the actual and predicted slip surfaces were almost the same. Moreover, the end position of the critical slip surface before the landslide—the headwater of the Jeneberang River—was a knick point at which the slope inclination became steeper. SSRFEM analysis may be useful for evaluating the slope stability of large-scale landslides, because the critical slip surface predicted by it was close to the actual surface, even though we assumed homogenous conditions without information on the degree of weathering or ground properties. As the knick point formed at the end of the critical slip surface and is equivalent to the end of the actual slip surface, we assume such topographic features to be a primary geomorphologic cause of the landslide.     

基于随机响应面法的边坡系统可靠度分析

提出了一套基于随机响应面法的边坡系统可靠度分析方法。该方法首先从大量潜在滑动面中筛选出代表性滑动面。针对每条代表性滑动面,采用Hermite多项式展开建立其安全系数与土体参数间的非线性显式函数关系（即随机响应面）。然后,采用直接蒙特卡洛模拟计算边坡系统失效概率。在蒙特卡罗模拟中,采用所有代表性滑动面的随机响应面计算每一组样本所对应的边坡最小安全系数。最后,以两个典型多层边坡系统可靠度问题为例验证了该方法的有效性。结果表明：文中提出的边坡系统可靠度分析方法能够有效地识别边坡代表性滑动面,具有较高的计算精度和效率,并且确定代表性滑动面时无需计算滑动面间的相关系数。同时该方法可以有效地计算低失效概率水平的边坡系统可靠度,为含相关非正态参数的边坡系统可靠度问题提供了一条有效的分析途径。此外,多层边坡可能同时存在多条潜在滑动面,基于单一滑动面（如临界确定性滑动面）或者部分代表性滑动面进行边坡系统可靠度分析均会低估边坡失效概率。     

水位变化过程中边坡临界滑动场

库水位升降过程中,会引起边坡体内孔隙水压力分布的变化,对边坡的稳定性产生影响。对边坡临界滑动场进行改进,提出了可以考虑水位变化过程的边坡临界滑动场数值模拟方法,方便、快速地计算出边坡局部、整体安全系数在水位变化过程中的变化历程。通过对一个典型边坡在水位上升和下降过程中的稳定性分析,并和其他方法进行比较分析,结果表明,临界滑动场方法能搜索任意形状的最危险滑面,计算的安全系数是合理的。     

Strength reduction analysis for slope reinforced with one row of piles

In this paper, slope reinforced with one row of piles is considered. The critical slip surface is found to be divided into two parts when the pile spacing is small, and these two parts gradually get connected with the increase of pile spacing until a clear critical slip surface is formed. The critical slip surface of the piled slope is found to be usually shallower than the corresponding slope without pile, and this result is totally different from the previous results based on the maximum point of shear force where a very deep critical slip surface is obtained. The optimal pile position for slope reinforced with one row of piles is found to be located between the middle of the slope and the middle of the critical slip surface. An optimal design procedure for the pile spacing has also been suggested in the present paper, and the procedure will help to reduce the time required for a piled slope three-dimensional analysis.     

Three dimensional analysis of slope stability

A 3-D mathematical approach to slope stability, which is based on limiting equilibrium and variational analysis, is presented. In the initial formulation there are three unknown functions: the slip surface, the normal stress and the shear stress direction over this surface. The minimum factor of safety is sought through variational extremization. The analysis indicates that the factor of safety is independent of the normal stress distribution over the critical slip surface. It also indicates that the direction of the elementary shear force over the slip surface depends on the slip surface function, but not on the normal stress function. The analysis yields a non-linear first order partial differential equation, relating the slip surface and its first partial derivatives. By limiting the analysis to symmetrical problems an ordinary differential equation, governing the slip surface path on the plane of symmetry, is derived. This equation enables the development of a numerical procedure to determine the minimal factor of safety of symmetrical 3-D slopes. Two possible failure modes are determined for homogeneous slopes. One mode consists of finite 3-D sliding body and the second represents cylindrical failure. Numerical analyses for some simple cases, of homogencous slopes are presented.     

Mohr-Coulomb准则下基于滑动面深度的边坡参数反分析方法

利用上限定理与强度折减法,提出了一种新的边坡反分析方法。该方法基于边坡的滑动面形状,不需要位移、应力等监测数据。分析边坡的临界滑动面位置与强度参数的关系,可以发现,当边坡的几何形状、土体重度以及孔隙水压力分布确定,则临界滑动面的形状只与土体黏聚力c、内摩擦角φ的关系式c/tanφ有关,并在此基础上引入边坡的无量纲系数c,λcφ。利用极限分析上限定理与强度折减法,得到了边坡临界滑动面方程,提出了其深度的计算公式,建立了边坡c,λcφ值与滑动面深度的关系。通过已有边坡的反分析实例,验证了方法的正确性和有效性。     

The effects of method of generating circular slip surfaces on determining the critical slip surface by particle swarm optimization

The main objectives of slope stability analysis are evaluating factor of safety for a given slip surface and determining the critical slip surface for a given slope. Factor of safety is usually calculated by limit equilibrium method. The main steps to determine the critical slip surface are generating trial slip surfaces as probable solutions and searching among them to determine the one with the lowest factor of safety. Although the process of searching the critical slip surface received much attention between researchers, the significance of method of generating slip surfaces is seldom addressed in the literature. The authors believe that this ignorance can affect the accuracy of the results of slope stability analysis even in the simplest problems with circular slip surfaces. Consequently, this paper focused on the method of generating circular trial slip surfaces as the simplest mechanism of sliding and considered its effect on determining the critical slip surface. A new method of generating circular slip surface was presented, which is more efficient and less restricted than the conventional method. A computer program was also developed to determine the critical slip surface of slopes by using particle swarm optimization. The performances of the proposed method and developed computer program were verified during comparative studies and sensitivity analysis. Based on the results, the effect of method of generating circular slip surfaces on determining the critical slip surface was confirmed successfully. In all considered problems, the proposed method of generating circular slip surfaces led to the lower values of factor of safety compare with the conventional method.     

基于应力影响系数法的高边坡临界滑动面研究

临界滑动面一方面是与岩体结构面和坡面的组合关系有关,另一方面与荷载大小及位置有关。三维条件下,均质岩高边坡的破坏机理比较复杂,临界滑动面的判定较为繁琐。本文基于应力影响系数法对荷载作用下对均质岩高边坡三维临界滑动面形态及位置进行初步判别,研究各因素变化对临界滑动面分布形态的影响,结果表明:荷载作用下均质岩高边坡空间临界滑动面主要发生在基础前侧区域,约为一不规则楔体,楔体尺寸与桥基尺寸、坡度及桥基距离有关,其中桥基距离影响最大,荷载强度对滑动面形态无影响。桥基荷载作用下均质岩高边坡临界滑动线范围随桥基尺寸、桥基距离增加而增大,随坡度增加,临界滑动线长度减小,深度增大。本文提出桥基荷载作用下均质岩高边坡临界滑动面范围确定公式,即:均质岩高边坡破坏长度Pw=08444L-05098+05759S+693716,破坏深度Ph=01148L+0171B+11725+00733S-5719。     

Critical slip surface and landslide volume of a soil slope under random earthquake ground motions

Finding the critical slip surface and estimating the landslide volume are of primary importance for slope seismic design. However, this may be difficult due to the uncertainty of ground motions. To address this problem, a new method for calculating uncertainties is recommended in this paper, especially for the critical slip surface and landslide volume under random earthquake ground motions. Firstly, a series of intensity–frequency nonstationary random earthquake ground motions were generated based on an improved orthogonal expansion method. A given number of potential slip surfaces were set in a soil slope. Subsequently, the factor of safety (FOS) of each slip surface for all ground motions was calculated and the minimum FOS curves were obtained. It was found that the critical slip surfaces and failure times are uncertain under different earthquakes. The Monte Carlo method was used to verify the accuracy of probability density evolution method (PDEM), and the results of the PDEM and the Monte Carlo method are consistent, meaning that the PEDM has higher computational efficiency. Moreover, the distributions of earthquake-triggered landslide volume and landslide depth were analyzed by considering equivalent extreme events. Both landslide volume and depth exhibit a normal distribution for a homogeneous soil slope. The framework of this study is meaningful for slope seismic design in engineering, for example, the location of critical slip surface can be used for slope reinforcement, and the distribution of sliding volume can be used for disaster assessment.     

一种基于位移场分析的临界滑动面确定方法研究

安全系数的计算和临界滑动面的确定是边坡稳定性分析的两个组成部分,强度折减法可以得到合理的安全系数,但却不能准确刻画临界滑动面的位置。基于对滑动现象的认识,提出了一种基于位移场分析的临界滑动面确定方法,认为当边坡处于临界状态时,其潜在滑动面附近的位移等值线最为密集,潜在滑动面上的点往往是在深度方向上沿垂直滑面的位移变化率达到最大值的位置。将本文方法与Spencer法、岩质边坡模型试验进行了对比,验证了该方法在曲线和折线滑动面搜索方面的适用性和可靠性。同时,探讨了单元形状、疏密程度、离散点间距等因素对滑动面的影响。     

土质陡坡降雨侵蚀的数学模型及求解方法

建立了土质陡坡降雨侵蚀的数学模型,模型包括陡坡细沟间坡面流控制方程、变沟宽陡坡细沟流控制方程、陡坡细沟间坡面侵蚀泥沙连续方程及陡坡变沟宽细沟侵蚀泥沙连续方程。通过紊流冲击分布概率确定土壤剥蚀方程,确定了泥沙源/汇项的表达形式及细沟高程与宽度的表达式。通过有限元和有限差分相结合的方法求解该数学模型,利用有限单元法对数学模型的水动力方程、运动方程及泥沙连续方程进行空间上的离散,利用有限差分法解决时间域的问题,得出了顺序求解的数值计算公式及模型数值求解的具体方法和步骤。     

降雨过程中边坡临界滑动场

降雨入渗过程中孔隙水压力的升高与基质吸力的降低引起边坡稳定性的下降,是导致边坡滑塌的主要诱导因素。利用饱和-非饱和渗流有限元计算得到的孔隙水压力场,基于Fredlund提出的非饱和土抗剪强度理论,对边坡临界滑动场进行改进,提出可以考虑降雨过程的边坡临界滑动场数值模拟方法,能够方便、快速地计算出边坡局部、整体安全系数和相对应的临界滑动面在降雨过程中的变化历程。将该法用于一个典型均质边坡和一个非均质边坡在降雨过程中的稳定性计算,分析降雨持续时间、降雨强度和非饱和强度参数取值等因素对边坡稳定性的影响,并将计算结果与其他方法进行比较,结果表明临界滑动场方法能搜索任意形状最危险滑面,计算的安全系数合理。     

Two Methods for Predicting Reliability Index and Critical Probabilistic Slip Surface of Soil Slopes

Two methods of reliability analysis of soil slopes are studied, and the representative flow charts of both methods are illustrated. Method 1 can predict the reliability index and the critical probabilistic slip surface directly and it is computational efficient, but it needs the development of new codes for integrating the reliability analysis code and the slope stability code. Method 2 makes the reliability analysis code call the slope stability analysis code directly, and each code can be considered as an intact part. The main result of Method 2 is the reliability index of soil slope. Combined with the proposed method for locating the critical slip surface, Method 2 can also predict the probabilistic slip surface. Although Method 2 needs much more callings of the subprogram of slope stability analysis code, it needs not the developing of new computer program. Thus, Method 2 is easy to use and can be applied to different reliability analysis methods and slope stability analysis methods.     

A Simple Monte Carlo Method for Locating the Three-dimensional Critical Slip Surface of a Slope

Based on the assumption of the plain-strain problem, various optimization or random search methods have been developed for locating the critical slip surfaces in slope-stability analysis, but none of such methods is applicable to the 3D case. In this paper, a simple Monte Carlo random simulation method is proposed to identify the 3D critical slip surface. Assuming the initial slip to be the lower part of a slip ellipsoid, the 3D critical slip surface is located by means of a minimized 3D safety factor. A column-based 3D slope stability analysis model is used to calculate this factor. In this study, some practical cases of known minimum safety factors and critical slip surfaces in 2D analysis are extended to 3D slope problems to locate the critical slip surfaces. Compared with the 2D result, the resulting 3D critical slip surface has no apparent difference in terms of only cross section, but the associated 3D safety factor is definitely higher.     

基于最小势能法的三维临界滑裂面搜索方法

为了得到三维边坡的临界滑裂面,提出了6个参数控制的三维滑裂面构造方法,基于边坡三维最小势能稳定性分析方法建立目标函数,采用遗传算法实现了临界滑裂面的搜索,并开发了相应的搜索程序。为检验文中搜索方法的合理性,将其与其他方法得到的临界滑裂面以及最小安全系数进行比较,并且将室内模型试验得到的临界滑面与理论搜索的结果进行对比。研究表明：搜索方法可实现稳步收敛;针对算例进行多变量同时变化的搜索验算,得到了与极限平衡法较为接近的结果,表明提出的搜索方法是合理的;理论计算结果与室内模型试验坡面加载得到的边坡临界滑面较为接近,再次验证三维边坡临界滑面搜索方法是可行的。     

裂隙黄土边坡三维稳定性极限分析

黄土边坡中竖直裂隙的发育往往会对边坡稳定产生影响。相对于平面应变机制,建立三维破坏机制下边坡稳定性分析方法更能接近实际边坡失稳情况。基于塑性极限分析上限法,考虑预先存在竖直裂隙的三维黄土边坡不同破坏机制（坡面破坏、坡脚破坏和坡底破坏）,建立能量平衡方程及其无量纲临界高度值γH/c表达式,采用随机搜索法得到了临界高度的上限解。分析了约束宽度、边坡坡度、内摩擦角以及裂隙深度对三维竖直裂隙黄土边坡临界高度值的影响。结果表明：对于坡脚破坏机制,临界高度值随着裂隙深度的增加而减小,减小至临界裂隙深度 (δ /H)_min后,裂隙深度的增加不再影响临界高度值;临界裂隙深度随着坡度β 的增大而增大,随着内摩擦角φ 的增大而减小。当约束宽度B/H<0.8时,大多数破坏机制为坡面破坏。当约束宽度B/H=0.8、内摩擦角φ =10° 及约束宽度B/H=0.6、内摩擦角φ =15° 时,边坡的破坏从坡面破坏机制逐渐过渡到坡脚破坏机制。存在竖直裂隙的黄土边坡比完整边坡具有更小的临界高度,约束宽度及内摩擦角会对三维黄土边坡破坏机制产生影响。     

Determination of the critical slip surface in slope stability computations

The analysis of the stability of slopes using limiting equilibrium considerations necessitates the determination of the critical slip surface which yields the minimal factor of safety. The numerous methods currently available for slope stability analysis provide a procedure for assigning a factor of safety to a given slip surface, but do not consider the problem of identifying the critical conditions. This paper presents an effective minimization procedure based on dynamic programming by which the minimal factor of safety, and the corresponding surface, are determined simultaneously. This procedure SSDP (Slope Stability by Dynamic Programming), couples the minimization scheme with Spencer's method of slope stability analysis. It may be applied to slopes of any geometry, layering, pore pressure and external load distributions. No arbitrary restrictions are placed on the shape of the slip surfaces, and the analysis satisfies all equilibrium equations. Application of the procedure to slope stability problems reported in the literature shows that for a given slip surface the procedure yields factors of safety which are almost identical to those reported, but in every case a more critical slip surface, with a lower factor of safety, may be found.