mécanique du point matériel et problème à deux corps

(les vecteurs sont en caractères gras)

Sommaire

I - Cinématique
I - 1 - Vitesse, quantité de mouvement, accélération
I - 2 - Mouvement à accélération centrale
II - Changement de repère; formules de composition des vitesses et des accélérations; notion de
force d'inertie
III - Dynamique
III - 1- Lois de Newton
III - 2- Travail, puissance, énergie cinétique, potentielle, mécanique
IV - Dynamique des mouvements de rotation
IV - 1 - Moment d'une force et moment cinétique
IV - 2 - Théorème du moment cinétique
IV - 3 - Compléments: mécanique du solide en rotation
V - Le problème à deux corps
V - 1 - Moment cinétique et quantité de mouvement
V - 2 - Référentiel barycentrique: quantité de mouvement et moment cinétique
V - 3 - Principe fondamental de la dynamique dans le référentiel barycentrique
V - 4 - Mouvement fictif de la particule de masse réduite µ dans le repère barycentrique
V - 5 - Energie mécanique et condition pour avoir une ellipse, une parabole ou une
hyperbole
V - 6 - Particularités du mouvement elliptique et lois de Kepler

I - Cinématique

I - 1 - Vitesse, quantité de mouvement, accélération

- vitesse v d'une particule de masse m située au point M dans un repère d'origine O

v = dOM/dt (unité: m s-1)

v a pour composantes (dx/dt, dy/dt, dz/dt) si les coordonnées de M sont (x, y, z) et dépendent du
temps.

- quantité de mouvement ou impulsion p = m v (unité: kg m s-1)

- accélération a d'une particule de masse m située au point M dans un repère d'origine O

a = dv/dt = d²OM/dt² (unité: m s-2)

a a pour composantes (d²x/dt², d²y/dt², d²z/dt²) si les coordonnées de M sont (x, y, z) et dépendent
du temps.

- Dans le repère de Frênet lié à la masse m (repère t, n où a

t est le vecteur unitaire tangent à la trajectoire et n est
le vecteur unitaire normal en M), on a: n
t
v = v t (v vitesse vectorielle, v vitesse algébrique)
M
 dans ce repère, si s est l'abscisse curviligne du point M, on a ds/dt = v et dt/ds = n/R, d'où:

a = (dv/dt) t + v²/R n

R est le rayon de courbure local (varie le long de la trajectoire)

M
- cas du mouvement circulaire de rayon R fixé v
à la vitesse angulaire ω constante (unité: radian s-1) : a
θ
v=ωRt
a = ω² R n O R

le vecteur vitesse est v = ω R t est tangent au cercle

le vecteur accélération a est normal au cercle et centripète dirigé vers le centre O

- Dans le plan en coordonnées polaires (r, θ) dans le repère (er, eθ):

OM = r er
eθ
La dérivée de er par rapport à θ est eθ
de sorte que der/dt = (der/dθ) (dθ/dt) = dθ/dt eθ, er
M
La dérivée de eθ par rapport à θ est - er, r
de sorte que deθ/dt = (deθ/dθ) (dθ/dt) = - dθ/dt er θ
O
v = dOM/dt = dr/dt er + r dθ/dt eθ

a = dv/dt = d²OM/dt² = [d²r/dt² - r (dθ/dt)²] er + [2 dr/dt dθ/dt + r d²θ/dt²] eθ

a = [d²r/dt² - r (dθ/dt)²] er + (1/r) d/dt(r² dθ/dt) eθ

I - 2 - Mouvement à accélération centrale (colinéaire au rayon vecteur r er)

Un tel mouvement n'a pas d'accélération sur eθ et est donc tel que d/dt(r² dθ/dt) = 0,

soit r² dθ/dt = C = constante des aires (unité m² s-1)

L'aire dS balayée par le rayon vecteur en dt est dS = [1/2 r² dθ/dt] dt

Alors dS/dt = C / 2 = constante = vitesse aréolaire (unité m² s-1)

Dans ce cas, on peut exprimer la vitesse et l'accélération uniquement en fonction de r(t):

v = dr/dt er + C/r eθ
a = [d²r/dt² - C²/r3] er

Un tel mouvement se produit lorsque la particule est soumise à une force dirigée selon le rayon
vecteur OM = r er

C'est le cas de:

- l'attraction gravitationnelle (force en - K m m'/r² er où K = 6.67 10-11 SI),

- de la force électrostatique (force en q q' / 4πε0r² er où 1/4πε0 = 9 109 SI),
- ou encore de l'oscillateur harmonique soumis à une force de rappel de la forme - k r er (k =
constante de raideur en N m-1).

II - Changement de repère; formules de composition des vitesses et des accélérations; notion

de force d'inertie

Considérons deux repères R et R', de même origine O, R' étant en rotation à la vitesse angulaire Ω
par rapport à R, et Ω étant le vecteur rotation. La direction du vecteur Ω est celle de l'axe de
rotation. La norme ||Ω|| est la vitesse angulaire de rotation de R' par rapport à R en rd/s.

On montre que si A est un vecteur, (dA/dt)R = (dA/dt)R' + Ω Λ A

Cette relation permet de changer de repère. Ainsi, pour un point M:

(dOM/dt)R = (dOM/dt)R' + Ω Λ OM

Posons va = (dOM/dt)R vitesse absolue dans R et vr = (dOM/dt)R' vitesse relative dans R'

alors va = vr + Ω Λ OM constitue la formule de composition des vitesses.

Dérivons par rapport au temps dans le référentiel R:

(dva/dt)R = (dvr/dt)R + [d(Ω Λ OM)/dt]R

Appliquons la formule de changement de repère (dA/dt)R = (dA/dt)R' + Ω Λ A

(dva/dt)R = (dvr/dt)R' + Ω Λ vr + dΩ/dt Λ OM + Ω Λ [(dOM/dt)R' + Ω Λ OM] Ω
M
= (dvr/dt)R' + 2 Ω Λ vr + dΩ/dt Λ OM + Ω Λ (Ω Λ OM) H

Or Ω Λ OM = Ω Λ HM où H est la projection de M sur l'axe de rotation

(OM = OH + HM); comme OH est colinéaire à Ω, O
celà implique Ω Λ OH = 0.

