Télécharger la plaquette Cast3M


fr en

Télécharger appfib.eso

Retour à la liste

Numérotation des lignes : 

appfib

  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