Télécharger la plaquette Cast3M


fr en

Télécharger lcpcdi.eso

Retour à la liste

Numérotation des lignes : 

lcpcdi

  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