MENUCours de Mécanique classique

Ce chapitre présente une application importante des lois vues jusqu'ici : les mouvements à force centrale. Après quelques généralités, on aborde le problème de l'interaction newtonienne et notamment celui du mouvement des planètes qui fit le succès de la théorie de Newton.

Lois de conservation

On considère un point matériel M de masse \(m\) soumis à une force centrale conservative de centre O, point fixe d'un référentiel galiléen. Nous allons voir que le problème se résout grâce à deux relations de conservation.

Généralités

Rappelons qu'une force est dite centrale quand la droite support de cette force passe constamment par un point fixe O. Si l'on repère la position de M à l'aide d'un système de coordonnées sphériques d'origine O, on a \[\overrightarrow{f}=f(r,\theta,\varphi)\,\overrightarrow{u_r}\] La force est attractive quand \(f<0\), répulsive dans l'autre cas. Si la force centrale est conservative alors \(f(r,\theta,\varphi)\) ne peut dépendre que de \(r\). En effet, son travail élémentaire vaut \[\eth W=\overrightarrow{f}\cdot\overrightarrow{\mathrm{d}\ell}=f(r,\theta,\varphi)\,\mathrm{d} r\] forme différentielle qui doit être différentielle totale exacte. Ceci n'est possible que si \(f(r,\theta,\varphi)\) ne dépend que de \(r\). Dans ce cas, l'énergie potentielle associée ne dépend que de \(r\) et vérifie

Interaction gravitationnelle

La force de gravitation entre un astre fixe (massif) situé en O, à symétrie sphérique de masse \(m_{1}\) et un astre mobile à symétrie sphérique de masse \(m_{2}\) vaut, d'après la loi de gravitation universelle \[\overrightarrow{f_{12}}=-\frac{\mathcal{{G}}m_{1}m_{2}}{r^{2}}\,\overrightarrow{u_r} \quad\Longrightarrow\quad \mathrm{d} \mathcal{E}_\textrm{p}=\frac{\mathcal{{G}}m_{1}m_{2}}{r^{2}}dr\] d'où l'on tire l'énergie potentielle

Interaction coulombienne

La force électrostatique entre une charge ponctuelle fixe (charge électrique \(q_{1}\)) et une charge ponctuelle mobile (charge \(q_{2}\)) est une force centrale et s'écrit \[\overrightarrow{f_{12}}=\frac{q_{1}q_{2}}{4\pi\epsilon_{0}r^{2}}\,\overrightarrow{u_r}\] Lorsque \(q_{1}\) et \(q_{2}\) sont de même signe, la force est répulsive. C'est aussi une force conservative, d'énergie potentielle

Conservation du moment cinétique

Comme nous l'avons vu dans le chapitre précédent, un point matériel soumis à une force centrale dans un référentiel galiléen voit son moment cinétique \(\overrightarrow{L}_{\textrm{O}}\) se conserver : \[\frac{\mathrm{d} \overrightarrow{L_{\textrm{O}}}(\text{M})}{\mathrm{d} t} = r\,\overrightarrow{u_r}\wedge f(r)\,\overrightarrow{u_r}=\overrightarrow{0} \quad\Longrightarrow\quad \overrightarrow{L}_{\textrm{O}}(\text{M})=\overrightarrow{\mathrm{C^{te}}}\] ce qui a trois conséquences :

1. le mouvement est plan ;
2. en coordonnées polaires, la conservation du moment cinétique se traduit par la relation

   Loi des aires

   \begin{equation} r^{2}\dot{\theta}=C \end{equation} où \(C\) est la constante des aires déterminée par les conditions initiales ;

3. l'aire balayée par le vecteur \(\overrightarrow{\text{OM}}\) par unité de temps est constante et égale à \(C/2\) (vitesse aréolaire) : \[\frac{\mathrm{d} \mathcal{A}}{\mathrm{d} t}=\frac{C}{2}\]

Conservation de l'énergie mécanique

La force étant conservative, l'énergie mécanique se conserve : \[\frac{1}{2}mv^{2}+\mathcal{E}_\textrm{p}(r)=\mathcal{E}_\textrm{m}=\mathrm{C^{te}}\] En coordonnées polaires, le carré de la vitesse s'écrit, en utilisant l'équation (1), \[v^{2}=\dot{r}^{2}+(r\dot{\theta})^{2}=\dot{r}^{2}+\frac{C^{2}}{r^{2}}\] ce qui donne une intégrale première

