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
C   Parameteridentifikation wichtig, Temperaturabhängigkeit über Interpolation
C 
C	Kriechgleichung: DECRA(1)=(A0*sig**Am0*(time(2)**(An0-1))*An0+A1*sig+A2*sig**
C     Am2*time(2)**(An2-1)*An2)*dtime	
C	Decra(1) = deviatorische Vergleichsdehnungsinkrement
C
C	Parameter bei T=125°C
      A0125=3.0
      A1125=0.0008
      A2125=0.0
      Am0125=2.0
      Am2125=0.0
      An0125=0.3
      An2125=0.0
C
C	Parameter bei T=200°C
      A0200=3.0
      A1200=0.0008
      A2200=0.0
      Am0200=2.0
      Am2200=0.0
      An0200=0.3
      An2200=0.0
C
C	Parameter bei T=300°C
      A0300=3.0
      A1300=0.0008
      A2300=0.0
      Am0300=2.0
      Am2300=0.0
      An0300=0.3
      An2300=0.0
C
C	Parameter bei T=400°C
      A0400=300.0
      A1400=0.0008
      A2400=0.0
      Am0400=2.0
      Am2400=0.0
      An0400=0.3
      An2400=0.0
C
C	Parameter bei T=450°C
      A0450=300.0
      A1450=0.0008
      A2450=0.0
      Am0450=2.0
      Am2450=0.0
      An0450=0.3
      An2450=0.0
C
C	Parameter bei T=500°C
      A0500=300.0
      A1500=0.0008
      A2500=0.0
      Am0500=2.0
      Am2500=0.0
      An0500=0.3
      An2500=0.0
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) auch für das Ende bestimmt wird
C
      LEXIMP=1
      LEND=1
C
C	Rückwirkend im Inkrement geltende Temperatur
C
      TEMPER = TEMP+DTEMP 
C	
C
C  
C	Wenn T kleiner 125°C erfolgt kein Kriechen
C
      IF(TEMPER.LT.125) THEN
      DECRA(1)=0	  
      DECRA(5)=0
      END IF
C	Damit der Rest des Programms nicht ausgeführt wird, ein Sprung ans Ende	  
      IF (TEMPER.LT.125) THEN
      GOTO 100
      ELSE
      CONTINUE
C	EVTL CONTINUE ÜBERFLUESSIG
      END IF
C	  
C
C
C	Wenn T zwischen 125°C und 200°C ist, wird die Dehnung interpoliert:  
C   
      IF (TEMPER.GE.125 .AND. TEMPER.LT.200) THEN
C	Dehnung bei 125°C
      TERM1=(A0125*sig**Am0125*(time(2)**(An0125-1))*An0125+A1125*sig+
     1 A2125*sig**Am2125*time(2)**(An2125-1)*An2125)*dtime
C	Dehnung bei 200°C
      TERM2=(A0200*sig**Am0200*(time(2)**(An0200-1))*An0200+A1200*sig+
     1 A2200*sig**Am2200*time(2)**(An2200-1)*An2)*dtime
C	Interpolation der deviatorischen Vergleichsinkrementdehnung
      DECRA(1)=TERM1+(TERM2-TERM1)/(200-125)*(TEMPER-125)
C
C Ableitung der deviatorischen Vergleichsdehnungsinkrement
C nach der von-Mises-Vergleichspannung, notwendig wenn
C implizit gerechnet wird
C
      IF(LEXIMP.EQ.1) THEN
C   Ableitung der deviatorische Vergleichsinkrementdehnung bei 125°C
      TERM3=(A0125*time(2)**(An0125-1)*An0125*sig**(Am0125-1)*Am0125+
     1 A1125+A2125*time(2)**(An2125-1)*An2125*sig**(Am2125-1)*Am2125)
     2 *dtime	 
C	Ableitung der deviatorische Vergleichsinkrementdehnung bei 200°C
      TERM4=(A0200*time(2)**(An0200-1)*An0200*sig**(Am0200-1)*Am0200+A1200+A2200
     1 A1200+A2200*time(2)**(An2200-1)*An2200*sig**(Am2200-1)*Am2200)
     2 *dtime	 	 
C	Interpolation der Ableitung der deviatorischen Vergleichsinkrementdehnung	 
      DECRA(5)=TERM3+(TERM4-TERM3)/(200-125)*(TEMPER-125) 
      END IF
      END IF  
C
C	Damit der Rest des Programms nicht ausgeführt wird, ein Sprung ans Ende	  
      IF (TEMPER.LT.200) THEN
      GOTO 100
      ELSE
      CONTINUE
C	EVTL CONTINUE ÜBERFLUESSIG
      END IF
C	
C
C
C	Wenn T zwischen 200°C und 300°C ist, wird die Dehnung interpoliert:  
C   
      IF (TEMPER.GE.200 .AND. TEMPER.LT.300) THEN
