Télécharger app3.eso

Retour à la liste

Numérotation des lignes :

  1. C APP3 SOURCE CHAT 10/12/16 21:15:12 6819
  2. C
  3. C====&===1=========2=========3=========4=========5=========6=========7==
  4. C Commentaires : Subroutine permettant de définir la relation entre
  5. C le taux de corrosion macroscopique TC et le taux d'éner-
  6. C gie additionel à prendre en compte pour modifier l'endom-
  7. C magement dans le cas du modèle RICJOI3D
  8. C avec interface imparfaite
  9. C Auteurs : B. Richard
  10. C Date : Février 2008
  11. C====&===1=========2=========3=========4=========5=========6=========7==
  12. SUBROUTINE APP3(X,Y,NDIM)
  13. C
  14. C-----DECLARATION GENERALE
  15. C
  16. IMPLICIT INTEGER(I-N)
  17. IMPLICIT REAL*8(A-H,O-Z)
  18. C
  19. C-----DECLARATION DES TYPE DE VARIABLES EN ARGUMENTS
  20. C
  21. REAL*8 XTC(20),XYR(20)
  22. C
  23. C-----VARIABLES LOCALES
  24. C
  25. INTEGER I
  26. C
  27. C-----DEFINITION DES TABLEAUX DE REFERENCE
  28. C
  29. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  30. XTC( 1) = 0.00D0
  31. XTC( 2) = 1.62D0
  32. XTC( 3) = 5.61D0
  33. XTC( 4) = 8.58D0
  34. XTC( 5) = 9.45D0
  35. XTC( 6) = 11.58D0
  36. XTC( 7) = 12.35D0
  37. XTC( 8) = 12.66D0
  38. XTC( 9) = 16.54D0
  39. XTC(10) = 17.09D0
  40. XTC(11) = 18.06D0
  41. XTC(12) = 20.45D0
  42. XTC(13) = 22.36D0
  43. XTC(14) = 22.99D0
  44. XTC(15) = 25.07D0
  45.  
  46. XYR( 1) = 0.00D0
  47. XYR( 2) = 0.00D0
  48. XYR( 3) = 1.40D4
  49. XYR( 4) = 1.45D4
  50. XYR( 5) = 1.65D4
  51. XYR( 6) = 1.60D4
  52. XYR( 7) = 1.30D4
  53. XYR( 8) = 1.25D4
  54. XYR( 9) = 2.00D3
  55. XYR(10) = -1.30D3
  56. XYR(11) = -1.00D4
  57. XYR(12) = -3.50D4
  58. XYR(13) = -8.00D4
  59. XYR(14) = -3.00D5
  60. XYR(15) = -6.00D5
  61. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  62.  
  63. C
  64. C-----TEST SI LE DEGRE DE CORROSION MAXI IDENTIFIE
  65. C N EST PAS DEPASSE
  66. C
  67. IF (X.GT.XTC(NDIM)) THEN
  68. Y = XYR(NDIM)
  69. RETURN
  70. END IF
  71. C
  72. C-----INTERPOLATION LINEAIRE PAR MORCEAU
  73. C
  74. DO I=1,(NDIM-1)
  75. IF ((XTC(I).LE.X).AND.(X.LE.XTC(I+1))) THEN
  76. Y=(XYR(I+1)-XYR(I))/(XTC(I+1)-XTC(I))*X+
  77. & (XTC(I+1)*XYR(I)-XTC(I)*XYR(I+1))/(XTC(I+1)-XTC(I))
  78. RETURN
  79. END IF
  80. END DO
  81.  
  82. RETURN
  83. END
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  

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