A High Order Symplectic Algorithm Based on Weighted Residual Method

LU Kelang,FU Minghui,LI Weihua,LI Renfei
DOI: https://doi.org/10.13471/j.cnki.acta.snus.2015.04.002
2015-01-01
Abstract:A new way to construct high order symplectic algorithms is proposed based on weighted residu-al method.Firstly,in the time subdomain,the corresponding integral equation of Galerkin method for Hamilton dual equation based on the idea of weighted residual method is proposed,then the generalized displacement and momentum are approximated by the same Lagrange interpolation within the time subdo-main,which are substituted into the corresponding integral equation.By numerical integration,the origi-nal initial value problem of dynamics is expressed as algebraic equations with displacement and momen-tum at the interpolation points as unknown variables.For nonlinear dynamic systems,a simple scheme of choosing initial values,which can significantly improve the computational efficiency for Newton-Raphson method,is presented.Finally,the symplecticity and performance of the proposed algorithms are dis-cussed in detail.Compared with the same order symplectic Runge-Kutta methods,the accuracy of the two methods are almost the same,but the proposed algorithms are much simpler and less computational ex-pense.The numerical results illustrate that the proposed algorithms show good performance in accuracy and efficiency.
What problem does this paper attempt to address?