Finalement la résolution de l'équation différentielle (2) donne \(r(t)\) puis \(\theta(t)\) en utilisant l'équation (1). On peut donc - en principe - trouver la trajectoire de M, soit de façon analytique quand l'équation différentielle est soluble, soit de façon numérique. Cependant, sans résoudre complètement le problème, on peut faire une analyse qualitative à partir du profil de \(\mathcal{E}_\textrm{p, eff}(r)\). En effet, la forme de l'énergie potentielle effective permet de savoir si la particule restera confinée autour du centre O ou s'en éloignera au bout d'un certain temps.

Par exemple, supposons que l'énergie potentielle effective ait l'allure représentée sur la figure ci-dessous.

Deux cas se présentent alors.

1. Si les conditions initiales sont représentées par un point P\(_{1}\) d'énergie \(\mathcal{E}_\textrm{m}\) situé dans un puits, alors \(r(t)\) oscille entre deux valeurs (cf. Oscillateurs mécaniques). Et comme \(r^{2}\dot{\theta}=C\), \(\theta(t)\) augmente ou diminue suivant le signe de \(C\). Ainsi, le point matériel décrit une orbite plus ou moins complexe autour du centre d'attraction. On dit que la particule est dans un état lié. Notons que si P\(_{1}\) possède une énergie mécanique correspondant à la valeur du fond du puits de potentiel, \(r\) reste constant au cours du temps ainsi que \(\dot \theta\). Le mouvement est alors circulaire uniforme.
2. Si les conditions initiales sont représentées par un point P\(_{2}\) situé en dehors du puits de potentiel, alors, après un éventuel rapprochement du centre d'attraction jusqu'à une distance d'approche minimale \(r_{\textrm{min}}\), le corps va s'éloigner indéfiniment du centre de force ; on dit que la particule est dans un état de diffusion ou état non lié.

Théorème de Bertrand

Lorsque le profil de \(\mathcal{E}_\textrm{p, eff}(r)\) présente un puits de potentiel, il existe une valeur \(r_{0}\) de \(r\) pour laquelle la trajectoire est circulaire. Cette valeur est donnée par le minimum de \(\mathcal{E}_\textrm{p, eff}(r)\) obtenue en annulant la fonction dérivée  :

Le mouvement est uniforme puisque \(r_{0}^{2}\dot\theta=C\) et la période orbitale vaut \[T_{\textrm{orb}}=\frac{2\pi}{\dot\theta}=\frac{2\pi r_{0}^{2}}{C}\] Imaginons maintenant que le corps en mouvement reçoive une petite quantité d'énergie de telle sorte que l'orbite circulaire est perturbée : la distance \(r\) va osciller autour de \(r_{0}\) avec une période d'oscillation \(T_{\textrm{osc}}\) donnée par la formule (cf. Oscillateurs mécaniques) :

Or, en général, \(T_{\textrm{orb}}\) et \(T_{\textrm{osc}}\) sont dans des rapports quelconques de sorte que la trajectoire, bien que liée, ne se referme pas sur elle même. La figure ci-dessous montre quelques exemples d'orbites.

Il existe cependant des cas où, quelles que soient les conditions initiales (en restant dans le cas où la trajectoire est bornée), ce rapport est commensurable : \(T_{\textrm{osc}}/T_{\textrm{orb}}=m/n\) avec \((m,n)\in\mathbb{N}^{2}\). En d'autres termes, après \(n\) révolutions, \(r(t)\) oscille exactement \(m\) fois ; la trajectoire se referme alors parfaitement. Analysons par exemple le cas de la force képlerienne \(f=-k/r^{2}\) de potentiel \(\mathcal{E}_\textrm{p}=-k/r\). L'équation (3) donne \[r_{0}=\frac{mC^{2}}{k} \quad\text{et}\quad T_{\textrm{orb}}=2\pi\frac{m^{2}C^{3}}{k^{2}}\] De la même manière, l'équation (4) donne \[T_{\textrm{osc}}=2\pi\sqrt{\frac{m}{\mathcal{E}''_\textrm{p, eff}}}= 2\pi\sqrt{\frac{m}{3m\,C^{2}/r_{0}^{4}-2k/r_{0}^{3}}}=2\pi\frac{m^{2}C^{3}}{k^{2}}\] Ainsi, on trouve \(T_{\textrm{osc}}=T_{\textrm{orb}}\). En d'autres termes, faiblement perturbée, l'orbite circulaire devient une orbite qui se referme après une révolution ; nous verrons qu'il s'agit d'une ellipse.

