Three-dimensional stability calculation method for high and large composite slopes formed by mining stope and inner dump in adjacent open pits

The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.

Three-dimensional stability of landslides based on local safety factor

Unlike the limit equilibrium method(LEM), with which only the global safety factor of the landslide can be calculated, a local safety factor(LSF) method is proposed to evaluate the stability of different sections of a landslide in this paper. Based on three-dimensional(3D) numerical simulation results, the local safety factor is defined as the ratio of the shear strength of the soil at an element on the slip zone to the shear stress parallel to the sliding direction at that element. The global safety factor of the landslide is defined as the weighted average of all local safety factors based on the area of the slip surface. Some example analyses show that the results computed by the LSF method agree well with those calculated by the General Limit Equilibrium(GLE) method in two-dimensional(2D) models and the distribution of the LSF in the 3D slip zone is consistent with that indicated by the observed deformation pattern of an actual landslide in China.

Limit analysis for the seismic stability of three-dimensional rock slopes using the generalized Hoek-Brown criterion

The parameters that influence slope stability and their criteria of failure are fairly understood but over-conservative design approaches are often preferred,which can result in excessive overburden removal that may jeopardize profitability in the context of open pit mining.Numerical methods such as finite element and discrete element modelling are instrumental to identify specific zones of stability,but they remain approximate and do not pinpoint the critical factors that influence stability without extensive parametric studies.A large number of degrees of freedom and input parameters may make the outcome of numerical modelling insufficient compared to analytical solutions.Existing analytical approaches have not tackled the stability of slopes using non-linear plasticity criteria and threedimensional failure mechanisms.This paper bridges this gap by using the yield design theory and the Hoek-Brown criterion.Moreover,the proposed model includes the effect of seismic forces,which are not always taken into account in slope stability analyses.The results are presented in the form of rigorous mathematical expressions and stability charts involving the loading conditions and the rock mass properties emanating from the plasticity criterion.

Probabilistic seismic stability of three-dimensional slopes by pseudo-dynamic approach

Probabilistic analysis is a rational approach for engineering design because it provides more insight than traditional deterministic analysis. Probabilistic evaluation on seismic stability of three dimensional (3D) slopes is studied in this paper. The slope safety factor is computed by combining the kinematic approach of limit analysis using a three-dimensional rotational failure mechanism with the pseudo-dynamic approach. The variability of input parameters, including six pseudo-dynamic parameters and two soil shear strength parameters, are taken into account by means of Monte-Carlo Simulations (MCS) method. The influences of pseudo-dynamic input variables on the computed failure probabilities are investigated and discussed. It is shown that the obtained failure probabilities increase with the pseudo-dynamic input variables and the pseudo-dynamic approach gives more conservative failure probability estimates compared with the pseudo-static approach.

Criteria for the Nonlinear Stability of Three-Dimensional Quasi-Geostrophic Motions

Some more proper criteria for the nonlinear stability of three-dimensional quasi-geostrophic motions are given by combining variational principle with a prior estimates method. The criteria are suitable for perturbations of initial condition as well as parameters in the model. The basic flow can be steady or unsteady. Particularly the difficulty due to the nonlinear boundary condition is completely overcome by the use of our method.

Undrained Stability Analysis of Three-Dimensional Rectangular Trapdoor in Clay

A new collapse model of the trapdoors,three-dimensional rectangular trapdoor(3DRT),is presented for ground surface collapse.Undrained stability of 3DRT is examined with the upper bound method of plasticity limit analysis theory.The soil where the trapdoors are located is assumed to be a perfectly plastic model with a Tresca yield criterion.Block analysis technique is employed to investigate the collapse of 3DRT.The model is divided into five different block types and added up to ten rigid blocks.According to the law of conservation of energy,the critical stability ratios of 3DRT are obtained through a search proceeding.The results of upper bound solution for 3DRT are given,and three trapdoor models with depth various are discussed during the application in the stability analysis of square trapdoors.The critical stability ratios can be used in the design of underground excavation and support force.

Three-dimensional critical slip surface locating and slope stability assessment for lava lobe of Unzen volcano

