Esope source files

   1 : C ECRISR    SOURCE    PV        11/03/07    21:16:37     6885           
   2 : c     Sous-Programme du module MISTRAL
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE ECRISR (TT,EPSEQ,PSI, IDP,NDPI,PECRI, RR)
   5 : C     --------------------------------------------------------------------------
   6 : C     Calcul, pour les deformations plastiques de nature IDP,
   7 : C     de la contrainte seuil RR
   8 : C     en fonction :
   9 : C     - de la temperature TT,
  10 : C     - de la deformation equivalente EPSEQ,
  11 : C     - de la variable de durcissement du a l'irradiation PSI.
  12 : C     --------------------------------------------------------------------------
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8 (A-H, O-Z)
  15 :       DIMENSION PECRI(1:*)
  16 :       NPECRI = 29*(IDP+NDPI-1)
  17 :       KRR = nint(PECRI(NPECRI+1))
  18 : C     ------------------------------------------------------------------
  19 :         IF (KRR.EQ.0) THEN
  20 :          RR = 0.
  21 :          RETURN
  22 :         END IF
  23 : C     ------------------------------------------------------------------
  24 :       EPSE0 = PECRI(NPECRI+2)
  25 :       G1 = PECRI(NPECRI+3)
  26 :       G2 = PECRI(NPECRI+4)
  27 :       LRRI = nint(PECRI(NPECRI+5))
  28 :       R10 = PECRI(NPECRI+6)
  29 :         IF (LRRI.EQ.1) THEN
  30 :        R11 = PECRI(NPECRI+7)
  31 :        R12 = PECRI(NPECRI+8)
  32 :         ELSE
  33 :        R1M1 = PECRI(NPECRI+7)
  34 :        R11 = PECRI(NPECRI+8)
  35 :        R12 = PECRI(NPECRI+9)
  36 :         END IF
  37 :       EN10 = PECRI(NPECRI+10)
  38 :       EN11 = PECRI(NPECRI+11)
  39 :       EN12 = PECRI(NPECRI+12)
  40 :       R20 = PECRI(NPECRI+13)
  41 :         IF (LRRI.EQ.1) THEN
  42 :        R21 = PECRI(NPECRI+14)
  43 :        R22 = PECRI(NPECRI+15)
  44 :         ELSE
  45 :        R2M1 = PECRI(NPECRI+14)
  46 :        R21 = PECRI(NPECRI+15)
  47 :        R22 = PECRI(NPECRI+16)
  48 :         END IF
  49 :       EN20 = PECRI(NPECRI+17)
  50 :       EN21 = PECRI(NPECRI+18)
  51 :       EN22 = PECRI(NPECRI+19)
  52 :       LRPSI = nint(PECRI(NPECRI+20))
  53 :       R00 = PECRI(NPECRI+21)
  54 :       RPSI1 = PECRI(NPECRI+22)
  55 :       RPSI2 = PECRI(NPECRI+23)
  56 :         IF (LRPSI.EQ.2) THEN
  57 :        RPSI3 = PECRI(NPECRI+24)
  58 :        RPSI4 = PECRI(NPECRI+25)
  59 :        RPSI5 = PECRI(NPECRI+26)
  60 :         END IF
  61 : C     ------------------------------------------------------------------
  62 :         IF (LRRI.EQ.1) THEN
  63 :          RR1 = R10+(R11+R12*TT)*TT
  64 :          RR2 = R20+(R21+R22*TT)*TT
  65 :         ELSE
  66 :          RR1 = R10*EXP(R1M1/TT+(R11+R12*TT)*TT)
  67 :          RR2 = R20*EXP(R2M1/TT+(R21+R22*TT)*TT)
  68 :         END IF
  69 :         EN1 = EN10+(EN11+EN12*TT)*TT
  70 :         EN2 = EN20+(EN21+EN22*TT)*TT
  71 :         RR = R00 + RR1*(EPSEQ+EPSE0)**EN1*EXP(-G1*EPSEQ)
  72 :      &         + RR2*(EPSEQ+EPSE0)**EN2*(1.-EXP(-G2*EPSEQ))
  73 :         IF (LRPSI.EQ.1) THEN
  74 :          RRPSI = RPSI1*LOG(1.+RPSI2*PSI)
  75 :         ELSE
  76 :          FT = 1./(1.+EXP(RPSI3*(TT/RPSI4-1.)))
  77 :          RRPSI = (RPSI1*FT+RPSI2*(1.-FT))*(1.-EXP(-RPSI5*PSI))
  78 :         END IF
  79 :         RR = RR+RRPSI
  80 :         RETURN
  81 :         END
  82 :  
  83 :  
  84 :