Télécharger la plaquette Cast3M


fr en

Télécharger formation_pasapas_2_initial.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : formation_pasapas_2_initial.dgibi
  2. *
  3. *----------------------------------------------------------------------*
  4. *              FORMATION AVANCEE SUR LA PROCEDURE PASAPAS              *
  5. *                                                                      *
  6. * Rupture d'une plaque trouee en traction                              *
  7. * Comportement elastique lineaire                                      *
  8. *                                                                      *
  9. * Ce fichier constitue la mise donnee initiale du probleme et realise  *
  10. * - le maillage                                                        *
  11. * - la definition du modele et des parametres materiaux                *
  12. * - la definition du chargement                                        *
  13. * - la resolution avec PASAPAS                                         *
  14. * - le post traitement (contraintes, courbe force/deplacement)         *
  15. *----------------------------------------------------------------------*
  16.  
  17.  
  18.  
  19.  
  20.  
  21. ** Options de trace / Plotting options
  22. itrac = FAUX ;
  23. ***OPTI 'TRAC' 'PSC' ;
  24. ipsc  = EGA (VALE 'TRAC') 'PSC' ;
  25.  
  26. ** Options generales / General options
  27. OPTI 'DIME' 2 'ELEM' 'QUA8' 'MODE' 'PLAN' 'DEFO' ;
  28.  
  29. ** Parametres geometriques / Geometrical parameters
  30. l1    = 1. ;
  31. r1    = 0.15 ;
  32.  
  33. ** Maillage / Mesh
  34. dens1 = R1 / 15. ;
  35. dens2 = l1 / 15. ;
  36. p0    = 0. 0. ;
  37. p1    = r1 0. ;
  38. p2    = 0. r1 ;
  39. p3    = l1 0. ;
  40. p4    = l1 l1 ;
  41. p5    = 0. l1 ;
  42. CER   = CERC 20 p1 p0 p2 ;
  43. lid   = DROI 10 p3 p4 ;
  44. lih   = DROI 10 p4 p5 ;
  45. s1    = REGL cer (lid ET lih) 'DINI' dens1 'DFIN' dens2 ;
  46. lig   = s1 COTE 2 ;
  47. lib   = s1 COTE 4 ;
  48. cs1   = CONT s1 ;
  49.  
  50. ** Modele et materiau / Model and material properties
  51. mo    = MODE s1 'MECANIQUE' 'ELASTIQUE' ;
  52. ma    = MATE mo 'YOUN' 200.E9 'NU' 0.3 ;
  53.  
  54. ** Blocages mecaniques
  55. bl1   = BLOQ 'UY' lib ;
  56. bl2   = BLOQ 'UX' lig ;
  57. bl3   = BLOQ 'UY' lih ;
  58. bl    = bl1 ET bl2 ET bl3 ;
  59.  
  60. ** Chargement en deplacement impose / Load with imposed displacement
  61. umax  = l1 / 20. ;
  62. f1    = DEPI bl3 umax ;
  63. ev1   = EVOL 'MANU' (PROG 0. 1.) (PROG 0. 1.) ;
  64. cha   = CHAR 'DIMP' f1 ev1 ;
  65. SI ipsc ;
  66.   OPTI 'FTRA' 'exo_2_chargement.ps' ;
  67. FINSI ;
  68. SI itrac ;
  69.   TRAC (VECT f1 ' ' 'FLX' 'VERT') ((RIGI mo ma) ET bl1 ET bl2) ;
  70. FINSI ;
  71.  
  72. ** Resolution
  73. tab1                         = TABL ;
  74. tab1 . 'MODELE'              = mo ;
  75. tab1 . 'CARACTERISTIQUES'    = ma ;
  76. tab1 . 'BLOCAGES_MECANIQUES' = bl ;
  77. tab1 . 'CHARGEMENT'          = cha ;
  78. tab1 . 'TEMPS_CALCULES'      = PROG 0.1 'PAS' 0.1 1. ;
  79. PASAPAS tab1 ;
  80.  
  81. ** Post-traitement
  82. liso  = PROG 0. 'PAS' 4.E9 4.E10 ;
  83. box   = (0. -0.2) ET (1. 1.4) ;
  84. nt1   = DIME (tab1 . 'TEMPS') ;
  85. SI ipsc ;
  86.   OPTI 'FTRA' 'exo_2_sigma.ps' ;
  87. FINSI ;
  88. REPE b1 nt1 ;
  89. * contraintes / stresses
  90.   sg1   = tab1 . 'CONTRAINTES' . (&b1 - 1) ;
  91.   sp1   = EXCO 'SI11' (PRIN sg1 mo) ;
  92. * deplacements / displacements
  93.   u1    = tab1 . 'DEPLACEMENTS' . (&b1 - 1) ;
  94. * reactions
  95.   r1    = tab1 . 'REACTIONS' . (&b1 - 1) ;
  96.   vec1  = VECT r1 5.E-10 'FORC' 'VIOL' ;
  97. * deformees / deformed shape
  98.   def1  = DEFO s1 u1 2. vec1 sp1 mo ;
  99. * traces / plot
  100.   SI itrac ;
  101.     TRAC def1 cs1 liso 'BOIT' box ;
  102.   FINSI ;
  103. FIN b1 ;
  104. * courbe force vs deplacement / force vs displacement curve
  105. evu   = EVOL 'TEMP' tab1 'DEPLACEMENTS' 'UY' p5 ;
  106. evf   = -1. * (@GLOBAL tab1 bl1 evu 'FY') ;
  107. evfu  = EVOL 'ROUG' 'MANU' 'Deplac. UY'  (EXTR evu 'ORDO')
  108.                            'Reaction FY' (EXTR evf 'ORDO') ;
  109. SI ipsc ;
  110.   OPTI 'FTRA' 'exo_2_evol.ps' ;
  111. FINSI ;
  112. SI itrac ;
  113.   DESS evfu ;
  114. FINSI ;
  115.  
  116.  
  117. FIN ;
  118.  
  119.  
  120.  

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