Télécharger la plaquette Cast3M


fr en

Télécharger calcp.procedur

Retour à la liste

Numérotation des lignes : 

  1. * CALCP     PROCEDUR  MAGN      99/12/09    21:15:07     3728           
  2. DEBPROC CALCP ;
  3. 'ARGU' T/'CHPOINT'
  4.        TF/'FLOTTANT'
  5.        TL/'LISTREEL' ;
  6. *-------------------------------------------------*
  7. * calcul de la capacite calorifique a pression    *
  8. * constante pour les incondensables en fonction   *
  9. * de la temperature  (J/kg/K)                     *
  10. *                                                 *
  11. * Entree T en Celsius                             *
  12. * Sortie Cp en J/kg/K                             *
  13. *                                                 *
  14. * PS : le resultat est du meme type que l'entree  *
  15. *-------------------------------------------------*
  16. A0 = 14.346 ;
  17. A1 = 8.16D-4 ;
  18. A2 = 4.35D-7 ;
  19. B0 = 5.193 ;
  20. B1 = 0.0 ;
  21. B2 = 0.0 ;
  22. C0 = 0.917 ;
  23. C1 = 1.404D-4 ;
  24. C2 = 3.0D-8 ;
  25. D0 = 1.0309 ;
  26. D1 = 2.044D-4 ;
  27. D2 = -3.3D-8 ;
  28. E0 = 0.836  ;
  29. E1 = 7.432D-4 ;
  30. E2 = -2.511D-7 ;
  31. F0 = 1.0302 ;
  32. F1 = 2.396D-4 ;
  33. F2 = -4.950D-8 ;
  34.  
  35. 'SI' ( 'EXIST' T ) ;
  36.     Cph2  = ( ((A2 * T) * T ) + ( T * A1 ) + A0 ) * 1.D3        ;
  37.     Cphe  = ( ((B2 * T) * T ) + ( T * B1 ) + B0 ) * 1.D3        ;
  38.     Cpo2  = ( ((C2 * T) * T ) + ( T * C1 ) + C0 ) * 1.D3        ;
  39.     Cpn2  = ( ((D2 * T) * T ) + ( T * D1 ) + D0 ) * 1.D3        ;
  40.     Cpco2 = ( ((E2 * T) * T ) + ( T * E1 ) + E0 ) * 1.D3        ;
  41.     Cpco  = ( ((F2 * T) * T ) + ( T * F1 ) + F0 ) * 1.D3        ;
  42.     Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  43.     'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  44.     'QUITTER' CALCP ;
  45. 'FINSI' ;
  46. 'SI' ( 'EXIST' TF) ;
  47.     Cph2  = ( ((A2 * TF) * TF ) + ( TF * A1 ) + A0 ) * 1.D3        ;
  48.     Cphe  = ( ((B2 * TF) * TF ) + ( TF * B1 ) + B0 ) * 1.D3        ;
  49.     Cpo2  = ( ((C2 * TF) * TF ) + ( TF * C1 ) + C0 ) * 1.D3        ;
  50.     Cpn2  = ( ((D2 * TF) * TF ) + ( TF * D1 ) + D0 ) * 1.D3        ;
  51.     Cpco2 = ( ((E2 * TF) * TF ) + ( TF * E1 ) + E0 ) * 1.D3        ;
  52.     Cpco  = ( ((F2 * TF) * TF ) + ( TF * F1 ) + F0 ) * 1.D3        ;
  53.     Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  54.     'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  55.     'QUITTER' CALCP ;
  56. 'FINSI' ;
  57. 'SI' ( 'EXIST' TL) ;
  58.     T00 = TL / TL ;
  59.     A00 = A0 * T00 ;
  60.     A10 = A1 * T00 ;
  61.     A20 = A2 * T00 ;
  62.     B00 = B0 * T00 ;
  63.     B10 = B1 * T00 ;
  64.     B20 = B2 * T00 ;
  65.     C00 = C0 * T00 ;
  66.     C10 = C1 * T00 ;
  67.     C20 = C2 * T00 ;
  68.     D00 = D0 * T00 ;
  69.     D10 = D1 * T00 ;
  70.     D20 = D2 * T00 ;
  71.     E00 = E0 * T00 ;
  72.     E10 = E1 * T00 ;
  73.     E20 = E2 * T00 ;
  74.     F00 = F0 * T00 ;
  75.     F10 = F1 * T00 ;
  76.     F20 = F2 * T00 ;
  77.  
  78.     Cph2  = ( ((A20* TL) * TL ) + ( TL * A10) + A00) * 1.D3        ;
  79.     Cphe  = ( ((B20* TL) * TL ) + ( TL * B10) + B00) * 1.D3        ;
  80.     Cpo2  = ( ((C20* TL) * TL ) + ( TL * C10) + C00) * 1.D3        ;
  81.     Cpn2  = ( ((D20* TL) * TL ) + ( TL * D10) + D00) * 1.D3        ;
  82.     Cpco2 = ( ((E20* TL) * TL ) + ( TL * E10) + E00) * 1.D3        ;
  83.     Cpco  = ( ((F20* TL) * TL ) + ( TL * F10) + F00) * 1.D3        ;
  84.     Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  85.     'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  86.     'QUITTER' CALCP ;
  87.  
  88. 'FINSI' ;
  89. 'ERRE' 5 ;
  90. FINPROC                ;
  91.  
  92.  
  93.  
  94.  

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