Ω Λ (Ω Λ HM) = (Ω.HM) Ω - Ω² HM (formule du double produit vectoriel, dans laquelle

Ω.HM = 0 en raison de leur orthogonalité).

Finalement, (d²va/dt²)R = (dvr/dt)R' + 2 Ω Λ vr + dΩ/dt Λ OM - Ω² HM

Posons aa = (d²OM/dt²)R accélération absolue dans R et ar = (d²OM/dt²)R' accélération relative

dans R', et supposons que Ω soit un vecteur constant; on obtient la formule de composition des
accélérations:

aa = ar + 2 Ω Λ vr - Ω² HM

- Ω² HM est l'accélération d'entraînement centripète

2 Ω Λ vr est l'accélération complémentaire de Coriolis (elle est nulle si vr = 0)
 Dans le repère R', une particule de masse m à la vitesse v subit alors deux forces d'inertie:
m Ω² HM est la force d'inertie centrifuge
- 2 m Ω Λ v est la force d'inertie de Coriolis (nulle si la vitesse v est nulle dans R')

Par exemple, la force d'inertie centrifuge est celle que l'on ressent sur un manège ou sur une route
dans un virage: elle est dirigée vers l'extérieur. La force d'inertie de Coriolis oriente le mouvement
des masses d'air dans le sens trigonométrique autour d'une dépression, ou dans le sens horaire
autour d'un anticyclone, dans l'hémisphère Nord de la Terre (et vice versa dans l'hémisphère Sud).
Elle est également à l'origine de mouvements de type anticyclonique dans l'atmosphère solaire au
voisinage des taches, contribuant à expliquer le cycle solaire de 11 ans.

III - Dynamique

III - 1- Lois de Newton

Un référentiel dans lequel s'applique les lois de Newton est dit galiléen. Par exemple, le référentiel
de Kepler, ayant pour centre le Soleil et trois axes othogonaux pointant vers trois étoiles de la
Galaxie peut être considéré comme galiléen pour décrire le mouvement des planètes et des comètes.

- 1ère loi de Newton: dans un référentiel galiléen, un point matériel isolé (donc ne subissant aucune
force) est en translation rectiligne uniforme.

- 2ème loi de Newton ou ou Principe Fondamental de la Dynamique (PFD):

dp/dt = m dv/dt = m a = ∑ F (somme des forces appliquées au point matériel, unité N ou Newton)

Si ∑ F = 0, alors v = constante (et on retrouve la 1ère loi de Newton)

- 3ème loi de Newton ou principe de l'action et de la réaction (ou des actions mutuelles) entre deux
points matériels A et B:

FA→B + FB→A = 0
où FA→B est la force exercée par A sur B et FB→A est la force exercée par B sur A.

Remarque: Le PFD s'applique à un solide de masse m sous la forme m dvG/dt = m aG = ∑ F

où vG et aG sont respectivement la vitesse et l'accélération du centre de masse ou de gravité.

III - 2- Travail, puissance, énergie cinétique, potentielle, mécanique

- Travail d'une force

dW = F.dl constitue le travail élémentaire dans un déplacement dl; dW = 0

lors d'un déplacement du point A vers le point B B
dW > 0
B
WAB = ∫ F.dl (unité: J ou Joule) dW < 0
A

dl
c'est un produit scalaire, le travail est moteur si dW > 0 A
résistant si dW < 0
 dl est tangent à la trajectoire: c'est le déplacement élémentaire de la masse m qui subit la force

exemple: lors d'un déplacement horizontal, le poids m g ne travaille pas car il est orthogonal
au déplacement; lors d'un déplacement vertical vers le haut d'une hauteur h, il vaut - m g h (il est
résistant); pour un déplacement vers le bas, vaut m g h et est moteur.

- Puissance instantanée d'une force

P = F.v (unité: W ou Watt)

(produit scalaire de la force avec la vitesse du point matériel, peut être > 0 ou < 0)

- Energie cinétique

Ec = 1/2 m v² (unité: J ou Joule)

pour un solide de masse m, Ec = 1/2 m vG² où vG est la vitesse du centre de masse ou de gravité

- théorème de l'énergie cinétique de la position A vers B:

EcB - EcA = WAB

démonstration simple à partir du PFD: m dv/dt = F

effectuons le produit scalaire avec le vecteur vitesse v = dOM/dt :
m v.dv/dt = F.v
donc dEc/dt = F.dOM/dt
dEc = F.dOM
B B
et de A vers B: ∫ dEc = ∫ F.dOM c'est à dire Ec(B) - Ec(A) = WAB
A A

- Force dérivant d'une énergie potentielle, dite conservative

Une force est conservative si elle dérive d'un gradient, soit si F = - grad Ep

Ep est la fonction énergie potentielle (unité: J ou Joule); elle dépend des variables de position.

La force est conservative parce que sa circulation sur un chemin fermé est nulle.

En effet, ∫ F.dl = - ∫ grad Ep . dl = - ∫ dEp = 0 car on revient au point de départ.

exemples de forces et d'énergie potentielle associée telle que F = - grad Ep

* poids m g, Ep(z) = m g z (z = hauteur de la masse m)

* force de gravitation entre deux masses m et m' distantes de r

loi de Newton F = - K m m'/r² u, Ep(r) = - Km m' / r

* force électrique entre deux charges q et q' distantes de r

loi de Coulomb F = q q'/4πε0r² u, Ep(r) = (1/4πε0) q q' / r

* force de rappel d'un ressort d'allongement x et de raideur k:

 F = - k x, Ep(x) = 1/2 k x²
et plus généralement, si F = - k OM, Ep(x) = 1/2 k ||OM||²

* charge q et de masse m dans un champ électrique dérivant du potentiel V:

E = - grad V, Ep = q V

- Principe de conservation de l'énergie mécanique

lorsque la force dérive d'une énergie potentielle:

on a vu ci dessus dans la démonstration du théorème de l'énergie cinétique que:
dEc = F.dOM
Pour une force dérivant d'une énergie potentielle, F.dOM = - grad(Ep).dOM = - dEp
d'où dEc + dEp = 0

donc Ec + Ep = constante (unité: J ou Joule)

Il n'y a pas conservation de l'énergie mécanique en cas de dissipation par frottement, parce que les
forces de frottement ne dérivent pas d'un potentiel et ne peuvent donc être conservatives.

quelques exemples:

masse m dans le champ de pesanteur:

1/2 m v² + m g z = constante

charge électrique en mouvement dans un champ électrique dérivant du potentiel V:

1/2 m v² + q V = constante

planète de masse m gravitant autour du Soleil de masse M à la distance r:

1/2 m v² - K m M/r = constante

électron de masse m et de charge -e, autour du noyau de charge Ze (Z nombre de protons du

noyau), à la distance r, et décrit par la mécanique classique:
1/2 m v² - Ze² / 4πε0r = constante

ressort de raideur k:
1/2 m v² + 1/2 k x² = constante

- position d'équilibre stable et instable; notion de particule liée/libre

Lorsqu'une particule de masse m possède une énergie potentielle Ep(x,y,z), les positions d'équilibre
correspondent aux extrema de la fonction Ep(x,y,z). Les positions des maxima sont des points
d'équilibre instable, tandis que les minima sont des points d'équilibre stable.

Soit une particule en mouvement sur l'axe Ox de vitesse v = dx/dt et une fonction énergie
potentielle Ep(x) du type de celle présentée sur la figure ci dessous. La position xeq est une position
d'équilibre stable: c'est un minimum d'énergie potentielle.

Mathématiquement, les positions d'équilibre sont données par dEp(x)/dx = 0

Les positions stables sont telles qu'au point d'équilibre, d²Ep(x)/dx² > 0 (concavité vers le haut)

Les positions instables sont telles qu'au point d'équilibre, d²Ep(x)/dx² > 0 (concavité vers le bas)
Si la particule possède l'énergie totale E constante, deux cas sont possibles:
- à gauche, la particule est liée car son espace possible est situé dans l'intervalle [x1, x2]
- à droite, la particule est libre car son espace possible est situé dans l'intervalle [x1, +∞[

Ep(x) 1/2 mv² + Ep = E

particule
x1 xeq x2 x1 xeq
x x
zone
interdite zone interdite zone libre
interdite
1/2 mv² + Ep = E
Ep(x)
particule
liée

Dans un mouvement à accélération centrale en coordonnées polaires (r, θ), on a vu que le vecteur
vitesse est v = dr/dt er + C/r eθ où C est la constante des aires.

L'énergie cinétique s'exprime par Ec = 1/2 m [ (dr/dt)² + C²/r² ]

à laquelle on doit ajouter l'énergie potentielle Ep(r) des forces conservatives (on suppose aucune
perte par frottement) pour avoir l'énergie totale E; le principe de conservation de l'énergie s'écrit:

