General stability of two-dimensional slopes based on Sarma's method

Sarma's method of taking a fictitious accerelation as a measure of safety is used to study a two-dimensional slope. The slope is divided into arbitrary slices. The relations among forces acting on the slices and Sarma's acceleration are assumed to be linear. Consequently, a general analytical expression of Sarma's acceleration is derived by means of Cramer's rule. Furthermore, for four commonly used slice methods with simple relations among the forces on each slice and Sarma's acceleration, a general closed-form solution of Sarma's acceleration is given. An example is calculated and the results agree well with those of Hoek. Copyright © 1999 John Wiley & Sons, Ltd.     

膨胀土边坡稳定性非连续变形分析新方法

许多膨胀土边坡在长期的地质和干湿循环的自然气候作用下裂隙逐渐发展,以致被裂隙分割成块体系统。将膨胀土边坡看成是一个具有非饱和土本构关系、块体间存在摩擦力和膨胀力和具有牵引式滑动及时空效应特点的块体系统,建立了一种新的边坡稳定性非连续变形分析方法,并应用于某高速公路段膨胀土边坡稳定性的分析中,获到了很好的效果。     

节理岩质边坡变形的DDA模拟

对原非连续变形分析方法（DDA方法）中边界约束方法进行了推广，将其应用于节理岩质边坡变形规律的数值模拟，并利用模型实验结果验证了推广方法的正确性。     

考虑降雨入渗的非连续性岩体边坡稳定分析

由于断层和节理切割,边坡岩体一般具有非连续特性,降雨入渗是影响其稳定的重要因素之一。应用饱和与非饱和渗流理论建立了岩体边坡的非稳定饱和与非饱和渗流分析模型,并采用块体系统的非连续变形分析建立了降雨入渗非饱和渗流场影响的岩体边坡稳定分析和评价模型。作为应用实例,对某高速公路岩体边坡在降雨入渗条件下的非稳定饱和与非饱和渗流场及其对边坡稳定性的影响进行详细分析,并对边坡的稳定性进行评价。     

参数反演在岩质高边坡变形与稳定分析中的应用

由于边坡岩体本身的非均质性和赋存环境的复杂性,以及边坡岩体的力学特性参数及赋存环境总是处在动态的演化过程之中,因此,决定参数取值是数值计算中的一个难点,它直接关系到计算成果的准确性与可靠性。在材料参数反演分析的基础上,对三峡永久船闸边坡的开挖效应进行了分析,并进一步分析了三峡船闸的应力变形状态,旨在对船闸运行期的变形稳定分析提供依据。     

A development of the discontinuous deformation analysis for rock fall analysis

Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd.     

岩体基本参数对岩质边坡变形的影响

文章利用有限差分程序FLAC,定量化讨论了边坡岩体(软岩)的变形模量、泊松比、密度、粘聚力、内摩擦角、边坡高度等参数对软岩边坡变形的影响。研究表明,边坡位移随变形模量E的增大而减小;泊松比μ的增加,则会引起x方向的位移的增加,y方向的位移减小;而岩体密度对岩石边坡位移的影响甚弱。边坡岩体的粘聚力c对边坡变形的影响程度与岩体的其他参数间存在很大关系。内摩擦角对边坡位移的影响大小,则取决于粘聚力c的取值。此外,随着边坡高度的增加边坡变形则会随之增加。其中x方向位移变化对边坡高度的变化更为敏感;同样,在边坡岩体其他参数不变的情况下,边坡倾角越大,边坡变形也越大。     

计算边坡安全系数的坡向离心法

边坡稳定分析方法发展最引人注目,它是经典土力学最早试图解决而至今仍未圆满解决的课题。在常用边坡分析方法的基础上,从边坡失稳机制出发,提出一种更方便的安全系数计算方法--坡向离心法。该法通过不断增大水平加速度,直至边坡失稳为止,依据水平加速度与重力加速度失稳影响机制求得安全系数。通过算例与传统极限平衡法和有限元强度折减法相比较,并对各物理参数进行敏感性分析。结果表明,坡向离心法在边坡工程稳定分析中的应用是切实可行的,其弹性模量和泊松比对该法所求安全系数影响不大。     

基于模拟退火算法的边坡稳定分析方法

岩石边坡开挖中存在大量的随机结构面,结构面与岩桥的组合破坏是边坡破坏的一种重要形式。如何确定抗剪强度最小的滑裂面是边坡工程中的重要研究方向。在岩体结构面调查统计基础上,运用随机网络模拟的方法得到边坡的裂隙网络图像,在此基础上应用模拟退火算法,搜索出边坡由结构面与岩桥组合形成的潜在最危险滑裂面的位置,将结构面组合滑动面与Sarma分析方法结合,形成一种改进的分析岩石边坡稳定性的方法--Sarma-SA分析法,通过一工程实例验证了该方法的合理性。     

