平面剪切岩质边坡滑裂面的确定及稳定性分析

由于传统搜索方法对岩质边坡滑裂面的确定无法兼顾效率与精度, 如何迅速准确确定潜在滑裂面仍然是个难题。极限平衡法在岩质边坡稳定性分析中备受认可, 采用岩质边坡平面剪切滑动模型, 以滑裂面的倾角来表征潜在滑裂面的位置; 基于极值法, 推导了极限平衡条件下平面剪切破坏型岩质边坡潜在滑裂面的解析解, 并结合香港秀茂坪路边坡对其准确性进行了验证, 进一步对四川宜宾打营盘山公路多级边坡进行了整体稳定性分析。结果表明: 香港秀茂坪边坡采用本文方法确定的边坡潜在滑裂面倾角与实际滑坡倾角基本一致。实际工程应用中, 采用Slide软件中布谷鸟搜索法和模拟退火法两种搜索方法得到的滑裂面倾角分别为38.0°和37.0°, 本解析法所得倾角为34.8°; 选用Janbu法、Morgenstern-Price法和Sarma法分别计算对应的稳定系数, 结果均为1.04左右, 本文所得稳定系数为1.15, 可见本文方法所得结果基本准确。通过参数敏感性分析发现, 随着黏聚力的增加, 边坡滑裂面倾角越来越小, 稳定系数也随之增加; 而当内摩擦角增大时, 边坡滑裂面倾角和稳定系数也随之增大。      

均质土坡潜在滑面的搜索

均质土坡潜在滑面的确定 ,是对均质土坡进行稳定性计算的必要前提。传统的以瑞典圆弧法为代表的刚体极限平衡方法 ,虽然能够确定出均质土坡的滑面 ,但工作量较大 ,使用不方便。因此 ,在塑性极限分析方法的基础上 ,提出了一种新的搜索滑面的方法     

太原市杏花岭区小返乡后李家山村边坡稳定性分析

通过对太原市杏花岭区小返乡后李家山村的两处边坡进行勘察,主要查明岩土类型、成因、工程特性以及岩土的透水性和地下水的出露情况,查明坡体岩土的物理力学性质和软弱结构面的抗剪强度,对边坡的稳定性进行了分析,稳定性分析采用圆弧滑动法进行,两处边坡均处于不稳定状态,在降水条件诱发下极易发生滑动,必须尽快进行综合整治.     

济南高陡土质边坡稳定性分析及治理方法探讨——以济南南部山区山间平原为例

在济南山前及山间倾斜平原区常见土质边坡存在较大安全隐患,本文结合实例,对此类边坡开展稳定性评价及治理方法研究,对不稳定边坡的成因做了详细分析,并采用毕肖普条分法对边坡稳定性进行了定量计算,提出清坡+格构+钢筋混凝土挡板+截水沟的方案解决土质边坡失稳问题。治理后边坡滑动安全系数为1.673,实现了土质边坡的有效防治,为此类土质边坡失稳问题提供了新的治理思路。     

免棱镜全站仪在边坡岩体结构面测量中的应用

很多岩体工程中都涉及到岩体边坡的施工,经常会遇到岩体稳定性的问题,而岩体结构面的几何信息如产状在岩质边坡稳定性评价中发挥着重要作用,所以能否快速、准确地采集岩体结构面的几何信息成为了岩体稳定性评价成败中的关键一环。     

图表法在边坡抗滑桩设计中的应用

抗滑桩作为加固滑坡体的一种有效措施,与其它抗滑工程如抗滑挡墙、锚杆等相比,具有抗滑能力强、施工安全简便,并能进一步核实地质条件等突出优点,现已广泛应用于边坡工程中,但其设计方法还不成熟.运用图表法分析简化了某土质边坡抗滑桩设计计算,这种方法对以后的抗滑桩设计以及进一步深入研究具有一定的参考意义.     

