C	Creep-Subroutine nach Norton Power-Law und alle drei
C	Kriechgebiete
C
      SUBROUTINE CREEP(DECRA,DESWA,STATEV,SERD,EC,ESW,P,QTILD,TEMP,DTEMP,
     1 PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND,COORDS,NSTATV,NOEL,NPT,
     2 LAYER,KSPT,KSTEP,KINC)
C
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
C
      DIMENSION DECRA(5),DESWA(5),STATEV(*),PREDEF(*),DPRED(*),TIME(2),
     1 COORDS(*),EC(2),ESW(2)
C
C Parameteridentifikation wichtig, Temperaturabhängigkeit
C Interpolation!
C
      A   = 3.0
      A1  = 0.0008
      A2  = 0.0
      Am0 = 2.0
      Am2 = 0.0
      An0 = 0.3
      An2 = 0.0
C
C
C
C Definieren der Berechnungsmethode
C LEXIMP=1, bedeutet implizite Kriechintegration
C LEND=1, damit wird die Routine erst am Ende des Inkrements
C aufgerufen.
C Notwendig, da die Kriechdehnung Decra(1) über die Gesamtzeit
C time(2)
C auch für das Ende bestimmt wird
C
      LEXIMP=1
      LEND=1
C
C deviatorische Vergleichsinkrementsdehnung
C
      DECRA(1)=(A0*sig**Am0*(time(2)**(An0-1))*An0+A1*sig+A2*sig**
     Am2*time(2)**(An2-1)*An2)*dtime
C
C Ableitung der deviatorischen Vergleichsinkrementsdehnung
C nach der von-Mises-Vergleichspannung, notwendig wenn
C implizit gerechnet wird
C
      IF(LEXIMP.EQ.1) THEN
C
      DECRA(5)=(A0*time(2)**(An0-1)*An0*sig**(Am0-1)+A1+A2*time(2)**
     3 (An2-1)*An2*sig**(Am2-1))*dtime
C
      END IF
C
      RETURN
      END