Esope source files

   1 : C ECOUV     SOURCE    STRU      08/12/22    21:15:13     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE ECOUV (SEQ,RR,EPSEQ,TT,FI,PSI, DSESDS, IDP, PECOU,
   5 :      &                  NPECOU, VEPSP,VEPSE)
   6 : C     --------------------------------------------------------------------------
   7 : C     Calcul, pour les deformations viscoplastiques de nature IDP :
   8 : C     - de la vitesse de deformation equivalente VEPSE,
   9 : C     - des vitesses de deformation VEPSP(IJ,IDP),
  10 : C     en fonction :
  11 : C     - de la contrainte effective equivalente SEQ,
  12 : C     - de la contrainte seuil RR,
  13 : C     - de la deformation equivalente EPSEQ,
  14 : C     - de la temperature TT,
  15 : C     - du flux de neutrons rapides FI,
  16 : C     - de la variable de durcissement du a l'irradiation PSI,
  17 : C     - des derivees partielles DSESDS(IJ) de SEQ par rapport aux contraintes
  18 : C     "effectives" S(IJ) = SIG(IJ)-XX(IJ,IDP),
  19 : C     - des parametres de la loi d'ecoulement transmis dans PECOU,
  20 : C     NPECOU (incremente ici) etant en entree le nombre de ceux
  21 : C     correspondant aux deformations de natures IDPP < IDP.
  22 : C     --------------------------------------------------------------------------
  23 :       IMPLICIT INTEGER (I-N)
  24 :       IMPLICIT REAL*8 (A-H, O-Z)
  25 :       PARAMETER ( NPTTM = 5 , NPPSIM = 6 )
  26 :       PARAMETER ( NPEDIM = 13+6*NPTTM+NPPSIM )
  27 : c     (NPEDIM = 49)
  28 :       PARAMETER ( NVM=2 , NPATDM = 14 , NPFIM = 9 )
  29 :       PARAMETER ( NPERVM = 11+NVM*(NPEDIM+2*NPATDM+NPFIM+11) )
  30 : c     (NPERVM = 205)
  31 :       DIMENSION DSESDS(1:*), PECOU(1:*)
  32 :       DIMENSION PERV(1:NPERVM)
  33 :         DIMENSION VEPSP(1:6,0:*)
  34 :       NPECOU = NPECOU+1
  35 :       NPERV = NINT(PECOU(NPECOU))
  36 :         IF (SEQ.LE.RR) THEN
  37 :          DO IJ = 1,6
  38 :           VEPSP(IJ,IDP) = 0.
  39 :          END DO
  40 :        VEPSE = 0.
  41 :         ELSE
  42 :        DO IP = 1,NPERV
  43 :         PERV(IP) = PECOU(NPECOU+IP)
  44 :        END DO
  44 :        CALL ECOREV (SEQ,RR,EPSEQ,TT,FI,PSI, PERV, VEPSE)
  45 :        DO IJ = 1,6
  46 :           VEPSP(IJ,IDP) = DSESDS(IJ)*VEPSE
  47 :          END DO
  48 :         END IF
  49 :       NPECOU = NPECOU+NPERV
  50 :         RETURN
  51 :         END
  52 :