Télécharger la plaquette Cast3M


fr en

Télécharger app2.eso

Retour à la liste

Numérotation des lignes : 

app2

  1. C APP2      SOURCE    CHAT      10/12/16    21:15:02     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 RICJOI2D
  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 APP2(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) =  2.80D4
  49.         XYR( 4) =  3.20D4
  50.         XYR( 5) =  4.00D4
  51.         XYR( 6) =  3.50D4
  52.         XYR( 7) =  2.50D4
  53.         XYR( 8) =  2.25D4
  54.         XYR( 9) = -7.00D3
  55.         XYR(10) = -2.60D3
  56.         XYR(11) = -6.00D4
  57.         XYR(12) = -1.50D5
  58.         XYR(13) = -3.50D5
  59.         XYR(14) = -1.50D6
  60.         XYR(15) = -3.00D6
  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