flèche de retour

Licence 1 > Mécanique 2 > Cours 22 : Forces centrales

M22 : Forces centrales

Les notions de ce chapitre sont disponibles sous la forme de 4 vidéos que vous retrouverez dans la playlist adéquate sur la chaine youtube de physagreg :

Voici également les quatre séquence de façon distincte :

  1. Définition d'une force centrale, constance du moment cinétique et conséquences
  2. Energie mécanique et énergie potentielle effective
  3. Energie effective (force attractive), équation polaire de la trajectoire, vitesse de libération
  4. Trajectoire circulaire et elliptique

Introduction

Ce chapitre va être l’occasion de revoir deux forces que l’on connaît bien, la force gravitationnelle (dite de Newton) et la force électrostatique (dite de Coulomb). En effet, nous l’avons déjà dit, ces forces présentent des similitudes, notamment leur variation en \(\dfrac{1}{r^{2}}\).

Dans ce chapitre, nous verrons les forces centrales conservatives, dont la force de Newton et celle de Coulomb font parties, et leurs caractéristiques ; puis nous étudierons le mouvement d’un point M soumis à une force centrale en remarquant la constance de certaines grandeurs.

Forces centrales conservatives

Définition

Une force centrale est une force qui s’écrit \(\overrightarrow{F} = F(r) \overrightarrow{u}_r\) en coordonnées sphériques.
Cela signifie :

  • que sa valeur ne dépend pas du temps ;

  • que a valeur de dépend que de r, la distance de M (point qui subit la force) à O (point appelé centre de force) ;

  • que sa droite d’action a la même direction que le vecteur \(\overrightarrow{OM}\).

Force centrale
M subit une force centrale de centre O

Cette force est conservative (le calcul de son travail ne dépend pas du chemin suivi), elle dérive donc d’une énergie potentielle :

\begin{equation}\boxed{\overrightarrow{F} = -\overrightarrow{\mathrm{grad}}\, E_P \quad \text{et ainsi} \quad F(r) = -\dfrac{\mathrm{d}E_P}{\mathrm{d}r}}\end{equation}

Exemples

Mouvement général d’un point M soumis à une force centrale conservative

Nous allons voir que cette notion de force centrale a des conséquences quant à la conservation de certaines grandeurs physiques, que l’on peut traduire en terme de mouvement.

Moment cinétique

Nous allons montrer que le fait que le point M ne soit soumis qu’à une force centrale rend son moment cinétique constant.
Appliquons le théorème du moment cinétique en O dans le référentiel galiléen (O,\(\overrightarrow{u}_x\),\(\overrightarrow{u}_y\),\(\overrightarrow{u}_z\)) :

\begin{equation}\dfrac{d\overrightarrow{L_O}(M)}{dt} = \overrightarrow{\mathcal{M}_O}(\overrightarrow{F}) = \overrightarrow{OM}\wedge\overrightarrow{F} = r\overrightarrow{u}_r\wedge F(r)\overrightarrow{u}_r = \overrightarrow{0}\end{equation}

Le fait que le moment cinétique soit constant à deux conséquences :

➲ Trouvons tout d’abord l’expression de la constante des aires $C$, liée au moment cinétique, en exprimant le moment cinétique en coordonnées cylindriques :

\begin{equation}\overrightarrow{L_O}(M)= m\,\overrightarrow{OM} \wedge \overrightarrow{v}(M) = m\,r\,\overrightarrow{u}_r \wedge (\dot{r}\overrightarrow{u}_r + r \dot{\theta} \overrightarrow{u}_{\theta})=mr^2\dot{\theta}\,\overrightarrow{u}_z\end{equation}

On note généralement \(\overrightarrow{L_O}(M) = m\,C\,\overrightarrow{u}_z\) avec \(\boxed{C=r^2\dot{\theta}}\,\).

➲ Exprimons maintenant l’aire balayée par le rayon vecteur pendant un temps dt :

\begin{equation}\mathrm{d}\mathcal{A} = \frac{1}{2}\,OM\times v\,\mathrm{d}t = \frac{1}{2}\,r\times r\dot{\theta}\,dt\end{equation}

car on peut voir cette portion infinitésimale d’aire comme un triangle de hauteur \(r\) et de base \(v\,\mathrm{d}t\).

Ainsi :

\begin{equation}\dfrac{\mathrm{d}\mathcal{A}}{\mathrm{d}t} = \dfrac{r^2\dot{\theta}}{2} = \dfrac{C}{2}\end{equation}

