Calcul d’une contrainte autour d’un cylindre FreeFEM
Outil premium pour estimer le champ de contraintes autour d’un cylindre ou d’un trou circulaire à partir de la solution analytique de Kirsch, utile pour vérifier un modèle FreeFEM en élasticité plane.
Résultats
Saisissez vos données puis cliquez sur le bouton pour afficher σrr, σθθ, τrθ et σvon Mises.
Guide expert : comment réaliser le calcul d’une contrainte autour d’un cylindre avec FreeFEM
Le calcul d’une contrainte autour d’un cylindre est un cas d’école fondamental en mécanique des milieux continus, en géomécanique, en résistance des matériaux et en simulation numérique par éléments finis. Dans la pratique, on l’utilise pour modéliser un trou dans une plaque, un puits dans une roche, une inclusion circulaire, un alésage mécanique, ou encore le voisinage d’un conduit soumis à des contraintes lointaines. Si vous utilisez FreeFEM pour résoudre ce problème, vous cherchez généralement à obtenir un champ de déplacements, puis à en déduire les composantes du tenseur des contraintes au voisinage de la paroi circulaire.
Ce sujet est important parce que la présence d’un cylindre ou d’un trou perturbe très fortement la distribution des contraintes. Même dans un matériau linéaire, homogène et isotrope, une traction uniforme appliquée au loin peut créer une concentration de contraintes significative au bord du cylindre. C’est précisément ce phénomène qui explique l’amorçage de fissures, la plastification locale ou les erreurs d’interprétation si le maillage FreeFEM n’est pas assez fin autour de la cavité.
1. Contexte physique du problème
On considère en général une plaque très grande par rapport au rayon du cylindre, de sorte que les frontières extérieures n’influencent que faiblement le voisinage du trou. Dans cette hypothèse dite de plaque quasi infinie, les équations analytiques donnent directement les composantes de contrainte en coordonnées polaires. En FreeFEM, on reconstruit numériquement ce même problème en 2D, en élasticité plane, avec un domaine annulaire ou un grand rectangle percé d’un trou.
Les grandeurs les plus importantes sont :
- le rayon a du cylindre ou du trou ;
- la distance radiale r au centre, avec la contrainte évaluée pour r ≥ a ;
- l’angle θ, qui permet d’identifier la position autour de la paroi ;
- les contraintes lointaines σx, σy et éventuellement τxy ;
- une éventuelle pression interne pi sur la paroi du cylindre.
Dans la plupart des validations FreeFEM, on commence par le cas le plus simple : traction uniaxiale, pas de cisaillement, pas de pression interne. On vérifie ensuite que les résultats numériques reproduisent bien la variation angulaire de la contrainte circonférentielle et la concentration maximale en bord de trou.
2. Pourquoi utiliser FreeFEM pour ce calcul
FreeFEM est particulièrement efficace pour ce type de problème car il permet de créer rapidement des géométries 2D, d’imposer des conditions aux limites mécaniques, d’affiner localement le maillage et d’extraire les champs de contraintes. Dans un flux de travail d’ingénierie, FreeFEM sert souvent à :
- tester un cas analytique simple avant un problème réel plus complexe ;
- analyser l’effet de la taille du domaine extérieur ;
- étudier l’influence du maillage près du cylindre ;
- ajouter ensuite anisotropie, non-linéarités, couplages thermo-mécaniques ou géométries non circulaires.
Le grand avantage est de pouvoir comparer la solution éléments finis à une référence théorique connue. Lorsque l’écart est faible, vous gagnez en confiance avant de passer à un modèle non analytique. Lorsque l’écart est fort, le problème vient souvent d’un bord extérieur trop proche, d’une conversion erronée entre déplacements et contraintes, d’un maillage trop grossier ou d’une mauvaise interprétation des conventions de signe.
3. Formulation de base utilisée par le calculateur
Le calculateur ci-dessus applique une forme usuelle de la solution de Kirsch en coordonnées polaires pour une plaque infinie contenant un trou circulaire. Il estime :
- la contrainte radiale σrr ;
- la contrainte circonférentielle σθθ ;
- le cisaillement polaire τrθ ;
- une contrainte équivalente de von Mises en hypothèse de contrainte plane.
Cette approche est très utile pour comparer vos sorties FreeFEM à des valeurs attendues à un point particulier. Si, par exemple, vous calculez en FreeFEM un point situé sur la frontière du trou à θ = 90° sous traction uniaxiale σx, vous devriez retrouver une contrainte circonférentielle proche du pic théorique. Dans le cas standard d’un trou dans une plaque infinie sous traction uniaxiale, la concentration maximale au bord du trou vaut en théorie Kt = 3.
4. Interprétation des résultats autour de la paroi
L’erreur la plus fréquente est de regarder seulement une valeur isolée au lieu de comprendre la distribution complète autour du cylindre. La contrainte varie avec l’angle. Sur certaines directions, la paroi est fortement sollicitée ; sur d’autres, la contrainte peut devenir nettement plus faible. Cette anisotropie géométrique est normale même si le matériau reste isotrope.
Le graphique généré dans le calculateur représente la variation de σθθ en fonction de l’angle. Cette courbe est précieuse pour une validation FreeFEM car elle vous permet de comparer une tendance globale, et pas uniquement un point. Si votre simulation reproduit bien l’amplitude, la position des maxima et minima, ainsi que la symétrie attendue, votre modèle est probablement cohérent.
5. Paramètres matériaux généralement associés à l’étude
Le champ de contraintes autour d’un cylindre dépend de la géométrie et du chargement, mais l’analyse complète en éléments finis nécessite aussi des propriétés mécaniques telles que le module d’Young E et le coefficient de Poisson ν. Les ordres de grandeur suivants sont couramment utilisés en simulation structurelle et géomécanique initiale :
| Matériau | Module d’Young E | Coefficient de Poisson ν | Contrainte admissible ou repère usuel |
|---|---|---|---|
| Acier de construction | 200 à 210 GPa | 0,27 à 0,30 | Limite d’élasticité typique 235 à 355 MPa |
| Aluminium 6061-T6 | 68 à 69 GPa | 0,33 | Limite d’élasticité typique autour de 276 MPa |
| Béton ordinaire | 25 à 35 GPa | 0,15 à 0,22 | Résistance en compression souvent 20 à 40 MPa |
| Roche dense | 30 à 80 GPa | 0,15 à 0,30 | Très variable selon texture, porosité et confinement |
Ces valeurs ne changent pas la structure analytique de Kirsch, mais elles deviennent essentielles dès que vous passez de la formule fermée à la simulation FreeFEM complète, puisque les déplacements, l’énergie élastique et la conversion déplacement-contrainte y dépendent explicitement.
6. Effet du maillage dans FreeFEM
Le maillage est déterminant. Près du bord circulaire, le gradient de contrainte peut être élevé ; un maillage trop grossier sous-estime alors le pic de contrainte et lisse artificiellement la courbe angulaire. En pratique, on recommande un raffinement local au voisinage du cylindre et plusieurs couches d’éléments réguliers dans la zone où la concentration de contraintes est la plus forte.
Le tableau suivant donne des ordres de grandeur de convergence fréquemment observés lors d’une validation numérique par rapport à la solution analytique, pour un cas linéaire bien posé :
| Qualité de maillage près du trou | Nombre d’éléments sur la circonférence | Erreur typique sur σθθ max | Usage recommandé |
|---|---|---|---|
| Grossier | 16 à 24 | 5 % à 15 % | Pré-étude rapide |
| Intermédiaire | 32 à 64 | 1 % à 5 % | Validation courante |
| Fin | 96 à 160 | souvent < 1 % | Benchmark et publication technique |
Il ne s’agit pas d’une loi universelle, mais d’une règle d’ingénierie raisonnable. L’erreur dépend également de la formulation des éléments, de la taille du domaine extérieur, du post-traitement des contraintes et de la manière dont vous interpolez la solution sur la frontière.
7. Méthodologie de validation d’un modèle FreeFEM
- Construire une géométrie avec un trou de rayon connu dans un domaine suffisamment grand.
- Imposer les contraintes ou déplacements à distance de façon cohérente avec le cas analytique choisi.
- Raffiner le maillage autour du cylindre.
- Résoudre le problème d’élasticité plane.
- Extraire les contraintes sur la paroi ou à différents rayons.
- Comparer point par point avec la solution analytique du calculateur.
- Vérifier la convergence en affinant progressivement le maillage.
Cette démarche doit être répétée tant que la différence entre FreeFEM et l’analytique n’est pas stabilisée. En ingénierie, on ne valide jamais sérieusement un modèle éléments finis sur un seul niveau de maillage.
8. Erreurs classiques à éviter
- Confondre le rayon a et la distance d’évaluation r.
- Évaluer la solution pour r < a, ce qui n’a pas de sens dans le cas d’un trou vide.
- Utiliser une convention de signe inversée entre pression et traction.
- Oublier que les contraintes en coordonnées cartésiennes et polaires ne sont pas identiques.
- Employer un domaine extérieur trop petit, ce qui modifie le champ lointain.
- Comparer une contrainte nodale lissée à une contrainte analytique ponctuelle sans précaution.
9. Quand le modèle analytique n’est plus suffisant
La solution analytique reste une référence remarquable, mais elle a ses limites. Dès que vous introduisez une forme non circulaire, des matériaux anisotropes, un comportement plastique, de grandes déformations, des fissures, un couplage thermique ou des conditions aux limites complexes, FreeFEM redevient indispensable. Dans ces cas, l’analytique sert surtout de base de validation sur un sous-problème simplifié.
Autrement dit, le bon usage n’est pas d’opposer solution analytique et simulation numérique, mais de les combiner intelligemment. L’analytique sécurise la phase de vérification ; le calcul éléments finis traite ensuite la complexité réelle du projet.
10. Lecture rapide des résultats du calculateur
Après calcul, vous obtenez :
- σrr : contrainte radiale, souvent égale à la pression opposée sur la paroi ;
- σθθ : contrainte circonférentielle, généralement la plus surveillée autour du cylindre ;
- τrθ : cisaillement polaire ;
- σvon Mises : indicateur utile pour une première vérification vis-à-vis d’un critère de plasticité.
Si vous êtes au bord du trou, la contrainte circonférentielle est souvent la grandeur la plus critique. En vérification FreeFEM, c’est aussi celle que l’on trace le plus souvent sur toute la circonférence pour mettre en évidence le facteur de concentration et la bonne symétrie de la réponse numérique.
11. Ressources d’autorité pour approfondir
MIT OpenCourseWare – Finite Element Analysis of Solids and Fluids
NASA – Finite Element Analysis Reference
NIST – Références et données scientifiques pour les matériaux et la modélisation
12. Conclusion pratique
Pour réussir un calcul d’une contrainte autour d’un cylindre avec FreeFEM, la meilleure stratégie consiste à partir d’un cas analytique propre, à comparer systématiquement vos résultats à la solution de Kirsch, puis à augmenter progressivement le niveau de complexité. Cette discipline de vérification est ce qui distingue un modèle numérique convaincant d’un simple calcul visuellement plausible. Le calculateur proposé ici peut ainsi servir de point de contrôle immédiat pour vos cas de traction lointaine, de cisaillement et de pression interne dans un cadre d’élasticité plane.