用有限元强度折减法评价节理岩质边坡稳定性

以某边坡为例,论述了在ABAQUS软件中实现节理岩质边坡稳定性评价的建模过程。采用具有一定厚度的、低强度实体单元模拟节理组,按连续介质处理,节理组与岩块之间创建相互接触作用,通过不断折减结构面强度参数,以模拟达到失稳状态。通过分析坡顶典型质点场变量与X方向的位移关系,以数值计算不收敛作为边坡失稳评价标准的稳定性安全系数为1.41,以坡顶质点出现明显位移作为失稳标准的稳定性安全系数为1.32,两种判断标准均表明边坡目前处于稳定状态。通过积分等效塑性应变区演化图,模拟得到了强度折减状态下边坡失稳时的潜在滑动面,滑动面直观可视且物理意义明确,无需像极限平衡法一样要事先假定滑动面的位置和形状。     

基于3DEC的岩质边坡安全系数影响因素敏感性研究

为研究岩质边坡中结构面几何参数及力学参数对边坡安全稳定性的影响,文中采用离散元软件3DEC进行数值模拟试验;以安全系数为指标,结合响应面分析法,研究边坡结构面几何参数倾向、倾角及力学参数法向刚度(Kn)、切向刚度(Ks)、粘聚力(c)、内摩擦角(φ)对影响边坡稳定性敏感程度;结果显示,结构面的几何参数是边坡的稳定性的前提条件,在具备滑塌条件的情况下,结构面的力学参数,将进一步影响边坡的稳定性;其中结构面的切向刚度和内摩擦角对边坡稳定性有高敏感度,按照敏感度由大到小排序,为切向刚度内摩擦角粘聚力法向刚度;同时发现内摩擦角和切向刚度与边坡安全性因素之间呈相互促进作用。     

京广铁路K1219路基土质边坡深层滑移失稳机制与整治对策

为研究土质边坡深层滑移失稳机制, 以京广铁路下行线K1219+000处路基边坡失稳为例, 通过现场调查测绘、工程地质钻探、原位试验和室内试验、深部位移监测和数值模拟等手段, 详细研究了该土质边坡变形破坏特征、地质力学过程和失稳模式。结果表明: 路基边坡表面裂缝宽度及深度呈坡顶至坡脚逐渐变浅变窄, 变形具有一定的旋转性, 牵引式特征明显, 属于深层滑移拉裂式失稳; 土质边坡经历了因坡脚开挖、抽水引起的应力场和渗流场重新分布阶段、雨水入渗软化导致下滑力不断增大而滑面逐渐迁移扩大加深阶段以及支挡结构抗力失效阶段3个地质力学过程, 其失稳模式包括浅层滑移、浅层滑面向深层迁移、动荷载触发深层滑移失稳3个阶段。在此基础上, 综合确定了滑面位置, 并通过反演方法确定了滑面力学参数。研究采用了刚架式双排抗滑桩的整治方案, 通过理论计算和数值分析, 边坡变形与抗滑桩变形基本一致, 且与监测结果整体吻合, 这表明土质边坡深层失稳理论分析准确且计算的力学参数科学, 整治对策稳妥、可靠。      

基于内外功率比的非直线坡面型边坡的稳定系数研究

