An Extension of 2D Janbu’s Generalized Procedure of Slices for 3D Slope Stability Analysis II—Numerical Method and Applications

This paper provides a numerical approach on achieving the limit equilibrium method for 3D slope stability analysis proposed in the theoretical part of the previous paper. Some programming techniques are presented to ensure the maneuverability of the method. Three examples are introduced to illustrate the use of this method. The results are given in detail such as the local factor of safety and local potential sliding direction for a slope. As the method is an extension of 2D Janbu's generalized procedure of slices (GPS), the results obtained by GPS for the longitudinal sections of a slope are also given for comparison with the 3D results. A practical landslide in Yunyang, the Three Gorges, of China, is also analyzed by the present method. Moreover, the proposed method has the advantages and disadvantages of GPS. The problem frequently encountered in calculation process is still about the convergency, especially in analyzing the stability of a cutting corner. Some advice on discretization is given to ensure convergence when the present method is used. However, the problem about convergency still needs to be further explored based on the rigorous theoretical background.

Analysis of critical excavation depth for a jointed rock slope using a face-to-face discrete element method

The critical excavation depth of a jointed rock slope is an important problem in rock engineering. This paper studies the critical excavation depth for two idealized jointed rock slopes by employing a face-to-face discrete element method (DEM). The DEM is based on the discontinuity analysis which can consider anisotropic and discontinuous deformations due to joints and their orientations. It uses four lump-points at each surface of rock blocks to describe their interactions. The relationship between the critical excavation depth D-s and the natural slope angle alpha, the joint inclination angle theta as well as the strength parameters of the joints c(r) ,phi(r) is analyzed, and the critical excavation depth obtained with this DEM and the limit equilibrium method (LEM) is compared. Furthermore, effects of joints on the failure modes are compared between DEM simulations and experimental observations. It is found that the DEM predicts a lower critical excavation depth than the LEM if the joint structures in the rock mass are not ignored.

边坡稳定性分析极限平衡法的简化条件

边坡稳定性分析极限平衡法在求解安全系数时,需要对条块间内力引入假定条件或某些平衡条件在求解时得不到满足.分析这些不同的假设条件对于计算结果的影响,建立一个满足极限平衡原理和合理性条件的数值优化模型.引入载荷系数,利用载荷系数与安全系数之间的单调关系,将模型的求解转化为求解一系列线性规划问题,算法中不存在数值收敛问题,且可以求出安全系数的全局最优解.通过对模型各个约束条件的变化来模拟引入的假定条件,以此探讨各种假定条件或静力平衡条件的缺失对计算结果影响及某些常用方法引入的假定条件可能带来的误差.

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Effects of rainfall infiltration on deep slope failure

With the finite element method and the limit equilibrium method, a numerical model has been estab-lished for examining the effects of rainfall infiltration on the stability of slopes. This model is able to reflect the variations in pore water pressure field in slopes, dead weight of the soil, and soil softening caused by rainfall infiltration. As a case study, an actual landslide located at the Nongji Jixiao in Chongqing was studied to analyze the effects of rainfall infiltration on the seepage field and slope sta-bility. The simulated results showed that a deep slope failure is prone to occur when rainfall infiltration leads to a remarkable variation in the seepage field, especially when the pore water pressure in slopes increases in a large range.

地理信息系统支持下的滑坡灾害空间分析预测-以云南小江流域为例

As a typical geological and environmental hazard, landslide has been causing more and more property and life losses. However, to predict its accurate occurring time is very difficult or even impossible due to landslide's complex nature. It has been realized that it is not a good solution to spend a lot of money to treat with and prevent landslide. The research trend is to study landslide's spatial distribution and predict its potential hazard zone under certain region and certain conditions. GIS(Geographical Information System) is a power tools for data management, spatial analysis based on reasonable spatial models and visualization. It is new and potential study field to do landslide hazard analysis and prediction based on GIS. This paper systematically studies the theory and methods for GIS based landslide hazard analysis. On the basis of project "Mountainous hazard study-landslide and debris flows" supported by Chinese Academy of Sciences and the former study foundation, this paper carries out model research, application, verification and model result analysis. The occurrence of landslide has its triggering factors. Landslide has its special landform and topographical feature which can be identify from field work and remote sensing image (aerial photo). Historical record of landslide is the key to predict the future behaviors of landslide. These are bases for landslide spatial data base construction. Based on the plenty of literatures reviews, the concept framework of model integration and unit combinations is formed. Two types of model, CF multiple regression model and landslide stability and hydrological distribution coupled model are bought forward. CF multiple regression model comes form statistics and possibility theory based on data. Data itself contains the uncertainty and random nature of landslide hazard, so it can be seen as a good method to study and understand landslide's complex feature and mechanics. CF multiple regression model integrates CF (landslide Certainty Factor) and multiple regression prediction model. CF can easily treat with the problems of data quantifying and combination of heteroecious data types. The combination of CF can assist to determine key landslide triggering factors which are then inputted into multiple regression model. CF regression model can provide better prediction results than traditional model. The process of landslide can be described and modeled by suitable physical and mechanical model. Landslide stability and hydrological distribution coupled model is such a physical deterministic model that can be easily used for landslide hazard analysis and prediction. It couples the general limit equilibrium method and hydrological distribution model based on DEM, and can be used as a effective approach to predict the occurrence of landslide under different precipitation conditions as well as landslide mechanics research. It can not only explain pre-existed landslides, but also predict the potential hazard region with environmental conditions changes. Finally, this paper carries out landslide hazard analysis and prediction in Yunnan Xiaojiang watershed, including landslide hazard sensitivity analysis and regression prediction model based on selected key factors, determining the relationship between landslide occurrence possibility and triggering factors. The result of landslide hazard analysis and prediction by coupled model is discussed in details. On the basis of model verification and validation, the modeling results are showing high accuracy and good applying potential in landslide research.

