My simulation fails with the system error code 5. I appreciate if anyone could help me to solve the issue.
The Wear is simulated in a 3D steel- steel contact pair by UMESHMOTION subroutine using Archard model in Abaqus. The whole model acts as a moving cylindrical pin on a fixed analytical rigid counterpart. ALE adaptive meshing is also used. The motion of nodes is independent of underlying material with the magnitude of 0.0004657 in x direction.
The model follows Penalty friction formulation.
Here is the part of UMESHMOTION subroutine:
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)
DIMENSION VEL(6)
CPRESS = 0.0D0
VELOCITY = 0.0D0
LTRN = 0
CALL GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN)
LOCNUM = 0
PARTNAME = ' '
CALL GETPARTINFO(NODE,0,PARTNAME,LOCNUM,JRCD)
NELEMS = NELEMMAX
CALL GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,
$ JRCD,JGVBLOCK)
JTYP = 0
CALL GETVRMAVGATNODE(NODE,JTYP,'CSTRESS',ARRAY,JRCD,
$ JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
CPRESS = ARRAY(1)
VELOCITY = VEL(1)
WEARRATE = 0.000035377 * VELOCITY * CPRESS
ULOCAL(NDIM) = ULOCAL(NDIM) - WEARRATE
RETURN
END