Télécharger la plaquette Cast3M


fr en

Télécharger bpeldi.eso

Retour à la liste

Numérotation des lignes : 

bpeldi

  1. C BPELDI    SOURCE    CB215821  16/04/21    21:15:23     8920
  2.        SUBROUTINE BPELDI(RM,RH,FC28,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 TABTP1(RM,RH,FC28,T0,TINI,TFIN,TEMPS,JDIM)
  49.         CALL FLURE1(RM,RH,FC28,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.         DO 41 J=0,4
  63.           TABGAM(INDICE,J)=GAMMA(J)
  64. 41      CONTINUE
  65. 40    CONTINUE
  66. *
  67. *     on recupere la liste des temps dans TEMPS
  68. *
  69.       T0=TEMP0(0)
  70.       TINI=TEMP0(0)
  71.       CALL TABTP1(RM,RH,FC28,T0,TINI,TFIN,TEMPS,JDIM)
  72. *
  73. *     on recupere la  liste des gamma(t)
  74. *
  75.       DO 50 I=0,MAILLE
  76.         EVO(0,I)=0.
  77.         DO 51 J=1,JDIM
  78.           EVO(J,I)=TABGAM(I,0)
  79.      1           +TABGAM(I,1)*EXP(-TR1(1)*TEMPS(J))
  80.      2           +TABGAM(I,2)*EXP(-TR1(2)*TEMPS(J))
  81.      3           +TABGAM(I,3)*EXP(-TR1(3)*TEMPS(J))
  82.      4           +TABGAM(I,4)*EXP(-TR1(4)*TEMPS(J))
  83.  
  84. 51      CONTINUE
  85. 50    CONTINUE
  86. *
  87. *     on recupere le module d'young
  88. *
  89.       YOUNG(0)=0.
  90.       DO 60 J=1,JDIM
  91.         YOUNG(J)=0.
  92.         DO 61 I=0,MAILLE
  93.           YOUNG(J)=YOUNG(J)+EVO(J,I)
  94. 61      CONTINUE
  95. 60    CONTINUE
  96.       END
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  

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