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 par une source de chaleur variable *
  5. * 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 variables *
  10. * - la definition du chargement (convection + source variable) *
  11. * - la resolution avec PASAPAS *
  12. * - le post traitement *
  13. ************************************************************************
  14.  
  15.  
  16. ** Options generales
  17. OPTI 'DIME' 2 'ELEM' 'QUA8' 'ECHO' 0 ;
  18. ITRAC = FAUX ;
  19.  
  20.  
  21. ** Parametres geometriques, materiaux, conditions aux limites
  22. TPSMAX = 1.E5 ;
  23. L = 1. ;
  24. R = 0.25 ;
  25. CP = 128. ;
  26. RHO = 19300. ;
  27. LT = PROG 0. 'PAS' 10. 2000. ;
  28. LL = 200. + (0.3 * LT) ;
  29. LAMBDA = EVOL 'MANU' 'T' LT 'K' LL ;
  30. LQ = 4.E6 * (EXP (-1. * (((LT - 1000.) / 700.) ** 2))) ;
  31. VSOUR = EVOL 'MANU' 'T' LT 'Q' LQ ;
  32. HCONV = EVOL 'MANU' 'TEMP' (PROG 0. TPSMAX) 'H' (PROG 10. 200.) ;
  33. T_INI = 20. ;
  34. T_EXT = 20. ;
  35.  
  36.  
  37. ** Maillage
  38. OPTI 'DENS' (R/5.) ;
  39. P0 = 0. 0. ;
  40. P1 = L 0. ;
  41. LBAS = DROI P0 P1 ;
  42. S1 = LBAS TRAN (0. L) ;
  43. CHX = S1 COOR 1 ;
  44. PTX = CHX POIN 'EGINFE' (R + 1.E-8) ;
  45. EL1 = S1 ELEM 'APPU' 'STRI' PTX ;
  46. CHY = EL1 COOR 2 ;
  47. PTY = CHY POIN 'EGINFE' (R + 1.E-8) ;
  48. MAIL1 = EL1 ELEM 'APPU' 'STRI' PTY ;
  49. MAIL2 = DIFF S1 MAIL1 ;
  50. MAILT = MAIL1 ET MAIL2 ;
  51. LDROITE = S1 COTE 2 ;
  52. LHAUT = S1 COTE 3 ;
  53. LEXT = LDROITE ET LHAUT ;
  54. P2 = MAILT 'POIN' 'PROC' (L L) ;
  55. SI ITRAC ;
  56. TRAC 'FACE' ((MAIL1 COUL 'ROUG') ET (MAIL2 COUL 'BLEU')) 'TITR' 'Maillage (zone source en rouge)' ;
  57. FINSI ;
  58.  
  59.  
  60. ** Modeles et materiaux
  61. MOD_THER = MODE MAILT 'THERMIQUE' 'CONDUCTION' ;
  62. MAT_THER = MATE MOD_THER 'K' LAMBDA 'C' CP 'RHO' RHO ;
  63. * sous-modele reduit au maillage de la source (pour l'operateur SOUR)
  64. MOD1 = REDU MOD_THER MAIL1 ;
  65.  
  66. MOD_CONV = MODE LEXT 'THERMIQUE' 'CONVECTION' ;
  67. MAT_CONV = MATE MOD_CONV 'H' HCONV ;
  68.  
  69.  
  70. ** Chargements thermiques
  71. EV1 = EVOL 'MANU' (PROG 0. TPSMAX) (PROG 1. 1.) ;
  72. CHTCONV = MANU 'CHPO' LEXT 1 'T' T_EXT ;
  73. CHA1 = CHAR 'TECO' CHTCONV EV1 ;
  74.  
  75.  
  76. ** Champ de temperature initial
  77. TEM0 = MANU 'CHPO' MAILT 1 'T' T_INI ;
  78.  
  79.  
  80. **Resolution avec PASAPAS
  81. TAB1 = TABL ;
  82. TAB1 . 'MODELE' = MOD_THER ET MOD_CONV ;
  83. TAB1 . 'CARACTERISTIQUES' = MAT_THER ET MAT_CONV ;
  84. TAB1 . 'TEMPERATURES' = TABL ;
  85. TAB1 . 'TEMPERATURES' . 0 = TEM0 ;
  86. TAB1 . 'CHARGEMENT' = CHA1 ;
  87. TAB1 . 'TEMPS_CALCULES' = PROG 0. 'PAS' (TPSMAX / 50.) TPSMAX ;
  88. TAB1 . 'PROCEDURE_CHARTHER' = VRAI ;
  89. TAB1 . 'AMOI' = TABL ;
  90. TAB1 . 'AMOI' . 'MAIL' = MAIL1 ;
  91. TAB1 . 'AMOI' . 'MODE' = MOD1 ;
  92. TAB1 . 'AMOI' . 'EVOL' = VSOUR ;
  93.  
  94.  
  95. ** Procedure CHARTHER
  96. DEBP CHARTHER T1*'TABLE' TPS1*'FLOTTANT' ;
  97. MA1 = T1 . 'AMOI' . 'MAIL' ;
  98. MO1 = T1 . 'AMOI' . 'MODE' ;
  99. EV1 = T1 . 'AMOI' . 'EVOL' ;
  100. * reduction du champ de temperatures au maillage de la source
  101. *** CHT1 = T1 . 'ESTIMATION' . 'TEMPERATURES' ;
  102. CHT1 = T1 . 'WTABLE' . 'THER_COURANT' ;
  103. CHT2 = REDU CHT1 MA1 ;
  104. * interpolation du champ de source a partir du champ de temperature
  105. CHS = IPOL CHT2 EV1 ;
  106. * creation du second membre
  107. CHFL1 = SOUR MO1 CHS ;
  108. * sortie du second membre
  109. T2 = TABL ;
  110. T2 . 'ADDI_SECOND' = CHFL1 ;
  111. FINP T2 ;
  112.  
  113.  
  114. ** Appel à la procedure PASAPAS
  115. PASAPAS TAB1 ;
  116.  
  117.  
  118. ** Post-traitement
  119. NBOU1 = DIME (TAB1 . 'TEMPS_CALCULES') ;
  120. T_MAX = 1800. ;
  121. LIS1 = PROG T_INI 'PAS' ((T_MAX - T_INI) / 20.) T_MAX ;
  122. REPE BOU1 NBOU1 ;
  123. I = &BOU1 - 1 ;
  124. * temps
  125. TEMPS_I = TAB1 . 'TEMPS' . I ;
  126. P100 = ENTI (100. * TEMPS_I / TPSMAX) ;
  127. * temperatures au temps i
  128. TEMP_I = TAB1 . 'TEMPERATURES' . I ;
  129. SI ITRAC ;
  130. TRAC TEMP_I MAILT LIS1 'TITR' (CHAI 'Temperature au temps' TEMPS_I ' (' P100 '%)') ;
  131. FINSI ;
  132. FIN BOU1 ;
  133. * courbes d'evolution temporelles des temperatures aux coins de la piece
  134. EVTP0 = EVOL 'ROUG' 'TEMP' TAB1 'TEMPERATURES' 'T' P0 ;
  135. EVTP2 = EVOL 'BLEU' 'TEMP' TAB1 'TEMPERATURES' 'T' P2 ;
  136. TL = TABL ;
  137. TL . 'TITRE' = TABL ;
  138. TL . 'TITRE' . 1 = 'Coin inf. gauche' ;
  139. TL . 'TITRE' . 2 = 'Coin sup. droite' ;
  140. SI ITRAC ;
  141. DESS (EVTP0 ET EVTP2) 'LEGE' 'SE' TL 'TITR' 'Temperature aux coins vs Temps' 'YBOR' 0. 2000. ;
  142. FINSI ;
  143.  
  144. FIN ;
  145.  
  146.  
  147.  

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