Télécharger la plaquette Cast3M


fr en

Télécharger fatdi.eso

Retour à la liste

Numérotation des lignes : 

fatdi

  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.  
  69.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales