TPN 4 - Introduction À La Résolution Approchée D'Edos: 1 Modélisation: Le Pendule Simple

Faculté des sciences et ingénierie (Toulouse III) Année universitaire
Département de mathématiques – M1 MAF 2012-2013

Introduction à la modélisation
TP no 4 – Introduction à la résolution approchée d’EDOs
Créer dans votre répertoire M1_Modélisation un répertoire TP04, dans lequel seront
placés les fichiers relatifs à ce TP.
1 Modélisation : le pendule simple

On s’intéresse dans ce texte au mouvement d’un pendule simple. On considère ainsi
une bille de masse m assimilée à un point matériel, reliée à un fil rigide, inextensible
de longueur ` et de masse négligeable. Dans un premier temps, on néglige les forces de
frottements subies par la masse au cours de son déplacement.
On note θ l’angle entre la verticale et la bille et on se place dans le repère mobile (e~r , e~θ )
lié à la bille.
Le principe fondamental de la dynamique appliqué à la bille s’écrit :

m~γ = P~ + T~ (1)
1
où P~ et T~ désignent respectivement le poids de la bille et la tension exercée par le fil sur
la bille. Dans le repère mobile lié à la bille, l’accélération ~γ a pour expression
~γ = `θ̈e~θ − `θ̇2 e~r .
Les forces extérieures appliquées à la bille ont pour composante tangentielle respective
P~ · e~θ = −mg sin θ et T~ · e~θ = 0.
En projetant l’équation (1) selon l’axe tangentiel e~θ , on obtient l’équation différentielle
du second ordre
g
θ̈ = − sin θ. (2)
`
Pour résoudre cette équation et obtenir l’évolution de θ en fonction du temps, on se
donne θ(0) = θ0 et θ̇(0) = θ1 les valeurs initiales respectives de l’angle et de la vitesse
angulaire.
2 Pistes pour l’étude théorique du pendule simple

1. Résoudre explicitement l’équation (2) (avec les conditions initiales données) dans
l’approximation “petits angles”, c’est-à-dire en supposant que sin θ peut être ap-
proché par θ.
2. Montrer que l’équation différentielle ordinaire du second ordre (2) peut être mise
sous la forme d’un système autonome d’équations différentielles ordinaires du pre-
mier ordre : !
θ(t)
Y 0 (t) = F (Y (t)), t ≥ 0, où Y (t) = (3)
θ0 (t)
et F : R2 → R2 est à déterminer.
En déduire un résultat d’exitence (d’unicité ?) des solutions de (3).
3. On définit l’énergie cinétique et l’énergie potentielle du système par
1
Ec (t) = m`2 θ̇(t)2 , Ep (t) = mg`(1 − cos θ(t)),
2
et l’énergie mécanique totale du système par Em (t) = Ec (t) + Ep (t).
Montrer que l’énergie mécanique de la bille est conservée au cours du temps.
4. On ajoute pour raffiner un peu le modèle un amortissement qui modélise une force
de frottement f~frott sous la forme
f~frott = −Λ~v , (4)
~v représentant la vitesse de la bille exprimée dans le repère mobile par ~v = `θ̇e~θ .

Le principe fondamental de la dynamique en projection sur l’axe tangentiel au
mouvement e~θ devient alors :
g
θ̈(t) + λθ̇(t) + sin θ(t) = 0. (5)
`
2
– Peut-on résoudre explicitement l’équation dans le cas de l’approximation “petits
angles” ?
– En mettant l’équation (5) sous forme d’un système d’équations différentielles
d’ordre 1, démontrer un résultat d’existence (d’unicité ?) des solutions à (5).
– L’énergie mécanique est-elle toujours conservée au cours du temps ?
3 Illustration avec Matlab de l’étude théorique du pendule

