[论文解读] On the simulation of multibody systems with holonomic constraints

Jukka Tuomela, Teijo Arponen+1 more|ArXiv.org|Nov 3, 2006

Dynamics and Control of Mechanical Systems参考文献 14被引用 3

一句话总结

本文提出一种基于喷射丛（jet space）的拉格朗日计算模型，用于模拟具有完整约束的多体系统，通过将约束投影到约束流形上，避免了漂移和虚假振荡。通过在喷射丛中嵌入动力学并几何地施加约束，该方法使显式龙格-库塔时间积分结合自适应步长控制，能够实现长期的物理一致模拟。

ABSTRACT

We use Lagrangian formalism and jet spaces to derive a computational model to simulate multibody dynamics with holonomic constraints. Our approach avoids the traditional problems of drift-off and spurious oscillations. Hence even long simulations remain physically relevant. We illustrate our method with several numerical examples.

研究动机与目标

• 解决多体系统在完整约束下长期模拟中持续存在的约束漂移与虚假振荡问题。
• 开发一种基于微分几何与拉格朗日力学的几何计算模型，以确保长时间模拟中的物理一致性。
• 将守恒定律（如能量不变性）自然地嵌入约束流形结构中，以保持物理不变量。
• 通过消除传统稳定化技术引入的刚性，使显式龙格-库塔方法能够实现高效且稳定的数值积分。
• 在典型多体系统（包括平面四边形与曲柄机构）上验证该方法的鲁棒性与效率。

提出的方法

• 在喷射丛框架内使用拉格朗日形式化表述动力学，实现对高阶导数与约束的几何处理。
• 将配置流形定义为喷射丛的子流形，通过隐式流形嵌入限制动力学，使其满足完整约束与不变量。
• 在每个积分步骤中应用正交投影以强制满足约束，防止漂移并保持物理一致性。
• 采用准正交投影策略加速约束投影过程中的收敛速度，相比标准正交投影显著降低计算成本。
• 基于 Dormand-Prince 龙格-库塔方法实现自适应步长控制，误差容差设定为（Atol=Rtol=10⁻⁷），以平衡精度与效率。
• 利用自动微分高效计算约束函数与不变量函数的梯度及高阶导数。

实验结果

研究问题

• RQ1 基于喷射丛的几何化多体动力学公式能否有效消除长期模拟中的约束漂移与虚假振荡？
• RQ2 在 CPU 时间与收敛性方面，正交投影的计算成本与传统稳定化方法相比如何？
• RQ3 在不损害数值效率的前提下，守恒定律（如能量不变性）能在多大程度上自然地嵌入约束流形中？
• RQ4 与标准正交投影相比，准正交投影策略在速度与迭代次数方面如何提升性能？
• RQ5 当采用基于投影的约束强制方法时，显式龙格-库塔方法能否在约束多体模拟中有效应用？

主要发现

• 该方法成功防止了漂移与虚假振荡，实现了具有完整约束的多体系统在长期模拟中的物理一致性。
• 在曲柄机构模拟中，准正交投影使投影的 CPU 时间减少了 75%（从 3198.90 秒降至 757.65 秒）。
• 约束梯度（dg）的函数评估次数显著减少（45,054 vs. 87,018），表明效率得到提升。
• 在曲柄机构模拟中，准正交投影下每步的平均牛顿迭代次数为 0.325，而标准投影为 1.42，表明收敛速度更快。
• 将能量守恒作为不变量引入后，对计算成本影响极小，但显著改善了能量保持性，如图 6.8 所示。
• 在曲柄机构案例中，使用准正交投影使总模拟时间减少 90%以上（从 3292.00 秒降至 848.67 秒），展现出显著的性能提升。