Télécharger formation_pasapas_4_solution_bis.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * FORMATION AVANCEE SUR LA PROCEDURE PASAPAS *
  3. * *
  4. * Variation d'un jeu sous l'action d'une sollicitation thermique en *
  5. * regime transitoire *
  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. * - l'utilisation d'un modele de convection pour representer l'echange *
  11. * thermique, avec un coefficient d'echange variable dependant de la *
  12. * pression de contact *
  13. * - la definition du chargement *
  14. * - la mise a jour du coefficient d'echange selon la pression *
  15. * (avec la procedure REEV_MEC) *
  16. * - la resolution avec PASAPAS *
  17. * - le post traitement *
  18. ************************************************************************
  19.  
  20.  
  21. ** Options generales
  22. OPTI 'DIME' 2 'ELEM' 'QUA4' 'MODE' 'PLAN' 'DEFO' 'ECHO' 0 ;
  23. ITRAC = FAUX ;
  24.  
  25.  
  26. ** Parametres geometriques, materiaux et CL
  27. L = 1. ;
  28. EPAI = 0.1 ;
  29. JEU_INI = 0.01*L ;
  30. T_INI = 20. ;
  31. T1 = 500. ;
  32. T4 = 20. ;
  33. TEMPS_F = 5.E5 ;
  34. XPAS = 5.E3 ;
  35.  
  36.  
  37. ** Maillage
  38. P1 = 0. 0. ;
  39. P2 = L 0. ;
  40. L0 = DROI 10 P1 P2 ;
  41. S1 = L0 TRAN 1 (0. EPAI) ;
  42. L1 = S1 COTE 4 ;
  43. L2 = S1 COTE 2 ;
  44. S2 = S1 PLUS ((L + JEU_INI) 0.) ;
  45. L3 = S2 COTE 4 ;
  46. L4 = S2 COTE 2 ;
  47. P3 = S2 POIN 'PROC' ((L + JEU_INI) 0.) ;
  48. P4 = S2 POIN 'PROC' (((2. * L) + JEU_INI) 0.) ;
  49. MAILT = S1 ET S2 ;
  50. SI ITRAC ;
  51. TRAC MAILT 'TITR' 'Maillage' ;
  52. FINSI ;
  53.  
  54.  
  55. ** Modeles et materiaux
  56. MODM1 = MODE S1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  57. MATM1 = MATE MODM1 'ALPHA' 12.E-6 'YOUN' 200.E9 'NU' 0.3 ;
  58. MODM2 = MODE S2 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  59. MATM2 = MATE MODM2 'ALPHA' 12.E-6 'YOUN' 200.E9 'NU' 0.3 ;
  60. MODT1 = MODE S1 'THERMIQUE' 'ISOTROPE' ;
  61. MATT1 = MATE MODT1 'K' 60. 'C' 800. 'RHO' 7800. ;
  62. MODT2 = MODE S2 'THERMIQUE' 'ISOTROPE' ;
  63. MATT2 = MATE MODT2 'K' 80. 'C' 100. 'RHO' 7800. ;
  64.  
  65. MRACC = RACC (1.1 * JEU_INI) L2 L3 ;
  66. MODRACC = MODE MRACC 'THERMIQUE' 'CONVECTION' ;
  67. HVSP = EVOL 'MANU' 'PC' (PROG 0. 5.E8 5.1E8) 'H' (PROG 0. 1.E4 1.E4) ;
  68. MATRACC = MATE MODRACC 'H' HVSP ;
  69.  
  70. MOD_TH = MODT1 ET MODT2 ET MODRACC ;
  71. MAT_TH = MATT1 ET MATT2 ET MATRACC ;
  72. MOD_ME = MODM1 ET MODM2 ;
  73. MAT_ME = MATM1 ET MATM2 ;
  74.  
  75.  
  76. ** Blocages mecaniques
  77. CL1 = BLOQ 'UY' MAILT ;
  78. CL2 = BLOQ 'UX' L1 ;
  79. CL3 = BLOQ 'UX' L4 ;
  80. RE1 = RELA 'MAXI' 'UX' L2 - 'UX' L3 ;
  81. CL_ME = CL1 ET CL2 ET CL3 ET RE1 ;
  82.  
  83.  
  84. ** Blocages thermiques
  85. CL4 = BLOQ 'T' L1 ;
  86. CL5 = BLOQ 'T' L4 ;
  87. CL_TH = CL4 ET CL5 ;
  88.  
  89.  
  90. ** Chargements mecaniques et thermiques
  91. EV1 = EVOL 'MANU' (PROG 0. TEMPS_F) (PROG 1. 1.) ;
  92. F1 = DEPI RE1 JEU_INI ;
  93. CHA1 = CHAR 'DIMP' F1 EV1 ;
  94. TEM1 = DEPI CL4 T1 ;
  95. TEM2 = DEPI CL5 T4 ;
  96. CHA2 = CHAR 'TIMP' (TEM1 ET TEM2) EV1 ;
  97.  
  98.  
  99. ** Chargement decrivant le parametre PC (pression de contact)
  100. CH_PC = MANU 'CHPO' MRACC 'PC' 0. ;
  101. CHA3 = CHAR 'PC' CH_PC EV1 ;
  102.  
  103.  
  104. ** Champ de temperature initiale
  105. TEM0 = MANU 'CHPO' MAILT 1 'T' T_INI ;
  106.  
  107.  
  108. ** Procedure REEV_MEC pour actualiser le chargement decrivant la
  109. * pression de contact
  110. DEBP REEV_MEC T1*'TABLE' ;
  111. R1 = T1 . 'ESTIMATION' . 'REACTIONS' ;
  112. WT = T1 . 'WTABLE' ;
  113. TAM = T1 . 'AMOI' ;
  114. L3 = TAM . 'LIGNE_3' ;
  115. MR = TAM . 'MAIL' ;
  116. SI ((T1 . 'ESTIMATION' . 'TEMPS') NEG 0.) ;
  117. * calcul de la pression de contact
  118. VAL_PC = (MAXI 'ABS' (RESU (REDU R1 L3))) / (MESU L3) ;
  119. * chargment decrivant le parametre 'PC'
  120. CH_PC = MANU 'CHPO' MR 'PC' VAL_PC ;
  121. EV1 = TAM . 'EVOL' ;
  122. CHA3 = CHAR 'PC' CH_PC EV1 ;
  123. * on ecrase le chargement global de WTABLE
  124. WT . 'CHARGEMENT' = (TAM . 'CHAR_0') ET CHA3 ;
  125. FINS ;
  126. FINP ;
  127.  
  128.  
  129. ** Resolution avec PASAPAS
  130. TAB1 = TABL ;
  131. TAB1 . 'MODELE' = MOD_ME ET MOD_TH ;
  132. TAB1 . 'CARACTERISTIQUES' = MAT_ME ET MAT_TH ;
  133. TAB1 . 'BLOCAGES_MECANIQUES' = CL_ME ;
  134. TAB1 . 'BLOCAGES_THERMIQUES' = CL_TH ;
  135. TAB1 . 'TEMPERATURES' = TABL ;
  136. TAB1 . 'TEMPERATURES' . 0 = TEM0 ;
  137. TAB1 . 'CHARGEMENT' = CHA1 ET CHA2 ET CHA3 ;
  138. TAB1 . 'TEMPS_CALCULES' = PROG XPAS 'PAS' XPAS TEMPS_F ;
  139. TAB1 . 'CONVERGENCE_MEC_THE' = VRAI ;
  140. TAB1 . 'PROCEDURE_REEV_MEC' = VRAI ;
  141. TAB1 . 'AMOI' = TABL ;
  142. TAB1 . 'AMOI' . 'LIGNE_3' = L3 ;
  143. TAB1 . 'AMOI' . 'MAIL' = MRACC ;
  144. TAB1 . 'AMOI' . 'EVOL' = EV1 ;
  145. TAB1 . 'AMOI' . 'CHAR_0' = CHA1 ET CHA2 ;
  146. PASAPAS TAB1 ;
  147.  
  148.  
  149. ** Post-traitement
  150. LTPS = PROG ;
  151. LJEU = PROG ;
  152. LT2 = PROG ;
  153. LT3 = PROG ;
  154. LPC = PROG ;
  155. NB1 = DIME (TAB1 . 'TEMPS') ;
  156. MESS ' Temps Jeu T2 T3' ;
  157. * listage du jeu et des temperatures a chaque pas de temps
  158. CONF0 = FORM ;
  159. REPE B1 NB1 ;
  160. TPS1 = TAB1 . 'TEMPS' . (&B1 - 1) ;
  161. LTPS = LTPS ET TPS1 ;
  162. TE1 = TAB1 . 'TEMPERATURES' . (&B1 - 1) ;
  163. U1 = TAB1 . 'DEPLACEMENTS' . (&B1 - 1) ;
  164. CONF1 = FORM U1 ;
  165. J1 = NORM (P3 MOIN P2) ;
  166. T2 = EXTR TE1 'T' P2 ;
  167. T3 = EXTR TE1 'T' P3 ;
  168. MESS TPS1 J1 T2 T3 ;
  169. LJEU = LJEU ET J1 ;
  170. LT2 = LT2 ET T2 ;
  171. LT3 = LT3 ET T3 ;
  172. R1 = TAB1 . 'REACTIONS' . (&B1 - 1) ;
  173. PC = (MAXI 'ABS' (RESU (REDU R1 L3))) / (MESU L3) ;
  174. LPC = LPC ET PC ;
  175. SI ITRAC ;
  176. TRAC TE1 MAILT (PROG T4 PAS ((T1 - T4) / 20.) T1) 'TITR' (CHAI 'Temperature au temps' TPS1) ;
  177. FINSI ;
  178. FORM CONF0 ;
  179. FIN B1 ;
  180. * jeu P2 <--> P3 VS temps
  181. JEUVST = EVOL 'ROUG' 'MANU' 'Temps' LTPS 'Jeu' LJEU ;
  182. * temperatures en P2 et P3 VS temps
  183. T2VST = EVOL 'ROUG' 'MANU' 'Temps' LTPS 'Temperature' LT2 ;
  184. T3VST = EVOL 'BLEU' 'MANU' 'Temps' LTPS 'Temperature' LT3 ;
  185. * pression de contact VS temps
  186. PCVST = EVOL 'ROUG' 'MANU' 'Temps' LTPS 'Pression contact' LPC ;
  187. TLEG = TABL ;
  188. TLEG . 'TITRE' = TABL ;
  189. TLEG . 'TITRE' . 1 = 'Point P2';
  190. TLEG . 'TITRE' . 2 = 'Point P3' ;
  191. SI ITRAC ;
  192. DESS JEUVST 'TITR' 'Jeu VS temps' ;
  193. DESS (T2VST ET T3VST) 'TITR' 'Temperature VS temps' 'LEGE' TLEG ;
  194. DESS PCVST 'TITR' 'Pression de contact VS temps' ;
  195. FINSI ;
  196.  
  197. FIN ;
  198.  
  199.  
  200.  
  201.  

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