C C Isotropic hardening for UHARDMSS.inp C SUBROUTINE UHARD(SYIELD,HARD,EQPLAS,EQPLASRT,TIME,DTIME,TEMP, $ DTEMP,NOEL,NPT,LAYER,KSPT,KSTEP,KINC, $ CMNAME,NSTATV,STATEV,NUMFIELDV, $ PREDEF,DPREDEF,NUMPROPS,PROPS) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 CMNAME C DIMENSION HARD(3),STATEV(NSTATV),TIME(*), $ PREDEF(NUMFIELDV),DPREDEF(*),PROPS(*) C IF (EQPLAS.LE.PROPS(2)) THEN SYIELD = PROPS(1) ELSE IF (EQPLAS.LE.PROPS(4)) THEN SYIELD = PROPS(1) + (PROPS(3)-PROPS(1)) * $ (EQPLAS - PROPS(2)) / (PROPS(4) - PROPS(2)) HARD(1) = (PROPS(3)-PROPS(1)) / (PROPS(4) - PROPS(2)) ELSE IF (EQPLAS.LE.PROPS(6)) THEN SYIELD = PROPS(3) + (PROPS(5)-PROPS(3)) * $ (EQPLAS - PROPS(4)) / (PROPS(6) - PROPS(4)) HARD(1) = (PROPS(5)-PROPS(3)) / (PROPS(6) - PROPS(4)) ELSE SYIELD = PROPS(5) ENDIF STATEV(1) = HARD(1) SUM = PREDEF(NUMFIELDV) + (TIME(1) + TIME(2))/2.0d0 SYIELD = (SUM/20.0d0) * SYIELD SUM = DPREDEF(NUMFIELDV) + DTIME SYIELD = SUM + SYIELD C RETURN END