Télécharger drvdt.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : drvdt.dgibi
  2. * -------------- VARI option DRVDT -------------------------------
  3. * Test de l'opérateur VARI DRVDT(P,T)
  4. * dervivee partielle de la densite de la vapeur d'eau par
  5. * rapport a la temperature
  6. * Les données sont un FLOTTANT, un LISTREEL ou un CHPO
  7. * notation : DT**** derivee partielle
  8. * notation : DC**** derivee calculee par accroissement
  9.  
  10. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  11. 'OPTI' 'ECHO' 0 ;
  12. DELTAT=0.1D0 ;
  13. DELTAP=1000D0;
  14. *---------------- FLOTTANT----------------------------------------
  15. P1 = 1.0133D5 ;
  16. T1=373.15 ;
  17. P2 = 5.021D5 ;
  18. T2=425.15;
  19. P3 = 10.027D5;
  20. T3=453.15 ;
  21. DTPOINT1 = VARI DRVDT P1 T1 ;
  22. DTPOINT2 = VARI DRVDT P2 T2 ;
  23. DTPOINT3 = VARI DRVDT P3 T3 ;
  24. DCPOINT4=((VARI ROVAP (P1) (T1+DELTAT))-(VARI ROVAP P1 T1))/DELTAT;
  25. DCPOINT5=((VARI ROVAP (P2) (T2+DELTAT))-(VARI ROVAP P2 T2))/DELTAT ;
  26. DCPOINT6=((VARI ROVAP (P3) (T3+DELTAT))-(VARI ROVAP P3 T3))/DELTAT;
  27. *----------------- LISTREEL----------------------------------------
  28. T1 = 'PROG' 373.15D0 425.15D0 453.15D0 ;
  29. P1 = 'PROG' 1.0133D5 5.021D5 10.027D5 ;
  30. DTLIST1 = VARI DRVDT P1 T1 ;
  31. DELTAT2 = 'PROG' 0.1D0 0.1D0 0.1D0 ;
  32. DELTAP2 = 'PROG' 100D0 100D0 100D0 ;
  33. DCLIST1= ((VARI ROVAP (P1) (T1+DELTAT2))-(VARI ROVAP P1 T1))/DELTAT2;
  34. *-----------------CHAMPOINT-----------------------------------------
  35. P11 = 1.0133D5 0. ;
  36. P22 = 5.021D5 0. ;
  37. P33 = 10.027D5 0. ;
  38. P11P33 = P11 'DROI' 1 P22 'DROI' 1 P33 ;
  39. P11P33_1 = CHAN 'POI1' P11P33 ;
  40. CT1= MANU 'CHPO' P11P33_1 'REAL' T1 ;
  41. CP1= MANU 'CHPO' P11P33_1 'REAL' P1;
  42. CT2 = MANU 'CHPO' P11P33_1 'REAL' (T1+DELTAT2);
  43. DCCHPO1= ((VARI ROVAP CP1 CT2)-(VARI ROVAP CP1 CT1))/DELTAT;
  44. DTCHPO1 = VARI DRVDT CP1 CT1 ;
  45. *--------------------------------------------------------------------
  46. *----RESULTATS ET STOCKAGES DES DONNEES------------------
  47. *----------extraction des listreal-----------------------
  48. LPOINT1= 'EXTR' DTLIST1 1 ;
  49. LPOINT2 = 'EXTR' DTLIST1 2 ;
  50. LPOINT3= 'EXTR' DTLIST1 3 ;
  51. LPOINT4 = 'EXTR' DCLIST1 1 ;
  52. LPOINT5 = 'EXTR' DCLIST1 2 ;
  53. LPOINT6= 'EXTR' DCLIST1 3 ;
  54. *----------extraction des champoints---------------------
  55. CPOINT4 = 'EXTR' DCCHPO1 'REAL' P11 ;
  56. CPOINT5 = 'EXTR' DCCHPO1 'REAL' P22 ;
  57. CPOINT6 = 'EXTR' DCCHPO1 'REAL' P33 ;
  58. CPOINT1 = 'EXTR' DTCHPO1 'REAL' P11 ;
  59. CPOINT2 = 'EXTR' DTCHPO1 'REAL' P22 ;
  60. CPOINT3= 'EXTR' DTCHPO1 'REAL' P33 ;
  61. *-----------------------------------------------------------------------
  62. *----CALCULS VALEURS REFERENCES XREFi CASTEM----------------------------
  63. *-----------------------------------------------------------------------
  64. XREF1=(DCPOINT4+LPOINT4+CPOINT4)/3;
  65. XREF2=(DCPOINT5+LPOINT5+CPOINT5)/3;
  66. XREF3=(DCPOINT6+LPOINT6+CPOINT6)/3;
  67. *------------CALCULS DES ERREURS---------------------------------------
  68. CTRL1 = DTPOINT1+ LPOINT1 + CPOINT1 / 3.d0 - DTPOINT1 / DTPOINT1 ;
  69. CTRL2 = DTPOINT2 + LPOINT2 + CPOINT2 / 3.d0 - DTPOINT2 / DTPOINT2 ;
  70. CTRL3 = DTPOINT3 + LPOINT3 + CPOINT3 / 3.d0 - DTPOINT3 / DTPOINT2 ;
  71. CTRL4 = CTRL1 + CTRL2 + CTRL3 ;
  72. CTRL5 = DTPOINT1 + LPOINT1 + CPOINT1 / 3.d0 - XREF1 / XREF1 ;
  73. CTRL6 = DTPOINT2 + LPOINT2 + CPOINT2 / 3.d0 - XREF2 / XREF2 ;
  74. CTRL7 = DTPOINT3 + LPOINT3 + CPOINT3 / 3.d0 - XREF3 / XREF3 ;
  75. CTRL8 = CTRL5 + CTRL6 + CTRL7 ;
  76. *-------------------------------------------> Affichage
  77. 'MESS' ' ' ;
  78. 'MESS' ' ' ;
  79. 'MESS' 'Pt 1 VARI DRVDT ----->' DTPOINT1 LPOINT1 CPOINT1 XREF1 ;
  80. 'MESS' 'Pt 2 VARI DRVDT ----->' DTPOINT2 LPOINT2 CPOINT2 XREF2 ;
  81. 'MESS' 'Pt 3 VARI DRVDT ----->' DTPOINT3 LPOINT3 CPOINT3 XREF3 ;
  82. 'MESS' ' ' ;
  83. 'MESS' 'Comparaison calculs rel.----->' CTRL4 ;
  84. 'MESS' 'Erreur rel. VDI ----->' CTRL8 ;
  85. 'MESS' ' ' ;
  86. 'MESS' ' ' ;
  87. 'MESS' ' ' ;
  88. 'MESS' ' ' ;
  89. 'MESS' ' ' ;
  90. *------------------------------------------> Compte-rendu et sortie
  91. EPS4 = 1.E-14 ;
  92. EPS8 = 2.E-3 ;
  93. CTRL4 = ABS CTRL4 ;
  94. CTRL8 = ABS CTRL8 ;
  95. LOG4 = CTRL4 > EPS4 ;
  96. LOG8 = CTRL8 > EPS8 ;
  97. L0 = LOG4 'OU' LOG8 ;
  98. 'SI' L0 ;
  99. 'ERREUR' 5 ;
  100. 'SINON' ;
  101. 'ERREUR' 0 ;
  102. 'FINSI' ;
  103. 'FIN' ;
  104.  
  105.  
  106.  
  107.  
  108.  

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