Télécharger formation_pasapas_3_solution.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * FORMATION AVANCEE SUR LA PROCEDURE PASAPAS *
  3. * *
  4. * Chauffage d'une plaque avec changementde phase par une source de *
  5. * chaleur variable, dependante de la temperature) *
  6. * *
  7. * Ce fichier constitue la mise donnee solution du probleme et realise *
  8. * - le maillage *
  9. * - la definition du modele et des parametres materiaux *
  10. * - la definition du chargement (convection + source variable avec la *
  11. * procedure CHARTHER) *
  12. * - la resolution avec PASAPAS *
  13. * - le post traitement *
  14. ************************************************************************
  15.  
  16.  
  17. ** Options generales
  18. OPTI 'DIME' 2 'ELEM' 'QUA4' 'ECHO' 0 ;
  19. ITRAC = FAUX ;
  20.  
  21.  
  22. ** Parametres geometriques, materiaux, conditions aux limites
  23. L = 1. ;
  24. R = 0.25 ;
  25. LAMBDA = EVOL 'MANU' 'T' (PROG 0. 27. 727.)
  26. 'K' (PROG 401. 398. 357.) ;
  27. CP = 385. ;
  28. RHO = EVOL 'MANU' 'T' (PROG 20. 1083. 1100. 1200. 1300.)
  29. 'RHO' (PROG 8930. 7940. 7924. 7846. 7764.) ;
  30. TFUS = 1083. ;
  31. QLATFUS = 204.8E3 ;
  32. HCONV = 150. ;
  33. VSOUR = EVOL 'MANU' 'T' (PROG 20. 150. 300. 500. 700. 1000. 1500.)
  34. 'Q' (1.E6*(PROG 5. 4. 3.7 3.5 3.3 3.1 3.)) ;
  35. T_INI = 20. ;
  36. T_EXT = 20. ;
  37.  
  38.  
  39. ** Maillage
  40. OPTI 'DENS' (R/5.) ;
  41. P0 = 0. 0. ;
  42. P1 = L 0. ;
  43. LBAS = DROI P0 P1 ;
  44. S1 = LBAS TRAN (0. L) ;
  45. CHX = S1 COOR 1 ;
  46. PTX = CHX POIN 'EGINFE' (R + 1.E-8) ;
  47. EL1 = S1 ELEM 'APPU' 'STRI' PTX ;
  48. CHY = EL1 COOR 2 ;
  49. PTY = CHY POIN 'EGINFE' (R + 1.E-8) ;
  50. MAIL1 = EL1 ELEM 'APPU' 'STRI' PTY ;
  51. MAIL2 = DIFF S1 MAIL1 ;
  52. MAILT = MAIL1 ET MAIL2 ;
  53. LDROITE = S1 COTE 2 ;
  54. LHAUT = S1 COTE 3 ;
  55. LEXT = LDROITE ET LHAUT ;
  56. SI ITRAC ;
  57. TRAC 'FACE' ((MAIL1 COUL 'ROUG') ET (MAIL2 COUL 'BLEU'))
  58. 'TITR' 'Maillage, zone source en rouge' ;
  59. FINSI ;
  60.  
  61.  
  62. ** Modeles et materiaux
  63. MOD_THER = MODE MAILT 'THERMIQUE' 'PHASE' ;
  64. MAT_THER = MATE MOD_THER 'K' LAMBDA 'C' CP 'RHO' RHO
  65. 'TPHA' TFUS 'QLAT' QLATFUS ;
  66. * sous-modele reduit au maillage de la source (pour l'operateur SOUR)
  67. MOD1 = REDU MOD_THER MAIL1 ;
  68.  
  69. MOD_CONV = MODE LEXT 'THERMIQUE' 'CONVECTION' ;
  70. MAT_CONV = MATE MOD_CONV 'H' HCONV ;
  71.  
  72.  
  73. ** Chargements thermiques
  74. EV1 = EVOL 'MANU' (PROG 0. 1.E5) (PROG 1. 1.) ;
  75. CHTCONV = MANU 'CHPO' LEXT 1 'T' T_EXT ;
  76. CHA1 = CHAR 'TECO' CHTCONV EV1 ;
  77. *CHFL1 = SOUR MOD1 VSOUR MAIL1 ;
  78. *CHA2 = CHAR 'Q' CHFL1 EV1 ;
  79.  
  80.  
  81. ** Champ de temperature initial
  82. TEM0 = MANU 'CHPO' MAILT 1 'T' T_INI ;
  83.  
  84.  
  85. **Resolution avec PASAPAS
  86. LTPSCAL1 = PROG 0. 'PAS' 0.01E5 0.2E5 'PAS' 0.05E5 1.E5 ;
  87. TAB1 = TABL ;
  88. TAB1 . 'MODELE' = MOD_THER ET MOD_CONV ;
  89. TAB1 . 'CARACTERISTIQUES' = MAT_THER ET MAT_CONV ;
  90. TAB1 . 'TEMPERATURES' = TABL ;
  91. TAB1 . 'TEMPERATURES' . 0 = TEM0 ;
  92. TAB1 . 'CHARGEMENT' = CHA1 ;
  93. TAB1 . 'TEMPS_CALCULES' = LTPSCAL1 ;
  94. TAB1 . 'PROCEDURE_CHARTHER' = VRAI ;
  95. TAB1 . 'AMOI' = TABL ;
  96. TAB1 . 'AMOI' . 'MAIL' = MAIL1 ;
  97. TAB1 . 'AMOI' . 'MODE' = MOD1 ;
  98. TAB1 . 'AMOI' . 'EVOL' = VSOUR ;
  99.  
  100.  
  101. ** Procedure CHARTHER
  102. DEBP CHARTHER T1*'TABLE' TPS1*'FLOTTANT' ;
  103. MA1 = T1 . 'AMOI' . 'MAIL' ;
  104. MO1 = T1 . 'AMOI' . 'MODE' ;
  105. EV1 = T1 . 'AMOI' . 'EVOL' ;
  106. * reduction du champ de temperatures au maillage de la source
  107. *** CHT1 = T1 . 'ESTIMATION' . 'TEMPERATURES' ;
  108. CHT1 = T1 . 'WTABLE' . 'THER_COURANT' ;
  109. CHT2 = REDU CHT1 MA1 ;
  110. * interpolation du champ de source a partir du champ de temperature
  111. CHS = IPOL CHT2 EV1 ;
  112. * creation du second membre
  113. CHFL1 = SOUR MO1 CHS ;
  114. * sortie du second membre
  115. T2 = TABL ;
  116. T2 . 'ADDI_SECOND' = CHFL1 ;
  117. FINP T2 ;
  118.  
  119.  
  120. ** Appel à la procedure PASAPAS
  121. PASAPAS TAB1 ;
  122.  
  123.  
  124. ** Post-traitement
  125. NBOU1 = DIME (TAB1 . 'TEMPS_CALCULES') ;
  126. T_MAX = MAXI (TAB1 . 'TEMPERATURES' . (NBOU1-1)) ;
  127. LIS1 = PROG T_INI 'PAS' ((T_MAX-T_INI)/20.) T_MAX ;
  128. LTPS = PROG ;
  129. LPHA = PROG ;
  130. TL = TABL ;
  131. TL . 'TITRE' = TABL ;
  132. TL . 'TITRE' . 1 = 'Coin inf. gauche' ;
  133. TL . 'TITRE' . 2 = 'Coin sup. droite' ;
  134. REPE BOU1 NBOU1 ;
  135. I = &BOU1 - 1 ;
  136. TEMPS_I = TAB1 . 'TEMPS' . I ;
  137. LTPS = LTPS ET TEMPS_I ;
  138. * temperatures au temps i
  139. TEMP_I = TAB1 . 'TEMPERATURES' . I ;
  140. SI ITRAC ;
  141. * TRAC TEMP_I MAILT LIS1 'TITR' (CHAI 'Temperature au temps' TEMPS_I) ;
  142. FINSI ;
  143. * proportions de phase au temps i
  144. PHASE_I = TAB1 . 'PROPORTIONS_PHASE' . I ;
  145. LPHA = LPHA ET (INTG MOD_THER PHASE_I) ;
  146. SI ITRAC ;
  147. TRAC PHASE_I MOD_THER 'TITR'
  148. (CHAI 'Proportions de phase au temps' TEMPS_I) ;
  149. FINSI ;
  150. FIN BOU1 ;
  151. EVPHA = EVOL 'ROUG' 'MANU' 'Temps' LTPS 'Vol Phase' LPHA ;
  152. SI ITRAC ;
  153. DESS EVPHA 'TITR' 'Volume de phase liquide vs Temps' ;
  154. FINSI ;
  155. EVTP0 = EVOL 'ROUG' 'TEMP' TAB1 'TEMPERATURES' 'T' P0 ;
  156. P2 = MAILT 'POIN' 'PROC' (L L) ;
  157. EVTP2 = EVOL 'VERT' 'TEMP' TAB1 'TEMPERATURES' 'T' P2 ;
  158. SI ITRAC ;
  159. DESS (EVTP0 ET EVTP2) 'LEGE' TL
  160. 'TITR' 'Temperature aux coins vs Temps' ;
  161. FINSI ;
  162.  
  163. FIN ;
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  

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