Drake项目中立方体水弹性接触不稳定性问题分析

问题现象

在Drake多体动力学仿真系统中，当使用离散接触近似处理立方体与地面接触时，观察到了两种异常现象：

1. 振动现象：立方体会在接触面上产生持续的微小振动，这种现象在使用零接触余量或推荐的1e-4m余量时尤为明显。
2. 倾斜现象：当增大接触余量至1cm时，虽然振动消失，但立方体会逐渐向某个角落倾斜，这与预期行为不符。

有趣的是，当立方体的某个边长被扰动（即使只有1e-4m的微小变化）时，这些问题就会消失。这表明问题可能特别出现在几何形状接近完美立方体的情况下。

问题根源

通过深入分析，我们发现问题的本质在于：

1. 网格粗糙度：当前的水弹性立方体网格为了准确重现压力场而采用了最简化的网格结构，这在离散接触模式下会导致问题。
2. 接触力施加方式：系统在每个多边形（通常是两个大三角形）的几何中心施加单个接触力。在平衡状态下，质心位于中间位置，但任何微小偏离都会导致压力中心前后移动。
3. 显式处理限制：几何形状在SAP（Semi-Analytical Primal）接触解析阶段被"冻结"，导致两个接触点来回振荡。

这相当于将立方体平衡放置在一条对角线的刀刃上，显然是非常不稳定的配置。

解决方案探索

我们尝试了三种不同的解决方案：

1. 压力中心(COP)修正

我们实现了基于张量形式的压力中心修正算法。虽然这种方法使振荡幅度减小、频率增高（因为COP对角变化的反应比质心更快），但由于接触点位置的显式处理，振荡问题仍然存在。

结论：不推荐此方案，因为虽然有所改善，但无法根本解决问题，且会增加代码复杂度和计算成本。

2. 自定义离散接触网格

我们为离散接触设计了一个特殊的立方体网格，将每个立方体面分割为四个三角形。这种方法通过隐式解析水弹性产生的恢复力矩，有效消除了不稳定性。

优势：

• 保持计算效率
• 不增加代码复杂度
• 有效解决问题

效果：仿真显示系统非常稳定，只有初始阶段的短暂瞬态过程。

最终解决方案

基于上述分析，我们确定自定义离散接触网格是最佳解决方案。这种方法在复杂性、计算成本和解决效果之间取得了良好平衡。具体实现将采用每个立方体面分割为四个三角形的网格结构，这足以隐式解析水弹性产生的恢复力矩，从根本上消除不稳定性。

技术意义

这个问题的解决不仅改善了立方体接触仿真的稳定性，也为Drake系统中处理类似几何形状的接触问题提供了参考方案。它展示了在物理仿真中，几何离散化方式对数值稳定性的重要影响，以及如何通过针对性的网格设计来解决特定问题。