在边坡稳定性上极限分析中, 由于考虑了材料的理想弹塑性本构关系与相关流动法则, 相比极限平衡法更符合岩土材料的特征。在以往的二维边坡极限分析上限法中, 要求坡面形态为规则的直线, 而无论是天然边坡还是人工边坡, 边坡的坡面形态往往并非规则的直线。此外, 以往的边坡极限上限分析中得到的稳定数N_s=c/γH主要针对土坡的临界高度计算, 且未考虑孔压等外力对边坡施加的外功率。不同于传统极限平衡法采用静力学的平衡条件, 本研究针对非直线型坡面边坡稳定性问题, 假定滑动面为对数螺旋线, 基于极限分析上限定理和虚功原理, 提出了一种边坡发生旋转破坏时旋转中心的确定方法, 推导出了非直线型坡面边坡的重力虚功功率和能量平衡方程的解析解, 并提出了基于内、外功率之比的稳定系数K用以评价边坡的稳定性。通过算例比较了不同形态天然边坡的稳定性和人工削坡对边坡稳定性的影响, 并分析了边坡的坡度(β)、土体的内摩擦角(φ)、黏聚力(c)以及孔压系数(r_u)对稳定系数K的影响规律。对于坡度较大的边坡, 通过削坡改变坡面形态提高了边坡的稳定系数K。稳定系数K随黏聚力的增加而非线性增大, 随孔压系数的增加而降低。当黏聚力相对于孔压系数对边坡稳定性影响更大时, 稳定系数K随内摩擦角的增加而增大; 反之, 稳定系数K则随内摩擦角的增加而减小。以上结果符合对边坡稳定性分析的普遍认知, 验证了模型的合理性。另外, 通过将该方法与传统的Bishop法的计算结果进行对比, 发现安全系数F_s=1与稳定系数K=0的临界状态物理意义相同, 稳定系数K随黏聚力非线性增加, 更符合边坡的渐进破坏过程。      

Seismic stability analysis of slopes with pre-existing slip surfaces

In analyzing seismic stability of a slope with upper bound limit analysis method, the slip surface is often assumed as a log-spiral or plane slip surface. However, due to the presence of a weak layer and unfavorable geological structural surface or a bedrock interface with overlying soft strata, the preexisting slip surface of the slope may be irregular and composed of a series of planes rather than strictly log-spiral or plane shape. A computational model is developed for analyzing the seismic stability of slopes with pre-existing slip surfaces. This model is based on the upper bound limit analysis method and can consider the effect of anchor bolts. The soil or rock is deemed to follow the Mohr-Coulomb yield criterion. The slope is divided into multiple block elements along the slip surface. According to the displacement compatibility and the associated flow rule, a kinematic velocity field of the slope can be obtained computationally. The proposed model allows not only calculation of the rate of external work owing to the combined effect of self-weight and seismic loading, but also that of the energy dissipation rate caused by the slip surface, interfaces of block elements and anchorage effect of the anchors. Considering a direct relationship between the rate of external work and the energy dissipation rate, the expressions of yield acceleration and permanent displacement of anchored slopes can be derived. Finally, the validity of this proposed model is illustrated by analysis on three typical slopes. The results showed that the proposed model is more easily formulated and does not need to solve complex equations or time consuming iterations compared with previous methods based on the conditions of force equilibrium.     

A method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes

This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any point in a slope is determined using the Coulomb’s strength principle and the extremum principle based on the ratio of the shear strength to the shear stress at that point. The ratio, which is considered as an analysis index, can be computed once the stress field of the soil slope is obtained. The critical slip direction at any point in the slope must be the tangential direction of a potential slip surface passing through the point. Therefore, starting from a point on the top of the slope surface or on the horizontal segment outside the slope toe, the increment with a small distance into the slope is used to choose another point and the corresponding slip direction at the point is computed. Connecting all the points used in the computation forms a potential slip surface exiting at the starting point. Then the factor of safety for any potential slip surface can be computed using limit equilibrium method like Spencer method. After factors of safety for all the potential slip surfaces are obtained, the minimum one is the factor of safety for the slope and the corresponding potential slip surface is the critical slip surface of the slope. The proposed method does not need to pre-assume the shape of potential slip surfaces. Thus it is suitable for any shape of slip surfaces. Moreover the method is very simple to be applied. Examples are presented in this paper to illustrate the feasibility of the proposed method programmed in ANSYS software by macro commands.     

含优势渗流层边坡降雨入渗下的可靠度分析

