Influence of joint spacing and rock characteristics on the toppling stability of cut rock slope through a simplified limit equilibrium method

Toppling failure of rock mass/soil slope is an important geological and environmental problem.Clarifying its failure mechanism under different conditions has great significance in engineering.The toppling failure of a cutting slope occurred in a hydropower station in Kyushu,Japan illustrates that the joint characteristic played a significant role in the occurrence of rock slope tipping failure.Thus,in order to consider the mechanical properties of jointed rock mass and the influence of geometric conditions,a simplified analytical approach based on the limit equilibrium method for modeling the flexural toppling of cut rock slopes is proposed to consider the influence of the mechanical properties and geometry condition of jointed rock mass.The theoretical solution is compared with the numerical solution taking Kyushu Hydropower Station in Japan as one case,and it is found that the theoretical solution obtained by the simplified analysis method is consistent with the numerical analytical solution,thus verifying the accuracy of the simplified method.Meanwhile,the Goodman-Bray approach conventionally used in engineering practice is improved according to the analytical results.The results show that the allowable slope angle may be obtained by the improved Goodman-Bray approach considering the joint spacing,the joint frictional angle and the tensile strength of rock mass together.

Influence of confined water on the limit support pressure of tunnel face in weakly water-rich strata

In the process of shield tunneling through soft soil layers,the presence of confined water ahead poses a significant threat to the stability of the tunnel face.Therefore,it is crucial to consider the impact of confined water on the limit support pressure of the tunnel face.This study employed the finite element method(FEM)to analyze the limit support pressure of shield tunnel face instability within a pressurized water-containing layer.Subsequently,a multiple linear regression approach was applied to derive a concise solution formula for the limit support pressure,incorporating various influencing factors.The analysis yields the following conclusions:1)The influence of confined water on the instability mode of the tunnel face in soft soil layers makes the displacement response of the strata not significant when the face is unstable;2)The limit support pressure increases approximately linearly with the pressure head,shield tunnel diameter,and tunnel burial depth.And inversely proportional to the thickness of the impermeable layer,soil cohesion and internal friction angle;3)Through an engineering case study analysis,the results align well with those obtained from traditional theoretical methods,thereby validating the rationality of the equations proposed in this paper.Furthermore,the proposed equations overcome the limitation of traditional theoretical approaches considering the influence of changes in impermeable layer thickness.It can accurately depict the dynamic variation in the required limit support pressure to maintain the stability of the tunnel face during shield tunneling,thus better reflecting engineering reality.

Probabilistic Analysis of Slope Using Finite Element Approach and Limit Equilibrium Approach around Amalpata Landslide of West Central, Nepal

The stability study of the ongoing and recurring Amalpata landslide in Baglung in Nepal’s Gandaki Province is presented in this research. The impacted slope is around 200 meters high, with two terraces that have different slope inclinations. The lower bench, located above the basement, consistently fails and sets others up for failure. The fluctuating water level of the slope, which travels down the slope masses, exacerbates the slide problem. The majority of these rocks are Amalpata landslide area experiences several structural disruptions. The area’s stability must be evaluated in order to prevent and control more harm from occurring to the nearby agricultural land and people living along the slope. The slopes’ failures increase the damages of house existing in nearby area and the erosion of the slope. Two modeling techniques the finite element approach and the limit equilibrium method were used to simulate the slope. The findings show that, in every case, the terrace above the basement is where the majority of the stress is concentrated, with a safety factor of near unity. Using probabilistic slope stability analysis, the failure probability was predicted to be between 98.90% and 100%.

Overhanging rock slope by design:An integrated approach using rock mass strength characterisation,large-scale numerical modelling and limit equilibrium methods

Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.

Lateral Bearing Capacity of Modified Suction Caissons Determined by Using the Limit Equilibrium Method

The modified suction caisson(MSC) adds a short-skirted structure around the regular suction caissons to increase the lateral bearing capacity and limit the deflection. The MSC is suitable for acting as the offshore wind turbine foundation subjected to larger lateral loads compared with the imposed vertical loads. Determination of the lateral bearing capacity is a key issue for the MSC design. The formula estimating the lateral bearing capacity of the MSC was proposed in terms of the limit equilibrium method and was verified by the test results. Parametric studies on the lateral bearing capacity were also carried out. It was found that the lateral bearing capacity of the MSC increases with the increasing length and radius of the external skirt, and the lateral bearing capacity increases linearly with the increasing coefficient of subgrade reaction. The maximum lateral bearing capacity of the MSC is attained when the ratio of the radii of the internal compartment to the external skirt equals 0.82 and the ratio of the lengths of the external skirt to the internal compartment equals 0.48, provided that the steel usage of the MSC is kept constant.