Even Unzen volcano has been declared to be in a state of relative dormancy,the latest formed lava lobe No.11 now represents a potential slope failure mass based on the latest research.This paper concentrates on the stability of the lava lobe No.11 and its possible critical sliding mass.It proposes geographic information systems (GIS) based three-dimensional (3D) slope stability analysis models.It uses a 3D locating approach to identify the 3D critical slip surface and to analyze the 3D stability of the lava lobe No.11.At the same time,the new 3D approach shows the effectiveness in selecting the range of the Monte Carlo random variables and locating the critical slip surface in different parts of the lava lobe No.11.The results are very valuable for judging the stability of the lava lobe and assigning the monitoring equipments.

STABILITY OF 3-D DISTURBANCE WAVES ON BOUNDARY LAYER WITH DISTRIBUTED LOCAL ROUGHNESS

The numerical solution of the stable basic flow on a 3-D boundary layer is obtained by using local ejection, local suction, and combination of local ejection and suction to simulate the local rough wall. The evolution of 3-D disturbance T-S wave is studied in the spatial processes, and the effects of form and distribution structure of local roughness on the growth rate of the 3-D disturbance wave and the flow stability are discussed. Numerical results show that the growth of the disturbance wave and the form of vortices are accelerated by the 3-D local roughness. The modification of basic flow owing to the evolvement of the finite amplitude disturbance wave and the existence of spanwise velocity induced by the 3-D local roughness affects the stability of boundary layer. Propagation direction and phase of the disturbance wave shift obviously for the 3-D local roughness of the wall. The flow stability characteristics change if the form of the 2-D local roughness varies.

STABILITY AND LOCAL BIFURCATION IN A SIMPLY-SUPPORTED BEAM CARRYING A MOVING MASS

The stability and local bifurcation of a simply-supported flexible beam （Bernoulli- Euler type） carrying a moving mass and subjected to harmonic axial excitation are investigated. In the theoretical analysis, the partial differential equation of motion with the fifth-order nonlinear term is solved using the method of multiple scales （a perturbation technique）. The stability and local bifurcation of the beam are analyzed for 1/2 sub harmonic resonance. The results show that some of the parameters, especially the velocity of moving mass and external excitation, affect the local bifurcation significantly. Therefore, these parameters play important roles in the system stability.

Three-dimensional numerical analysis of plant-soil hydraulic interactions on pore water pressure of vegetated slope under different rainfall patterns

Understanding the pore water pressure distribution in unsaturated soil is crucial in predicting shallow landslides triggered by rainfall,mainly when dealing with different temporal patterns of rainfall intensity.However,the hydrological response of vegetated slopes,especially three-dimensional(3D)slopes covered with shrubs,under different rainfall patterns remains unclear and requires further investigation.To address this issue,this study adopts a novel 3D numerical model for simulating hydraulic interactions between the root system of the shrub and the surrounding soil.Three series of numerical parametric studies are conducted to investigate the influences of slope inclination,rainfall pattern and rainfall duration.Four rainfall patterns(advanced,bimodal,delayed,and uniform)and two rainfall durations(4-h intense and 168-h mild rainfall)are considered to study the hydrological response of the slope.The computed results show that 17%higher transpiration-induced suction is found for a steeper slope,which remains even after a short,intense rainfall with a 100-year return period.The extreme rainfalls with advanced(PA),bimodal(PB)and uniform(PU)rainfall patterns need to be considered for the short rainfall duration(4 h),while the delayed(PD)and uniform(PU)rainfall patterns are highly recommended for long rainfall durations(168 h).The presence of plants can improve slope stability markedly under extreme rainfall with a short duration(4 h).For the long duration(168 h),the benefit of the plant in preserving pore-water pressure(PWP)and slope stability may not be sufficient.

Global stability coefficient of large underground caverns under static loading and earthquake wave condition

