Projet 1

Comportement élastostatique d’un cylindre soumis à son poids propre

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

I) Partie théorique

1) Contexte et problématique

Cette étude porte sur le comportement mécanique d’un cylindre vertical soumis à l’action de son propre poids. L’objectif est de déterminer les champs de contraintes et de déplacements induits par la gravité au sein de la structure. Le système est un cylindre de rayon R et de hauteur H, modélisé comme un milieu élastique, homogène et isotrope. L’analyse s’inscrit dans le cadre de l’élastostatique linéarisée en supposant l’hypothèse des petites déformations. Physiquement, ce problème permet de comprendre comment une structure massive se déforme sous sa propre charge, un aspect fondamental pour la conception d’ouvrages de génie civil ou de composants mécaniques de grande dimension. L’enjeu est ici de valider la solution théorique par une approche numérique afin de garantir la fiabilité des modèles de simulation utilisés en ingénierie.

2) Approche analytique

La résolution analytique repose sur une approche en contraintes [1.1]. La démarche débute par l’établissement des équations d’équilibre local, où la divergence du tenseur des contraintes doit compenser la force volumique de pesanteur. Pour simplifier la résolution, le problème initialement encastré est transformé en un problème de type 2, où la base inférieure est soumise à une distribution de forces surfaciques équivalente au poids total du cylindre. En utilisant un champ de contraintes unidirectionnel selon l’axe vertical, nous intégrons les équations d’équilibre en tenant compte des conditions aux limites : la surface latérale et la face supérieure sont libres d’efforts. Une fois le champ de contraintes déterminé, la loi de comportement de Hooke permet de déduire le tenseur des déformations. Enfin, l’intégration de ces déformations, en satisfaisant les équations de compatibilité géométrique, fournit le champ de déplacement global. Cette méthode révèle que la face supérieure adopte une courbure parabolique due à l’effet du coefficient de Poisson.

II) Partie numérique (FEM)

1) Méthode numérique employée

L’étude numérique est réalisée à l’aide de la méthode des éléments finis (FEM), mise en œuvre via la bibliothèque FEniCSx. Cette approche repose sur la formulation variationnelle des équations de l’élasticité linéaire. Le domaine cylindrique est discrétisé en un maillage non structuré composé d’éléments tétraédriques, dont la finesse peut être ajustée pour tester la convergence des résultats. Contrairement à l’approche analytique simplifiée, le modèle numérique impose une condition d’encastrement sur la base inférieure, se traduisant par un déplacement nul. Les propriétés matériau correspondent à celles d’un acier standard, avec un module de Young de 210 GPa et un coefficient de Poisson de 0,3. L’intérêt de la méthode numérique réside dans sa capacité à traiter des conditions aux limites complexes et à fournir une solution sur l’ensemble du volume discrétisé.

2) Résultats numériques

Les simulations numériques permettent de visualiser précisément la distribution des déplacements et des contraintes au sein du cylindre. Le résultat majeur est le déplacement vertical, noté uz, qui atteint sa valeur maximale au centre de la face supérieure. Pour un cylindre de 50 cm de hauteur, les ordres de grandeur obtenus pour le déplacement maximal sont de l’ordre de 4,51e-08 mètres. Les visualisations confirment le gradient de déplacement le long de l’axe vertical, avec une base fixe et une compression progressive du matériau sous l’effet de la masse située au-dessus de chaque section. Le modèle numérique reproduit fidèlement la tendance physique attendue d’un tassement gravitaire et la déformée de la surface supérieure.

III) Comparaison analytique – numérique

La comparaison entre les deux approches montre une excellente cohérence globale, confirmant la validité du modèle. Pour le déplacement au centre de la face supérieure, la valeur théorique calculée est de -0,458 nm (normalisée), tandis que la simulation par éléments finis fournit une valeur de -0,4519 nm pour le maillage le plus fin. L’erreur relative observée est extrêmement faible, s’établissant à 1,36 % pour un maillage de 73 586 éléments. On observe que l’écart diminue progressivement à mesure que le maillage est raffiné, illustrant le phénomène de convergence numérique. La légère différence résiduelle s’explique par la nature discrète de la méthode des éléments finis et les différences de conditions aux limites (encastrement strict vs efforts imposés). Cette confrontation valide la précision de l’outil numérique pour ce type d’analyse mécanique.

Conclusion orientée validation

Cette étude démontre la robustesse de la simulation numérique par éléments finis pour traiter des problèmes d’élastostatique fondamentale. La convergence de l’erreur relative vers des valeurs proches de 1 % atteste de la précision du modèle face à une solution analytique exacte. Sur le plan pédagogique, cet exercice souligne l’importance de la vérification croisée entre théorie et numérique pour assurer la fiabilité des calculs de structure. La concordance des résultats prouve que la méthode est capable de capturer des phénomènes physiques de très faible amplitude. Enfin, cette démarche illustre comment un moteur d’éléments finis moderne permet de reproduire et d’automatiser l’analyse de structures complexes sous chargement gravitaire, facilitant la validation de modèles d’ingénierie plus élaborés.