Donc :

\begin{equation}\boxed{\mathcal{A}(t) = \dfrac{C}{2}\,t + \mathrm{cste}}\end{equation}
Loi des aires
Aire balayée par OM pendant $\mathrm{d}t$

Remarque
La grandeur \(\dfrac{\mathrm{d}\mathcal{A}}{\mathrm{d}t}\) se nomme parfois vitesse aréolaire, vitesse de balayage d’une aire.

Énergie mécanique

Conservation de l’énergie mécanique

Le fait que la force centrale soit conservative implique que l’énergie mécanique du point M est conservée au cours du mouvement.
En effet, d’après le théorème de l’énergie cinétique, pour le point M qui se déplace entre la position A et la position B :

\begin{equation}E_C(B) - E_C(A) = W_{AB}(\overrightarrow{F})\end{equation}

Or comme la force \(\overrightarrow{F}\) est conservative, on peut définir une énergie potentielle telle que :

\begin{equation}W_{AB}(\overrightarrow{F}) = E_P(A) - E_P(B)\end{equation}

En effet :

\begin{equation} W_{AB}(\overrightarrow{F}) = \int_A^B F(r)\,\mathrm{d}r = \int_A^B \dfrac{K}{r^2}\,\mathrm{d}r = \left[-\dfrac{K}{r}\right]_A^B = E_P(A) - E_P(B) \end{equation}

Donc en deux positions quelconques A et B du point M :

\begin{equation}E_C(A) + E_P(A) = E_C(B) + E_P(B) = \mathrm{cste}\end{equation}

On peut donc écrire :

\begin{equation}E_\mathrm{M} = E_C + E_P = \mathrm{cste}\end{equation}

Si cette énergie est constante c’est qu’elle a à tout instant la valeur qu’elle avait dans l’instant initial : on peut donc déterminer sa valeur à partir des conditions initiales.

Définition d’une énergie potentielle effective

On peut exprimer cette énergie mécanique en fonction de l’unique variable r. On a :

\begin{equation}E_\mathrm{M} = \frac{1}{2}\,m\,v^2 + E_P(r)\end{equation}

Car nous avons dit que l’énergie potentielle ne dépendait que de r \(\left(E_P = \dfrac{K}{r} + \mathrm{cste}\right)\).

On peut exprimer la vitesse en coordonnées polaires :

\begin{equation} E_\mathrm{M} = \frac{1}{2}m(\dot{r}^2+r^2\dot{\theta}^2) + E_P(r) \end{equation}

Car \(\overrightarrow{v}=\dot{r}\overrightarrow{u}_r + r\dot{\theta} \overrightarrow{u}_{\theta}\), si on porte cette expression au carré, le terme "$2ab$" fait apparaître le produit scalaire \(\overrightarrow{u}_r\cdot\overrightarrow{u}_{\theta}\) qui est nul.

Enfin on peut faire apparaître la constante des aires et ainsi définir une nouvelle énergie potentielle :

\begin{equation}\boxed{E_\mathrm{M} = \frac{1}{2}m\dot{r}^2 +\dfrac{m\,C^2}{2\,r^2} + E_P(r) = \frac{1}{2}\,m\,\dot{r}^2 + E_{\mathrm{Peff}}(r)}\end{equation}

\(E_{\mathrm{Peff}}(r)= E_P(r) + \dfrac{m\,C^2}{2\,r^2}\) est appelée énergie potentielle effective, elle comprend l’énergie potentielle et une partie de l’énergie cinétique du point M.

Le problème se résume alors à l’étude d’un point matériel de masse \(m\) dont la position est décrite par un seul degré de liberté, \(r\) ; et soumis à une force conservative dont l’énergie potentielle est \(E_\mathrm{p\mathrm{eff}}\).

Pourquoi définir une énergie potentielle effective ?

Cette énergie, qui n’a pas réellement de sens physique, va permettre par son étude, de trouver les formes de mouvements possibles pour le point M en fonction du signe de $K$ et de la valeur de la constante \(E_\mathrm{M}\).

Étude des mouvements possibles et conditions d’existence

Nous allons utiliser l’énergie potentielle effective définie précédemment pour identifier les mouvements possibles. Selon l’allure de \(E_{\mathrm{Peff}}=f(r)\) et la valeur de l’énergie mécanique de M, nous distinguons plusieurs cas :

Cas d’une force répulsive (\(K>0\))

La fonction \(E_{\mathrm{Peff}}=f(r)\) a l’allure suivante :

