C	Creep-Subroutine nach Norton Power-Law und alle drei 
C	Kriechgebiete

SUBROUTINE CREEP(DECRA,DESWA,STATEV,SERD,EC,ESW,P,QTILD, TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND, COORDS,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) 

C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
C
	DIMENSION DECRA(5),DESWA(5),STATEV(*),PREDEF(*),DPRED(*),TIME(2),1COORDS(*),EC(2),ESW(2)
C

C Parameteridentifikation wichtig, Temperaturabhängigkeit 
C Interpolation!
	
	A=2.1E-005
	An=3.3
	Am=0
	

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	

	LEXIMP=1 
	LEND=1

C deviatorische Vergleichsinkrementsdehnung

	DECRA(1)=A*(sig**An)*(time(2)**Am)*dtime

C Ableitung der deviatorischen Vergleichsinkrementsdehnung 
C nach der von-Mises-Vergleichspannung, notwendig wenn 
C implizit gerechnet wird
	
	IF(LEXIMP.EQ.1) THEN

	DECRA(5)=A*An*(sig**(An-1))*(time(2)**Am)*dtime

	END IF

      RETURN
      END