SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK) C INCLUDE 'ABA_PARAM.INC' C DIMENSION ULOCAL(NDIM),JELEMLIST(100) DIMENSION ALOCAL(NDIM,*),TIME(1) DIMENSION JMATYP(*),JGVBLOCK(*) CHARACTER*80 PARTNAME DIMENSION ARRAY(3) DIMENSION VEL(6) CPRESS=0.0D0 VELOCITY=0.0D0 if (NODE.EQ.3)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.6)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.9)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.12)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.16)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.18)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.21)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.24)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.27)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.30)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.33)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.36)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.39)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.42)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.45)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.49)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.51)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.54)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.57)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else if (NODE.EQ.60)then ULOCAL(NDIM)=0 call GETVRN(NODE,'V',VEL,JRCD,JGVBLOCK,LTRN) else call GETPARTINFO(NODE,0,PARTNAME,LOCNUM,JRCD) call GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,JRCD,JGVBLOCK) call GETVRMAVGATNODE(NODE,'CSTRESS',ARRAY,JRCD,JELEMLIST,NELEMS,JMATYP,JGVBLOCK) end if CPRESS = ARRAY(1) VELOCITY = VEL(1) WEARRATE=0.000035377*VELOCITY*CPRESS ULOCAL(NDIM)=ULOCAL(NDIM)-WEARRATE c write(7,*),WEARRATE RETURN END