C	Dehnung bei 200°C
      TERM1=(A0200*sig**Am0200*(time(2)**(An0200-1))*An0200+A1200*sig+
     1 A2200*sig**Am2200*time(2)**(An2200-1)*An2200)*dtime
C	Dehnung bei 300°C
      TERM2=(A0300*sig**Am0300*(time(2)**(An0300-1))*An0300+A1300*sig+
     1 A2300*sig**Am2300*time(2)**(An2300-1)*An2)*dtime
C	Interpolation der deviatorischen Vergleichsinkrementdehnung
      DECRA(1)=TERM1+(TERM2-TERM1)/(300-200)*(TEMPER-200)
C
C Ableitung der deviatorischen Vergleichsinkrementdehnung
C nach der von-Mises-Vergleichspannung, notwendig wenn
C implizit gerechnet wird
C
      IF(LEXIMP.EQ.1) THEN
C   Ableitung der deviatorische Vergleichsinkrementdehnung bei 200°C
      TERM3=(A0200*time(2)**(An0200-1)*An0200*sig**(Am0200-1)*Am0200+
     1 A1200+A2200*time(2)**(An2200-1)*An2200*sig**(Am2200-1)*Am2200)
     2 *dtime	 
C	Ableitung der deviatorische Vergleichsinkrementdehnung bei 300°C
      TERM4=(A0300*time(2)**(An0300-1)*An0300*sig**(Am0300-1)*Am0300+
     1 A1300+A2300*time(2)**(An2300-1)*An2300*sig**(Am2300-1)*Am2300)
     2 *dtime	 	 
C	Interpolation der Ableitung der deviatorischen Vergleichsinkrementdehnung	 
      DECRA(5)=TERM3+(TERM4-TERM3)/(300-200)*(TEMPER-200) 
      END IF
      END IF  
C
C	Damit der Rest des Programms nicht ausgeführt wird, ein Sprung ans Ende	  
      IF (TEMPER.LT.300) THEN
      GOTO 100
      ELSE
      CONTINUE
C	EVTL CONTINUE ÜBERFLUESSIG
      END IF
C	
C
C
C	Wenn T zwischen 300°C und 400°C ist, wird die Dehnung interpoliert:  
C   
      IF (TEMPER.GE.300 .AND. TEMPER.LT.400) THEN
C	Dehnung bei 300°C
      TERM1=(A0300*sig**Am0300*(time(2)**(An0300-1))*An0300+A1300*sig+
     1 A2300*sig**Am2300*time(2)**(An2300-1)*An2300)*dtime
C	Dehnung bei 400°C
      TERM2=(A0400*sig**Am0400*(time(2)**(An0400-1))*An0400+A1400*sig+
     1 A2400*sig**Am2400*time(2)**(An2400-1)*An2)*dtime
C	Interpolation der deviatorischen Vergleichsinkrementdehnung
      DECRA(1)=TERM1+(TERM2-TERM1)/(400-300)*(TEMPER-300)
C
C Ableitung der deviatorischen Vergleichsinkrementdehnung
C nach der von-Mises-Vergleichspannung, notwendig wenn
C implizit gerechnet wird
C
      IF(LEXIMP.EQ.1) THEN
C   Ableitung der deviatorische Vergleichsinkrementdehnung bei 300°C
      TERM3=(A0300*time(2)**(An0300-1)*An0300*sig**(Am0300-1)*Am0300+
     1 A1300+A2300*time(2)**(An2300-1)*An2300*sig**(Am2300-1)*Am2300)
     2 *dtime	 
C	Ableitung der deviatorische Vergleichsinkrementdehnung bei 400°C
      TERM4=(A0400*time(2)**(An0400-1)*An0400*sig**(Am0400-1)*Am0300+
     1 A1400+A2400*time(2)**(An2400-1)*An2400*sig**(Am2400-1)*Am2400)
     2 *dtime	 	 
C	Interpolation der Ableitung der deviatorischen Vergleichsinkrementdehnung	 
      DECRA(5)=TERM3+(TERM4-TERM3)/(400-300)*(TEMPER-300) 
      END IF
      END IF  
C
C	Damit der Rest des Programms nicht ausgeführt wird, ein Sprung ans Ende	  
      IF (TEMPER.LT.400) THEN
      GOTO 100
      ELSE
      CONTINUE
C	EVTL CONTINUE ÜBERFLUESSIG
      END IF
C	
C
C
C	Wenn T zwischen 400°C und 450°C ist, wird die Dehnung interpoliert:  
C   
      IF (TEMPER.GE.400 .AND. TEMPER.LT.450) THEN
C	Dehnung bei 400°C
      TERM1=(A0400*sig**Am0400*(time(2)**(An0400-1))*An0400+A1400*sig+
     1 A2400*sig**Am2400*time(2)**(An2400-1)*An2400)*dtime