含优势渗流层边坡在降雨入渗的作用下其渗流场往往具有较高的不确定性,这给边坡的稳定性评价带来困难,通常采用概率的方法解决此类问题。针对含优势渗流层边坡降雨入渗下的可靠度问题,通过将应力分析中的点估计-有限元法引入到边坡渗流-稳定性分析,提出了考虑优势渗流层渗透特性不确定性的渗流概率分析和边坡可靠度分析方法;其次以广西某含碎石夹层土坡为例,分析了降雨入渗下碎石夹层的优势渗流效应及渗流概率,并基于此开展了该边坡降雨入渗下的可靠度分析。结果表明:①含优势渗流层边坡雨水沿优势渗流层渗入坡体内部的深度显著高于沿坡面渗入的深度;优势渗流层渗透特性的不确定性对渗流结果的影响较大,使得边坡稳定性分析具有较强的不确定性;②随着雨水入渗持时的增加,含优势渗流层边坡不同滑动面的失效概率总体呈现增加趋势,最危险滑动面的位置不断向边坡下部演化;依托工程滑动面位置的预测结果与工程实际吻合;③提出的概率分析方法适用于分析含优势渗流层边坡降雨入渗影响下的稳定性问题,而且具有计算量小的优势,可作这类边坡可靠度分析的一种新方法。      

Evaluation of mountain slope stability considering the impact of geological interfaces using discrete fractures model

The stability of rock slope is often controlled by the existing discontinuous surfaces, such as discrete fractures, which are ubiquitously distributing in a geological medium. In contrast with the traditional approaches used in soil slope with a continuous assumption, the simulation methods of jointed rock slope are different from that of in soil slope. This paper presents a study on jointed rock slope stability using the proposed discontinuous approach, which considers the effects of discrete fractures. Comparing with traditional methods to model fractures in an implicit way, the presented approach provides a method to simulate fractures in an explicit way, where grids between rock matrix and fractures are independent. To complete geometric components generation and mesh partition for the model, the corresponding algorithms were devised. To evaluate the stability state of rock slope quantitatively, the strength reduction method was integrated into our analysis framework. A benchmark example was used to verify the validation of the approach. A jointed rock slope, which contains natural fractures, was selected as a case study and was simulated regarding the workflow of our framework. It was set up in the light of the geological condition of the site. Slope stability was evaluated under different loading conditions with various fracture patterns. Numerical results show that fractures have significant contributions to slope stability, and different fracture patterns would lead to different shapes of the slip surface. The devised method has the ability to calculate a non-circular slip surface, which is different from a circular slip surface obtained by classical methods.     

A generalized limit equilibrium method for the solution of active earth pressure on a retaining wall

In this paper, a generalized limit equilibrium method of solving the active earth pressure problem behind a retaining wall is proposed. Differing from other limit equilibrium methods, an arbitrary slip surface shape without any assumptions of pre-defined shapes is needed in the current framework, which is verified to find the most probable failure slip surface. Based on the current computational framework, numerical comparisons with experiment, discrete element method and other methods are carried out. In addition, the influences of the inclination of the wall, the soil cohesion, the angle of the internal friction of the soil, the slope inclination of the backfill soil on the critical pressure coefficient of the soil, the point of application of the resultant earth pressure and the shape of the slip surface are also carefully investigated. The results demonstrate that limit equilibrium solution from predefined slip plane assumption, including Coulomb solution, is a special case of current computational framework. It is well illustrated that the current method is feasible to evaluate the characteristics of earth pressure problem.     

基于最短折减路径法的边坡安全系数研究

