Hits:

Indexed by:期刊论文

Date of Publication:2012-01-27

Journal:INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING

Included Journals:SCIE、EI、Scopus

Volume:89

Issue:4

Page Number:438-508

ISSN No.:0029-5981

Key Words:symplectic; Hamiltonian system; principle of least action; generating function; dual

Abstract:In this paper, the principle of least action and generating functions are used to construct symplectic numerical algorithms for finite dimensional autonomous Hamiltonian systems. The approximate action is obtained by approximating the generalized coordinates and momentums by Lagrange polynomials and performing Gaussian quadrature. Based on the principle of least action and the requirements of a canonical transformation, different types of symplectic algorithms have been constructed by choosing different types of independent variables at two ends of the time step. The symmetric property of the four types of symplectic algorithms proposed in this paper is discussed, and the exact linear stability domain for small m, n and g is discussed. The linear stability and precision of different types of symplectic algorithms are tested using numerical examples. Copyright (c) 2011 John Wiley & Sons, Ltd.