Seismic-induced surficial failure of cohesive slopes using three-dimensional limit analysis:A case study of the Wangjiayan landslide in Beichuan, China

A seismic-induced landslide is a common geological catastrophe that occurs in nature.The Wangjiayan landslide,which was triggered by the Wenchuan earthquake,is a typical case in point.The Wanjiayan landslide caused many casualties and resulted in enormous property loss.This study constructs a simple surficial failure model based on the upper bound approach of three-dimensional(3D)limit analysis to evaluate the slope stability of the Wangjiayan case,while a traditional two-dimensional(2D)analysis is also conducted as a reference for comparison with the results of the 3D analysis.A quasi-static calculation is used to study the effect of the earthquake in terms of horizontal ground acceleration,while a parametric study is conducted to evaluate the critical cohesion of slopes.Rather than employing a 3D analysis,using the 2D analysis yields an underestimation regarding the safety factor.In the Wangjiayan landslide,the difference in the factors of safety between the 3D and 2D analyses can reach 20%.The sliding surface morphology as determined by the 3D method is similar to actual morphology,and the parameters of both are also compared to analyze the reliability of the proposed 3D method.

An approach for determination of lateral limit angle in kinematic planar sliding analysis for rock slopes

Planar sliding is one of the frequently observed types of failure in rock slopes.Kinematic analysis is a classic and widely used method to examine the potential failure modes in rock masses.The accuracy of planar sliding kinematic analysis is significantly influenced by the value assigned to the lateral limit angleγlim.However,the assignment ofγlim is currently used generally based on an empirical criterion.This study aims to propose an approach for determining the value ofγlim in deterministic and probabilistic kinematic planar sliding analysis.A new perspective is presented to reveal thatγlim essentially influences the probability of forming a potential planar sliding block.The procedure to calculate this probability is introduced using the block theory method.It is found that the probability is correlated with the number of discontinuity sets presented in rock masses.Thus,different values ofγlim for rock masses with different sets of discontinuities are recommended in both probabilistic and deterministic planar sliding kinematic analyses;whereas a fixed value ofγlim is commonly assigned to different types of rock masses in traditional method.Finally,an engineering case was used to compare the proposed and traditional kinematic analysis methods.The error rates of the traditional method vary from 45%to 119%,while that of the proposed method ranges between 1%and 17%.Therefore,it is likely that the proposed method is superior to the traditional one.

Analysis of Low-Frequency Vibrational Modes and Particle Rearrangements in Marginally Jammed Amorphous Solid under Quasi-Static Shear

We present the numerical simulation results of a model granular assembly formed by spherical particles with tIertzian interaction subjected to a simple shear in the athermal quasi-static limit. The stress-strain curve is shown to separate into smooth, elastic branches followed by a subsequent plastic event. Mode analysis shows that the lowest-frequency vibrational mode is more localized, and eigenvalues and participation ratios of low- frequency modes exhibit similar power-law behavior as the system approaches plastic instability, indicating that the nature of plastic events in the granular system is also a saddle node bifurcation. The analysis of projection and spatial structure shows that over 75% contributions to the non-affine displacement field at a plastic instability come from the lowest-frequency mode, and the lowest-frequency mode is strongly spatially correlated with local plastic rearrangements, inferring that the lowest-frequency mode could be used as a predictor for future plastic rearrangements in the disordered system jammed marginally.

A Modified Quasi-Static Method for CALM System Analysis

The quasi-static analysis method introduced by API RP 2P is well known and accepted as a very useful mooring analysis method. In the early design stage, this method is widely used for preliminary analysis and mooring parameter selection. However, the quasi-static method of API RP 2P is developed for single-floating-body condition, i. e., only one floating body is considered in the computation procedure. Difficulties arise when it is used for the analysis of a CALM system, which is comprised of two floating bodies (tanker and buoy). This paper presents an analysis procedure for a two-floating-body system based on the quasi-static procedure of API RP 2P with some modifications reflecting special characteristics of the CALM system. Finally, the analysis results of a CALM system are given to illustrate the use of this procedure.

FEM and EFG Quasi-Static Explicit Buckling Analysis for Thin-Walled Members

The quasi-static explicit finite element method (FEM) and element free Galerkin (EFG) method are applied to trace the post-buckling equilibrium path of thin-walled members in this paper. The factors that primarily control the explicit buckling solutions, such as the computation time, loading function and dynamic relaxation, are investigated and suggested for the buckling analysis of thin-walled members. Three examples of different buckling modes, namely snap-through, overall and local buckling, are studied based on the implicit FEM, quasi-static explicit FEM and EFG method via the commercial software LS-DYNA. The convergence rate and accuracy of the explicit methods are compared with the conventional implicit arc-length method. It is drawn that EFG quasi-static explicit buckling analysis presents the same accurate results as implicit finite element solution, but is without convergence problem and of less-consumption of computing time than FEM.

Archetype Analysis of Limits-to-growth of Rural Tourism Development in China and Countermeasures

