Télécharger formation_pasapas_1_solution.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * FORMATION AVANCEE SUR LA PROCEDURE PASAPAS *
  3. * *
  4. * Flexion d'une poutre en grands deplacements avec chargemet suiveur *
  5. * *
  6. * Ce fichier constitue la mise donnee solution du probleme et realise *
  7. * - le maillage *
  8. * - la definition du modele et des parametres materiaux *
  9. * - la definition du chargement sur la configuration deformee avec la *
  10. * procedure CHARMECA *
  11. * - la resolution avec PASAPAS en hypothese grands deplacements *
  12. * - le post traitement du deplacement de la poutre *
  13. * *
  14. ************************************************************************
  15.  
  16.  
  17. ** Options generales
  18. OPTI 'DIME' 3 'ELEM' 'SEG2' 'ECHO' 0 'COSC' 'BLAN' ;
  19. ITRAC = FAUX ;
  20. OEIL1 = 0. (-1.E6) 0. ;
  21. XPAS = 0.05 ;
  22.  
  23.  
  24. ** Parametres geometriques
  25. LONG = 1. ;
  26. SEC = 1.E-2 ;
  27. INR1 = 1.E-4 ;
  28. TOR1 = 2.E-4 ;
  29.  
  30.  
  31. ** Maillage
  32. P1 = 0. 0. 0. ;
  33. P2 = 0. 0. LONG ;
  34. L1 = DROI 25 P1 P2 ;
  35.  
  36.  
  37. ** Modele et materiau (elastique)
  38. MO1 = MODE L1 'MECANIQUE' 'ELASTIQUE' 'POUT' ;
  39. MA1 = MATE MO1 'YOUN' 50.E6 'NU' 0.3
  40. 'SECT' SEC 'INRY' INR1 'INRZ' INR1 'TORS' TOR1 ;
  41.  
  42.  
  43. ** Blocages mecaniques (encastrement de la base)
  44. BLO1 = BLOQ 'DEPL' 'ROTA' P1 ;
  45.  
  46.  
  47. ** Chargement en effort impose
  48. F1 = FORC (1. 0. 0.) P2 ;
  49. EV1 = EVOL 'MANU' 'Temps' (PROG 0. 1.) 'Coef' (PROG 0. 2.E4) ;
  50. CHAR1 = CHAR 'MECA' F1 EV1 ;
  51.  
  52.  
  53. ** Procedure CHARMECA pour appliquer un effort perpendiculaire a la
  54. * poutre au cours de son deplacement
  55. DEBP CHARMECA T1*'TABLE' TPS1*'FLOTTANT' ;
  56. * donnees utiles
  57. TAM = T1 . 'AMOI' ;
  58. PF = TAM . 'POINT' ;
  59. EV1 = TAM . 'EVOL' ;
  60. U1 = T1 . 'ESTIMATION' . 'DEPLACEMENTS' ;
  61. * angle de rotation de la section droite
  62. RYS = EXTR U1 'RY' PF ;
  63. RYSD = RYS * 180. / PI ;
  64. * direction de l'effort
  65. DIR1 = (COS RYSD) 0. (-1. * (SIN RYSD)) ;
  66. * creation du second membre
  67. VAL1 = IPOL TPS1 EV1 ;
  68. F1 = FORC (VAL1 * DIR1) PF ;
  69. T2 = TABLE ;
  70. T2 . 'ADDI_SECOND' = F1 ;
  71. FINP T2 ;
  72.  
  73.  
  74. ** Resolution
  75. TAB1 = TABL ;
  76. TAB1 . 'MODELE' = MO1 ;
  77. TAB1 . 'CARACTERISTIQUES' = MA1 ;
  78. TAB1 . 'BLOCAGES_MECANIQUES' = BLO1 ;
  79. TAB1 . 'TEMPS_CALCULES' = PROG XPAS 'PAS' XPAS 1. ;
  80. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  81. TAB1 . 'PREDICTEUR' = 'HPP' ;
  82. TAB1 . 'PROCEDURE_CHAR_MECA' = VRAI ;
  83. * donnees supplementaires
  84. TAB1 . 'AMOI' = TABL ;
  85. TAB1 . 'AMOI' . 'EVOL' = EV1 ;
  86. TAB1 . 'AMOI' . 'POINT' = P2 ;
  87. PASAPAS TAB1 ;
  88.  
  89.  
  90. ** Post-traitement
  91. TT = TAB1 . 'TEMPS' ;
  92. TU1 = TAB1 . 'DEPLACEMENTS' ;
  93. TR1 = TAB1 . 'REACTIONS' ;
  94. DEF0 = DEFO L1 (TU1 . 0) 0. ;
  95. CONF0 = FORM ;
  96. * trace des deformees successives
  97. REPE B1 (DIME TT) ;
  98. I1 = &B1 - 1 ;
  99. TPS1 = TT . I1 ;
  100. U1 = TU1 . I1 ;
  101. R1 = TR1 . I1 ;
  102. V1 = VECT R1 5.E-5 'FORC' 'ROUG' ;
  103. DEF1 = DEFO L1 U1 1. V1 'BLEU' ;
  104. SI ITRAC ;
  105. TRAC OEIL1 (DEF0 ET DEF1) 'TITR' (CHAI 'Deformee au temps ' TPS1) ;
  106. FINSI ;
  107. FIN B1 ;
  108. * evolution des deplacements du sommet au cours du temps
  109. TLEG = TABL ;
  110. TLEG . 'TITRE' = TABL ;
  111. TLEG . 'TITRE' . 1 = 'Ux' ;
  112. TLEG . 'TITRE' . 2 = 'Uz' ;
  113. EVUX1 = EVOL 'ROUG' 'TEMP' TAB1 'DEPLACEMENTS' 'UX' P2 ;
  114. EVUZ1 = EVOL 'BLEU' 'TEMP' TAB1 'DEPLACEMENTS' 'UZ' P2 ;
  115. SI ITRAC ;
  116. DESS (EVUX1 ET EVUZ1) 'TITR' 'Deplacements vs Temps' 'LEGE' TLEG ;
  117. FINSI ;
  118.  
  119. FIN ;
  120.  
  121.  
  122.  
  123.  

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