A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path...A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path-planning process into two layers is proposed by designing the velocityobstacle strategy for satisfying timeliness and effectiveness.The upper-level layer focuses on creating an efficient Dubins initial path considering the dynamic constraints of the fixed wing.Subsequently,the lower-level layer detects potential collisions and adjusts its flight paths to avoid collisions by using the threedimensional velocity obstacle method,which describes the maneuvering space of collision avoidance as the intersection space of half space.To further handle the dynamic and collisionavoidance constraints,a priority mechanism is designed to ensure that the adjusted path is still feasible for fixed-wing UAVs.Simulation experiments demonstrate the effectiveness of the proposed method.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
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 ...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.展开更多
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 irregula...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.展开更多
The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation met...The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation methods, a new approach based on the shadow area percentage （SAP） that can be used to quantify surface roughness is proposed in this article. Firstly, by the help of laser scanning technique, the three-dimensional model of the surface of rock discontinuity was established. Secondly, a light source was simulated, and there would be some shadows produced on the model surface. Thirdly, to obtain the value of SAP of each specimen, the shadow detection technique was introduced for use. Fourthly, compared with the result from direct shear testing and based on statistics, an empirical for- mula was found among SAP, normal stress, and shear strength. Data of Yujian （~ River were used as an example, and the following conclusions have been made. （1） In the case of equal normal stress, the peak shear stress is positively proportional to the SAP. （2） The formula for estimating was derived, and the predictions of peak-shear strength made with this equation well agreed with the experimental re- suits obtained in laboratory tests.展开更多
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 met...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 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 elem...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.展开更多
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 su...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.展开更多
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 tim...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.展开更多
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 forwar...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.展开更多
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 t...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.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
The singular hybrid boundary node method （SHBNM） is proposed for solving three-dimensional problems in linear elasticity. The SHBNM represents a coupling between the hybrid displacement variational formulations and ...The singular hybrid boundary node method （SHBNM） is proposed for solving three-dimensional problems in linear elasticity. The SHBNM represents a coupling between the hybrid displacement variational formulations and moving least squares （MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the later. The rigid movement method was employed to solve the hyper-singular integrations. The ＇boundary layer effect＇, which is the main drawback of the original Hybrid BNM, was overcome by an adaptive integration scheme. The source points of the fundamental solution were arranged directly on the boundary. Thus the uncertain scale factor taken in the regular hybrid boundary node method （RHBNM） can be avoided. Numerical examples for some 3D elastic problems were given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method were studied through the numerical examples.展开更多
A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the ...A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.展开更多
This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite differenc...This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.展开更多
In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through min...In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.展开更多
In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination...In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.展开更多
In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method wit...In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method with the k-ε turbulence mathematical model. Then, the water hydrodynamic characteristics and transport processes for BOD_5, NH_(3^-)N and TP were analyzed. The results showed that the water surface of convex bank was a little lower than that of concave bank due to the centrifugal force near the bend, and most concentrations were inferior to the type Ⅴ standard indexes of surface water environmental quality. The model validation indicated that the errors between the simulated and monitored values were comparatively small, satisfying the application demands and providing scientific basis and decision support for the restoration and protection of water quality.展开更多
There are relatively few studies on large rotation or deformation by means of the three-dimensional(3D)numerical manifold method(NMM).A new modified symmetric and antisymmetric decomposition(MSAD)theory is developed a...There are relatively few studies on large rotation or deformation by means of the three-dimensional(3D)numerical manifold method(NMM).A new modified symmetric and antisymmetric decomposition(MSAD)theory is developed and implemented into the 3D NMM,eliminating the false-volume expansion and false-rotation strain/stress problems.The Jaumann rate is used to measure the material rotation,and the geometric stiffness built on the Jaumann rate is deduced.The incremental formulas of the MSAD-based 3D NMM and a practical guide on the implementation of the MSAD theory are given in detail and exemplified.The new theory and formulas can be applied to analyze both large rotation and large deformation problems.Based on the hypoelasto-plasticity theory and the unified strength theory,the unified yield criterion with associated flow rule is implemented into the MSAD-based 3D NMM.Several typical examples are studied,showing the advantage and potential of the new MSAD theory and the MSAD-based 3D NMM.展开更多
A appropriate size with three-dimension(3 D) channels for lithium diffusion plays an important role in constructing highperforming LiNi_(0.5)Mn_(1.5)O_4(LNMO) cathode materials, as it can not only reduce the transport...A appropriate size with three-dimension(3 D) channels for lithium diffusion plays an important role in constructing highperforming LiNi_(0.5)Mn_(1.5)O_4(LNMO) cathode materials, as it can not only reduce the transport path of lithium ions and electrons, but also reduce the side effects and withstand the structural strain in the process of repetitive Li~+ intercalation/deintercalation. In this work, an e fficient method for designing the hollow LNMO microsphere with 3 D channels structure by using polyethylene oxide(PEO) as soft template agent assisted solvothermal method is proposed. Experimental results indicate that PEO can make the reagents mingle evenly and nucleate slowly in the solvothermal process, thus obtaining a homogeneous distribution of carbonate precursors. In the final LNMO products, the hollow 3 D channels structure obtained by the decomposition of PEO and carbonate precursor in the calcination can provide abundant electroactive zones and electron/ion transport paths during the charge/discharge process, which benefits to improve the cycling performance and rate capability. The LNMO prepared by adding 1 g PEO possesses the most outstanding electrochemical performance, which presented an excellent discharge capacity of 143.1 mAh g~(-1) at 0.1 C and with a capacity retention of 92.2% after 100 cycles at 1 C. The superior performance attributed to the 3 D channels structure of hollow microspheres, which provide uninterrupted conductive systems and therefore achieve the stable transfer for electron/ion.展开更多
基金supported by the National Science Fund for Distinguished Young Scholars(52425211)BIT Research Fund Program for Young Scholars(XSQD-202201005).
文摘A three-dimensional path-planning approach has been developed to coordinate multiple fixed-wing unmanned aerial vehicles(UAVs)while avoiding collisions.The hierarchical path-planning architecture that divides the path-planning process into two layers is proposed by designing the velocityobstacle strategy for satisfying timeliness and effectiveness.The upper-level layer focuses on creating an efficient Dubins initial path considering the dynamic constraints of the fixed wing.Subsequently,the lower-level layer detects potential collisions and adjusts its flight paths to avoid collisions by using the threedimensional velocity obstacle method,which describes the maneuvering space of collision avoidance as the intersection space of half space.To further handle the dynamic and collisionavoidance constraints,a priority mechanism is designed to ensure that the adjusted path is still feasible for fixed-wing UAVs.Simulation experiments demonstrate the effectiveness of the proposed method.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
文摘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.
基金The authors wish to acknowledge financial support from the National Natural Science Foundation of China(51822407 and 51774327)Natural Science Foundation of Hunan Province in China(2018JJ1037)Innovation Driven project of Central South University(2020CX014).
文摘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.
基金supported by the China Geological Survey (No.1212011014030)the Major State Basic Research Development Program of China (973 Program) (No.2011CB710600)
文摘The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation methods, a new approach based on the shadow area percentage （SAP） that can be used to quantify surface roughness is proposed in this article. Firstly, by the help of laser scanning technique, the three-dimensional model of the surface of rock discontinuity was established. Secondly, a light source was simulated, and there would be some shadows produced on the model surface. Thirdly, to obtain the value of SAP of each specimen, the shadow detection technique was introduced for use. Fourthly, compared with the result from direct shear testing and based on statistics, an empirical for- mula was found among SAP, normal stress, and shear strength. Data of Yujian （~ River were used as an example, and the following conclusions have been made. （1） In the case of equal normal stress, the peak shear stress is positively proportional to the SAP. （2） The formula for estimating was derived, and the predictions of peak-shear strength made with this equation well agreed with the experimental re- suits obtained in laboratory tests.
基金Project supported by the National Key Laboratory Foundation,China(Grant No.9140C530103110C5301)
文摘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.
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘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.
基金Sponsored by National Natural Science Foundation of China(50175095)Provincial Natural Science Foundation of Hebei of China(502173)
文摘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.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘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.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘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.
基金This work was supported by the National Natural Science Foundation of China （No. 50275094）.
文摘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.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
基金Project supported by the Program of the Key Laboratory of Rock and Soil Mechanics of Chinese Academy of Sciences (No.Z110507)
文摘The singular hybrid boundary node method （SHBNM） is proposed for solving three-dimensional problems in linear elasticity. The SHBNM represents a coupling between the hybrid displacement variational formulations and moving least squares （MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the later. The rigid movement method was employed to solve the hyper-singular integrations. The ＇boundary layer effect＇, which is the main drawback of the original Hybrid BNM, was overcome by an adaptive integration scheme. The source points of the fundamental solution were arranged directly on the boundary. Thus the uncertain scale factor taken in the regular hybrid boundary node method （RHBNM） can be avoided. Numerical examples for some 3D elastic problems were given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method were studied through the numerical examples.
文摘A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.
文摘This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.
基金supported by the National Key Research and Development Program of China(Grant No.2019YFC1510400)the National Natural Science Foundation of China(Grant Nos.41975054 and 41930967)the Special Fund for Forecasters of China Meteorological Administration(Grant No.CMAYBY2018-040)。
文摘In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.
基金Project(2017YFC0404802)supported by the National Key R&D Program of ChinaProjects(U1965206,51979143)supported by the National Natural Science Foundation of China。
文摘In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.
基金Supported by the Innovative Research Groups of National Natural Science Foundation of China(No.51321065)the Major Science and Technology Program for Water Pollution Control and Treatment(2012ZX07101-008)the National Natural Science Foundation of China(No.51439005)
文摘In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method with the k-ε turbulence mathematical model. Then, the water hydrodynamic characteristics and transport processes for BOD_5, NH_(3^-)N and TP were analyzed. The results showed that the water surface of convex bank was a little lower than that of concave bank due to the centrifugal force near the bend, and most concentrations were inferior to the type Ⅴ standard indexes of surface water environmental quality. The model validation indicated that the errors between the simulated and monitored values were comparatively small, satisfying the application demands and providing scientific basis and decision support for the restoration and protection of water quality.
基金This research is supported by the National Basic Research Program of China(973 Program,Grant No.2014CB047100)the National Natural Science Foundation of China(Grant Nos.41472289,51179185 and 41807275).
文摘There are relatively few studies on large rotation or deformation by means of the three-dimensional(3D)numerical manifold method(NMM).A new modified symmetric and antisymmetric decomposition(MSAD)theory is developed and implemented into the 3D NMM,eliminating the false-volume expansion and false-rotation strain/stress problems.The Jaumann rate is used to measure the material rotation,and the geometric stiffness built on the Jaumann rate is deduced.The incremental formulas of the MSAD-based 3D NMM and a practical guide on the implementation of the MSAD theory are given in detail and exemplified.The new theory and formulas can be applied to analyze both large rotation and large deformation problems.Based on the hypoelasto-plasticity theory and the unified strength theory,the unified yield criterion with associated flow rule is implemented into the MSAD-based 3D NMM.Several typical examples are studied,showing the advantage and potential of the new MSAD theory and the MSAD-based 3D NMM.
基金funded by the National Natural Science Foundation of China(No.21776051)the Natural Science Foundations of Guangdong(No.2018A030313423)。
文摘A appropriate size with three-dimension(3 D) channels for lithium diffusion plays an important role in constructing highperforming LiNi_(0.5)Mn_(1.5)O_4(LNMO) cathode materials, as it can not only reduce the transport path of lithium ions and electrons, but also reduce the side effects and withstand the structural strain in the process of repetitive Li~+ intercalation/deintercalation. In this work, an e fficient method for designing the hollow LNMO microsphere with 3 D channels structure by using polyethylene oxide(PEO) as soft template agent assisted solvothermal method is proposed. Experimental results indicate that PEO can make the reagents mingle evenly and nucleate slowly in the solvothermal process, thus obtaining a homogeneous distribution of carbonate precursors. In the final LNMO products, the hollow 3 D channels structure obtained by the decomposition of PEO and carbonate precursor in the calcination can provide abundant electroactive zones and electron/ion transport paths during the charge/discharge process, which benefits to improve the cycling performance and rate capability. The LNMO prepared by adding 1 g PEO possesses the most outstanding electrochemical performance, which presented an excellent discharge capacity of 143.1 mAh g~(-1) at 0.1 C and with a capacity retention of 92.2% after 100 cycles at 1 C. The superior performance attributed to the 3 D channels structure of hollow microspheres, which provide uninterrupted conductive systems and therefore achieve the stable transfer for electron/ion.
