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

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