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