Analyse d'une poutre console : Euler-Bernoulli face aux Éléments Finis de FEniCs-Link
I) Partie théorique
1) Contexte et problématique
L’étude de la poutre en console, ou cantilever, constitue l'un des problèmes les plus emblématiques de la mécanique des structures. Ses origines remontent à 1638, lorsque Galilée proposa, dans son ouvrage Discours concernant deux sciences nouvelles, une première analyse d’une poutre encastrée soumise à une charge à son extrémité. Bien que son modèle initial ait été approximatif, il introduisit le concept fondamental du levier avec une zone critique située à l’encastrement.
Aujourd'hui, ce problème qui consiste en une poutre de longueur L, parfaitement encastrée à une extrémité et soumise à un chargement ponctuel à son extrémité libre est traité dans le cadre de la théorie d’Euler-Bernoulli, qui suppose un matériau homogène, isotrope et un comportement élastique linéaire. L’enjeu physique est de prédire la flèche et les contraintes internes pour garantir l'intégrité de structures en porte-à-faux. Cette étude vise à confronter la rigueur du calcul différentiel aux capacités de la simulation numérique moderne pour valider la précision des modèles de conception utilisés par les ingénieurs
2) Approche analytique
La résolution analytique s'appuie sur les équations fondamentales de l’équilibre statique pour déterminer les réactions d'appui à l'encastrement, incluant une force verticale et un moment de rappel. La méthode des coupures est ensuite employée, en isolant une section droite à l'abscisse x, pour définir le torseur des efforts internes le long de la poutre, identifiant ainsi l'effort tranchant V(x) et le moment fléchissant M(x). En vertu du principe de Saint-Venant, ces efforts permettent de décrire localement l'état de sollicitation de la structure.
L'équation différentielle de la ligne moyenne, EIy′′(x)=M(x), reliant le moment fléchissant à la courbure, est alors intégrée deux fois pour obtenir l'expression de la flèche y(x). Cette approche permet également de calculer les contraintes normales de flexion, qui atteignent leur paroxysme aux fibres extrêmes de la section au niveau du support. Cette solution théorique exacte sert de référence absolue pour l'évaluation de la performance du solveur numérique.
II) Partie numérique (FEM)
1) Méthode numérique employée
L'analyse numérique est réalisée via la méthode des éléments finis (FEM) intégrée dans un environnement de calcul reproductible : FEniCS-Link.
Un webinaire sera prochainement organisé pour présenter une démonstration complète et une prise en main approfondie.
Accédez au projet FEniCs-Link en open source dès maintenant 👉
alfenec
Cette méthode repose sur une formulation variationnelle des équations de l'élasticité, permettant de transformer le problème différentiel continu en un système algébrique discret. Contrairement à la théorie 1D simplifiée, la simulation FEM traite ici la poutre comme un volume continu discrétisé par des éléments de type Lagrange, permettant de capturer les effets de déformation locaux. La formulation variationnelle du problème est mise en œuvre pour relier le champ de contraintes de Hooke aux déplacements nodaux. Les conditions de Dirichlet sont appliquées pour modéliser l'encastrement strict sur une face complète du volume, tandis qu'une charge ponctuelle est simulée sur la face opposée, reproduisant fidèlement le cas théorique.
L'intérêt de cette approche numérique, au-delà de la simple validation, réside dans sa capacité à explorer rapidement des paramètres physiques et à automatiser des workflows complexes.
2) Résultats numériques
Les simulations numériques fournissent une vision globale du comportement de la poutre sous charge. La déformée obtenue montre un déplacement vertical progressif dont la valeur maximale se situe, comme attendu, à l'extrémité libre de la console. Pour la configuration étudiée, la flèche maximale simulée atteint environ 0,012 mm.
Notebook Poutre Console
Notebook Jupyter complet avec code source pour modéliser une poutre console (simulation FEM FEniCS + données FEniCS-Link).
Parallèlement, la visualisation tridimensionnelle du champ de contraintes de von Mises met en lumière la distribution spatiale des sollicitations au sein du matériau. Les valeurs de contraintes les plus élevées sont clairement identifiées au voisinage immédiat de l'encastrement, atteignant environ 0,68 MPa sur les faces supérieures et inférieures.
Ces ordres de grandeur confirment que le modèle numérique capture fidèlement la localisation des zones de traction et de compression maximales prédites par la théorie classique de Navier-Bernoulli.
💡 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.
🔒 L’accès à ce webinaire est réservé aux membres MecaNews Prime.
👉 Inscrivez-vous gratuitement pour recevoir le lien d’accès et être informé de la date.
