Télécharger appfib.eso

Retour à la liste

Numérotation des lignes :

  1. C APPFIB SOURCE CHAT 10/12/16 21:15:13 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 INTIMP
  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 APPFIB(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) = 5.00D3
  49. XYR( 4) = 5.30D3
  50. XYR( 5) = 6.50D3
  51. XYR( 6) = 6.00D3
  52. XYR( 7) = 5.00D3
  53. XYR( 8) = 4.80D3
  54. XYR( 9) = -2.00D3
  55. XYR(10) = -4.00D3
  56. XYR(11) = -1.00D4
  57. XYR(12) = -2.50D4
  58. XYR(13) = -7.00D4
  59. XYR(14) = -5.00D5
  60. XYR(15) = -1.00D6
  61. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  62. C
  63. C-----TEST SI LE DEGRE DE CORROSION MAXI IDENTIFIE
  64. C N EST PAS DEPASSE
  65. C
  66. IF (X.GT.XTC(NDIM)) THEN
  67. Y = XYR(NDIM)
  68. RETURN
  69. END IF
  70. C
  71. C-----INTERPOLATION LINEAIRE PAR MORCEAU
  72. C
  73. DO I=1,(NDIM-1)
  74. IF ((XTC(I).LE.X).AND.(X.LE.XTC(I+1))) THEN
  75. Y=(XYR(I+1)-XYR(I))/(XTC(I+1)-XTC(I))*X+
  76. & (XTC(I+1)*XYR(I)-XTC(I)*XYR(I+1))/(XTC(I+1)-XTC(I))
  77. RETURN
  78. END IF
  79. END DO
  80.  
  81. RETURN
  82. END
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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