We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on ex...We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for two non-canonical Hamiltonian systems. Numerical tests show that the proposed methods exhibit good numerical performance in preserving the phase orbit and the energy of the system over long time, whereas higher order Runge–Kutta methods do not preserve these properties. Numerical tests also show that the K-symplectic methods exhibit better efficiency than that of the same order implicit symplectic, explicit and implicit symplectic methods for the original nonseparable non-canonical systems. On the other hand, the fourth order K-symplectic method is more efficient than the fourth order Yoshida’s method, the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om explicit K-symplectic methods for the extended phase space Hamiltonians, but less efficient than the the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om extended phase space symplectic-like methods with the midpoint permutation.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this...Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this paper, for the sake of illustrating in detail how dynamic explicit finite element method is applied to the numerical simulation of the autobody panel forming process,an example of optimization of stamping process pain meters of an inner door panel is presented. Using dynamic explicit finite element code Ls-DYNA3D, the inner door panel has been optimized by adapting pa- rameters such as the initial blank geometry and position, blank-holder forces and the location of drawbeads, and satisfied results are obtained.展开更多
A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a...A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow pheno...Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.展开更多
It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration method's concom...It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration method's concomitant explicit method was used as an example and is described in detail in this paper.The relationship between the implicit method and explicit method is defined,which provides some guidance about how to create a new explicit method that has high precision and computational efficiency.展开更多
We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusi...We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.展开更多
In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on ...In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.展开更多
An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are as...An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.展开更多
The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation too...The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation tool for analyzing complex three-dimensional sheet metal forming problems. The theory and features of the dynamic explicit finite element methods are introduced and the available various commercial finite element method codes used for sheet metal forming simulation in the world are discussed,and the civil and international status quo of automobile panel simulation as well. The front door outer panel of one certain new automobile is regarded as one example that the dynamic explicit FEM code Dynaform is used for the simulation of the front door outer panel forming process. Process defects such as ruptures are predicted. The improving methods can be given according to the simulation results. Foreground of sheet metal forming simulation is outlined.展开更多
In this paper,we introduce a new class of explicit numerical methods called the tamed stochastic Runge-Kutta-Chebyshev(t-SRKC)methods,which apply the idea of taming to the stochastic Runge-Kutta-Chebyshev(SRKC)methods...In this paper,we introduce a new class of explicit numerical methods called the tamed stochastic Runge-Kutta-Chebyshev(t-SRKC)methods,which apply the idea of taming to the stochastic Runge-Kutta-Chebyshev(SRKC)methods.The key advantage of our explicit methods is that they can be suitable for stochastic differential equations with non-globally Lipschitz coefficients and stiffness.Under certain non-globally Lipschitz conditions,we study the strong convergence of our methods and prove that the order of strong convergence is 1/2.To show the advantages of our methods,we compare them with some existing explicit methods(including the Euler-Maruyama method,balanced Euler-Maruyama method and two types of SRKC methods)through several numerical examples.The numerical results show that our t-SRKC methods are efficient,especially for stiff stochastic differential equations.展开更多
We study a temporal step size control of explicit Runge-Kutta(RK)methods for com-pressible computational fuid dynamics(CFD),including the Navier-Stokes equations and hyperbolic systems of conservation laws such as the...We study a temporal step size control of explicit Runge-Kutta(RK)methods for com-pressible computational fuid dynamics(CFD),including the Navier-Stokes equations and hyperbolic systems of conservation laws such as the Euler equations.We demonstrate that error-based approaches are convenient in a wide range of applications and compare them to more classical step size control based on a Courant-Friedrichs-Lewy(CFL)number.Our numerical examples show that the error-based step size control is easy to use,robust,and efcient,e.g.,for(initial)transient periods,complex geometries,nonlinear shock captur-ing approaches,and schemes that use nonlinear entropy projections.We demonstrate these properties for problems ranging from well-understood academic test cases to industrially relevant large-scale computations with two disjoint code bases,the open source Julia pack-ages Trixi.jl with OrdinaryDiffEq.jl and the C/Fortran code SSDC based on PETSc.展开更多
In this paper,numerical methods for the time-changed stochastic differential equations of the form dY(t)=a(Y(t))dt+b(Y(t))dE(t)+s(Y(t))dB(E(t))are investigated,where all the coefficients a(·),b(·)and s(·...In this paper,numerical methods for the time-changed stochastic differential equations of the form dY(t)=a(Y(t))dt+b(Y(t))dE(t)+s(Y(t))dB(E(t))are investigated,where all the coefficients a(·),b(·)and s(·)are allowed to contain some super-linearly growing terms.An explicit method is proposed by using the idea of truncating terms that grow too fast.Strong convergence in the finite time of the proposed method is proved and the convergence rate is obtained.The proposed method is also proved to be able to reproduce the asymptotic stability of the underlying equation in the almost sure sense.Simulations are provided to demonstrate the theoretical results.展开更多
Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to dra...Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to draw out governing equation of beam. The interactions between the dam and the soil were modeled by using a translational spring and a rotational spring. A Sommerfeld＇s radiation condition at the infinity boundary of the fluid domain was adopted. The effects of the reservoir bottom absorption and surface waves on the dam-reservoir-foundation interaction due to the earthquake were studied. To avoid the instability of solution, a semi-implicit scheme was used for the discretization of the governing equation of dam and an explicit scheme was used for the discretization of the governing equation of fluid. The results show that as the slope of upstream dam increases, the hydrodynamic pressure on the dam is reduced. Moreover, when the Timoshenko beam theory is used, the system response increases.展开更多
The stowing and deploying experiment was conducted for three 700 mm long thin-walled tubes,and the structural behavior characteristics parameters were measured clearly,including strain,deformation and wrapping moment....The stowing and deploying experiment was conducted for three 700 mm long thin-walled tubes,and the structural behavior characteristics parameters were measured clearly,including strain,deformation and wrapping moment.3D finite element models(FEM)were built subsequently and explicit dynamic method was used to simulate the stowing and deploying of the lenticular carbon fiber reinforced polymer(CFRP)thin-walled tubular space boom,which was designed as four-ply(45°/-45°/45°/-45°)lay-up.The stress and energy during the wrapping process were got and compared with different wrapping angular velocity,the reasonable wrapping angular velocity and effective method were conformed,and structural behavior characteristics were obtained.The results were compared and discussed as well,and the results show that the numerical results by 0.628 rad/s velocity agree well with the measured values.In this paper,the numerical procedure and experimental results are valuable to the optimization design of CFRP thin-walled tubular space boom and future research.展开更多
Broken gap is an extremely dangerous state in the service of high-speed rails,and the violent wheel–rail impact forces will be intensified when a vehicle passes the gap at high speeds,which may cause a secondary frac...Broken gap is an extremely dangerous state in the service of high-speed rails,and the violent wheel–rail impact forces will be intensified when a vehicle passes the gap at high speeds,which may cause a secondary fracture to rail and threaten the running safety of the vehicle.To recognize the damage tolerance of rail fracture length,the implicit–explicit sequential approach is adopted to simulate the wheel–rail high-frequency impact,which considers the factors such as the coupling effect between frictional contact and structural vibration,nonlinear material and real geometric profile.The results demonstrate that the plastic deformation and stress are distributed in crescent shape during the impact at the back rail end,increasing with the rail fracture length.The axle box acceleration in the frequency domain displays two characteristic modes with frequencies around 1,637 and 404 Hz.The limit of the rail fracture length is 60 mm for high-speed railway at a speed of 250 km/h.展开更多
By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integ...By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integrate each part by calculating its phase flow,and develop explicit symplectic integrators of different orders for the original Hamiltonian by composing the phase flows.A 2nd-order reversible constructed symplectic scheme is employed to simulate solitons motion and invariants behavior of the NLSE.The simulation results are compared with a 3rd-order non-symplectic implicit Runge-Kutta method,and the convergence of the formal energy of this symplectic integrator is also verified.The numerical results indicate that the explicit symplectic scheme obtained via L-L-N splitting is an effective numerical tool for solving the NLSE.展开更多
In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complic...In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.展开更多
In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum alo...In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 11901564 and 12171466)。
文摘We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for two non-canonical Hamiltonian systems. Numerical tests show that the proposed methods exhibit good numerical performance in preserving the phase orbit and the energy of the system over long time, whereas higher order Runge–Kutta methods do not preserve these properties. Numerical tests also show that the K-symplectic methods exhibit better efficiency than that of the same order implicit symplectic, explicit and implicit symplectic methods for the original nonseparable non-canonical systems. On the other hand, the fourth order K-symplectic method is more efficient than the fourth order Yoshida’s method, the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om explicit K-symplectic methods for the extended phase space Hamiltonians, but less efficient than the the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om extended phase space symplectic-like methods with the midpoint permutation.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
文摘Taking CPU time cost and analysis accuracy into account, dynamic explicit finite ele- ment method is adopted to optimize the forming process of autobody panels that often have large sizes and complex geometry. In this paper, for the sake of illustrating in detail how dynamic explicit finite element method is applied to the numerical simulation of the autobody panel forming process,an example of optimization of stamping process pain meters of an inner door panel is presented. Using dynamic explicit finite element code Ls-DYNA3D, the inner door panel has been optimized by adapting pa- rameters such as the initial blank geometry and position, blank-holder forces and the location of drawbeads, and satisfied results are obtained.
基金Science Council, Chinese Taipei Under Grant No. NSC-95-2221-E-027-099
文摘A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金King Mongkut’s University of Technology North Bangkok (KMUTNB)the Office of the Higher Education Commission (OHEC)the National Metal and Materials Technology Center (MTEC) for supporting this research work
文摘Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.
基金Fundamental Research Funds for the Central Universities
文摘It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration method's concomitant explicit method was used as an example and is described in detail in this paper.The relationship between the implicit method and explicit method is defined,which provides some guidance about how to create a new explicit method that has high precision and computational efficiency.
基金Open Access funding enabled and organized by Projekt DEAL.
文摘We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.
基金National Natural Scienccs Foundation of China (50178005).
文摘In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.
基金Fundamental Research Funds for the Central Universities under Grant No.2682022CX072the Research and Development Plan in Key Areas of Guangdong Province under Grant No.2020B0202010008。
文摘An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.
文摘The paper starts with a brief overview to the necessity of sheet metal forming simulation and the complexity of automobile panel forming, then leads to finite element analysis (FEA) which is a powerful simulation tool for analyzing complex three-dimensional sheet metal forming problems. The theory and features of the dynamic explicit finite element methods are introduced and the available various commercial finite element method codes used for sheet metal forming simulation in the world are discussed,and the civil and international status quo of automobile panel simulation as well. The front door outer panel of one certain new automobile is regarded as one example that the dynamic explicit FEM code Dynaform is used for the simulation of the front door outer panel forming process. Process defects such as ruptures are predicted. The improving methods can be given according to the simulation results. Foreground of sheet metal forming simulation is outlined.
基金supported by the National Natural Science Foundation of China(Grant Nos.12101525,12071403)by the Natural Science Foundation of Hunan Province of China(Grant No.2023JJ40615)+1 种基金by the Research Foundation of Education Department of Hunan Province of China(Grant No.21A0108)by the Research Initiation Fund Project of Xiangtan University(Grant No.21QDZ16).
文摘In this paper,we introduce a new class of explicit numerical methods called the tamed stochastic Runge-Kutta-Chebyshev(t-SRKC)methods,which apply the idea of taming to the stochastic Runge-Kutta-Chebyshev(SRKC)methods.The key advantage of our explicit methods is that they can be suitable for stochastic differential equations with non-globally Lipschitz coefficients and stiffness.Under certain non-globally Lipschitz conditions,we study the strong convergence of our methods and prove that the order of strong convergence is 1/2.To show the advantages of our methods,we compare them with some existing explicit methods(including the Euler-Maruyama method,balanced Euler-Maruyama method and two types of SRKC methods)through several numerical examples.The numerical results show that our t-SRKC methods are efficient,especially for stiff stochastic differential equations.
基金Open Access funding enabled and organized by Projekt DEAL.Andrew Winters was funded through Vetenskapsrådet,Sweden Grant Agreement 2020-03642 VR.Some computations were enabled by resources provided by the Swedish National Infrastructure for Computing(SNIC)at Tetralith,par-tially funded by the Swedish Research Council under Grant Agreement No.2018-05973Hugo Guillermo Castro was funded through the award P2021-0004 of King Abdullah University of Science and Technol-ogy.Some of the simulations were enabled by the Supercomputing Laboratory and the Extreme Comput-ing Research Center at King Abdullah University of Science and Technology.Gregor Gassner acknowl-edges funding through the Klaus-Tschira Stiftung via the project“HiFiLab”.Gregor Gassner and Michael Schlottke-Lakemper acknowledge funding from the Deutsche Forschungsgemeinschaft through the research unit“SNuBIC”(DFG-FOR5409).
文摘We study a temporal step size control of explicit Runge-Kutta(RK)methods for com-pressible computational fuid dynamics(CFD),including the Navier-Stokes equations and hyperbolic systems of conservation laws such as the Euler equations.We demonstrate that error-based approaches are convenient in a wide range of applications and compare them to more classical step size control based on a Courant-Friedrichs-Lewy(CFL)number.Our numerical examples show that the error-based step size control is easy to use,robust,and efcient,e.g.,for(initial)transient periods,complex geometries,nonlinear shock captur-ing approaches,and schemes that use nonlinear entropy projections.We demonstrate these properties for problems ranging from well-understood academic test cases to industrially relevant large-scale computations with two disjoint code bases,the open source Julia pack-ages Trixi.jl with OrdinaryDiffEq.jl and the C/Fortran code SSDC based on PETSc.
基金Wei Liu would like to thank Shanghai Rising-Star Program(Grant No.22QA1406900)Science and Technology Innovation Plan of Shanghai(Grant No.20JC1414200)the National Natural Science Foundation of China(Grant Nos.11871343,11971316 and 12271368)for their financial support.
文摘In this paper,numerical methods for the time-changed stochastic differential equations of the form dY(t)=a(Y(t))dt+b(Y(t))dE(t)+s(Y(t))dB(E(t))are investigated,where all the coefficients a(·),b(·)and s(·)are allowed to contain some super-linearly growing terms.An explicit method is proposed by using the idea of truncating terms that grow too fast.Strong convergence in the finite time of the proposed method is proved and the convergence rate is obtained.The proposed method is also proved to be able to reproduce the asymptotic stability of the underlying equation in the almost sure sense.Simulations are provided to demonstrate the theoretical results.
文摘Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to draw out governing equation of beam. The interactions between the dam and the soil were modeled by using a translational spring and a rotational spring. A Sommerfeld＇s radiation condition at the infinity boundary of the fluid domain was adopted. The effects of the reservoir bottom absorption and surface waves on the dam-reservoir-foundation interaction due to the earthquake were studied. To avoid the instability of solution, a semi-implicit scheme was used for the discretization of the governing equation of dam and an explicit scheme was used for the discretization of the governing equation of fluid. The results show that as the slope of upstream dam increases, the hydrodynamic pressure on the dam is reduced. Moreover, when the Timoshenko beam theory is used, the system response increases.
基金the National Natural Science Foundation of China(No.50878128)the Shanghai Aerospace Foundation(No.HTJ10-15)
文摘The stowing and deploying experiment was conducted for three 700 mm long thin-walled tubes,and the structural behavior characteristics parameters were measured clearly,including strain,deformation and wrapping moment.3D finite element models(FEM)were built subsequently and explicit dynamic method was used to simulate the stowing and deploying of the lenticular carbon fiber reinforced polymer(CFRP)thin-walled tubular space boom,which was designed as four-ply(45°/-45°/45°/-45°)lay-up.The stress and energy during the wrapping process were got and compared with different wrapping angular velocity,the reasonable wrapping angular velocity and effective method were conformed,and structural behavior characteristics were obtained.The results were compared and discussed as well,and the results show that the numerical results by 0.628 rad/s velocity agree well with the measured values.In this paper,the numerical procedure and experimental results are valuable to the optimization design of CFRP thin-walled tubular space boom and future research.
基金The work is supported by the National Natural Science Foundation of China(Nos.51608459,51778542 and U1734207)Fundamental Research Funds for the Central Universities(No.2682018CX01)Cultivation Program for the Excellent Doctoral Dissertation of Southwest Jiaotong University.
文摘Broken gap is an extremely dangerous state in the service of high-speed rails,and the violent wheel–rail impact forces will be intensified when a vehicle passes the gap at high speeds,which may cause a secondary fracture to rail and threaten the running safety of the vehicle.To recognize the damage tolerance of rail fracture length,the implicit–explicit sequential approach is adopted to simulate the wheel–rail high-frequency impact,which considers the factors such as the coupling effect between frictional contact and structural vibration,nonlinear material and real geometric profile.The results demonstrate that the plastic deformation and stress are distributed in crescent shape during the impact at the back rail end,increasing with the rail fracture length.The axle box acceleration in the frequency domain displays two characteristic modes with frequencies around 1,637 and 404 Hz.The limit of the rail fracture length is 60 mm for high-speed railway at a speed of 250 km/h.
基金This research is partially supported by the Informatization Construction of Knowledge Innovation Projects of the Chinese Academy of Sciences“Supercomputing En-vironment Construction and Application”(INF105-SCE)National Natural Science Foundation of China(Grant Nos.10471145 and 10672143).
文摘By performing a particular spatial discretization to the nonlinear Schrodinger equation(NLSE),we obtain a non-integrable Hamiltonian system which can be decomposed into three integrable parts(L-L-N splitting).We integrate each part by calculating its phase flow,and develop explicit symplectic integrators of different orders for the original Hamiltonian by composing the phase flows.A 2nd-order reversible constructed symplectic scheme is employed to simulate solitons motion and invariants behavior of the NLSE.The simulation results are compared with a 3rd-order non-symplectic implicit Runge-Kutta method,and the convergence of the formal energy of this symplectic integrator is also verified.The numerical results indicate that the explicit symplectic scheme obtained via L-L-N splitting is an effective numerical tool for solving the NLSE.
基金National Natural Science Foundation of China Under Grant No. 50178065
文摘In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy.
文摘In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.
