Yaya Toure L2Modelisation Projet EDO Simulation Numerique Du Systeme Oscillant

Projet EDO: simulation
numérique de
Étudiant en Licence 2
déplacement d’un Modélisation
Mathématique, Analyse et
système oscillant(masse Simulation Numérique
+ ressort)
Yaya Touré / 76-593-03-47

Responsable du cours: Dr O Diop
Centre National de Calculs Scientifiques | Diamniadio
Description du probleme:
Le déplacement x (t) d’un systeme oscillant composé d’une masse
et d’un ressort, soumis à une force de frottement proportionnelle à
la vitesse, est d’ecrit par l’equation differentielle du second ordre
:mx ” + λx ′ + kx = 0 (∗) avec m : masse du chariot et λ > 0 :
modelise les frotements et k : constante de raideur
1.Mettons mx ” + λx ′ + kx = 0 sous la forme de X ′ = AX
Pour transformer l’equation differentielle du second ordre sous la
forme de X ′ = AX , il suffit juste de poser x ′ = y relation (1) et de
trouver l’expression de y ′ dans (∗). En effet remplaçons
x ′ = y et x ” = y ′ dans (∗) pour obtenir l’expression de y ′ , donc
my ′ + λy + kx = 0 =⇒ my ′ = −kx − λy =⇒ y ′ = − m k
x−m λ
y
relation (2). Nous allons combiner la relation (1) et la relation (2)
pour en fair un systeme d’equation, Et de cette dernier on peut
facilement recuperer la matrice A
Systeme obtenu a l’aide des relations (1) et (2)
Rappels des relations : x ′ = y = 0x + y relation(1) et
y′ = −m k
x−mλ
y relation (2) ; finalement on obtient le systeme
suivant :
 ′
 x = 0x + y
k λ (1)
 y′ = − x− y
m m
De cette systeme (1) on en deduit aisement la relation X ′ = AX
! ! !
x′ 0 1 x
= −k −λ
y′ m m y
X ′ = AX CQFD
2.Determinons en fonction de α et ω02 , les valeurs propres de la
matrice A
Pour trouver les valeurs propres de la matrice A,nous devons
proceder comme suit :
!
1 0
✓ on definit la matrice identite (2X 2) : I = =⇒
0 1
! !
γ 0 0 1
γI = et rappelons la matrice A = −k −λ
0 γ m m
✓ Evaluons: ! ! !
0 1 γ 0 −γ 1
A − γI = −k −λ - = −k λ
m m 0 γ m −( m + γ)
✓ Calculons: det(A − γI) = γ(γ + λ
m) - ( −k
m ) = γ2 + λ
mγ + k
m
✓ Resolvons det(A − γI) = 0 pour trouver les valeurs propres,on
λ k
pose α = 2m ⇐⇒ λ = 2αm et ω02 = m il suffit juste de
2 2 λ k
remplacer λ et ω0 dans γ + m γ + m = 0
⇐⇒ γ 2 + 2αm 2 2
m γ + ω0 = 0 ⇐⇒ γ + 2αγ + ω0 = 0 (∗∗)
2
Calculons ∆ = b 2 − 4ac dans l’equation γ 2 + 2αγ + ω02 = 0 :
On voit bien que notre equation depend bien des coefficients α et
ω02 , donc les valeurs propres de la matrice dependeront bien
d’eux.Enqeffet : ∆ = 4α2 − 4ω02 = 4(α2 − ω02 ) > 0 =⇒
√
∆ = 2 α2 − ω02 , nous avons des valeurs reels:
√
−2α+2 α2 −ω02
q
γ1 = 2 = −α + α2 − ω02 et
√
−2α−2 α2 −ω02
q
γ2 = 2 = −α − α2 − ω02
Finalement nous obtenons les valeurs propres de la matrice A en
fonction de α et ω02
q q
γ1 = −α + α2 − ω02 et γ2 = −α − α2 − ω02 CQFD
On en deduit la soltion general:
X (t) = c1 e γ1 t + c2 e γ2 t
3. Resolvons le probleme lineair par la methode de
Runge-Kutta d’ordre 4 et commenter les resultats pour les
scenarios suivants :
Rappels domaines de travail pour la resolution numerique :
X(0) = 0 , X’(0) = 1 et [0; 5]
(a) Frottements importants, i.e : α2 > ω02 On montrera que

les courbes X (t) solutions de l’equation diff erentielle
presentent alors une asymptote pour t −→ ∞
✓ Si α2 > ω02 ⇐⇒ γ1 > γ2 et cela est verifier quelque soit la
condition, donc prenons γ1 = −1 et γ2 = −3 alors la solution
general devient : X (t) = c1 e −t + c2 e −3t =⇒
V (t) = X ′ (t) = −c1 e −t − 3c2 e −3t
✓ Evaluons X (0) = 0 ⇐⇒ c1 + c2 = 0 =⇒ c1 = −c2
✓ Evaluons V (0) = X ′ (0) ⇐⇒ −c1 − 3c2 = 1 ⇐⇒
c2 − 3c2 = 1 ⇐⇒ −2c2 = 1 =⇒ c2 = − 12 avec c1 = -c2 =
−(− 12 ) = 12 donc c1 = 12 . Finalement nous avons:
X (t) = 12 e −t − 12 e −3t et V (t) = X ′ (t) = − 12 e −t − 32 e −3t
Cas α2 − ω02 > 0 (voirs les codes maxima et scilab ):
Graph Maxima
"On obtient une asymptote horizontal lorsque X(t) tend vers infini"
"Ce qui fait que le résultat de la limitte donne — 0—-" " "et c’est
la raison pour laquelle, les courbes tendent vers une asymptote
horizontal (y=0) "
Graph scilab avec la méthode RK4
Cette constante ce justifie par la presence d’asymptote horizontal

