C EVOLTH    SOURCE    CHAT      05/01/12    23:47:28     5004
C  EVOLTH    SOURCE   NECH   99/22/11
      SUBROUTINE EVOLTH(T0,TF,BETJEF,NECH0,NECH3)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*      COMMON /DBETJEF/AA,BETA,RB0,ALFA0,YOUN0,XNU,GFC0,GFT0,CAR,ETA,
*     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
*     &                ITER,ISIM,IBB,IGAU1,IZON
*      COMMON /DBETNECH0/DT,DC,ALFG,S0,ENDO
*      COMMON /DBETNECH3/RB,ALFA,YOUN,GFC,GFT,ALPH,AT,AC,BC
      SEGMENT BETJEF
       REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
     &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
       INTEGER ICT,ICC,IMOD,IVISS,ITER,
     &                ISIM,IBB,IGAU1,IZON
      ENDSEGMENT
      SEGMENT NECH0
       REAL*8 DT,DC,ALFG,S0,ENDO
      ENDSEGMENT
      SEGMENT NECH3
       REAL*8 RBT,ALFAT,YOUNT,GFCT,GFTT,ALPH
      ENDSEGMENT
C
*--------------------------------------------------------------------*
*                                                                    *
*    CALCUL DE L'EVOLUTION DES CARECTERISTIQUES THERMO-MECANIQUE     *
*                 AVEC LA TEMPERATURE                                *
*                                                                    *
*--------------------------------------------------------------------*
C
      DIMENSION TE(6), E(6), FC(6), FT(6), GFTD(6), FAC(6), FBC(6)
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       1: DEFINITION DU TABLEAU DE TEMPERATURE
C
      TE(1)=20.D0
      TE(2)=150.D0
      TE(3)=350.D0
      TE(4)=450.D0
      TE(5)=550.D0
      TE(6)=750.D0
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       2: DEFINITION DE L'EVOLUTION DU MODULE DE YOUNG
C
      E(1)=1.D0
      E(2)=0.828D0
      E(3)=0.496D0
      E(4)=0.344D0
      E(5)=0.217D0
      E(6)=0.063D0
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       3: DEFINITION DE L'EVOLUTION DE LA RESISTANCE
C                EN COMPRESSION
C
      FC(1)=1.D0
      FC(2)=0.877D0
      FC(3)=0.836D0
      FC(4)=0.771D0
      FC(5)=0.495D0
      FC(6)=0.253D0
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       4: DEFINITION DE L'EVOLUTION DE LA RESISTANCE
C                EN TRACTION
C
      FT(1)=1.D0
      FT(2)=0.885D0
      FT(3)=0.592D0
      FT(4)=0.448D0
      FT(5)=0.322D0
      FT(6)=0.124D0
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       5: DEFINITION DE L'EVOLUTION DE L'ENERGIE DE
C                FISSURATION
C
      GFTD(1)=1.D0
      GFTD(2)=0.862D0
      GFTD(3)=1.207D0
      GFTD(4)=1.086D0
      GFTD(5)=0.879D0
      GFTD(6)=0.793D0
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C       Calcul des differents parametres avec la temperature
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
      DDT=TF-T0
      IF(DDT.GE.0)THEN
       CALL INTERE(TF,TE,E,EV1)
       CALL INTERE(TF,TE,FC,EV2)
       CALL INTERE(TF,TE,FT,EV3)
       CALL INTERE(TF,TE,GFTD,EV4)
C
       YOUN= YOUNT*EV1
       RB  = RBT*EV2
       ALFA= ALFAT*(EV3/EV2)
       GFT = GFTT*EV4
       GFC = GFCT*EV4
C       ALPH= ALPH
      ELSE
       RB  = RBT
       ALFA= ALFAT
       YOUN= YOUNT
       GFT = GFTT
       GFC = GFCT
C       ALPH= ALPH
      ENDIF
C
C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C
        RETURN
        END