看完就没有不会的磨损仿真：Abaqus磨损从原理到实战指南
1. 磨损建模理论基础
1.1 基本原理
零件在接触过程中，因摩擦、挤压等因素会导致表面材料损耗，即磨损。Abaqus采用改良版的Archard磨损模型进行磨损量计算。其核心特点是：
模拟时仅在接触层面考虑由磨损引起的节点位置偏移。 磨损不影响零件内部单元的计算。 对于中小程度的磨损量，无需启用ALE（任意拉格朗日-欧拉）自适应网格技术即可保证精度。
对于大磨损量（如磨损厚度接近网格单元尺寸）的模拟，可在 Abaqus/Explicit 中结合ALE自适应网格技术，使网格随磨损后的形状动态调整，从而保持计算精度。此外，软件也支持通过特定方法高效模拟如轴承旋转等数百万次循环的磨损过程。
1.2 Archard磨损模型及其形式
Abaqus提供了4种Archard磨损模型，主要分为两大类：
第一类：与摩擦系数相关
适用于磨损主要由滑动摩擦主导的场景（如零件滑动摩擦）。 其物理本质与摩擦消耗的能量相关。 磨损速率与摩擦能密度耗散率成正比。
第二类：与摩擦系数无关
适用于低摩擦或无摩擦场景（如纯挤压导致的磨损）。 磨损主要与接触压力相关。
磨损系数可定义为磨损量、接触压力、温度及场变量的函数，从而实现复杂行为的模拟。例如，可模拟带防护涂层的零件：当涂层被磨掉后，通过设置磨损系数增大来反映基材更易磨损的特性。
Archard磨损方程的不同形式
| 与摩擦系数无关 | ||
| 显式依赖摩擦系数 |
参数说明：
：瞬时磨损速率（单位：长度/时间，如 m/s） ：无量纲磨损系数 ：有量纲集总磨损系数（单位：长度/(应力·时间)，如 m/(Pa·s)） ：摩擦系数 ：法向接触压力（单位：应力，如 Pa） ：滑移速率模（单位：1/时间，如 1/s） ：摩擦应力 ：材料硬度（单位：应力，如 Pa）
当滑移速率不为零且摩擦应力存在时，依赖摩擦系数的磨损速率可写为： 。其中，摩擦能密度耗散率为 ，因此 ，表明磨损速率与摩擦能密度耗散率成正比。
1.3 磨损计算注意事项
在Abaqus/Standard中，仅当接触公式支持接触面局部区域接触时，才会在接触表面累积磨损。因此，通用接触中的补充接触公式（边-面、边-边、点-面）不参与磨损累积。 线性摄动分析步不产生磨损累积。
2. 磨损仿真适用范围与设置
2.1 求解器与接触类型支持
| 支持的接触形式 | ||
| 结合ALE自适应网格 | ||
| 结合分析步循环 | ||
| 分析步支持 | *MATERIAL MIGRATION和稳态传输分析步中的特定情况外，常规分析步均支持 | |
| 单元类型限制 | 仅支持实体单元 |
2.2 磨损属性的定义与分配
定义方法
磨损属性需作为表面属性进行定义。
输入文件用法： *SURFACE PROPERTY, NAME=<表面磨损属性名称>
*WEAR SURFACE PROPERTIES, UNITLESS WEAR COEF=YES/NO, FRIC COEF DEPENDENT=YES/NO, REFERENCE STRESS=<硬度值>
<磨损系数数据行>Abaqus/CAE操作： 进入 相互作用模块→创建相互作用属性→ 类型选择“磨损”。
分配方式
对于接触对：需将命名表面磨损属性分别分配至接触对的两个接触面。
*SURFACE, NAME=<接触面名称>, TYPE=NODE/ELEMENT, PROPERTY=<表面磨损属性名称>
*CONTACT PAIR, INTERACTION=<相互作用属性名称>
<从面名称>,<主面名称>
对于通用接触：通过*SURFACE PROPERTY ASSIGNMENT选项进行分配。
*SURFACE PROPERTY ASSIGNMENT, PROPERTY=WEAR
<接触面A名称>, <磨损属性名称1>
<接触面B名称>, <磨损属性名称2>
注意：如果接触面存在重叠区域，重叠部分将采用最后分配的磨损属性。
3. Abaqus/CAE中的磨损属性设置详解
在Abaqus/CAE的磨损属性设置窗口中，各项参数含义如下：
磨损属性设置主界面
| Name/Type | Wear。 |
| Friction coefficient dependent | 是否让磨损模型依赖摩擦系数 • 勾选：使用“显式依赖摩擦系数”的模型。 • 不勾选（如图）：使用“与摩擦系数无关”的模型。 |
| Unitless wear coefficient | 是否使用无量纲磨损系数 • 勾选：使用无量纲系数 k，需在下方Reference stress中填写材料硬度H。• 不勾选（如图）：使用有量纲集总磨损系数 \hat{k}。 |
| Reference stress | Unitless wear coefficient时激活，用于填写材料硬度H。 |
| Wear Coefficient (Data区域) | k或hat{k}。单位需与模型一致。 |
| Use surface temperature dependent data | |
| Use surface wear distance dependent data | |
| Use contact pressure dependent data | |
| Number of field variables |
3.1 INP文件导出注意事项
从Abaqus/CAE导出INP文件时，磨损属性部分可能出现格式错误，例如**结束符与数据行错位，导致计算报错。错误示例：
*Surface Property, name=WEAR1
*WEAR SURFACE PROPERTIES, FRIC COEF DEPENDENT=NO, UNITLESS WEAR COEF=NO
7.5e-05, , , ,
** <-- 此星号可能与上一行合并
解决方法：手动打开INP文件，确保**独占一行。
4. 高级建模策略：大量循环磨损
4.1 方法一：缩放磨损系数 (Abaqus/Explicit)
对于需要模拟大量物理循环（如数百万次）的场景，可将N个物理循环合并为单个模拟循环。
原理：在磨损系数中引入缩放因子 N。即，在仿真中使用的磨损系数为实际系数乘以N。要求：需合理控制 N的大小，以确保单个模拟循环内局部接触应力分布变化不大，从而保证磨损积分精度。
4.2 方法二：分析步循环 (Abaqus/Standard)
利用Abaqus/Standard的分析步循环功能，重复执行某一分析步或分析步序列，直至满足终止条件（如达到目标磨损量）。
需要在INP文件中手动添加关键字。 关键参数： 首个模拟循环等效的物理循环数。 单个模拟循环允许的最大增量磨损量阈值。 磨损终止条件（如最大总磨损量）。 后续循环的等效循环数会根据前一步的磨损量自动调整。
输入文件示例：
*STEP CYCLING, START, NAME=Sample cycle repeat
10, ! 最大重复次数
10000,25000 ! 初始等效循环数，最大等效循环数
*STEP,NAME=WearStep-1
... ! 分析步定义
*END STEP
*STEP, NAME=WearStep-2
*STEP CYCLING CONTROL, ACTION=END CYCLING, CRITERION=WEAR-BASED, NAME=SCC-1
NodeSet1, AbsMin, 0.004, 0.002 ! 节点集1，判断类型，终止阈值，单循环目标值
NodeSet2, AbsMax, 0.012, 0.006 ! 节点集2，判断类型，终止阈值，单循环目标值
... ! 分析步定义
*END STEP
*STEP CYCLING, END
5. 特殊工况：稳态传输分析
在稳态传输分析步中，材料在固定网格内对流。
磨损计算特点：网格节点处的磨损量，等于某一时刻占据该节点的材料点在旋转坐标系转动一周、沿其流线运动过程中所累积的磨损量。 计算方法：在稳态假设下，材料点的磨损速率历程，通过沿其流线的局部磨损速率空间分布进行平均得到。磨损仅发生在流线与接触面相作用的弧段。 注意：过大的磨损偏移量可能导致接触穿透，影响计算精度。
6. 启用ALE自适应网格 (Abaqus/Explicit)
当预测磨损量与单元尺寸相当时，建议启用ALE自适应网格技术以保持精度。
工作原理： 节点位移主要由磨损量决定（作为空间网格约束）。 在ALE网格调整增量步之间（默认每10个增量步执行一次），磨损量以节点矢量形式参与接触计算。 当触发ALE增量步时，累积的磨损矢量将转化为网格节点的实际位移。
7. 结果输出
磨损仿真结果通过接触磨损变量 CWEAR 输出，它是接触面节点的累积磨损矢量。
Abaqus/Standard：支持场输出和历程输出。 Abaqus/Explicit：仅支持场输出。 注意：输出 CWEAR需要在INP文件中手动请求。
输入文件用法：
请求场输出： *OUTPUT, FIELD
*CONTACT OUTPUT
CWEAR请求历程输出（仅Standard）： *OUTPUT, HISTORY
*CONTACT OUTPUT, NSET=<节点集名称>
CWEAR
