Tek-Tips is the largest IT community on the Internet today!

Members share and learn making Tek-Tips Forums the best source of peer-reviewed technical information on the Internet!

  • Congratulations cowski on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

using UMESH MOTION subroutine

Status
Not open for further replies.

Ajith Kurian Baby

Mechanical
Joined
Mar 28, 2020
Messages
53
Location
IN
I want to use UMESH MOTION subroutine for a reciprocating wear problem. I don't have much knowledge of writing subroutine. I went through Abaqus documentation but was not that helpful for me. I had written a code but its always showing compilation error. kindly help me to solve this issue.
I don't find any tutorial for this subroutine. pls support me

SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,
* LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,
* KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK,LSMOOTH)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION ULOCAL(NDIM),JELEMLIST(*)
DIMENSION ALOCAL(NDIM,*),TIME(2)
DIMENSION JMATYP(*),JGVBLOCK(*)

CHARACTER*80 PARTNAME
DIMENSION ARRAY(3)
DIMENSION VEL(6)

CPRESS=0.0D0
VELOCITY=0.0D0
if (NODE.EQ.8)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)
CPRESS = ARRAY(1)
VELOCITY = VEL(1)
WEARRATE=0.000035377*VELOCITY*CPRESS
ULOCAL(NDIM)=ULOCAL(NDIM)-WEARRATE
RETURN
END
 
Status
Not open for further replies.

Part and Inventory Search

Sponsor

Back
Top