[论文解读] 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$ 上求解离散轨迹规划问题，展示该方法的精度与几何保真度。
• 将该方法推广至加速度以外的更高阶导数，以拓展其在光滑插值与最优控制中的广泛应用。

提出的方法

• 该方法在连续时间中采用高阶哈密顿-庞特里亚金变分原理，并通过拉格朗日乘子强制实施二阶运动学约束。
• 通过将龙格-库塔-芒特-卡斯（RKMK）方法与变分积分原理相结合，构建离散时间动力学，以保持连续系统的几何结构。
• 离散流映射源自李群高阶切从丛上的变分原理，确保辛结构与动量守恒。
• 该格式整合了离散伴随运动与动量映射，确保在节点处实现精确动量守恒，并在中间阶段保持连续性。
• 采用打靶法求解离散轨迹规划问题，通过在边界条件与节点约束下最小化代价泛函实现。
• 实现采用初始动量空间上的梯度下降法，代价泛函基于最小化平均平方协变加速度。

实验结果

研究问题

• RQ1 如何将变分积分方法从一阶系统扩展至李群上的高阶机械系统？
• RQ2 何种离散数值格式可在具有加速度依赖性拉格朗日量的高阶力学中同时保持辛结构与动量守恒？
• RQ3 哈密顿-庞特里亚金框架能否被适配至高阶系统，以实现在 $S^2$ 等流形上的结构保持轨迹规划？
• RQ4 黎曼立方的几何特性（如光滑性与曲率）如何在李群上的离散积分器中体现？
• RQ5 该高阶积分方法在求解具有指定中间姿态的插值问题时，其精度与稳定性如何？

主要发现

• 所提出的高阶哈密顿-庞特里亚金积分方法在离散时间中成功保持了辛形式与动量映射，确保长期几何保真度。
• 该方法在节点索引处实现精确动量守恒，并在中间步骤中保持动量连续性，如 $S^2$ 插值示例所验证。
• 数值结果表明，$S^2$ 上的离散插值器以 $\sigma = 0.025$ 的容差精确通过预设目标点。
• 动量范数 $\|\mu_k\|$ 与 $\|\nu_k\|$ 在节点处实现精确守恒，并在区间间保持连续，证实了理论上的结构保持性。
• 离散轨迹以颜色编码速度进行可视化，显示平滑过渡且无虚假振荡，表明插值质量优异。
• 该方法可直接推广至加速度以外的更高阶导数，为未来实现加加加速度（jerk）或加加加加速度（snap）最小化轨迹提供支持。