Joseph Bertrand se posa la question de la fermeture des trajectoires dans un cadre beaucoup moins restrictif que le notre (notre raisonnement n'est valable que pour de petites perturbations au voisinage de \(r_{0}\)) et trouva un résultat remarquable qu'il envoya à l'académie des sciences. Ce court article[1] fut publié en 1873 et montre le théorème suivant :

Le problème de Kepler

Le problème de Kepler fait référence à l'étude du mouvement d'une planète sphérique de masse \(m\) soumise à l'attraction gravitationnelle de la part d'une étoile de masse \(m_{\star}\) considérée fixe et à symétrie sphériqueTout écart à la sphéricité induit une force non centrale.. La seule force est centrale, newtonienne et attractive : \[\overrightarrow{f}=-\frac{\mathcal{G}m\,m_{\star}}{r^{2}}\,\overrightarrow{u_r}\] Bien entendu, tout corps de masse \(m\) lié gravitationnellement à un astre fixe de masse \(m'\) est soumis aux lois de Kepler ; il suffit alors de remplacer \(m\,m_{\star}\) par \(m\,m'\).

Lois de Kepler

Les trois lois sur le mouvement des astres du système solaire, dites lois de Kepler, s'énoncent ainsi :

1. Les planètes parcourent des orbites planes, elliptiques. Le Soleil occupe l'un des foyers de l'ellipse.
2. En des durées égales, les planètes balayent des aires égales.
3. Le rapport du carré de la période de rotation au cube du demi-grand axe est identique pour toutes les planètes du système solaire.

Ces trois lois permettront à Newton de confirmer sa théorie du mouvement des corps (Philosophiæ Naturalis Principia Mathematica 1687). Nous avons déjà vu comment la deuxième loi est une conséquence de la conservation du moment cinétique, voyons comment dériver les deux autres lois à partir des lois de Newton.

Cas du mouvement circulaire

Le profil de l'énergie potentielle effective \[\mathcal{E}_\textrm{p, eff}=\frac{mC^{2}}{2r^{2}}-\frac{\mathcal{G}mm_{\star}}{r}\] montre que les états liés sont possibles.

Si l'énergie mécanique correspond à la valeur minimale de \(\mathcal{E}_\textrm{p, eff}\), alors \(r\) reste constant au cours du temps et l'orbite est circulaire. On voit de plus que cette orbite circulaire est stable. Supposons donc que M décrit une orbite circulaire de rayon \(r_{0}\). En vertu de la conservation du moment cinétique, on a \[{r_{0}}^{2}\dot{\theta}=C\] ce qui implique que le mouvement est circulaire uniforme (vitesse angulaire constante). La relation fondamentale de la dynamique appliquée à M donne \[ma=m\frac{v^{2}}{r_{0}}=\frac{\mathcal{G}mm_{\star}}{r_{0}^{2}}\] Il vient alors \[ v=\sqrt{\frac{\mathcal{G}m_{\star}}{r_{0}}} \quad \text{et} \quad \mathcal{E}_\textrm{c}=-\frac{1}{2}\mathcal{E}_\textrm{p} \] la vitesse décroît quand \(r_{0}\) croît.

La période de révolution est reliée au rayon orbital par \(2\pi r_{0}=v\,T\) puisque le mouvement est uniforme. Il vient, en élevant au carré : \[ \frac{r_{0}^{3}}{T^{2}}=\frac{\mathcal{G}m_{\star}}{4\pi^{2}} \] relation qui donne une version simplifiée de la troisième loi de Kepler.

Enfin l'énergie mécanique est constante et égale à \[\mathcal{E}_\textrm{m}=\mathcal{E}_\textrm{c}+\mathcal{E}_\textrm{p}=\frac{\mathcal{E}_\textrm{p}}{2}=-\frac{\mathcal{G}mm_{\star}}{2r_{0}}\]

Ensemble des trajectoires solutions

L'orbite circulaire est donc une solution stable particulière. L'ensemble des trajectoires possibles s'obtient à partir des deux relations de conservation : \[\left\{\begin{array}{rcl} \dfrac{1}{2}m\dot{r}^{2}+\dfrac{mC^{2}}{2r^{2}}-\dfrac{\mathcal{G}mm_{\star}}{r} &=& \mathcal{E}_\textrm{m}\\[4mm] r^{2}\dot{\theta} &=& C \end{array}\right.\] Il est alors judicieux de procéder au changement de variable \(u=\frac{1}{r}\) afin d'obtenir l'équation polaire de la trajectoire, à savoir la relation \(r(\theta)\). On a \[\dot{r}=\frac{\mathrm{d} \left(1/u\right)}{\mathrm{d} t}= -\frac{1}{u^{2}}\frac{\mathrm{d} u}{\mathrm{d} t}=-\frac{1}{u^{2}}\frac{\mathrm{d} u}{\mathrm{d} \theta}\dot{\theta}= -C\frac{\mathrm{d} u}{\mathrm{d} \theta}\] Le système d'équations devient : \[\left\{\begin{array}{rcl} \dfrac{1}{2}mC^{2}\left(\dfrac{\mathrm{d}u}{\mathrm{d}\theta}\right)^{2}+\dfrac{mC^{2}u^{2}}{2}-\mathcal{G}mm_{\star}u &=&\mathcal{E}_\textrm{m}\\ \dot{\theta} &=& Cu^{2} \end{array}\right.\] Si l'on dérive la première relation par rapport à \(\theta\) il vient \[\frac{\mathrm{d}^{2}u}{\mathrm{d} \theta^{2}}+u=\frac{\mathcal{G}m_{\star}}{C^{2}}\] équation différentielle linéaire avec un second membre constant dont la solution s'écrit \[u(\theta)=\frac{\mathcal{G}m_{\star}}{C^{2}}+A\cos(\theta-\theta_{0})\] c'est-à-dire

On reconnaît l'équation polaire \(r(\theta)\) d'une conique de paramètres \(p\) et \(e\) dont l'origine O est l'un des foyers et dont l'axe focal est la droite \(\theta=\theta_{0}\) (cf. les coniques). La trajectoire est donc une conique de foyer le centre d'attraction et dont la forme dépend des conditions initiales puisque \(p\) et \(e\) en dépendent. Suivant la valeur de \(e\), on aura des états liés (l'ellipse pour \(e<1\)) ou des états de diffusion (la parabole pour \(e=1\) ou l'hyperbole pour \(e > 1\)). Les planètes, confinées autour du Soleil, décrivent ainsi une ellipse de foyer le centre du Soleil.

On peut noter que l'équation (5) ne permet pas d'avoir l'équation horaire de l'astre. Il faudrait pour cela connaître l'évolution de l'angle \(\theta\) au cours du temps. Cependant, la loi des aires \(\dot{\theta}=C/r^{2}\) se met sous la forme \(\dot{\theta}=f(\theta)\) — équation différentielle ordinaire du premier ordre — qu'il est toujours possible de résoudre par une approche numérique.

Troisième loi de Képler

La troisième loi est une conséquence des deux premières et ne s'applique que pour les mouvements elliptiques. Tout d'abord la loi des aires se traduit par \[\frac{\mathrm{d} \mathcal{A}}{\mathrm{d} t}=\frac{C}{2} \quad\Longrightarrow\quad \frac{\mathcal{A}}{T}=\frac{C}{2}\] où \(T\) est la période orbitale et \(\mathcal{A}\) l'aire de l'ellipse. Les mathématiques donnent \(\mathcal{A}=\pi ab\) où \(a\) et \(b\) sont le grand et le petit-axe de l'ellipse. De plus, \(b=\sqrt{ap}\) (cf. les coniques) de sorte que si l'on élève au carré la loi des aires, on obtient \[\frac{\pi^{2}a^{3}p}{T^{2}}=\frac{C^{2}}{4}\] Or, selon la solution générale (5) \(p=C^{2}/\mathcal{G}m_{\star}\) d'où l'on déduit la fameuse formule

En d'autres termes, le rapport \(a^{3}/T^{2}\) est identique pour toutes les planètes du système solaire comme l'avait constaté Kepler. La contribution de Newton est de relier cette loi à la constante de gravitation universelle et la masse du Soleil.

Energie

Allons au delà des lois de Kepler et montrons que l'énergie mécanique s'exprime simplement en fonction du grand-axe de la conique.

Reprenons les expressions de \(r\) et de \(\dot{r}\) (par souci de simplification on choisit l'axe O\(x\) de façon à ce que \(\theta_{0}=0\)). \[\left\{\begin{array}{rcl} r &= &\dfrac{p}{1+e\cos\theta} \\[3mm] \dot{r} &= &-C\dfrac{\mathrm{d} u}{\mathrm{d} \theta} = \dfrac{Ce\sin\theta}{p} \end{array}\right.\] pour les substituer dans l'énergie mécanique. Il vient alors \[\begin{array}{rcl} \mathcal{E}_\textrm{m} &= &\dfrac{1}{2}m\dot{r}^{2}+\dfrac{mC^{2}}{2r^{2}}-\dfrac{\mathcal{G}mm_{\star}}{r}\\[3mm] &= &\dfrac{1}{2}mC^{2}\left[\dfrac{(1+e\cos\theta)^{2}}{p^{2}}+\dfrac{e^{2}\sin^{2}\theta}{p^{2}}\right]-\dfrac{\mathcal{G}mm_{\star}(1+e\cos\theta)}{p}\\[3mm] \mathcal{E}_\textrm{m} &= & \dfrac{mC^{2}}{2p^{2}}\left(1+e^{2}+2e\cos\theta\right)-\dfrac{\mathcal{G}mm_{\star}(1+e\cos\theta)}{p}\\ \end{array}\] Or, \(C^{2}=\mathcal{G}m_{\star}p\) d'où \[ \mathcal{E}_\textrm{m}=-\frac{\mathcal{G}mm_{\star}}{2p}\left(1-e^{2}\right) \] On peut alors distinguer trois cas.

1. Pour une hyperbole, \(e>1\) et \(\mathcal{E}_\textrm{m}>0\). Les états ne sont pas liés et le corps s'éloigne indéfiniment du centre d'attraction avec une énergie cinétique non nulle lorsque \(r\to\infty\). Le grand axe d'une hyperbole vaut \(a=p/(e^{2}-1)\) d'où \[\mathcal{E}_\textrm{m}=\frac{\mathcal{G}mm_{\star}}{2a}\quad\text{[hyperbole]}\]
2. Pour une parabole, \(e=1\) et \(\mathcal{E}_\textrm{m}=0\), ce qui signifie que le corps va s'éloigner du centre avec une vitesse qui tend vers 0 lorsque \(r\to\infty\).
3. Pour une ellipse, \(e<1\) et \(\mathcal{E}_\textrm{m}<0\) : le corps est lié à l'astre central. Pour une ellipse, le grand axe vaut \(a=p/(1-e^{2})\) ce qui donne \[\mathcal{E}_\textrm{m}=-\frac{\mathcal{G}mm_{\star}}{2a}\quad\text{[ellipse]}\]

On pourra retenir que dans le cadre du problème de Kepler, l'énergie mécanique d'une planète en orbite autour d'une étoile de masse \(m_{\star}\) vaut \[ \fcolorbox{#FF9D00}{#FEF5EB}{\quad \(\displaystyle \mathcal{E}_\textrm{m}=-\frac{\mathcal{G}m\,m_{\star}}{2a} \)\quad} \;\color{#FF9900}{\heartsuit} \]

Vitesses cosmiques

Supposons que l'on veuille envoyer un satellite S de masse \(m\) en orbite autour de la Terre. Pour cela, on utilise un lanceur qui, dans une phase propulsive, donne l'énergie suffisante au satellite. Lorsque le satellite se détache du lanceur à la hauteur \(h\), il entre dans la phase balistique avec une certaine vitesse \(v\) et une énergie \(\mathcal{E}_\textrm{m}\). Si l'on assimile la Terre à un astre à symétrie sphérique de masse \(m_{\textrm{T}}\) et de rayon \(R_{\textrm{T}}\), on a (L'altitude \(h\) est suffisamment basse pour avoir \(h\ll R\) et suffisamment haute pour négliger les frottements dus à une atmosphère très raréfiée) : \[\mathcal{E}_\textrm{m}=\frac{1}{2}mv^{2}-\frac{\mathcal{G}m\,m_{\textrm{T}}}{R_{\textrm{T}}+h}\simeq\frac{1}{2}mv^{2}-\frac{\mathcal{G}m\,m_{\textrm{T}}}{R_{\textrm{T}}}\] Deux cas limites se présentent.

• \(\mathcal{E}_\textrm{m}\geq 0\) : la trajectoire n'est pas liée et le satellite s'éloigne indéfiniment de la Terre Plus exactement, elle quitte l'attraction de la Terre pour rejoindre la «sphère d'influence» d'un autre astre.. Dans ce cas, on a \[v\geq v_{\textrm{lib}}=\sqrt{2gR_{\textrm{T}}}\approx 11\;\mathrm{km.s^{-1}}\] Ce qui définit la vitesse de libération \(v_{\textrm{lib}}\).
• \(\mathcal{E}_\textrm{m}<0\) : la trajectoire est une ellipse de foyer le centre de la Terre. Pour éviter une collision avec la terre, le périgée de l'orbite doit se situer à une distance \(r_{\textrm{min}}=a(1-e)>R_{\textrm{T}}\). Il faut donc communiquer une énergie \[\mathcal{E}_\textrm{m}=-\frac{\mathcal{G}m\,m_{\textrm{T}}}{2a}> -\frac{\mathcal{G}m\,m_{\textrm{T}}}{2R_{\textrm{T}}}(1-e)\] Dans le cas d'une orbite circulaire - cas où l'énergie minimale est la plus faible - on a \[\mathcal{E}_\textrm{m} >-\frac{\mathcal{G}m\,m_{\textrm{T}}}{2R_{\textrm{T}}}\] Condition qui, traduit en termes de vitesse, donne \[v\geq v_{\textrm{sat}}=\sqrt{gR_{\textrm{T}}}\approx 8\;\mathrm{km.s^{-1}}\] ce qui définit la vitesse de satellisation \(v_{\textrm{sat}}\).

Notez que ces vitesses ne dépendent pas de la masse. Elles concernent aussi bien les astres que les molécules.

Equation horaire

Le problème de Kepler n'est pas complètement résolu au sens où nous n'avons trouvé que l'ensemble des trajectoires. Il nous reste à préciser le mouvement de M le long de cette trajectoire, c'est-à-dire à trouver la relation entre \(\theta\) et le temps \(t\). Il est commode de fixer l'origine des temps lorsque M se situe au péricentre (\(\theta(0)=\theta_0\)).

Dans la cas d'une orbite fermée de période \(T\) on définit la vitesse angulaire moyenne \[\overline{\omega}=\frac{2\pi}{T}\] reliée à la loi de Kepler par \(\overline{\omega}^2 a^{3}=\mathcal{G}m_{\star}\).

Cas du cercle

Lorsque l'orbite est circulaire de rayon \(r_{0}\) nous avons montré que M décrit l'orbite à une vitesse angulaire constante \(\dot{\theta}=\overline{\omega}\). Ainsi, le mouvement de l'astre est régi par l'équation horaire \[\left\{\begin{array}{rcl} r &=& r_{0} \\ \theta &=& \overline{\omega} t+\theta_0 \end{array}\right.\]

Cas des faibles excentricités

Comme le montre le tableau, les planètes du système solaire ont une trajectoire elliptique de faible excentricité. On peut dans ce cas, approcher l'équation horaire par un développement d'ordre 1 en \(e\). Pour trouver la relation \(\theta(t)\), il faut résoudre l'équation différentielle du premier ordre \[r^{2}\dot{\theta}=C \quad\text{avec}\quad r=\frac{p}{1+e\cos\theta} \quad\text{et}\quad p=\frac{C^{2}}{\mathcal{G}m_{\star}}\] où on a pris soin de choisir l'axe O\(x\) de façon à avoir \(\theta_{0}=0\). L'équation différentielle se résout par la méthode de séparation des variables :

Intégrale que l'on approche à l'aide d'un développement limité à l'ordre 1 de l'intégrant : \[(1+e\cos\theta')^{-2}\simeq 1-2e\cos\theta' \qquad\text{d'où}\qquad p^{2}(\theta-2e\sin\theta)\simeq C\,t\] De plus, \[C^{2}=p\,\mathcal{G}m_{\star}=p\,\overline{\omega}^{2}a^{3}=\frac{\overline{\omega}^{2}p^{4}}{(1-e^{2})^{3}}\] Finalement, en négligeant \(e^{2}\) devant 1, on trouve \(C/p^{2}\simeq \overline{\omega}\) et l'équation horaire s'écrit \[ \theta-2e\sin\theta=\overline{\omega}t \] Ainsi, le mouvement des planètes du système solaire peut être décrit par le système d'équations \[\left\{\begin{array}{rcl} r &=& \dfrac{p}{1+e\cos\theta} \\[3mm] \theta-2e\sin\theta &=&\overline{\omega}t \end{array}\right.\]

