An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions ...An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.展开更多
As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely use...As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.展开更多
Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy...Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.展开更多
The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numerical integration of an orbit on the unstable manifold of a periodic solution.This algorithm is matrix-free and emp...The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numerical integration of an orbit on the unstable manifold of a periodic solution.This algorithm is matrix-free and employs a combination of the Newton-Raphson method and the Krylov subspace method.Moreover,the algorithm adopts a multiple shooting method to address the problem of orbital instability due to long-term numerical integration.The algorithm is described through computing the extension of unstable manifold of a recomputed Nagata′s lowerbranch steady solution of plane Couette flow,which is an example of an exact coherent state that has recently been studied in subcritical transition to turbulence.展开更多
With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and re...With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.展开更多
This paper deals with the parallel information-based complexity of numerical integration on Sobolev class. We obtain tight bounds on the complexity, considered as a function of two variables simultaneously:the number ...This paper deals with the parallel information-based complexity of numerical integration on Sobolev class. We obtain tight bounds on the complexity, considered as a function of two variables simultaneously:the number of processors, the rquired precision. This result seems to be new even in serial case.展开更多
In this paper,we introduce a type of tensor neural network.For the first time,we propose its numerical integration scheme and prove the computational complexity to be the polynomial scale of the dimension.Based on the...In this paper,we introduce a type of tensor neural network.For the first time,we propose its numerical integration scheme and prove the computational complexity to be the polynomial scale of the dimension.Based on the tensor product structure,we develop an efficient numerical integration method by using fixed quadrature points for the functions of the tensor neural network.The corresponding machine learning method is also introduced for solving high-dimensional problems.Some numerical examples are also provided to validate the theoretical results and the numerical algorithm.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
Numerical integration over the implicitly defined domains is challenging due to topological variances of implicit functions.In this paper,we use interval arithmetic to identify the boundary of the integration domain e...Numerical integration over the implicitly defined domains is challenging due to topological variances of implicit functions.In this paper,we use interval arithmetic to identify the boundary of the integration domain exactly,thus getting the correct topology of the domain.Furthermore,a geometry-based local error estimate is explored to guide the hierarchical subdivision and save the computation cost.Numerical experiments are presented to demonstrate the accuracy and the potential of the proposed method.展开更多
Several numerical integration schemes for the evaluation of matrix elements in density functional theory calculations have been studied and compared by computational practice. The best scheme was found to be the combi...Several numerical integration schemes for the evaluation of matrix elements in density functional theory calculations have been studied and compared by computational practice. The best scheme was found to be the combination of the atomic partition function proposed by Becke with the scaled generalized Gauss-Laguerre quadrature formula for radial integration suggested by Yang, which achieve the highest convergence rate to the numerical integration. With the same number of integration points, the accuracy of the calculated results by this scheme is higher by 1 to 2 orders of magnitudes than that by other schemes. The reason for achieving higher accuracy by this scheme has been proposed preliminarily.展开更多
The numerical evaluation of an integral is a frequently encountered problem in antenna analysis. A special Gauss Christoffel quadrature formula for nonclassical weight function is constructed for solving the pseu...The numerical evaluation of an integral is a frequently encountered problem in antenna analysis. A special Gauss Christoffel quadrature formula for nonclassical weight function is constructed for solving the pseudo singular integration problem arising from the analysis of thin wire antennas. High integration accuracy is obtained at comparable low computation cost by the quadrature formula constructed. This integration method can be also used in other electromagnetic integral equation problems.展开更多
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.展开更多
The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise w...The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.展开更多
An edge wave is a kind of surface gravity wave basically travelling along a shoaling beach. Based on the periodic assumption in the longshore direction, a second order ordinary differential equation is obtained for nu...An edge wave is a kind of surface gravity wave basically travelling along a shoaling beach. Based on the periodic assumption in the longshore direction, a second order ordinary differential equation is obtained for numerical simulation of the cross-shore surface elevation. Given parameters at the shoreline, a cross-shore elevation profile is obtained through integration with fourth-order Runge Kutta technique. For a compound slope, a longshore wavenumber is obtained by following a geometrical approach and solving a transcendental equation with an asymptotic method. Numerical results on uniform and compound sloping beaches with different wave periods, slope angles, modes and turning point positions are presented. Some special scenarios, which cannot be predicted by analytical models are also discussed.展开更多
Hybrid simulation is a powerful test method for evaluating the seismic performance of structural systems. This method makes it feasible that only critical components of a structure be experimentally tested. This paper...Hybrid simulation is a powerful test method for evaluating the seismic performance of structural systems. This method makes it feasible that only critical components of a structure be experimentally tested. This paper presents a newly proposed integration algorithm for seismic hybrid simulation which is aimed to extend its capabilities to a wide range of systems where existing methods encounter some limitations. In the proposed method, which is termed the variable time step （VTS） integration method, an implicit scheme is employed for hybrid simulation by eliminating the iterative phase on experimental element, the phase which is necessary in regular implicit applications. In order to study the effectiveness of the VTS method, a series of numerical investigations are conducted which show the successfulness of the VTS method in obtaining accurate, stable and converged responses. Then, in a comparative approach, the improved accuracy of the VTS method over commonly used integration methods is demonstrated. The stability of the VTS method is also studied and the results show that it provides conditional stability; however, its stability limit is well beyond the accuracy limit. The effect of time delay on the VTS method results is also investigated and it is shown that the VTS method is quite successful in handling this experimental error.展开更多
A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can ...A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can easily be got by directly solving linear matrixequations rather than structure motion differential equations. Moreover, in order to solve thecorresponding linear matrix equations, the numerical integration fast algorithm is presented. Thenaccording to the results, dynamic design and life-span estimation can be done. Besides, the newalgorithm can solve non-proportion damp structure response.展开更多
This paper focuses on studying a new energy-work relationship numericM integration scheme of nonholonomic Hamiltonian systems. The signal-stage numerical, multi-stage and parallel composition numerical integration sch...This paper focuses on studying a new energy-work relationship numericM integration scheme of nonholonomic Hamiltonian systems. The signal-stage numerical, multi-stage and parallel composition numerical integration schemes are presented. The high-order energy-work relation scheme of the system is constructed by a parallel connection of n multi-stage schemes of order 2, its order of accuracy is 2n. The connection, which is discrete analogue of usual case, between the change of energy and work of nonholonomic constraint forces is obtained for nonholonomie Hamiltonian systems. This paper also gives that there is smaller error of the scheme when taking a large number of stages than a less one. Finally, an applied example is discussed to illustrate these results.展开更多
If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restric...If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.展开更多
A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with...A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with explicit and predict-correct and self-starting and fourth-order accuracy to integrate the integral equations is given. Theoretical analysis and numerical examples shaw that DIM-IM described in this paper suitable for strong nonlinear and non-conservative system have higher accuracy than central difference, Houbolt, Newmark and Wilson-Theta methods.展开更多
文摘An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.
基金supported by National Key Basic Research Program (973 Program, Grant No. 2011CB706804)National Natural Science Foundation of China (Grant No. 50805093)Science & Technology Commission of Shanghai Municipality, China (Grant No. 09QH1401500)
文摘As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.
基金National Natural Science Foundation of China under Grant Nos.51639006 and 51725901
文摘Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.
文摘The recent result of an orbit continuation algorithm has provided a rigorous method for long-term numerical integration of an orbit on the unstable manifold of a periodic solution.This algorithm is matrix-free and employs a combination of the Newton-Raphson method and the Krylov subspace method.Moreover,the algorithm adopts a multiple shooting method to address the problem of orbital instability due to long-term numerical integration.The algorithm is described through computing the extension of unstable manifold of a recomputed Nagata′s lowerbranch steady solution of plane Couette flow,which is an example of an exact coherent state that has recently been studied in subcritical transition to turbulence.
文摘With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.
基金this work was supported by china State Major Key Project for Basic Researchers
文摘This paper deals with the parallel information-based complexity of numerical integration on Sobolev class. We obtain tight bounds on the complexity, considered as a function of two variables simultaneously:the number of processors, the rquired precision. This result seems to be new even in serial case.
基金supported in part by the National Key Research and Development Program of China(Grant No.2019YFA0709601)by the National Center for Mathematics and Interdisciplinary Science,CAS.
文摘In this paper,we introduce a type of tensor neural network.For the first time,we propose its numerical integration scheme and prove the computational complexity to be the polynomial scale of the dimension.Based on the tensor product structure,we develop an efficient numerical integration method by using fixed quadrature points for the functions of the tensor neural network.The corresponding machine learning method is also introduced for solving high-dimensional problems.Some numerical examples are also provided to validate the theoretical results and the numerical algorithm.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
基金The work is supported by the National Natural Science Foundation of China(No.11771420).
文摘Numerical integration over the implicitly defined domains is challenging due to topological variances of implicit functions.In this paper,we use interval arithmetic to identify the boundary of the integration domain exactly,thus getting the correct topology of the domain.Furthermore,a geometry-based local error estimate is explored to guide the hierarchical subdivision and save the computation cost.Numerical experiments are presented to demonstrate the accuracy and the potential of the proposed method.
基金Supported by State Major Key Project for Basic Researches and the National Natural Science Foundation of China.
文摘Several numerical integration schemes for the evaluation of matrix elements in density functional theory calculations have been studied and compared by computational practice. The best scheme was found to be the combination of the atomic partition function proposed by Becke with the scaled generalized Gauss-Laguerre quadrature formula for radial integration suggested by Yang, which achieve the highest convergence rate to the numerical integration. With the same number of integration points, the accuracy of the calculated results by this scheme is higher by 1 to 2 orders of magnitudes than that by other schemes. The reason for achieving higher accuracy by this scheme has been proposed preliminarily.
文摘The numerical evaluation of an integral is a frequently encountered problem in antenna analysis. A special Gauss Christoffel quadrature formula for nonclassical weight function is constructed for solving the pseudo singular integration problem arising from the analysis of thin wire antennas. High integration accuracy is obtained at comparable low computation cost by the quadrature formula constructed. This integration method can be also used in other electromagnetic integral equation problems.
基金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.
文摘The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.
基金financially supported by the National Natural Science Foundation of China(Grant No.51279055)the Fundamental Research Funds for the Central Universities(Grant No.2015B35114)the Open Fund of Jiangsu Key Laboratory of Coast Ocean Resources Development and Environment Security of Hohai University(Grant No.201506)
文摘An edge wave is a kind of surface gravity wave basically travelling along a shoaling beach. Based on the periodic assumption in the longshore direction, a second order ordinary differential equation is obtained for numerical simulation of the cross-shore surface elevation. Given parameters at the shoreline, a cross-shore elevation profile is obtained through integration with fourth-order Runge Kutta technique. For a compound slope, a longshore wavenumber is obtained by following a geometrical approach and solving a transcendental equation with an asymptotic method. Numerical results on uniform and compound sloping beaches with different wave periods, slope angles, modes and turning point positions are presented. Some special scenarios, which cannot be predicted by analytical models are also discussed.
文摘Hybrid simulation is a powerful test method for evaluating the seismic performance of structural systems. This method makes it feasible that only critical components of a structure be experimentally tested. This paper presents a newly proposed integration algorithm for seismic hybrid simulation which is aimed to extend its capabilities to a wide range of systems where existing methods encounter some limitations. In the proposed method, which is termed the variable time step （VTS） integration method, an implicit scheme is employed for hybrid simulation by eliminating the iterative phase on experimental element, the phase which is necessary in regular implicit applications. In order to study the effectiveness of the VTS method, a series of numerical investigations are conducted which show the successfulness of the VTS method in obtaining accurate, stable and converged responses. Then, in a comparative approach, the improved accuracy of the VTS method over commonly used integration methods is demonstrated. The stability of the VTS method is also studied and the results show that it provides conditional stability; however, its stability limit is well beyond the accuracy limit. The effect of time delay on the VTS method results is also investigated and it is shown that the VTS method is quite successful in handling this experimental error.
基金This project is supported by National Natural Science Foundation of China (No.59805001)
文摘A new algorithm of structure random response numerical characteristics, namedas matrix algebra algorithm of structure analysis is presented. Using the algorithm, structurerandom response numerical characteristics can easily be got by directly solving linear matrixequations rather than structure motion differential equations. Moreover, in order to solve thecorresponding linear matrix equations, the numerical integration fast algorithm is presented. Thenaccording to the results, dynamic design and life-span estimation can be done. Besides, the newalgorithm can solve non-proportion damp structure response.
基金supported by National Natural Science Foundation of China under Grant No.10672143the Natural Science Foundation of Henan Province under Grant No.0511022200
文摘This paper focuses on studying a new energy-work relationship numericM integration scheme of nonholonomic Hamiltonian systems. The signal-stage numerical, multi-stage and parallel composition numerical integration schemes are presented. The high-order energy-work relation scheme of the system is constructed by a parallel connection of n multi-stage schemes of order 2, its order of accuracy is 2n. The connection, which is discrete analogue of usual case, between the change of energy and work of nonholonomic constraint forces is obtained for nonholonomie Hamiltonian systems. This paper also gives that there is smaller error of the scheme when taking a large number of stages than a less one. Finally, an applied example is discussed to illustrate these results.
基金National Natural Science Foundation of China (50178065), 973 Program (2002CB412706), National Social Com-monweal Research Foundation (2002DIB30076) and Joint Seismological Science Foundation (101066).
文摘If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.
文摘A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with explicit and predict-correct and self-starting and fourth-order accuracy to integrate the integral equations is given. Theoretical analysis and numerical examples shaw that DIM-IM described in this paper suitable for strong nonlinear and non-conservative system have higher accuracy than central difference, Houbolt, Newmark and Wilson-Theta methods.
