Télécharger formation_pasapas_2_solution.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * FORMATION AVANCEE SUR LA PROCEDURE PASAPAS *
  3. * *
  4. * Rupture d'une poutre en flexion 3 points *
  5. * Comportement elasto-plastique a ecrouissage isotrope *
  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 *
  11. * - la resolution avec PASAPAS *
  12. * - la suppression d'elements du maillage selon un critere de rupture *
  13. * base sur la 1ere contrainte principale (avec la procedure PERSO1) *
  14. * - le post traitement (courbe effort/fleche) *
  15. * *
  16. ************************************************************************
  17.  
  18.  
  19. ** Options generales
  20. OPTI 'DIME' 2 'ELEM' 'QUA4' 'MODE' 'PLAN' 'DEFO' 'ECHO' 0 ;
  21. ITRAC = FAUX ;
  22. XPAS = 0.05 ;
  23.  
  24.  
  25. ** Parametres geometriques
  26. L = 40. ;
  27. EP = 10. ;
  28. ENT = 1. ;
  29.  
  30.  
  31. ** Maillage
  32. NE1 = 10 ;
  33. DENS1 = EP / NE1 ;
  34. OPTI 'DENS' DENS1 ;
  35. P1 = 0. 0. ;
  36. P2 = 0. EP ;
  37. P3 = ((L / 2.) - (4. * ENT)) 0. ;
  38. P4 = ((L / 2.) - (4. * ENT)) EP ;
  39. L12 = DROI P1 P2 ;
  40. L34 = DROI P3 P4 ;
  41. S1 = REGL L12 L34 ;
  42. DENS2 = ENT / NE1 ;
  43. OPTI 'DENS' DENS2 ;
  44. P5 = ((L / 2.) - ENT) 0. ;
  45. P6 = (L / 2.) ENT ;
  46. P7 = ((L / 2.) + ENT) 0. ;
  47. OPTI 'DENS' DENS1 ;
  48. P8 = ((L / 2.) + ENT) EP ;
  49. P9 = (L / 2.) EP ;
  50. P10 = ((L / 2.) - ENT) EP ;
  51. CE = DROI P5 P6 DROI P7 DROI P8 DROI P9 DROI P10 DROI P5 ;
  52. S2 = SURF CE ;
  53. L410 = DROI P4 P10 ;
  54. L105 = CE ELEM 'COMP' P10 P5 ;
  55. L53 = DROI P5 P3 ;
  56. S3 = SURF (L34 ET L410 ET L105 ET L53) ;
  57. S4 = (S1 ET S3) SYME 'DROI' P6 P9 ;
  58. S1 = S1 ET S2 ET S3 ET S4 ;
  59. ELIM S1 1.E-9 ;
  60. CS1 = CONT S1 ;
  61. P11 = CS1 POIN 'PROC' (L 0.) ;
  62. SI ITRAC ;
  63. TRAC S1 'TITR' 'Maillage' ;
  64. FINSI ;
  65.  
  66.  
  67. ** Modele et materiau
  68. MO = MODE S1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'CINEMATIQUE' ;
  69. MA = MATE MO 'YOUN' 40000.E6 'NU' 0.3 'SIGY' 400.E6 'H' 5000.E6 ;
  70.  
  71.  
  72. ** Blocages mecaniques
  73. CL1 = BLOQ 'UY' (P1 ET P11) ;
  74. CL2 = BLOQ 'UX' P9 ;
  75.  
  76.  
  77. ** Chargement en effort impose
  78. FMAX = 1.E6 ;
  79. F1 = FORC (0. -1.) P9 ;
  80. EV1 = EVOL 'MANU' (PROG 0. 1.) (PROG 0. FMAX) ;
  81. CHA1 = CHAR 'MECA' F1 EV1 ;
  82.  
  83.  
  84. ** Procedure PERSO1 pour supprimer les elements selon un critere sur
  85. * les contraintes principales
  86. DEBP PERSO1 T1*'TABLE' ;
  87. ES = T1 . 'ESTIMATION' ;
  88. WT = T1 . 'WTABLE' ;
  89. MO1 = WT . 'MO_TOT' ;
  90. MAIL0 = EXTR MO1 'MAIL' ;
  91. SG1 = ES . 'CONTRAINTES' ;
  92. SGP1 = PRIN SG1 MO1 ;
  93. SG11 = CHAN 'GRAVITE' (EXCO 'SI11' SGP1) MO1 ;
  94. * maillage sain a conserver
  95. MAIL1 = SG11 ELEM 'INFERIEUR' 1.E6 ;
  96. * maillage complementaire (les elements retires)
  97. MAIL2 = DIFF MAIL0 MAIL1 ;
  98. * on supprime aussi les elements "pendants"
  99. SI ((NBEL MAIL2) > 0) ;
  100. MAIL3 = MAIL1 ELEM 'APPUYE' 'STRICTEMENT' (CONT MAIL2) 'NOVERIF' ;
  101. SI ((NBEL MAIL3) > 0) ;
  102. MAIL1 = DIFF MAIL1 MAIL3 ;
  103. FINSI ;
  104. FINSI ;
  105. * suppression des elements du modele
  106. MO2 = REDU MO1 MAIL1 ;
  107. WT . 'MO_TOT' = MO2 ;
  108. WT . 'MO_TOTAL' = MO2 ;
  109. WT . 'MOD_MEC' = MO2 ;
  110. FINP ;
  111.  
  112.  
  113. ** Resolution
  114. TAB1 = TABL ;
  115. TAB1 . 'MODELE' = MO ;
  116. TAB1 . 'CARACTERISTIQUES' = MA ;
  117. TAB1 . 'BLOCAGES_MECANIQUES' = CL1 ET CL2 ;
  118. TAB1 . 'CHARGEMENT' = CHA1 ;
  119. TAB1 . 'TEMPS_CALCULES' = PROG XPAS 'PAS' XPAS 1. ;
  120. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  121. TAB1 . 'PROCEDURE_PERSO1' = VRAI ;
  122. PASAPAS TAB1 ;
  123.  
  124.  
  125. ** Post-traitement (trace des deformees et des contraintes)
  126. TT = TAB1 . 'TEMPS' ;
  127. TU = TAB1 . 'DEPLACEMENTS' ;
  128. TS = TAB1 . 'CONTRAINTES' ;
  129. DEF0 = DEFO CS1 (TU . 0) 0. 'VERT' ;
  130. LTPS = PROG 0. ;
  131. NT1 = (DIME TT) - 1 ;
  132. REPE B1 NT1 ;
  133. TPS1 = TT . &B1 ;
  134. LTPS = LTPS ET TPS1 ;
  135. SG1 = TS . &B1 ;
  136. MAIL1 = EXTR SG1 'MAIL' ;
  137. CONT1 = CONT MAIL1 ;
  138. MO1 = REDU MO MAIL1 ;
  139. U1 = REDU (TU . &B1) MAIL1 ;
  140. DEF1 = DEFO CONT1 U1 4.E3 'ROUG' ;
  141. TIT1 = CHAI 'Deformee au temps ' TPS1 ;
  142. SI ITRAC ;
  143. * TRAC (DEF0 ET DEF1) 'TITR' TIT1 ;
  144. FINSI ;
  145. DEF1 = DEFO MAIL1 U1 4.E3 ;
  146. SGP1 = PRIN SG1 MO1 ;
  147. TIT1 = CHAI 'Contraintes pinc. au temps ' TPS1 ;
  148. SI ITRAC ;
  149. TRAC SGP1 MO1 DEF1 CONT1 'TITR' TIT1 ;
  150. FINSI ;
  151. FIN B1 ;
  152. * courbe force vs fleche
  153. EVU = EVOL 'TEMP' TAB1 'DEPLACEMENTS' 'UY' P9 ;
  154. LU = ABS (EXTR EVU 'ORDO') ;
  155. LF = IPOL LTPS EV1 ;
  156. EVFU = EVOL 'ROUG' 'MANU' 'Deplacement' LU 'Force' LF ;
  157. SI ITRAC ;
  158. DESS EVFU 'TITR' 'Force vs Fleche' ;
  159. FINSI ;
  160.  
  161. FIN ;
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  

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