ABAQUS----UMESHMOTION子程序，磨损分析

3月前浏览18326

本文摘要（由AI生成）：

文章主要介绍了UMESHMOTION子程序在磨损和烧蚀方面的应用，并提供了一些相关资料和案例。UMESHMOTION子程序可用于实现磨损和烧蚀的仿真，但网上关于该子程序的资料较少。文章作者整理了一些相关资料，包括ABAQUS用户手册、子程序讲义、磨损书籍、外文资料、百度文库中的例子以及一个简单的材料烧蚀帖子。这些资料可以帮助读者更好地理解和应用UMESHMOTION子程序。

UMESHMOTION子程序可用于磨损和烧蚀，网上关于UMESHMOTION子程序的资料较少，笔者找了一些资料整理如下。由于本人并非专业人士，无法对质量进行过多评价，希望能对需要的人员有所帮助：

1. ABAQUS用户手册——轮胎磨损经典例子，考虑轮胎纹理边角的磨损，涉及到不同磨损区域的判定和局部磨损方向的定义。（本人对其进行了逐行的详细解释）

http://wufengyun.com:888/books/exa/default.htm

2. ABAQUS子程序讲义，内含该程序的基本概念解释。

https://fdocuments.in/document/user-sub.html

3. 一本介绍磨损的书籍，内含一些ABAQUS磨损案例，可通过这些经典案例的曲线形状判断子程序结果准确与否。

书名：Numerical Modeling of the Effect of Fretting Wear on Fretting Fatigue

http://eprints.nottingham.ac.uk/10681/1/ThesisFinal.pdf

4. 一个外文资料，基于Archard理论和实验进行对比，并对相关理论公式进行说明。

https://www.scirp.org/html/1-2190164_96936.htm#f3

5. 百度文库中的一个例子，简单二维磨损案例有配图，没有基于磨损理论进行子程序编写，也没有说明磨损过程。

https://wenku.baidu.com/view/305572e86aec0975f46527d3240c844769eaa0c4.html

6. 一个简单的材料烧蚀帖子及启发（附案例）。

https://www.smart3de.com/forum.php?mod=viewthread&tid=1270&extra=page%3D2

7. 李玲等考虑到Archard模型的局限性，基于能量模型：《不同加载条件下柱面 /平面微动磨损有限元分析》,及其引用文章。（康乐学位论文：《栓接结合部的微动磨损动态仿真研究》）

文章检索号： DOI: 10．13433/ j．cnki．1003-8728．20180132

8. 实现轴孔过盈配合磨损的相关描述:《微动磨损对过盈配合结构微动参量的影响》,及其引用文章。

文章检索号： DOI: 10．16078/ j．tribology．2015．04.018

知网上还有很多UMESHMOTION子程序的文章和涉及的硕士/博士论文，这里就不一一列举了。

————————————————————————————————————————————————

UMESHMOTION子程序接口的相关解释在上述资料中，该子程序的核心是定义：化学反应速率/磨损速率——通过ULOCAL进行传递，同时需要一些功能函数辅助，摘引介绍如下图：

1 - 副本.png

根据实验标定、编写程序并调试其相关参数（k值）是一件十分繁琐的事，笔者并未结合磨损加载实际工况和参数选取，仅参照参考文献4的相关理论，完成了一个简单的三维模型的磨损实现案例，磨损处加密网格，目的是实现子程序和ALE自适应网格的仿真结合，仅做为案例操作演示CAE操作用。

1 00_00_00-00_00_30~7.gif

付费内容为上述动图案例的CAE(2020版本)、.inp、.for文件。以及本人针对帮助文档中轮胎磨损子程序逐字句的中文解释（解释+程序337行），读懂这个最复杂的简单的平面磨损也就迎刃而解了。ABAQUS帮助文档中的磨损本构公式

子程序截取片段如下：

      NELEMS = NELEMMAX
c     前文定义(NELEMMAX=100)
      CALL GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,
     $     JRCD,JGVBLOCK)
      luseendisp = 1
c     KINC：增量步数量 kmeshsweep：网格扫略数量 lvalidinc：初始赋值标记（初值为-1）
      IF (KINC.EQ.1.AND.kmeshsweep.eq.0.and.lvalidinc.eq.-1) then
C       first time in this routine altogether
        lvalidinc = 2
        lvalidsweep = 1
        luseendisp = 0
      else if (kinc .eq. lvalidinc) then
C       first time in this routine this new increment/ 第一次循环
C       Move back copy of energy to front copy