A Generalized Limit Equilibrium Method for the Solution of Active Earth Pressure on a Retaining Wall

In this paper, a generalized limit equilibrium method of solving the active earth pressure problem behind a retaining wall is proposed.Differing from other limit equilibrium methods, an arbitrary slip surface shape without any assumptions of pre-defined shapes is needed in the current framework, which is verified to find the most probable failure slip surface. Based on the current computational framework, numerical comparisons with experiment, discrete element method and other methods are carried out. In addition, the influences of the inclination of the wall, the soil cohesion, the angle of the internal friction of the soil, the slope inclination of the backfill soil on the critical pressure coefficient of the soil, the point of application of the resultant earth pressure and the shape of the slip surface are also carefully investigated. The results demonstrate that limit equilibrium solution from predefined slip plane assumption, including Coulomb solution, is a special case of current computational framework. It is well illustrated that the current method is feasible to evaluate the characteristics of earth pressure problem.

Comprehensive analysis of slope stability and determination of stable slopes in the Chador-Malu iron ore mine using numerical and limit equilibrium methods

One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was considered vital to perform a comprehensive slope stability analysis. At first, we divided the existing rock hosting pit into six zones and a geotechnical map was prepared. Then, the value of MRMR （Mining Rock Mass Rating） was determined for each zone. Owing to the fact that the Chador-Malu iron ore mine is located in a highly tectonic area and the rock mass completely crushed, the Hoek-Brown failure criterion was found suitable to estimate geo-mechanical parameters. After that, the value of cohesion （c） and friction angle （tp） were calculated for different geotechnical zones and relative graphs and equations were derived as a function of slope height. The stability analyses using numerical and limit equilibrium methods showed that some instability problems might occur by increasing the slope height. Therefore, stable slopes for each geotechnical zone and prepared sections were calculated and presented as a function of slope height.

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.

THE APPLICATION OF WEINSTEIN-CHIEN′S METHOD——THE UPPER AND LOWER LIMITS OF FUNDAMENTAL FREQUENCY OF RECTANGULAR PLATES WITH EDGES ARE THE MIXTURE OF SIMPLY SUPPORTED PORTIONS AND CLAMPED PORTIONS

In this paper, the method of relaxed boundary conditions is applied to rectangular plates with edges which are a sort of the mixture of simply supported portions and clamped portions, so that the lower limit of fundamental frequency of such plates is evaluated. A kind of polynomial satisfying the displacement boundary conditions is designed, os that it is enabled to evaluate the upper limit of fundamental frequency by Ritz' method. The practical calculation examples solved by these methods have given satisfactory results. At the end of this paper, it is pointed out that the socalled exact solution of such plates usually evaluated by the force superposition method is essentially a kind of lower limit of solution, if the truncated error of series which occurs in actual calculation is considered.

Stability analysis of tunnel face reinforced with face bolts

Face bolting has been widely utilized to enhance the stability of tunnel face,particularly in soft soil tunnels.However,the influence of bolt reinforcement and its layout on tunnel face stability has not been systematically studied.Based on the theory of linear elastic mechanics,this study delved into the specific mechanisms of bolt reinforcement on the tunnel face in both horizontal and vertical dimensions.It also identified the primary failure types of bolts.Additionally,a design approach for tunnel face bolts that incorporates spatial layout was established using the limit equilibrium method to enhance the conventional wedge-prism model.The proposed model was subsequently validated through various means,and the specific influence of relevant bolt design parameters on tunnel face stability was analyzed.Furthermore,design principles for tunnel face bolts under different geological conditions were presented.The findings indicate that bolt failure can be categorized into three stages:tensile failure,pullout failure,and comprehensive failure.Increasing cohesion,internal friction angle,bolt density,and overlap length can effectively enhance tunnel face stability.Due to significant variations in stratum conditions,tailored design approaches based on specific failure stages are necessary for bolt design.

A new method for the stability analysis of geosynthetic-reinforced slopes

