Télécharger la plaquette Cast3M


fr en

Télécharger drx_impact_anneau.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  drx_impact_anneau.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * chute et rebond d'un anneau dans un cone rigide
  5. * calcul drexus explicite avec impact ligne-ligne
  6. *
  7. * hypotheses : - materiau elastique lineaire
  8. *              - grandes deformations
  9. *
  10. *
  11. *
  12. *                *******
  13. *               *       *
  14. *         *     *       *     *
  15. *          *    *       *    *
  16. *           *    *******    *
  17. *            *      |      *
  18. *             *     |     *
  19. *              *    v    *
  20. *               *       *
  21. *                *     *
  22. *                 *   *
  23. *                  * *
  24. *                   *
  25. *
  26. *
  27. graph= faux;
  28. opti dime 2 elem seg2 mode plan cont;
  29.  
  30. * maillage
  31.  
  32. pc = 0. 10. ; p1 = 0 0 ;
  33. l1 = lign rota 60 pc p1 360. ;
  34. elim l1 .1 ;
  35.  
  36. tg20 = (sin 20.) / (cos 20.) ;
  37. e = 10 - (10.1/(sin 20.)) ;
  38. p2 = 0  e ; p3 = ((20 - e) * tg20) 20 ;
  39. p4 = (-1*(20 - e) * tg20) 20 ;
  40. l2 = p3 d 1 p2 d 1 p4;
  41. tout = l1 et l2 ;
  42.  
  43. * donnees drexus
  44.  
  45. mod1 = mode l1 mecanique coq2 ;
  46. mat1 = mate mod1 youn 100 nu 0. rho .01 epai 1. dim3 1. ;
  47. vit0 = manu chpo l1 3 ux 0 uy -4 rz 0.;
  48.  
  49. etab = table;
  50. etab . modele = mod1 ;
  51. etab . grandes_deformations = vrai ;
  52. etab . caracteristiques = mat1;
  53. etab . vitesse_initiale = vit0;
  54. etab . frequence_sortie = 20 ;
  55. etab . pas_temps = 2e-3 ;
  56. etab . npasmax = 2000;
  57. etab . impact = table;
  58. etab . impact . maitre = l2 ;
  59. etab . impact . esclave = l1 ;
  60.  
  61. drexus etab ;
  62.  
  63. * post-traitement
  64.  
  65. i = 1 ;
  66. lt = prog 0.;
  67. lz = prog 0. ;
  68.  
  69. repeter bou1 ( (dime etab . deplacements ) - 1) ;
  70.   lt = lt et ( prog etab . temps . i ) ;
  71.   dep1 = etab . deplacements . i;
  72.   def1 = defo tout dep1 1;
  73.   si (ega i 1);
  74.     def = def1;
  75.   sinon;
  76.     def = def et def1;
  77.   finsi;
  78.   lz1=extr dep1 p1 UY;
  79.   lz = lz et ( prog lz1 );
  80.   i = i + 1 ;
  81. fin bou1 ;
  82.  
  83. evdz = evol bleu  manu 'Temps' lt 'Dy' lz  ;
  84. si  graph ;
  85. dess evdz;
  86.  
  87. trac def anime;
  88. finsi;
  89.  
  90. * test de fonctionnement
  91.  
  92. err1 = 1. + (lz1 / 15.902) ;
  93. mess 'Erreur sur la force ' (err1 * 100. ) '%' ;
  94. lerr1 = err1 >eg 0.05 ;
  95.  
  96. si lerr1 ;
  97.   mess 'Erreur dans le cas test impact_anneau' ;
  98.   erreur 5 ; 
  99. finsi ;
  100.  
  101.  
  102. fin;
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  

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