simple
Comme on ne connaît pas de solution explicite des équations (3) et (5), la seule
illustration que l’on peut proposer pour cette partie théorique est la dynamique du
pendule dans l’approximation “petits angles”.
Partant donc d’une position initiale dans cette approximation, on propose d’illustrer
la dynamique du pendule. Pour cela, il y a plusieurs choses que l’on peut représenter :
– l’évolution de l’angle θ en fonction du temps en partant d’un petit angle, avec ou
sans vitesse angulaire initiale. Il s’agira donc, du point de vue de la programmation,
du tracé de la fonction t 7→ θ(t), ce qui correspond à un graphe d’une fonction
connue explicitement.
Exercice 1. Programmer une fonction petitsangles dont les entrées-sorties sont :
entrées : θ0 , θ1
sortie : graphe t 7→ θ(t)
– la dynamique du pendule lui-même en représentant la trajectoire de la bille au
cours du temps. On représentera donc la bille par un point matériel (x, y) dont on
va illustrer la dynamique au cours du temps, partant du fait que, dans le repère
immobile (Oxy), on a :
∀t ≥ 0, x(t) = ` sin θ(t), y(t) = −` cos θ(t).
Du point de vue de la programmation, on va représenter la trajectoire d’un point

2D en fonction du temps. Il s’agira donc d’un “film”, où, à chaque pas de temps,
on place le point représentant la bille dans le plan (Oxy).
Exercice 2. Programmer une animation qui permet de visualiser l’évolution de
la position dans le plan de la bille au cours du temps. Pour cela, on rappelle qu’il
faut faire une boucle en temps (parcourir des temps tn , n = 0, . . . , N , placer à
chaque pas de temps la position (x(θ(t)), y(θ(t))) en superposant les graphiques.
entrées : θ0 , θ1
sortie : film (x(θ(tn )), y(θ(tn ))) n = 0, . . . , N
3
4 Approximation numérique des solutions de (2) et (5)
Pour étudier et illustrer la dynamique du pendule simple, du fait de la non-linéarité
en sin, on est amené à faire une étude approchée des solutions de ce problème sur un
intervalle de temps [0, T ]. On définit pour cela une discrétisation temporelle uniforme de
pas h 0 = t0 < t1 < · · · < tN = T avec tn = nh. Pour résoudre numériquement l’équation
(3), on va programmer les méthodes d’Euler explicite et de Kunge-Kutta d’ordre 4 (notée
RK4).
Mise en place des méthodes numériques On rappelle que, pour l’approximation

numérique d’une EDO autonome
y 0 (t) = f (y(t)), t ≥ 0, y(0) = y0 ∈ Rd , f : Rd → Rd ,
les schémas numériques d’Euler explicite et de Runge-Kutta d’ordre 4 sont définis en

calculant l’approximation yn de y(tn ) par :
Euler explicite y0 = y0 , yn+1 = yn + hf (yn ).

RK4 y0 = y0 , yn+1 = yn + h2 [k1n + 2k2n + 2k3n + k4n ],
k1n = f (yn ),
k2n = f (yn + h2 k1n ),
k3n = f (yn + h2 k2n ),
k4n = f (yn + h2 k3n ).
Précision des méthodes : ordre de convergence Une étude théorique des deux
méthodes ci dessus (on définit pour cela les notions d’ordre de consistance et de stabilité
des méthodes numériques) permet de démontrer que, si la fonction f qui régit l’EDO
est suffisamment régulière, alors
max |yn − y(tn )| ≤ Chp sup y (p+1) (t) . (6)

0≤n≤N 0≤t≤T
Le plus petit tel entier p est appelé ordre de convergence de la méthode et l’on peut
démontrer que la méthode d’Euler explicite est d’ordre 1 et la méthode RK4 d’ordre 4.
Implémentation des méthodes numériques : résolution approchée du pro-

blème du pendule Pour simplifier votre première implémentation de méthode nu-
mérique pour les EDOs, on vous propose un programme pour implémenter la méthode
d’Euler explicite pour la résolution approchée du problème du pendule avec une donnée
initiale y0 = (θ0 , θ1 )T , un pas de temps h et un temps final T . On considérera pour
simplifier que g` = 1.
Créer un fichier eulerE .m et y reproduire la fonction eulerE suivante :
4
f u n c t i o n y=e u l e r E ( y0 , h , T)
y =[ y0 ] ; yaux=y0 ;
N=f l o o r (T/h ) ;
f o r i =1:N
yaux=yaux+h∗ f ( yaux ) ;
y =[ y yaux ] ;
end
qui nécessite bien sûr la création de la fonction f dans un fichier f .m :
f u n c t i o n y=f ( x )
y =[ x (2); − s i n ( x ( 1 ) ) ] ;
end
Exercice 3. Tester la fonction eulerE pour les données numériques suivantes : θ0 =
0, θ1 = 0.3, T = 10, h = 0.1. Que renvoie-t-elle ? Comment peut-on l’utiliser (s’inspirer
de l’exercice 2 pour illustrer la dynamique du pendule) ?
Exercice 4. Programmer une fonction RK4 qui implémente la méthode de RK4 appli-
quée au problème du pendule pour une donnée initiale y0 = (θ0 , θ1 ), un pas de temps h
et un temps final T . L’utiliser pour illustrer, comme dans l’exercice 2 la dynamique du
pendule.
Exercice 5. Tracer l’évolution de l’énergie mécanique numérique Em = 21 θ̇n2 + 1 − cos θn
en fonction du temps. Qu’observe-t-on ?
Exercice 6. Adapter les schémas numériques utilisés pour résoudre de manière appro-
chée le problème du pendule amorti (5).
5 Diagrammes d’erreur d’approximation pour l’équation

y 0 = −y
Pour illustrer la précision des méthodes utilisées, i.e. l’estimation (6), il est nécessaire
de s’attaquer à un problème dont on connaît la solution exacte (pour pouvoir tracer
|yn − y(tn )| où y(tn ) est la valeur de la solution exacte au problème au temps tn ). On
s’intéresse donc dans cette partie au problème très simple
∀t ≥ 0, y 0 (t) = −y(t), y(0) = y0 , (7)
dont la solution bien connue est bien sûr y(t) = y0 e−t .
Pour mettre en évidence l’ordre
de convergence des méthodes numériques implémentées, on va chercher à tracer l’erreur
Eh = max |yn − y(tn )| en fonction de h,
0≤n≤N
et on va chercher à illustrer le fait que, quand h devient petit, cette erreur se comporte
comme une puissance de h. Pour cela, on calcule Eh pour un certain nombre de pas h,
puis on trace log(Eh ) en fonction de | log h| et on s’attend à trouver une droite, dont on
calcule la pente (c’est-à-dire la puissance de h dans l’équation (6)) avec la commande
polyfit .

TPN 4 - Introduction À La Résolution Approchée D'Edos: 1 Modélisation: Le Pendule Simple

Transféré par

Droits d'auteur :

Formats disponibles

TPN 4 - Introduction À La Résolution Approchée D'Edos: 1 Modélisation: Le Pendule Simple

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

TPN 4 - Introduction À La Résolution Approchée D'Edos: 1 Modélisation: Le Pendule Simple

Transféré par

Droits d'auteur :

Formats disponibles

Faculté des sciences et ingénierie (Toulouse III) Année universitaire

Département de mathématiques – M1 MAF 2012-2013

1 Modélisation : le pendule simple

Le principe fondamental de la dynamique appliqué à la bille s’écrit :

~γ = `θ̈e~θ − `θ̇2 e~r .

P~ · e~θ = −mg sin θ et T~ · e~θ = 0.

2 Pistes pour l’étude théorique du pendule simple

f~frott = −Λ~v , (4)

~v représentant la vitesse de la bille exprimée dans le repère mobile par ~v = `θ̇e~θ .

3 Illustration avec Matlab de l’étude théorique du pendule

∀t ≥ 0, x(t) = ` sin θ(t), y(t) = −` cos θ(t).

Du point de vue de la programmation, on va représenter la trajectoire d’un point

Mise en place des méthodes numériques On rappelle que, pour l’approximation

y 0 (t) = f (y(t)), t ≥ 0, y(0) = y0 ∈ Rd , f : Rd → Rd ,

les schémas numériques d’Euler explicite et de Runge-Kutta d’ordre 4 sont définis en

Euler explicite y0 = y0 , yn+1 = yn + hf (yn ).

max |yn − y(tn )| ≤ Chp sup y (p+1) (t) . (6)

Implémentation des méthodes numériques : résolution approchée du pro-

5 Diagrammes d’erreur d’approximation pour l’équation

Vous aimerez peut-être aussi