Télécharger lcpcdi.eso

Retour à la liste

Numérotation des lignes :

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

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