L'inégalité de Sundman (1912) concerne le Problème à N corps en mécanique céleste :
![]() |
elle constitue une des clefs d'entrée dans l'étude du destin du problème à 3 corps.

La moitié de sa dérivée temporelle
où est 2.E(cinétique) toujours positive ; et , l'énergie gravitationnelle des N points, toujours négative.
Bien sûr, le Lagrangien est T- V.
L'Hamiltonien T+V = H(x, y)= Ho constante ( {y} est l'ensemble des quantités de mouvement)
L'inégalité conduit à :
, le deuxième membre s'appelle la barrière centrifuge (Leibniz 1689).
On peut utiliser la notation ABC pour le triangle, et les formules usuelles des distances mutuelles ; on pose D = max(a, b, c) et d = min(a, b, c).
On déclare P(x) == (x-m1)(x-m2)(x-m3) == x^3- M.x^2+ P.x -M3^3
Il s'ensuit que :
Donc
et
.
Soit :
Un raisonnement similaire avec l'énergie potentielle donne :
Et par conséquent jouera vraisemblablement un rôle important dans l'analyse de la situation.
Pour un système stable (existent cercles périgée inf d(t) et apogée sup D(t)), Eo est négative. démonstration : I" =4T + 2 V ⇒ viriel 2
et si Eo>0, le système est ouvert : D→infty : en effet I" >0 démonstration : on utilise la fonction de Sundman (voir plus loin)
Attention : si Eo<0, cela n'implique pas que le système soit stable ! En effet, deux étoiles peuvent indéfiniment se rapprocher, et la troisième partir à l'infini doucement.
Il n'y a de collision triple stricte que si L = 0 demonstration : pour avoir I(t0) =0 , V(t0)= -infty, donc d'après le viriel I"(t) >0 près de to, donc I(t) décroit de manière monotone I'(t)<0 : alors l'inégalité de Sundman donne 2L^2 Ln[I(t1)/I(t)] < cste; donc I(t) reste borné sauf si L=0
On introduit par généralisation du problème à 2 corps les notations :
, appelé demi-grand axe, constant.
, appelé paramètre, constant.
, en gros la taille, variable.
, environ la petite taille, variable.
avec les inégalités usuelles de convexité :
(soit
)
Le viriel devient :
La fonction de Sundman S(t) est :
L'inégalité de Sundman devient :
⇒
,
donc R varie comme S(t).
L'inégalité triangulaire s'écrit : 
d'où
, soit
|
|
D'autre part on en déduit , ce qui semble naturel, mais fait intervenir la décomposition d'Eckart des vitesses :
avec ,
donc en fait, ce n'est pas si simple à démontrer (les chimistes appellent ce référentiel tournant le référentiel d'Eckart ; en astronomie c'est le SAM de Saari).
le cas intéressant est celui dit de la binaire de Hill : si p>> a, le système tourne très vite à |énergie| faible, alors le système restera toujours une binaire serrée autour de laquelle tourne la troisième étoile (qui peut ou non s'échapper).
L'habitude veut que l'on représente les iso (R/r) qui ressemblent vaguement à la surface de Hill du problème restreint( m3<< m1 et m2), dans le référentiel où AB est pris comme axe des x, et C joue le rôle de troisième étoile ; mais cette représentation est un peu fallacieuse, car elle ne met pas en exergue le rôle symétrique des trois masses.