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.

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.

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.

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.

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.

Three-dimensional upper bound limit analysis of underground cavities using nonlinear Baker failure criterion

A generalized nonlinear Baker failure criterion is employed with the upper bound limit analysis to study the surrounding rock stability of underground cavities. A three-dimensional(3D) failure mode is established by extending the two-dimensional(2D) failure mode, which offers an upper bound expression of the surrounding rock pressure. This method is validated with a series of examples before the influence of four parameters of scale parameter, curvature parameter, shift parameter and lateral pressure coefficient, on the surrounding rock pressure is analyzed. According to these results, failure ranges of the underground cavities are determined. The following conclusions are reached:(1) the proposed approach is more accurate to predict surrounding rock pressure than the Mohr-Coulomb failure criterion;(2) the surrounding rock with large scale parameter, curvature parameter, shift parameter, and lateral pressure coefficient can lead to a more stable underground cavity;(3) the failure range in 3D mode can be predicted according to the upper bound solutions.

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.

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.

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 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.

Limit analysis of ultimate uplift capacity and failure mechanism ofshallow plate anchors in multi-layered soils

Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of shallow strips and circular plate anchors in multi-layered soils.The nonlinear strength criterion and non-associated flow rule of geotechnical materials were introduced to investigate the influence of nonuniformity on the pullout performance and failure mechanism of shallow plate anchors.The expressions of the detaching curves or surfaces were obtained to reflect the failure mechanism,which can be used to figure out the ultimate uplift capacity and failure range.The results are generally in agreement with the numerical simulations and previous research.The effects of various parameters on the ultimate uplift capacity and failure mechanism of plate anchors in multi-layered soils were investigated,and it is found that the ultimate uplift capacity and failure range of shallow anchors increase with the increase of initial cohesion and dilatancy coefficient,but decrease with the unit weight,axial tensile stress and nonlinear coefficient.

The limit analysis in soil and rock:a mature discipline of geomechanics

The solution of a slope stability problem can be approached by its least upper-bound and maximum lower-bound with high accuracy. The limit equilibrium methods that employ vertical slices imply a lower bound of the factor of safety. It has been successfully extended to the area of active earth pressure analysis that accounts for different input of locations of earth pressure applications. Those methods that employ slices with inclined interfaces give an upper-bound approach to the stability analysis. It enjoys a sound mechanical background and is able to provide accurate solutions of soil plasticity. It has been successfully extended to the area of bearing capacity analysis in which various empirical coefficients are no longer necessary. The 3D upper- and lower-bound methods under this framework have been made possible and show great potential for solving various engineering problems.

LOWER BOUND LIMIT ANALYSIS OF THREE-DIMENSIONAL ELASTOPLASTIC STRUCTURES BY BOUNDARY ELEMENT METHOD

Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.

Upper bound limit analysis of roof collapse in shallow tunnels with arbitrary cross sections under condition of seepage force

The analytical solutions for predicting the exact shape of collapse mechanisms in shallow tunnels with arbitrary excavation profiles were obtained by virtue of the upper bound theorem of limit analysis and variation principle according to Hoek-Brown failure criterion. The seepage force was included in the upper bound limit analysis, and it was computed from the gradient of excess pore pressure distribution. The seepage was regarded as a work rate of external force. The numerical results of roof collapse in square and circular tunnels with different rock parameters were derived and discussed, which proves to be valid in comparison with the previous work. The influences of different parameters on the shape of collapsing blocks were also discussed.

Limit analysis method for active earth pressure on laggings between stabilizing piles

Stabilizing pile is a kind of earth shoring structure frequently used in slope engineering. When the piles have cantilever segments above the ground,laggings are usually installed to avoid collapse of soil between piles. Evaluating the earth pressure acting on laggings is of great importance in design process.Since laggings are usually less stiff than piles,the lateral pressure on lagging is much closer to active earth pressure. In order to estimate the lateral earth pressure on lagging more accurately,first,a model test of cantilever stabilizing pile and lagging systems was carried out. Then,basing the experimental results a three-dimensional sliding wedge model was established. Last,the calculation process of the total active force on lagging is presented based on the kinematic approach of limit analysis. A comparison is made between the total active force on lagging calculated by the formula presented in this study and the force on a same-size rigid retaining wall obtained from Rankine's theory. It is found that the proposed method fits well with the experimental results.Parametric studies show that the total active force on lagging increases with the growth of the lagging height and the lagging clear span; while decreases asthe soil internal friction angle and soil cohesion increase.

3D limit analysis of rock slopes based on equivalent linear failure criterion with tension cut-off

Hoek-Brown (HB) failure criterion is widely used to predict the strength of intact or heavily jointed rock mass. For stability analysis of rock slopes governed by the HB failure criterion, the equivalent linearity to Mohr-Coulomb (MC) criterion is often adopted, leading to the well-known equivalent Mohr-Coulomb method (EMCM). Existing studies on EMCM analysis mainly consider the shear strength of rock material, while consideration of the tensile strength is rare. This contradicts the fact that the underlying tensile strength of rock mass has considerable impact on the rock slope stability in real world. In this regard, this paper proposes a limit analysis-based approach that can account for tension in the three-dimensional (3D) stability analysis of HB rock slope. This approach is established on the equivalent linearity of the HB criterion with consideration of tensile strength, known as the equivalent tension cut-off MC method (ETMCM), and using a horn-like 3D mechanism of limit analysis. The safety factor solutions given by the proposed approach are validated by previous studies and numerical results. Parametric studies are conducted to investigate the effect of rock tensile strength on slope stability. Results show that the consideration of tension leads to a more conservative safety factor and a sharper curvature of the failure surface, and these impacts tend to be more obvious with the increases in slope inclination and slope width. Finally, the stability of the HB rock slope under seepage conditions is studied using the proposed approach. The results indicate that the effect of tensile strength is highly remarkable in seepage circumstances.

NONSMOOTH MODEL FOR PLASTIC LIMIT ANALYSIS AND ITS SMOOTHING ALGORITHM

By means of Lagrange duality of Hill＇s maximum plastic work principle theory of the convex program, a dual problem under Mises＇ yield condition has been derived and whereby a non-differentiable convex optimization model for the limit analysis is developed. With this model, it is not necessary to linearize the yield condition and its discrete form becomes a minimization problem of the sum of Euclidean norms subject to linear constraints. Aimed at resolving the non-differentiability of Euclidean norms, a smoothing algorithm for the limit analysis of perfect-plastic continuum media is proposed. Its efficiency is demonstrated by computing the limit load factor and the collapse state for some plane stress and plain strain problems.

Limit Analysis for the Mechanical Structure of the ITER Neutron Shielding Block

The ITER neutron shielding blocks are located between the inner shell and the outer shell of the vacuum vessel （VV） with the main function of providing neutron shielding. Conskicring the combined loads of the shielding blocks during the plasma operation of the ITER, limit analysis for one typical ferromagnetic （FM） shielding block has been performed and the structural design has bccn evaluated based on the American Society of Mechanical Engineers （ASME） criterion and European standards. Results show that the collapse load of this shielding block is three times the specified load, which is much higher than the design requirement of 1.25. The structure of this neutron shielding block has a sufficient safety margin.

A MATHEMATICAL PROGRAMMING ALGORITHM FOR LIMIT ANALYSIS

This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound load multiplier has been established,and an iteration algorithm proposed accordingly.In this algorithm the plastic and rigid zones are distinguished for every iteration step,and the goal function is modified gradually.The difficulties caused by the nonsmoothness of the goal function are over- come.Some examples solved by this algorithm are presented.