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 $𝐫 = 𝐎 𝐏 (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:

𝐯 = \frac{d 𝐫}{d t}

;

𝐚 = \frac{d 𝐯}{d t} = \frac{d^{2} 𝐫}{d t^{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 $𝐞_{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 $𝐞_{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 $𝐞_{b}$ , perpendiculaire au plan osculateur donné par $(𝐞_{t}, 𝐞_{n})$ , et orienté de telle sorte que les trois vecteur unitaires $(𝐞_{t}, 𝐞_{n}, 𝐞_{b})$ forment un trièdre direct, il suffit de poser: $𝐞_{b} = 𝐞_{t} \times 𝐞_{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 𝐞_{t}}{d t} = \frac{d 𝐞_{t}}{d s} \frac{d s}{d t} = v \frac{d 𝐞_{t}}{d s}$ . 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:

𝐞'_{t} = \frac{d 𝐞_{t}}{d s} = γ 𝐞_{n}

, où

γ = \frac{1}{ρ}

est la courbure de la trajectoire au point P, inverse du rayon de courbure

ρ

, qui lui correspond au rayon du "cercle tangent instantané" à la trajectoire dans le plan osculateur;

𝐞'_{n} = \frac{d 𝐞_{n}}{d s} = - γ 𝐞_{t} + τ 𝐞_{b}

, où

τ

est la torsion de la courbe, qui traduit la tendance du plan osculateur à s'incliner si la trajectoire n'est pas plane (sinon

τ = 0

).

𝐞'_{b} = \frac{d 𝐞_{b}}{d s} = - τ 𝐞_{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:

𝐯 = \frac{d s}{d t} 𝐞_{t} = v 𝐞_{t}

;

𝐚 = \frac{d v}{d t} 𝐞_{t} + \frac{v^{2}}{ρ} 𝐞_{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

Repérage du mouvement du point matériel dans l'espace.

Dans sa publication originale^[1] Siacci propose une décomposition du vecteur accélération selon la direction du vecteur position, appelée direction radiale, notée $𝐞_{r}$ et la direction tangentielle $𝐞_{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:

𝐚 = - \frac{γ v^{2} r}{p} 𝐞_{r} + \frac{h}{p^{2}} \frac{d h}{d s} 𝐞_{t}

,

ou encore, en faisant apparaître plutôt le rayon de courbure:

𝐚 = - \frac{v^{2} r}{ρ p} 𝐞_{r} + \frac{h}{p^{2}} \frac{d h}{d s} 𝐞_{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 = p v$ , qui correspond au moment cinétique par unité de masse du point matériel par rapport à O^[2]Modèle:,^[3].

Démonstration

La démonstration du théorème de Siacci est assez facile^[4] et consiste à utiliser l'expression du vecteur position dans le trièdre de Serret-Frenet:

𝐫 = q 𝐞_{t} - p 𝐞_{n}

,

pour introduire le vecteur unitaire de la direction radiale $𝐞_{r} = \frac{𝐫}{r}$ . Il est alors facile d'exprimer le vecteur unitaire normal en fonction de $𝐞_{t}$ et de $𝐞_{r}$ , il vient:

𝐞_{n} = - \frac{r}{p} 𝐞_{r} + \frac{q}{p} 𝐞_{t}

, expression valable pour Modèle:Nobr.

La substitution dans l'expression de Frenet de l'accélération permet d'obtenir sans difficulté l'expression:

𝐚 = - \frac{γ r v^{2}}{p} 𝐞_{r} + (\frac{γ q v^{2}}{p} + v \frac{d v}{d s}) 𝐞_{t}

.

Afin de faire apparaître la quantité h, il convient d'utiliser le fait que $p^{'} = \frac{d p}{d s} = γ q$ ^[5].

Il vient alors facilement $\frac{d h}{d s} = \frac{d (p v)}{d s} = p^{'} v + p \frac{d v}{d s} = γ v q + p \frac{d v}{d s} = \frac{p}{v} (\frac{γ q v^{2}}{p} + v \frac{d v}{d s}) = \frac{p^{2}}{h} (\frac{γ q v^{2}}{p} + v \frac{d v}{d s})$ , 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)^[6].

Dans ce cas l'expression précédente est modifiée en ce qui concerne la composante tangentielle et devient:

𝐚 = - \frac{γ v^{2} r}{p} 𝐞_{r} + \frac{h}{p^{2}} (\frac{d h}{d s} + τ b v) 𝐞_{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:

𝐚 = - \frac{γ r v^{2}}{p} 𝐞_{r} + (\frac{γ q v^{2}}{p} + v \frac{d v}{d s}) 𝐞_{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^{'} = γ q - τ 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{d h}{d s} = \frac{1}{2 p^{2}} d (h^{2}) d s = 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{γ v^{2} r}{p} = \frac{γ h^{2} r}{p^{3}} = - f (r)

, or il est possible de montrer^[4] que

γ r = \frac{d p}{d r}

, par suite il vient en intégrant:

\frac{h^{2}}{2 p^{2}} + \int f (r) d r = c t e

, relation qui correspond en fait à la conservation de l'énergie.

Notes et références

Modèle:Références

Bibliographie

Voir aussi

Articles connexes

Liens externes

Modèle:Portail

↑ Modèle:Harvsp
↑ En effet, il est possible d'écrire pour le vecteur position $𝐫 = q 𝐞_{t} - p 𝐞_{n}$ , la distance q étant le projeté orthogonal du vecteur position sur la tangente à la trajectoire à la date considérée (cf. figure). Il vient alors pour le moment cinétique par rapport à l'origine du point matériel (de masse notée m):
$𝐋_{O} = m 𝐫 \times 𝐯 = m (q 𝐞_{t} - p 𝐞_{n}) \times (v 𝐞_{t}) = m p v 𝐞_{b} = m h 𝐞_{b}$ ,
où il a été fait usage de l'expression du vecteur vitesse dans le trièdre de Frenet et du fait que $𝐞_{b} = 𝐞_{t} \times 𝐞_{n}$ .
↑ Dans ce contexte de la mécanique céleste, il s'agit en fait de la constante des aires généralement notée C: cf. article mouvement képlerien
↑ ^4,0 et ^4,1 Modèle:Harvsp
↑ Ceci peut être facilement démontrer replaçant l'expression du vecteur position $𝐫 = q 𝐞_{t} - p 𝐞_{n}$ dans la définition du vecteur unitaire tangent $𝐞_{t} = \frac{d 𝐫}{d s}$ , il vient facilement $𝐞_{t} = (q^{'} + γ p) 𝐞_{t} + (γ q - p^{'}) 𝐞_{n}$ , il est évident que cette égalité n'est vérifiée que si $q^{'} = 1 - γ p$ et $p^{'} = γ p$ .
↑ Modèle:Harvsp

[Siacci1879a-1] Modèle:Harvsp

[2] En effet, il est possible d'écrire pour le vecteur position $𝐫 = q 𝐞_{t} - p 𝐞_{n}$ , la distance q étant le projeté orthogonal du vecteur position sur la tangente à la trajectoire à la date considérée (cf. figure). Il vient alors pour le moment cinétique par rapport à l'origine du point matériel (de masse notée m):
$𝐋_{O} = m 𝐫 \times 𝐯 = m (q 𝐞_{t} - p 𝐞_{n}) \times (v 𝐞_{t}) = m p v 𝐞_{b} = m h 𝐞_{b}$ ,
où il a été fait usage de l'expression du vecteur vitesse dans le trièdre de Frenet et du fait que $𝐞_{b} = 𝐞_{t} \times 𝐞_{n}$ .

[3] Dans ce contexte de la mécanique céleste, il s'agit en fait de la constante des aires généralement notée C: cf. article mouvement képlerien

[Casey2011-4] 4,0 et ^4,1 Modèle:Harvsp

[5] Ceci peut être facilement démontrer replaçant l'expression du vecteur position $𝐫 = q 𝐞_{t} - p 𝐞_{n}$ dans la définition du vecteur unitaire tangent $𝐞_{t} = \frac{d 𝐫}{d s}$ , il vient facilement $𝐞_{t} = (q^{'} + γ p) 𝐞_{t} + (γ q - p^{'}) 𝐞_{n}$ , il est évident que cette égalité n'est vérifiée que si $q^{'} = 1 - γ p$ et $p^{'} = γ p$ .

[Siacci1879b-6] Modèle:Harvsp

[1]

[2]

[3]

[4]

[5]

[6]

Théorème de Siacci

Sommaire

Notions préliminaires : repère de Serret-Frenet

Définitions générales

Repère de Frenet

Expressions des vecteurs vitesse et accélération

Théorème de Siacci

Cas d'une trajectoire plane

Enoncé

Démonstration

Cas d'une trajectoire non-plane

Applications

Notes et références

Bibliographie

Voir aussi

Articles connexes

Liens externes

Menu de navigation

Théorème de Siacci

Notions préliminaires : repère de Serret-Frenet

Définitions générales

Repère de Frenet

Expressions des vecteurs vitesse et accélération

Théorème de Siacci

Cas d'une trajectoire plane

Enoncé

Démonstration

Cas d'une trajectoire non-plane

Applications

Notes et références

Bibliographie

Voir aussi

Articles connexes

Liens externes

Menu de navigation

Rechercher