Projet 2

Validation numérique de l’élastostatique d’un cylindre en torsion par comparaison analytique et éléments finis

Accès immédiat + email envoyé après paiement

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.

III) Comparaison analytique – numérique

La confrontation des deux approches révèle une convergence remarquable entre les modèles. L’analyse des résultats montre que l’angle de torsion théorique calculé à 0,010140 rad est très proche de la valeur simulée de 0,010098 rad. L’erreur relative sur la composante de rotation principale est limitée à 0,41 %, et tombe même à 0,04 % lorsqu’on considère la norme globale de la rotation. Concernant les contraintes, la valeur principale maximale théorique de 101,40 MPa concorde avec les 99,90 MPa obtenus par le solveur éléments finis. Ces légers écarts s’expliquent par la discrétisation du maillage et l’approximation numérique inhérente à la méthode FEM. Néanmoins, la cohérence des ordres de grandeur et la distribution spatiale des efforts valident sans ambiguïté le modèle numérique. Cette étape de validation confirme que le simulateur capture avec précision les phénomènes de cisaillement et de torsion.

Conclusion orientée validation

Cette étude démontre l’efficacité de la méthode des éléments finis pour reproduire des solutions analytiques de référence en élasticité linéaire. La comparaison rigoureuse entre les calculs théoriques et les résultats issus de FEniCSx prouve que l’erreur numérique peut être maintenue en dessous de 1 % avec un maillage et un ordre d’interpolation adaptés. Sur le plan pédagogique, ce travail souligne l’importance de maîtriser les directions principales de contraintes pour interpréter les faciès de rupture réels. L’intérêt majeur réside dans la validation du modèle, qui peut désormais être étendu à des géométries plus complexes, comme des sections elliptiques, où les solutions analytiques sont plus ardues. L’utilisation d’un moteur d’éléments finis moderne permet ainsi d’automatiser ces analyses de structure avec une fiabilité garantie pour l’ingénierie.