Télécharger formation_pasapas_2_initial.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 initiale 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. * - le post traitement (courbe effort/fleche) *
  13. * *
  14. ************************************************************************
  15.  
  16.  
  17. ** Options generales
  18. OPTI 'DIME' 2 'ELEM' 'QUA4' 'MODE' 'PLAN' 'DEFO' 'ECHO' 0 ;
  19. ITRAC = FAUX ;
  20. XPAS = 0.05 ;
  21.  
  22.  
  23. ** Parametres geometriques
  24. L = 40. ;
  25. EP = 10. ;
  26. ENT = 1. ;
  27.  
  28.  
  29. ** Maillage
  30. NE1 = 10 ;
  31. DENS1 = EP / NE1 ;
  32. OPTI 'DENS' DENS1 ;
  33. P1 = 0. 0. ;
  34. P2 = 0. EP ;
  35. P3 = ((L / 2.) - (4. * ENT)) 0. ;
  36. P4 = ((L / 2.) - (4. * ENT)) EP ;
  37. L12 = DROI P1 P2 ;
  38. L34 = DROI P3 P4 ;
  39. S1 = REGL L12 L34 ;
  40. DENS2 = ENT / NE1 ;
  41. OPTI 'DENS' DENS2 ;
  42. P5 = ((L / 2.) - ENT) 0. ;
  43. P6 = (L / 2.) ENT ;
  44. P7 = ((L / 2.) + ENT) 0. ;
  45. OPTI 'DENS' DENS1 ;
  46. P8 = ((L / 2.) + ENT) EP ;
  47. P9 = (L / 2.) EP ;
  48. P10 = ((L / 2.) - ENT) EP ;
  49. CE = DROI P5 P6 DROI P7 DROI P8 DROI P9 DROI P10 DROI P5 ;
  50. S2 = SURF CE ;
  51. L410 = DROI P4 P10 ;
  52. L105 = CE ELEM 'COMP' P10 P5 ;
  53. L53 = DROI P5 P3 ;
  54. S3 = SURF (L34 ET L410 ET L105 ET L53) ;
  55. S4 = (S1 ET S3) SYME 'DROI' P6 P9 ;
  56. S1 = S1 ET S2 ET S3 ET S4 ;
  57. ELIM S1 1.E-9 ;
  58. CS1 = CONT S1 ;
  59. P11 = CS1 POIN 'PROC' (L 0.) ;
  60. SI ITRAC ;
  61. TRAC S1 'TITR' 'Maillage' ;
  62. FINSI ;
  63.  
  64.  
  65. ** Modele et materiau
  66. MO = MODE S1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'CINEMATIQUE' ;
  67. MA = MATE MO 'YOUN' 40000.E6 'NU' 0.3 'SIGY' 400.E6 'H' 5000.E6 ;
  68.  
  69.  
  70. ** Blocages mecaniques
  71. CL1 = BLOQ 'UY' (P1 ET P11) ;
  72. CL2 = BLOQ 'UX' P9 ;
  73.  
  74.  
  75. ** Chargement en effort impose
  76. FMAX = 1.E6 ;
  77. F1 = FORC (0. -1.) P9 ;
  78. EV1 = EVOL 'MANU' (PROG 0. 1.) (PROG 0. FMAX) ;
  79. CHA1 = CHAR 'MECA' F1 EV1 ;
  80.  
  81.  
  82. ** Resolution
  83. TAB1 = TABL ;
  84. TAB1 . 'MODELE' = MO ;
  85. TAB1 . 'CARACTERISTIQUES' = MA ;
  86. TAB1 . 'BLOCAGES_MECANIQUES' = CL1 ET CL2 ;
  87. TAB1 . 'CHARGEMENT' = CHA1 ;
  88. TAB1 . 'TEMPS_CALCULES' = PROG XPAS 'PAS' XPAS 1. ;
  89. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  90. PASAPAS TAB1 ;
  91.  
  92.  
  93. ** Post-traitement (trace des deformees et des contraintes)
  94. TT = TAB1 . 'TEMPS' ;
  95. TU = TAB1 . 'DEPLACEMENTS' ;
  96. TS = TAB1 . 'CONTRAINTES' ;
  97. DEF0 = DEFO CS1 (TU . 0) 0. 'VERT' ;
  98. LTPS = PROG 0. ;
  99. NT1 = (DIME TT) - 1 ;
  100. REPE B1 NT1 ;
  101. TPS1 = TT . &B1 ;
  102. LTPS = LTPS ET TPS1 ;
  103. SG1 = TS . &B1 ;
  104. MAIL1 = EXTR SG1 'MAIL' ;
  105. CONT1 = CONT MAIL1 ;
  106. MO1 = REDU MO MAIL1 ;
  107. U1 = REDU (TU . &B1) MAIL1 ;
  108. DEF1 = DEFO CONT1 U1 4.E3 'ROUG' ;
  109. TIT1 = CHAI 'Deformee au temps ' TPS1 ;
  110. SI ITRAC ;
  111. * TRAC (DEF0 ET DEF1) 'TITR' TIT1 ;
  112. FINSI ;
  113. DEF1 = DEFO MAIL1 U1 4.E3 ;
  114. SGP1 = PRIN SG1 MO1 ;
  115. TIT1 = CHAI 'Contraintes pinc. au temps ' TPS1 ;
  116. SI ITRAC ;
  117. TRAC SGP1 MO1 DEF1 CONT1 'TITR' TIT1 ;
  118. FINSI ;
  119. FIN B1 ;
  120. * courbe force vs fleche
  121. EVU = EVOL 'TEMP' TAB1 'DEPLACEMENTS' 'UY' P9 ;
  122. LU = ABS (EXTR EVU 'ORDO') ;
  123. LF = IPOL LTPS EV1 ;
  124. EVFU = EVOL 'ROUG' 'MANU' 'Deplacement' LU 'Force' LF ;
  125. SI ITRAC ;
  126. DESS EVFU 'TITR' 'Force vs Fleche' ;
  127. FINSI ;
  128.  
  129. FIN ;
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  

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