y = 0 quand X (t) −→ ∞
(b) Frotement faibles : 0 < α2 < ω02 :

En effet c’est le cas ou notre ∆ < 0 ce qui fait que nous allons
travailler dans l’ensemble C , alors il suffit de voir dans les deux
cotes et apres de long calculs on prend α = 0, 2 et ω02 = 1, 04
(b) Frotement faibles : 0 < α2 < ω02 : (voir les codes Maxima
et scilab )
que nous allons remplacer dans notre equation precedente:
γ 2 + 2αγ + ω02 = 0 (∗∗) ⇐⇒ γ 2 + 2(0, 2)γ + 1, 04 = 0 (∗∗)
⇐⇒ (∗∗) ⇐⇒ γ 2 + (0, 4)γ + 1, 04 = 0 (∗∗)
✓ Le calcule de ∆ = −4 < 0 =⇒ γ1 = −0, 2 + i et γ2 = −0, 2 − i
=⇒ la solution general est :
X (t) = c1 e −0,2t cos(t) + c2 e −0,2t sin(t)
✓ En appliquant les conditions initials,on trouve:

X (t) = e −0,2t sin(t) et V (t) = X ′ (t) = e −0,2t (cos(t) − 0, 2sin(t))
Cas 0 < α2 < ω02 : (voir les codes Maxima et scilab ):
Graph Maxima | Il
s’agit d’un cas avec des frotement faible,donc la période des
oscillations T = 2π qui est l’intervalle de temps qui sépare 2
passages consécutifs de même nature reste contante pour les deux
courbes. A noté que ces 2 courbes continuent leurs chemin jusqu’à
l’amortissement, c’est le domaine de travail [0; 5](imposé par le
concepteur du projet) qui freine l’activité des 2 courbes
Graph scilab avec RK4
(c) Mouvement sans frottements, i.e :α = 0 ;on montrera que
les solutions X (t) sont des ellipses.
Pour montrer que les solutions sont des ellipses, il suffit juste de
rappeler notre ∆ = 4α2 − 4ω02 = 4(α2 − ω02 ). Etant donné que
nous travaillons dans le cas d’un mouvement sans frottement, i.e
α = 0, donc l’expression de ∆ devient √
∆ = 4(0 − ω02 ) = −4ω02 = i 2 22 ω02 =⇒ ∆ = i2ω0 Naturellement
nous sommes dans C , posons donc le systeme suivant

2 2
 x + y = 2ω0

x2 − y2 = 0 (2)

 2xy = 2ω
0
✓ en évaluant ce système ont trouve :

• si x = ω0 alors y = 1 −→ γ1 = ω0 + i
• si x = −ω0 alors y = −1 −→ γ2 = −ω0 − i
Par suite notre solution elliptique s’écrit :
x 2 (t) y 2 (t)
(ω0 +i)2
+ (−ω0 −i)2
=1
Memes questions dans le cas intermediaire α2 = ω02 .
Montrons en particulier que les solutions X (t) presentent une
asymptote pour t −→ ∞.
En effet : α2 = ω02 ⇐⇒ α2 − ω02 = 0 ce qui revient à dire que
notre ∆ = 0. Rappelons l ′ expression : γ 2 + 2αγ + ω02 = 0 En
prenant α = 2 et ω02 = 4 on a : γ 2 + 2(2)γ + 4 = 0 ⇐⇒
γ 2 + 4γ + 4 = 0 et en pour suivant les calculs on trouve la solution
générale du probléme initial :
x (t) = c1 e −2t + c2 te −2t
✓ en utilisant les conditions initiales on trouve :

X (t) = te −2t et V (t) = X ′ (t) = (1 − 2t)e −2t
Voir dans le graphe l’allure de X (t) et V (t) qui resemble avec
le cas γ 2 > ω02
Maxima Graphe
On voit bien que pour t −→ ∞ les solutions X (t) et X ′ (t)
admetent l’asypmtote horizontal d’ quation y = 0
Cas α2 = ω02 | Graphe scilab avec la méthode de RK4 :
Cette constante ce justifie par le fait que quand x −→ ∞, alors

y =0
Veuillez consulter les codes sources du projet et pour cela vous

devez installer les logiciels Maxima et scilab dans votre ordi
FIN De Projet EDO

Yaya Toure L2Modelisation Projet EDO Simulation Numerique Du Systeme Oscillant

Transféré par

Droits d'auteur :

Formats disponibles

Yaya Toure L2Modelisation Projet EDO Simulation Numerique Du Systeme Oscillant

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Yaya Toure L2Modelisation Projet EDO Simulation Numerique Du Systeme Oscillant

Transféré par

Droits d'auteur :

Formats disponibles

Projet EDO: simulation

Yaya Touré / 76-593-03-47

On en deduit la soltion general:

(a) Frottements importants, i.e : α2 > ω02 On montrera que

Cette constante ce justifie par la presence d’asymptote horizontal

(b) Frotement faibles : 0 < α2 < ω02 :

X (t) = c1 e −0,2t cos(t) + c2 e −0,2t sin(t)

✓ En appliquant les conditions initials,on trouve:

✓ en évaluant ce système ont trouve :

x (t) = c1 e −2t + c2 te −2t

✓ en utilisant les conditions initiales on trouve :

Cette constante ce justifie par le fait que quand x −→ ∞, alors

Veuillez consulter les codes sources du projet et pour cela vous

FIN De Projet EDO

Vous aimerez peut-être aussi