Profil d'énergie potentielle effective, cas du K>0
Allure de \(E_{\mathrm{Peff}} = f(r)\) dans le cas d’une force répulsive

Rappelons la relation donnant l’énergie mécanique : \(E_\mathrm{M} = \frac{1}{2}\,m\,\dot{r}^2 + E_{\mathrm{Peff}}(r)\)
Ainsi, comme l’énergie cinétique radiale est nécessairement positive, pour qu’il y ait mouvement il faut que \(E_\mathrm{M} \geq E_{\mathrm{Peff}}\).

Le mouvement ne peut donc se faire qu’entre \(r = r_{\mathrm{min}}\) et l’\(\infty\), ce mouvement est non borné :
Lorsque la force centrale est répulsive, on parle d’un état de diffusion.

Pour trouver la valeur de \(r_{\mathrm{min}}\), il faut résoudre l’équation :

\begin{equation}E_\mathrm{M} = E_{\mathrm{Peff}} \Longleftrightarrow E_\mathrm{M} - \dfrac{m\,C^2}{2\,r^2} - \dfrac{K}{r} = 0\end{equation}

Ce qui équivaut à l’équation du second degré suivante :

\begin{equation}E_\mathrm{M}\,r^2 - K\,r - \dfrac{m\,C^2}{2} = 0\end{equation}
Etat de diffusion quand il y a répulsion
Etat de diffusion

Cette équation admet deux solutions, la solution positive donne l’expression de \(r_{\mathrm{min}}\) :

\begin{equation}r_{\mathrm{min}} = \dfrac{1}{2E_\mathrm{M}}\left(K+\sqrt{K^2+2\,m\,C^2E_\mathrm{M}}\right)\end{equation}

Cas d’une force attractive (\(K<0\))

Regardons une nouvelle fois la forme de la courbe \(E_{\mathrm{Peff}} = f(r)\) :

Profil d'énergie potentielle effective, cas du K<0
Allure de \(E_{\mathrm{Peff}} = f(r)\) dans le cas d’une force attractive

Cette fois, plusieurs cas sont possibles selon le signe de l’énergie mécanique du point M :

  • Si \(E_\mathrm{M} > 0\), on se retrouve dans la même configuration que lorsque \(K > 0\), C’est à dire que le seul mouvement possible s’effectue entre \(r_1\) et l’\(\infty\), on a encore à faire à un état de diffusion.
Etat de diffusion dans le cas attractif
Etat de diffusion lorsque \(K <0\)
  • Si \(E_\mathrm{M} <0\), le mouvement est borné entre \(r_{\mathrm{min}}\) et \(r_{\mathrm{max}}\), on parle alors d’un état lié. le mouvement est dans ce cas elliptique.
    \(r_{\mathrm{min}}\) et \(r_{\mathrm{max}}\) sont solutions de l’équation du second degré :
    \begin{equation}E_\mathrm{M}\,r^2 - K\,r - \dfrac{m\,C^2}{2} = 0\end{equation}
    Il existe un cas particulier, où l’équation précédente admet une solution double (discriminant du polynôme nul) ; dans ce cas l’état est toujours lié, avec un \(r = r_0 = \mathrm{cste}\) : le mouvement est circulaire autour de O.
Etat lié
État lié lorsque \(K <0\)

Équation polaire de la trajectoire

En repartant de l’expression de l’énergie mécanique, il est possible, à l’aide d’un changement de variable et de quelques astuces, de trouver l’équation \(r=f(\theta)\) des trajectoires évoquées ci-dessus.
On utilise notamment le changement de variable \(u = \frac{1}{r}\) qui nous permet d’obtenir une équation différentielle en u que nous savons résoudre.

Nous nous contenterons de donner ici les équations polaires finalement obtenues :

Force attractive \(K<0\)

\begin{equation}\text{Si } K<0 \text{, } r = \dfrac{p}{1+e\cos \theta} \text{ avec } p=\left|\dfrac{m\,C^2}{K}\right| \text{ et } e = \left|\dfrac{AmC^2}{K}\right| (A = \mathrm{cste})\end{equation}

Cette équation polaire est celle d’une conique, la valeur de l’excentricité e donne la forme de la conique :

Force répulsive \(K>0\)

\begin{equation}\text{Si } K>0 \text{, } r = \dfrac{p}{e\cos \theta - 1} \text{ avec } p \text{ et } e \text{ positifs}\end{equation}

On obtient donc une hyperbole (état de diffusion, voir figure 4).

