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.

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.

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.

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.

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.

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

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.

Lower Bound Limit Analysis of Anisotropic Soils

Previous approaches can only tackle anisotropic problems with cohesion varying with direction.A novel linearization of the Mohr-Coulomb yield criterion associated with plane strain problem has been achieved by simulating the Mohr’s circle with orientation lines inσ-τspace,which allows for lower bound solution of soils with cohesion and friction coefficient varying with direction.The finite element lower limit analysis formulation using the modified anisotropic yield criterion is then developed.Several examples are given to illustrate the capability and effectiveness of the proposed numerical procedure for computing rigorous lower bounds for anisotropic soils.

Upper Bound Limit Analysis of Anisotropic Soils

In this paper,a novel discretization method inσ-τspace is developed to calculate the upper bound limit loads and failure modes of anisotropic Mohr-Coulomb materials.To achieve this objective,the Mohr-Coulomb yield criterion is linearized inσ-τspace,which allows for upper bound solution of soils whose cohesion and friction coefficient varying with direction.The finite element upper limit analysis formulation using the modified anisotropic yield criterion is then developed.Several examples are given to illustrate the capability and effectiveness of the proposed numerical procedure for computing rigorous upper bounds for anisotropic soils.

PLASTIC LIMIT ANALYSIS OF INCOMPATIBLE FINITE ELEMENT METHOD

This paper describes an incompatible finite element model satisfying the consistency condition of energy to solve the numerical precision problem of finite element solution in perfectly plastic analysis. In this paper the reason and criterion of the application of the model to plastic limit analysis are discussed, and an algorithm of computing plastic limit load is given.

Upper-Bound Limit Analysis of the Multi-Layer Slope Stability and Failure Mode Based on Generalized Horizontal Slice Method

Multi-layer slopes are widely found in clay residue receiving fields.A generalized horizontal slice method(GHSM)for assessing the stability of multi-layer slopes that considers the energy dissipation between adjacent horizontal slices is presented.In view of the upper-bound limit analysis theory,the energy equation is derived and the ultimate failure mode is generated by comparing the sliding surface passing through the slope toe(mode A)with that below(mode B).In addition,the influence of the number of slices on the stability coefficients in the GHSM is studied and the stable value is obtained.Compared to the original method(Chen’s method),the GHSM can acquire more precise results,which takes into account the energy dissipation in the inner sliding soil mass.Moreover,the GHSM,limit equilibrium method(LEM)and numerical simulation method(NSM)are applied to analyze the stability of a multi-layer slope with different slope angles and the results of the safety factor and failure mode are very close in each case.The ultimate failure modes are shown to be mode B when the slope angle is not more than 28°.It illustrates that the determination of the ultimate sliding surface requires comparison of multiple failure modes,not only mode A.

Characterizing model uncertainty of upper-bound limit analysis on slopes using 3D rotational failure mechanism

The slope stability assessment is a classical problem in geotechnical engineering.This topic have attracted many researcher’s attention and various theoretical models for predicting critical slope heights or safety factors in the light of the limit equilibrium(LE)method and the kinematical approach of limit analysis(LA)method.Meanwhile,a large number of experimental studies have been conducted to check the slope stability.Using centrifuge testing results,this paper aims to employ Bayesian method to characterize the model uncertainties of the classical three-dimensional rotational failure mechanism proposed by Michalowski and Drescher(2009)to predict critical slope heights in frictional soils,by incorporating the test uncertainties and parameter uncertainties.The obtained results show that the LA three-dimensional rotational failure mechanism overestimates the critical slope height compared with the LE method,and the experimental observational uncertainty has negligible influences on the posterior statistics of model uncertainty.

Inverse reliability analysis and design for tunnel face stability considering soil spatial variability

The traditional deterministic analysis for tunnel face stability neglects the uncertainties of geotechnical parameters,while the simplified reliability analysis which models the potential uncertainties by means of random variables usually fails to account for soil spatial variability.To overcome these limitations,this study proposes an efficient framework for conducting reliability analysis and reliability-based design(RBD)of tunnel face stability in spatially variable soil strata.The three-dimensional(3D)rotational failure mechanism of the tunnel face is extended to account for the soil spatial variability,and a probabilistic framework is established by coupling the extended mechanism with the improved Hasofer-Lind-Rackwits-Fiessler recursive algorithm(iHLRF)as well as its inverse analysis formulation.The proposed framework allows for rapid and precise reliability analysis and RBD of tunnel face stability.To demonstrate the feasibility and efficacy of the proposed framework,an illustrative case of tunnelling in frictional soils is presented,where the soil's cohesion and friction angle are modelled as two anisotropic cross-correlated lognormal random fields.The results show that the proposed method can accurately estimate the failure probability(or reliability index)regarding the tunnel face stability and can efficiently determine the required supporting pressure for a target reliability index with soil spatial variability being taken into account.Furthermore,this study reveals the impact of various factors on the support pressure,including coefficient of variation,cross-correlation between cohesion and friction angle,as well as autocorrelation distance of spatially variable soil strata.The results also demonstrate the feasibility of using the forward and/or inverse first-order reliability method(FORM)in high-dimensional stochastic problems.It is hoped that this study may provide a practical and reliable framework for determining the stability of tunnels in complex soil strata.

Limit load and failure mechanisms of a vertical Hoek-Brown rock slope

The problem considered in this short note is the limit load determination of a vertical rock slope.The classical limit theorem is employed with the use of adaptive finite elements and nonlinear programming to determine upper and lower bound limit loads of a Hoek-Brown vertical rock slope.The objective function of the mathematical programming problem is such as to optimize a boundary load,which is known as the limit load,resembling the ultimate bearing capacity of a strip footing.While focusing on the vertical slope,parametric studies are carried out for several dimensionless ratios such as the dimensionless footing distance ratio,the dimensionless height ratio,and the dimensionless rock strength ratio.A comprehensive set of design charts is presented,and failure envelopes shown with the results explained in terms of three identified failure mechanisms,i.e.the face,the toe,and the Prandtl-type failures.These novel results can be used with great confidence in design practice,in particularly noting that the current industry-based design procedures for the presented problem are rarely found.

Limit analysis on seismic stability of anisotropic and nonhomogeneous slopes with anti-slide piles

This study employs the limit analysis method to evaluate the seismic stability of anisotropic and nonhomogeneous slopes stabilized with anti-slide piles. The pseudo-static approach is used to simplify the earthquake load. The yield seismic acceleration factor is obtained from the optimization procedure and the results are verified with the published data. Then, the seismically-unstable slope is reinforced with anti-slide piles, and the seismic stability of the reinforced slope is explored. The results show that the anisotropy and nonhomogeneity of soils have significant effects on the stabilizing force required from the anti-slide piles and the optimal location of the pile is near the toe of the slope.