Ec = 1/2 m (dr/dt)² + [ 1/2 m C²/r² + Ep(r) ] = constante

La stabilité du mouvement nécessite dans ce cas l'analyse de la fonction 1/2 m C²/r² + Ep(r) en
fonction de r.

IV - Dynamique des mouvements de rotation

Imaginons qu'une particule de masse m en M effectue un mouvement de rotation autour du point O

sous l'action d'une force F.

IV - 1 - Moment d'une force et moment cinétique

- Moment de la force au point O: MO = OM Λ F F

(unité N m)
O
M
La force F est appliquée au point M.
En norme, le moment est maximal lorsque
OM et F sont orthogonaux (il est nul s'ils sont colinéaires).

- Notion de couple de forces

F
Il s'agit de la somme des moments en O A
de deux forces F et - F égales B
en norme mais opposées en direction.
-F
 MO = OA Λ (- F) + OB Λ F = (AO + OB) Λ F = AB Λ F
est indépendant de O (unité N m); on écrira alors que le couple de forces est C = AB Λ F

exemple: action d'un tournevis sur la tête d'une vis.

- Moment cinétique de la particule de masse m située en M, en un point O: KO = OM Λ m v

C'est le moment en O de la quantité de mouvement p = m v (unité: m² kg s-1)

Le moment cinétique est une notion utile pour décrire les mouvements de rotation.

exemple1: dans un mouvement circulaire uniforme de rayon R à la vitesse angulaire ω, on a en

norme: v = ω R, p = m v = m ω R, KO = R p = m ω R²

exemple2: le moment cinétique de l'électron de l'atome d'Hydrogène est quantifié par la relation
de Bohr: KO = n (h/2π) où n est un nombre entier positif. La constante de Planck réduite (h/2π)
apparaît donc comme un quantum de moment cinétique.

- Moment cinétique en coordonnées polaires (r, θ) ; cas des mouvements à accélération centrale

On a vu que v = dr/dt er + C/r eθ

Alors KO = OM Λ m v = r er Λ m (dr/dt er + C/r eθ) = m C ez

Dans un mouvement à accélération centrale, la valeur du moment cinétique est lié

à la constante des aires C par la loi KO = m C = constante.

IV - 2 - Théorème moment cinétique

- théorème du moment cinétique: dKO/dt = OM Λ F

c'est l'analogue du PFD, pour les mouvements de rotation.

démonstration simple:
KO = OM Λ p
dérivons par rapport au temps:
dKO/dt = dOM/dt Λ p + OM Λ dp/dt
or dOM/dt = v et p = m v sont colinéaires (donc leur produit vectiriel nul) et d'après le PFD, on a:
dp/dt = F d'où
dKO/dt = OM Λ F

exemple: mouvement du pendule de masse m, de longueur l dans le champ de pesanteur g, en

négligeant le poids de la tige ou du fil.
Soit ez le vecteur unitaire orthogonal au plan de la figure. O
On a un mouvement de rotation circulaire à vitesse angulaire
dθ/dt non constante. l
θ
KO = m l² dθ/dt ez M, masse m

→ dKO/dt= m l² d²θ/dt² ez mg
le moment du poids par rapport à O est: MO = - m g l sinθ ez

→ m l² d²θ/dt² = - m g l sinθ

et pour les petits mouvements (θ << 1 implique sinθ ≈ θ), d²θ/dt² + (g/l) θ = 0

le mouvement est périodique de pulsation ω = (g/l)1/2 et de période T = 2π/ω = 2π (l/g)1/2

IV - 3 - Compléments: mécanique du solide en rotation

- complément 1: moment cinétique d'un solide en rotation autour d'un axe ∆

K∆ = ∫∫∫ HM Λ v(M) dm
Ω (axe ∆)
Le moment cinétique du solide par rapport à un axe
de rotation ∆ est la somme des moments cinétiques
élémentaires, en affectant à chaque point M du solide
H M
une masse élémentaire dm = ρ dv occupant le
volume élémentaire dv (ρ masse volumique locale).
La vitesse du point M est v(M) = Ω Λ HM
où H est la projection de M sur l'axe de rotation ∆ O
et Ω le vecteur rotation autour de ∆.

K∆ = ∫∫∫ HM Λ (Ω Λ HM) dm

= ∫∫∫ [HM² Ω - (HM.Ω) HM)] dm

