Mouvement central discret

Un article de Wikipédia, l'encyclopédie libre.

Vous avez de nouveaux messages (diff^?).

Cet article a une forme trop académique (extrait de cours) et nécessite d'être réécrit afin de davantage correspondre aux standards de l'encyclopédie.

Cet article, bien que respectant dans une certaine mesure la neutralité de point de vue, est jugé trop incomplet dans son développement ou dans l'expression des concepts et des idées. Son contenu est donc à considérer avec précaution. Pour toute information complémentaire, veuillez consulter sa page de discussion et la liste des articles incomplets.

Pour des personnes n'ayant aucun bagage en calcul différentiel et intégral, la mécanique peut quand même se comprendre en utilisant le calcul ombrial : c'est-à-dire transposer les règles différentielles en règles de différence, mais ajustées de manière adroite.

On peut ainsi faire comprendre le Mouvement à force centrale

[modifier] Mouvement newtonien discret

Le mouvement est plan. Le centre de force est le Soleil placé à l'origine O . Les coordonnées de la planète $P\,$ sont appelées $(x , y)\,$ . Le théorème de Pythagore donne $OP = r =\sqrt{x^2 + y^2}\,$ .

On notera les positions successives de la planète $P$ aux temps (ndt. $d t$ étant une durée appelé pas temporel discret du calcul : plus il est petit et meilleure est la trajectoire, mais plus les calculs sont longs), $P(n)\,$ de coordonnées $(x(n),y(n))\,$ .

On appelle vitesse de la planète selon l'axe $O_x\,$ , $u\,$ et selon l'axe $O_y\,$ , $v\,$ :

Le théorème de Pythagore donne $V =\sqrt{u^2+v^2}\,$ .

Mais attention, on appelle vitesse discrète $(u , v)\,$ entre le point $P(n)\,$ et le point $P(n+1)\,$ :

et idem pour $O_y\,$ :

L'accélération $(a, b)\,$ sera $a(n) = \frac{u(n+1)-u(n)}{dt}\,$ et selon $O_y\,$ : $b(n) = \frac{v(n+1)-v(n)}{dt}\,$ .

Le principe fondamental de la dynamique énoncé par Newton en 1687 consiste à dire : étant deux points initiaux $P(0)\,$ et $P(1)\,$ et la règle générale qui à partir de deux points permet de calculer le suivant , alors :

On peut en déduire $P(2)\,$
Puis à partir de $P(1)\,$ et $P(2)\,$ , calculer $P(3)\,$
Et continuer ainsi pour calculer tous les points de la trajectoire. Les ordinateurs de la Nasa ou de l'ESA ne procèdent pas autrement.

[modifier] Loi de force newtonienne discrète

Dans le cas du mouvement des planètes, cas discret, la loi est légèrement modifiée : la composante de l'accélération selon $O_x\,$ est :

et idem selon $O_y\,$ :

Grâce à cette loi de force, on peut démontrer par un calcul dit télescopique (les termes s'éliminent progressivement) que :

l'énergie cinétique s'écrit $T(n) = \frac{1}{2}\cdot m \cdot\frac{u(n)^2+v(n)^2}\,$
le travail de la force s'écrit $w(dt,n)= x(n+1) \cdot F_x(n)+y(n+1) \cdot F_y(n)\,$ , et le travail total $W(n)\,$ est $\sum_{k=0}^{n-1}w(k)\,$ .
On trouve alors la conservation de l'énergie mécanique :

En faisant tous les calculs, on trouve $W(n)=\frac{-k}{r(n)}\,$ tout simplement.

[modifier] Méthode d'itération

Là se pose le vrai problème :

En effet , il faudra résoudre le système de 4 équations (non linéaires) à 4 inconnues x(2),y(2),u(2) et v(2), qui est assez compliqué , mais que tout système type maple est capable de résoudre.

Il suffit alors de réitérer avec x(1),y(1),u(1),v(1) comme point de départ.

On vérifiera qu'au cours du calcul, l'énergie reste bien constante aux arrondis de calculs près. Les trajectoires sont au moins aussi belles que celles données par une méthode de Runge-Kutta d'ordre 4.

Remarque : comme pour tout calcul numérique, on n'opère jamais brutalement avec les valeurs réelles en S.I. qui n'a évidemment aucun sens pour un ordinateur. On opère donc en unité réduite, bâtie à partir du système d'unités naturelles des 4 équations récurrentes à étudier.