Télécharger la plaquette Cast3M


fr en

Télécharger Test_celsius.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  Test_celsius.dgibi
  2.  
  3. *******************************************
  4. *******************************************
  5.  
  6. *******************************************************
  7. *           Test Test_celsius.dgibi: Jeux de données  *
  8. *           ---------------------------------         *
  9. *                                                     *
  10. *  - Suite a une anomalie, ce cas-test permet de      *
  11. *    tester que les reactions rendues par PASAPAS sont*
  12. *    correctes en THERMIQUE suivant si l'indice       *
  13. *   'CESLIUS' est defini a VRAI ou FAUX dans la TABLE *
  14. *******************************************************
  15.  
  16. OPTI DIME 2 ELEM QUA8 ;
  17. OPTI TRAC PSC;
  18.  
  19. T_IMP = 100.;
  20. T_INI_AIR = 25.;
  21. TPSFIN = 40000.;
  22. N_PAS = 100;
  23.  
  24.  
  25. P0 = 0. 0.;
  26. P1 = 0.1 0.;
  27.  
  28.  
  29.  
  30. CERC1 = CERC 50 'ROTA' 360. P1 P0;
  31. ELIM CERC1 1E-10;
  32.  
  33.  
  34. SURF_CERC = SURF CERC1;
  35.  
  36. MOD_COND = MODE SURF_CERC 'THERMIQUE';
  37. MAT_COND = MATE MOD_COND 'K' 10. 'C' 1000. 'RHO' 1000. 'TINI' T_INI_AIR;
  38.  
  39. *On impose la temperature au point de chauffe
  40. PIMP = SURF_CERC POIN 'PROC' P0;
  41.  
  42. BL_PANIER = BLOQ PIMP 'T';
  43. FL_PANIER = DEPI BL_PANIER T_IMP;
  44.  
  45.  
  46.  
  47. LISTE1  = PROG 0. TPSFIN;
  48. LISTE2  = PROG 1. 1.;
  49. EVT1    = EVOL 'MANU' LISTE1 LISTE2;
  50. CHATIMP = CHAR 'TIMP' FL_PANIER EVT1;
  51.  
  52. TPS_CALC = PROG 0. PAS (TPSFIN/N_PAS) TPSFIN ;
  53.  
  54.  
  55.  
  56. TRAC SURF_CERC;
  57.  
  58.  
  59. TAPAS                      = TABL;
  60. TAPAS.CELSIUS              = VRAI ;
  61. TAPAS.PRECISION            = 1E-5;
  62. TAPAS.MODELE               = MOD_COND;
  63. TAPAS.CARACTERISTIQUES     = MAT_COND;
  64. TAPAS.BLOCAGES_THERMIQUES  = BL_PANIER ;
  65. TAPAS.CHARGEMENT           = CHATIMP ;
  66. TAPAS.TEMPS_CALCULES       = TPS_CALC ;
  67.  
  68. PASAPAS TAPAS;
  69.  
  70.  
  71. MAIL_CAL ='EXTR' (TAPAS.MODELE) 'MAIL';
  72.  
  73.  
  74.  
  75. NB_PAS   ='DIME' TAPAS.'TEMPS';
  76.  
  77.  
  78.  
  79. CTK1 = (TAPAS.TEMPERATURES. 0)  ;
  80. CTC1 = CTK1  ;
  81. titre  'Temperature rayonnement t=0' ;
  82. TRAC CTC1 MAIL_CAL ;
  83. CTK1 = (TAPAS.TEMPERATURES.(NB_PAS - 1))   ;
  84. CTC1 = CTK1  ;
  85. titre  'Temperature rayonnement t=t_{final}' ;
  86. TRAC CTC1 MAIL_CAL ;
  87.  
  88.  
  89. OPTI EPTR 5;
  90.  
  91. * suivi de la temperature
  92. TEMP1 = PROG ;TEMP2 = PROG ;LISTPUI = PROG;
  93. ITER  = 0 ;
  94. NBT1  = NB_PAS ;
  95. REPETER BCL1  NBT1 ;
  96.   TP1    = REDU (EXCO 'T' (TAPAS.'TEMPERATURES'.ITER ) 'T' ) SURF_CERC   ;
  97.   T1     = MINI ( TP1 ) ;
  98.   TEMP1  = TEMP1 'ET' T1 ;
  99.   TEMP2  = TEMP2 ET T_IMP;
  100.   CH_PUI = REAC BL_PANIER (TAPAS.'TEMPERATURES'.ITER );
  101.   VAL_PUI = MAXI (RESU CH_PUI);
  102.   LISTPUI = LISTPUI ET VAL_PUI;
  103. ITER = ITER + 1;
  104. FIN BCL1 ;
  105.  
  106. TITRE 'EVOLUTION DES TEMPERATURES EN CELSIUS' ;
  107. EV1  = EVOL VERT MANU 'TEMPS (s)' TPS_CALC 'T (C)' TEMP1 ;
  108. EV2  = EVOL ROUG MANU 'TEMPS (s)' TPS_CALC 'T (C)' TEMP2 ;
  109.  
  110.  
  111. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  112. TABDES.'TITRE'. 1  = 'Temperature minimale du cercle' ;
  113. TABDES.'TITRE'. 2  = 'Temperature imposee' ;
  114. DESS (EV1 ET EV2) TABDES 'MIMA' 'LEGE';
  115.  
  116.  
  117. TITRE 'Evolution de la puissance avec temperature en celsius' ;
  118. EV3 = EVOL ROUG MANU 'TEMPS (s)' TPS_CALC 'Q (W)' LISTPUI ;
  119. EV4 = EVOL ROUG MANU 'TEMPS (s)' TPS_CALC 'E (J)' LISTPUI ;
  120.  
  121.  
  122. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  123. TABDES.'TITRE'. 1 = 'Puissance de chauffe' ;
  124.  
  125. DESS (EV3) TABDES 'MIMA' 'LEGE' ;
  126.  
  127.  
  128. EV5 = PRIM EV4;
  129. TITRE 'Evolution de energie avec temperature en celsius' ;
  130. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  131. TABDES.'TITRE'. 1 = 'Energie totale du systeme' ;
  132. DESS(EV5) TABDES 'MIMA' 'LEGE';
  133.  
  134.  
  135. E_four_vrai = INTG EV3;
  136. LIST E_four_vrai;
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147. *** Température en kelvin
  148.  
  149. TAPAS                      = TABL;
  150. TAPAS.CELSIUS              = FAUX ;
  151. TAPAS.PRECISION            = 1E-5;
  152. TAPAS.MODELE               = MOD_COND;
  153. TAPAS.CARACTERISTIQUES     = MAT_COND;
  154. TAPAS.BLOCAGES_THERMIQUES  = BL_PANIER ;
  155. TAPAS.CHARGEMENT           = CHATIMP ;
  156. TAPAS.TEMPS_CALCULES       = TPS_CALC ;
  157.  
  158. PASAPAS TAPAS;
  159.  
  160.  
  161. MAIL_CAL ='EXTR' (TAPAS.MODELE) 'MAIL';
  162.  
  163.  
  164.  
  165. NB_PAS   ='DIME' TAPAS.'TEMPS';
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172. * suivi de la temperature
  173. TEMP1 = PROG ;TEMP2 = PROG ;LISTPUI = PROG;
  174. ITER  = 0 ;
  175. NBT1  = NB_PAS ;
  176. REPETER BCL1  NBT1 ;
  177.   TP1    = REDU (EXCO 'T' (TAPAS.'TEMPERATURES'.ITER ) 'T' ) SURF_CERC   ;
  178.   T1     = MINI ( TP1 ) ;
  179.   TEMP1  = TEMP1 'ET' T1 ;
  180.   TEMP2  = TEMP2 ET T_IMP;
  181.   CH_PUI = REAC BL_PANIER (TAPAS.'TEMPERATURES'.ITER );
  182.   VAL_PUI = MAXI (RESU CH_PUI);
  183.   LISTPUI = LISTPUI ET VAL_PUI;
  184. ITER = ITER + 1;
  185. FIN BCL1 ;
  186.  
  187. TITRE ' EVOLUTION DES TEMPERATURES EN KELVIN' ;
  188. EV11  = EVOL VERT MANU 'TEMPS (s)' TPS_CALC 'T (K)' TEMP1 ;
  189. EV21  = EVOL ROUG MANU 'TEMPS (s)' TPS_CALC 'T (K)' TEMP2 ;
  190.  
  191.  
  192. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  193. TABDES.'TITRE'. 1  = 'Temperature minimale du cercle' ;
  194. TABDES.'TITRE'. 2  = 'Temperature imposee' ;
  195. DESS (EV11 ET EV21) TABDES 'MIMA' 'LEGE';
  196.  
  197.  
  198. TITRE 'Evolution de la puissance en fonction du temps avec temperature en kelvin' ;
  199. EV31 = EVOL ROUG MANU 'TEMPS (s)' TPS_CALC 'Q (W)' LISTPUI ;
  200. EV41 = EVOL ROUG MANU 'TEMPS (s)' TPS_CALC 'E (J)' LISTPUI ;
  201.  
  202.  
  203. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  204. TABDES.'TITRE'. 1 = 'Puissance de chauffe' ;
  205.  
  206. DESS (EV31) TABDES 'MIMA' 'LEGE' ;
  207.  
  208.  
  209. EV51 = PRIM EV41;
  210. TITRE 'Evolution de energie en fonction du temps avec temperature en kelvin' ;
  211. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  212. TABDES.'TITRE'. 1 = 'Energie du totale du systeme' ;
  213. DESS(EV51) TABDES 'MIMA' 'LEGE';
  214.  
  215.  
  216. E_four_faux = INTG EV31;
  217. LIST E_four_faux;
  218.  
  219. E_four_diff = E_four_faux - E_four_vrai;
  220. EV_pui_diff = EV31 - EV3;
  221.  
  222. TITRE 'Evolution de la difference de puissance entre les cas temperature en kelvin et en celsius' ;
  223. TABDES = TABLE; TABDES.'TITRE' = TABLE ;
  224. TABDES.'TITRE'. 1 = 'Difference de puissance' ;
  225. DESS EV_pui_diff TABDES 'MIMA' 'LEGE';
  226.  
  227. Critere='MAXI' 'ABS' ('EXTR' EV_pui_diff 'ORDO' 1);
  228.  
  229.   MESS 'Critere :' Critere ;
  230.  
  231. SI ((ABS Critere) >EG 1E-10 );
  232.   MESS 'Erreur de puissance';
  233.   ERRE 21 ;
  234. FINSI;
  235.  
  236.  
  237.  
  238. FIN; 
  239.  
  240.  

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