Energie mécanique et trajectoires

En manipulant l’énergie mécanique, il est possible de l’exprimer en fonction des paramètres déjà utilisés pour décrire l’équation polaire de la trajectoire. On obtient :

\begin{equation}E_\mathrm{M}=-\dfrac{|K|}{2\,p}(1-e^2)= \dfrac{K^2}{2\,m\,C^2}(e^2 - 1) \end{equation}

Etant donné que les valeur de $E_M$ dépendent directement des valeurs de $e$, la valeur de $E_M$ nous indique directement la nature de la trajectoire :

Études de trajectoires particulières

Trajectoire parabolique et vitesse de libération

On parle de vitesse de libération lorsqu’un corps, soumis à l’attraction gravitationnelle d’un autre corps et distant de \(r\), a une vitesse suffisante pour s’"échapper" de cette attraction.
On se place donc dans le cas d’une force Newtonienne attractive (\(K < 0\)).

Cette vitesse de libération est atteinte si le corps est dans un état de diffusion, donc sur une trajectoire parabolique.
Ainsi, \(e = 1\) et d’après l’expression de l’énergie mécanique vue précédemment, \(E_\mathrm{M} = 0\).
Ce qui donne :

\begin{equation}\frac{1}{2}m\,v_l^2 + \dfrac{K}{r} = 0 \Longleftrightarrow v_l=\sqrt{\dfrac{-2\,K}{m\,r}}\end{equation}

Plaçons nous dans le cas de la Terre et d’un satellite de masse m, sa vitesse de libération depuis la surface terrestre est :

\begin{equation}v_l=\sqrt{\dfrac{2\,G\,m_T\,m}{m\,R_T}} = \sqrt{\dfrac{2\,G\,m_T}{R_T}} = 11\,\mathrm{km.s^{-1}}\end{equation}

Trajectoire elliptique et lois de Kepler

Cette trajectoire est importante : les lois de Kepler et en particulier la première explique que chaque planète du système solaire décrit une orbite elliptique autour du soleil qui constitue un des foyers de l’ellipse.
D’après ce qui a été vu, on a \(r=\dfrac{p}{1+e\,\cos \theta}\). Ainsi :

Expression de la vitesse sur la trajectoire

  • On peut tout d’abord exprimer l’énergie mécanique en fonction du demi-grand axe de l’ellipse : d’après le schéma ci-contre :

    \begin{equation}\begin{aligned} 2a &= r_p + r_a \\ \Longleftrightarrow 2a &= \dfrac{p}{1+e} + \dfrac{p}{1-e} = \dfrac{2p}{1-e^2}\end{aligned}\end{equation}

    D’où :

    \begin{equation}p = a(1-e^2)\end{equation}
  • Remplaçons ceci dans l’expression de l’énergie mécanique (equation ) :

    \begin{equation}E_\mathrm{M} = -\dfrac{|K|}{2a} = \dfrac{K}{2a}\end{equation}
Ellipse
Trajectoire elliptique

Ce résultat est important : dans une trajectoire elliptique, l'énergie mécanique ne dépend que du demi grand-axe de l'ellipse.

Troisième loi de Kepler

Cherchons à retrouver l’expression de la troisième loi de Kepler : \(\dfrac{T^2}{a^3}=\dfrac{4\,\pi^2}{G\,m_O}\).

Cette loi de Kepler est valable dans le cas de la trajectoire circulaire, on remplace alors \(a\) par \(r_0\), rayon de la trajectoire circulaire.

Vous trouverez à cette adresse une animation/vidéo du CEA qui résume ce que sont les lois de Kepler

Références

Annexe 1 : établissement de l’équation de la trajectoire

L’énergie mécanique s’écrit :

\begin{equation}E_\mathrm{m} = E_\mathrm{C} + E_\mathrm{P} = \dfrac{1}{2}\,m\,\dot{r}^2 + \dfrac{m\,C^2}{2\,r^2} + \dfrac{K}{r}\end{equation}

On cherche l’expression de \(r(\theta)\), donc il faut transformer le \(\dot{r}\).

On a : \(\dot{r} = \dfrac{\mathrm{d}r}{\mathrm{d}t}=\dfrac{\mathrm{d}r}{\mathrm{d}\theta}\times \dfrac{\mathrm{d}\theta}{\mathrm{d}t} = \dot{\theta}\dfrac{\mathrm{d}r}{\mathrm{d}\theta} = \dfrac{C}{r^2}\dfrac{\mathrm{d}r}{\mathrm{d}\theta}\).