Search For Critical Slip Surface In Slope Stability Analysis By Spline-Based Ga Method

The stability of a soil slope is usually analyzed by limit equilibrium methods, in which the identification of the critical slip surface is of principal importance. In this study the spline curve in conjunction with a genetic algorithm is used to search the critical slip surface, and Spencer's method is employed to calculate the factor of safety. Three examples are presented to illustrate the reliability and efficiency of the method. Slip surfaces defined by a series of straight lines are compared with those defined by spline curves, and the results indicate that use of spline curves renders better results for a given number of slip surface nodal points comparing with the approximation using straight line segments.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

三峡库区堆积层斜坡稳定性分级与应用研究

Considering the lacking of standard for the classification of Accumulative slopes so far, research working was conducted based on the results of geological investigation, data analysis and experiment carried out in Wanzhou. By mean of statistical method and grey system, the author studied in detail inflationary factors to Accumulative slopes. In order to study the mechanism of Rock-Soil Aggregate (RSA), numerical testing method was used. Coordinates in the two and three dimensional space and its corresponding rock fragments in the sample were generated randomly by VB and Particle flow code. After being built the models of RSA with different rock content, uniaxial and triaxial numerical simulation tests were carried out respectively. In order to study the effect of rainfall in Accumulative slopes, in situ infiltration testing had been conducted on site in Wanzhou, Three Gorges Area. Relationship between the infiltration rate and amount of precipitation has been obtained. Eleven factors are considered in the classification of Accumulative slopes in this paper.(1)On the basis of four basic factors and four inducing factors, sum-and-difference method for the classification system has been built. (2)After weight of factors being determined by analytic hierarchy process and membership function of Accumulative slopes stability being built in virtue of fuzzy mathematics, AHP-FM model of Accumulative slopes stability has been completed. In the end of this paper, having been applied on stability of Accumulative slopes in Three Gorge area and compared with result by limit equilibrium, classification system has good effect.

Release of Sulfur and Chlorine during Cofiring RDF and Coal in an Internally Circulating Fluidized Bed

An internally circulating fluidized bed (ICFB) was applied to investigate the behavior of chlorine and sulfur during cofiring RDF and coal. The pollutant emissions in the flue gas were measured by Fourier transform infrared (FTIR) spectrometry (Gasmet DX-3000). In the tests, the concentrations of the species CO, CO2, HCl, and SO2 were measured online. Results indicated when cofiring RDF and char, due to the higher content of chlorine in RDF, the formation of HCl significantly increases. The concentration of SO2 is relatively low because alkaline metal in the fuel ash can absorb SO2. The concentration of CO emission during firing pure RDF is relatively higher and fluctuates sharply. With the CaO addition, the sulfur absorption by calcium quickly increases, and the desulfuration ratio is bigger than the dechlorination ratio. The chemical equilibrium method is applied to predict the behavior of chlorine. Results show that gaseous HCl emission increases with increasing RDF fraction, and gaseous KCl and NaCl formation might occur.

Release of Sulfur and Chlorine During Cofiring Rdf and Coal In An Internally Circulating Fluidized Bed

An internally circulating fluidized bed (ICFB) was applied to investigate the behavior of chlorine and sulfur during cofiring RDF and coal. The pollutant emissions in the flue gas were measured by Fourier transform infrared (FTIR) spectrometry (Gasmet DX-3000). In the tests, the concentrations of the species CO, CO2, HCl, and SO2 were measured online. Results indicated when cofiring RDF and char, due to the higher content of chlorine in RDF, the formation of HCl significantly increases. The concentration Of SO2 is relatively low because alkaline metal in the fuel ash can absorb SO2. The concentration of CO emission during firing pure RDF is relatively higher and fluctuates sharply. With the CaO addition, the sulfur absorption by calcium quickly increases, and the desulfuration ratio is bigger than the dechlorination ratio. The chemical equilibrium method is applied to predict the behavior of chlorine. Results show that gaseous HCl emission increases with increasing RDF fraction, and gaseous KCl and NaCl formation might occur.

Prediction of soil organic carbon partition coefficients by soil column liquid chromatography

To avoid the limitation of the widely used prediction methods of soil organic carbon partition coefficients (K-OC) from hydrophobic parameters, e.g., the n-octanol/water partition coefficients (K-OW) and the reversed phase high performance liquid chromatographic (RP-HPLC) retention factors, the soil column liquid chromatographic (SCLC) method was developed for K-OC prediction. The real soils were used as the packing materials of RP-HPLC columns, and the correlations between the retention factors of organic compounds on soil columns (k(soil)) and K-OC measured by batch equilibrium method were studied. Good correlations were achieved between k(soil) and K-OC for three types of soils with different properties. All the square of the correlation coefficients (R-2) of the linear regression between log k(soi) and log K-OC were higher than 0.89 with standard deviations of less than 0.21. In addition, the prediction of K-OC from K-OW and the RP-HPLC retention factors on cyanopropyl (CN) stationary phase (k(CN)) was comparatively evaluated for the three types of soils. The results show that the prediction of K-OC from k(CN) and K-OW is only applicable to some specific types of soils. The results obtained in the present study proved that the SCLC method is appropriate for the K-OC prediction for different types of soils, however the applicability of using hydrophobic parameters to predict K-OC largely depends on the properties of soil concerned. (C) 2004 Elsevier B.V. All rights reserved.

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.