C
C     Isotropic hardening for UHARDCTS.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     
      SIG0 = 200.0D0
      Q = 0.0D0
      B = 0.25D0
      CURRENTTEMP = TEMP + DTEMP
      CURRENTFV = PREDEF(1) + DPREDEF(1)
      IF (CURRENTFV.LT.1D-5) THEN
         Q = PROPS(1) - (PROPS(1)-PROPS(2)) *
     $        (CURRENTTEMP / 100.0D0)
      ELSE
         Q = PROPS(2) - (PROPS(2)-PROPS(3)) *
     $        (CURRENTFV / 100.0D0)
      ENDIF
      SYIELD = SIG0 + Q * (1.0D0 - EXP(-B*EQPLAS))
      HARD(1) = B*(Q-(SYIELD-SIG0))
      RETURN
      END