Or le produit scalaire HM.Ω est nul; il reste alors: K∆ = ∫∫∫ HM² Ω dm = Ω ∫∫∫ HM² dm, soit

K∆ = J Ω où J = ∫∫∫ HM² dm désigne le moment d'inertie du solide par rapport à l'axe ∆

Le théorème du moment cinétique s'écrit dans ce cas dK∆/dt = ∫∫∫ HM Λ dF(M)

où dF(M) est la force élémentaire appliquée à la masse dm = ρ dv au point M.

- complément 2: quelque exemples de moments d'inertie pour des solides homogènes de masse m
en rotation autour d'un axe ∆

sphère pleine de rayon R: J = 2/5 m R²

disque ou cylindre plein de rayon R: J = 1/2 m R²
anneau fin de rayon R: J = m R²
masse ponctuelle à la distance l: J = m l²
tige pleine de longueur l par rapport à une extrémité: J = 1/3 m l²
tige pleine de longueur l par rapport à son milieu: J = 1/12 m l²

- complément 3: énergie cinétique d'un solide en rotation autour d'un axe ∆

Ec = ∫∫∫ 1/2 v² dm où v(M) = Ω Λ HM

= ∫∫∫ 1/2 (Ω Λ HM).(Ω Λ HM) dm

On a un produit mixte dans l'intégrale que l'on peut transformer:

(Ω Λ HM).(Ω Λ HM) = Ω. [HM Λ (Ω Λ HM)]

= Ω. [ HM² Ω - (Ω.HM) HM]
= HM² Ω² car Ω et HM sont orthogonaux

donc Ec = ∫∫∫ 1/2 HM² Ω² dm = 1/2 Ω² ∫∫∫ 1/2 HM² dm,

Ec = 1/2 J Ω² où J est le moment d'inertie du solide par rapport à l'axe ∆

Si le solide de masse m et de moment d'inertie J possède en plus un mouvement de translation,

Ec = 1/2 m vg² + 1/2 J Ω² où vg est la vitesse du centre de gravité du solide

Rappelons que le centre de gravité G d'un solide de masse m est défini par la relation vectorielle:

∫∫∫ GM dm = 0 ou encore OG = (1/m) ∫∫∫ OM dm si O est l'origine d'un repère lié au solide.

Le point M décrit l'ensemble des points du solide; on lui affecte une masse élémentaire dm.

La vitesse du centre de gravité est définie par vG = (1/m) ∫∫∫ v(M) dm

Un corps comme la Terre possède deux mouvements rotatoires: la révolution autour du Soleil
d'énergie cinétique Ec1 = 1/2 m vg² = 1/2 m l² ω² où l est la distance Soleil Terre et ω la vitesse
angulaire de révolution (T = 2π/ω = 365.25 jours); la rotation diurne autour de l'axe de pôles avec
Ec2 = 1/2 J Ω² où Ω est la vitesse angulaire de rotation (T = 2π/Ω = 23 h 56 min).

Application numérique: la Terre se déplace dans l'espace à la vitesse vg = 30 km/s et tourne sur elle
même à la vitesse angulaire Ω = 3 10-4 rd/s; sa masse est m = 6 1024 kg; pour une sphère homogène
de masse m et de rayon de rayon R, J = 2/5 m R²; avec R = 6400 km, J = 1038 kg m².
Ec1 = 1/2 m vg² = 2.51033 J
Ec2 = 1/2 J Ω² = 4.5 1030 J
L'énergie cinétique de révolution autour du Soleil est donc 1000 fois plus grande que celle due à la
rotation autour de l'axe des pôles.

V - Le problème à deux corps

V - 1 - Moment cinétique et quantité de mouvement er

Comment décrire le mouvement de deux corps
isolés et en interaction gravitationnelle mutuelle ?
Tel est l'objet de cette section. F = F2
M2 , m 2
On considère deux masses ponctuelles m1 et m2
de positions M1 et M2 et de vitesses vectorielles
respectives V1 et V2 dans un référentiel galiléen, F1
dans lequel on pourra appliquer les lois de Newton.
Il peut par exemple s'agir du référentiel de Kepler
pour décrire le mouvement d'une planète ou d'une
comète autour du Soleil (en faisant abstraction M1 , m 1
des autres corps du système solaire, le problème
étant à deux corps).
L'attraction gravitationnelle entre les masses m1 et m2 est donnée par la loi de Newton:

F2 = - K m1 m2 / r² er et F1 = - F2 ; on posera F2 = - F1 = F

- Moment cinétique total

L'application du théorème du moment cinétique au système à deux corps donne, avec

KO = OM1 Λ m1 V1 + OM2 Λ m2 V2 moment cinétique total où O est l'origine du référentiel:

dKO/dt = OM1 Λ F1 + OM2 Λ F2 = (OM2 - OM1) Λ F = M1M2 Λ F = 0

car M1M2 et F sont colinéaires.

Le moment cinétique total KO est donc une constante vectorielle, et le mouvement est plan (KO est
un vecteur constant orthogonal à ce plan).

- Quantité de mouvement totale

L'application du principe fondamental de la dynamique au système à deux corps donne, avec

P = m1 V1 + m2 V2 = P1 + P2 quantité de mouvement totale:

dP/dt = F1 + F2 = 0

La quantité de mouvement totale P = P1 + P2 est donc une constante vectorielle.

V - 2 - Référentiel barycentrique: quantité de mouvement et moment cinétique

Introduisons maintenant le barycentre G du système des deux corps. G est défini par:

OG = (m1 OM1 + m2 OM2) / (m1 + m2)

La vitesse du barycentre (ou centre de gravité, ou centre de masse) est donnée par:

