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 ;
  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. * angle de rotation de la section droite
  61. U1 = T1 . 'ESTIMATION' . 'DEPLACEMENTS' ;
  62. RYS = EXTR U1 'RY' PF ;
  63. **** une autre possibilite : mettre a jour la force sur la config. fin de pas
  64. *** DRYS = 0. ;
  65. *** SI (EGA (TYPE ZDEPT) 'CHPOINT') ;
  66. *** DRYS = EXTR ZDEPT 'RY' PF ;
  67. *** FINSI ;
  68. *** RYS = (EXTR (T1 . 'CONTINUATION' . 'DEPLACEMENTS') 'RY' PF) + DRYS ;
  69. RYSD = RYS * 180. / PI ;
  70. * direction de l'effort
  71. DIR1 = (COS RYSD) 0. (-1. * (SIN RYSD)) ;
  72. * creation du second membre
  73. VAL1 = IPOL TPS1 EV1 ;
  74. F1 = FORC (VAL1 * DIR1) PF ;
  75. T2 = TABLE ;
  76. T2 . 'ADDI_SECOND' = F1 ;
  77. FINP T2 ;
  78.  
  79.  
  80. ** Resolution
  81. TAB1 = TABL ;
  82. TAB1 . 'MODELE' = MO1 ;
  83. TAB1 . 'CARACTERISTIQUES' = MA1 ;
  84. TAB1 . 'BLOCAGES_MECANIQUES' = BLO1 ;
  85. *TAB1 . 'CHARGEMENT' = CHAR1 ;
  86. TAB1 . 'TEMPS_CALCULES' = PROG XPAS 'PAS' XPAS 1. ;
  87. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  88. TAB1 . 'PREDICTEUR' = 'HPP' ;
  89. TAB1 . 'PROCEDURE_CHARMECA' = VRAI ;
  90. TAB1 . 'AMOI' = TABL ;
  91. TAB1 . 'AMOI' . 'EVOL' = EV1 ;
  92. TAB1 . 'AMOI' . 'POINT' = P2 ;
  93. PASAPAS TAB1 ;
  94.  
  95.  
  96. ** Post-traitement
  97. TT = TAB1 . 'TEMPS' ;
  98. TU1 = TAB1 . 'DEPLACEMENTS' ;
  99. TR1 = TAB1 . 'REACTIONS' ;
  100. DEF0 = DEFO L1 (TU1 . 0) 0. ;
  101. CONF0 = FORM ;
  102. * trace des deformees successives
  103. REPE B1 (DIME TT) ;
  104. I1 = &B1 - 1 ;
  105. TPS1 = TT . I1 ;
  106. U1 = TU1 . I1 ;
  107. R1 = TR1 . I1 ;
  108. V1 = VECT R1 5.E-5 'FORC' 'ROUG' ;
  109. DEF1 = DEFO L1 U1 1. V1 'BLEU' ;
  110. SI ITRAC ;
  111. TRAC OEIL1 (DEF0 ET DEF1) 'TITR' (CHAI 'Deformee au temps ' TPS1) ;
  112. FINSI ;
  113. FIN B1 ;
  114. * evolution des deplacements du sommet au cours du temps
  115. TLEG = TABL ;
  116. TLEG . 'TITRE' = TABL ;
  117. TLEG . 'TITRE' . 1 = 'Ux' ;
  118. TLEG . 'TITRE' . 2 = 'Uz' ;
  119. EVUX1 = EVOL 'ROUG' 'TEMP' TAB1 'DEPLACEMENTS' 'UX' P2 ;
  120. EVUZ1 = EVOL 'BLEU' 'TEMP' TAB1 'DEPLACEMENTS' 'UZ' P2 ;
  121. SI ITRAC ;
  122. DESS (EVUX1 ET EVUZ1) 'TITR' 'Deplacements vs Temps' 'LEGE' TLEG ;
  123. FINSI ;
  124.  
  125. FIN ;
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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