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 plaque trouee en traction *
  5. * Comportement elastique lineaire *
  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 (contraintes, courbe force/deplacement) *
  13. * *
  14. ************************************************************************
  15.  
  16.  
  17. ** Options generales
  18. OPTI 'DIME' 2 'ELEM' 'QUA8' 'MODE' 'PLAN' 'DEFO' 'ECHO' 0 ;
  19. ITRAC = FAUX ;
  20.  
  21.  
  22. ** Parametres geometriques
  23. L1 = 1. ;
  24. R1 = 0.15 ;
  25.  
  26.  
  27. ** Maillage
  28. DENS1 = R1 / 15. ;
  29. DENS2 = L1 / 15. ;
  30. P0 = 0. 0. ;
  31. P1 = R1 0. ;
  32. P2 = 0. R1 ;
  33. P3 = L1 0. ;
  34. P4 = L1 L1 ;
  35. P5 = 0. L1 ;
  36. CER = CERC 20 P1 P0 P2 ;
  37. LID = DROI 10 P3 P4 ;
  38. LIH = DROI 10 P4 P5 ;
  39. S1 = REGL CER (LID ET LIH) 'DINI' DENS1 'DFIN' DENS2 ;
  40. LIG = S1 COTE 2 ;
  41. LIB = S1 COTE 4 ;
  42. CS1 = CONT S1 ;
  43.  
  44.  
  45. ** Modele et materiau
  46. MO = MODE S1 'MECANIQUE' 'ELASTIQUE' ;
  47. MA = MATE MO 'YOUN' 200.E9 'NU' 0.3 ;
  48.  
  49.  
  50. ** Blocages mecaniques
  51. BL1 = BLOQ 'UY' LIB ;
  52. BL2 = BLOQ 'UX' LIG ;
  53. BL3 = BLOQ 'UY' LIH ;
  54. BL = BL1 ET BL2 ET BL3 ;
  55.  
  56.  
  57. ** Chargement en deplacement impose
  58. UMAX = L1 / 20. ;
  59. F1 = DEPI BL3 UMAX ;
  60. EV1 = EVOL 'MANU' (PROG 0. 1.) (PROG 0. 1.) ;
  61. CH = CHAR 'DIMP' F1 EV1 ;
  62.  
  63.  
  64. ** Resolution
  65. TAB1 = TABL ;
  66. TAB1 . 'MODELE' = MO ;
  67. TAB1 . 'CARACTERISTIQUES' = MA ;
  68. TAB1 . 'BLOCAGES_MECANIQUES' = BL ;
  69. TAB1 . 'CHARGEMENT' = CH ;
  70. TAB1 . 'TEMPS_CALCULES' = PROG 0.1 'PAS' 0.1 1. ;
  71. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  72. PASAPAS TAB1 ;
  73.  
  74.  
  75. ** Post-traitement (deformees, contraintes, reaction aux appuis, courbe force/deplacement)
  76. TT = TAB1 . 'TEMPS' ;
  77. TU = TAB1 . 'DEPLACEMENTS' ;
  78. TS = TAB1 . 'CONTRAINTES' ;
  79. TR = TAB1 . 'REACTIONS' ;
  80. LISO = PROG 0. 'PAS' 4.E9 4.E10 ;
  81. DEF0 = DEFO CS1 (TU . 0) 0. 'ROUG' ;
  82. LTPS = PROG 0. ;
  83. LDEP = PROG 0. ;
  84. LFOR = PROG 0. ;
  85. NT1 = (DIME TT) - 1 ;
  86. REPE B1 NT1 ;
  87. * instant
  88. TPS1 = TT . &B1 ;
  89. LTPS = LTPS ET TPS1 ;
  90. * contraintes
  91. SG1 = TS . &B1 ;
  92. SP1 = EXCO 'SI11' (PRIN SG1 MO) ;
  93. * deplacements
  94. U1 = TU . &B1 ;
  95. UY1 = EXTR U1 'UY' P5 ;
  96. LDEP = LDEP ET UY1 ;
  97. * reactions
  98. R1 = TR . &B1 ;
  99. VEC1 = VECT R1 5.E-10 'FORC' 'ROSE' ;
  100. RES1 = MAXI 'ABS' (RESU (REDU R1 LIB)) ;
  101. LFOR = LFOR ET RES1 ;
  102. * deformees
  103. DEF1 = DEFO S1 U1 1. VEC1 SP1 MO ;
  104. * traces
  105. TIT1 = CHAI 'Contraintes au temps ' TPS1 ;
  106. SI ITRAC ;
  107. TRAC (DEF1 ET DEF0) 'TITR' TIT1 LISO ;
  108. FINSI ;
  109. FIN B1 ;
  110. * courbe force vs deplacement
  111. EVFC = EVOL 'MANU' 'Deplacement' LDEP 'Force' LFOR ;
  112. SI ITRAC ;
  113. DESS EVFC 'TITR' 'Force vs Deplacement' ;
  114. FINSI ;
  115.  
  116.  
  117. FIN ;
  118.  
  119.  
  120.  

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