Based on archetype analysis of system dynamics,the limits-to-growth archetype of Peter Senge(one of his eight archetypes) was borrowed,and various influence factors of rural tourism development combined to construct the limits-to-growth archetype which consists of limits-to-growth archetype of scientific planning balance,and limits-to-growth archetype of infrastructure input balance,and limits-to-growth archetype of high-quality staff balance.Through the join of limits-to-growth sub-archetypes,limits-to-growth archetype of rural tourism development in overseas countries were obtained,on the basis of which further researches were carried out to analyze management countermeasures to eliminate limits to growth via rural tourism development.It was stressed that government should play the guiding role,unified planning,scientific management and policy support further enhanced;more efforts should be devoted to the infrastructure construction,and more financial support given by the government;more professional talents should be introduced in multiple ways,and training for present rural tourism staff enhanced.

Limit analysis of extended reach drilling in South China Sea

Extended reach wells （ERWs）, especially horizontal extended reach well with a high HD （horizontal displacement） to TVD （true vertical depth） ratio, represent a frontier technology and challenge the drilling limitations. Oil and gas reservoir in beaches or lakes and offshore can be effectively exploited by using extended reach drilling （ERD） technology. This paper focuses on the difficult technological problems encountered during exploiting the Liuhua 11-1 oil field in the South China Sea, China. Emphasis is on investigating the key subjects including prediction and control of open hole limit extension in offshore ERD, prediction of casing wear and its prevention and torque reduction, φ244.5mm casing running with floating collars to control drag force, and steerable drilling modes. The basic concept of limit extension in ERD is presented and the prediction method for open hole limit extension is given in this paper. A set of advanced drilling mechanics and control technology has been established and its practical results are verified by field cases. All those efforts may be significant for further investigating and practicing ERD limit theory and control technology in the future.

Combined influence of nonlinearity and dilation on slope stability evaluated by upper-bound limit analysis

The combined influence of nonlinearity and dilation on slope stability was evaluated using the upper-bound limit analysis theorem.The mechanism of slope collapse was analyzed by dividing it into arbitrary discrete soil blocks with the nonlinear Mohr–Coulomb failure criterion and nonassociated flow rule.The multipoint tangent(multi-tangent) technique was used to analyze the slope stability by linearizing the nonlinear failure criterion.A general expression for the slope safety factor was derived based on the virtual work principle and the strength reduction technique,and the global slope safety factor can be obtained by the optimization method of nonlinear sequential quadratic programming.The results show better agreement with previous research result when the nonlinear failure criterion reduces to a linear failure criterion or the non-associated flow rule reduces to an associated flow rule,which demonstrates the rationality of the presented method.Slope safety factors calculated by the multi-tangent inclined-slices technique were smaller than those obtained by the traditional single-tangent inclined-slices technique.The results show that the multi-tangent inclined-slices technique is a safe and effective method of slope stability limit analysis.The combined effect of nonlinearity and dilation on slope stability was analyzed,and the parameter analysis indicates that nonlinearity and dilation have significant influence on the result of slope stability analysis.

RIGID FINITE ELEMENT AND LIMIT ANALYSIS

According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.

Undrained stability of a spherical cavity in cohesive soils using finite element limit analysis

A parametric study of undrained stability of a spherical cavity in clays is investigated by finite elementlimit analysis with an axisymmetric condition. Influences of cover depth ratio of cavity and dimensionlessoverburden factor on predicted failure mechanisms and dimensionless load factor are examined.It is found that a previously recommended and up-to-date lower bound solution to the problem wassignificantly inaccurate for practice use. Thus, an accurate approximate solution to the problem is proposedfrom nonlinear regression analysis of the computed average bound solutions. New cavity stabilityfactors for the soil cohesion and soil unit weight are proposed. New findings are revealed for the threedimensionaleffect of the cavity shape on these factors between the axisymmetric and plane strainconditions, and their applications to the undrained stability evaluation of cavity problems in practice aredescribed.

PLASTIC LIMIT ANALYSIS OF DUCTILE COMPOSITE STRUCTURES FROM MICRO-TO MACRO-MECHANICAL ANALYSIS

The load-bearing capacity of ductile composite structures comprised of periodic composites is studied by a combined micro/macromechanicai approach. Firstly, on the microscopic level, a representative volume element （RVE） is selected to reflect the microstructures of the composite materials and the constituents are assumed to be elastic perfectly-plastic. Based on the homogenization theory and the static limit theorem, an optimization formulation to directly calculate the macroscopic strength domain of the RVE is obtained. The finite element modeling of the static limit analysis is formulated as a nonlinear mathematical programming and solved by the sequential quadratic programming method, where the temperature parameter method is used to construct the self-stress field. Secondly, Hill＇s yield criterion is adopted to connect the micromechanicai and macromechanical analyses. And the limit loads of composite structures are worked out on the macroscopic scale. Finally, some examples and comparisons are shown.

A pseudo-coupled numerical approach for stability analysis of frozen soil slopes based on finite element limit analysis method

