Télécharger la plaquette Cast3M


fr en

Télécharger reacdyna.dgibi

Retour à la liste

Numérotation des lignes : 

  1. *fichier : reacdyna.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * calcul d'un ressort avec une masse au bout 
  6. * par la procedur dynamic
  7. * k=2000 m=100 v0=100 
  8. * 
  9. graph = 'N';
  10. opti dime 2 elem seg2;
  11. pa = 0 0 ; pb = 100 0 ;
  12. liab = pa droi 1 pb;
  13. bly = bloq UY liab;
  14. blx = bloq UX pa;
  15. mod = model liab mecanique elastique barre;
  16. mat = mate mod young 20000 nu 0.3 rho 1.e-15 sect 10.;
  17. ri = rigi mod mat;
  18. ma = masse mod mat;
  19. masad= manu rigidite type masse ( manu poi1 pb)
  20. ( MOTS 'UX') (prog 100);
  21. matot = ma et masad;
  22. x = (coor 1 liab)NOMC UX;list x;
  23. ta = table;
  24. ff = force  ( 0 0 ) pb;
  25. evt = evol manu 'TEMPS' (prog 0. 10000) 'coef' (prog 1. 1.);
  26. ch = char meca ff evt;
  27. ta = table;
  28. ta.'CHAR' = ch;
  29. ta.'DEPL'= manu chpo liab 2 UX 0. UY 0.;
  30. ta.'RIGI' = ri et bly et blx;
  31. ta.'MASS'= matot;
  32. list matot;
  33. dt = 0.005;
  34. ta.'INST' = prog 0. pas dt 0.6;
  35. *ta.'INST' = prog 0. pas dt 0.015;
  36. ta.'SAUV' = vrai;
  37. ta.'VITE' = x;
  38. ta . 'FREQ' = 0.25 / dt;
  39. aa = dynamic ta;
  40. te= prog 0.;
  41. de= prog 0.;
  42. re = prog 0.;
  43. repe tot ( (dime aa) - 1);
  44. tt= aa. ( &tot + 1) . 'TEMP';
  45. te = te et ( prog tt);
  46. bd = aa . ( &tot + 1) . 'DEPL';
  47. md = maxi ( (exco bd 'UX') redu (manu poi1 pb));
  48. de = de et ( prog md);
  49. rr= reac bd blx;
  50. rm= maxi rr;
  51. re = re et ( prog rm);
  52. fin tot;
  53. re = re /1000.;
  54. evdep = evol manu 'temps ' te ' depl ' de;
  55. evrea= evol manu 'temps' te ' reac ' re;
  56. *titr ' deplacement point pb';dess evdep;
  57.  
  58. *titr ' reaction au point pa';dess evrea;
  59. titr 'reaction/1000 et deplacement (k=2000 m=100 V0=100)';
  60. si ( ega graph 'O');
  61. dess ( evdep et evrea);
  62. finsi;
  63. * 0.5mv2 = 0.5 kx2 >>>> 
  64. m=100.; v=100.; k = 2000.;
  65. det = (( 0.5 * m * v * v ) /(0.5 * k))**0.5;
  66. ret = det * k;
  67. mess ' deplacement max théorique' det;
  68. mess ' reaction max théorique ' ret; 
  69. mre= (maxi abs re) * 1000.;
  70. mde= maxi abs de;
  71. ermre = abs (( mre - ret) / ret);
  72. mess ' erreur sur maxi des reactions ' ermre '%';
  73. ermde= abs ((mde - det) / det);
  74. mess ' erreur sur maxi des deplacements ' ermde '%';
  75.  
  76.  
  77. si(( ermre + ermde) > 4.e-5);
  78. erreur 5;
  79. finsi;
  80. fin;
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  

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