[论文解读] Dual Quaternion Variational Integrator for Rigid Body Dynamic Simulation

Jiafeng Xu, Karl Henning Halse|arXiv (Cornell University)|Nov 2, 2016

Numerical methods for differential equations参考文献 19被引用 9

一句话总结

本文提出了双四元数变分积分器（DQVI），一种保持辛结构的单步李群积分器，利用双四元数统一表示旋转与平动，实现对刚体在六自由度运动中的模拟。通过保持系统的几何结构、能量与动量，DQVI能够在任意定义的参考点和完整的6×6惯性矩阵下，实现高精度的实时模拟，且时间步长较大。

ABSTRACT

We introduce a symplectic dual quaternion variational integrator(DQVI) for simulating single rigid body motion in all six degrees of freedom. Dual quaternion is used to represent rigid body kinematics and one-step Lie group variational integrator is used to conserve the geometric structure, energy and momentum of the system during the simulation. The combination of these two becomes the first Lie group variational integrator for rigid body simulation without decoupling translations and rotations. Newton-Raphson method is used to solve the recursive dynamic equation. This method is suitable for real-time rigid body simulations with high precision under large time step. DQVI respects the symplectic structure of the system with excellent long-term conservation of geometry structure, momentum and energy. It also allows the reference point and 6-by-6 inertia matrix to be arbitrarily defined, which is very convenient for a variety of engineering problems.

研究动机与目标

• 开发一种几何积分器，以在模拟过程中保持刚体系统的辛结构、能量与动量。
• 通过特殊欧氏群SE(3)统一处理刚体动力学中的平动与旋转，避免其解耦。
• 在保持高精度与数值稳定性的同时，实现大时间步长的实时模拟。
• 允许任意定义参考点与完整的6×6惯性矩阵，以增强其在工程应用中的广泛适用性。
• 为多步积分与未来扩展至含势场及非保守力的系统提供基础。

提出的方法

• 使用双四元数在SE(3)中表示刚体位姿，提供无奇点、紧凑且统一的旋转与平动表示。
• 采用单步李群变分积分器，确保系统在李群上自然演化，无需重新投影或施加约束。
• 通过六个独立变量（Φ与Ψ）对李代数se(3)的指数映射进行参数化，确保配置始终位于SE(3)上。
• 基于变分原理推导离散动力学方程，保持连续系统的几何结构。
• 使用牛顿-拉夫森法求解所得非线性方程，并以增量位姿变化作为初始估计。
• 采用残差方程的解析雅可比表达式，以确保快速且鲁棒的收敛。

实验结果

研究问题

• RQ1 能否构建一种变分积分器，将完整的6自由度刚体运动作为单一李群处理，而无需解耦平动与旋转？
• RQ2 在李群变分积分器中使用双四元数，是否相比传统方法能更好地保持能量、动量与几何结构的守恒？
• RQ3 该积分器是否能在大时间步长下仍保持高精度与稳定性，同时支持任意参考点与6×6惯性矩阵？
• RQ4 所提出的初始估计策略下，牛顿-拉夫森求解器在实现快速收敛方面的有效性如何？
• RQ5 该方法未来能否扩展至包含非保守力与势场的系统？

主要发现

• DQVI方法保持了系统的辛结构，从而实现了能量与动量的优异长期守恒。
• 该方法允许任意定义参考点与完整的6×6惯性矩阵，显著增强了其在工程问题（如水面航行器的附加质量和阻尼）中的适用性。
• 当初始值设置得当时，牛顿-拉夫森求解器可在仅三步迭代内收敛至机器精度，从而实现高效的实时模拟。
• 基于增量位姿变化的初始估计确保了参数化的有效性（|Φ| < 1），从而将时间步长限制在旋转不超过180度的范围内。
• 当仅发生平动且参考点位于质心时，牛顿-拉夫森求解器退化为单步前向欧拉法，验证了其在极限情况下的自洽性。
• DQVI是首个将SE(3)作为整体处理的刚体动力学李群变分积分器，未对旋转与平动动力学进行解耦。