Télécharger la plaquette Cast3M


fr en

Télécharger @SATURAT.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @SATURAT  PROCEDUR  PASCAL    08/10/06    21:15:05     6176           
  2. * PROCEDURE @SATURAT
  3. *---------------------------------------------------------------------*
  4. * NOM         : @SATURAT                                              *
  5. *                                                                     *
  6. * DESCRIPTION :                                                       *
  7. *                                                                     *
  8. * SYNTAXE     :                                                       *
  9. *                                                                     *
  10. *                                                                     *
  11. * LANGAGE     : GIBIANE-CAST3M                                        *
  12. * AUTEUR      : S. PASCAL (CEA/DEN/DM2S/SEMT/LM2S)                    *
  13. * COURRIEL    : serge.pascal@cea.fr                                   *
  14. *---------------------------------------------------------------------*
  15. * VERSION    : v1, xx/xx/xx08, version initiale                       *
  16. * HISTORIQUE : v1, xx/xx/xx08, creation                               *
  17. * HISTORIQUE :                                                        *
  18. * HISTORIQUE :                                                        *
  19. *---------------------------------------------------------------------*
  20. * Priere de PRENDRE LE TEMPS de completer les commentaires            *
  21. * en cas de modification de ce sous-programme afin de faciliter       *
  22. * la maintenance !                                                    *
  23. *---------------------------------------------------------------------*
  24. 'DEBP' @SATURAT ;
  25. 'ARGU' CHTK1*'CHPOINT' CHPC1*'CHPOINT' ;
  26. * Ok !
  27. * ENTREES                                                             *
  28. *  - Temperature (K)                                                  *
  29. *  - Pression capillaire (Pa)                                         *
  30. * SORTIE                                                              *
  31. *  - Saturation (-) : SW0                                             *
  32. * PARAMETRES                                                          *
  33. * |-------------|---------|---------|---------|-----------|
  34. * |coefficient  | BO      |   BHP   |   MO    |     MH    |
  35. * ||------------|---------|---------|---------|-----------|
  36. * |a[MPa]       |18.6237  |46.9364  |37.5479  |96.2837    |
  37. * |-------------|---------|---------|---------|-----------|
  38. * |b[-]         |2.2748   |2.0601   |2.1684   |1.9540     |
  39. * |-------------|---------|---------|---------|-----------|
  40. * Q2       : parametre de calcul du parametre "a", relatif a T>Tcrit  *
  41. * Q3       : parametre de calcul du parametre "a", relatif a T=20°C   *
  42. * PB1      : exposant "b"                                             *
  43. * GN1      : exposant "N" de E(theta)                                 *
  44. * PZ1      : parametre "z" de E(theta) qd T > TCRI1                   *
  45. * E0       : valeur inconnue => 1. !                                  *
  46. * TREF1    : Temperature (K) "ambiante" (initiale ?)                  *
  47. * TSAT1    : Temperature (K) au-dessus de laquelle le param. "a" vari *
  48. * TCRI1    : Temperature (K) critique                                 *
  49. Q2       = 7.E6 ;
  50. Q3       = 18.62E6 ;
  51. PB1      = 2.27 ;
  52. GN1      = 1.2 ;
  53. PZ1      = 0.5 ;
  54. E0       = 1. ;
  55. TREF1    = 293.15 ;
  56. TSAT1    = 373.15 ; 
  57. TCRI1    = 647.15 ;
  58. *                                                                     *
  59. * Calcul du parametre "a" : PA0                                       *
  60. Q0       = (CHTK1 - TSAT1) / (TCRI1 - TSAT1) ;
  61. Q02      = -3. * (Q0 ** 2) ;
  62. Q03      =  2. * (Q0 ** 3) ;
  63. Q0       = (Q3 - Q2) * (1. + Q02 + Q03) ;
  64. MASQI1   = CHTK1 'MASQ' 'EGIN' TSAT1 ;
  65. MASQS1   = CHTK1 'MASQ' 'SUPE' TSAT1 ;
  66. PA0      = (Q3 * MASQI1) + ((Q0 + Q2) * MASQS1) ;
  67. *                                                                     *
  68. * Calcul de E(theta) : ETHETA0                                        *
  69. MASQI2   = CHTK1 'MASQ' 'INFE' TCRI1 ;
  70. MASQS2   = CHTK1 'MASQ' 'EGSU' TCRI1 ;
  71. CHTKSC1  = MASQI2 * CHTK1 ;
  72. ETHETA1  = (((TCRI1 - CHTKSC1) ** -1) * (TCRI1 - TREF1)) ** GN1 ;
  73. ETHETA2  = GN1 / PZ1 * (CHTK1 - TCRI1 + PZ1) ;
  74. ETHETA2  = E0 * (1. + ETHETA2) ;
  75. * ATTENTION : jeux de donnees SdP ne distinguent pas T>TSAT1...       *
  76. *ETHETA0  = MASQI1 + (ETHETA1 * MASQI2 * MASQS1) + (ETHETA2 * MASQS2) ;
  77. ETHETA0  = ETHETA1 ;
  78. *                                                                     *
  79. * La saturation :                                                     *
  80. M1SPB1   = -1. / PB1 ;
  81. PA0      = PA0 ** -1 ;
  82. SW0      = (ETHETA0 * PA0 * CHPC1) ** (1. / (1. + M1SPB1)) + 1. ;
  83. SW0      = SW0 ** M1SPB1 ;
  84. *                                                                     *
  85. 'RESP' SW0 ;
  86. *                                                                     *
  87. 'FINP' ;
  88. *---------------------------------------------------------------------*
  89. * FIN PROCEDURE @SATURAT
  90.  

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