VG = (m1 V1 + m2 V2) / (m1 + m2) = P / (m1 + m2) = constante vectorielle

On en conclut que le référentiel barycentrique (ayant pour origine G) est en translation rectiligne
uniforme par rapport au référentiel supposé galiléen de centre O de départ; le référentiel
barycentrique est donc galiléen, et on peut étudier le mouvement des deux corps dans ce repère.

Vitesses et quantités de mouvement dans le référentiel barycentrique

posons v1 = V1 - VG et v2 = V2 - VG

les vitesses des corps de masse m1 et m2 dans le référentiel barycentrique.

v1 = V1 - (m1 V1 + m2 V2) / (m1 + m2) = [m2 / (m1 + m2)] (V1 - V2)

v2 = V2 - (m1 V1 + m2 V2) / (m1 + m2) = [m1 / (m1 + m2)] (V2 - V1)

Posons v = V2 - V1 vitesse relative de M2 par rapport à M1.

v1 = - [m2 / (m1 + m2)] v et v2 = [m1 / (m1 + m2)] v

Dans le repère barycentrique, les quantités de mouvement s'écrivent:

p1 = m1 v1 = - [m1 m2 / (m1 + m2)] v et p2 = m2 v2 = [m1 m2 / (m1 + m2)] v

Appelons la quantité µ= m1 m2 / (m1 + m2) masse réduite du système à deux corps; alors:

p1 = m1 v1 = - µ v et p2 = m2 v2 = µ v avec p1 + p2 = 0

La quantité de mouvement totale est nulle dans le repère barycentrique.

- Moment cinétique dans le référentiel barycentrique

KG = GM1 Λ m1 v1 + GM2 Λ m2 v2 = GM1 Λ (- µ v) + GM2 Λ µ v = M1M2 Λ µ v

Soit M le point défini par GM = M1M2

alors KG = GM Λ µ v

est le moment cinétique d'une particule fictive de masse µ = m1 m2 / (m1 + m2) située au point
fictif M tel que GM = M1M2 = r er et soumise à la force gravitationnelle F = - K m1 m2 / r² er

Ce moment cinétique est constant, car d'après le théorème du moment cinétique,

dKG/dt = GM1 Λ (- F) + GM2 Λ F = M1M2 Λ F = 0 car M1M2 et F sont colinéaires.

La vitesse du point M fictif est v = V2 - V1 puisque

dGM/dt = dM1M2/dt = d(OM2 - OM1)/dt = dOM2/dt - dOM1/dt = V2 - V1

V - 3 - Principe fondamental de la dynamique dans le référentiel barycentrique

Appliquons le PFD à la masse m2 dans le référentiel barycentrique

m2 dv2/dt = F

or m2 v2 = µ v d'où l'on déduit µ dv/dt = F = - K m1 m2 / r² er

Ainsi, le PFD permet de déterminer le mouvement de la particule fictive M de masse

µ = m1 m2 / (m1 + m2) dans le référentiel barycentrique de centre G, avec r = GM.

Une fois connu ce mouvement fictif, il sera aisé d'en déduire le mouvement réel des masses m1 et
m2 aux points respectifs M1 et M2 par la relation barycentrique m1 GM1 + m2 GM2 = 0.

En effet, m1 GM1 + m2 (GM1 + M1M2) = 0 et on a posé GM = M1M2. Il vient les relations:

GM1 = - m2 GM / (m1 + m2) et GM2 = m1 GM / (m1 + m2)

v1 = - m2 v / (m1 + m2) et v2 = m1 v / (m1 + m2)
On remarque que GM1 = - (m2/m1) GM2 et v1 = - (m2/m1) v2

Ainsi, les positions et les vitesses sont dans le rapport des masses.

Que se passe t-il si m1 >> m2 ? Dans ce cas, µ ≈ m2, GM1 ≈ 0 et GM2 ≈ GM

C'est presque le cas dans le système solaire si M1 représente le Soleil de masse m1.

V - 4 - Mouvement fictif de la particule M de masse réduite µ dans le repère barycentrique

Posons GM = r er et utilisons les coordonnées polaires (r, θ) dans le plan du mouvement.

Le mouvement étant à accélération centrale, on a vu que:

v = dr/dt er + C/r eθ , a = [d²r/dt² - C²/r3] er où C est la constante des aires.

Le PFD s'écrit alors µ [d²r/dt² - C²/r3] = - K m1 m2 / r²

soit encore d²r/dt² - C²/r3 = - K (m1 + m2) / r²

Pour résoudre cette équation, on fait le changement de variable suivant: u = 1/r

et sachant que:
dθ/dt = C u² (loi des aires), on transforme l'équation en r(t) contre une équation en u(θ).

dr/dt = d(1/u)/dθ dθ/dt = - (1/u²) du/dθ C u² = - C du/dθ

d²r/dt² = d(- C du/dθ)/dt = - C d²u/dθ2 dθ/dt = - C² u² d²u/dθ2

Ainsi, a = - C² u² (d²u/dθ² + u) se nomme formule de Binet pour l'accélération

et le PDF s'écrit: - C² u² (d²u/dθ² + u) = - K (m1 + m2) u²

et finalement, d²u/dθ² + u = K (m1 + m2) / C²

dont la solution est de la forme u = K (m1 + m2) / C² + A cosθ = 1/r (avec A constante), d'où

r(θ) = [C²/ K (m1 + m2)] / (1+ e cosθ) = p / (1+ e cosθ) où p = C²/ [K (m1 + m2)] (en km)

qui constitue l'équation d'une conique dont la nature dépend de la valeur de l'excentricité e.

