Bonjour,
je cherche a résoudre un système d'oscillation simple (masse + ressort), mais j'ai un petit doute avec mon algorithme (methode euler-chromer)
Les oscillations obéissent a la loi suivante :
mx''̈ + bx' + kx = f (t) (x'=dx/dt, x''=dx'/dt)
que l'on peut séparer en deux équations en introduisant les constantes
v' + γx' + ω0^2 x = f (t)/m et v(t) = x' (γ=b/m, ω0^2=k/m)
et avec la methode euler chromer on a :
vn+1 = vn − (γ*vn + ω0^2* xn )∆t + (1/m)*fn* ∆t
et xn+1 = xn + vn+1* ∆t
Mon problème et au niveau de la définition de f(t)=A*cos(Ωt),
fn correspondant sera simplement fn=A*cos(Ω*∆t) ?
Merci de votre aide.
̇
Partager