Télécharger evolth.eso

Retour à la liste

Numérotation des lignes :

  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