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.

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.

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.

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.

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.

A Bayesian multi-model inference methodology for imprecise momentindependent global sensitivity analysis of rock structures

Traditional global sensitivity analysis(GSA)neglects the epistemic uncertainties associated with the probabilistic characteristics(i.e.type of distribution type and its parameters)of input rock properties emanating due to the small size of datasets while mapping the relative importance of properties to the model response.This paper proposes an augmented Bayesian multi-model inference(BMMI)coupled with GSA methodology(BMMI-GSA)to address this issue by estimating the imprecision in the momentindependent sensitivity indices of rock structures arising from the small size of input data.The methodology employs BMMI to quantify the epistemic uncertainties associated with model type and parameters of input properties.The estimated uncertainties are propagated in estimating imprecision in moment-independent Borgonovo’s indices by employing a reweighting approach on candidate probabilistic models.The proposed methodology is showcased for a rock slope prone to stress-controlled failure in the Himalayan region of India.The proposed methodology was superior to the conventional GSA(neglects all epistemic uncertainties)and Bayesian coupled GSA(B-GSA)(neglects model uncertainty)due to its capability to incorporate the uncertainties in both model type and parameters of properties.Imprecise Borgonovo’s indices estimated via proposed methodology provide the confidence intervals of the sensitivity indices instead of their fixed-point estimates,which makes the user more informed in the data collection efforts.Analyses performed with the varying sample sizes suggested that the uncertainties in sensitivity indices reduce significantly with the increasing sample sizes.The accurate importance ranking of properties was only possible via samples of large sizes.Further,the impact of the prior knowledge in terms of prior ranges and distributions was significant;hence,any related assumption should be made carefully.

Change Point Analysis to Detect the Effect of Pruning Severity on Tree Growth

The effect of pruning severity on tree growth was analyzed by change point detection using segmented regression. The present study applied this analysis to a well-known published data set including diameter growth response, tree age, pruning severity and pretreatment crown size. First, multiple regression analysis was performed to assess the effect of tree age, pruning severity and pretreatment crown size on diameter growth response. Next, segmented regression analysis was performed to assess the effect of pruning severity on diameter growth response. The results of the multiple regression showed that diameter growth response was significantly influenced by pruning severity and pretreatment crown size. The results of the segmented regression showed that in the whole data set, an abrupt change toward a decrease in diameter growth response was detected at 25% of the live crown removed. However, in the group of fully crowned and open-grown, diameter growth response continuously decreased with increasing pruning severity with no significant abrupt change, whereas in the group of 70% - 90% live crown, diameter growth response did not significantly decrease up to the break point (53% crown removed) and then abruptly decreased. This may be the first study to show the numerical evaluation of the effect of pruning severity on tree growth by change point analysis.

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.

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.

Multivariate adaptive regression splines analysis for 3D slope stability in anisotropic and heterogenous clay

Little research can be found in relation to the stability of anisotropic and heterogenous soils in three dimensions.In this paper,we propose a study on the three-dimensional(3D)undrained slopes in anisotropic and heterogenous clay using advanced upper and lower bounds finite element limit analysis(FELA).The obtained stability solutions are normalized,and presented by a stability number that is a function of three geometrical ratios and two material ratios,i.e.depth ratio,length ratio,slope angle,shear strength gradient ratio and anisotropic strength ratio.Numerical results are compared with experimental data in the literature,and charts are presented to cover a wide range of design parameters.Using the multivariate adaptive regression splines(MARS)analysis,the respective influence and sensitivity of each design parameter on the stability number and the failure mechanism are investigated.An empirical equation is also developed to effectively estimate the stability number.

Two-dimensional face stability analysis in rock masses governed by the Hoek-Brown strength criterion with a new multi-horn mechanism

The face stability problem is a major concern for tunnels excavated in rock masses governed by the Hoek-Brown strength criterion.To provide an accurate prediction for the theoretical solution of the critical face pressure,this study adopts the piecewise linear method(PLM)to account for the nonlinearity of the strength envelope and proposes a new multi-horn rotational mechanism based on the Hoek-Brown strength criterion and the associative flow rule.The analytical solution of critical support pressure is derived from the energy-work balance equation in the framework of the plastic limit theorem;it is formulated as a multivariable nonlinear optimization problem relying on 2m dependent variables(m is the number of segments).Meanwhile,two classic linearized measures,the generalized tangential technique(GTT)and equivalent Mohr-Coulomb parameters method(EMM),are incorporated into the analysis for comparison.Surprisingly,the parametric study indicates a significant improvement in support pressure by up to 13%compared with the GTT,and as expected,the stability of the tunnel face is greatly influenced by the rock strength parameters.The stress distribution on the rupture surface is calculated to gain an intuitive understanding of the failure at the limit state.Although the limit analysis is incapable of calculating the true stress distribution in rock masses,a rough approximation of the stress vector on the rupture surface is permitted.In the end,sets of normalized face pressure are provided in the form of charts for a quick assessment of face stability in rock masses.

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.

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.

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.