Hi,
Hi,
I urgently need some with my subroutine in USDFLD. I am modelling a stiffness dependency on mean stress (E= f(s11+s22+s33)). I enter in a table in the input file, the values of Young's modulus and the related field variable (pressure). I have tried many things but somehow the subroutine gives incorrect results for the strains. The stresses are evaluated correctly! PLease help.
This is the subroutine.
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME,ORNAMECHARACTER*3 FLGRAY(15)DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),1 T(3,3),TIME(2)DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
C
C ABSOLUTE VALUE OF MEAN PRESSURE
call GETVRM ('SINV',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,MATLAYO,1 LACCFLA)
PRESS=ABS(ARRAY(3))
C USE ABSOLUTE VALUE OF MEAN PRESSURE AS A FIELD VARIABLE
FIELD(1)=PRESS
C
WRITE(6,*) FIELD(1)WRITE(6,*) TIME(2),TIME(1)
C
IF ERROR COMMENT ON .DAT FILEIF (JRCD.NE.0) THENWRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER', NOEL,
1 'INTERGATION POINT NUMBER', NPT
C
ENDIF
RETURN
END