To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on the computed thermal-hydro field, the stability of frozen soil slope is evaluated. Although the shear strength for frozen soil is very complicated and is usually represented by a nonlinear MC failure criterion, a simple linear MC yield criterion is utilized. In this method, the internal friction angle is expressed as a function of volumetric ice content and the cohesion is fitted as a simple bilinear expression of Tand volumetric water content. To assess slope stability, the limit analysis is employed in conjunction with the recently developed a-section search algorithm. A frozen soil slope example is used to examine the proposed pseudo-coupled numerical approach, and numerical studies validate its effectiveness. Based on numerical results, it is seen that slope stability may be remarkably influenced by warming air （or grotmd surface） temperature. With increasing ground surface temperature, slope stability indicated by FOS may reduce to 1.0, implying that wanning air temperature could be a trigger of frozen soil slope failure.

Three-dimensional stability analysis of anisotropic and non-homogeneous slopes using limit analysis

A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modification and extension of the solutions. The influences of friction angle, anisotropy factor, nonhomogeneous factor, slope angle, ratio of width to depth, and load on the slope crest were investigated. The results show that solutions are suitable to deal with the purely cohesive soils and frictional/cohesive soils, isotropic and anisotropic, homogeneous and nonhomogeneous, loaded and unloaded cases.

Safe range analysis of clear distance of twin shallow tunnels based on limit analysis and reliability theory

A method of analyzing the stability of twin shallow tunnels was presented using both limit analysis with nonlinear failure criterion and reliability theory.In the condition of nonlinear failure criterion,the critical clear distancesof twin shallow tunnels were obtained by analyzing the change of surrounding pressure.A reliability model was established based on limit state equation,and the failure probability was solved by virtue of Monte Carlo method.Safety factor and corresponding clear distance of different safety levels were obtained by introducing a target reliability index.The scope of clear distance for different safety levels is described,which can be used as a supplement and improvement to the design codes of tunnels.

Force analysis of pile foundation in rock slope based on upper-bound theorem of limit

Based on the characteristic that the potential sliding surfaces of rock slope are commonly in the shape of either line or fold line,analysis thought of conventional pile foundation in the flat ground under complex load condition was applied and the upper-bound theorem of limit analysis was used to compute thrust of rock layers with all possible distribution shapes. The interaction of slope and pile was considered design load in terms of slope thrust,and the finite difference method was derived to calculate inner-force and displacement of bridge pile foundation in rock slope under complex load condition. The result of example shows that the distribution model of slope thrust has certain impact on displacement and inner-force of bridge pile foundation. The maximum displacement growth rate reaches 54% and the maximum moment and shear growth rates reach only 15% and 20%,respectively,but the trends of inner-force and displacement of bridge pile foundation are basically the same as those of the conventional pile foundation in the flat ground. When the piles bear the same level lateral thrust,the distribution shapes of slope thrust have different influence on inner-force of pile foundation,especially the rectangle distribution,and the triangle thrust has the smallest displacement and inner-force of pile foundation.

Predicting seismic permanent displacement of soil walls under surcharge based on limit analysis approach

Seismic permanent displacement of the soil walls plays an important role in design of these structures. Due to the increase in growth of urban areas and the limitations in use of flat grounds, many structures are built near slopes and retaining walls. During earthquakes, these structures can apply an additional surcharge on the wall. The intensity and location of the surcharge is of considerable importance on the seismic displacements of the soil wall. In this study, by using the limit analysis and upper bound theorem, seismic permanent displacement of the soil wall under surcharge has been analyzed. Thus, a formulation is presented for calculating the yield acceleration and seismic displacement for different surcharge conditions. The effect of seismic acceleration, surcharge intensity, its location and soil properties is investigated. A parameter called the ＂displacement coefficient＂ is proposed, and is a potential modification for Newmark’s sliding-block method.

A LOWER BOUND LIMIT ANALYSIS OF DUCTILE COMPOSITE MATERIALS

The plastic load-bearing capacity of ductile composites such as metal matrix composites is studied with an insight into the microstructures. The macroscopic strength of a composite is obtained by combining the homogenization theory with static limit analysis, where the temperature parameter method is used to construct the self-equilibrium stress field. An interface failure model is proposed to account for the effects of the interface on the failure of composites. The static limit analysis with the finite-element method is then formulated as a constrained nonlinear programming problem, which is solved by the Sequential Quadratic Programming （SQP） method. Finally, the macroscopic transverse strength of perforated materials, the macroscopic transverse and off-axis strength of fiber-reinforced composites are obtained through numerical calculation. The computational results are in good agreement with the experimental data.

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead ofthe finite element method is used t perform lower bound limit andshakedown analysis of structures. The self-equilibrium stress fieldsare constructed by a linear combination of several basicself-equilibrium stress fields with parameters to be determined.These basic self-equilibrium stress fields are expressed as elasticresponses of the body to im- posed permanent strains and obtainedthrough elastic-plastic incremental analysis.

Upper-bound limit analysis based on the natural element method

The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic the- orem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nons- moothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.

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.