Télécharger la plaquette Cast3M


fr en

Télécharger htc_per.procedur

Retour à la liste

Numérotation des lignes : 

  1. * HTC_PER   PROCEDUR  AM        11/09/13    21:15:02     7126           
  2. *-----------------------------------------------------------------------
  3. *23456789012345678901234567890123456789012345678901234567890123456789012
  4. *-----------------------------------------------------------------------
  5. *------------------------------------------------
  6. * HTC_PER Permeability of heated concrete
  7. * (after BAZANT and THONGUTHAI)
  8. * Date: 16-03-98
  9. *------------------------------------------------
  10. 'DEBP' HTC_PER T0*CHPOINT HR*CHPOINT TAV1*TABLE ;
  11. *------------------------------------------------
  12. *
  13.  FLG1= TAV1.'FLG1';
  14. 'SI' FLG1; 'MESS' 'PROCEDURE-HTC_PER'; 'FINSI'; 
  15.  TREF=25.;
  16.  GR=9.80665;
  17.  HC=.75;
  18.  A0=TAV1. 'A0';  
  19.  TRK = 273.15 + TREF;
  20.  F2T95='EXP' (2700. * ((368.15 - TRK)/(TRK * 368.15))); 
  21.  A0P=A0*F2T95;
  22. * 
  23.  RESO1= T0 'MASQUE' 'EGINFE'    95.;    
  24.  RESO2= T0 'MASQUE' 'SUPERIEUR' 95.;
  25.  RESO3= HR 'MASQUE' 'EGSUPE'    1. ;
  26.  RESO4= HR 'MASQUE' 'INFERIEUR' 1. ;
  27.  RESO5= T0 'MASQUE' 'EGSUPE'    TREF;    
  28. *
  29. 'SI' (('MINI' RESO1) > .9);
  30. *  T0 TUTTO UGUALE SOTTO 95
  31. *  ----> f2t          
  32.   TTT = 273.15 + T0 ;                                             
  33.   F2T = 'EXP' (RESO5 * 2700. * ((1. / TRK)-(TTT**-1)));
  34. *  CASO MISTO T SOTTO 95 E  H SOTTO E SOPRA 1
  35. *f1h dans le cas h >= 1
  36.   F1HS  = RESO3;
  37. *f1h dans le cas h <  1
  38.   ALFAT = (1. + ((19./(95.-TREF))*(95. - T0)))**-1;  
  39.   F1HI  = RESO4 * (ALFAT + ((1.-ALFAT)/(1.+(((1.- HR)/(1.-HC))**4))));
  40.   F1H   = F1HS + F1HI;
  41.   AGDG=A0*F1H*F2T/GR; 
  42. *
  43.  'DETR' TTT  ; 'DETR' ALFAT; 'DETR' F2T  ;  
  44.  'DETR' F1H  ; 'DETR' F1HS ; 'DETR' F1HI ; 
  45.  'DETR' RESO1; 'DETR' RESO2; 'DETR' RESO3;
  46.  'DETR' RESO4; 'DETR' RESO5;
  47. 'SINON';
  48.  'SI' (('MINI' RESO2) > .9);
  49. * T0 TUTTO SOPRA 95
  50.    T095=  T0 - 95. ;
  51.    F3T = 'EXP' (T095/(0.881 + (0.214*T095)));
  52.    AGDG=A0P*F3T/GR; 
  53.   'DETR' T095; 'DETR' F3T;
  54.   'DETR' RESO1; 'DETR' RESO2; 'DETR' RESO3;
  55.   'DETR' RESO4; 'DETR' RESO5;
  56.  'SINON'; 
  57. *T0 SOPRA E SOTTO 95
  58.    T095= RESO2 * ( T0 - 95.);
  59.    F3T = 'EXP' (T095/(0.881 + (0.214*T095)));
  60. *f1h dans le cas h >= 1
  61.    F1HS  = RESO3;
  62. *f1h dans le cas h <  1
  63.    ALFAT= (1. + ((19./(95.-TREF))*(95. - T0)))**-1;  
  64.    F1HI =RESO4 * (ALFAT + ((1.-ALFAT)/(1.+(((1.- HR)/(1.-HC))**4))));
  65.    F1H  = F1HS + F1HI;
  66. *  ----> f2t          
  67.    TTT = 273.15 + T0 ;                                 
  68.    F2T = 'EXP' (RESO5 * 2700. * ((1. / TRK)-(TTT**-1)));
  69.    AMA95=RESO2*(A0P*F3T);
  70.    AIE95=RESO1*(A0*F1H*F2T);
  71.    AGDG = (AMA95 + AIE95)/ GR;  
  72.   'DETR' T095; 'DETR' F3T ;
  73.   'DETR' TTT  ; 'DETR' ALFAT; 'DETR' F2T  ;  
  74.   'DETR' F1H  ; 'DETR' F1HS ; 'DETR' F1HI ;
  75.   'DETR' AMA95; 'DETR' AIE95;
  76.   'DETR' RESO1; 'DETR' RESO2; 'DETR' RESO3;
  77.   'DETR' RESO4; 'DETR' RESO5;
  78.  'FINSI';
  79. 'FINSI';                                        
  80. *-----------------------------------------------------------------------
  81. 'FINPROC' AGDG;
  82.  
  83.  

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