Télécharger la plaquette Cast3M


fr en

Télécharger newmark1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  newmark1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Ce cas test verifie que le bilan energetique en dynamique est correct
  6. *
  7. * Verification en comportement et en choc
  8. *
  9. optio dime 2 elem qua4;dens 1;
  10. optio mode plan cont;
  11. *
  12. *   test comportement
  13. *
  14. p1=(0 0); p2= 1 0;
  15. l1 = manu seg2 p1 p2;
  16. mo = model l1  MECANIQUE ELASTIQUE PLASTIQUE parfait barr;
  17. ma = mater mo young 1 nu 0.0 rho 0 sigy 1 sect 1;
  18. mas = masse 'DEPL' 1 (manu poi1 p2);
  19. bl = bloq depl (p1 );
  20. lis1=prog 0 pas 0.1   3;
  21. fo=force (manu poi1 p2 ) (2.5 0);
  22. li1 = prog 0  1000;
  23. li2 = prog 1 1;
  24. ev = evol manu t li1 f(t) li2;
  25. cha1 = char meca fo ev;
  26. tab1 = table;
  27. tab1.'BLOCAGES_MECANIQUES'=bl;
  28. tab1.'MODELE'=mo;
  29. tab1.'CARACTERISTIQUES'=ma;
  30. tab1.'CHARGEMENT'=cha1;
  31. tab1.'DYNAMIQUE'=vrai;
  32. tab1.'MASSE_CONSTANTE'=mas;
  33. tab1.'TEMPS_CALCULES'=lis1;
  34. tab1.'SOUS_INCREMENT'=100;
  35. tab1.'MAXDEFOR'=1e-1;
  36. tab1.'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  37. pasapas tab1;
  38. vxx = tab1.vitesses.30 extr p2 ux;
  39. uxx = tab1.deplacements.30 extr p2 ux;
  40. ecin = vxx*vxx / 2.;
  41. eabs = 1.0*(uxx - 1) + 0.5;
  42. eext = 2.5* uxx;
  43. optio echo 0;
  44. mess 'position            ' uxx   ;
  45. mess 'vitesse             ' vxx   ;
  46. mess 'energie cinetique   ' ecin ;
  47. mess 'energie absorbee    ' eabs ;
  48. mess 'travail             ' eext ;
  49. bilan = eext - ecin - eabs;
  50. mess 'bilan               ' bilan;
  51. optio echo 1;
  52. si ((ABS bilan) > 5E-3) ; erreur 5; finsi;
  53.  
  54.  
  55. *
  56. * test choc
  57. *
  58. ma = mater mo young 1 nu 0.0 rho 0 sigy 1000 sect 1;
  59. blu = bloq 'MAXI' 'UX' p2;
  60. lis2=prog 0 pas 0.1 1;
  61. di = depi blu 0.5;
  62. cha2 = cha1 et (char dimp di ev);
  63. tab2 = table;
  64. tab2.'BLOCAGES_MECANIQUES'=(bl et blu);
  65. tab2.'MODELE'=mo;
  66. tab2.'CARACTERISTIQUES'=ma;
  67. tab2.'CHARGEMENT'=cha2;
  68. tab2.'DYNAMIQUE'=vrai;
  69. tab2.'MASSE_CONSTANTE'=mas;
  70. tab2.'TEMPS_CALCULES'=lis2;
  71. tab2.'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  72. pasapas tab2;
  73. vxx = tab2.vitesses.10 extr p2 'UX';
  74. uxx = tab2.deplacements.10 extr p2 'UX';
  75. ecin = vxx*vxx / 2.;
  76. eabs = uxx*uxx/2.;
  77. eext = 2.5* uxx;
  78. optio echo 0;
  79. mess 'position            ' uxx   ;
  80. mess 'vitesse             ' vxx   ;
  81. mess 'energie cinetique   ' ecin ;
  82. mess 'energie absorbee    ' eabs ;
  83. mess 'travail             ' eext ;
  84. bilan = eext - ecin - eabs;
  85. mess 'bilan               ' bilan;
  86. optio echo 1;
  87. si ((ABS bilan) > 5E-3) ; erreur 5; finsi;
  88. fin;
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  

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