4. 接触属性与摩擦模型
在 Abaqus 的接触分析中，法向行为决定“能不能压进去”，切向行为决定“会不会滑、滑多远、滑得顺不顺”。切向接触属性一旦设置不当，轻则收敛困难，重则得到完全错误的装配载荷路径。本节围绕库仑摩擦、各向异性摩擦、摩擦衰减、粘滑行为、摩擦生热以及参数标定展开，重点说明这些模型在 Abaqus 中如何落地，以及工程上如何判断参数是否合理。
一、库仑摩擦：切向接触的基本模型
Abaqus/Standard 默认的切向接触是“无摩擦”，即切向不传递剪应力。要引入摩擦，需要在 Interaction Property 的 Tangential Behavior 中选择 Penalty 或 Lagrange 乘子法，并指定摩擦系数 μ。库仑摩擦的核心判据是：当等效剪应力 τ_eq 小于 μp（p 为接触压力）时处于粘着状态；达到临界值后开始滑动。Penalty 法允许弹性滑移，滑移量由罚刚度 k_s 控制，k_s 越大越接近理想粘着，但收敛越困难。
工程上常用“弹性滑移容差”来间接控制罚刚度。例如在 Abaqus 中设置 Slip Tolerance 为特征单元尺寸的 0.5%，意味着当切向位移达到该量级时剪应力接近 μp。对于一般结构钢-钢干摩擦，μ 可取 0.10～0.15；带氧化皮或加工表面可取 0.2～0.3；润滑表面可低至 0.05。这些数值必须结合表面状态说明，不能直接当作材料常数。
二、各向异性摩擦：方向相关的切向阻力
当接触面存在明显纹理、拉拔方向或纤维铺层方向时，各方向摩擦阻力并不相同。Abaqus 通过 Anisotropic Friction 定义两个正交方向的摩擦系数 μ1、μ2，并指定一个局部方向（通常用局部坐标系或用户子程序定义）。其滑动判据可写为：
f = sqrt( (τ1/μ1)² + (τ2/μ2)² ) - p ≤ 0当 f < 0 时粘着，f = 0 时沿椭圆势面方向滑动。典型场景包括：复合材料层合板沿纤维方向与垂直纤维方向摩擦不同；冷轧钢板沿轧制方向与横向摩擦不同；橡胶密封件沿唇口方向与轴向摩擦不同。设置时需注意局部方向必须随接触面转动而更新，否则大滑移后会失去物理意义。
三、摩擦衰减与粘滑行为
摩擦系数并非恒定。摩擦衰减描述的是 μ 随滑移速度或滑移距离增大而下降的现象。Abaqus 中可通过 Friction Coefficient 的 Velocity 依赖或 Field 依赖来定义 μ(v) 曲线。例如制动器摩擦片在低速时 μ 可达 0.4，高速时降至 0.25，这种负斜率会诱发系统自激振动。
粘滑行为则是粘着与滑动交替出现的状态。在 Abaqus/Standard 中，粘滑通常表现为收敛困难、切向力锯齿状波动；在 Abaqus/Explicit 中则表现为接触力高频振荡。判断粘滑是否发生，可观察切向力-位移曲线是否出现“上升—突降—再上升”的循环。若需要稳定求解，可适当降低罚刚度、引入阻尼或改用显式动力学。
四、摩擦生热基础
当接触面发生相对滑动时，摩擦功转化为热量。Abaqus 中通过 Interaction Property 的 Thermal Conductance 与 Heat Generation 选项启用摩擦生热。关键参数是生热分数 η：η 表示摩擦耗散功率中转化为热量的比例，通常取 0.9～1.0；剩余部分以声能、材料损伤等形式耗散。热流在上、下表面的分配由材料热导率与热容决定，Abaqus 默认按接触对两侧的热材料属性自动分配。
摩擦生热分析必须使用耦合温度-位移分析步（Coupled Temperature-Displacement），并定义材料的热传导率、比热容和密度。若忽略温度对 μ 的影响，可在高温工况下产生显著偏差。例如制动盘表面温度从 20°C 升至 300°C 时，μ 可能下降 20%～30%，直接改变制动力矩分配。
五、工程参数标定方法
摩擦参数不能仅凭手册选取，应通过试验标定。常用方法包括：
- 斜面法：将试件置于可调倾角斜面，记录开始滑动时的临界角 θ，则 μ = tanθ。适用于低速、干摩擦的初步估计。
- 销-盘试验：在恒定法向载荷下测量切向力，得到 μ-速度曲线，用于拟合摩擦衰减。
- 环-块试验：更接近面接触工况，可同时测量温升与摩擦系数变化。
- 反演标定：用试验测得的力-位移曲线，在 Abaqus 中调整罚刚度与滑移容差，使仿真曲线与试验曲线吻合。
标定顺序建议为：先确定室温、低速下的基准 μ₀；再通过变速度试验拟合 μ(v)；若有温度影响，补充变温度试验拟合 μ(T)；最后用耦合热分析验证温升与摩擦系数的耦合响应。标定完成后，应在 Abaqus 中检查切向力-位移曲线是否出现异常振荡，并确认能量平衡中摩擦耗散与生热总量一致。
综上，接触切向行为的设置是一条从“界面受力”到“参数选择”再到“状态判据”和“热-力反馈”的完整链条。库仑摩擦提供基础框架，各向异性摩擦处理方向差异，摩擦衰减与粘滑行为描述动态特性，摩擦生热则将力学与热学耦合起来。工程参数标定必须依托试验，并在 Abaqus 中通过力-位移曲线、能量平衡和温度场进行交叉验证，才能保证复杂装配结构仿真结果的可靠性。