Télécharger cohill.eso

Retour à la liste

Numérotation des lignes :

  1. C COHILL SOURCE STRU 08/12/22 21:15:05 6246
  2. c Sous-Programme du module MISTRAL0
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE COHILL (TT,PSI, ID, PCOHI, NP, H,DHDTT,DHDPSI)
  5. C --------------------------------------------------------------------------
  6. C Calcul des coefficients de Hill H(I),
  7. C et de leur derivees partielles DHDTT(I) et DHDPSI(I) si ID = 1,
  8. C en fonction de :
  9. C - la temperature TT,
  10. C - la variable de durcissement du a l'irradiation PSI,
  11. C - des parametres transmis dans PCOHI,
  12. C NP (incremente ici) etant en entree le nombre
  13. C de ceux correspondant aux deformations de natures IDPP < IDP.
  14. C --------------------------------------------------------------------------
  15. IMPLICIT INTEGER (I-N)
  16. IMPLICIT REAL*8 (A-H, O-Z)
  17. DIMENSION PCOHI(1:*)
  18. c (NPCHM = 6*10, NPCOHM = NPCHM*NDP)
  19. DIMENSION H(1:*),DHDTT(1:*),DHDPSI(1:*)
  20. DO I = 1,6
  21. NP = NP+1
  22. KH = NINT(PCOHI(NP))
  23. IF (KH.EQ.1) THEN
  24. NP = NP+1
  25. H(I) = PCOHI(NP)
  26. IF (ID.EQ.1) THEN
  27. DHDTT(I) = 0.D0
  28. DHDPSI(I) = 0.D0
  29. END IF
  30. ELSE IF (KH.EQ.2) THEN
  31. H(I) = PCOHI(NP+1)+(PCOHI(NP+2)+PCOHI(NP+3)*TT)*TT
  32. IF (ID.EQ.1) THEN
  33. DHDTT(I) = PCOHI(NP+2)+2.*PCOHI(NP+3)*TT
  34. DHDPSI(I) = 0.D0
  35. END IF
  36. NP = NP+3
  37. ELSE IF (KH.EQ.3) THEN
  38. HI0 = PCOHI(NP+1)
  39. HI1 = PCOHI(NP+2)
  40. HI2 = PCOHI(NP+3)
  41. HI3 = PCOHI(NP+4)
  42. HI4 = PCOHI(NP+5)
  43. HI5 = PCOHI(NP+6)
  44. HI6 = PCOHI(NP+7)
  45. HI7 = PCOHI(NP+8)
  46. HI8 = PCOHI(NP+9)
  47. EMHI7F = EXP(-HI7*PSI)
  48. EMHI8F = EXP(-HI8*PSI)
  49. HIUN = HI1+(HI2+HI3*TT)*TT
  50. HIDE = HI4+(HI5+HI6*TT)*TT
  51. H(I) = HI0+HIUN*EMHI7F+HIDE*(1.-EMHI8F)
  52. IF (ID.EQ.1) THEN
  53. DHDTT(I) = (HI2+2.*HI3*TT)*EMHI7F+(HI5+2.*HI6*TT)*(1.-EMHI8F)
  54. DHDPSI(I) = -HI7*HIUN*EMHI7F+HI8*HIDE*EMHI8F
  55. END IF
  56. NP = NP+9
  57. ELSE IF (KH.EQ.4) THEN
  58. HI0 = PCOHI(NP+1)
  59. HI1 = PCOHI(NP+2)
  60. HI2 = PCOHI(NP+3)
  61. HI3 = PCOHI(NP+4)
  62. HI4 = PCOHI(NP+5)
  63. HI5 = PCOHI(NP+6)
  64. HI6 = PCOHI(NP+7)
  65. HI7 = PCOHI(NP+8)
  66. HI8 = PCOHI(NP+9)
  67. EMHPSI = EXP(-HI7*PSI**HI8)
  68. HIUN = HI1+(HI2+HI3*TT)*TT
  69. HIDE = HI4+(HI5+HI6*TT)*TT
  70. H(I) = HI0+HIUN*EMHPSI+HIDE*(1.-EMHPSI)
  71. IF (ID.EQ.1) THEN
  72. DHDTT(I) = (HI2+2.*HI3*TT)*EMHPSI+(HI5+2.*HI6*TT)*(1.-EMHPSI)
  73. DHDPSI(I) = -(HIUN-HIDE)*HI7*HI8*PSI**(HI8-1.)*EMHPSI
  74. END IF
  75. NP = NP+9
  76. END IF
  77. END DO
  78. RETURN
  79. END
  80.  
  81.  

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