e > 1: hyperbole d'asymptotes cosθ = -1/e et de foyer G (la particule part à l'infini)
e = 1: parabole de foyer G (la particule part à l'infini)
e < 1: ellipse de foyer G
e = 0: cercle de centre G

Le périgée de la conique est obtenu pour θ = 0 et vaut rp = p / (1 + e)

L'apogée n'existe que pour une ellipse, dans ce cas, θ = π et vaut ra = p / (1 - e)

V - 5 - Energie mécanique et condition pour avoir une ellipse, une parabole ou une hyperbole

L'énergie mécanique du mouvement est conservée et donnée par E = Ec + Ep, soit

E = 1/2 µ v² - K m1 m2 / r

où v² = (dr/dt)² + (r dθ/dt)² = [(dr/dθ)² + r²] (dθ/dt)² = [(dr/dθ)² + r²] C²/r4

Or dr/dθ = p e sinθ / (1+ e cosθ)² = r² e sinθ / p

donc v² = C² e² sin²θ / p² + C² (1 + e cosθ)² / p² = (C²/p²) (e² + 1 + 2 e cosθ)

E = 1/2 µ (C²/p²) (e² + 1 + 2 e cosθ) - (1 + e cosθ) K m1 m2 / p

Remarquant que C² = K p (m1 + m2) , on a µ C² = K p m1 m2

E = 1/2 K m1 m2 (e² + 1 + 2 e cosθ) / p - (1 + e cosθ) K m1 m2 / p

E = 1/2 K m1 m2 (e² - 1) / p = constante, qui amène à la conclusion suivante:

e > 1 → E > 0; trajectoire hyperbolique

e = 1 → E = 0; trajectoire parabolique
e < 1 → E < 0; trajectoire elliptique

- Energie en fonction de la position et de la vitesse initiale

C = r0 v0 v0
On suppose qu'en t = 0, r = r0 et v = v0, le vecteur vitesse v0 v0
étant orthogonal au rayon vecteur GM0. r0
G M0
Alors E = 1/2 µ v0² - K m1 m2 / r0 = 1/2 K m1 m2 (e² - 1) / p

Introduisons la vitesse de libération vL obtenue lorsque e = 1, E = 0 (trajectoire parabolique):

1/2 µ vL² = K m1 m2 / r0 soit vL = [2 K (m1 + m2) / r0]1/2

E = 1/2 µ (v0² - vL²) = 1/4 µ vL² (e² - 1) r0/p

d'où v0 = vL [ 1 + (e² - 1) r0/2p ]1/2

Cas de l'hyperbole ou de la parabole (e ≥ 1)

r0 = p / (1 + e) est le périgée de la trajectoire, alors v0 = vL [ (1 + e)/2 ]1/2 ≥ vL

Cas de l'ellipse (0 < e < 1)

si r0 = p / (1 + e) est le périgée de la trajectoire, alors v0 = vL [ (1 + e)/2 ]1/2 et vL/√2 < v0 < vL
si r0 = p / (1 - e) est l'apogée de la trajectoire, alors v0 = vL [ (1 - e)/2 ]1/2 et 0 < v0 < vL/√2

Cas du cercle (e = 0)
r0 = p et v0 = vL/√2

En résumé, la vitesse initiale conditionne la nature de la trajectoire:

0 < v0 < vL/√2 ellipse dont la position initiale est l'apogée

v0 = vL/√2 cercle vL = [2 K (m1 + m2) / r0]1/2
vL/√2 < v0 < vL ellipse dont la position initiale est le périgée vitesse de libération
v0 = vL parabole
v0 > vL hyperbole
 Hyperbole

équation cartésienne:
x²/a² - y²/b² = 1
asymptotes y = ± (b/a) x
M
r foyer OG = c = (a² + b²)1/2
périgée OP = a
excentricité e = c/a > 1
G r0 P O b = a (e² - 1)1/2

équation polaire:
r = p/(1 + e cosθ) = GM
p = a (e² - 1)
périgée rp = GP = p/(1 + e)

Parabole

M équation cartésienne:
r x = - y² / 2p
foyer OG = c = p/2
excentricité e = 1
G r0 O
équation polaire:
r = p/(1 + cosθ) = GM
périgée rp = GO = p/2

Ellipse
B
équation cartésienne:
M x²/a² + y²/b² = 1
OA = a et OB = b
r
foyer OG = c = (a² - b²)1/2
GB = a
C O G r0 A
excentricité e = c/a < 1
b = a (1 - e²)1/2

équation polaire:
r = p/(1 + e cosθ) = GM
p = a (1 - e²)
périgée ra = GA = a (1 - e)
apogée rp = GC = a (1 + e)
Dans chaque cas, p = C² / K (m1 + m2) et C = r0 v0
Trajectoires possibles d'un objet de masse réduite µ = m1 m2 / (m1 + m2) autour du centre de gravité
ou foyer noté F. La vitesse initiale v0 au point de départ D est supposée orthoradiale (donc
orthogonale à FD). La nature de la trajectoire dépend du paramètre x = v0/vL où vL est la vitesse de
libération. L'excentricité e dépend aussi du rapport x = v0/vL. Si m1 >> m2, la masse m1 est en F et le
corps de masse m2 gravite autour. La vitesse de libération est dans ce cas vL = [2 K m1 / r0]1/2.
Trajectoires partant à l'infini:
x > 1 : hyperbole (et ses asymptotes) en violet
x = 1 : parabole en vert
Trajectoires fermées:
1/√2 < x < 1: ellipse dont D est le périgée en noir
x = 1/√2 : cercle en rouge (e = 0, r = constante)
x < 1/√2 : ellipse dont D est l'apogée en noir
Cinématique des positions d'un objet de masse réduite µ = m1 m2 / (m1 + m2) autour du centre de
gravité ou foyer F. La vitesse initiale v0 au point de départ D est supposée orthoradiale (donc
orthogonale à FD). Entre deux positions successives (représentées par les points), l'intervalle de
temps est constant. Plus les points sont proches, et plus la vitesse est faible. Avec x = v0/vL,
Trajectoires partant à l'infini:
x > 1 : hyperboles en violet, vitesse décroissante vers une valeur asymptote non nulle
x = 1 : parabole en vert, vitesse décroissante vers zéro
Trajectoires fermées:
1/√2 < x < 1: ellipse dont D est le périgée en noir (vitesse maximale au périgée D)
x = 1/√2 : cercle en rouge (e = 0, vitesse constante)
x < 1/√2 : ellipse dont D est l'apogée en noir (vitesse minimale à l'apogée D)

