Télécharger la plaquette Cast3M


fr en

Télécharger devflo.eso

Retour à la liste

Numérotation des lignes : 

devflo

  1. C DEVFLO    SOURCE    PV        22/04/21    08:32:51     11344          
  2.        SUBROUTINE DEVFLO(INPLAS,SIG,EPSV,VAR,
  3.      * XMAT,NCOMAT,NSTRS,NVARI,EVP1,VARP1,TAU)
  4.       IMPLICIT INTEGER(I-N)
  5.        IMPLICIT REAL*8 (A-H,O-Z)
  6.        DIMENSION SIG(*),VAR(*),XMAT(*),EVP1(*),VARP1(*)
  7.        DIMENSION SS(6),EPS0(6),epsv(8)
  8. C_____________________________________________________________________
  9. C
  10. C STEP-2 : COMPUTATION OF SEQUV
  11. C ------------------------------
  12.        SIGM=0.333333333333333D0*(SIG(1)+SIG(2)+SIG(3))
  13.        DO 40 I=1,NSTRS
  14.        A        =0.0D0
  15.        IF(I.LE.3) A=1.0D0
  16.        SS(I)=SIG(I)-A*SIGM
  17.  40    CONTINUE
  18.        SEQUV = PROCON(SS,SS,NSTRS)
  19.        SEQUV = SQRT(1.5D0*SEQUV)
  20.        EPSE  = VAR(1)
  21.        IF (SEQUV.LT.10.) THEN
  22.           SEQUV = 10.
  23.        ENDIF
  24.        DEPSCP = 0.
  25.        DEPSCS = 0.
  26.        VAR(12) = 0.
  27.  
  28.        IF (INPLAS.EQ.85) THEN
  29. C
  30. C   CCPL
  31. C
  32.          COE1 = XMAT(6)  * (SEQUV ** XMAT(7))  * XMAT(8)
  33.          COE2 = XMAT(9)  * (SEQUV ** XMAT(10)) * XMAT(11)
  34.          COE3 = XMAT(12) * (SEQUV ** XMAT(13)) * XMAT(14)
  35.  
  36.          IF (COE1.GT.COE2) THEN
  37.            COEI = COE1
  38.          ELSE
  39.            COEI = COE2
  40.          ENDIF
  41.          IF (COEI.LT.COE3) THEN
  42.            COEFF = COEI
  43.          ELSE
  44.            COEFF = COE3
  45.          ENDIF
  46.          DEP11 = COEFF * (1 + XMAT(19) * XMAT(18))
  47.          DEP12 = XMAT(15) * SEQUV ** XMAT(16) * XMAT(17) *
  48.      &         XMAT(18)
  49.          DEPSCP = DEP11 + DEP12
  50.  
  51.        ELSE IF (INPLAS.EQ.102) THEN
  52. C
  53. C   SODERBERG
  54. C
  55.          EPSE1 = VAR(2)
  56.          EPSE2 = VAR(3)
  57. C        terme secondaire
  58.  
  59.          COE1 = XMAT(6)  * (SEQUV ** XMAT(7))  * XMAT(8)
  60.          COE2 = XMAT(9)  * (SEQUV ** XMAT(10)) * XMAT(11)
  61.          COE3 = XMAT(12) * (SEQUV ** XMAT(13)) * XMAT(14)
  62.  
  63.          IF (COE1.GT.COE2) THEN
  64.            COEI = COE1
  65.          ELSE
  66.            COEI = COE2
  67.          ENDIF
  68.          IF (COEI.LT.COE3) THEN
  69.            COEFF = COEI
  70.          ELSE
  71.            COEFF = COE3
  72.          ENDIF
  73.          DEP11 = COEFF * (1 + XMAT(19) * XMAT(18))
  74.          DEP12 = XMAT(15) * SEQUV ** XMAT(16) * XMAT(17) *
  75.      &         XMAT(18)
  76.          DEPSCS = DEP11 + DEP12
  77.  
  78. C        terme primaire de type E1 * (1-exp(-Rt))
  79.          E1_AMP = (XMAT(20) * EXP (XMAT(21) * SEQUV)) + XMAT(22)
  80.          RTOT   = (XMAT(23) + (XMAT(24) * SEQUV)) ** XMAT(25)
  81.          DEPSCP = RTOT * (E1_AMP - EPSE1)
  82.  
  83.          VARP1(2) = DEPSCP
  84.          VARP1(3) = DEPSCS
  85.  
  86.        ELSE IF (INPLAS.EQ.86) THEN
  87. C
  88. C   X11
  89. C
  90.          E0  = XMAT(15) * TANH(XMAT(16) * SEQUV)
  91.          VP0 = XMAT(13) * SINH(XMAT(14) * SEQUV)
  92.          VS0 = XMAT(11) * SINH(XMAT(12) * SEQUV)
  93.          VPF = XMAT(6)  * EXP(-XMAT(7) * XMAT(17))+
  94.      &         (1. - XMAT(6)) * EXP(-XMAT(8) * XMAT(17))
  95.          VSF = (1. - XMAT(9)) * EXP(-XMAT(10)*XMAT(17))
  96.      &         + XMAT(9)
  97.          DEPSCP = (VS0 * VSF) + ((VP0 * VPF) - (VS0 * VSF))
  98.      &            * EXP(-EPSE/E0)
  99.  
  100.        ELSE IF (INPLAS.EQ.84) THEN
  101. C
  102. C   COMETE
  103. C
  104.          EPSE1 = VAR(2)
  105.          EPSE2 = VAR(3)
  106.  
  107.          IF (EPSE1.LT.1E-40) THEN
  108.            EPFIN1 = XMAT(6)  * (SEQUV ** XMAT(7)) *
  109.      &              (TAU**XMAT(8))
  110.            COE12  = EPFIN1/TAU
  111.  
  112.            VAR(12) = 2.
  113.          ELSE
  114.            PUI1  = (XMAT(8) - 1.)/XMAT(8)
  115.            COE11 = XMAT(6)  * (SEQUV ** XMAT(7))
  116.            COE12 = XMAT(8)  * (COE11 ** (1/XMAT(8))) *
  117.      &           (EPSE1 ** PUI1)
  118.          ENDIF
  119.  
  120.          IF (EPSE2.LT.1E-40) THEN
  121.            EPFIN2 = XMAT(9)  * (SEQUV ** XMAT(10)) *
  122.      &              (TAU**XMAT(12)) * (XMAT(13) ** XMAT(11))
  123.            COE22  = EPFIN2/TAU
  124.  
  125.            VAR(12) = 3.
  126.          ELSE
  127.            PUI2  = (XMAT(12) - 1.)/XMAT(12)
  128.            COE21 = XMAT(9)  * (SEQUV ** XMAT(10))
  129.            COE22 = XMAT(12) * (COE21 ** (1/XMAT(12))) *
  130.      &           (EPSE2 ** PUI2) * (XMAT(13) ** XMAT(11))
  131.          ENDIF
  132.  
  133.          DEPSCP = COE12
  134.          DEPSCS = COE22
  135.          VARP1(2) = DEPSCP
  136.          VARP1(3) = DEPSCS
  137.  
  138.        ENDIF
  139.  
  140.  
  141.        DEPST = DEPSCP + DEPSCS
  142.        CONS=1.5D0*DEPST/SEQUV
  143.  
  144.        DO 70 I=1,NSTRS
  145.    70  EPS0(I)=CONS*SS(I)
  146.  
  147.        DO 11 I=1,NSTRS
  148.          EVP1(I)=EPS0(I)
  149.  11    CONTINUE
  150.        VARP1(1) = DEPST
  151. C        print*,VARP1(1),EVP1(1),EVP1(2),EVP1(3)
  152.  
  153.  
  154.  
  155.         RETURN
  156.        END
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  

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