Accueil > Cours de Physique > 6ième année > Chapitres 10 et 14 Oscillateurs mécaniques et modélisation des (...) > 6A Oscillateurs mécaniques et méthode d’Euler

6A Oscillateurs mécaniques et méthode d’Euler

mercredi 30 janvier 2013 par Luc Tartier

En mécanique, l’équation différentielle vient de l’application de la seconde loi de Newton :
 
Pour un pendule élastique horizontal, le poids est annulé par la réaction de la barre et si les frottements étaient nuls on aurait : m.a = m.dv/dt = m.d²x/dt² = -k.x (projection de la relation vectorielle sur un axe)
 
 
 
La méthode d’Euler (1707-1783) est facile à comprendre et à appliquer avec un tableur-grapheur. On prendra comme exemple une masse m = 0,200 kg attachée à un ressort de raideur k = 50N.m-1 :
  1. On part des conditions initiales v0 = 0 m/s et x0 = Xmax = 0,05 m (On écarte la masse m de 5cm de sa position d’équilibre et on la lâche sans vitesse initiale.).
  2. On laisse le temps s’écouler pas à pas par intervalle δt = 0,01s très petits devant la période attendue. On crée dans le tableau une première colonne pour le temps qui commence à t0 = 0s et qui est incrémentée de δt jusqu’à t = 0,5s.
  3. L’accélération est donnée par la seconde loi de Newton. Ici ai = -k.xi/m et dépend de l’allongement du ressort à un instant donné. On crée une seconde colonne pour l’accélération qui commence t0 = 0s) à a0 = -k.x0/m 
  4. Par définition l’accélération est la dérivée de la vitesse par rapport au temps : a = dv/dt. Cela devient, à l’instant de date ti, ai = (vi+1 - vi)/δt = -k.xi/m avec la deuxième loi de Newton. On en déduit vi+1 = vi -k.xi/m.δt. On crée une troisième colonne pour la vitesse qui commence à v0 = 0 m/s. Pour calculer les prochaines valeurs de vi, on applique la formule à tout le reste de la colonne.
  5. Par définition, la vitesse est la dérivée de la position par rapport au temps v = dx/dt. On retrouve donc le même raisonnement : à l’instant de date ti, vi = (xi+1 - xi)/δt que l’on transforme en vi.δt = xi+1 - xi et xi+1 = xi + vi.δt. On crée une quatrième colonne pour la position qui commence à x0 = Xmax = 0,05 m et on applique la formule pour calculer les prochaines valeurs de la position.
  6. On trace les graphiques a(t) v(t) et x(t) pour faire apparaître l’évolution de la masse qui devrait être sinusoïdale si l’intervalle de temps δt est suffisamment petit.
 
A) Ouvrir le fichier joint "TP_OscillateurHarmoniqueEuler6A.xls" pour remplir les colonnes préparées à gauche pour l’oscillateur idéal et voir les graphiques s’afficher.
 
B) Refaire la réflexion pour une oscillateur amorti par des frottements fluides dont l’intensité est proportionnelle à la vitesse f = b.v et remplir le tableau de droite.

Documents joints


Contact | | Plan du site | Mentions légales | Suivre la vie du site RSS 2.0 | Haut de page | SPIP |
Physique et Chimie au lycée Jan Neruda de Prague (Lycée Bilingue Jan Neruda)