使用双参数折减方法分析边坡稳定性的研究较多, 如何把两个折减系数定义为单一的综合安全系数是目前研究的一项重要内容。Isakov提出的最短折减路径法能够保证在不同工况下得到最小安全系数, 但是该方法的缺点在于计算复杂, 不适合工程应用。通过有限元数值模拟, 利用最短折减路径方法计算不同强度黏土构成的不同坡度均质土坡的最小安全系数和对应的折减系数, 探索了最小安全系数与土的初始黏聚力、内摩擦角以及边坡坡度的关系, 分析了初始强度对折减系数的影响。结果表明, 相同坡度下不同强度的黏土边坡在失稳时, 最小安全系数对应的临界破坏强度相同。临界破坏强度与坡度近似成线性正相关关系。由此基于最短折减路径法提出了一种新的计算最小安全系数的方法, 该方法得到的安全系数与目前常用的极限平衡方法所得结果相近, 并且计算简单, 因此可以用于边坡稳定性分析。      

基于一阶线性应变软化理论的边坡稳定性研究

在持续降雨或开挖卸荷作用下, 土体的强度指标会发生劣化, 但现今采用的边坡稳定性计算大多直接将其视为一个常数。为接近真实的边坡失稳破坏模式, 基于瑞典条分法以及一阶线性应变软化机制, 提出了一种边坡渐进破坏分析新方法, 推导出应变软化型边坡极限平衡表达式, 并获得了各破坏进度下边坡的安全系数。通过对模拟算例进行分析, 证明了条分-软化法的可靠性, 且计算结果表明渐进破坏过程中安全系数不仅取决于边坡的破坏方式与强度参数, 还与岩土体的软化模量密切相关。同时, 通过与滑坡实际案例的对比验算, 证实其强度指标存在不同的衰减系数, 即黏聚力的衰减系数大于摩擦角。从理论到应用, 最终获得的条分-软化法, 不仅考虑了岩土体的强度劣化效应以及滑动面的渐进发展, 还能有效地服务于实际工程背景下边坡的稳定性分析, 可以为滑坡的预防与治理提供指导建议。      

Analysis on velocity distribution and displacement along the profile of a slope using both empirical and analytical methods

Assessing the slope deformation is significant for landslide prediction. Many researchers have studied the slope displacement based on field data from the inclinometer in combination with complicated numerical analysis. They found that there was a shear zone above the slip surface, and they usually focused on the distribution of velocity and displacement within the shear zone. In this paper, two simple methods are proposed to analyze the distribution of displacement and velocity along the whole profile of a slope from the slip surface to the slope surface during slow movement. In the empirical method, the slope soil above the shear zone is assumed as a rigid body. Dual or triple piecewise fitting functions are empirically proposed for the distribution of velocity along the profile of a slope. In the analytical method, the slope soil is not assumed as a rigid body but as a deformable material. Continuous functions of the velocity and displacement along the profile of a slope are directly obtained by solving the Newton's equation of motion associated with the Bingham model. Using the two proposed methods respectively, the displacement and velocity along the slope profiles of three slopes are determined. A reasonable agreement between the measured data and the calculated results of the two proposed methods has been reached. In comparison with the empirical method, the analytical method would be more beneficial for slope deformation analysis in slope engineering, because the parameters are material constants in the analytical solution independent of time t, and the nonlinear viscosity of the soil can be considered.     

Evaluation of the possible slip surface of a highly heterogeneous rock slope using dynamic reduction method

A new method, the dynamic reduction method (DRM) combined with the strain-softening method, was applied to evaluate the possible slip surface of a highly heterogeneous rock slope of the Dagangshan hydropower station in Southwest China. In DRM, only the strength of the failure elements is reduced and the softening reduction factor K is adopted to calculate the strength parameters. The simulation results calculated by DRM show that the further slip surface on the right slope of the Dagangshan hydropower station is limited in the middle part of the slope, while both SRM (strength reduction method) and LEM (limit equilibrium method) predict a failure surface which extends upper and longer. The observations and analysis from the three recorded sliding events indicate that the failure mode predicted by DRM is more likely the scenario. The results in this study illustrate that for highly heterogeneous slopes with geological discontinuities in different length scales, the proposed DRM can provide a reliable prediction of the location of the slip surface.