Télécharger phase_01.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : phase_01.dgibi
  2. optio norm auto;
  3. ************************************************************************
  4. ************************************************************************
  5.  
  6. ************************************************************************
  7. * Cas test : phase01.dgibi
  8. * --------
  9. *
  10. * Categorie : Verification & Validation (Formule Analytique)
  11. * ---------
  12. *
  13. * Description :
  14. * -----------
  15. * plaque 2D sur laquelle on applique 3 modeles
  16. * 1- Diffusion de la chaleur
  17. * 2- Changement de phase THERMIQUE à 'T' = 0°C
  18. * 3- Changement de phase THERMIQUE 'T' =100°C
  19. *
  20. * Le point P1 est soumis alternativement a un flux impose positif puis
  21. * negatif.
  22. *
  23. * Les 2 phases apparaissent au chauffage et la 2eme phase disparait
  24. * au refroidissement
  25. *
  26. * Validation :
  27. * ----------
  28. * L'integrale de l'energie injectee par P1 au cours du temps doit etre
  29. * egale (a la precision machine pres) a la chaleur reguliere (RHO.CP)
  30. * plus la chaleur singuliere de changement de phase (Chaleur latente)
  31. *
  32. ************************************************************************
  33.  
  34.  
  35. GRAPH = FAUX;
  36. OPTI TRAC 'PSC';
  37. OPTI DIME 2 ELEM 'QUA4' MODE plan ;
  38.  
  39. *************** Temperature de changement de phase ************************
  40. T_ch1 = 0. ;
  41. lat1 = 1000. * 334.00D3 ;
  42. T_ch2 = 100. ;
  43. lat2 = 1000. * 2264.76D3 ;
  44. C_ch3 = 98. ;
  45. lat3 = 2.e9 ;
  46.  
  47. *************** Geometrie *************************************************
  48. P1 = 0. 0. ;
  49. P2 = 5.e-3 0. ;
  50. P3 = 5.e-3 10.e-3 ;
  51. P4 = 0. 10.e-3 ;
  52.  
  53. n = 20 ;
  54. L23 = P2 'DROI' n P3 ;
  55. L41 = P4 'DROI' n P1 ;
  56. L14 ='INVE' L41 ;
  57.  
  58. PA ='POIN' L41 'PROC' (0. 2.5e-3) ;
  59.  
  60. n = 10 ;
  61. MQUA4= ('INVE' L41) 'REGL' n L23 ;
  62. MTRI3='CHAN' 'TRI3' (L23 'REGL' n (L23 'PLUS' P2));
  63.  
  64. su1 = MQUA4 'ET' MTRI3;
  65.  
  66. MPOI1='CHAN' 'POI1' su1 ;
  67. MPOI1='POIN' MPOI1 'DROI' P1 P2 1.D-10 ;
  68. L_BAS='ELEM' ('CONT' su1) 'APPU' 'STRI' MPOI1;
  69.  
  70. **************** Modeles et caracteristiques *********************************
  71. modth1 ='MODE' su1 'THERMIQUE' 'CONDUCTION' 'CONS' 'CON1' ;
  72. modph1 ='MODE' su1 'CHANGEMENT_PHASE' 'PARFAIT' 'INCO' 'T' 'Q' 'CONS' 'CHP1' ;
  73. modph2 ='MODE' su1 'CHANGEMENT_PHASE' 'PARFAIT' 'INCO' 'T' 'Q' 'CONS' 'CHP2' ;
  74.  
  75. modtot = modth1 'ET' modph1 'ET' modph2 ;
  76. modcal = modth1 'ET' modph1 'ET' modph2 ;
  77.  
  78. rho1 = 4500. ;
  79. cp1 = 900. ;
  80. cd1 = rho1*cp1 ;
  81.  
  82. matth1 ='MATE' modth1 'C' cp1 'K' 50. 'RHO' rho1 ;
  83. matph1 ='MATE' modph1 'PRIM' T_ch1 'DUAL' lat1 ;
  84. matph2 ='MATE' modph2 'PRIM' T_ch2 'DUAL' lat2 ;
  85.  
  86. mattot = matth1 'ET' matph1 'ET' matph2 ;
  87.  
  88.  
  89. *************** Pas de temps "grossier" et fin ********************************
  90. * Durée totale
  91. TT1 = 1.2 ;
  92. TT2 = 2.0 ;
  93. dt1 = 0.1 ;
  94. *
  95. * pas de temps du début "grossier" jusqu'à Temps1 et ensuite fin jusqu'à Temps2 ;
  96. *
  97. Temps1 = 0.350e-3 ;
  98. Temps2 = 0.378e-3 ;
  99. Temps3 = 0.385e-3 ;
  100. listT = prog 0. 'PAS' dt1 TT1 'PAS' dt1 TT2;
  101. n1 = DIME listT ;
  102.  
  103. K1 = 1.D0 ;
  104. K2 =-1.D0 ;
  105. EvS = EVOL MANU 'TEMP' (PROG 0. TT1 (TT1+1.D-12) TT2) 'SOUR' (PROG K1 K1 K2 K2);
  106. si graph;
  107. dess EvS TITR 'evolution chargement';
  108. finsi;
  109.  
  110. T_ini ='MANU' 'CHPO' su1 1 'T' -20. 'NATURE' 'DIFFUS' ;
  111. P_ini ='ZERO' modtot 'VARINTER' ;
  112.  
  113. si graph ;
  114. TRAC su1 T_ini 'TITR' ' Temperature initiale' ;
  115. finsi;
  116.  
  117. * CHARGEMENT
  118. srcth1 ='MANU' 'CHPO' ('MANU' 'POI1' P1) 1 'Q' 5.D4 ;
  119. charth1 ='CHAR' 'Q' srcth1 EvS ;
  120.  
  121. chartot = charth1 ;
  122.  
  123. * PASAPAS
  124. xtab = TABL ;
  125. xtab.'TEMPERATURES' = table ;
  126. xtab.'TEMPERATURES' . 0 = T_ini ;
  127. *xtab.'CELSIUS' = VRAI ;
  128. *xtab.'PROPORTIONS_PHASE' = TABL ;
  129. *xtab.'PROPORTIONS_PHASE'. 0= P_ini ;
  130. xtab.'MODELE' = modcal ;
  131. xtab.'CARACTERISTIQUES' = mattot ;
  132. xtab.'CHARGEMENT' = chartot ;
  133. xtab.'TEMPS_CALCULES' = listT ;
  134. xtab.'PROCESSEURS' ='MONOPROCESSEUR' ;
  135.  
  136. PASAPAS xtab ;
  137.  
  138. * POST-TRAITEMENT
  139. *'OPTI' 'EPTR' 10 ;
  140. csu1 ='CONT' su1 ;
  141. MCHLAT ='EXCO' mattot 'DUAL' 'PPHA' ;
  142.  
  143. NSAUV ='DIME' xtab.'TEMPERATURES' ;
  144.  
  145. LTEMP=VIDE 'LISTREEL' ;
  146. LTPS =VIDE 'LISTREEL' ;
  147. Pmoy1=VIDE 'LISTREEL' ;
  148. Pmoy2=VIDE 'LISTREEL' ;
  149.  
  150. LISOV=PROG 0. 'PAS' 0.01 2. ;
  151.  
  152. REPE SURI NSAUV;
  153. II =&SURI ;
  154. TPSi = xtab.'TEMPS' . (II-1) ;
  155. TEMPi= xtab.'TEMPERATURES' . (II-1) ;
  156. PROPi= xtab.'PROPORTIONS_PHASE'. (II-1) ;
  157.  
  158. TEMPi='EXCO' TEMPi 'T' 'T' ;
  159. T1 ='EXTR' TEMPi 'T' PA ;
  160.  
  161. PMi1 ='INTG' modph1 PROPi ;
  162. PMi2 ='INTG' modph2 PROPi ;
  163.  
  164. LTPS =LTPS 'ET' TPSi ;
  165. LTEMP=LTEMP 'ET' T1 ;
  166. Pmoy1=Pmoy1 'ET' PMi1 ;
  167. Pmoy2=Pmoy2 'ET' PMi2 ;
  168.  
  169. PROP_SOMME =('REDU' PROPi modph1) + ('CHAN' ('REDU' PROPi modph2) 'CONS' 'CHP1');
  170. 'TRAC' PROP_SOMME modph1 'TITR' ('CHAI' 'Proportion sommee :' TPSi 's') LISOV ;
  171. *'TRAC' PROPi modph1 'TITR' ('CHAI' 'Proportion de PHASE 1:' TPSi 's') LISOV;
  172. *'TRAC' PROPi modph2 'TITR' ('CHAI' 'Proportion de PHASE 2:' TPSi 's') LISOV;
  173. *'DESS' ('EVOL' 'CHPO' TEMPi 'T' L14) ('CHAI' 'Temperature sur la ligne L1' TPSi 's');
  174. FIN SURI;
  175.  
  176. *EVO2 = 'EVOL' 'ROUG' 'MANU' 'Tps' LTPS 'TPA ' LTEMP ;
  177. *EVO3 =('EVOL' 'BLEU' 'MANU' 'Tps' LTPS 'PMO1' Pmoy1) 'ET'
  178. * ('EVOL' 'BOUT' 'MANU' 'Tps' LTPS 'PMO2' Pmoy2) ;
  179. *'DESS' EVO2 ('CHAI' 'Temperature au point PA au cours du temps');
  180. *'DESS' EVO3 ('CHAI' 'Moyennes des proportions de phase au cours du temps');
  181.  
  182. * Bilan d'energie :
  183. E_ext ='MAXI' ('RESU' ((K1*srcth1*TT1) + (K2*srcth1*(TT2-TT1)))) ;
  184. E_reg ='INTG' modth1 ('CHAN' 'CHAM' modth1 (rho1 * cp1 * (TEMPi - T_ini))) ;
  185. MCHLA1 ='REDU' modph1 MCHLAT ;
  186. DELTA1 ='REDU' modph1 (PROPi - P_ini) ;
  187. E_sing1='INTG' modph1 (MCHLA1 * DELTA1) 'PPHA' ;
  188. MCHLA2 ='REDU' modph2 MCHLAT ;
  189. DELTA2 ='REDU' modph2 (PROPi - P_ini) ;
  190. E_sing2='INTG' modph2 (MCHLA2 * DELTA2) 'PPHA' ;
  191.  
  192.  
  193. 'SI'(E_ext 'NEG' 0.D0);
  194. ERR1 =(E_ext - (E_reg + E_sing1 + E_sing2)) / E_ext ;
  195. 'SINO';
  196. ERR1 = E_ext - (E_reg + E_sing1 + E_sing2) ;
  197. 'FINS';
  198.  
  199. CHAI1 ='CHAI' E_ext ' =?' (E_reg + E_sing1 + E_sing2) ' |' ERR1 ;
  200. 'MESS' CHAI1 ;
  201. 'SI' (('ABS' ERR1) >EG 1.D-14);
  202. 'ERRE' 5 ;
  203. 'FINS';
  204.  
  205. 'FIN';
  206.  
  207.  
  208.  
  209.  
  210.  

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