Télécharger eurodi.eso

Retour à la liste

Numérotation des lignes :

  1. C EURODI SOURCE CB215821 16/04/21 21:16:42 8920
  2. SUBROUTINE EURODI(HO,ROH,FCM,S,TMAX,MAILLE,TEMPS
  3. 1 ,TREL,EVO,YOUNG,JDIM)
  4.  
  5. *
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. -INC CCOPTIO
  9. *
  10. REAL*8 TEMPS(0:1000),REL(0:1000),YOUNG(0:1000)
  11. REAL*8 TEMP0(0:100),TR1(0:4),GAMMA(0:4),TR(0:8)
  12. REAL*8 MODULE(0:10),TABMOD(0:100,0:10)
  13. REAL*8 TABGAM(0:10,0:4),TREL(0:10)
  14. REAL*8 EVO(0:1000,0:10)
  15. *
  16. TINI=1.
  17. TFIN=TMAX
  18. DO 10 I=0,100
  19. DO 11 J=0,10
  20. TABMOD(I,J)=0.
  21. 11 CONTINUE
  22. 10 CONTINUE
  23. DO 12 I=0,10
  24. DO 13 J=0,4
  25. TABGAM(I,J)=0.
  26. 13 CONTINUE
  27. 12 CONTINUE
  28. *
  29. * liste des t0 a calculer
  30. *
  31. TEMP0(0)=TINI
  32. TMP=LOG((TFIN-TINI)/.25)/.5
  33. NMAX=INT(TMP)
  34. IF (NMAX.GT.100) THEN
  35. NMAX=100
  36. ENDIF
  37. DO 20 I=1,NMAX
  38. TEMP0(I)=.25*EXP(.5*FLOAT(I))+TINI
  39. 20 CONTINUE
  40. *
  41. * boucle sur les t0
  42. *
  43. DO 30 L=0,NMAX
  44. T0=TEMP0(L)
  45. TINI=TEMP0(L)
  46. CALL TABTMP(HO,ROH,FCM,S,T0,TINI,TFIN,TEMPS,JDIM)
  47. CALL FLUREL(HO,ROH,FCM,S,T0,TEMPS,JDIM,REL)
  48. *
  49. CALL CALMOD(TEMPS,JDIM,REL,MAILLE,MODULE,TR)
  50. IF(IERR.NE.0) RETURN
  51. *
  52. DO 31 I=0,MAILLE
  53. TABMOD(L,I)=MODULE(I)
  54. TREL(I)=TR(I)
  55. 31 CONTINUE
  56. 30 CONTINUE
  57. *
  58. DO 40 INDICE=0,MAILLE
  59. CALL CALEVO(TEMP0,NMAX,TABMOD,INDICE,GAMMA,TR1)
  60. IF(IERR.NE.0) RETURN
  61. DO 41 J=0,4
  62. TABGAM(INDICE,J)=GAMMA(J)
  63. 41 CONTINUE
  64. 40 CONTINUE
  65. *
  66. * on recupere la liste des temps dans TEMPS
  67. *
  68. T0=TEMP0(0)
  69. TINI=TEMP0(0)
  70. CALL TABTMP(HO,ROH,FCM,S,T0,TINI,TFIN,TEMPS,JDIM)
  71. *
  72. * on recupere la liste des gamma(t)
  73. *
  74. DO 50 I=0,MAILLE
  75. EVO(0,I)=0.
  76. DO 51 J=0,JDIM
  77. EVO(J,I)=TABGAM(I,0)
  78. 1 +TABGAM(I,1)*EXP(-TR1(1)*TEMPS(J))
  79. 2 +TABGAM(I,2)*EXP(-TR1(2)*TEMPS(J))
  80. 3 +TABGAM(I,3)*EXP(-TR1(3)*TEMPS(J))
  81. 4 +TABGAM(I,4)*EXP(-TR1(4)*TEMPS(J))
  82.  
  83. 51 CONTINUE
  84. 50 CONTINUE
  85. *
  86. * on recupere le module d'young
  87. *
  88. YOUNG(0)=0.
  89. DO 60 J=0,JDIM
  90. YOUNG(J)=0.
  91. DO 61 I=0,MAILLE
  92. YOUNG(J)=YOUNG(J)+EVO(J,I)
  93. 61 CONTINUE
  94. 60 CONTINUE
  95.  
  96.  
  97. END
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  

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