abaqus磨损分析

417182775

模型可简化为两个方块压在一起的周期性的相互滑动（尺寸不一样），材料可设为钢，压紧力自己给。
目的：1 磨损率；2 预测较小滑块在一定滑动一定周期后的磨损量，参数不全可以自己给主要是作为参考。
目前建模中的难点：将archard损伤模型写入自适应子程序
提示：帮助文档有个轮胎模型的自适应子程序，本模型的子程序算法应该比帮助文档容易的多，因为只有平东，磨损节点一直处于连续接触过程中，没有转动时节点接触状态不连续的情况——可惜本人才，abaqus接触时间短，对fortran更是完全陌生。

关于奖励：除了仿真币外，另外可以追加200——500RMB的报酬。
邮箱yun417182775@yeah.net

417182775

模型可简化为两个方块压在一起的周期性的相互滑动（尺寸不一样），材料可设为钢，压紧力自己给。
目的：1 磨损率；2 预测较小滑块在一定滑动一定周期后的磨损量，参数不全可以自己给，主要是作为参考。
目前建模中的难点：将archard损伤模型写入自适应子程序
提示：帮助文档有个轮胎模型的自适应子程序，本模型的子程序算法应该比帮助文档容易的多，因为只有平动，磨损节点一直处于连续接触过程中，没有转动时节点接触状态不连续的情况——可惜本人不才，abaqus接触时间短，对fortran更是完全陌生。

关于奖励：除了仿真币外，另外可以追加200——500RMB的报酬。

邮箱yun417182775@yeah.net
稍后给出我的INP文件。

417182775

附件INP文件，跪求ABAQUS先知指导

Lai_weiqing

417182775 发表于 2012-3-25 22:31
附件INP文件，跪求ABAQUS先知指导

哥们儿弄出来了吗？我现在也是做这个

llm890727

钢丝刷洗属于这个问题么

good1990

，路过路过

余温花知晓

帮顶。。。。。。。。

menglinishimeng

，路过路过

791698650

高大上啊

ytao123

在线等大神

suyin

路过,楼主成功了吗，可以分享一下程序吗

DWY1

这是Fortran语言的一个磨损的例子，涉及到的archard公式比较简单，我是做模具磨损；你可以参考一下
C     USER INPUT FOR ADAPTIVE MESH CONSTRAINT
C
      SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,
     $     LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,
     $     KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK,LSMOOTH)
C
      include 'ABA_PARAM.INC'

      CHARACTER*80 PARTNAME
      DIMENSION ARRAY(3), JPOS(15)
      DIMENSION ULOCAL(NDIM)
      DIMENSION ALOCAL(NDIM,*), TIME(2)
      DIMENSION JMATYP(*), JGVBLOCK(*)
      PARAMETER (NELEMMAX=100)
      DIMENSION JELEMLIST(NELEMMAX), JELEMTYPE(NELEMMAX)
      PARAMETER (k = 0.00000251)

      CALL GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,
     $     JRCD,JGVBLOCK)     
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc      
      CALL GETVRMAVGATNODE(LOCNUM,JTYP,'CDISP',ARRAY,JRCD,
     $     JELEMLIST,NELEMS,JMATYP,JGVBLOCK)     
      CSLIP = SQRT(ARRAY(2)**2+ARRAY(3)**2)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc        
      CALL GETVRMAVGATNODE(LOCNUM,JTYP,'CSTRESS',ARRAY,JRCD,
     $     JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
      CPRESS = ARRAY(1)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc      
      if (CPRESS.GT.0) then
      ULOCAL(NDIM) = -1 * k * CSLIP * CPRESS
      end if

      RETURN
      END

DWY1

另外，如果往复式（反复）磨损，需要考虑磨损之后的表面节点的变化，这就需要用到python提取表面节点坐标，然后修改模型的inp文件里的节点坐标，重构表面，再进行第二次磨损分析，依次下去。。。我也正在做这个

hqu00

同求啊，有偿的