磨损失效是机械设备最常见的失效模式没有之一。齿轮、轴承、刹车盘、模具、人工关节几乎你能想到的接触类部件最终都逃不过磨损问题。但搞过结构分析的工程师都清楚常规有限元软件里默认是没有磨损这个选项的——Ansys Mechanical能算接触应力、算变形、算疲劳但材料在一点点消失这件事它默认不帮你做。我刚接触这个需求时也头疼了好一阵客户只给了一个接触副模型问模拟跑十万次循环后能磨掉多少软件里找了一圈没找到磨损模块上来就要写APDL命令流。这篇文章就把我在Ansys Mechanical里折腾磨损建模的完整经验写下来包括Archard磨损模型怎么落到软件里、接触计算怎么跟磨损迭代耦合、网格畸变了怎么处理、结果怎么验证希望对做接触磨损仿真的同行有点帮助。1. 磨损仿真不是玄学先把Archard模型吃透1.1 为什么常规有限元跑不出磨损你打开Ansys Mechanical的Contact工具栏能找到摩擦应力、接触压力、滑移距离、间隙、渗透量但找不到磨损深度。原因很简单磨损在物理上是一个渐进的材料迁移过程每一步磨损量极小但累积效果显著而有限元本质是在固定几何上做力与位移求解。几何不变磨损就无从谈起——这是个底层矛盾。所以在Mechanical中对磨损建模这件事本质上不是在软件里找某个按钮而是要自己搭建一个求解→提取接触参量→计算磨损量→更新几何→再求解的闭环。Ansys真正提供的是一个开放的内核接触压力、滑移量、节点坐标底层都能用APDL命令读出来、改回去。这就是它能做磨损仿真的根本原因。1.2 Archard磨损模型的工程化表达磨损模型有很多种Archard模型是工程界用得最广的经典模型原因是它足够简单且参数容易标定。基本表达式为[ V K \frac{F \times S}{H} ]其中V是磨损体积F是法向载荷S是滑移距离H是材料硬度K是磨损系数。但这个式子算的是宏观总量到了有限元里我们要把它局部化。对接触面上的任意一个节点其磨损深度增量可以写成[ \Delta h_i k_w \times p_i \times \Delta s_i ]这里( k_w K / H )即磨损率系数常用单位是m²/N或mm²/N( p_i )是该节点的接触压力( \Delta s_i )是该节点在当前增量步内的滑移距离增量。这个局部化的推导过程不复杂把体积增量分配到单元接触面积上深度增量就等于K乘以压力乘以滑移再除以硬度。关键点在于这个公式默认了磨损只发生在当前接触状态持续期间所以实际求解时每一步的滑移距离增量和接触压力要用同一个荷载步内完成收敛后的结果这一点对接下来的APDL实现至关重要。使用Archard时我建议注意三个坑。第一这个模型不区分磨粒磨损、粘着磨损还是氧化磨损它只给一个总磨损量如果你需要区分磨损机制得单独做机制图判断。第二K值在仿真里通常被当成常数但真实材料在跑合期和稳定磨损期的K值能差一个数量级如果实验数据充足建议做成分段函数。第三单位务必统一工程上经常有人把mm和m混着用最后算出来的磨损深度差三个数量级我后文会给单位自检表。1.3 Ansys Mechanical磨损建模的技术路线对比用户可能搜到过各种磨损仿真的方案这里把主流的几条路线做一个对比方便选择适合自己的技术路线实现方式优点缺点适用场景直接计算法忽略几何更新提取接触压力和滑移量一次性估算总磨损深度速度快、稳定不考虑磨损后接触面积的改变误差大磨损量少于初始尺寸5%的粗估USRCLS/间隙法用单元生死模拟材料去除相对自动大量单元生死导致收敛困难网格依赖性强磨粒磨损凹坑模拟APDL节点偏移法每步求解后按磨损量修改接触面节点坐标精度可控、几何实时更新需要写命令流磨损量过大时网格畸变中低磨损量下的工程模拟网格重构法每N步重新划分网格插值传递状态变量适合大磨损量实现复杂接触状态重映射容易丢深度磨损/磨损量接近尺寸量级我最常用的是第三条APDL节点偏移法。它在Workbench里可以通过Command对象嵌入既能利用Workbench方便的前后处理接口又能直接操作底层的节点和接触结果。后面的内容都是以这条路线为核心展开的。如果你只是粗略估计寿命方案一行也能出个数给领导交差但如果是做减摩结构优化、磨损形貌预测、或者跟实验对标磨痕宽度老老实实走方案三。2. 开工前的关键准备材料参数、接触设置与单位防呆2.1 磨损系数的获取与适用范围很多同学拿到项目第一步就想建模但我强烈建议先花时间确定磨损系数这个系数的正确性决定了整个仿真有没有意义。查文献要注意Archard模型的K值来源五花八门有销盘试验测的、有实际工件反推的、有不同润滑条件下测的随便拿来用很可能跟你的工况相差一个数量级。我自己的处理流程是先检索同材料副、同润滑条件、同载荷量级下的销盘磨损实验数据记录所有文献的K值范围如果条件允许自己做一组简单的销盘实验标定K值哪怕只做两个载荷点也能验证线性假设是否成立实在没条件做实验最终仿真报告里必须把K值作为修正因子讨论并给出不同K值下的磨损结果范围而不是拍脑袋给一个数另外提醒一下磨损系数不是纯材料属性它跟接触面的粗糙度、温度、润滑膜厚度都有关。对标实验数据时确保摩擦副材料、表面处理状态、润滑状态尽量一致。2.2 接触算法选择对磨损计算的影响Ansys Mechanics里接触算法有Pure Penalty纯罚函数、Augmented Lagrange增广拉格朗日、MPC、Normal Lagrange等。磨损仿真的接触设置跟普通强度计算有明显的区别。我的经验是磨损计算优先选Augmented Lagrange不要选Pure Penalty。原因有二一是增广拉格朗日法对接触压力精度更高而磨损深度直接跟接触压力成正比压力差了5%累计几十个循环后磨损量就差得离谱二是它迭代收敛性比Normal Lagrange好磨损计算本身就是多步循环收敛慢会让你痛不欲生。接触刚度设置也要调整。默认的接触刚度为了收敛往往会偏软体现在结果上就是存在较小的接触穿透。普通强度分析那点穿透量无伤大雅但磨损计算里接触刚度偏软会让接触面积偏大中心压力偏小最终磨痕形貌跟实际对不上。建议把Normal Stiffness系数从默认值调高到1.0到3.0之间同时关注接触穿透量是否控制在特征尺寸的1%以内。接触探测点的设置根据CZM、表面投影位置等选择不展开细说但要注意你后续计算磨损深度的节点应该与接触单元的目标面节点一一对应否则提取数据和更新坐标时会有“错位”问题。2.3 一个典型的销-盘磨损模型构建为了方便后文讲解我们用一个经典的销-盘模型作为载体销是Arcas2的圆柱体底面是半径为5 mm的圆面压在旋转的盘上载荷为100 N盘转速设定为使其表面相对滑移速度0.1 m/s。在Workbench里搭建这个模型时注意以下几点几何处理销和盘都用六面体网格接触区域做局部加密接触区网格尺寸建议不大于0.2 mm远离区可以放粗到1 mm。磨痕宽度通常比接触宽度大一些网格太粗会导致磨痕形貌变成锯齿状。接触设置摩擦系数设0.3干摩擦接触公式用Augmented Lagrange接触面选销底面目标面选盘顶面。这里的目标面选择很重要后文节点更新主要作用在目标面上把盘做成目标面更便于模拟盘的磨损取决于你想磨谁。边界条件约束盘底面所有自由度销顶面施加100 N压力盘顶面施加旋转位移或速度。如果是二维轴对称简化模型用Cyclic Symmetry能减少计算量但磨损本身会破坏轴对称性磨损深度随径向分布不均匀所以我的经验是别滥用轴对称简化要做全模型或至少半模型。这个模型我会在后面第3章和第5章反复引用你可以把自己的模型对应到相同概念上。2.4 单位防呆检查表Ansys本身不强制统一单位这就导致磨损仿真特别容易翻车。我列一个自检表每次建模前过一遍物理量SIm-kg-smm-tonne-s常见错误长度mmm把mm当m用力NN质量单位不对压力/应力PaMPaPa和MPa混用磨损系数k_wm²/Nmm²/N数量级错3个量级密度kg/m³tonne/mm³把kg/m³当tonne/mm³输入滑移距离mmm转动位移单位不一致我最惨的一次教训是磨损系数用了mm²/N但几何用的是m结果磨损深度算出来比实际大1000倍当时差点怀疑Archard模型本身的正确性。后来所有参数统一从表格勾一遍这类低级错误再也没犯过。3. APDL命令流驱动磨损迭代核心操作全解3.1 磨损循环的整体程序设计要把Ansys Mechanical变成磨损仿真器需要在求解循环中嵌入一个APDL命令对象。整体逻辑可以用下面这段流程描述第1步施加当前循环的载荷与运动边界求解接触状态第2步进入后处理提取接触面上每个节点的接触压力和滑移增量第3步按局部Archard公式计算每个节点的磨损深度增量第4步按磨损深度增量修改接触面节点的坐标即向内偏移第5步进入下一循环基于更新后的几何再次求解这里有一个重要的设计决策一个载荷步跟多少个磨损循环对应。因为实际工况往往运行几十万次滑动循环如果你一个载荷步对应一次真实循环那要算几十万步完全不现实。实际做法是做一个时间片“压缩”让一个载荷步对应N次真实循环那么这个载荷步内的滑移距离就应该是真实单次滑移距离的N倍。这个方案有效的前提是在一个载荷步时间尺度内接触压力和滑移分布变化不大。磨损是个慢过程这个前提在大多数情况下是成立的。在Ansys里实现循环我用的是Workbench的Command对象串联多个载荷步或者直接写/SOLU下的DO循环。但DO循环在Workbench里要小心两个问题一是GUI交互下循环可能不会自动更新结果文件需要OUTRES,ALL,ALL强制写出二是接触表面节点坐标更新后必须重写接触单元否则接触关系还停留在旧几何上导致下一载荷步接触探测出错。下面对每一步详细展开。3.2 提取接触压力与滑移距离提取接触压力与滑移距离需要进入后处理利用ETABLE命令从接触单元中读取。假设接触单元类型是CONTA174目标单元是TARGE170我们先选中接触单元! 进入后处理 /POST1 ! 读取当前载荷步结果 SET,LAST ! 选择接触单元并定义单元表 ESEL,S,ENAME,,174 ETABLE,CPRESS,CONT,PRES ! 接触压力 ETABLE,SLIDE,CONT,SLID ! 累积滑移距离 ETABLE,STAT,CONT,STAT ! 接触状态在Ansys的接触单元输出中CONT,PRES是节点接触压力CONT,SLID是累计滑移距离。要注意这个滑移距离是从初始状态开始累积的不是当前载荷步的增量。为了得到当前载荷步滑移增量需要记录上一载荷步结束时的滑移量两者相减! 定义磨损增量数组 *DIM,WEAR_DEPTH_INCREMENT,ARRAY,NODE_COUNT,1 ! 遍历接触节点 *DO,I,1,NODE_COUNT NODEI WEAR_NODES(I) PRESI GET_ETABLE_VALUE(NODEI,CPRESS) SLIDE_TOTAL GET_ETABLE_VALUE(NODEI,SLIDE) SLIDE_INCREMENT SLIDE_TOTAL - PREV_SLIDE(NODEI) ! Archard局部模型 WEAR_DEPTH_INCREMENT(I) KW * PRESI * SLIDE_INCREMENT *ENDDO这里说明一下GET_ETABLE_VALUE是自定义宏通过*GET读取节点关联的第一个单元的表值。一个接触节点可能对应多个接触单元建议取平均值或最大值视工程判断而定。我通常对同一节点多个单元的压力和滑移取面积加权平均这样更贴合物理。获取当前滑移值的替代方案是使用CONT,SLID对应的节点分量。如果接触设置允许parting注意滑移量会在脱离接触后又重新累计所以出现负增量时直接归零处理。3.3 节点坐标更新与磨损深度修正拿到每个接触节点的磨损深度增量后更新坐标是至关重要的一步。对于盘上的磨损磨损方向通常认为是目标面法向向内也就是把目标面节点坐标沿接触面法向减去磨损深度! 将磨损深度增量转化为节点坐标偏移量 *DO,I,1,NODE_COUNT NODEI WEAR_NODES(I) NX NX(NODEI) NY NY(NODEI) NZ NZ(NODEI) ! 拿到该节点的接触法向 NORMAL_X GET_NORMAL_COMPONENT(NODEI, 1) NORMAL_Y GET_NORMAL_COMPONENT(NODEI, 2) NORMAL_Z GET_NORMAL_COMPONENT(NODEI, 3) ! 更新节点坐标 NX(NODEI) NX - WEAR_DEPTH_INCREMENT(I) * NORMAL_X NY(NODEI) NY - WEAR_DEPTH_INCREMENT(I) * NORMAL_Y NZ(NODEI) NZ - WEAR_DEPTH_INCREMENT(I) * NORMAL_Z *ENDDO ! 通过UPGEOM将坐标更新到有限元模型中 UPGEOM,FLAG实际用NX()函数是不能直接修改节点坐标的这里我更正一下正确的做法是用NMODIF,NODEI,NX,NZ,NZ。顺便提醒如果节点坐标更新后不更新接触单元下一载荷步求解会把当前接触压力误判为间隙或穿透可能导致接触力突变甚至发散。所以更新节点坐标后必须执行! 重新生成接触对 ESEL,S,ENAME,,174 ESEL,A,ENAME,,170 EALIVE,ALL更稳妥的做法是直接把这一载荷步的接触单元删除然后在下一载荷步开始前重新建立。APDL里对应的是GCGEN或重新定义REAL合集上的接触对。在Workbench环境里如果接触对是通过Automatic Detection自动生成的你更新坐标后建议在Command对象里显式执行ALLSEL ! 基于更新后的几何重新生成目标单元和接触单元 GCGEN,,,TRGT_NAME,CONT_NAMEGCGEN是生成接触对的高级命令用法比较绕不展开。实际项目里我通常是直接修改INP文件在每次循环里重新定义接触对或者干脆在Workbench里把磨损更新放到外部脚本中逐次提交求解。3.4 多磨损周期加速技巧实际工程中不可能一个循环一个循环跑也不建议一个载荷步模拟所有循环我通常的做法是磨损分段加速。把总循环次数划分为若干阶段每个阶段用一个载荷步代表一定次数的循环载荷步内滑移增量乘以循环次数倍数前期磨损变化快、接触面积和压力分布迅速变化建议每个载荷步代表较少的循环次数如1000次一个载荷步后期磨损进入稳定阶段可以用较大的步长比如5000到10000次一个载荷步但是修改滑移增量时有一个矛盾接触计算的滑移量是由几何和运动边界决定的你不能简单地把运动速度放大N倍因为Ansys瞬态求解的动力效应会改变接触状态。更合适的做法是保持运动边界不变但把材料磨损系数乘以循环倍数这也不对因为磨损系数只是材料属性。这里我给出的工程简化方案是如果一个载荷步模拟N次循环那么在第4步坐标更新时磨损深度增量应乘以N而接触压力和单次滑移增量仍保持原来的物理量WEAR_DEPTH_STEP KW * PRESI * SLIDE_INCREMENT * CYCLE_PER_STEP这个方案的关键假设是N次循环的接触压力和每次的滑移距离都近似相同所以总磨损量只是单次磨损量的N倍。这个假设是否成立需要你在结果里观察接触压力分布是否发生显著变化。如果压力变化超过10%就该缩短载荷步对应的循环数。我实际项目里销盘实验跑10万次循环一般分20个载荷步每步5000次循环。拿一个几万节点的模型单核求解每个载荷步约几分钟总共也就一两个小时效率上完全可以接受。4. 收敛控制与网格自适应的工程经验4.1 磨损步进量与子步数的配合磨损计算里最常见的不收敛原因是坐标更新量过大导致接触状态发生突变求解器来回震荡或者直接发散。原因是接触压力对几何形貌极其敏感尤其当接触区域被磨损出一个小凹坑时压力分布可能重新分配。我的经验法则是每个载荷步内的最大磨损深度增量应控制在接触区网格尺寸的10%以内。如果网格是0.2 mm那每一步磨损深度就别超过0.02 mm。如果实际磨损量较大就增加载荷步数而不是增大单步磨损量。同时每个载荷步内部要使用自动时间步长AUTOTS,ON。磨损更新后的第一个子步接触状态经常变化剧烈自动步长可以帮助平滑过渡。我通常会设置子步数为5到20最大子步40最小子步1初始子步大小可以是载荷步时间的1/10。4.2 接触穿透报警的处理用过Ansys做接触的人都知道输出窗口偶尔会出现接触穿透警告但磨损建模时这个警告出现的频率明显更高。原因在于接触几何更新后接触参数尤其是接触刚度没有跟上几何变化导致原本的接触闭合量不再合适。处理方式优先级如下把接触刚度Normal Stiffness再提高一个量级检查目标面和接触面的节点是否完全重合有没有间隙或重叠更新接触单元时检查KEYOPT(5)和KEYOPT(7)设置确认是自动接触闭合还是要求初始闭合量若还不行把接触算法临时改成Pure Penalty收敛后再改回Augmented Lagrange有一个容易忽略的点如果磨损导致接触面出现非常浅的凹陷而接触面间隙值ICONT设置得不够小求解器可能认为目标面上的磨损凹坑只是粗糙度直接把它过滤掉不重新建立接触。所以磨损模型中建议把ICONT设置得非常小比如特征尺寸的0.01%或者设置成基于当前最大间隙的自动值让接触探测不遗漏微小几何变化。4.3 网格畸变的识别与重划分时机磨损量积累后网格畸变几乎是必然的。主要表现为接触面附近的单元被压扁或拉长、雅可比矩阵退化、求解出现负特征值警告。网格畸变的判断可以通过后处理检查单元形态参数Ansys里可以用CHECK命令或者查看SHPP结果。我的经验是当接触区单元的最大长宽比超过10:1时结果精度已经明显下降该重画网格了。Ansys Mechanical本身有网格重划分功能但在磨损过程中自动重划分仍然需要手动介入较多。我的APDL循环里加了一个畸变判定! 检查接触表面单元的JACOBIAN值 ETABLE,JACOBIAN,GEOM,JACOBIAN *GET,JACOMIN,JACOBIAN,1,DEF,MININE *IF,JACOMIN,LT,0.2,THEN ! 触发重新划分 /PREP7 ETDEL,ALL ! 重新生成几何和网格 ! ... *ENDIF实际项目中更靠谱的做法是设定固定的重划分间隔比如每5个载荷步重新划分一次网格并把磨损表面作为新的几何边界从当前磨损深度轮廓重建几何。在Workbench里你可以把这个命令流放在Model树的某个分支中或者写成外部Python脚本驱动Mechanical重划分前处理。但我必须承认Ansys里的网格重划分还不是完全自动化的对复杂模型手工操作量很大。如果磨损量真的很大比如超过初始尺寸的20%我个人建议直接换用专门的磨损分析工具或者考虑做自适应网格加密而不是全局重划分。5. 结果可信度检验从跑出来到敢用5.1 磨损形貌的定性校验软件跑出结果与结论正确是两回事。第一关是磨损形貌的定性检查。以销盘模型为例磨损区域应该是一个环带靠近摩擦副边缘的磨损深度通常大于中心部位——因为销边缘的接触压力尖峰会随着磨损而平缓化。如果你算出来的磨损形状是一个规则的圆盘状凹坑也合理但如果出现锯齿状、间断的磨损带那多半是网格问题或接触探测问题而不是真实的磨损形貌。检验方法很简单在后处理中显示目标面的位移云图磨损量表现为接触面法向向内位移观察磨损区是否连续、光滑、过渡自然。个别节点出现孤立高磨痕或低磨痕通常是压力提取或坐标更新时节点编号错位。5.2 数据后处理与磨损曲线输出定量分析磨损结果要把磨损深度随循环次数变化、磨损体积随时间变化的曲线导出来。APDL后处理里可以这样输出! 输出磨损深度和磨损体积 *CFOPEN,WEAR_RESULT,TXT *VWRITE,CYCLE_NUMBER,WEAR_MAXDEPTH,WEAR_VOLUME (F10.0,3E14.5) *CFCLOS这里磨损体积可以通过对每个接触节点的磨损深度乘以该节点所属的单元面积并求和得到注意节点面积不是单元面积要平均分摊。横轴循环次数、纵轴磨损深度/磨损体积的曲线可以直观看出磨损是否进入稳定阶段。如果曲线一直保持线性增长说明接触压力分布没有发生显著集中或重新分布如果后期斜率变缓说明磨损面积增大、压力降低是典型的跑合期到稳定期过渡。5.3 实战案例销-盘摩擦磨损实验对标我用一个实际的例子来说明结果可信度的验证逻辑。某次项目里客户要求验证某销-盘摩擦副在设计载荷下运行10万次循环后的磨损深度。材料磨损系数来自同工况销盘实验标定数值为2.3×10⁻⁶ mm²/N。仿真设置销半径5 mm盘直径60 mm法向载荷200 N滑移速度0.1 m/s总共分20个载荷步每步对应5000次循环。网格方面接触区尺寸0.15 mm整体约5.4万节点在普通工作站上单核跑了约3小时。仿真结果最大磨损深度0.032 mm磨损体积8.4×10⁻³ mm³磨损带宽度1.8 mm。实验对标相同条件下销盘试验机实际测得的磨痕宽度约2.0 mm最大磨损深度0.029 mm磨损深度相对偏差约10%。考虑到磨损系数本身的标定误差和网格离散误差这个结果完全可以用于工程判断。那一次做下来的体会是Archard模型的框架下仿真精度上限往往由磨损系数的准确度决定而不是由求解器设置决定。如果磨损系数有30%误差你花再多精力调网格也白搭。所以做磨损仿真实验标定K值应该占项目周期的一半以上。6. 踩坑记录与效率优化建议6.1 最常见的三个翻车现场这几年做磨损仿真我见过不少同行踩坑自己也踩过分享三个典型问题翻车一单位没统一导致的天价磨损。有人用mm建模、N作为力单位、而磨损系数却查文献用的m²/N结果每个循环磨损深度算出来是几百毫米整个接触面直接磨穿了。这种问题从物理常识就能发现但很多人不核对就往下走了。所以我在第2.4节的单位自检表建议打印出来贴在工位前。翻车二磨损深度更新在了错误的节点上。有次我提取接触压力时用了接触单元Contact Element上的节点坐标而坐标更新时又用了目标单元Target Element的节点集两套节点在接触面上不完全重合导致盘面出现周期性波动噪声。后来统一改成目标面节点为基准提取压力和滑移都用目标面节点上的接触值更新坐标也用同一组节点问题才消除。翻车三接触状态不收敛但结果不报错。Ansys的接触计算有一个特点接触压力结果即使不够精确也不会显式报错只在输出窗口打印一条警告信息。很多人不看警告就直接取压力值导致磨损量偏差大。我的建议是每个载荷步计算后检查接触状态总和确认穿透量在允许范围内再进入下一步。6.2 让磨损仿真提速的实用配置磨损仿真由于是多载荷步反复求解速度很重要。以下是我验证有效的提速策略使用半模型或四分之一模型在磨损不破坏对称性时有效。但注意一旦磨损形貌偏离对称比如边缘局部磨损这种方法就失效了。接触区网格局部加密、非接触区粗化磨损主要发生在接触带上远离区网格对结果影响很小。用静态/准静态求解如果惯性效应不明显尽量用静态而不是瞬态动力学求解速度快很多还避免动力噪声干扰接触压力结果。并行求解Ansys Mechanical的分布式求解对接触问题的并行加速效果一般但多核跑稀疏矩阵还是有提升的。建议在求解设置里开8到16核。外部脚本驱动批量提交如果循环次数很多把载荷步管理放到Python外部脚本里每次只提交一个载荷步并读取结果、更新几何提交下一步。这比在APDL内部循环灵活得多也方便挂高算集群。我最近做一个模具磨损项目用外部Python脚本驱动Mechanical逐个载荷步提交总循环20步每步约2分钟40分钟就跑完了。如果再配合HPC效率还可以翻倍。6.3 磨损仿真往后的扩展思路写完这篇文章我最想说的是Ansys Mechanical里的磨损建模本质上是用APDL命令流把Archard模型嫁接到接触分析框架上。它不完美网格畸变、自动重划分、收敛跳动都困扰过我们但它在工程上足够用——尤其是做方案对比、参数敏感性分析和磨损形貌预测。如果你后续想往深处走可以关注几个方向一是考虑温度影响的磨损模型结合热-结构耦合分析二是磨损表面粗糙度演化的多尺度模拟那就需要引入统计模型或者离散元了三是跟优化算法结合把磨损分析封装成模块后做耐磨结构的参数优化。我个人的体会是磨损失效分析这个方向越做越发现真实工程问题远比Archard模型复杂很多情况下磨损机制不是单一常数能描述的但作为工程判断工具这套方法已经能帮我们解决大量实际问题。希望大家能在我的经验基础上做出更精细的磨损仿真模型。如果你在操作过程中遇到跟这篇文章里类似的坑欢迎回来交流尤其是网格畸变和接触收敛这两块每个人遇到的问题都值得记录和分享。