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. * ajout de option epsilon lineaire pour la precision des test!
  12. OPTION epsilon lineaire;
  13. *
  14. * test comportement
  15. *
  16. p1=(0 0); p2= 1 0;
  17. l1 = manu seg2 p1 p2;
  18. mo = model l1 MECANIQUE ELASTIQUE PLASTIQUE parfait barr;
  19. ma = mater mo young 1 nu 0.0 rho 0 sigy 1 sect 1;
  20. mas = masse 'DEPL' 1 (manu poi1 p2);
  21. bl = bloq depl (p1 );
  22. lis1=prog 0 pas 0.1 3;
  23. fo=force (manu poi1 p2 ) (2.5 0);
  24. li1 = prog 0 1000;
  25. li2 = prog 1 1;
  26. ev = evol manu t li1 f(t) li2;
  27. cha1 = char meca fo ev;
  28. tab1 = table;
  29. tab1.'BLOCAGES_MECANIQUES'=bl;
  30. tab1.'MODELE'=mo;
  31. tab1.'CARACTERISTIQUES'=ma;
  32. tab1.'CHARGEMENT'=cha1;
  33. tab1.'DYNAMIQUE'=vrai;
  34. tab1.'MASSE_CONSTANTE'=mas;
  35. tab1.'TEMPS_CALCULES'=lis1;
  36. tab1.'SOUS_INCREMENT'=100;
  37. tab1.'MAXDEFOR'=1e-1;
  38. pasapas tab1;
  39. vxx = tab1.vitesses.30 extr p2 ux;
  40. uxx = tab1.deplacements.30 extr p2 ux;
  41. ecin = vxx*vxx / 2.;
  42. eabs = 1.0*(uxx - 1) + 0.5;
  43. eext = 2.5* uxx;
  44. optio echo 0;
  45. mess 'position ' uxx ;
  46. mess 'vitesse ' vxx ;
  47. mess 'energie cinetique ' ecin ;
  48. mess 'energie absorbee ' eabs ;
  49. mess 'travail ' eext ;
  50. bilan = eext - ecin - eabs;
  51. mess 'bilan ' bilan;
  52. optio echo 1;
  53. si ((ABS bilan) > 5E-3) ; erreur 5; finsi;
  54.  
  55.  
  56. *
  57. * test choc
  58. *
  59. ma = mater mo young 1 nu 0.0 rho 0 sigy 1000 sect 1;
  60. blu = bloq 'MAXI' 'UX' p2;
  61. lis2=prog 0 pas 0.1 1;
  62. di = depi blu 0.5;
  63. cha2 = cha1 et (char dimp di ev);
  64. tab2 = table;
  65. tab2.'BLOCAGES_MECANIQUES'=(bl et blu);
  66. tab2.'MODELE'=mo;
  67. tab2.'CARACTERISTIQUES'=ma;
  68. tab2.'CHARGEMENT'=cha2;
  69. tab2.'DYNAMIQUE'=vrai;
  70. tab2.'MASSE_CONSTANTE'=mas;
  71. tab2.'TEMPS_CALCULES'=lis2;
  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.  

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