This paper is concerned with the stability analysis of reinforced slopes.A new approach based on the limit equilibrium principle is proposed to evaluate the stability of the reinforced slopes.The effect of reinforcement is modeled as an equivalent restoring force acting the bottom of the slice and added into the general limit equilibrium(GLE) method.The equations of force and moment equilibrium of the slice are derived and corresponding iterative solution methods are provided.The new method can satisfy both the force and the moment equilibrium and be applicable to the critical failure surface of arbitrary form.Furthermore,the results predicted by the proposed method are compared with the calculation examples of other researchers and the centrifuge model test results to validate its correctness and effectiveness.

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.

Hybrid method of limit equilibrium and finite element internal force for analysis of arch dam stability against sliding

A hybrid method of limit equilibrium and finite element internal force for analysis of arch dam stability against sliding is presented.The finite element internal force method(FEIFM) is used to provide more accurate thrust forces acting on the faces of a slip body,and the limit equilibrium method(LEM) is employed to evaluate the factor of safety of the slip body.The method presented can deal with a slip body with large amount of geometrically complex slip faces.In addition,compared with the traditional LEM,it can meet the balance condition of the forces in the slip faces.An example shows that the factor of safety obtained by the method presented agrees well with the theoretical solution.A practical example is also presented to demonstrate the application of the method in the stability analysis of an arch dam project.The results from the examples show that the method is promising in analysis of arch dam stability against sliding.

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.

Deformation mechanism and limit support pressure of cutting steel plate during connection between pipes in large spacing using pipe curtain structure method

The pipe curtain structure method(PSM)is a novel construction method to control ground deformation strictly.Compared with the traditional pipe-roofing and pipe jacking method,the connection between pipes in large spacings using PSM is widely acknowledged as a unique construction procedure.Further study on this connection procedure is needed to resolve similar cases in that the pipes are inevitably constructed on both sides of existing piles.Cutting the steel plate during the connection procedure is the first step,which is crucial to control the safety and stability of the surrounding environment and existing structures.The deformation mechanism and limit support pressure of the cutting steel plate during the connection between pipes in large spacings are studied in this paper,relying on the undercrossing Yifeng gate tower project of Jianning West Road River Crossing Channel in Nanjing,China.A modified 3D wedge-prism failure model is proposed using the 3D discrete element method.Combined with Terzaghi loose earth pressure theory and the limit equilibrium theory,the analytical solutions for the limit support pressure of the excavation face of the cutting steel plate are derived.The modified 3D wedge-prism failure model and corresponding analytical solutions are categorised into two cases:(a)unilateral cutting scheme,and(b)bilateral cutting scheme.The analytical solutions for the two cases are verified from the numerical simulation and in-situ data and compared with the previous solutions.The comparative analysis between the unilateral and bilateral cutting schemes indicates that the bilateral cutting scheme can be adopted as a priority.The bilateral cutting scheme saves more time and induces less ground deformation than the unilateral one due to the resistance generated from the superimposed wedge.In addition,the parametric sensitivity analysis is carried out using an orthogonal experimental design.The main influencing factors arranged from high to low are the pipe spacing,the cutting size,and the pipe burial depth.The ground deformation increases with the increased cutting size and pipe spacing.The pipe burial depth slightly affects the ground deformation if the other two factors are minor.Cutting steel plates in small sizes,excavating soil under low disturbance,and supporting pipes for high frequency can effectively reduce the ground surface subsidence.

Graphical method based on modified maximum force criterion to indicate forming limit curves of 22MnB5 boron steel sheets at elevated temperatures

A new approach for predicting forming limit curves(FLCs)at elevated temperatures was proposed herein.FLCs are often used to predict failure and determine the optimal forming parameters of automotive parts.First,a graphical method based on a modified maximum force criterion was applied to estimate the FLCs of 22MnB5 boron steel sheets at room temperature using various hardening laws.Subsequently,the predicted FLC data at room temperature were compared with corresponding data obtained from Nakazima's tests to obtain the best prediction.To estimate the FLC at elevated temperatures,tensile tests were conducted at various temperatures to determine the ratios of equivalent fracture strains between the corresponding elevated temperatures and room temperature.FLCs at elevated temperatures could be established based on obtained ratios.However,the predicted FLCs at elevated temperatures did not agree well with the corresponding FLC experimental data of Zhou et al.A new method was proposed herein to improve the prediction of FLCs at elevated temperatures.An FLC calculated at room tem-perature was utilized to predict the failure of Nakazima's samples via finite element simulation.Based on the simulation results at room temperature,the mathematical relationships between the equivalent ductile fracture strain versus stress triaxiality and strain ratio were established and then combined with ratios between elevated and room temperatures to calculate the FLCs at different temperatures.The predicted FLCs at elevated temperatures agree well with the corresponding experimental FLC data.

