Esope source files

   1 : C FATDI     SOURCE    STRU      08/12/22    21:15:19     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE FATDI (TT,S,PSI, PAR, V)
   5 : C     --------------------------------------------------------------------------
   6 : C     Calcul de V = FTT(TT,S,PSI)*FPSI(PSI) en fonction :
   7 : C     - de la temperature TT,
   8 : C     - d'une contrainte S,
   9 : C     - de la variable de durcissement du a l'irradiation PSI,
  10 : C     et des parametres transmis dans PAR.
  11 : C     --------------------------------------------------------------------------
  12 :       IMPLICIT INTEGER (I-N)
  13 :       IMPLICIT REAL*8 (A-H, O-Z)
  14 :       DIMENSION PAR(1:*)
  15 : c     (NPATDM = 14)
  16 : C     ------------------------------------------------------------------
  17 :       L = NINT(PAR(1))
  18 :       IP = 1
  19 :         IF (L.EQ.-1) THEN
  20 :        IP = IP+1
  21 :          FTT = PAR(IP)
  22 :         ELSE IF (L.EQ.1) THEN
  23 :          FTT = PAR(IP+1)*EXP(-(PAR(IP+2)-PAR(IP+3)*S)/TT)
  24 :        IP = IP+3
  25 :         ELSE IF (L.EQ.2) THEN
  26 :        FT = 1./(1.+EXP(PAR(IP+3)*(TT/PAR(IP+4)-1.)))
  27 :        FTT = PAR(IP+1)*FT+PAR(IP+2)*(1.-FT)
  28 :        IP = IP+4
  29 :         ELSE IF (L.EQ.3) THEN
  30 :          FTT = PAR(IP+1)
  31 :      &  *EXP(-(PAR(IP+2)+PAR(IP+3)*EXP(PAR(IP+4)*PSI**PAR(IP+5)))/TT)
  32 :        IP = IP+5
  33 :         ELSE IF (L.EQ.4) THEN
  34 :          FTT = PAR(IP+1)*EXP(PAR(IP+2)*TT)
  35 :        IP = IP+2
  36 :         END IF
  37 : C     ------------------------------------------------------------------
  38 :       IP = IP+1
  39 :       L = NINT(PAR(IP))
  40 :         IF (L.EQ.0) THEN
  41 :          V = FTT
  42 :        RETURN
  43 : C     ------------------------------------------------------------------
  44 :         ELSE IF (L.EQ.1) THEN
  45 :          FPSI = (PSI+PAR(IP+1))**PAR(IP+2)
  46 :         ELSE IF (L.EQ.2) THEN
  47 :          FPSI = PAR(IP+1)*EXP(PAR(IP+2)*PSI**PAR(IP+3))
  48 :      & + (1.-PAR(IP+1))*EXP(PAR(IP+4)*PSI**PAR(IP+5))
  49 :         ELSE IF (L.EQ.3) THEN
  50 :          PSIPL = PSI+PAR(IP+2)
  51 :          FPSI = PAR(IP+1)*PSIPL**PAR(IP+3)*EXP(-PAR(IP+4)*PSI)
  52 :      &      + PAR(IP+5)*PSIPL**PAR(IP+6)*(1.-EXP(-PAR(IP+7)*PSI))
  53 :         ELSE IF (L.EQ.4) THEN
  54 :          FPSI = EXP(PAR(IP+1)*PSI**PAR(IP+2)+PAR(IP+3)*PSI**PAR(IP+4))
  55 :         ELSE IF (L.EQ.5) THEN
  56 :          FPSI = PAR(IP+1) + (1.-PAR(IP+1))*EXP(PAR(IP+2)*PSI
  57 :      &      +PAR(IP+3)*PSI**PAR(IP+5)+PAR(IP+4)*PSI**PAR(IP+6))
  58 :         ELSE IF (L.EQ.6) THEN
  59 :          FPSI = PAR(IP+1) + PAR(IP+2)*PSI + PAR(IP+3)*PSI**2
  60 :         ELSE IF (L.EQ.7) THEN
  61 :          FPSI = (PAR(IP+1)*PSI+PAR(IP+2)) / (PAR(IP+3)*PSI+PAR(IP+4))
  62 :          FPSI = FPSI**PAR(IP+5)
  63 :         END IF
  64 : C     ------------------------------------------------------------------
  65 :         V = FTT*FPSI
  66 :         RETURN
  67 :         END
  68 :