C ECRISR    SOURCE    PV        11/03/07    21:16:37     6885           
c     Sous-Programme du module MISTRAL
C     --------------------------------------------------------------------------
      SUBROUTINE ECRISR (TT,EPSEQ,PSI, IDP,NDPI,PECRI, RR)
C     --------------------------------------------------------------------------
C     Calcul, pour les deformations plastiques de nature IDP,
C     de la contrainte seuil RR
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     --------------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H, O-Z)
      DIMENSION PECRI(1:*)
      NPECRI = 29*(IDP+NDPI-1)
      KRR = nint(PECRI(NPECRI+1))
C     ------------------------------------------------------------------
        IF (KRR.EQ.0) THEN
         RR = 0.
         RETURN
        END IF
C     ------------------------------------------------------------------
      EPSE0 = PECRI(NPECRI+2)
      G1 = PECRI(NPECRI+3)
      G2 = PECRI(NPECRI+4)
      LRRI = nint(PECRI(NPECRI+5))
      R10 = PECRI(NPECRI+6)
        IF (LRRI.EQ.1) THEN
       R11 = PECRI(NPECRI+7)
       R12 = PECRI(NPECRI+8)
        ELSE
       R1M1 = PECRI(NPECRI+7)
       R11 = PECRI(NPECRI+8)
       R12 = PECRI(NPECRI+9)
        END IF
      EN10 = PECRI(NPECRI+10)
      EN11 = PECRI(NPECRI+11)
      EN12 = PECRI(NPECRI+12)
      R20 = PECRI(NPECRI+13)
        IF (LRRI.EQ.1) THEN
       R21 = PECRI(NPECRI+14)
       R22 = PECRI(NPECRI+15)
        ELSE
       R2M1 = PECRI(NPECRI+14)
       R21 = PECRI(NPECRI+15)
       R22 = PECRI(NPECRI+16)
        END IF
      EN20 = PECRI(NPECRI+17)
      EN21 = PECRI(NPECRI+18)
      EN22 = PECRI(NPECRI+19)
      LRPSI = nint(PECRI(NPECRI+20))
      R00 = PECRI(NPECRI+21)
      RPSI1 = PECRI(NPECRI+22)
      RPSI2 = PECRI(NPECRI+23)
        IF (LRPSI.EQ.2) THEN
       RPSI3 = PECRI(NPECRI+24)
       RPSI4 = PECRI(NPECRI+25)
       RPSI5 = PECRI(NPECRI+26)
        END IF
C     ------------------------------------------------------------------
        IF (LRRI.EQ.1) THEN
         RR1 = R10+(R11+R12*TT)*TT
         RR2 = R20+(R21+R22*TT)*TT
        ELSE
         RR1 = R10*EXP(R1M1/TT+(R11+R12*TT)*TT)
         RR2 = R20*EXP(R2M1/TT+(R21+R22*TT)*TT)
        END IF
        EN1 = EN10+(EN11+EN12*TT)*TT
        EN2 = EN20+(EN21+EN22*TT)*TT
        RR = R00 + RR1*(EPSEQ+EPSE0)**EN1*EXP(-G1*EPSEQ)
     &         + RR2*(EPSEQ+EPSE0)**EN2*(1.-EXP(-G2*EPSEQ))
        IF (LRPSI.EQ.1) THEN
         RRPSI = RPSI1*LOG(1.+RPSI2*PSI)
        ELSE
         FT = 1./(1.+EXP(RPSI3*(TT/RPSI4-1.)))
         RRPSI = (RPSI1*FT+RPSI2*(1.-FT))*(1.-EXP(-RPSI5*PSI))
        END IF
        RR = RR+RRPSI
        RETURN
        END