Underground energy and resource development,deep underground energy storage and other projects involve the global stability of multiple interconnected cavern groups under internal and external dynamic disturbances.An evaluation method of the global stability coefficient of underground caverns based on static overload and dynamic overload was proposed.Firstly,the global failure criterion for caverns was defined based on its band connection of plastic-strain between multi-caverns.Then,overloading calculation of the boundary geostress and seismic intensity on the caverns model was carried out,and the critical unstable state of multi-caverns can be identified,if the plastic-strain band appeared between caverns during these overloading processes.Thus,the global stability coefficient for the multi-caverns under static loading and earthquake was obtained based on the corresponding overloading coefficient.Practical analysis for the Yingliangbao(YLB)hydraulic caverns indicated that this method can not only effectively obtain the global stability coefficient of caverns under static and dynamic earthquake conditions,but also identify the caverns’high-risk zone of local instability through localized plastic strain of surrounding rock.This study can provide some reference for the layout design and seismic optimization of underground cavern group.

Study of gas shielding stability for underwater local dry welding based on FLUENT

The research of gas shielding stability technology, playing a key role in underwater local dry welding, was introduced in this paper. The study includes shielding cup design, gas flow simulation, draining and welding experiments. The commercial computational fluid dynamics （ CFD ） software FLUENT was applied to get and compare the speed and pressure contour images of gas in the three kinds of cup with different intake mode. The computed results show that the cup with stilling chamber on the top has the best shielding performance. The underwater welding experiment of 15 meters proves that the shielding cup with stilling chamber can offer a good dry space, the welding arc can burn stably in it and the weld quality is perfect. The research will facilitate the application of underwater local dry welding technology in the maintenance of nuclear power station widely.

STABILIZATION EFFECT OF FRICTIONS FOR TRANSONIC SHOCKS IN STEADY COMPRESSIBLE EULER FLOWS PASSING THREE-DIMENSIONAL DUCTS

Transonic shocks play a pivotal role in designation of supersonic inlets and ramjets.For the three-dimensional steady non-isentropic compressible Euler system with frictions,we constructe a family of transonic shock solutions in rectilinear ducts with square cross-sections.In this article,we are devoted to proving rigorously that a large class of these transonic shock solutions are stable,under multidimensional small perturbations of the upcoming supersonic flows and back pressures at the exits of ducts in suitable function spaces.This manifests that frictions have a stabilization effect on transonic shocks in ducts,in consideration of previous works which shown that transonic shocks in purely steady Euler flows are not stable in such ducts.Except its implications to applications,because frictions lead to a stronger coupling between the elliptic and hyperbolic parts of the three-dimensional steady subsonic Euler system,we develop the framework established in previous works to study more complex and interesting Venttsel problems of nonlocal elliptic equations.

Three-Dimensional Cooperative Localization via Space-Air-Ground Integrated Networks

The space-air-ground integrated network(SAGIN)combines the superiority of the satellite,aerial,and ground communications,which is envisioned to provide high-precision positioning ability as well as seamless connectivity in the 5G and Beyond 5G(B5G)systems.In this paper,we propose a three-dimensional SAGIN localization scheme for ground agents utilizing multi-source information from satellites,base stations and unmanned aerial vehicles(UAVs).Based on the designed scheme,we derive the positioning performance bound and establish a distributed maximum likelihood algorithm to jointly estimate the positions and clock offsets of ground agents.Simulation results demonstrate the validity of the SAGIN localization scheme and reveal the effects of the number of satellites,the number of base stations,the number of UAVs and clock noise on positioning performance.

High-precision three-dimensional Rydberg atom localization in a four-level atomic system

Rydberg atoms have been widely investigated due to their large size,long radiative lifetime,huge polarizability and strong dipole-dipole interactions.The position information of Rydberg atoms provides more possibilities for quantum optics research,which can be obtained under the localization method.We study the behavior of three-dimensional(3D)Rydberg atom localization in a four-level configuration with the measurement of the spatial optical absorption.The atomic localization precision depends strongly on the detuning and Rabi frequency of the involved laser fields.A 100%probability of finding the Rydberg atom at a specific 3D position is achieved with precision of~0.031λ.This work demonstrates the possibility for achieving the 3D atom localization of the Rydberg atom in the experiment.

Stability analysis of slope in strain-softening soils using local arc-length solution scheme