Donc l’expression de l’énergie mécanique devient :

\begin{equation}E_\mathrm{m}= \dfrac{m\,C^2}{2\,r^4}\left(\dfrac{\mathrm{d}r}{\mathrm{d}\theta}\right)^2 + \dfrac{m\,C^2}{2\,r^2} + \dfrac{K}{r}\end{equation}

Utilisons un changement de variable : posons \(u=\dfrac{1}{r}\) (\(r = \dfrac{1}{u}\)), on a alors \(\dfrac{\mathrm{d}r}{\mathrm{d}\theta}=-\dfrac{1}{u^2}\dfrac{\mathrm{d}u}{\mathrm{d}\theta}\). Alors :

\begin{equation}\begin{aligned} E_\mathrm{m}&= \dfrac{m\,C^2\,u^4}{2}\times \dfrac{1}{u^4} \times \left(\dfrac{\mathrm{d}u}{\mathrm{d}\theta}\right)^2 + \dfrac{m\,C^2\,u^2}{2} + K\,u \\ &= \dfrac{m\,C^2}{2} \left(\left(\dfrac{\mathrm{d}u}{\mathrm{d}\theta}\right)^2+u^2\right) + K\,u\end{aligned}\end{equation}

Utilisons à présent la constance de l’énergie mécanique :

\begin{equation}\begin{aligned} E_\mathrm{m} = \mathrm{cste} \Longleftrightarrow \dfrac{\mathrm{d}E_\mathrm{m}}{\mathrm{d}\theta} = 0 & \Longleftrightarrow m\,C^2 \left(\dfrac{\mathrm{d}u}{\mathrm{d}\theta}\dfrac{\mathrm{d}^2u}{\mathrm{d}\theta^2}+u\,\dfrac{\mathrm{d}u}{\mathrm{d}\theta}\right)+ K\,\dfrac{\mathrm{d}u}{\mathrm{d}\theta} = 0 \\ & \Longleftrightarrow \dfrac{\mathrm{d}u}{\mathrm{d}\theta} \left(m\,C^2\,\left(\dfrac{\mathrm{d}^2u}{\mathrm{d}\theta^2} + u\right)+K\right) = 0\end{aligned}\end{equation}

Il y a donc deux possibilités :

Annexe 2 : paramètres d’une ellipse

Demi grand-axe \(a\)

On a établi que \(a = \dfrac{p}{1-e^2}\)

Distance focale \(c\)

C’est la distance entre le centre de l’ellipse et un des foyers (le centre de force) :

\begin{equation}c = a - OP = a - \dfrac{p}{1+e} = \dfrac{p}{1-e^2} - \dfrac{p}{1+e} = \dfrac{p\,e}{1-e^2}\end{equation}

Excentricité \(e\)

L’excentricité de l’ellipse peut alors s’exprimer de la façon suivante : \(e = \dfrac{c}{a}\).

Demi-petit axe \(b\)

On cherche ici à relier \(b\) à \(a\) et \(p\). On connaît une relation générale des coniques qui donne : \(a^2 = b^2 + c^2\). On peut alors écrire :

\begin{equation}b^2 = a^2 - c^2 = \dfrac{p^2}{(1-e^2)^2} - \dfrac{p^2e^2}{(1-e^2)^2} = \dfrac{p^2}{1-e^2}\end{equation}

D’où :

\begin{equation}b = \dfrac{p}{\sqrt{1-e^2}} = a\,\sqrt{1-e^2} = \sqrt{a\,p} \quad \text{ou} \quad b^2 = a\,p\end{equation}

Derniers ajouts

Physique à l'ENSCR

physique àl'ENSCR

Retrouver, entre autres, des contenus de travaux pratiques, produits par l'équipe de physique de l'ENSCR


écrire :

\begin{equation}b^2 = a^2 - c^2 = \dfrac{p^2}{(1-e^2)^2} - \dfrac{p^2e^2}{(1-e^2)^2} = \dfrac{p^2}{1-e^2}\end{equation}

D’où :

\begin{equation}b = \dfrac{p}{\sqrt{1-e^2}} = a\,\sqrt{1-e^2} = \sqrt{a\,p} \quad \text{ou} \quad b^2 = a\,p\end{equation}

Documents de référence

Derniers ajouts

Liens

Physique à l'ENSCR

physique àl'ENSCR

Retrouver, entre autres, des contenus de travaux pratiques, produits par l'équipe de physique de l'ENSCR