Esope source files

   1 : C ECDIRR    SOURCE    STRU      08/12/22    21:15:10     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE ECDIRR (TT,EPSEQ,PSI, ID, PEDIR, NPEDIR, RR,
   5 :      &                   DRRDTT,DRDEPS,DRDPSI)
   6 : C     --------------------------------------------------------------------------
   7 : C     Calcul, pour les deformations plastiques de nature IDP,
   8 : C     de la contrainte seuil RR
   9 : C     et de ses derivees partielles DRRDTT, DRDEPS et DRDPSI si ID = 1,
  10 : C     en fonction :
  11 : C     - de la temperature TT,
  12 : C     - de la deformation equivalente EPSEQ,
  13 : C     - de la variable de durcissement du a l'irradiation PSI,
  14 : C     et des parametres transmis dans PEDIR,
  15 : C     NPEDIR (incremente ici) etant en entree le nombre de ceux
  16 : C     correspondant aux deformations de natures IDPP < IDP.
  17 : C     --------------------------------------------------------------------------
  18 :       IMPLICIT INTEGER (I-N)
  19 :       IMPLICIT REAL*8 (A-H, O-Z)
  20 :       PARAMETER ( NPTTM = 5 , NPPSIM = 6 )
  21 :       PARAMETER ( NPEDIM = 13+6*NPTTM+NPPSIM )
  22 : c     (NPEDIM = 49)
  23 :       DIMENSION PEDIR(1:*)
  24 :       DIMENSION PECDI(1:NPEDIM)
  25 :       NPEDIR = NPEDIR+1
  26 :       KRR = NINT(PEDIR(NPEDIR+1))
  27 :         IF (KRR.EQ.0) THEN
  28 :        NPEDIR = NPEDIR+1
  29 :          RR = 0.D0
  30 :        IF (ID.EQ.1) THEN
  31 :         DRRDTT = 0.D0
  32 :         DRDEPS = 0.D0
  33 :         DRDPSI = 0.D0
  34 :        END IF
  35 :         ELSE IF (KRR.EQ.-1) THEN
  36 :        NPEDIR = NPEDIR+2
  37 :        RR = PEDIR(NPEDIR)
  38 :        IF (ID.EQ.1) THEN
  39 :         DRRDTT = 0.D0
  40 :         DRDEPS = 0.D0
  41 :         DRDPSI = 0.D0
  42 :        END IF
  43 :         ELSE
  44 :        NPECDI = NINT(PEDIR(NPEDIR))
  45 :        DO IP = 1,NPECDI
  46 :         PECDI(IP) = PEDIR(NPEDIR+IP)
  47 :        END DO
  48 :        NPEDIR = NPEDIR+NPECDI
  48 :        CALL FECRDI (TT,EPSEQ,PSI, ID, PECDI, RR,DRRDTT,DRDEPS,DRDPSI)
  49 :         END IF
  50 :         RETURN
  51 :         END
  52 :