Sources Esope | Cast3M

ecdirr
C ECDIRR    SOURCE    STRU      08/12/22    21:15:10     6246
c     Sous-Programme du module MISTRAL0
C     --------------------------------------------------------------------------
      SUBROUTINE ECDIRR (TT,EPSEQ,PSI, ID, PEDIR, NPEDIR, RR,
     &                   DRRDTT,DRDEPS,DRDPSI)
C     --------------------------------------------------------------------------
C     Calcul, pour les deformations plastiques de nature IDP,
C     de la contrainte seuil RR
C     et de ses derivees partielles DRRDTT, DRDEPS et DRDPSI si ID = 1,
C     en fonction :
C     - de la temperature TT,
C     - de la deformation equivalente EPSEQ,
C     - de la variable de durcissement du a l'irradiation PSI,
C     et des parametres transmis dans PEDIR,
C     NPEDIR (incremente ici) etant en entree le nombre de ceux
C     correspondant aux deformations de natures IDPP &lt; IDP.
C     --------------------------------------------------------------------------
      IMPLICIT INTEGER (I-N)
      IMPLICIT REAL*8 (A-H, O-Z)
      PARAMETER ( NPTTM = 5 , NPPSIM = 6 )
      PARAMETER ( NPEDIM = 13+6*NPTTM+NPPSIM )
c     (NPEDIM = 49)
      DIMENSION PEDIR(1:*)
      DIMENSION PECDI(1:NPEDIM)
      NPEDIR = NPEDIR+1
      KRR = NINT(PEDIR(NPEDIR+1))
        IF (KRR.EQ.0) THEN
       NPEDIR = NPEDIR+1
         RR = 0.D0
       IF (ID.EQ.1) THEN
        DRRDTT = 0.D0
        DRDEPS = 0.D0
        DRDPSI = 0.D0
       END IF
        ELSE IF (KRR.EQ.-1) THEN
       NPEDIR = NPEDIR+2
       RR = PEDIR(NPEDIR)
       IF (ID.EQ.1) THEN
        DRRDTT = 0.D0
        DRDEPS = 0.D0
        DRDPSI = 0.D0
       END IF
        ELSE
       NPECDI = NINT(PEDIR(NPEDIR))
       DO IP = 1,NPECDI
        PECDI(IP) = PEDIR(NPEDIR+IP)
       END DO
       NPEDIR = NPEDIR+NPECDI
       CALL FECRDI (TT,EPSEQ,PSI, ID, PECDI, RR,DRRDTT,DRDEPS,DRDPSI)
        END IF
        RETURN
        END