Esope source files

   1 : C EVOLTH    SOURCE    CHAT      05/01/12    23:47:28     5004           
   2 : C  EVOLTH    SOURCE   NECH   99/22/11
   3 :       SUBROUTINE EVOLTH(T0,TF,BETJEF,NECH0,NECH3)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : *      COMMON /DBETJEF/AA,BETA,RB0,ALFA0,YOUN0,XNU,GFC0,GFT0,CAR,ETA,      
   8 : *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
   9 : *     &                ITER,ISIM,IBB,IGAU1,IZON
  10 : *      COMMON /DBETNECH0/DT,DC,ALFG,S0,ENDO
  11 : *      COMMON /DBETNECH3/RB,ALFA,YOUN,GFC,GFT,ALPH,AT,AC,BC  
  12 :       SEGMENT BETJEF
  13 :        REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  14 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  15 :        INTEGER ICT,ICC,IMOD,IVISS,ITER,
  16 :      &                ISIM,IBB,IGAU1,IZON
  17 :       ENDSEGMENT
  18 :       SEGMENT NECH0
  19 :        REAL*8 DT,DC,ALFG,S0,ENDO
  20 :       ENDSEGMENT
  21 :       SEGMENT NECH3
  22 :        REAL*8 RBT,ALFAT,YOUNT,GFCT,GFTT,ALPH
  23 :       ENDSEGMENT
  24 : C
  25 : *--------------------------------------------------------------------*
  26 : *                                                                    *
  27 : *    CALCUL DE L'EVOLUTION DES CARECTERISTIQUES THERMO-MECANIQUE     * 
  28 : *                 AVEC LA TEMPERATURE                                *
  29 : *                                                                    *
  30 : *--------------------------------------------------------------------*
  31 : C
  32 :       DIMENSION TE(6), E(6), FC(6), FT(6), GFTD(6), FAC(6), FBC(6)
  33 : C
  34 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  35 : C       1: DEFINITION DU TABLEAU DE TEMPERATURE
  36 : C
  37 :       TE(1)=20.D0
  38 :       TE(2)=150.D0
  39 :       TE(3)=350.D0
  40 :       TE(4)=450.D0
  41 :       TE(5)=550.D0
  42 :       TE(6)=750.D0
  43 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  44 : C       2: DEFINITION DE L'EVOLUTION DU MODULE DE YOUNG
  45 : C
  46 :       E(1)=1.D0
  47 :       E(2)=0.828D0
  48 :       E(3)=0.496D0
  49 :       E(4)=0.344D0
  50 :       E(5)=0.217D0
  51 :       E(6)=0.063D0 
  52 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  53 : C       3: DEFINITION DE L'EVOLUTION DE LA RESISTANCE 
  54 : C                EN COMPRESSION
  55 : C
  56 :       FC(1)=1.D0
  57 :       FC(2)=0.877D0
  58 :       FC(3)=0.836D0
  59 :       FC(4)=0.771D0
  60 :       FC(5)=0.495D0
  61 :       FC(6)=0.253D0 
  62 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  63 : C       4: DEFINITION DE L'EVOLUTION DE LA RESISTANCE 
  64 : C                EN TRACTION
  65 : C
  66 :       FT(1)=1.D0
  67 :       FT(2)=0.885D0
  68 :       FT(3)=0.592D0
  69 :       FT(4)=0.448D0
  70 :       FT(5)=0.322D0
  71 :       FT(6)=0.124D0   
  72 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  73 : C       5: DEFINITION DE L'EVOLUTION DE L'ENERGIE DE  
  74 : C                FISSURATION
  75 : C
  76 :       GFTD(1)=1.D0
  77 :       GFTD(2)=0.862D0
  78 :       GFTD(3)=1.207D0
  79 :       GFTD(4)=1.086D0
  80 :       GFTD(5)=0.879D0
  81 :       GFTD(6)=0.793D0  
  82 : C
  83 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  84 : C       Calcul des differents parametres avec la temperature
  85 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
  86 : C                                      
  87 :       DDT=TF-T0
  88 :       IF(DDT.GE.0)THEN
  88 :        CALL INTERE(TF,TE,E,EV1)
  88 :        CALL INTERE(TF,TE,FC,EV2)
  88 :        CALL INTERE(TF,TE,FT,EV3)
  88 :        CALL INTERE(TF,TE,GFTD,EV4)
  89 : C                    
  90 :        YOUN= YOUNT*EV1
  91 :        RB  = RBT*EV2
  92 :        ALFA= ALFAT*(EV3/EV2)
  93 :        GFT = GFTT*EV4
  94 :        GFC = GFCT*EV4
  95 : C       ALPH= ALPH
  96 :       ELSE
  97 :        RB  = RBT
  98 :        ALFA= ALFAT
  99 :        YOUN= YOUNT
 100 :        GFT = GFTT
 101 :        GFC = GFCT
 102 : C       ALPH= ALPH     
 103 :       ENDIF
 104 : C
 105 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 106 : C                       
 107 :         RETURN
 108 :         END 
 109 :  
 110 :