Télécharger la plaquette Cast3M


fr en

Télécharger snap_non_associe.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  snap.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *************************************************************************       
  5. *  TEST SNAP
  6. *                                                                               
  7. *      EXEMPLE D UTILISATION DE LA PROCEDURE PASAPAS ( ET INCREME)               
  8. *                                                                               
  9. *      PROBLEME DE GRANDS DEPLACEMENTS
  10. *                                                                               
  11. *      PROBLEME DU SNAP une seule barre est maillee et calculee
  12. *
  13. *      Utilisation d'une relation non associee entre le deplacement du point de controle
  14. *      et la force de chargement
  15. *
  16. *      Utilisation de K_SIGMA pour eviter la raideur tangente nulle au declanchement du snap
  17. *
  18. *
  19. *         ||
  20. *         ||         
  21. *         \/ F
  22. *         /\
  23. *        /  \
  24. *     | /    \ |        
  25. *     |/      \|
  26. *     |        |                                                       
  27. *************************************************************************       
  28. *                                                                               
  29. *                                                                               
  30. TITRE ' SNAP ' ;
  31. *
  32. *TEMPS ;
  33. *                                                                               
  34. *------------- geometrie   ligne ST formee d' 1  SEG2 -------------------
  35. *              PAR SYMETRIE, ON N'ETUDIE QUE LA MOITIE
  36. * 
  37. OPTI ECHO 0 ;                                                                      
  38. OPTI DIME 2 ELEM SEG2 MODE PLAN CONT ;DENS 1. ;
  39. opti trac psc;
  40. P1 = 0.    1.    ; P2 = 10. 0.  ;
  41. ST = P1 D 1  P2 ;
  42. *                                                                               
  43. *------------ calcul mecanique .  ------------------
  44. *                                                                               
  45.  MO = MODE ST  MECANIQUE  BARR     ;
  46. *
  47. *   MATERIAU ET CARACTERISTIQUES
  48. *
  49. *                                                                               
  50. MA1 = MATE MO YOUN 2.1E11  NU 0.3 ;
  51. CAR1 = CARA MO SECT 0.05 ;
  52. MACA= MA1 ET CAR1;
  53. *                                                                               
  54. *----------- calcul de la rigidite --------------------------------------       
  55. *                                                                               
  56. RI1 =    RIGI MACA  MO ;
  57. *                                                                               
  58. *----------- definition des conditions aux limites ----------------------       
  59. *
  60. CL1 = BLOQ UX P1 ;
  61. CL2 = BLOQ UX P2 ;
  62. CL3 = BLOQ UY P2 ;
  63. CL4 = BLOQ UY p1;
  64. CL = CL1 ET CL2 ET CL3 ;  
  65. *                                                                               
  66. *----------- definition du chargement -----------------------------------       
  67. *                                                                         
  68. FP11 = FORCE (  0 -12.5e5 ) p1;
  69. LIX1 = PROG 0. 40. ;                                                            
  70. LIY1 = PROG 0. 40. ;                                                            
  71. EV1 = EVOL MANU T LIX1 F(T) LIY1 ;                                                                                  
  72. *                                                                               
  73. *----------- resolution par la procedure NONLIN -------------------------       
  74. *            CALCUL EN GRANDS DEPLACEMENTS
  75. *                                                                               
  76. TAB2 = TABLE ; 
  77. **CHA1 = CHAR MECA FP11   EV1 ;                                                                 
  78. TAB2.'GRANDS_DEPLACEMENTS'=VRAI;
  79. *
  80.  
  81. * relation non associee pour avoir un chargement en deplacement impose
  82. chs = manu chpo p1 'UY' -1;
  83. cls = rela chs dual (fp11 * 1d-2);
  84. dis = depi cls 1.10 ;
  85.  
  86. TAB2.'BLOCAGES_MECANIQUES' = CL et cls;
  87. TAB2.'MODELE' = MO;
  88. TAB2.'CARACTERISTIQUES' = MACA;
  89. TAB2.'CHARGEMENT' =  char dimp dis ev1;
  90. LIS11 =   PROG 0. pas 0.067  2. ;
  91. TAB2.'TEMPS_CALCULES' = LIS11;
  92. TAB2.'PRECISION' = 1d-8; 
  93. TAB2.'K_SIGMA' = vrai;
  94. PASAPAS TAB2 ;  
  95. *                                                                              
  96. *----------- resultats --------------------------------------------------       
  97. *        courbe de snap through : montee descente montee                                                                   
  98. PGX = PROG 0.;
  99. PGY = PROG 0.;
  100. NDIM = (DIME ( TAB2 . DEPLACEMENTS )) - 1 ;
  101.  
  102. REPETER TBOU2 NDIM ;
  103.  
  104. LEDEP = TAB2 . DEPLACEMENTS.  (&TBOU2);
  105. REA1 = REAC CL3 LEDEP ;
  106. V   = EXTR LEDEP UY P1 ;
  107. PGX = PGX ET ( PROG V ) ;
  108. VV  = EXTR REA1  FY P2 ;
  109. PGY = PGY ET ( PROG VV ) ;
  110. FIN TBOU2 ;
  111. *
  112. pgx = pgx * -1.;
  113. strou = mini pgy;
  114. sref = -2.e+6;
  115. smp= sref * 1.05;smb = sref * 0.95;
  116. mess ' borne min ' smp 'valeur trouvée' strou 'borne max' smb; 
  117. *
  118. si (( strou > smb ) ou ( strou < smp) ) ;
  119.    erreur (5);
  120. finsi;
  121. *
  122. EV5 = EVOL TURQ MANU  'Deplacement' PGX 'Force de reaction' PGY ;
  123. DESS EV5 ;
  124. fin;
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  

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