This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order...This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a ＂major compo- nent＂ and a ＂high order small quantity＂ and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.展开更多
Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved thro...Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.展开更多
This paper presents an efficient time-integration method for obtaining reliable solutions to the second-order nonlinear dynamic problems in structural engineering. This method employs both the backward-acceleration di...This paper presents an efficient time-integration method for obtaining reliable solutions to the second-order nonlinear dynamic problems in structural engineering. This method employs both the backward-acceleration differentiation formula and the trapezoidal rule, resulting in a self-starting, single step, second-order accurate algorithm. With the same computational effort as the trapezoidal rule, the proposed method remains stable in large deformation and long time range solutions even when the trapezoidal rule fails. Meanwhile, the proposed method has the following characteristics: (1) it is applicable to linear as well as general nonlinear analyses; (2) it does not involve additional variables (e.g. Lagrange multipliers) and artificial parameters; (3) it is a single-solver algorithm at the discrete time points with symmetric effective stiffness matrix and effective load vectors; and (4) it is easy to implement in an existing computational software. Some numerical results indicate that the proposed method is a powerful tool with some notable features for practical nonlinear dynamic analyses.展开更多
基金supported by the National Natural Science Foun-dation of China (11172334)
文摘This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a ＂major compo- nent＂ and a ＂high order small quantity＂ and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.
基金Project supported by the National Basic Research Program of China (973 program) (No.G1999032804)
文摘Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.
基金sponsored by the Scientific Foundation for Returned Oversea Scholars of China (Grant No.20101020044)the State Key Laboratory of Hydro–Science and Engineering (Grant Nos. 2008Z6 and 2009-TC-2)
文摘This paper presents an efficient time-integration method for obtaining reliable solutions to the second-order nonlinear dynamic problems in structural engineering. This method employs both the backward-acceleration differentiation formula and the trapezoidal rule, resulting in a self-starting, single step, second-order accurate algorithm. With the same computational effort as the trapezoidal rule, the proposed method remains stable in large deformation and long time range solutions even when the trapezoidal rule fails. Meanwhile, the proposed method has the following characteristics: (1) it is applicable to linear as well as general nonlinear analyses; (2) it does not involve additional variables (e.g. Lagrange multipliers) and artificial parameters; (3) it is a single-solver algorithm at the discrete time points with symmetric effective stiffness matrix and effective load vectors; and (4) it is easy to implement in an existing computational software. Some numerical results indicate that the proposed method is a powerful tool with some notable features for practical nonlinear dynamic analyses.
