[論文レビュー] Geometric integrators for higher-order mechanics on Lie groups

Christopher L. Burnett, Darryl D. Holm+1 more|arXiv (Cornell University)|Dec 28, 2011

Numerical methods for differential equations参考文献 47被引用数 13

ひとこと要約

本稿では、加速に依存するラグランジアンを有する高次力学系に対する、リー群上での構造保存型数値積分法を提示する。この手法は、ハミルトン＝ポントリャーギン変分積分法フレームワークを高次力学系へ拡張したものであり、ルンゲ＝クッタ＝ムンテ＝カース法と離散変分原理を組み合わせることで、シンプレクティック性と運動量保存を保持する。この手法は、$S^2$ 上のリーマン多様体立方曲線の補間問題に成功裏に適用され、運動量の正確な保存と高精度な軌道再構成が達成された。

ABSTRACT

This paper develops a structure-preserving numerical integration scheme for a class of higher-order mechanical systems. The dynamics of these systems are governed by invariant variational principles defined on higher-order tangent bundles of Lie groups. The variational principles admit Lagrangians that depend on acceleration, for example. The symmetry reduction method used in the Hamilton--Pontryagin approach for developing variational integrators of first-order mechanics is extended here to higher order. The paper discusses the general approach and then focuses on the primary example of Riemannian cubics. Higher-order variational integrators are developed both for the discrete-time integration of the initial value problem and for a particular type of trajectory-planning problem. The solution of the discrete trajectory-planning problem for higher-order interpolation among points on the sphere illustrates the approach.

研究の動機と目的

• リー群上での高次力学系のための変分積分法を構築し、シンプレクティック性や運動量保存といった幾何的構造を保存すること。
• ハミルトン＝ポントリャーギン変分形式を一次力学系から高次力学系へ拡張し、特に加速に依存するラグランジアンを有する系に適用すること。
• 得られた積分法を、$SO(3)$ 上での回転の補間問題といった軌道計画問題に応用し、計算解剖学やロボット工学への応用を図ること。
• 球面 $S^2$ 上での離散的軌道計画問題の数値解法を通じて、本手法の精度と幾何的忠実性を実証すること。
• 加速度を超える高次導関数に対しても一般化可能であり、滑らかな補間や最適制御への広範な応用を可能にする。

提案手法

• 連続時間において、2階の運動学的制約をラグランジュ乗数を用いて強制する高次ハミルトン＝ポントリャーギン変分原理を採用する。
• 離散時間における力学系は、ルンゲ＝クッタ＝ムンテ＝カース（RKMK）法と変分積分法原理を組み合わせることで、連続系の幾何的構造を保存する。
• 離散的な流れ写像は、リー群の高次接バンドル上での変分原理から導出され、シンプレクティック性と運動量保存を保証する。
• 離散的随伴運動と運動量写像を組み込み、ノード点における正確な運動量保存と中間段階における連続性を実現する。
• 境界条件およびノード制約を満たすようにコスト関数を最小化するためのシューティング法を用いて離散的軌道計画問題を解く。
• 初期運動量空間上で勾配降下法を用いて実装し、平均二乗共変微分加速度を最小化するコスト関数を基盤とする。

実験結果

リサーチクエスチョン

• RQ1 どのようにして、リー群上での一次力学系から高次力学系への変分積分法を拡張できるか？
• RQ2 加速に依存するラグランジアンを有する高次力学系において、シンプレクティック性と運動量保存の両方を保持する離散数値スキームは何か？
• RQ3 ハミルトン＝ポントリャーギンフレームワークを高次系に適応させることで、$S^2$ のような多様体上での構造保存型軌道計画が可能になるか？
• RQ4 リーマン多様体立方曲線の幾何的性質（滑らかさ、曲率など）は、リー群上での離散積分法にどのように反映されるか？
• RQ5 所定の中間回転姿勢を満たす補間問題を解く際、得られた高次積分法の精度と安定性はどの程度か？

主な発見

• 提案された高次ハミルトン＝ポントリャーギン積分法は、離散時間においてシンプレクティック形式と運動量写像を正確に保存し、長期的な幾何的忠実性を確保する。
• $S^2$ 上の補間例を通じて、ノードインデックスにおける正確な運動量保存と中間ステップにおける運動量の連続性が確認された。
• 数値結果から、$S^2$ 上の離散的補間曲線が、$\sigma = 0.025$ の許容誤差内で所定の目標点を正確に通過することが示された。
• 運動量ノルム $\|\mu_k\|$ と $\|\nu_k\|$ はノード点で正確に保存され、区間間で連続的であることが確認され、理論的な構造保存性が裏付けられた。
• 離散的軌道は、色分けされた速度で可視化され、滑らかな遷移と不自然な振動の欠如が確認され、高品質な補間が達成された。
• 本手法は加速度を超える高次導関数に対しても容易に一般化可能であり、今後のジャーブやスナップ最小化軌道への応用が期待できる。