Télécharger fatdi.eso

Retour à la liste

Numérotation des lignes :

  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