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