磨耗模拟
磨耗理论基础
磨耗模拟是什么
老师，用FEM可以模拟磨耗（磨损）吗？
可以。基于Archard磨耗规律计算接触面的材料去除，并更新形状。
$V$ 是磨耗体积，$s$ 是滑动距离，$K$ 是磨耗系数，$F_n$ 是法向力，$H$ 是硬度。
按滑动距离成比例磨耗。这是一个很简单的模型呢。
Archard公式是最简单的，但对许多磨耗问题能得到实用精度。
FEM中的实现
磨耗模拟的步骤（迭代式）：
1. 接触分析 — 计算接触压和滑动量
2. 磨耗量计算 — 用Archard公式计算磨耗深度
3. 网格更新 — 沿磨耗量移动节点
4. 重新接触分析 — 用更新的形状
5. 迭代至收敛
每次都更新网格似乎很困难呢。
Abaqus的UMESHMOTION（自适应网格）或LS-DYNA的*MAT_WEAR可自动化。
总结
要点：
- Archard磨耗规律 — $dV/ds = K F_n / H$
- 接触→磨耗量→网格更新的迭代 — 磨耗时间演化
- Abaqus的UMESHMOTION — 自适应网格更新磨耗面
- 轴承、齿轮、制动衬垫 — 主要应用
Archard磨耗规律1953年
磨耗模拟基础的Archard磨耗规律源自J.F. Archard和W. Hirst在1953年发表的论文。磨耗量V = k × W × s / H（k:磨耗系数、W:荷载、s:滑动距离、H:硬度）虽然是简单比例规律，但60多年后的今天仍被工业磨耗预测作为标准模型使用。该论文先于Wear期刊创刊（1957年），对摩擦学（tribology）学术分支确立做出了巨大贡献。
磨耗数值计算方法
磨耗的FEM实现
```
*ADAPTIVE MESH, TYPE=CONTACT SURFACE
wear_surface
*ADAPTIVE MESH CONTROLS, MESHING FREQUENCY=1
*UMESHMOTION $ 用户子程序计算磨耗量
```
```
*MAT_WEAR
K, H $ Archard系数和硬度
```
LS-DYNA有专用的磨耗材料模型呢。
*MAT_WEAR自动将接触面的单元变薄来表示磨耗。不需要重新网格化，磨耗会自动进行。
总结
形状更新的数值方案
磨耗FEM的数值核心在于"磨耗量→节点移动→网格更新→再接触"的迭代处理。Podra & Andersson（1997年）提出的欧拉方法（Eulerian approach）中，计算一个周期的磨耗量用Archard规律，之后沿法向后退边界节点，再用平滑算法修正单元变形。该手法至今仍是ABAQUS UMESHMOTION用户子程序实现的参考标准。
磨耗实务应用
磨耗的实务
实务检查清单
切削刀具磨耗预测
从2000年代起，将Archard磨耗规律融入FEM的金属切削刀具寿命预测逐渐普及。山特维克在2010年代公开的硬质合金立铣刀分析中，使用DEFORM代码在切削速度200m/min、进给量1mm条件下，成功以±18%的误差预测了工具侧面磨耗宽度VB=0.3mm时的加工长度与实验值的对比。该结果被用于刀具涂层优化的设计判断。
磨耗软件对比
磨耗工具
选型指南
DEFORM磨耗模块的历史
金属加工专用FEM代码DEFORM（Scientific Forming Technologies Corporation）从1989年建立以来一直专注切削、锻造分析，磨耗模块的商用提供始于DEFORM V6.1（2000年左右）。该版本实现了Archard规律和Usui规律（切削用）两种，工具寿命预测精度成为业界标准。现行DEFORM V13中增加了机器学习自动补正磨耗系数的自适应功能。
磨耗先端研究
磨耗的先端研究
数据驱动磨耗规律
2019年以来，将分子动力学模拟与PINNs（物理信息神经网络）结合，将纳米尺度磨耗机制上标度到宏观磨耗规律的研究活跃化。2023年IFToMM World Congress的报告中，由Fe-Fe接触MD模拟10万个案例学习的PINNs成功以15%以内的预测误差捕捉到Archard模型无法再现的低荷载域非线性磨耗行为。
磨耗故障对应
磨耗的故障
磨耗系数的波动与失败
磨耗模拟中最大的不确定因素是磨耗系数k的值。即使是相同钢材，由于润滑状态、表面粗糙度、温度的不同，k可变动3个数量级从10⁻⁶到10⁻³。2015年国内汽车零部件制造商经历的失败案例中，直接使用文献值k=5×10⁻⁵时，预测磨耗量竟是实测的7倍，导致设计判断失误。根本对策是用Pin-on-Disk摩擦计进行使用条件对应的k实测。
有用
详情
错误