A variational integrator of arbitrarily high-order on the special orthogonal group $SO(n)$ is constructed using the polar decomposition and the constrained Galerkin method. It has the advantage of avoiding the second-order derivative of the exponential map that arises in traditional Lie group variational methods. In addition, a reduced Lie--Poisson integrator is constructed and the resulting algorithms can naturally be implemented by fixed-point iteration. The proposed methods are validated by numerical simulations on $SO(3)$ which demonstrate that they are comparable to variational Runge--Kutta--Munthe-Kaas methods in terms of computational efficiency. However, the methods we have proposed preserve the Lie group structure much more accurately and and exhibit better near energy preservation.
翻译:使用极地分解法和受限制的Galerkin法,构建了特殊正统组(USSO(n))任意高顺序的变异集成体,其优点是避免了传统立伊组变异法中产生的指数图的二阶衍生物,此外,还建造了一个减低的利皮-普瓦松集成体,由此产生的算法自然可以通过固定点迭代法加以实施。提议的方法通过对USSO(3)美元的数值模拟加以验证,这些模拟表明,在计算效率方面,这些方法与变异的龙格-库塔-蒙特-卡亚斯方法相仿。然而,我们提议的方法更精确地维护利伊组结构,并在靠近能源保护的地方展示得更好。