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