刚体动力学方程的一个辛积分方法

摘要

针对四元数和对应广义动量表示的刚体定点动力学方程,利用一种位移格式的微分—代数方程积分方案,实现了非独立广义动量表示的拉格朗日方程的辛积分算法.数值实验显示该算法具有精度高和保持系统守恒量的特点.更为重要的是,广义动量表示的拉格朗日方程较之传统形式的拉格朗日方程在辛积分中表现出独特的优越性.