Vitesse v/v0 en fonction de la distance r/r0 au foyer F. D point de départ; v0 vitesse initiale.
Trajectoires partant à l'infini, r(t) > r0 , à vitesse décroissante v(t) < v0 :
x > 1 : hyperbole (violet, vitesse asymptote) et x = 1 : parabole (vert, vitesse nulle à l'infini)
Trajectoires fermées à vitesse v(t) < v0 : 1/√2 < x < 1, ellipse dont D est le périgée en noir, r(t) > r0
Trajectoires fermées à vitesse v(t) = v0 : x = 1/√2, cercle = point rouge, r(t) = r0
Trajectoires fermées à vitesse v(t) > v0 : x < 1/√2, ellipse dont D est l'apogée en noir, r(t) < r0
Cinématique du problème à deux corps de masse m1 et m2 (la masse du corps vert est double de
celle du rouge) gravitant autour de leur centre de gravité situé au foyer F. Les droites en trait
pointillé relient les positions respectives des deux corps.
Entre deux positions successives (représentées par les points), l'intervalle de temps est constant.
L'excentricité des ellipses dépend du paramètre x = v0/vL (v0 vitesse initiale du corps fictif de masse
réduite µ = m1 m2 /(m1 + m2), vL = [2K(m1 + m2)/r0]1/2 , r0 étant la distance initiale au centre de
gravité F). Cas où x = 0.95. Les vitesses et distances des corps sont dans le rapport des masses.

Même figure, mais avec

x = 0.75 (proche de la
valeur critique 0.707
pour laquelle on aurait
eu deux cercles
concentriques de centre
F (ou centre de gravité
du système). La masse
du corps vert est double
de celle du rouge. On a
les relations:

m1 FM1 + m2 FM2 = 0
m 1 v1 + m 2 v2 = 0
V - 6 - Particularités du mouvement elliptique

L'énergie E prend la forme suivante, sachant que p = a (1 - e²)

E = - 1/2 K m1 m2 / a

où a est le demi grand axe de l'ellipse.

Si e est son excentricité (0 < e < 1), le demi petit axe est b = a (1 - e²)1/2

e est déterminée par la condition initiale:

vL/√2 < v0 < vL , e = 2 (v0 /vL)² - 1

0 < v0 < vL/√2 , e = 1 - 2 (v0 /vL)²

r = p/(1 + e cosθ) fournit l'apogée ra = a (1 + e) et le périgée rp = a (1 - e)

On a ra + rp = 2 a

Si va et vp désignent la vitesse à l'apogée et au périgée, la loi des aires s'écrit:

C = ra va = rp vp = r0 v0

où r0 et v0 sont respectivement la position et la vitesse initiale (supposée orthoradiale).

Par ailleurs, la conservation de l'énergie entre l'apogée et le périgée donne une seconde relation:

E = 1/2 µ va² - K m1 m2 / ra = 1/2 µ vp² - K m1 m2 / rp où µ = m1 m2 / (m1 + m2) masse réduite

comme vp = (ra / rp) va

1/2 µ va² (1 - (ra/rp)²) = - K m1 m2 (1/rp -1/ra)

On a vu que ra = a (1 + e) et le périgée rp = a (1 - e), il vient alors:

1/2 µ va² [1 - (1 + e)²/(1 - e)²] = - K m1 m2 (2 a e) / a² (1 - e²)

1/2 va² (-4 e) / (1 - e)² = - K (m1 + m2) (2 e) / a (1 - e²)

ce qui donne finalement va = [K (m1 + m2) / a]1/2 [(1 - e) / (1 + e)]1/2

puis vp = [K (m1 + m2) / a]1/2 [(1 + e) / (1 - e)]1/2

- Période de révolution T

On peut la calculer à lartir de la loi des aires: dS/dt = C/2 = constante

ce qui donne en intégrant: S = surface de l'ellipse = π a b = C T/2 = ra va T/2

d'où T = 2π a b / ra va

Comme ra = a (1 + e) et va = [K (m1 + m2) / a]1/2 [(1 - e) / (1 + e)]1/2

En élevant T au carré, on obtient: T² = 4π² a² b² / K (m1 + m2) a (1 - e²)

Cependant, b² = a² (1 - e²), d'où T² = 4π² a3 / K (m1 + m2), loi qui s'énonce ainsi:

T²/a3 = 4π² / K (m1 + m2)

- Lois de Kepler

* 1ère loi: la trajectoire d'une planète est une ellipse dont le centre du Soleil est un foyer

* 2ème loi: les aires balayées par le rayon vecteur sont proportionnelles au temps mis à les parcourir
(loi des aires r² dθ/dt = C = constante des aires)

* 3ème loi: T²/a3 = 4π² / KMs

où T est la période de révolution, a le demi grand axe de l'ellipse, K la constante de gravitation
(6.67 10-11 SI) et Ms la masse du Soleil (2 1030 kg)

- Données sur le Soleil et les planètes du système solaire

Soleil: masse 2 1030 kg, rayon 700 000 km, densité moyenne 1.4.

Excentricités actuelles (elles ne sont pas constantes dans le temps !), masse en kg et demi grand axe
des planètes en UA (1 UA = 149.6 109 m)
e Masse (kg) demi grand axe a (UA)
23
Mercure 0.2056 la plus forte 3.3 10 0.39
Vénus 0.0067 la plus faible 4.9 1024 0.72
24
Terre 0.0167 6.0 10 1.00
Mars 0.0934 6.3 1023 1.52
Jupiter 0.0484 1.9 1027 la plus massive 5.20
26
Saturne 0.0542 5.7 10 9.55
Uranus 0.0472 8.7 1025 19.22
26
Neptune 0.0086 1.1 10 30.11

Période de révolution T autour du Soleil (jours), rayon équatorial R des planètes (km) et densité d
T (jours) R (km) d (sans unité)
Mercure 88.0 2439 5.5
Vénus 224.7 6052 5.1
Terre 365.25 6378 5.5
Mars 687.0 3397 3.9
Jupiter 4332.6 71400 1.3
Saturne 10759 60000 0.7
Uranus 30688 25400 1.3
Neptune 60181 24300 1.6
Exercice: le Soleil à un rayon de 700 000 km. En supposant que toute la masse du système solaire
est concentrée dans Jupiter (soit 1/1050 de la masse du Soleil), à quelle distance du centre du Soleil
(en km) se trouve le centre de gravité G ?
Donnée: distance Soleil Jupiter = 5.202 UA et 1 UA = 149.6 millions de km.

Réponse:
on a la relation vectorielle de définition du centre de gravité S:
MS SG + MJ JG = 0
et la relation de Chasles
JG = JS + SG
où J est le centre de Jupiter de masse MJ, et S le centre du Soleil de masse MS
d'où SG = SJ MJ / (MS + MJ) ≈ SJ MJ / MS
et en norme, SG ≈ SJ /1050 soit 4.95 10-3 UA = 741000 km soit environ le rayon solaire. Le centre
de gravité est proche de la surface solaire. La période de révolution autour du centre de gravité est
celle de Jupiter soit 12 ans environ.

Exercice: l'excentricité de l'orbite terrestre "e" autour du Soleil varie de 0 à 0.06 en cent mille ans
environ. On fait l'hypothèse que le demi grand axe "a" est invariable et égal à sa valeur actuelle de
149.6 106 km.
1) dans quel intervalle varie l'apogée ra = a (1 + e) quand e varie de 0 à 0.06 ?
2) dans quel intervalle varie le périgée rp = a (1 - e) quand e varie de 0 à 0.06 ?
3) la puissance P(t) reçue par mètre carré au niveau de l'orbite terrestre est P(t) = L/4πr(t)² (où r(t)
est la distance Soleil Terre). Que vaut P(t) quand e = 0 donc r = a au cours d'une révolution
circulaire (on donne L = 3.86 1026 W luminosité du Soleil) ?
4) La puissance instantanée P(t) reçue par la Terre par mètre carré est (L/4πC) dθ/dt où C est la
constante des aires. Elle varie en fonction du temps (car dθ/dt varie sur l'orbite). Exprimer
l'énergie totale E reçue par mètre carré sur une révolution en fonction de L et de C.
5) Sachant que C = [K Ms a (1-e²)] 1/2, où K est la constante de gravitation et Ms la masse du Soleil,
exprimer E en fonction de L, (K Ms a) et e.
6) sachant que (1-e²)-1/2 ≈ 1 - e²/2, que vaut la variation relative de E quand e varie de 0 à 0.06 ?
Cette variation à elle seule est-elle susceptible d'influencer fortement le climat ?
7) dans quel intervalle varie la puissance reçue par mètre carré P(t) = L/4πr(t)² entre apogée et
périgée quand e = 0.0167 (valeur actuelle) ?
8) dans quel intervalle varie la puissance reçue par mètre carré P(t) = L/4πr(t)² entre apogée et
périgée quand e = 0.06 (valeue maximale) ?

