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