Force‐stepping integrators in Lagrangian mechanics
Force‐stepping integrators in Lagrangian mechanics
复制标题
拉格朗日力学中的力步进积分器
DOI:
--
复制
发表时间:
2010
期刊:
影响因子:
--
通讯作者:
Magdalena Ortiz
中科院分区:
文献类型:
--
作者:
M. González;B. Schmidt;Magdalena Ortiz
We formulate an integration scheme for Lagrangian mechanics, referred to as the force‐stepping scheme, which is symplectic, energy conserving, time‐reversible, and convergent with automatic selection of the time‐step size. The scheme also conserves approximately all the momentum maps associated with the symmetries of the system. The exact conservation of momentum maps may additionally be achieved by recourse to the Lagrangian reduction. The force‐stepping scheme is obtained by replacing the potential energy by a piecewise affine approximation over a simplicial grid or regular triangulation. By taking triangulations of diminishing size, an approximating sequence of energies is generated. The trajectories of the resulting approximate Lagrangians can be characterized explicitly and consist of piecewise parabolic motion, or free fall. Selected numerical tests demonstrate the excellent long‐term behavior of force‐stepping, its automatic time‐step selection property, and the ease with which it deals with constraints, including contact problems. Copyright © 2010 John Wiley & Sons, Ltd.