Télécharger la plaquette Cast3M


fr en

Télécharger evolth.eso

Retour à la liste

Numérotation des lignes : 

evolth

  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
  89.        CALL INTERE(TF,TE,E,EV1)
  90.        CALL INTERE(TF,TE,FC,EV2)
  91.        CALL INTERE(TF,TE,FT,EV3)
  92.        CALL INTERE(TF,TE,GFTD,EV4)
  93. C
  94.        YOUN= YOUNT*EV1
  95.        RB  = RBT*EV2
  96.        ALFA= ALFAT*(EV3/EV2)
  97.        GFT = GFTT*EV4
  98.        GFC = GFCT*EV4
  99. C       ALPH= ALPH
  100.       ELSE
  101.        RB  = RBT
  102.        ALFA= ALFAT
  103.        YOUN= YOUNT
  104.        GFT = GFTT
  105.        GFC = GFCT
  106. C       ALPH= ALPH
  107.       ENDIF
  108. C
  109. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  110. C
  111.         RETURN
  112.         END
  113.  
  114.  
  115.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales