# Théorème de Siacci

Le théorème de Siacci est un théorème de cinématique qui permet d'exprimer l'accélération d'un point matériel en la projetant sur les directions du vecteur position (direction radiale) et du vecteur vitesse (direction tangentielle) de ce point à un instant donné. Il est important de souligner que ces deux directions ne sont en général pas orthogonales, même pour une trajectoire plane. Il est dû au mathématicien et balisticien Italien Francesco Siacci (1839–1907).
La décomposition de Siacci diffère donc sensiblement de l'expression du vecteur accélération du point selon les directions tangentielle et orthoradiale, donc dans un repère local orthonormé par construction, donnée dans le cas de la décomposition de Serret-Frenet.
L'intérêt du théorème de Siacci est surtout historique, mais il peut être intéressant dans le cas du mouvement à force centrale.

## Notions préliminaires : repère de Serret-Frenet
### Définitions générales

En cinématique du point, la trajectoire est la courbe géométrique décrite par le lieu des positions successives du point matériel (noté M) au cours du temps t, dans son déplacement par rapport à un référentiel donné. Il s'agit donc d'une courbe paramétrique, qui n'est pas nécessairement plane. Il est cependant toujours possible d'orienter positivement cette courbe selon le sens du mouvement, et de prendre un point quelconque de la trajectoire comme origine.
Dans ce contexte, il est possible de définir l'abscisse curviligne s(t) comme la distance parcourue le long de la trajectoire par le point matériel depuis l'origine sur la trajectoire, à une date t donnée. La distance est bien sûr évaluée selon le sens du mouvement, et donc avec cette convention l'abscisse curviligne est une grandeur positive.
Le vecteur position du point matériel P est le vecteur $\mathbf {r} =\mathbf {OP} (t)$, O étant un point quelconque de l'espace, donc pas nécessairement situé sur la trajectoire, pris pour origine.
Les vecteurs vitesse et accélération du point matériel P sont alors définis comme les deux dérivées successives par rapport au temps du vecteur position:
    $\mathbf {v} ={\frac {d\mathbf {r} }{dt}}$;
    $\mathbf {a} ={\frac {d\mathbf {v} }{dt}}={\frac {d^{2}\mathbf {r} }{dt^{2}}}$.
Afin de résoudre les problèmes, il faut exprimer ces grandeurs dans des repères appropriés.

### Repère de Frenet

Plutôt que d'exprimer les vecteurs position, vitesse et accélération dans un repère géométrique fixe lié à l'origine O, il est possible d'utiliser un repère mobile, d'origine au point P, défini par:
* le vecteur unitaire tangent $\mathbf {e} _{t}$, dont la direction est celle de la tangente à la trajectoire au point P et le sens est celui du mouvement;
* le vecteur unitaire orthonormal $\mathbf {e} _{n}$, dont la direction est celle de la perpendiculaire à la tangente en P dans le plan osculateur à la trajectoire, dirigé selon le centre de courbure;
* le vecteur unitaire binormal $\mathbf {e} _{b}$, perpendiculaire au plan osculateur donné par $(\mathbf {e} _{t},\mathbf {e} _{n})$, et orienté de telle sorte que les trois vecteur unitaires $(\mathbf {e} _{t},\mathbf {e} _{n},\mathbf {e} _{b})$ forment un trièdre direct, il suffit de poser: $\mathbf {e} _{b}=\mathbf {e} _{t}\times \mathbf {e} _{n}$.
Bien entendu l'expression des vecteurs vitesse et accélération nécessite de dériver les vecteurs de base du trièdre par rapport au temps, et il est qu'alors il vient par exemple ${\frac {d\mathbf {e} _{t}}{dt}}={\frac {d\mathbf {e} _{t}}{ds}}{\frac {ds}{dt}}=v{\frac {d\mathbf {e} _{t}}{ds}}$. Il est possible de montrer que la dérivation des vecteurs de base par rapport à l'abscisse curviligne peut se mettre sous la forme suivante:
    $\mathbf {e} '_{t}={\frac {d\mathbf {e} _{t}}{ds}}=\gamma \mathbf {e} _{n}$, où $\gamma ={\frac {1}{\rho }}$ est la courbure de la trajectoire au point P, inverse du rayon de courbure $\rho$, qui lui correspond au rayon du "cercle tangent instantané" à la trajectoire dans le plan osculateur;
    $\mathbf {e} '_{n}={\frac {d\mathbf {e} _{n}}{ds}}=-\gamma \mathbf {e} _{t}+\tau \mathbf {e} _{b}$, où $\tau$ est la torsion de la courbe, qui traduit la tendance du plan osculateur à s'incliner si la trajectoire n'est pas plane (sinon $\tau =0$).
    $\mathbf {e} '_{b}={\frac {d\mathbf {e} _{b}}{ds}}=-\tau \mathbf {e} _{n}$, bien entendu si la trajectoire est plane cette quantité est nulle et la binormale garde une direction fixe dans d'espace.

### Expressions des vecteurs vitesse et accélération

À partir des expressions précédentes il est possible de démontrer que les vecteurs vitesses et accélération se décomposent ainsi dans le repère de Frenet:
    $\mathbf {v} ={\frac {ds}{dt}}\mathbf {e} _{t}=v\mathbf {e} _{t}$;
    $\mathbf {a} ={\frac {dv}{dt}}\mathbf {e} _{t}+{\frac {v^{2}}{\rho }}\mathbf {e} _{n}$.
Comme le repère local de Frenet est orthonormal, ces expressions correspondent bien aux composantes des deux vecteurs dans ce repère, et il est facile d'exprimer la valeur de la vitesse et de l'accélération en fonction de celles-ci.

## Théorème de Siacci

