Télécharge Exercices de mathématique 9 - La méthode d’Euler et plus Exercices au format PDF de Méthodes Mathématiques sur Docsity uniquement!
La méthode d’Euler : lancer d’un projectile
A. Fiche élève Le modèle est celui du canon : un obus est tiré dans l’axe du canon, suivant un angle avec l’horizontale. L’obus sort de la bouche du canon avec une vitesse initiale v 0 et suivant la tangente à la trajectoire. A un moment donné, les forces agissant sur l’obus sont la résistance de l’air R , proportionnelle au carré de la vitesse, dirigée à l’opposé de cette dernière (la vitesse est un vecteur représenté par la tangente à la trajectoire) et le poids P orienté vers le bas. Nous regardons tout d’abord ce qui se passe lorsqu’on ne tient pas compte des frottements. On a dans le repère orthonormal ( O ; i , j ) les coordonnées ( x ( t ) ; y ( t )) du centre de gravité M de l’obus, la vitesse v x ( '( ) ; t y '( )) t et l’accélération a x ( ''( ) ; t y ''( )) t. Les conditions initiales sont
M 0 O ( (0) ; (0)) x y 0 et v 0 ( x '(0) ; y '(0)) ( v 0 cos ; v 0 sin ).
Les équations du mouvement sont en suivant les lois de Newton : ''( ) 0 ''( )
x t ma R P ma mg a g (^) y t g
^
Toutes les représentations se feront avec un angle de 60° et une vitesse v 0 de 100 m.s−1. Ces valeurs
seront stockées dans les cellules B1 et B2 du tableur de manière à pouvoir être modifiées.
- On pose X t ( ) x '( ) t et Y t ( ) y '( ) t.
a. Vérifier que l’approximation affine du système ( 1 ) est (^ )^ ( ) ( ) ( )
X t h X t Y t h Y t gh
b. Quelles conditions initiales donneriez-vous? c. Représentez les solutions X ( t ) et Y ( t ) à l’aide du tableur (500 valeurs). Quelle est l’interprétation physique de ces solutions?
d. Calculez pour toutes les valeurs précédemment obtenues la quantité ( ) t X^2 ( ) t Y^2 ( ) t. Interprétez
les résultats physiquement.
- a. Donnez la solution exacte du système
X t Y t g
^
avec les conditions initiales^0 0
(0) cos (0) sin
X v Y v
Que vaut la quantité ?
b. Evaluez l’écart entre les solutions exactes et celles obtenues par la méthode d’Euler. Que pensez-vous du résultat?
- On cherche la trajectoire de l’obus. a. Montrez que le problème revient à résoudre '( ) ( ) '( ) ( )
x t X t y t Y t
^
où X et Y sont les solutions trouvées précédemment.
b. Vérifiez que l’approximation affine du système est
x t h x t hX t y t h y t hY t
^ ^ ^
. Quelles sont les conditions initiales?
M
R
P
v
v 0
O x
y
c. Représentez les solutions x ( t ) et y ( t ) à l’aide du tableur (500 valeurs). Quel est le type de courbe obtenue? d. On veut atteindre avec l’obus une cible située à 1 km de distance. Quelle doit-être l’inclinaison du canon?
- Vérification des solutions et de la trajectoire
a. Montrez que les solutions exactes du système sont
0 (^2 )
( ) cos ( ) 1 sin 2
x t v t y t gt v t
^
b. Tracez ces solutions sur la même représentation qu’au 3. c. Calculez l’écart entre les solutions exactes et celles du 3. Conclusion.
d. Vérifiez que la courbe décrite par l’obus a pour équation 2 2 2 0
(1 tan ) tan 2
y g x x v
e. Quelle doit-être l’inclinaison du canon pour atteindre avec l’obus une cible située à 1 km de distance?
B. Commentaires Cette méthode fait appel à la notion d’approximation linéaire d’une fonction en un point : si f est
dérivable en x 0 alors lim h 0^ f x^ (^0^^ h )^^ f x (^0^ ) f '( x 0 ) h ^ ^ ; en ne passant pas à la limite, mais en prenant h
suffisamment petit, on peut écrire^ f x^ (^0^^ h )^^ f x (^0^ )^ f '( x 0 (^) ) f x ( 0 (^) h ) f x ( 0 (^) ) hf '( x 0 ) h
^ ^ .
On dispose en général d’une relation du type f '( ) x ( , x f x ( )) (une équation différentielle), ce qui permet lorsqu’on connaît x et f ( x ) de calculer f ’( x ). Il s’agit alors de construire la courbe de la fonction inconnue f en utilisant ces relations : on part d’un premier point (^) x 0 , d’ordonnée connue (^) f x ( 0 ), on en déduit le point (^) x 1 (^) x 0 h d’ordonnée
f x ( 1 ) f x ( 0 h ) f x ( 0 ) hf '( x 0 ) f x ( 0 ) h ( x 0 , f x ( 0 )).
Il reste à recommencer de nombreuses fois pour se faire une idée de la courbe de la fonction f. Evidemment tout ceci reste approximatif, mais dans de nombreuses situations il est impossible ou très difficile de faire mieux. Au niveau physique la question est traitée en cours de physique, mais sans rentrer dans les détails. On ne fait pas intervenir les frottements ici, mais lorsque c’est le cas la situation devient beaucoup trop complexe. Il y a quand même le cas de la chute verticale où une solution analytique est à la portée des élèves. Si vous êtes intéressés, voir http://promenadesmaths.free.fr/chute_des_corps/Chute_libre.pdf dont je peux fournir le texte en word.
C. Pédagogie Barême : 1 point par partie (autant faire simple). Les résolutions analytiques sont vraiment très simples, mais ce n’est pas l’objectif ici. Par contre le fait de se frotter à un système du second ordre est fort intéressant et montre la puissance des méthodes numériques. Il me semble qu’il faut insister lourdement sur le contrôle de l’erreur, qui est ici possible grâce à la résolution exacte. Si on fait intervenir les frottements, il n’y a plus de solution exacte et le contrôle est évidemment beaucoup plus délicat.
