Procedures | Cast3M

   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 :