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 coordinates test method with uncertain projectile proximity explosion position based on dynamic seven photoelectric detection screen

To objectively obtain the three-dimensional coordinates of the projectile fuze proximity explosion when projectile intersects the head of missile target, we propose a dynamic seven photoelectric detection screen test method, which is made up of six plane detection screens and a flash photoelectric dynamic detection screen. The three-dimensional coordinates calculation model of the projectile proximity explosion position based on seven plane detection screens with dynamic characteristics is established.According to the relation of the dynamic seven photoelectric detection screen planes and the time values,the analytical function of the projectile proximity explosion position parameters under non-linear motion is derived. The projectile signal filtering method based on discrete wavelet transform is explored in this work. Additionally, the projectile signal recognition algorithm using an improved particle swarm is proposed. Based on the characteristics of the time duration and the signal peak error for the projectile passing through the detection screen, the signals attribution of the same projectile passing through six detection screens are analyzed for obtaining precise time values of the same projectile passing through the detection screens. On the basis of the projectile fuze proximity explosion test, the linear motion model and the proposed non-linear motion model are used to calculate and compare the same group of projectiles proximity explosion position parameters. The comparison of test results verifies that the proposed test method and calculation model in this work accurately obtain the actual projectile proximity explosion position parameters.

Research on the Coupled and Coordinated Development of Tourism and Transport in Chengdu City

Based on a total of 16 indicators selected from the tourism and transport industries,an evaluation index system of the coupling and coordination development level of tourism and transport is constructed.The entropy value method and the coupling coordination degree model are used to conduct an empirical study on the development level and coupling coordination level of the transport and tourism industries in Chengdu City from 2011 to 2020.The results show that,on the whole,the coupling coordination degree of transport and tourism in Chengdu is poor and has been in a state of mild to moderate dysfunction.The development level of tourism lagged behind the development of transport from 2011 to 2012,and the two were in a state of mild dysfunction.However,from 2013 onwards,the development level of tourism was prioritized over the development level of transport.This shift caused the coupling coordination degree of the two industries to decline sharply to 0.23305 in 2013.The development level of the tourism industry increased again,reaching 0.34206 in 2019,which marked an improvement.Consequently,the coordination degree of the transport and tourism industries evolved from moderate dislocation to mild dislocation.Finally,the results of the empirical research are analyzed,and corresponding suggestions are put forward to promote the sustainable growth of the transport and tourism industries in Chengdu City.These suggestions aim to improve the coupled and coordinated development level of the two industries.

Three-dimensional acoustic propagation model for shallow waters based on an indirect boundary element method

This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary element method(BEM),which uses a half-space Green function instead of the more conventional free-space Green function.On the other hand,a numerical implementation is illustrated to explore the formulation in simple idealized cases,controlled by a few parameters,which provides necessary tests for the accuracy and performance of the model.The half-space Green's function,which has been previously used in scattering and diffraction,adds terms to the usual expressions of the integral operators without altering their continuity properties.Verifications against the wavenumber integration solution of the Pekeris waveguide suggest that the model allows an adequate prediction for the acoustic field.Likewise,numerical experiments in relation to the necessary mesh size for the description of the water-marine sediment interface lead to the conclusion that a transmission loss prediction with acceptable accuracy can be obtained with the use of a limited mesh around the desired evaluation region.

Free Vibration Analysis of Rectangular Plate with Cutouts under Elastic Boundary Conditions in Independent Coordinate Coupling Method

Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.

THE THREE-DIMENSIONAL FUNDAMENTAL SOLUTION TO STOKES FLOW IN THE OBLATE SPHEROIDAL COORDINATES WITH APPLICATIONS TO MULTIPLES SPHEROID PROBLEMS

A new three-dimensional fundamental solution to the Stokes flow was proposed by transforming the solid harmonic functions in Lamb's solution into expressions in terms Of the oblate spheroidal coordinates. These fundamental solutions are advantageous in treating flows past an arbitrary number of arbitrarily positioned and oriented oblate spheroids. The least squares technique was adopted herein so that the convergence difficulties often encountered in solving three-dimensional problems were completely avoided. The examples demonstrate that present approach is highly accurate, consistently stable and computationally efficient. The oblate spheroid may be used to model a variety of particle shapes between a circular disk and a sphere. For the first time, the effect of various geometric factors on the forces and torques exerted on two oblate spheroids were systematically studied by using the proposed fundamental solutions. The generality of this approach was illustrated by two problems of three spheroids.

Three-dimensional numerical simulation of the flow past a circular cylinder based on LES method

The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.

Evaluation of the coordinated development of China’s forest resources-economy-environment system