Dans sa publication originale Siacci propose une décomposition du vecteur accélération selon la direction du vecteur position, appelée direction radiale, notée $\mathbf {e} _{r}$ et la direction tangentielle $\mathbf {e} _{t}$ (cf. figure ci-contre).

### Cas d'une trajectoire plane
#### Enoncé

Dans le cas simple d'une trajectoire contenue à tout moment dans un plan, le théorème de Siacci s'écrit:
    $\mathbf {a} =-{\frac {\gamma v^{2}r}{p}}\mathbf {e} _{r}+{\frac {h}{p^{2}}}{\frac {dh}{ds}}\mathbf {e} _{t}$,
ou encore, en faisant apparaître plutôt le rayon de courbure:
    $\mathbf {a} =-{\frac {v^{2}r}{\rho p}}\mathbf {e} _{r}+{\frac {h}{p^{2}}}{\frac {dh}{ds}}\mathbf {e} _{t}$,
où p est la distance entre la projection orthogonale de l'origine O sur la direction de la tangente à la trajectoire en P, autrement dit la distance à la podaire par rapport à l'origine O, et $h=pv$, qui correspond au moment cinétique par unité de masse du point matériel par rapport à O$^,$.

#### Démonstration

La démonstration du théorème de Siacci est assez facile et consiste à utiliser l'expression du vecteur position dans le trièdre de Serret-Frenet:
    $\mathbf {r} =q\mathbf {e} _{t}-p\mathbf {e} _{n}$,
pour introduire le vecteur unitaire de la direction radiale $\mathbf {e} _{r}={\frac {\mathbf {r} }{r}}$. Il est alors facile d'exprimer le vecteur unitaire normal en fonction de $\mathbf {e} _{t}$ et de $\mathbf {e} _{r}$, il vient:
    $\mathbf {e} _{n}=-{\frac {r}{p}}\mathbf {e} _{r}+{\frac {q}{p}}\mathbf {e} _{t}$, expression valable pour p ≠ 0.
La substitution dans l'expression de Frenet de l'accélération permet d'obtenir sans difficulté l'expression:
    $\mathbf {a} =-{\frac {\gamma rv^{2}}{p}}\mathbf {e} _{r}+\left({\frac {\gamma qv^{2}}{p}}+v{\frac {dv}{ds}}\right)\mathbf {e} _{t}$.
Afin de faire apparaître la quantité h, il convient d'utiliser le fait que $p'={\frac {dp}{ds}}=\gamma q$.
Il vient alors facilement ${\frac {dh}{ds}}={\frac {d(pv)}{ds}}=p'v+p{\frac {dv}{ds}}=\gamma vq+p{\frac {dv}{ds}}={\frac {p}{v}}\left({\frac {\gamma qv^{2}}{p}}+v{\frac {dv}{ds}}\right)={\frac {p^{2}}{h}}\left({\frac {\gamma qv^{2}}{p}}+v{\frac {dv}{ds}}\right)$, et finalement l'expression donnée plus haut du théorème de Siacci dans le cas du mouvement plan.

### Cas d'une trajectoire non-plane

Le théorème de Siacci peut être généralisé au cas où la trajectoire du point matériel P dans le référentiel d'étude n'est pas plane (courbe gauche).
Dans ce cas l'expression précédente est modifiée en ce qui concerne la composante tangentielle et devient:
    $\mathbf {a} =-{\frac {\gamma v^{2}r}{p}}\mathbf {e} _{r}+{\frac {h}{p^{2}}}\left({\frac {dh}{ds}}+\tau bv\right)\mathbf {e} _{t}$,
où b est la distance à la date t entre le plan osculateur et l'origine du repère.
En réalité, l'expression établie lors de la démonstration dans le cas d'une courbe plane, sans faire intervenir le moment cinétique par unité de masse h, est toujours valable:
    $\mathbf {a} =-{\frac {\gamma rv^{2}}{p}}\mathbf {e} _{r}+\left({\frac {\gamma qv^{2}}{p}}+v{\frac {dv}{ds}}\right)\mathbf {e} _{t}$.
Toutefois, dans le cas d'une trajectoire gauche, l'expression de la dérivée curviligne de h est modifiée car dans ce cas $p'=\gamma q-\tau b$, ce qui conduit après réarrangements à l'expression proposée plus haut du théorème de Siacci.

## Applications

Il est important de souligner que le théorème de Siacci donne une décomposition de l'accélération d'un point matériel selon un repère non-orthogonal. Par suite en toute rigueur il est erroné de parler de "composantes" radiale et tangentielle de l'accélération avec cette décomposition.
Le théorème de Siacci peut être utilisé dans le cas d'un champ central de centre O. Dans ce cas la trajectoire est plane et la partie tangentielle de l'accélération sera par hypothèse nulle:
    ${\frac {h}{p^{2}}}{\frac {dh}{ds}}={\frac {1}{2p^{2}}}{d(h^{2})}{ds}=0$, par suite $h'=0$, le moment cinétique garde une valeur fixe tout le long de la trajectoire. En fait cela correspond à une vitesse aréolaire constante (loi des aires), puisque que le moment cinétique par unité de masse est dans ce cas proportionnel à cette vitesse.
Par ailleurs en notant f(r) pour la valeur de l'accélération radiale, il vient d'après le théorème de Siacci la relation:
    ${\frac {\gamma v^{2}r}{p}}={\frac {\gamma h^{2}r}{p^{3}}}=-f(r)$, or il est possible de montrer que $\gamma r={\frac {dp}{dr}}$, par suite il vient en intégrant:
    ${\frac {h^{2}}{2p^{2}}}+\int {f(r)dr}=cte$, relation qui correspond en fait à la conservation de l'énergie.