4. 接触力学建模
在灵巧手抓取规划中，接触力学模型是连接“几何位姿”与“力/力矩可控性”的桥梁。手指与物体之间的接触并非抽象约束，而是由材料、曲率、法向载荷与摩擦条件共同决定的物理界面。若接触模型选错，后续的力分配、抓取稳定性判据与滑移预测都会出现系统性偏差。本节围绕点接触、线接触、面接触、软指接触、赫兹接触以及接触力旋量表示展开，建立从微观变形到宏观力旋量的完整认知链。
一、接触类型与运动约束
接触类型决定了接触点处允许的相对运动自由度，也决定了可传递的力与力矩分量。最基础的是无摩擦点接触：接触面被视为一个理想点，只能传递沿法向的推力，切向力为零。若加入摩擦，则成为摩擦点接触，切向可传递不超过 μN 的摩擦力，但仍不能传递绕法向的力矩。当接触区域退化为一条线时，称为线接触，常见于圆柱指面与棱边或圆柱面之间的配合，其法向分布沿接触线变化，可传递绕接触线方向的力矩。面接触则具有有限接触面积，能够传递法向力、切向力以及绕法向的力矩，是包覆抓取和软指抓取中常见的接触形式。
软指接触是面接触的重要特例。由于指面材料（如硅胶、橡胶）弹性模量低，接触区域会随法向力增大而扩展，形成“面积可变”的接触。这种接触不仅传递力旋量，还能通过变形适应物体表面，提高抓取鲁棒性。在建模时，软指接触常被简化为多个点接触的集合，或采用弹性半空间模型描述压力分布。
二、赫兹接触：从曲率到压力分布
赫兹接触理论描述了两个弹性体在法向载荷作用下，接触区域内的压力分布与变形关系。对于两个球体或球-平面接触，接触区为圆形，压力呈半椭球分布，中心压力最大，边缘为零。其核心参数包括等效曲率半径 R、等效弹性模量 E*、法向载荷 N 与接触半径 a。关系可写为：
a = (3NR / 4E*)^(1/3)
p(r) = p0 * sqrt(1 - (r/a)^2)
p0 = 3N / (2πa^2)
δ = a^2 / R其中 δ 为法向趋近量，p0 为最大接触压力。赫兹模型的价值在于：它把“法向力”与“接触面积、压力峰值、变形量”定量联系起来。在灵巧手中，当指尖为球形且物体表面曲率已知时，可用赫兹接触估算接触半径，进而判断是否进入面接触状态。若接触半径远小于指面尺寸，点接触近似成立；若接触半径接近指面宽度，则必须按面接触或软指接触处理。
三、接触力旋量表示
接触力旋量是抓取力分析的统一语言。对于第 i 个接触点，设接触坐标系原点位于接触点，法向为 n，切向为 t、o。接触力旋量 w_i 由三维力 f_i 与三维力矩 τ_i 组成，可写为：
w_i = [ f_i ; τ_i ]
f_i = f_n * n + f_t * t + f_o * o
τ_i = τ_n * n + τ_t * t + τ_o * o不同接触类型对应不同的力旋量约束：无摩擦点接触仅 f_n ≥ 0，其余分量为零；摩擦点接触满足 sqrt(f_t^2 + f_o^2) ≤ μ f_n，且 τ = 0；软指接触则允许 τ_n 存在，并常以摩擦锥与力矩上限共同描述。将各接触点的力旋量映射到物体质心，可得到抓取矩阵 G，进而判断抓取是否满足力平衡与力矩平衡。
四、方法流程：从接触到抓取力旋量
在实际工程中，接触力学建模通常遵循以下流程：
- 识别接触几何：根据指尖形状与物体表面曲率，判断是点、线还是面接触。球形指尖与平面物体通常为点接触；圆柱指面与圆柱物体可能形成线接触；软指与任意曲面易形成面接触。
- 选择接触模型：若接触半径远小于指面尺寸，采用点接触加摩擦锥；若接触面积不可忽略，采用面接触或软指接触；若需计算压力分布与变形，引入赫兹接触。
- 建立接触坐标系：在每个接触点建立局部坐标系，确定法向 n 与切向 t、o，并计算接触力旋量分量。
- 构造抓取矩阵：将各接触力旋量映射到物体质心，形成抓取矩阵 G，用于力平衡与稳定性分析。
- 求解与校验：在摩擦锥约束下求解期望接触力，校验是否满足力封闭条件，必要时调整接触位置或指面材料。
五、案例：三指灵巧手抓取圆柱杯
以三指灵巧手抓取圆柱杯为例。指尖为半球形硅胶软指，杯壁为刚性圆柱面。初始接触时，接触区近似为圆形，半径约 2 mm，远小于指面直径 20 mm，可先按赫兹点接触估算。当法向力增大到 5 N 时，接触半径扩展至约 6 mm，此时接触区已不可忽略，需切换为软指面接触模型。每个接触点可传递法向力、切向摩擦力以及绕法向的力矩。将三个接触力旋量映射到杯体质心，构造 6×18 的抓取矩阵，求解满足力平衡的接触力分布。若某指接触力接近摩擦锥边界，则提示滑移风险，需增大法向力或调整抓取位形。
| 接触类型 | 可传递力分量 | 可传递力矩分量 | 典型场景 |
|---|---|---|---|
| 无摩擦点接触 | 法向力 fn | 无 | 理想化分析、形封闭 |
| 摩擦点接触 | fn, ft, fo | 无 | 硬指尖抓取、力封闭 |
| 线接触 | 沿接触线分布 | 绕接触线方向 | 圆柱指面与棱边配合 |
| 面接触 | fn, ft, fo | τn, τt, τo | 包覆抓取、宽指面 |
| 软指接触 | fn, ft, fo | τn 为主 | 硅胶指、易碎物体 |
| 赫兹接触 | 法向压力分布 | 无（理想） | 球形指尖、曲率接触 |
六、常见误区与复盘
误区一：所有接触都按点接触处理。当指面为软材料或法向力较大时，接触面积显著扩展，点接触假设会低估可传递力矩，导致抓取稳定性误判。
误区二：忽略赫兹接触的适用条件。赫兹理论要求小变形、无摩擦、各向同性。若指面有纹理或物体为塑料、橡胶等软材料，直接套用赫兹公式会产生较大误差。
误区三：力旋量只算法向力。切向摩擦力与绕法向力矩同样影响抓取平衡，尤其在面接触和软指接触中，力矩分量不可省略。
复盘时，建议按“接触几何→材料属性→载荷范围→模型选择→实验校验”的顺序逐项检查。对于灵巧手抓取规划，接触力学模型不是越复杂越好，而是要与任务精度、计算实时性和传感器能力匹配。点接触适合快速规划，面接触与软指接触适合精细操作，赫兹接触则用于需要压力分布信息的场景。只有把接触类型、压力分布与力旋量表示统一起来，才能为后续的抓取力优化与滑移控制提供可靠基础。