Traditional Chinese medicine (TCM) with its millenniaold wisdom rooted in the principles of holistic Yin-Yang balance and “Bianzheng Lunzhi”[辨证论治, Zhenghou(证候) differentiation and treatment], has long offered ...Traditional Chinese medicine (TCM) with its millenniaold wisdom rooted in the principles of holistic Yin-Yang balance and “Bianzheng Lunzhi”[辨证论治, Zhenghou(证候) differentiation and treatment], has long offered a unique lens to understand human health and disease.However, the modern scientific interpretation of TCM remains at the stage of “knowing that it works, but not knowing why it works”.展开更多
Based on protein-DNA complex crystal structural data in up-to-date Nucleic Acid Database,the related parameters of DNA Kinetic Structure were investigated by Monte-Carlo Multiple Integrals on the base of modified DNA ...Based on protein-DNA complex crystal structural data in up-to-date Nucleic Acid Database,the related parameters of DNA Kinetic Structure were investigated by Monte-Carlo Multiple Integrals on the base of modified DNA structure statistical mechanical model,and time complexity and precision were analyzed on the calculated results.展开更多
By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variable...By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variables. Firstly, the spatial space and temporal domain are discretized by FEM and precise integral algorithm respectively. Then, the high accuracy semi-analytical solution of direct problem can be got. Finally, based on the solution, the computing model of inverse problem and expression of sensitivity analysis are established. Single variable and variables combined identifications including thermal parameters, boundary conditions and source-related terms etc. are given to validate the approach proposed in 1-D and 2-D cases. The effects of noise data and initial guess on the results are investigated. The numerical examples show the effectiveness of this approach.展开更多
A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in t...A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.展开更多
In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DN...In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DNN)method.The NPIM is employed to modify the exponential matrix and loading vector based on the DNN integral method.This involves incorporating the basic assumption of the Newmark-βmethod into the dynamic equation and eliminating the acceleration term from the dynamic equilibrium equation.As a result,the equation is reduced to a first-order linear equation system.Subsequently,the PIM is applied to integrate the system step by step within the NPIM.The DNN method is adopted to solve the inhomogeneous term through fitting the integrand and the original function with a pair of neural networks,and the integral term is solved using the Newton–Leibniz formula.Numerical examples demonstrate that the proposed methodology significantly improves computing efficiency and provides sufficient precision compared to the DNN method.This is particularly evident when analyzing large-scale structures under blast loading conditions.展开更多
The vehicle-road coupling dynamics problem is a prominent issue in transportation,drawing significant attention in recent years.These dynamic equations are characterized by high-dimensionality,coupling,and time-varyin...The vehicle-road coupling dynamics problem is a prominent issue in transportation,drawing significant attention in recent years.These dynamic equations are characterized by high-dimensionality,coupling,and time-varying dynamics,making the exact solutions challenging to obtain.As a result,numerical integration methods are typically employed.However,conventional methods often suffer from low computational efficiency.To address this,this paper explores the application of the parameter freezing precise exponential integrator to vehicle-road coupling models.The model accounts for road roughness irregularities,incorporating all terms unrelated to the linear part into the algorithm's inhomogeneous vector.The general construction process of the algorithm is detailed.The validity of numerical results is verified through approximate analytical solutions(AASs),and the advantages of this method over traditional numerical integration methods are demonstrated.Multiple parameter freezing precise exponential integrator schemes are constructed based on the Runge-Kutta framework,with the fourth-order four-stage scheme identified as the optimal one.The study indicates that this method can quickly and accurately capture the dynamic system's vibration response,offering a new,efficient approach for numerical studies of high-dimensional vehicle-road coupling systems.展开更多
Seismic modeling is a useful tool for studying the propagation of seismic waves within complex structures. However, traditional methods of seismic simulation cannot meet the needs for studying seismic wavefields in th...Seismic modeling is a useful tool for studying the propagation of seismic waves within complex structures. However, traditional methods of seismic simulation cannot meet the needs for studying seismic wavefields in the complex geological structures found in seismic exploration of the mountainous area in Northwestern China. More powerful techniques of seismic modeling are demanded for this purpose. In this paper, two methods of finite element-finite difference method （FE-FDM） and arbitrary difference precise integration （ADPI） for seismic forward modeling have been developed and implemented to understand the behavior of seismic waves in complex geological subsurface structures and reservoirs. Two case studies show that the FE-FDM and ADPI techniques are well suited to modeling seismic wave propagation in complex geology.展开更多
The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration a...D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.展开更多
A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analys...This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analysis so that the parameter N, which has been takenas a constant and an independent pa- rameter without consideration ofthe problems in the original method, can be generated automaticallyby the algorithm itself.展开更多
Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix d...Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.展开更多
Earthquake is a kind of sudden and destructive random excitation in nature.It is significant to determine the probability distribution characteristics of the corresponding dynamic indicators to ensure the safety and t...Earthquake is a kind of sudden and destructive random excitation in nature.It is significant to determine the probability distribution characteristics of the corresponding dynamic indicators to ensure the safety and the stability of structures when the intensive seismic excitation,the intensity of which is larger than 7,acts in train-bridge system.Firstly,the motion equations of a two-dimensional train-bridge system under the vertical random excitation of track irregularity and the vertical seismic acceleration are established,where the train subsystem is composed of 8 mutually independent vehicle elements with 48 degrees of freedom,while the single-span simple supported bridge subsystem is composed of 102D beam elements with 20 degrees of freedom on beam and 2 large mass degrees of freedom at the support.Secondly,Monte Carlo method and pseudo excitation method are adopted to analyze the statistical parameters of the system.The power spectrum density of random excitation is used to define a series of non-stationary pseudo excitation in pseudo excitation method and the trigonometric series of random vibration history samples in Monte Carlo method,respectively solved by precise integral method and Newmark-βmethod through the inter-system iterative procedure.Finally,the results are compared with the case under the weak seismic excitation,and show that the samples of vertical acceleration response of bridge and the offload factor of train obeys the normal distribution.In a high probability,the intensive earthquakes pose a greater threat to the safety and stability of bridges and trains than the weak ones.展开更多
This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using prec...This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.展开更多
The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear...The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.展开更多
The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcom...The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.展开更多
The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was...The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was discussed in this study.To efficiently and accurately compute the exponential of the dynamics state matrix and the matrix functions due to external loads,an adaptively filtered precise integration method was proposed,which inherits the high precision of the precise integrationmethod,improves the computational efficiency and saves the memory required.Moreover,the perturbation method was introduced to avoid repeated computations of matrix exponential and terms due to external loads.Based on the filtering and perturbation techniques,an adaptively filtered precise integration method considering perturbation for stochastic dynamics problems was developed.Two numerical experiments,including a model of phononic crystal and a bridge model considering random parameters,were performed to test the performance of the proposed method in terms of accuracy and efficiency.Numerical results show that the accuracy and efficiency of the proposed method are better than those of the existing precise integration method,the Newmark-βmethod and the Wilson-θmethod.展开更多
The estimation of the disturbance input acting on a vehicle from its given responses is an inverse problem.To overcome some of the issues related to ill-posed inverse problems,this work proposes a method of reconstruc...The estimation of the disturbance input acting on a vehicle from its given responses is an inverse problem.To overcome some of the issues related to ill-posed inverse problems,this work proposes a method of reconstructing the road roughness based on the Kalman filter method.A half-car model that considers both the vehicle and equipment is established,and the joint input-state estimation method is used to identify the road profile.The capabilities of this methodology in the presence of noise are numerically demonstrated.Moreover,to reduce the influence of the driving speed on the estimation results,a method of choosing the calculation frequency is proposed.A road vibration test is conducted to benchmark the proposed method.展开更多
An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary diffe...The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.展开更多
文摘Traditional Chinese medicine (TCM) with its millenniaold wisdom rooted in the principles of holistic Yin-Yang balance and “Bianzheng Lunzhi”[辨证论治, Zhenghou(证候) differentiation and treatment], has long offered a unique lens to understand human health and disease.However, the modern scientific interpretation of TCM remains at the stage of “knowing that it works, but not knowing why it works”.
基金Supported by Inner Mongolia Natural Science Foundation(200711020112)Innovation Fundation of Inner Mongolia University of Science and Technology (2009NC064)~~
文摘Based on protein-DNA complex crystal structural data in up-to-date Nucleic Acid Database,the related parameters of DNA Kinetic Structure were investigated by Monte-Carlo Multiple Integrals on the base of modified DNA structure statistical mechanical model,and time complexity and precision were analyzed on the calculated results.
文摘By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variables. Firstly, the spatial space and temporal domain are discretized by FEM and precise integral algorithm respectively. Then, the high accuracy semi-analytical solution of direct problem can be got. Finally, based on the solution, the computing model of inverse problem and expression of sensitivity analysis are established. Single variable and variables combined identifications including thermal parameters, boundary conditions and source-related terms etc. are given to validate the approach proposed in 1-D and 2-D cases. The effects of noise data and initial guess on the results are investigated. The numerical examples show the effectiveness of this approach.
基金the National Natural Science Foundation of China(No.11672064)。
文摘A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.
基金supported by the National Natural Science Foundation of China(Grant Nos.12072288,U2241274,and 12272319).
文摘In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DNN)method.The NPIM is employed to modify the exponential matrix and loading vector based on the DNN integral method.This involves incorporating the basic assumption of the Newmark-βmethod into the dynamic equation and eliminating the acceleration term from the dynamic equilibrium equation.As a result,the equation is reduced to a first-order linear equation system.Subsequently,the PIM is applied to integrate the system step by step within the NPIM.The DNN method is adopted to solve the inhomogeneous term through fitting the integrand and the original function with a pair of neural networks,and the integral term is solved using the Newton–Leibniz formula.Numerical examples demonstrate that the proposed methodology significantly improves computing efficiency and provides sufficient precision compared to the DNN method.This is particularly evident when analyzing large-scale structures under blast loading conditions.
基金Supported by the National Natural Science Foundation of China(No.U22A20246)the Key Project of Natural Science Foundation of Hebei Province of China(Basic Research Base Project)(No.A2023210064)the Science and Technology Program of Hebei Province of China(Nos.246Z1904G and 225676162GH)。
文摘The vehicle-road coupling dynamics problem is a prominent issue in transportation,drawing significant attention in recent years.These dynamic equations are characterized by high-dimensionality,coupling,and time-varying dynamics,making the exact solutions challenging to obtain.As a result,numerical integration methods are typically employed.However,conventional methods often suffer from low computational efficiency.To address this,this paper explores the application of the parameter freezing precise exponential integrator to vehicle-road coupling models.The model accounts for road roughness irregularities,incorporating all terms unrelated to the linear part into the algorithm's inhomogeneous vector.The general construction process of the algorithm is detailed.The validity of numerical results is verified through approximate analytical solutions(AASs),and the advantages of this method over traditional numerical integration methods are demonstrated.Multiple parameter freezing precise exponential integrator schemes are constructed based on the Runge-Kutta framework,with the fourth-order four-stage scheme identified as the optimal one.The study indicates that this method can quickly and accurately capture the dynamic system's vibration response,offering a new,efficient approach for numerical studies of high-dimensional vehicle-road coupling systems.
基金supported by the Natural Science Foundation of China(Grant No.40574050,40821062)the National Basic Research Program of China(Grant No.2007CB209602)the Key Research Program of China National Petroleum Corporation(Grant No.06A10101)
文摘Seismic modeling is a useful tool for studying the propagation of seismic waves within complex structures. However, traditional methods of seismic simulation cannot meet the needs for studying seismic wavefields in the complex geological structures found in seismic exploration of the mountainous area in Northwestern China. More powerful techniques of seismic modeling are demanded for this purpose. In this paper, two methods of finite element-finite difference method （FE-FDM） and arbitrary difference precise integration （ADPI） for seismic forward modeling have been developed and implemented to understand the behavior of seismic waves in complex geological subsurface structures and reservoirs. Two case studies show that the FE-FDM and ADPI techniques are well suited to modeling seismic wave propagation in complex geology.
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.
基金This project is sponsored by the Specialized Prophasic Basic Research of the"973"Programme,contract No:2001cca02300
文摘D seismic modeling can be used to study the propagation of seismic wave exactly and it is also a tool of 3-D seismic data processing and interpretation. In this paper the arbitrary difference and precise integration are used to solve seismic wave equation, which means difference scheme for space domain and analytic integration for time domain. Both the principle and algorithm of this method are introduced in the paper. Based on the theory, the numerical examples prove that this hybrid method can lead to higher accuracy than the traditional finite difference method and the solution is very close to the exact one. Also the seismic modeling examples show the good performance of this method even in the case of complex surface conditions and complicated structures.
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
基金the National Natural Science Foundation of China (No.19872016,19872017)the National Key Basic Research Special Foundation (G1999032805)the Foundation for University Key Teachers by the Ministry of Education of China
文摘This paper presents an improved precise integration algorithm fortransient analysis of heat transfer and some other problems. Theoriginal precise integration method is improved by means of the inve-rse accuracy analysis so that the parameter N, which has been takenas a constant and an independent pa- rameter without consideration ofthe problems in the original method, can be generated automaticallyby the algorithm itself.
基金Project supported by the National Natural Science Foundation of China(Nos.60273048and60174023).
文摘Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.
基金Project(52178101) supported by the National Natural Science Foundation of China。
文摘Earthquake is a kind of sudden and destructive random excitation in nature.It is significant to determine the probability distribution characteristics of the corresponding dynamic indicators to ensure the safety and the stability of structures when the intensive seismic excitation,the intensity of which is larger than 7,acts in train-bridge system.Firstly,the motion equations of a two-dimensional train-bridge system under the vertical random excitation of track irregularity and the vertical seismic acceleration are established,where the train subsystem is composed of 8 mutually independent vehicle elements with 48 degrees of freedom,while the single-span simple supported bridge subsystem is composed of 102D beam elements with 20 degrees of freedom on beam and 2 large mass degrees of freedom at the support.Secondly,Monte Carlo method and pseudo excitation method are adopted to analyze the statistical parameters of the system.The power spectrum density of random excitation is used to define a series of non-stationary pseudo excitation in pseudo excitation method and the trigonometric series of random vibration history samples in Monte Carlo method,respectively solved by precise integral method and Newmark-βmethod through the inter-system iterative procedure.Finally,the results are compared with the case under the weak seismic excitation,and show that the samples of vertical acceleration response of bridge and the offload factor of train obeys the normal distribution.In a high probability,the intensive earthquakes pose a greater threat to the safety and stability of bridges and trains than the weak ones.
文摘This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.
基金supported by the National Natural Science Foundation of China (No.10662003)Educational Commission of Guangxi Province of China (No.200807MS109)
文摘The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.
文摘The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.
基金the support of the National Natural Science Foundation of China(Grant Nos.11472067 and 51609034)the Science Foundation of Liaoning Province of China(No.2021-MS-119)+1 种基金the Dalian Youth Science and Technology Star Project(No.2018RQ06)the Fundamental Research Funds for the Central Universities(Grant No.DUT20GJ216).
文摘The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was discussed in this study.To efficiently and accurately compute the exponential of the dynamics state matrix and the matrix functions due to external loads,an adaptively filtered precise integration method was proposed,which inherits the high precision of the precise integrationmethod,improves the computational efficiency and saves the memory required.Moreover,the perturbation method was introduced to avoid repeated computations of matrix exponential and terms due to external loads.Based on the filtering and perturbation techniques,an adaptively filtered precise integration method considering perturbation for stochastic dynamics problems was developed.Two numerical experiments,including a model of phononic crystal and a bridge model considering random parameters,were performed to test the performance of the proposed method in terms of accuracy and efficiency.Numerical results show that the accuracy and efficiency of the proposed method are better than those of the existing precise integration method,the Newmark-βmethod and the Wilson-θmethod.
基金This work was supported by the Natural Science Foundation of Shaanxi Province(Grant No.2021KW-25)the Astronautics Supporting Technology Foundation of China(Grant No.2019-HT-XG)the Fundamental Research Funds for the Central Universities(Grant No.3102018ZY015).
文摘The estimation of the disturbance input acting on a vehicle from its given responses is an inverse problem.To overcome some of the issues related to ill-posed inverse problems,this work proposes a method of reconstructing the road roughness based on the Kalman filter method.A half-car model that considers both the vehicle and equipment is established,and the joint input-state estimation method is used to identify the road profile.The capabilities of this methodology in the presence of noise are numerically demonstrated.Moreover,to reduce the influence of the driving speed on the estimation results,a method of choosing the calculation frequency is proposed.A road vibration test is conducted to benchmark the proposed method.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
文摘The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.
