Projet 5
Validation numérique de la concentration de contraintes d’une plaque trouée par comparaison analytique et éléments finis
Accès immédiat + email envoyé après paiementVoici le contenu du fichier article.md, rédigé avec rigueur scientifique en me basant exclusivement sur les sources théoriques et numériques fournies concernant la plaque trouée.
I) Partie théorique
1) Contexte et problématique
L’étude de la concentration de contraintes au voisinage d’une discontinuité géométrique est un enjeu majeur de la mécanique de la rupture et du dimensionnement des structures. Le problème considéré porte sur une plaque élastique, isotrope et linéaire, de grandes dimensions, percée en son centre d’un trou circulaire de rayon \(a\). La plaque est soumise à une contrainte de traction uniforme \(\sigma\) appliquée à l’infini. Physiquement, la présence du trou perturbe le passage des lignes de force, créant des surcontraintes locales susceptibles d’initier des fissures. L’objectif de ce travail est de quantifier cette perturbation en supposant un état de contraintes planes et l’absence de forces volumiques. Cette analyse, résolue historiquement par Kirsch en 1898, constitue une référence fondamentale pour valider la précision des outils de simulation numérique face à des singularités géométriques.
2) Approche analytique
La résolution analytique repose sur l’utilisation de la fonction d’Airy, notée \(\Phi\), une fonction scalaire dont les dérivées secondes permettent de déduire les composantes du tenseur des contraintes tout en vérifiant l’équilibre statique. Pour que le champ de contrainte soit admissible, cette fonction doit satisfaire l’équation biharmonique. Compte tenu de la géométrie du problème, le passage en coordonnées polaires \((r, \theta)\) est privilégié. La démarche consiste à chercher une solution sous la forme d’une somme d’une partie axisymétrique et d’une partie périodique. En imposant les conditions aux limites — à savoir une surface libre de contrainte au bord du trou (\(r = a\)) et une traction uniforme à l’infini — on détermine les constantes d’intégration. Cette méthode aboutit à la solution de Kirsch, qui démontre que la contrainte ortho-radiale au bord du trou est triplée par rapport à la contrainte appliquée au sommet de la perforation, établissant un facteur de concentration théorique \(K_t = 3\).
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) à l’aide de la bibliothèque FEniCSx. Le problème est traité par une formulation variationnelle qui traduit l’équilibre du système sous forme intégrale, reliant l’énergie de déformation virtuelle au travail des forces extérieures. La géométrie, une plaque carrée de côté \(L\), est discrétisée par un maillage généré avec Gmsh. Un raffinement local crucial est appliqué au bord du trou pour capturer précisément le fort gradient de contrainte prédit par la théorie. Pour l’espace fonctionnel, on utilise des éléments finis vectoriels de Lagrange d’ordre 2 (P2), garantissant une meilleure précision sur les déplacements et, par dérivation, sur les contraintes. Les conditions aux limites numériques bloquent les déplacements rigides (symétries) tandis qu’une traction de 10 MPa est imposée sur le bord supérieur. Le système linéaire résultant est résolu par une méthode de décomposition LU via le solveur PETSc. Enfin, les contraintes sont projetées sur un espace continu pour faciliter l’extraction des valeurs nodales au bord du trou.
2) Résultats numériques
Les simulations numériques permettent de cartographier le champ de contrainte ortho-radiale \(\sigma_{\theta\theta}\) sur l’ensemble du domaine. Les résultats confirment visuellement la localisation de la concentration de contrainte aux sommets du trou, là où la direction est perpendiculaire à l’axe de traction. Pour une traction appliquée de 10 MPa, le solveur identifie une valeur maximale proche de 30 MPa, tandis qu’une zone de compression de l’ordre de -10 MPa est observée à l’équateur du trou. Les visualisations obtenues via PyVista mettent en évidence l’efficacité du raffinement de maillage : la transition entre la zone perturbée et le champ uniforme lointain est fluide et cohérente avec la physique du problème. Ces ordres de grandeur numériques valident l’aptitude du modèle à reproduire le comportement singulier induit par la perforation circulaire.
III) Comparaison analytique – numérique
La confrontation entre la solution théorique de Kirsch et les résultats issus de FEniCSx révèle une excellente cohérence globale. Sur le plan quantitatif, la valeur maximale de la contrainte ortho-radiale obtenue par simulation présente une erreur relative inférieure à 2 % par rapport aux 30 MPa prédits par le calcul analytique. Le facteur de concentration de contrainte numérique \(K_t\) converge très précisément vers la valeur théorique de 3. L’analyse montre également que la perturbation thermique et mécanique s’estompe rapidement : à une distance radiale supérieure à trois fois le rayon du trou, l’écart entre la solution numérique et la traction nominale devient négligeable, conformément aux prévisions théoriques. Cette convergence démontre que le maillage raffiné et l’utilisation d’éléments d’ordre 2 permettent de surmonter les difficultés liées aux gradients élevés près des bords libres.
Conclusion orientée validation
Cette étude valide la robustesse de la simulation par éléments finis pour traiter des problèmes de concentration de contraintes complexes. La parfaite superposition des résultats numériques avec la solution exacte de Kirsch confirme la fiabilité de la formulation variationnelle et du processus de discrétisation employé. Sur le plan pédagogique, cet exercice souligne l’importance de coupler l’analyse théorique (fonction d’Airy) et la vérification numérique pour garantir l’intégrité des structures d’ingénierie. La concordance obtenue prouve que les outils modernes de calcul permettent de capturer des phénomènes locaux avec une grande précision. En conclusion, la maîtrise de ce cas de référence autorise l’extension de la méthode à des géométries plus complexes, telles que des entailles ou des fissures, où les solutions analytiques ne sont plus disponibles.