[論文レビュー] 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

ひとこと要約

本稿では、二重クォータニオンを用いて回転と並進を統合することで、6自由度の剛体運動をシミュレートするシンプレクティックな1ステップリー群積分法、Dual Quaternion Variational Integrator (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) 上の剛体ポーズを二重クォータニオンで表現し、特異点のない、コンactかつ統合的な回転と並進の表現を実現する。
• リーマン多様体上での自然な運動を保証する1ステップリー群変分積分法を採用し、再投影や制約なしにリー群上を自然に進化させる。
• リー代数 se(3) の指数写像を表す6つの独立変数（Φ と Ψ）を用いたパarametrization を適用し、配置が常に SE(3) 上に保たれるようにする。
• 変分原理を用いて離散的運動方程式を導出し、連続系の幾何的構造を保存する。
• 得られた非線形方程式を、インクリメンタルなポーズ変化から得られる初期推定値を用いてニュートン・ラプソン法で解く。
• 残差方程式のヤコビアンに対して解析的表現を用いることで、高速かつ安定した収束を実現する。

実験結果

リサーチクエスチョン

• RQ1 6自由度の剛体運動を回転と並進を分離せずに1つのリー群として扱う変分積分法を構築できるか？
• RQ2 リー群変分積分法に二重クォータニオンを用いることで、従来の手法と比較してエネルギー、運動量、幾何的構造の保存性が向上するか？
• RQ3 任意の基準点と6×6の慣性マトリクスを許容しながらも、大きな時間刻み幅でも高精度と安定性を維持できるか？
• RQ4 提案された初期推定戦略を用いたニュートン・ラプソン法は、高速収束を達成できるか？
• RQ5 今後の研究において、非保存力やポテンシャル場を含む系への拡張は可能か？

主な発見

• DQVI法は系のシンプレクティック構造を保存し、長期的にエネルギーと運動量の優れた保存性を実現する。
• 基準点と6×6の慣性マトリクスを任意に定義可能であり、海洋ビークルにおける追加質量や減衰の問題など、工学的応用への適用性が向上する。
• 適切に初期化された場合、ニュートン・ラプソン法は最悪でも3回の反復で機械精度に収束し、効率的なリアルタイムシミュレーションを可能にする。
• インクリメンタルポーズ変化に基づく初期推定により、パarametrization が有効（|Φ| < 1）を保ち、回転角が180度未満に制限されるように時間刻み幅が制御される。
• 質量中心に基準点を置いた純粋な並進運動の場合、ニュートン・ラプソン法は1ステップの前進オイラー法に帰着し、極限における一貫性が確認される。
• DQVIは、回転と並進の動的挙動を分離せずにSE(3)全体を一つの対象として扱う、剛体動力学における最初のリー群変分積分法である。