Télécharger la plaquette Cast3M


fr en

Télécharger @ch_dama.procedur

Retour à la liste

Numérotation des lignes : 

  1. * @CH_DAMA  PROCEDUR  FD218221  26/02/16    21:15:05     12474          
  2. *---------------------------------------------------------------------*
  3. * @CH_DAMA  PROCEDUR                                                  *
  4. * PROCEDURE @CH_DAMA                                                  *
  5. * APPELLE PAR : @BET_THM                                              *
  6. *---------------------------------------------------------------------*
  7. * NOM         : @CH_DAMA                                              *
  8. *                                                                     *
  9. * LANGAGE     : GIBIANE-CAST3M                                        *
  10. * AUTEUR      : G. Sciumè (I2M - University of Bordeaux)              *
  11. * COURRIEL    : giuseppe.sciume@u-bordeuax.fr                         *
  12. *---------------------------------------------------------------------*
  13. * COMMENTAIRE: Mise a jour du degre d'hydratation                     *
  14. *---------------------------------------------------------------------*
  15. * VERSION    : v1, 31/10/2024, version initiale                       *
  16. * HISTORIQUE : v2, XX/XX/20XX,                                        *
  17. *---------------------------------------------------------------------*
  18. DEBPROC @CH_DAMA WORKTAB*'TABLE' CH_GLOB*'CHPOINT';
  19.  
  20.  
  21.  
  22. tab_zone = EXTR (WORKTAB . 'MOD_THM') 'ZONE';
  23. nzone = (dime tab_zone)/2;
  24.  
  25.  
  26. *** BOUCLE SUR LES DIFFERENTS ZONES THERMOHYDRIQUES
  27. indz = 1;
  28. 'REPETER' BOUCLE_Z nzone;
  29.  
  30. MOD_THM  = tab_zone.indz;
  31. MAIL_THM = 'EXTR' MOD_THM 'MAIL';
  32. SOLU_W  = 'REDU' MAIL_THM  CH_GLOB;
  33. *HYDRA_1 = 'REDU' MAIL_THM (WORKTAB . 'HYD1');
  34.  
  35.  
  36. SI ('EXIS' PRECED . WTABLE 'VAR_THM1') ;
  37.         VAR_THM1  = 'REDU' (PRECED . WTABLE . 'VAR_THM1' ) MOD_THM;
  38. SINO;
  39.         VAR_THM1  = 'REDU' (PRECED . 'VARIABLES_THM' . 0) MOD_THM;
  40. FINS;
  41.  
  42. DCH_1   = 'EXCO' 'DCH' VAR_THM1;
  43.  
  44. THETA_W = WORKTAB . 'RELAXATION_THETA';
  45. DELTA_T = WORKTAB . 'DT';
  46.  
  47. T_W = REDU MAIL_THM ((EXCO 'T' SOLU_W) 'NOMC' 'SCAL');
  48.  
  49.  
  50. * ENTREE                                                              *
  51. *  - Temperature (Celcius)                                            *
  52. * SORTIES                                                             *
  53. *  - Taux de deshydratation (kg/s) : DCHW                             *
  54. * PARAMETRES (from Gawin et al. 2011)                                 *
  55. * a1      : 1.7151 x 1E-3                                             *
  56. * a2      : 4.0006 x 1E-7                                             *
  57. * a3      : 2.9507 x 1E-10                                            *
  58. *
  59. a1 =  1.7151E-3  ;                                          
  60. a2 = -4.0006E-7  ;                                            
  61. a3 = -2.9507E-10 ;                                                                                  
  62. *
  63. SI (WORKTAB . 'HT_SOL1' 'EGA' 0.);
  64. CHTKSC1 = ('EXCO' PRECED . 'CONTINUATION' . 'TEMPERATURES' 'T') 'NOMC' 'SCAL' ;
  65. CHTKSC2 = ('EXCO' PRECED . 'CONTINUATION' . 'TEMPERATURES' 'T') 'NOMC' 'SCAL' ;
  66. SINO;
  67. CHTKSC1 = ('EXCO' WORKTAB . 'HT_SOL1' 'T') 'NOMC' 'SCAL' ;
  68. CHTKSC2 = ('EXCO' WORKTAB . 'HT_SOL2' 'T') 'NOMC' 'SCAL' ;
  69. FINS;
  70.  
  71. TCMTD1  = BORN (CHTKSC1 - 378.15) 'SCAL' 'MINIMUM' 0.;
  72. TCMTD2  = BORN (CHTKSC2 - 378.15) 'SCAL' 'MINIMUM' 0.;
  73. FIRE1   = (a3*(TCMTD1**3.)) + (a2*(TCMTD1**2.)) +  (a1*TCMTD1);
  74. FIRE2   = (a3*(TCMTD2**3.)) + (a2*(TCMTD2**2.)) +  (a1*TCMTD2);
  75. DDCHi    = BORN (FIRE2 - FIRE1) 'SCAL' 'MINIMUM' 0.;
  76. *
  77. DCHi     = FIRE1 + (THETA_W * DDCHi);
  78. *--------------------------------------------------------------------*
  79.  
  80. DCH_Wi   =  'CHAN' 'CHAM' DCHi  MOD_THM;
  81. DDCH_Wi  =  'CHAN' 'CHAM' DDCHi  MOD_THM;
  82.  
  83. 'SI' (indz EGA 1);
  84.     DCH_W  = DCH_Wi;
  85.     DDCH_W = DDCH_Wi;
  86. 'SINON';
  87.     DCH_W  = DCH_W  + DCH_Wi;
  88.     DDCH_W = DDCH_W + DDCH_Wi;
  89. 'FINSI';
  90.  
  91.  indz = indz + 2;
  92. 'FIN' BOUCLE_Z;
  93.  
  94. FINPROC DDCH_W DCH_W;
  95. ***********************************************************************
  96. * FIN PROCEDURE @CH_DAMA
  97. ***********************************************************************
  98.  
  99.  

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