Forests provide enormous ecological, economic, and social benefits, as such forest development should be ori‐ented toward resource-economy-environment harmonization. This paper constructs a comprehensive evalua‐tion index system and uses the coupled coordination degree model to measure the coordinated development level of China’ s forest resources-economy-environment system. The results show that, across 2006‒2020, the level of coupled coordinated development of China’ s forest resources-economy-environment composite system fluctuates in an upward trend, thus gradually developing from an initial imbalance to a high degree of coordi‐nation;the level of coordinated development of each subsystem of the forest resources, economy, and environ‐ment also shows an upward trend. The factors influencing the coordinated development of the forest resource economy-environment system are, in order, the government’ s financial capacity, market environment, scien‐tific and technological innovation capacity, level of economic development, and strength of policy implemen‐tation. Therefore, this paper proposes some measures to improve the coordinated development.

Velocity-Free MS/AE Source Location Method for Three-Dimensional Hole-Containing Structures

Microseismic/acoustic emission(MS/AE)source localization method is crucial for predicting and controlling of potentially dangerous sources of complex structures.However,the locating errors induced by both the irregular structure and pre-measured velocity are poorly understood in existing methods.To meet the high-accuracy locating requirements in complex three-dimensional hole-containing structures,a velocity-free MS/AE source location method is developed in this paper.It avoids manual repetitive training by using equidistant grid points to search the path,which introduces A*search algorithm and uses grid points to accommodate complex structures with irregular holes.It also takes advantage of the velocity-free source location method.To verify the validity of the proposed method,lead-breaking tests were performed on a cubic concrete test specimen with a size of 10 cm10 cm10 cm.It was cut out into a cylindrical empty space with a size of/6cm10 cm.Based on the arrivals,the classical Geiger method and the proposed method are used to locate lead-breaking sources.Results show that the locating error of the proposed method is 1.20 cm,which is less than 2.02 cm of the Geiger method.Hence,the proposed method can effectively locate sources in the complex three-dimensional structure with holes and achieve higher precision requirements.

Three-dimensional simulation method of multipactor in microwave components for high-power space application

Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.

The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate

In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.

Microstructure modeling and virtual test of asphalt mixture based on three-dimensional discrete element method

The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.

Micromechanical modeling of asphalt concrete fracture using a user-defined three-dimensional discrete element method

A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.

Predicting the Dynamic Behavior of Asphalt Concrete Using Three-dimensional Discrete Element Method

A user-defined three-dimensional （3D） discrete element model was presented to predict the dynamic modulus and phase angle of asphalt concrete （AC）. The 3D discrete element method （DEM） model of AC was constructed employing a user-defined computer program developed using the ＂Fish＂ language in PFC3D. Important microstructural features of AC were modeled, including aggregate gradation, air voids and mastic. The irregular shape of aggregate particle was modeled using a clump of spheres. The developed model was validated through comparing with experimental measurements and then used to simulate the cyclic uniaxial compression test, based on which the dynamic modulus and phase angle were calculated from the output stress- strain relationship. The effects of air void content, aggregate stiffness and volumetric fraction on AC modulus were further investigated. The experimental results show that the 3D DEM model is able to accurately predict both dynamic modulus and phase angle of AC across a range of temperature and loading frequencies. The user- defined 3D model also demonstrated significant improvement over the general existing two-dimensional models.

Stream Surface Strip Element Method and Simulation of Three-Dimensional Deformation of Continuous Hot Rolled Strip

A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.

An RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in a Lagrangian coordinate

In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.

Three-dimensional frequency-domain full waveform inversion based on the nearly-analytic discrete method

The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.

Three-Dimensional Analytical Modeling of Axial-Flux Permanent Magnet Drivers

In this paper, the axial-flux permanent magnet driver is modeledand analyzed in a simple and novel way under three-dimensional cylindricalcoordinates. The inherent three-dimensional characteristics of the deviceare comprehensively considered, and the governing equations are solved bysimplifying the boundary conditions. The axial magnetization of the sectorshapedpermanent magnets is accurately described in an algebraic form bythe parameters, which makes the physical meaning more explicit than thepurely mathematical expression in general series forms. The parameters of theBessel function are determined simply and the magnetic field distribution ofpermanent magnets and the air-gap is solved. Furthermore, the field solutionsare completely analytical, which provides convenience and satisfactoryaccuracy for modeling a series of electromagnetic performance parameters,such as the axial electromagnetic force density, axial electromagnetic force,and electromagnetic torque. The correctness and accuracy of the analyticalmodels are fully verified by three-dimensional finite element simulations and a15 kW prototype and the results of calculations, simulations, and experimentsunder three methods are highly consistent. The influence of several designparameters on magnetic field distribution and performance is studied and discussed.The results indicate that the modeling method proposed in this papercan calculate the magnetic field distribution and performance accurately andrapidly, which affords an important reference for the design and optimizationof axial-flux permanent magnet drivers.

Three-dimensional forward modeling for magnetotelluric sounding by finite element method

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

ANALYSIS OF THREE-DIMENSIONAL UPSETTING PROCESS BY THE RIGID-PLASTIC REPRODUCING KERNEL PARTICLE METHOD

A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.