A new double reduction method for slope stability analysis

The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a detailed calculation procedure and a definition of FOS for slope stability was developed based on the understanding of SRM. When constructing the new definition of FOS, efforts were made to make sure that it has concise physical meanings and fully reflects the shear strength of the slope. Two examples, slopes A and B with the slope angles of 63° and 34° respectively, were given to verify the method presented. It is found that, for these two slopes, the FOSs from original strength reduction method are respectively 1.5% and 38% higher than those from double reduction method. It is also found that the double reduction method predicts a deeper potential slide line and a larger slide mass. These results show that on one hand, the double reduction method is comparative to the traditional methods and is reasonable, and on the other hand, the original strength reduction method may overestimate the safety of a slope. The method presented is advised to be considered as an additional option in the practical slope stability evaluations although more useful experience is required.

Seismic passive earth resistance using modified pseudo-dynamic method

In earthquake prone areas, understanding of the seismic passive earth resistance is very important for the design of different geotechnical earth retaining structures. In this study, the limit equilibrium method is used for estimation of critical seismic passive earth resistance for an inclined wall supporting horizontal cohesionless backfill. A composite failure surface is considered in the present analysis. Seismic forces are computed assuming the backfill soil as a viscoelastic material overlying a rigid stratum and the rigid stratum is subjected to a harmonic shaking. The present method satisfies the boundary conditions. The amplification of acceleration depends on the properties of the backfill soil and on the characteristics of the input motion. The acceleration distribution along the depth of the backfill is found to be nonlinear in nature. The present study shows that the horizontal and vertical acceleration distribution in the backfill soil is not always in-phase for the critical value of the seismic passive earth pressure coefficient. The effect of different parameters on the seismic passive earth pressure is studied in detail. A comparison of the present method with other theories is also presented, which shows the merits of the present study.

IMPROVEMENTS ON THE ARC-LENGTH-TYPE METHOD

Arc-length-type and energy-type methods are two main strategies used in structural nonlinear tracing analysis, but the former is widely used due to the explicitness and clarity in conception, as well as the convenience and reliability in calculation. It is very important to trace the complete load-deflection path in order to know comprehensively the characteristics of structures subjected to loads. Unfortunately, the nonlinear analysis techniques are only workable for tracing the limit-point-type equilibrium path. For the bifurcation-point-type path, most of them cannot secure a satisfactory result. In this paper, main arc-length-type methods are reviewed and compared, and the possible reasons of failures in tracing analysis are briefly discussed. Some improvements are proposed, a displacement perturbation method and a force perturbation method are presented for tracing the bifurcation-point-type paths. Finally, two examples are analyzed to verify the ideas, and some conclusions are drawn with respect to the arc-length-type methods.

The prediction of external flow field and hydrodynamic force with limited data using deep neural network

Predicting the external flow field with limited data or limited measurements has attracted long-time interests of researchers in many industrial applications.Physics informed neural network(PINN)provides a seamless framework for combining the measured data with the deep neural network,making the neural network capable of executing certain physical constraints.Unlike the data-driven model to learn the end-to-end mapping between the sensor data and high-dimensional flow field,PINN need no prior high-dimensional field as the training dataset and can construct the mapping from sensor data to high dimensional flow field directly.However,the extrapolation of the flow field in the temporal direction is limited due to the lack of training data.Therefore,we apply the long short-term memory(LSTM)network and physics-informed neural network(PINN)to predict the flow field and hydrodynamic force in the future temporal domain with limited data measured in the spatial domain.The physical constraints(conservation laws of fluid flow,e.g.,Navier-Stokes equations)are embedded into the loss function to enforce the trained neural network to capture some latent physical relation between the output fluid parameters and input tempo-spatial parameters.The sparsely measured points in this work are obtained from computational fluid dynamics(CFD)solver based on the local radial basis function(RBF)method.Different numbers of spatial measured points(4–35)downstream the cylinder are trained with/without the prior knowledge of Reynolds number to validate the availability and accuracy of the proposed approach.More practical applications of flow field prediction can compute the drag and lift force along with the cylinder,while different geometry shapes are taken into account.By comparing the flow field reconstruction and force prediction with CFD results,the proposed approach produces a comparable level of accuracy while significantly fewer data in the spatial domain is needed.The numerical results demonstrate that the proposed approach with a specific deep neural network configuration is of great potential for emerging cases where the measured data are often limited.