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 (ou bien PERSO1) *
  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 'SECT' SEC 'INRY' INR1 'INRZ' INR1 'TORS' TOR1 ;
  40.  
  41.  
  42. ** Blocages mecaniques (encastrement de la base)
  43. BLO1 = BLOQ 'DEPL' 'ROTA' P1 ;
  44.  
  45.  
  46. ** Chargement en effort impose
  47. F1 = FORC (1. 0. 0.) P2 ;
  48. EV1 = EVOL 'MANU' 'Temps' (PROG 0. 1.) 'Coef' (PROG 0. 2.E4) ;
  49. CHAR1 = CHAR 'MECA' F1 EV1 ;
  50.  
  51.  
  52. ** Procedure CHARMECA pour appliquer un effort perpendiculaire a la
  53. * poutre au cours de son deplacement
  54. DEBP CHARMECA T1*'TABLE' TPS1*'FLOTTANT' ;
  55. * donnees utiles
  56. TAM = T1 . 'AMOI' ;
  57. PF = TAM . 'POINT' ;
  58. EV1 = TAM . 'EVOL' ;
  59. * angle de rotation de la section droite
  60. U1 = T1 . 'ESTIMATION' . 'DEPLACEMENTS' ;
  61. RYS = EXTR U1 'RY' PF ;
  62. **** une autre possibilite : mettre a jour la force sur la config. fin de pas
  63. *** DRYS = 0. ;
  64. *** SI (EGA (TYPE ZDEPT) 'CHPOINT') ;
  65. *** DRYS = EXTR ZDEPT 'RY' PF ;
  66. *** FINSI ;
  67. *** RYS = (EXTR (T1 . 'CONTINUATION' . 'DEPLACEMENTS') 'RY' PF) + DRYS ;
  68. RYSD = RYS * 180. / PI ;
  69. * direction de l'effort
  70. DIR1 = (COS RYSD) 0. (-1. * (SIN RYSD)) ;
  71. * creation du second membre
  72. VAL1 = IPOL TPS1 EV1 ;
  73. F1 = FORC (VAL1 * DIR1) PF ;
  74. T2 = TABLE ;
  75. T2 . 'ADDI_SECOND' = F1 ;
  76. FINP T2 ;
  77.  
  78.  
  79. ** Il est aussi possible de modeliser la force suiveuse avec PERSO1
  80. * (mais les modifications des CL sont pises en compte au pas de temps suivant,
  81. * les resultats auront donc un pas de temps de retard !)
  82. DEBP PERSO1 T1*'TABLE' ;
  83. * donnees utiles
  84. TAM = T1 . 'AMOI' ;
  85. PF = TAM . 'POINT' ;
  86. EV1 = TAM . 'EVOL' ;
  87. * angle de rotation de la section droite
  88. U1 = T1 . 'ESTIMATION' . 'DEPLACEMENTS' ;
  89. RYS = EXTR U1 'RY' PF ;
  90. RYSD = RYS * 180. / PI ;
  91. * direction de l'effort
  92. DIR1 = (COS RYSD) 0. (-1. * (SIN RYSD)) ;
  93. * creation du second membre (pour le pas de temps suivant)
  94. TPS1 = (T1 . 'ESTIMATION' . 'TEMPS') + XPAS ;
  95. VAL1 = IPOL TPS1 EV1 ;
  96. F1 = FORC (VAL1 * DIR1) PF ;
  97. CHAR2 = CHAR 'MECA' F1 (EVOL 'MANU' (PROG 0. 1.) (PROG 1. 1.)) ;
  98. T1 . 'WTABLE' . 'CHARGEMENT' = CHAR2 ;
  99. FINP ;
  100.  
  101.  
  102. ** Resolution
  103. TAB1 = TABL ;
  104. TAB1 . 'MODELE' = MO1 ;
  105. TAB1 . 'CARACTERISTIQUES' = MA1 ;
  106. TAB1 . 'BLOCAGES_MECANIQUES' = BLO1 ;
  107. *TAB1 . 'CHARGEMENT' = CHAR1 ;
  108. TAB1 . 'TEMPS_CALCULES' = PROG XPAS 'PAS' XPAS 1. ;
  109. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  110. TAB1 . 'PREDICTEUR' = 'HPP' ;
  111. TAB1 . 'PROCEDURE_CHARMECA' = VRAI ;
  112. ***TAB1 . 'PROCEDURE_PERSO1' = VRAI ;
  113. TAB1 . 'AMOI' = TABL ;
  114. TAB1 . 'AMOI' . 'EVOL' = EV1 ;
  115. TAB1 . 'AMOI' . 'POINT' = P2 ;
  116. PASAPAS TAB1 ;
  117.  
  118.  
  119. ** Post-traitement
  120. TT = TAB1 . 'TEMPS' ;
  121. TU1 = TAB1 . 'DEPLACEMENTS' ;
  122. TR1 = TAB1 . 'REACTIONS' ;
  123. DEF0 = DEFO L1 (TU1 . 0) 0. ;
  124. CONF0 = FORM ;
  125. * trace des deformees successives
  126. REPE B1 (DIME TT) ;
  127. I1 = &B1 - 1 ;
  128. TPS1 = TT . I1 ;
  129. U1 = TU1 . I1 ;
  130. R1 = TR1 . I1 ;
  131. V1 = VECT R1 5.E-5 'FORC' 'ROUG' ;
  132. DEF1 = DEFO L1 U1 1. V1 'BLEU' ;
  133. SI ITRAC ;
  134. TRAC OEIL1 (DEF0 ET DEF1) 'TITR' (CHAI 'Deformee au temps ' TPS1) ;
  135. FINSI ;
  136. FIN B1 ;
  137. * evolution des deplacements du sommet au cours du temps
  138. TLEG = TABL ;
  139. TLEG . 'TITRE' = TABL ;
  140. TLEG . 'TITRE' . 1 = 'Ux' ;
  141. TLEG . 'TITRE' . 2 = 'Uz' ;
  142. EVUX1 = EVOL 'ROUG' 'TEMP' TAB1 'DEPLACEMENTS' 'UX' P2 ;
  143. EVUZ1 = EVOL 'BLEU' 'TEMP' TAB1 'DEPLACEMENTS' 'UZ' P2 ;
  144. SI ITRAC ;
  145. DESS (EVUX1 ET EVUZ1) 'TITR' 'Deplacements vs Temps' 'LEGE' TLEG ;
  146. FINSI ;
  147.  
  148. FIN ;
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  

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