Analyse élastostatique d'un cylindre en torsion
I) Partie théorique
1) Contexte et problématique
Cette étude porte sur l’analyse mécanique d’un arbre cylindrique de hauteur h et de rayon R, sollicité en torsion pure. L’objectif principal est de déterminer l’état de contraintes internes et les déformations résultantes afin de prévoir les zones de rupture potentielles, notamment pour les matériaux fragiles. Le système est modélisé dans le cadre de l’élastostatique linéarisée, en supposant un matériau homogène, isotrope, et en adoptant l’hypothèse des petites déformations. Le problème physique réside dans la transmission d'un couple moteur sans force résultante, une configuration classique en ingénierie mécanique pour les arbres de transmission. Comprendre la distribution des contraintes, maximales à la périphérie du cylindre, est essentiel pour garantir l'intégrité structurelle du composant sous charge. L’enjeu ici est de fournir une solution de référence permettant de valider les approches numériques plus complexes.
2) Approche analytique
La résolution analytique repose sur l'intégration des équations d'équilibre local en l'absence de forces volumiques. En partant d’un tenseur des contraintes proposé, on démontre que sa divergence est nulle, confirmant que seule la sollicitation aux limites régit le comportement du système. Les conditions aux limites imposent une surface latérale libre d'effort, tandis que les bases du cylindre reçoivent des distributions de contraintes de cisaillement. Le calcul des éléments de réduction montre que la résultante des efforts est nulle, alors que le moment résultant correspond à un couple de torsion pur agissant selon l'axe longitudinal. La démarche permet d’identifier que la contrainte principale maximale se situe sur la peau du cylindre. Cette analyse théorique prédit également une rupture hélicoïdale à 45 degrés pour les matériaux fragiles, une conséquence directe de l'orientation des directions principales des contraintes.
II) Partie numérique (FEM)
1) Méthode numérique employée
L'analyse numérique est conduite via la méthode des éléments finis (FEM) en utilisant la bibliothèque FEniCSx. Le problème est abordé par une formulation variationnelle qui traduit l’équilibre global du système sous forme intégrale, reliant le tenseur des contraintes de Hooke au tenseur des déformations. La géométrie est discrétisée par un maillage tétraédrique généré avec Gmsh, permettant une représentation précise du domaine cylindrique. Pour assurer une convergence optimale et une précision élevée dans le calcul des gradients, des éléments de Lagrange quadratiques d'ordre 2 sont sélectionnés pour l'espace fonctionnel. Les conditions aux limites numériques reproduisent fidèlement le montage physique : une base est parfaitement encastrée (déplacement nul), tandis qu'un champ de déplacement tangentiel de torsion est imposé sur la face opposée. Les propriétés du matériau sont intégrées via les paramètres de Lamé, calculés à partir du module d'Young et du coefficient de Poisson. L'intérêt de cette approche numérique réside dans sa capacité à extraire des champs de contraintes de von Mises et des rotations en tout point du volume.
2) Résultats numériques
Les simulations fournissent une cartographie complète des déplacements et des contraintes. Le résultat principal met en évidence une rotation maximale selon l'axe longitudinal au niveau de la face libre. Pour les paramètres d'application, la rotation maximale simulée atteint environ 0,0101 radian. En termes de sollicitations, les contraintes équivalentes de von Mises sont calculées pour identifier les risques de plastification, tandis que les contraintes principales sont extraites pour la validation théorique. L'ordre de grandeur de la contrainte principale maximale numérique est d'environ 99,9 MPa. La visualisation 3D des résultats confirme que les contraintes les plus élevées sont localisées sur la surface extérieure du cylindre, conformément aux prévisions physiques, et décroissent linéairement vers le centre de la section transversale.
Notebook Cylindre en torsion
Notebook Jupyter complet avec code source pour modéliser un cylindre en torsion (simulation FEM FEniCS + explications).
💡 Reproduire cette simulation
Cette étude est entièrement reproductible dans FEneC, un environnement de calcul en Python dédié à l'analyse et à la simulation par la méthode des éléments finis. Les notebooks complets sont inclus.