节理岩体边坡稳定性分析新方法

节理岩体边坡失稳破坏同时受控于节理与岩体抗剪强度。在对具有2组平行节理的岩体坡边失稳破坏机制研究基础上,研究了具有2组平行节理岩体边坡的极限平衡分析方法,并推导了相应的边坡稳定性分析计算公式,编制了基于潜在滑动面自动搜索的边坡稳定性研究程序。通过与文献中算例的对比研究,证明了所提出的节理岩体边坡稳定性分析方法的正确性,为具有2组平行节理岩体边坡的稳定性分析提供一条新的有效途径。     

再论岩质高边坡稳定性分析方法

岩质高边坡稳定性分析是岩土工程中经典的研究课题,其假定不同,分析方法就不同。从现有稳定性分析方法中存在的问题入手,指出其存在的优劣性,在刚体极限平衡和有限元等传统边坡稳定性分析方法的基础上,提出了用有限单元法直接求解边坡稳定性。该方法的思路是用潜在滑动面上在开挖、支护完成后尚未发生滑动的所有滑面单元模拟边坡在各工况下的真实应力场和变形场,根据这个真实应力场直接计算边坡安全系数。结合现场即时监测手段和计算机仿真技术可逐步模拟边坡在分步开挖期间的稳定性。算例分析表明,该方法完善了刚体极限平衡法的不足,克服了强度折减法强度参数c、?值折减比例的难题,弥补了超载法力学方向不明确的缺陷。计算所得的边坡安全系数和稳定状态与工程实际相当接近,表明提出的方法是可靠有效的新方法,为复杂岩质边坡稳定性分析提供了新的思路和途径。     

水力和超载条件下锚固岩石边坡动态稳定性拟静力分析

基于极限平衡理论,综合考虑水力条件、坡顶超载、地震荷载效应和锚固效应对岩石边坡进行了全面的稳定性分析。计算给出了多影响因素条件下岩石边坡稳定性安全系数的表达式,并重点分析了几种相关参数组合对岩石边坡稳定性的影响。分析表明,坡顶张拉裂缝积水、地下水渗流作用、滑面出流缝被堵塞、地震影响效应不利于岩石边坡抗滑稳定性,而锚索锚固效应则对提高边坡抗滑稳定性有积极作用;坡顶张拉裂缝积水、滑面出流缝被堵塞、水平向地震影响效应都不利于岩石边坡抗倾覆稳定性,但锚索锚固效应、坡顶超载、与竖直方向地震效应则对提高边坡抗倾覆稳定性有益。最后针对工程实际,提出了相应的工程建议。     

基于稳定性耦合分析法的余王扁边坡稳定性分析

边坡稳定性分析的极限平衡法与有限元法(FEM)的耦合分析法,首先利用有限元法分析获得边坡岩土体的整体应力场,在此基础上利用极限平衡法进行边坡的稳定系数求解。该方法既反映了边坡的稳定和变形之间的关系,又克服了极限平衡法与有限元法的不足,使二者的优点相互补充,获得的稳定系数基于极限平衡理论体系,可以同传统的稳定系数评价体系接轨。以西安市雁塔区余王扁削坡后边坡为例,用稳定性耦合分析法对其进行了稳定性分析,并把分析结果与各种极限平衡法的计算结果进行了比较,证明了耦合分析方法的可靠性和可行性。     

Modified predictor-corrector solution approach for efficient discontinuous deformation analysis of jointed rock masses

The high computational costs associated with the implicit formulation of discontinuous deformation analysis (DDA) have been one of the major obstacles for its implementation to engineering problems involving jointed rock masses with large numbers of blocks. In this paper, the Newmark-based predictor-corrector solution (NPC) approach was modified to improve the performance of the original DDA solution module in modeling discontinuous problems. The equation of motion for a discrete block system is first established with emphasis on the consideration of contact constraints. A family of modified Newmark-based predictor-corrector integration (MNPC) scheme is then proposed and implemented into a unified analysis framework. Comparisons are made between the proposed approach and the widely used constant acceleration (CA) integration approach and central difference (CD) approach, regarding the stability and numerical damping features for a single-degree-of-freedom model, where the implications of the proposed approach on open-close iteration are also discussed. The validity of the proposed approach is verified by several benchmarking examples, and it is then applied to two typical problems with different numbers of blocks. The results show that the original CA approach in DDA is efficient for the simulation of quasi-static deformation of jointed rock masses, while the proposed MNPC approach leads to improved computational efficiency for dynamic analysis of large-scale jointed rock masses. The MNPC approach therefore provides an additional option for efficient DDA of jointed rock masses.     