C	Dehnung bei 450°C
      TERM2=(A0450*sig**Am0450*(time(2)**(An0450-1))*An0450+A1450*sig+
     1 A2450*sig**Am2450*time(2)**(An2450-1)*An2)*dtime
C	Interpolation der deviatorischen Vergleichsinkrementdehnung
      DECRA(1)=TERM1+(TERM2-TERM1)/(450-400)*(TEMPER-400)
C
C Ableitung der deviatorischen Vergleichsinkrementdehnung
C nach der von-Mises-Vergleichspannung, notwendig wenn
C implizit gerechnet wird
C
      IF(LEXIMP.EQ.1) THEN
C   Ableitung der deviatorische Vergleichsinkrementdehnung bei 400°C
      TERM3=(A0400*time(2)**(An0400-1)*An0400*sig**(Am0400-1)*Am0400+
     1 A1400+A2400*time(2)**(An2400-1)*An2400*sig**(Am2400-1)*Am2400)
     2 *dtime	 
C	Ableitung der deviatorische Vergleichsinkrementdehnung bei 450°C
      TERM4=(A0450*time(2)**(An0450-1)*An0450*sig**(Am0450-1)*Am0450+
     1 A1450+A2450*time(2)**(An2450-1)*An2450*sig**(Am2450-1)*Am2450)
     2 *dtime	 	 
C	Interpolation der Ableitung der deviatorischen Vergleichsinkrementdehnung	 
      DECRA(5)=TERM3+(TERM4-TERM3)/(450-400)*(TEMPER-400) 
      END IF
      END IF  
C
C	Damit der Rest des Programms nicht ausgeführt wird, ein Sprung ans Ende	  
      IF (TEMPER.LT.450) THEN
      GOTO 100
      ELSE
      CONTINUE
C	EVTL CONTINUE ÜBERFLUESSIG
      END IF
C	
C
C
C	Wenn T zwischen 450°C und 500°C ist, wird die Dehnung interpoliert:  
C   
      IF (TEMPER.GE.450 .AND. TEMPER.LT.500) THEN
C	Dehnung bei 400°C
      TERM1=(A0450*sig**Am0450*(time(2)**(An0450-1))*An0450+A1450*sig+
     1 A2450*sig**Am2450*time(2)**(An2450-1)*An2450)*dtime
C	Dehnung bei 500°C
      TERM2=(A0500*sig**Am0500*(time(2)**(An0500-1))*An0500+A1500*sig+
     1 A2500*sig**Am2500*time(2)**(An2500-1)*An2)*dtime
C	Interpolation der deviatorischen Vergleichsinkrementdehnung
      DECRA(1)=TERM1+(TERM2-TERM1)/(500-450)*(TEMPER-450)
C
C Ableitung der deviatorischen Vergleichsinkrementdehnung
C nach der von-Mises-Vergleichspannung, notwendig wenn
C implizit gerechnet wird
C
      IF(LEXIMP.EQ.1) THEN
C   Ableitung der deviatorische Vergleichsinkrementdehnung bei 450°C
      TERM3=(A0450*time(2)**(An0450-1)*An0450*sig**(Am0450-1)*Am0450+
     1 A1450+A2450*time(2)**(An2450-1)*An2450*sig**(Am2450-1)*Am2450)
     2 *dtime	 
C	Ableitung der deviatorische Vergleichsinkrementdehnung bei 500°C
      TERM4=(A0500*time(2)**(An0500-1)*An0500*sig**(Am0500-1)*Am0500+
     1 A1500+A2500*time(2)**(An2500-1)*An2500*sig**(Am2500-1)*Am2500)
     2 *dtime	 	 
C	Interpolation der Ableitung der deviatorischen Vergleichsinkrementdehnung	 
      DECRA(5)=TERM3+(TERM4-TERM3)/(500-450)*(TEMPER-450) 
      END IF
      END IF  
C
C	Damit der Rest des Programms nicht ausgeführt wird, ein Sprung ans Ende	  
      IF (TEMPER.LT.500) THEN
      GOTO 100
      ELSE
      CONTINUE
C	EVTL CONTINUE ÜBERFLUESSIG
      END IF
C
C
C  
C	Wenn T größergleich 500 °C, wird aus Gründen der Extrapolationssicherheit die Kriechdehnung bei 500°C verwendet:
C
      IF (TEMPER.GE.500) THEN
C 
      DECRA(1)=(A0500*sig**Am0500*(time(2)**(An0500-1))*An0500+A1500*sig+
     1 A2500*sig**Am2500*time(2)**(An2500-1)*An2)*dtime
C	  
      IF(LEXIMP.EQ.1) THEN
      DECRA(5)=(A0500*time(2)**(An0500-1)*An0500*sig**(Am0500-1)*Am0500+
     1 A1500+A2500*time(2)**(An2500-1)*An2500*sig**(Am2500-1)*Am2500)
     2 *dtime
      END IF 
      END IF
C
C	  
100    RETURN
      END