本体感知如何在高速跳跃中表征颗粒地形？

仅靠深度的地形模型在冲击时失效

多数本体感知地形方法假设颗粒阻力主要随足部穿透深度增长。腿缓慢下沉时这一假设成立，但跳跃腿是以一定速度触地的。论文的受控跳跃实验显示，由此得到的力—深度曲线会系统性地偏离仅含深度的线性拟合。这种偏离不是随机散布，而是集中在步态的特定时刻。

实验使用一台垂直约束的直驱跳跃机器人，足部为直径 80 mm 的圆盘，在每次试验前用流化床方法制备的 300 µm 玻璃珠床面上跳跃。五个标称触地速度（0.22、0.51、0.83、1.01 和 1.22 m/s）与 2.50、3.75 和 5.00 N/cm 三种压缩刚度交叉组合，每种配置五次试验，共 60 次试验。每次只分析第一个步态，以免先前接触压实表面。足部安装的力传感器提供真值力，电机电流、编码器、IMU 和飞行时间传感器提供本体感知信号。

两个瞬态承担了大部分误差

在一个代表性步态中，足部并未像在刚性地面上那样触地即停；颗粒屈服使它在约 50 ms 内减速。足部 IMU 在冲击时记录到一次急剧向上的加速度瞬态，力传感器测得的力上升到一个较早的峰值，该峰值比刚性地面冲击更低、更宽。在压缩—伸展切换处，控制器提高虚拟腿刚度，又出现了第二次明显的加速度瞬态。

这两个区间——触地后不久和刚度切换之后——正是深度加速度力模型留下大残差的地方。论文把这种位置重合当作诊断线索：两个区间都伴随足部快速减速，指向一种与加速度而非深度或速度相关的力贡献。

被裹挟的颗粒表现为附加质量

论文把地形反作用力拆成三部分：随深度变化的阻力、随速度变化的二次阻力项，以及随加速度变化的附加质量项。附加质量的含义是，高速运动的足部会带动其下方一个锥形区域的颗粒，因此加速足部也同时加速了这部分颗粒质量。在匀速侵入时加速度项为零，这也是仅靠匀速试验无法暴露它的原因。

为检验这一点，作者以 0.022 至 1.1 m/s、步长 0.02 m/s 做了匀速侵入试验，并拟合出一个深度—速度力曲面。加入速度依赖后吻合度有所改善，但两个瞬态残差依然存在。估计出附加质量分布并乘以实测足部加速度后，所得项复现了这些残差的时间位置和形状，只是幅值有适度差异——作者把这一差异归因于梯度积分和整体模型拟合的不确定性。这是支持附加质量为主要贡献者的模型证据，而非对被裹挟颗粒的直接测量。

动量观测器消除惯性偏差

仅由电机电流估计的本体感知力必然包含腿自身的刚体惯性和重力，因此原始准静态估计与力传感器不一致。论文加入一个动量观测器来补偿这些项，从而由机载传感得到修正后的接触力估计。

对地形推断的影响以相对线性执行器真值颗粒阻力 k_GT 的误差来衡量。在固定压缩刚度下扫描触地速度、以及在约 1.0 m/s 触地速度下扫描刚度时，动量观测器修正后的估计都比准静态估计更接近 k_GT，使试验级 RMSE 从 4.84 N/cm 降至 3.47 N/cm。残差仍然存在，尤其是在加速度较大的高触地速度下——这正是下一步要处理的问题。

降低高加速度样本权重缩小差距

由于剩余误差与较大的足部加速度相关，作者改的是回归而不是力模型。每个样本获得一个逆方差权重，一旦足部加速度绝对值超过阈值 a0，sigmoid 函数就提高其被赋予的方差。高加速度样本被视为噪声更大的测量，但仍保留在拟合中，因此回归更依赖步态中较平稳的部分，而不丢弃瞬态。

采用这种加速度感知加权最小二乘后，估计阻力在触地速度扫描和刚度扫描中都更贴近 k_GT，试验级 RMSE 进一步从 3.47 N/cm 降至 2.35 N/cm。两项修正各司其职：动量观测器修正力信号，加权则修正把该信号转成阻力值的方式。论文认为两者都是在不同运动条件下获得一致估计所必需的。

结果适用的范围与边界

证据来自单一平台和单一介质：一台单自由度垂直跳跃机器人，平圆足，在干燥、均匀的玻璃珠上运动，且每次试验只分析第一个步态。作者明确把推广到其他腿形态、多自由度接触以及黏性或流体饱和基底列为未来工作，并指出届时可能需要依赖几何和方向的力表示。报告的 RMSE 是相对线性执行器真值的试验级误差，而非自然地形上的实地性能。

在这一范围内，论文的实践优先级很明确：当触地速度或刚度切换产生较大的足部加速度时，必须显式建模加速度相关的颗粒效应。真正重要的比较是三种处理的排序——准静态 4.84 N/cm、动量观测器修正 3.47 N/cm、动量观测器加加速度感知加权 2.35 N/cm。读者若在决定是否给地形估计器增加复杂度，应先看自己的步态是否会产生论文所指的加速度瞬态；如果会，基于深度的回归就是错误的基线。

本文引用的文献