Seismic stability analysis of a layered rock slope

In general, the determination of the factor of safety and the location of the critical slip surface are two major challenges in seismic slope stability analysis. In this paper, a new approach for determining the factor of safety and the corresponding critical slip surface of a layered rock slope subjected to seismic excitations is presented, through a case study based on the combination of the shear strength reduction technique and distinct element method. According to this proposed method, the seismic factor of safety and the critical slip surface of the slope are estimated and compared with those obtained by the pseudo-static approach, combined with the limit equilibrium method. It is found that the factor of safety obtained from the proposed method is slightly greater than that computed by the pseudo-static analysis, with a difference of 4.2%, and that the critical slip surface obtained from the two methods is identical, which confirms the reasonability and feasibility of the proposed method.     

Computation of the safety factor for slope stability using discontinuous deformation analysis and the vector sum method

This study introduces the vector sum method into discontinuum-based methods by considering the sliding vector and the stress state of the discrete block system. The sliding direction computation and force projection in the new approach are detailed, and the safety factor is solved by explicit equations. The vector sum method is implemented in the discontinuous deformation analysis (DDA) program and is used to compute the safety factors for two numerical examples. A comparison of the solutions obtained with the theoretical analysis and limit equilibrium analysis demonstrates that the new method is suitable for calculating the safety factor of a slope.     

边坡稳定分析的非线性有限元混合解法

在对有限元强度折减法和有限元迭代解法的优缺点分析的基础上,提出了边坡稳定分析的非线性有限元混合解法,该法将有限元强度折减法和有限元迭代解法联合运用于边坡稳定分析,充分利用了两种方法的优点。即由有限元强度折减法搜索边坡可能滑动面,将可能滑动面在网格中画出,由迭代解法逐步迭代求得其滑动面的安全系数。并以算例和工程实例说明了此种方法的正确性和合理性。     

Boundary setting method for the seismic dynamic response analysis of engineering rock mass structures using the discontinuous deformation analysis method

Large deformations and discontinuous problems can be calculated using the discontinuous deformation analysis (DDA) method by solving time steps, and this method is suitable for simulating the seismic dynamic response of engineering rock mass structures. However, the boundary setting must be carefully analyzed. In this paper, four boundary settings for the DDA method are investigated. First, the contributions to the DDA equations for nonreflecting boundaries (including the viscous boundary and the viscoelastic boundary) are deduced based on the Newmark method. Second, a free‐field boundary is introduced in the DDA method with boundary grid generation and coupling calculation algorithms to accurately simulate external source wave motion, such as earthquakes. Third, seismic input boundary treatments are intensively examined, and the force input method is introduced based on nonreflecting boundaries. Finally, the static‐dynamic unified boundary is implemented to ensure consistent boundary transformation. The boundary setting method in the DDA method is discussed, and the suggested treatments are used to analyze the seismic dynamic response of underground caverns. Copyright © 2015 John Wiley & Sons, Ltd.     

Multi-scale limit load analysis for discontinuous rock mass based on the homogenization method

A new method to evaluate the strength of rock mass structures is proposed and examined. The method is based on the collapse load analysis of elasto-perfectly plastic material along with the homogenization method, which enables the multi-scale analyses for heterogeneous media. The homogenization process replaces a rock mass with cracks by an equivalent continuum medium with macroscopic stiffness while the failure criterion for the rock mass is estimated in the localization process. It is shown that both the averaged stiffness and the macroscopic failure criterion of the discontinuous rock mass are numerically obtained via the finite element analyses. Thus, the failure strength of a rock mass structure is evaluated by the collapse load analysis in the form of Linear Programming with the macroscopic failure criterion. This is the first attempt to apply the homogenization method to the strength analysis of rock mass. Copyright © 1999 John Wiley & Sons, Ltd.     

顺层岩质边坡稳定性极限分析上限法

基于极限分析上限法理论,运用体积力增量法,考虑单层滑动面极限分析模型的缺陷,建立考虑含结构面的多岩层错动的任意块体模型,利用岩块体在外荷载作用下达到极限破坏时,外荷载做的功与岩体消耗的功相等为基础建立等式方程,由此推导得到此类边坡稳定系数的计算公式。在分析典型边坡算例的基础上,对比离散元数值模拟结果,验证了该方法的可行性和适用性。对云南省普宣高速公路某顺层岩质边坡进行计算,其稳定系数接近离散元数值模拟结果,数值稍有偏大,边坡稳定。其研究结果为此类边坡稳定性评价提供了一种新的计算方法。