Réponses:
1) ra = a (1 + e) varie de 149.6 à 158.6 106 km
2) rp = a (1 - e) varie de 149.6 à 140.6 106 km
3) la puissance reçue par la Terre L/4πr² est de 1372 W/m² et ne varie pas dans le temps (e = 0, ra =
rp = a = 149.6 106 km)
4) E = L / 2C
5) E = L / (2 [K Ms a (1-e²)]1/2 )
6) variation relative 0.06²/2 = 0.0018 soit seulement 0.18 % (faible influence climatique)
7) ra = 152.1 106 km et rp = 147.1 106 km; la puissance reçue par la Terre L/4πr² varie entre 1328
W/m² (apogée à e = 0.0167) à 1420 W/m² (périgée à e = 0.0167), cette variation de 7 % au cours de
la révolution est modérée.
8) la puissance reçue par la Terre L/4πr² varie entre 1221 W/m² (apogée à e = 0.06) à 1554 W/m²
(périgée à e = 0.06), cette variation de 25% au cours de la révolution est importante.
Nota: la luminosité L du Soleil varie avec le cycle solaire mais à raison de moins de 0.1 %, soit
nettement moins que les fluctuations d'origine orbitale.
Exercice: l'excentricité de l'orbite terrestre e autour du Soleil est de 0.0167. La distance moyenne
Soleil Terre est a = 149.6 106 km. Calculer le diamètre angulaire (donc sur le ciel pour un
observateur) en minutes de degré entre l'apogée (Juillet) et le périgée (Janvier) sachant que le
diamètre réel du Soleil est de 1.4 106 km.

Réponses:
apogée = a (1+e) = 152.1 106 km
périgée = a (1-e) = 147.1 106 km
diamètre apparent en radians = diamètre solaire / distance Soleil-Terre
1 radian = 360 degrés
1 degré = 60 minutes
donne une variation de 31.5' à 32.5' de l'apogée au périgée.

Exercice: la comète de Halley a une période T = 76 ans. Son excentricité e vaut 0.97. On demande
de calculer le demi grand axe a (en km puis en UA), l'apogée, le périgée (en km puis en UA), la
vitesse à l'apogée et la vitesse au périgée (en km/s). On donne la masse du Soleil: 2 1030 kg, et
1 UA = 149.6 106 km.

Réponses:
3ème loi de Kepler, a = 2.7 109 km = 18.0 UA
apogée = a (1+e) = 5.3 109 km = 35.4 UA
périgée = a (1-e) = 8.1 107 km = 0.55 UA
vitesse à l'apogée = 0.9 km/s
vitesse au périgée = 57 km/s

Exercice: vitesse de libération.

1) Calculer la vitesse de libération à la surface du Soleil, de la Terre et de la Lune en km/s.
Calculer ensuite la vitesse de libération du Soleil depuis l'orbite terrestre en km/s.
2) un satellite tourne autour de la Terre en orbite basse circulaire; donner en appliquant le principe
fondamental de la dynamique l'expression de sa vitesse de rotation v en fonction de K (constante de
gravitation), de M (masse de la Terre) et du rayon R de la Terre. Que vaut numériquement cette
vitesse en km/s ?
3) calculer de la même façon la vitesse de la Terre sur son orbite en la supposant circulaire (en
km/s).
Données:
Masses: Soleil, 2 1030 kg, Terre 6 1024 kg, Lune, 7.4 1022 kg
Rayons: Soleil, 700 000 km, Terre, 6400 km, Lune, 1700 km
Distance Soleil-Terre: 149 106 km

Réponses:
1) vitesse de libération aux surfaces:
du Soleil, 617 km/s
de la Terre, 11.2 km/s
de la Lune, 2.4 km/s
vitesse de libération du Soleil depuis l'orbite terrestre: 42.1 km/s
2) v = (K M / R)1/2
v = 7.9 km/s
3) v = 29.8 km/s