On peut poursuivre le développement à l'ordre 1 en \(e\), en remarquant que \(\theta=2e\sin\theta+\overline{\omega}t\) et donc que \[2e\sin\theta=2e\sin(2e\sin\theta+\overline{\omega}t)=2e\sin(2e\sin\theta)\cos(\overline{\omega}t)+2e\sin(\overline{\omega}t)\cos(2e\sin\theta)\] Si l'on ne garde que les termes d'ordre 1, on trouve \[2e\sin\theta=2e\sin(\overline{\omega}t)\] En conclusion, si l'on néglige les termes d'ordre supérieur ou égal à deux, on trouve \[\left\{\begin{array}{rcl} r &=& \dfrac{p}{1+e\cos\theta} \\[3mm] \theta &=&\overline{\omega}t+2e\sin(\overline{\omega}t) \end{array}\right.\]

Cas elliptique

Sans se restreindre aux petites excentricités, on peut obtenir l'équation horaire au prix de développements assez calculatoires. On obtient : \[\left\{\begin{array}{rcl} r &=& \dfrac{p}{1+e\cos\theta} \\[3mm] \tan\theta/2 &=& \sqrt{\dfrac{1+e}{1-e}}\tan E/2\\[3mm] E-e\sin E &=& \overline{\omega} t \end{array}\right.\] La variable \(E\) désigne ici l'anomalie excentrique. Elle se confond avec \(\theta\) lorsque \(e\to 0\). Pour obtenir la position du corps à chaque instant, il est nécessaire de résoudre l'équation transcendante \(E-e\sin E=\overline{\omega}t\) ce qui peut être réalisé grâce à une méthode numérique.

