Télécharger phase2d.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : phase2d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. graph= faux;
  5. *
  6. * Cas test : phase2d.dgibi
  7. *
  8. * barreau ayant une temperature variant de 0 à 250 jusqu'a son milieu
  9. * puis de 250 à 0°C à son extremité. On suppose une température de
  10. * fusion à 200°C, on laisse évoluer la température
  11. * tout devient solide.
  12. *
  13. *
  14.  
  15.  
  16. OPTI DIME 2 ELEM qua4 MODE plan;
  17. *
  18. *************** Temperature de changement de phase ************************
  19. Tfma = 200. ;
  20. *
  21. *************** Geometrie *************************************************
  22. *p1 = 0. 0. 0. ; p11 = (0. 0. 10.e-3); p12=(0. 0. -10.e-3);
  23. *p2 = 10.e-3 0. 0. ;
  24. *p3 = 10.e-3 10.e-3 0. ;
  25. *p4 = 0. 10.e-3 0. ;
  26. p1 = 0. 0. ; p11 = (0. 10.e-3); p12=(0. -10.e-3);
  27. p2 = 10.e-3 0. ;
  28. *p3 = 10.e-3 10.e-3 0. ;
  29. *p4 = 0. 10.e-3 0. ;
  30. n=1;
  31. l1 = p12 droi 100 p1 droi 100 p11;
  32. si graph;trac l1;finsi;
  33. su1= l1 trans n p2;
  34. *vol1= su1 volu trans n p4;
  35. *trac cach vol1;
  36. **************** Modeles et caracteristiques *********************************
  37. modth1 = MODE su1 'THERMIQUE' 'CONDUCTION' 'PHASE';
  38.  
  39. latmin= 16.25e19; latint = 16.85;latmax= 25.63e19;
  40.  
  41. LAT= 16.5e19;
  42.  
  43. LAT=1.e5;
  44. CAPCALO= 900.;
  45. matth1 = MATE modth1 'C' capcalo 'K' 100.
  46. 'RHO' 4500. 'TPHASE' Tfma 'QLAT' LAT;
  47. *
  48. *************** Pas de temps "grossier" et fin ********************************
  49. * Durée totale
  50. TT1 = 1.;
  51. * pas de temps "grossier"
  52. dt1 = 5.e-2 ;
  53. * pas de temps "fin"
  54. dt2 = 0.5e-6 ;
  55. * pas de temps "très fin" juste pour avant la disparition totale de la zone fondue
  56. dt3 = 0.5e-8 ;
  57. *
  58. * pas de temps du début "grossier" jusqu'à Temps1 et ensuite fin jusqu'à Temps2 ;
  59. *
  60. Temps1 = 0.350e-3 ;
  61. Temps2 = 0.378e-3 ;
  62. Temps3 = 0.385e-3 ;
  63. listT = PROG 0. PAS dT1 Temps1 PAS dt2 Temps2 PAS dt3 Temps3
  64. PAS dt1 TT1 ;
  65. listt= prog 0. pas dt1 tt1;
  66. n1 = DIME listT ;
  67. listS = PROG 0. pas 1. (n1-1) ;
  68. listS = (1./(n1-1.)) * listS ;
  69. EvS = EVOL MANU 'TEMP' listT 'SOUR' listS ;
  70. si graph; dess EvS TITR 'evolution chargement'; finsi;
  71. *
  72. **************** Définition du champ de temperatures initial ****************
  73. *xco = COOR 3 vol1 ;
  74. xco = COOR 2 su1 ;
  75. xco = ABS xco ;
  76. max1 = MAXI xco ;
  77. xco = max1 - xco ;
  78. xco = (1./max1) * xco ;
  79. *
  80. * eventuellement, pour un pic de temperature plus localisé...
  81. *xco = xco * xco ;
  82. *
  83. Txco = Tfma * 1.25 * xco ;
  84. Txco = CHAN COMP 'T' Txco;
  85. Txco = CHAN 'ATTRIBUT' Txco NATURE 'DIFFUS' ;
  86. *TRAC vol1 Txco 'TITR' ' Temperature initiale' ;
  87. si graph ;TRAC su1 Txco 'TITR' ' Temperature initiale' ; finsi;
  88. evch = evol chpo l1 Txco;
  89. si graph ; dess evch; finsi;
  90. ********************************************************************
  91. * chargement nul !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  92. ********************************************************************
  93. *sour0 = MANU CHPO vol1 1 'Q' 0. 'NATURE' 'DISC';
  94. sour0 = MANU CHPO su1 1 'Q' 0. 'NATURE' 'DISC';
  95. chthQ0 = CHAR 'Q' sour0 EvS ;
  96. ********************************************************************
  97. xtab = table ;
  98. *xtab.'CELSIUS'= VRAI ;
  99. xtab.'TEMPERATURES' = table ;
  100. xtab.'TEMPERATURES' . 0 = Txco ;
  101. xtab.'MODELE' = modth1 ;
  102. xtab.'CARACTERISTIQUES' = matth1 ;
  103. xtab.'CHARGEMENT' = chthQ0 ;
  104. xtab.'PROCEDURE_PERSO2' = VRAI ;
  105. DEBP PERSO2 ;
  106. chT = PECHE xtab.'TEMPERATURES' ;
  107. cht = exco cht 'T' 'T';
  108. * TRAC vol1 chT ;
  109. Tmax1 = MAXI chT ;
  110. titre ' profil de temperature le long du barreau';
  111. evch= evol chpo l1 cht ;
  112. si graph ; dess evch nclk; finsi;
  113. FINP ;
  114. xtab.'TEMPS_CALCULES' = listT ;
  115. *xtab.'TEMPS_SAUVES' = (PROG 0. TT1) ;
  116. *xtab.'PAS_AJUSTE' = vrai ;
  117. *xtab.'RELAXATION_THETA' = 1.;
  118. *xtab.'CRITERE' = 1.e-3 ;
  119. *
  120. PASAPAS xtab ;
  121.  
  122.  
  123.  
  124. * on compare energie fournie/dépensée par chaleur latente dql
  125. * à la différence d'énergie calorique % à la solution homogène 125°C
  126. * obtenue sans chaleur latente (dqt).
  127.  
  128. tder = peche xtab temperatures ;
  129. tder= tder exco 'T';
  130. evt = evol chpo tder l1;
  131. so =( ( INTG evt) extr 1 ) / 2.e-2;
  132. list so;
  133. evini= evol chpo txco l1;
  134. sini= ( ( INTG evini) extr 1 ) / 2.e-2;
  135.  
  136. dqt = capcalo * ( abs( so - sini));
  137.  
  138. *proportions de phase, initiale, prise dans le cacul
  139. pr0= xtab.proportions_phase . 0;
  140. zo= intg pr0 modth1;
  141. zz = manu chml modth1 'SCAL' 1.;
  142.  
  143. zotot= intg zz modth1;
  144. na = dime xtab.proportions_phase ;
  145. prfi= xtab . proportions_phase . (na - 1);
  146. zfi= intg prfi modth1;
  147. zmod = abs ( zfi - zo ) ;
  148.  
  149. propchzo=abs( zmod / zotot);
  150. dql = LAT * propchzo ;
  151.  
  152. mess ' dqt' dqt ' dql ' dql;
  153.  
  154. err= abs ( ( dql- dqt ) / dql);
  155. si ( err > 1.e-4) ;
  156. mess ' erreur dans phase2d.dgibi' ;
  157. erreur 5 ; finsi;
  158.  
  159. fin;
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  

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