Soils with strain-softening behavior — manifesting as a reduction of strength with increasing plastic strain — are commonly found in the natural environment. For slopes in these soils,a progressive failure mechanism can occur due to a reduction of strength with increasing strain. Finite element method based numerical approaches have been widely performed for simulating such failure mechanism,owning to their ability for tracing the formation and development of the localized shear strain. However,the reliability of the currently used approaches are often affected by poor convergence or significant mesh-dependency,and their applicability is limited by the use of complicated soil models. This paper aims to overcome these limitations by developing a finite element approach using a local arc-length controlled iterative algorithm as the solution strategy. In the proposed finite element approach,the soils are simulated with an elastoplastic constitutive model in conjunction with the Mohr-Coulomb yield function. The strain-softening behavior is represented by a piece-wise linearrelationship between the Mohr-Coulomb strength parameters and the deviatoric plastic strain. To assess the reliability of the proposed finite element approach,comparisons of the numerical solutions obtained by different finite element methods and meshes with various qualities are presented. Moreover,a landslide triggered by excavation in a real expressway construction project is analyzed by the presented finite element approach to demonstrate its applicability for practical engineering problems.

Preventing Pressure Oscillations Does Not Fix Local Linear Stability Issues of Entropy-Based Split-Form High-Order Schemes

Recently,it was discovered that the entropy-conserving/dissipative high-order split-form discontinuous Galerkin discretizations have robustness issues when trying to solve the sim-ple density wave propagation example for the compressible Euler equations.The issue is related to missing local linear stability,i.e.,the stability of the discretization towards per-turbations added to a stable base flow.This is strongly related to an anti-diffusion mech-anism,that is inherent in entropy-conserving two-point fluxes,which are a key ingredi-ent for the high-order discontinuous Galerkin extension.In this paper,we investigate if pressure equilibrium preservation is a remedy to these recently found local linear stability issues of entropy-conservative/dissipative high-order split-form discontinuous Galerkin methods for the compressible Euler equations.Pressure equilibrium preservation describes the property of a discretization to keep pressure and velocity constant for pure density wave propagation.We present the full theoretical derivation,analysis,and show corresponding numerical results to underline our findings.In addition,we characterize numerical fluxes for the Euler equations that are entropy-conservative,kinetic-energy-preserving,pressure-equilibrium-preserving,and have a density flux that does not depend on the pressure.The source code to reproduce all numerical experiments presented in this article is available online(https://doi.org/10.5281/zenodo.4054366).

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

Coexistence and local Mittag–Leffler stability of fractional-order recurrent neural networks with discontinuous activation functions

In this paper, coexistence and local Mittag–Leffler stability of fractional-order recurrent neural networks with discontinuous activation functions are addressed. Because of the discontinuity of the activation function, Filippov solution of the neural network is defined. Based on Brouwer's fixed point theorem and definition of Mittag–Leffler stability, sufficient criteria are established to ensure the existence of (2k + 3)~n (k ≥ 1) equilibrium points, among which (k + 2)~n equilibrium points are locally Mittag–Leffler stable. Compared with the existing results, the derived results cover local Mittag–Leffler stability of both fractional-order and integral-order recurrent neural networks. Meanwhile discontinuous networks might have higher storage capacity than the continuous ones. Two numerical examples are elaborated to substantiate the effective of the theoretical results.

LOCAL STABILITY OF TRAVELLING FRONTS FOR A DAMPED WAVE EQUATION

The paper is concerned with the long-time behaviour of the travelling fronts of the damped wave equation αutt ＋ut = uxx -V′（u） on R. The long-time asymptotics of the solutions of this equation are quite similar to those of the corresponding reaction-diffusion equation ut = uxx - V′（u）. Whereas a lot is known about the local stability of travelling fronts in parabolic systems, for the hyperbolic equations it is only briefly discussed when the potential V is of bistable type. However, for the combustion or monostable type of V, the problem is much more complicated. In this paper, a local stability result for travelling fronts of this equation with combustion type of nonlinearity is established. And then, the result is extended to the damped wave equation with a case of monostable pushed front.