Cas parabolique

Analysons le cas d'un astre sur une orbite parabolique arrivant au péricentre à la vitesse \(v_{\textrm{max}}\). Dans le cas, \(e=1\) et l'équation (6) donne \[\int_{0}^{\theta}\dfrac{p^{2}}{(1+\cos\theta')^{2}}\,\mathrm{d}\theta'= \int_{0}^{\theta}\dfrac{p^{2}}{4\cos^4\theta'/2}\,\mathrm{d} \theta'=Ct\] Il est facile de vérifier que la primitive de \(1/\cos^{4}x\) est \(\tan x+1/3\tan^{3}x\) d'où : \[\tan \theta/2+\frac{1}{3}\tan^{3}\theta/2=\frac{2C}{p^{2}}t\] Lorsque M atteint le péricentre, il se trouve à une distance minimum \(r_{\textrm{min}}=r(\theta=0)=p/2\) avec une vitesse maximum \(\overrightarrow{v}_{\textrm{max}}\) perpendiculaire au vecteur position d'où \[C=r_{\textrm{min}}v_{\textrm{max}}=\frac{p v_{\textrm{max}}}{2}\] Finalement l'équation horaire d'une orbite parabolique de paramètre \(p\) s'écrit \[\left\{\begin{array}{rcl} r &=& \dfrac{p}{1+\cos\theta} \\[3mm] \tan\theta/2+\dfrac{1}{3}\tan^{3}\theta/2 &=&\dfrac{v_{\textrm{max}}}{p} t \end{array}\right.\]

Interaction coulombienne

On se place maintenant dans le cas où une charge ponctuelle fixe \(q_{1}\) interagit avec une charge mobile \(q_{2}\) via une force centrale attractive ou répulsive selon le signe du produit des charges. \[\overrightarrow{f}=\frac{q_{1}q_{2}}{4\pi\epsilon_{0}r^{2}}\,\overrightarrow{u_r}\]

Cas attractif

Considérons dans un premier temps le cas ou les deux charges sont opposées. Par exemple le problème de l'électron lié électriquement à un proton 1836 fois plus lourd et donc quasiment fixe, peut être un point de départ pour modéliser l'atome d'hydrogène, si l'on accepte de le décrire dans un cadre newtonien Rigoureusement, la physique de l'atome d'hydrogène et des autres éléments de la classification périodique relève de la mécanique quantique..

Il est clair que le problème mathématique est formellement identique au problème de Kepler, il suffit simplement d'opérer le changement suivant : \[ -\mathcal{G}m\,m_{\star} \leftrightarrow \frac{q_{1}q_{2}}{4\pi\epsilon_{0}} \] La particule de charge \(q_{2}\) décrit donc une conique dont un des foyers coïncide avec la charge \(q_{1}\) et dont l'équation est donnée par \[r=\frac{p}{1+e\cos(\theta-\theta_{0})} \quad\text{avec}\quad \left\{\begin{array}{rcl} p&=&-\dfrac{4\pi\epsilon_{0}m\,C^{2}}{q_{1}q_{2}} \\ e&\geq &0 \\ \end{array}\right.\] et l'énergie vaut \[\mathcal{E}_\textrm{m}=\frac{q_{1}q_{2}}{8\pi\epsilon_{0}p}\left(1-e^{2}\right)\]

Cas répulsif

Plaçons nous dans le cas où les deux charges sont de même signe : \(q_{1}q_{2}>0\). L'énergie potentielle effective \[\mathcal{E}_\textrm{p, eff}=\frac{mC^{2}}{2r^{2}}+\frac{q_{1}q_{2}}{4\pi\epsilon_{0}r}\] étant une fonction monotone décroissante montre que les états liés sont impossibles : la particule chargée, après une éventuelle approche du centre répulsif, va nécessairement s'en éloigner définitivement. En reprenant la démarche utilisée pour le problème de Kepler, on trouve que la solution générale s'écrit \[r=\frac{p}{e\cos(\theta-\theta_{0})-1} \quad\text{avec}\quad \left\{\begin{array}{rcl} p &=& \dfrac{4\pi\epsilon_{0}m\,C^{2}}{q_{1}q_{2}} \\ e & > &1 \end{array}\right.\] Ici \(r>0\) implique \(e>1\) : la trajectoire est donc une branche d'hyperbole de paramètres \(p\), d'excentricité \(e\) et de foyer, le centre répulsif. \(\theta=\theta_{0}\) correspond à l'axe focal de la branche d'hyperbole. Si l'on fixe l'origine des angles de telle sorte que \(\theta_{0}=0\), on a \[r=\frac{p}{e\cos\theta-1}\] La différence essentielle par rapport au cas d'un potentiel newtonien attractif est que le point matériel décrit une branche d'hyperbole qui ne contourne pas le centre répulsif.

Ainsi, une particule chargée d'énergie cinétique initiale \(\mathcal{E}_\textrm{c}\) lancée en direction d'une charge fixe, va subir une déflexionOn parle de diffusion de Rutherford. du fait de la répulsion. La trajectoire de la particule subit une déviation angulaire \(\Theta_{\textrm{dev}}\) qu'il est possible d'exprimer en fonction de \(\mathcal{E}_\textrm{c}\), des charges en interaction et du demi-petit axe \(b\) que l'on appelle ici le paramètre d'impact. Tout d'abord, l'angle \(\Theta_{\textrm{dev}}\) est lié à l'angle \(\pm\theta_{1}\) que font les asymptotes. On a \[\Theta_{\textrm{dev}}=\pi-2\theta_{1} \quad\Longrightarrow\quad \tan\theta_{1}=\frac{1}{\tan(\Theta_{\textrm{dev}}/2)}\] Or, on sait[7] que \(\tan\theta_{1}=\sqrt{e^{2}-1}\). Par ailleurs, la constante des aires vaut \[C=\Vert\overrightarrow{r}\wedge\overrightarrow{v}\Vert_{r\to\infty}=b\,v_{\infty} \quad\Longrightarrow\quad p= \frac{4\pi\epsilon_{0}m\,b^{2}v_{\infty}^{2}}{q_{1}q_{2}}\] Enfin, \(p=a(e^{2}-1)\) et \(b=a\sqrt{e^{2}-1}\) de telle sorte que \[\tan\theta_{1}=\frac{p}{b}=\frac{4\pi\epsilon_{0}m\,b\,v_{\infty}^{2}}{q_{1}q_{2}}\] La loi de diffusion de Rutherford est donc donnée par \[ \tan\Theta_{\textrm{dev}}/2=\frac{q_{1}q_{2}}{8\pi\epsilon_{0}\mathcal{E}_\textrm{c}b} \] Les particules sont d'autant plus diffusées que le paramètre d'impact est faible. Lorsque \(b=0\), la particule arrive de façon frontale et perd toute son énergie cinétique jusqu'à atteindre la distance minimale d'approche \(d_{\textrm{min}}\) puis rebrousse chemin, l'énergie potentielle étant convertie en énergie cinétique. On a alors ce que l'on appelle une rétrodiffusion. La distance minimale d'approche est donnée par l'égalité entre l'énergie cinétique initiale et l'énergie potentielle en ce point : \[\mathcal{E}_\textrm{c}=\frac{q_{1}q_{2}}{4\pi\epsilon_{0}d_{\textrm{min}}} \quad\Longrightarrow\quad d_{\textrm{min}}=\frac{q_{1}q_{2}}{4\pi\epsilon_{0}\mathcal{E}_\textrm{c}}\]