Towards Automatic Stress Analysis using Scaled Boundary Finite Element Method with Quadtree Mesh of High-order Elements

Cet article présente une technique automatisée d'analyse des contraintes et de la fissuration utilisant la méthode des éléments finis à frontière mise à l'échelle (SBFEM) avec des maillages quadtree d'éléments d'ordre élevé, permettant de modéliser précisément les singularités de contrainte et les frontières courbes sans enrichissement asymptotique ni traitement spécial des nœuds pendus.

L'essentiel

Imaginez que vous essayez de prédire comment une structure (comme un pont, une pièce de machine ou même un réacteur nucléaire) va réagir lorsqu'on la pousse ou la tire. Pour le faire, les ingénieurs utilisent une méthode appelée Analyse par Éléments Finis. C'est comme découper l'objet en millions de petits morceaux (des "briques") pour calculer les contraintes sur chacune d'elles.

Le problème ? Découper un objet complexe en briques parfaites est un cauchemar pour les ordinateurs. Surtout si l'objet a des formes bizarres, des trous ronds ou des fissures.

Voici l'histoire de cette recherche, racontée simplement :

1. Le Problème : Le Puzzle Impossible

Traditionnellement, pour modéliser une forme ronde ou une fissure, il faut créer un maillage (un réseau de briques) très précis.

• Le dilemme : Si vous voulez que le contour d'un cercle soit lisse, vous devez utiliser des briques minuscules tout autour. Cela crée des millions de pièces, ce qui ralentit énormément l'ordinateur.
• Le casse-tête des "nœuds pendus" : Quand vous avez des bignes de tailles différentes (grosses d'un côté, petites de l'autre), elles ne s'emboîtent pas toujours parfaitement. Il reste des "nœuds pendus" (des coins qui flottent dans le vide sans être connectés). Les méthodes classiques doivent faire des calculs compliqués pour "coller" ces nœuds ensemble, comme essayer de faire tenir un puzzle avec des pièces de tailles différentes.

2. La Solution Magique : La Méthode "Squelette" (SBFEM)

Les auteurs de ce papier (Hou Man et son équipe) ont eu une idée brillante : arrêter de découper l'intérieur de l'objet en briques.

Au lieu de cela, ils utilisent une technique appelée Méthode des Éléments Finis à Frontière Échelle (SBFEM).

• L'analogie du projecteur : Imaginez que vous avez une lampe (le centre de l'objet) et que vous projetez la lumière sur les bords de la pièce. Vous n'avez besoin de modéliser que le contour (les murs), pas l'intérieur.
• La puissance : Cette méthode permet de calculer mathématiquement ce qui se passe à l'intérieur, même si la forme est bizarre ou si une fissure traverse la pièce. Elle gère les "points de rupture" (comme le bout d'une fissure) avec une précision incroyable, sans avoir besoin de micro-briques partout.

3. L'Outil d'Automatisation : Le "Quadtree" (L'Arbre Carré)

Pour rendre tout cela automatique, ils utilisent une structure appelée Quadtree.

• L'analogie du zoom : Imaginez une image numérique. Au début, c'est un seul carré géant. Si une zone est intéressante (par exemple, autour d'un trou), l'ordinateur divise ce carré en 4 petits carrés. Si l'un de ces 4 est encore trop gros, il le divise à nouveau en 4, et ainsi de suite. C'est comme un zoom automatique qui ne s'arrête que là où c'est nécessaire.
• Le résultat : Vous obtenez un maillage où les carrés sont gros dans les zones simples et minuscules près des détails complexes, le tout généré automatiquement.

4. La Grande Révolution : Combiner les deux

C'est ici que le papier devient génial. Ils ont combiné le Zoom automatique (Quadtree) avec la Méthode du projecteur (SBFEM).

• Gestion des "nœuds pendus" : Dans un maillage automatique, il y a souvent des carrés de tailles différentes qui se touchent. Normalement, cela crée des trous ou des nœuds pendus. Avec leur méthode, chaque carré (même celui qui a des voisins de tailles différentes) est traité comme une forme polygonale unique. Plus besoin de tricher pour coller les pièces ensemble ! L'ordinateur gère tout naturellement.
• Les bords courbes : Au lieu de découper des milliers de petits carrés pour imiter un cercle (ce qui donne un effet "escalier"), ils utilisent des éléments d'ordre élevé. Imaginez que le bord de votre carré n'est pas une ligne droite, mais une courbe mathématique lisse qui épouse parfaitement le cercle. Résultat : une précision parfaite avec très peu de pièces.
• Les fissures : Si une fissure traverse la pièce, l'ordinateur place simplement un carré autour de la pointe de la fissure. La méthode mathématique calcule automatiquement la "catastrophe" (la concentration de contrainte) à ce point précis, sans avoir besoin de raffiner le maillage autour.

5. Pourquoi c'est génial pour tout le monde ?

Imaginez que vous devez analyser un réacteur nucléaire fissuré ou une pièce de voiture avec des trous complexes.

• Avant : Un ingénieur passait des jours à dessiner manuellement le maillage, à ajuster les tailles des briques, à corriger les erreurs de connexion.
• Maintenant (avec cette méthode) : L'ingénieur donne juste la forme de l'objet et quelques points de contrôle. L'ordinateur génère le maillage en quelques secondes (moins d'une seconde pour certains exemples !).
• Précision : Les résultats sont aussi précis que les méthodes manuelles, mais beaucoup plus rapides.

En résumé :
Les auteurs ont créé un "robot maçon" qui sait construire des murs (le maillage) automatiquement, même sur des terrains accidentés. Grâce à une astuce mathématique (la méthode SBFEM), ce robot n'a pas besoin de poser des briques partout ; il suffit qu'il connaisse les contours pour prédire exactement comment la structure va réagir aux forces, même si elle se fissure. C'est plus rapide, plus simple et tout aussi précis.

Résumé technique

1. Problématique et Contexte

L'analyse par éléments finis (FEA) est un outil central en ingénierie, mais sa précision dépend fortement de la qualité du maillage. La génération de maillages de haute qualité pour des géométries complexes, comportant des frontières courbes ou des singularités (comme les fissures), reste une étape chronophage et difficile à automatiser.

Les méthodes traditionnelles basées sur les arbres de quadtree (maillages hiérarchiques carrés) souffrent de deux problèmes majeurs :

1. Les nœuds flottants (hanging nodes) : Lorsqu'une cellule est subdivisée, des nœuds apparaissent sur les arêtes des cellules adjacentes non subdivisées. Cela crée des incompatibilités de déplacement qui nécessitent des traitements spéciaux (sous-triangulation, fonctions de forme spéciales ou contraintes).
2. Approximation des frontières courbes : Les cellules quadtree étant carrées, elles ne suivent pas naturellement les frontières courbes. Une subdivision excessive (raffinement local) est souvent requise pour approximer la géométrie, ce qui augmente le nombre de degrés de liberté (DDL) et la complexité du calcul.

L'objectif de cet article est de proposer une technique automatisée pour l'analyse des contraintes et de la rupture qui surmonte ces limitations en combinant la Méthode des Éléments Finis à Frontière Mise à l'Échelle (SBFEM) avec des maillages quadtree utilisant des éléments d'ordre élevé.

2. Méthodologie

La méthode proposée intègre la SBFEM dans le cadre des cellules quadtree, transformant chaque cellule (carrée ou polygonale après découpe) en un polygone SBFEM.

A. Traitement des Nœuds Flottants et des Cellules

Contrairement aux méthodes FEM classiques, la SBFEM traite chaque cellule quadtree comme un polygone à nombre arbitraire de côtés.

• Absence de sous-triangulation : Les nœuds flottants sont traités comme des nœuds ordinaires (sommets du polygone). Aucune subdivision de cellule n'est nécessaire pour assurer la compatibilité des déplacements.
• Formulation semi-analytique : La solution à l'intérieur de la cellule est exprimée semi-analytiquement en fonction de la coordonnée radiale $\xi$ (depuis le centre de mise à l'échelle jusqu'à la frontière) et numériquement le long de la frontière via des éléments linéaires 1D.
• Éléments d'ordre élevé : Il est possible d'utiliser directement des éléments d'ordre élevé sur les arêtes des cellules quadtree, permettant une précision accrue sans augmenter excessivement le nombre de cellules.

B. Génération Automatique du Maillage

L'algorithme de maillage est entièrement automatique et repose sur :

• Fonction de distance signée : La géométrie est définie par une fonction de distance signée, permettant de gérer facilement les géométries complexes et les opérations booléennes.
• Points de semence (Seed points) : L'utilisateur définit uniquement des points de semence sur les frontières et les zones d'intérêt pour contrôler la densité du maillage.
• Règle 2:1 : Le maillage respecte la règle 2:1 (différence maximale de niveau de subdivision entre cellules adjacentes), limitant le nombre de types de cellules maîtresses à 6 (ou 24 avec rotation), ce qui optimise le calcul.
• Adaptation aux frontières courbes : Les cellules intersectées par une frontière courbe sont découpées en polygones. Au lieu de raffiner excessivement (h-raffinement), la méthode utilise des éléments d'ordre élevé (p-raffinement) pour capturer précisément la courbure avec peu d'éléments.

C. Analyse de Rupture et Singularités

• Modélisation des fissures : Le sommet d'une fissure est choisi comme centre de mise à l'échelle d'une cellule SBFEM.
• Singularité analytique : La variation des contraintes (y compris la singularité) le long de la direction radiale est résolue analytiquement. Aucune enrichissement asymptotique (comme dans la XFEM) ni raffinement local extrême n'est nécessaire autour de la pointe de fissure.
• Facteurs d'intensité de contrainte (SIF) : Ils sont calculés directement à partir des modes propres de la matrice hamiltonienne de la SBFEM.

3. Contributions Clés

1. Automatisation complète : Élimination des étapes manuelles complexes de maillage (comme la définition de zones de raffinement autour des fissures ou la gestion des nœuds flottants).
2. Gestion native des nœuds flottants : La formulation SBFEM permet d'utiliser les cellules quadtree telles quelles, sans fonctions de forme de transition complexes ni subdivision triangulaire.
3. Précision avec peu de DDL : L'utilisation d'éléments d'ordre élevé permet de modéliser des frontières courbes et des champs de contraintes complexes avec un nombre de degrés de liberté nettement inférieur aux méthodes FEM classiques.
4. Efficacité computationnelle : Une stratégie de construction de la matrice de rigidité globale est proposée. Seules les cellules maîtresses (types de base) doivent être calculées une fois ; les autres cellules réutilisent ces matrices par rotation ou translation, réduisant drastiquement le temps de calcul pour les grands maillages.
5. Robustesse géométrique : Capacité à traiter des géométries complexes, des frontières courbes et des configurations de fissures multiples sans dégradation de la qualité du maillage.

4. Résultats Numériques

Cinq exemples numériques valident la technique :

1. Plaque infinie avec trou circulaire : Comparaison avec la solution analytique. La méthode montre une convergence monotone et rapide. L'utilisation d'éléments d'ordre 4 ou 5 donne une précision supérieure avec moins de DDL que les éléments linéaires.
2. Plaque carrée avec 9 trous : Démontre la flexibilité de l'automatisation pour gérer les transitions de maillage entre des caractéristiques géométriques de tailles variées. Les résultats concordent bien avec ceux d'ANSYS, mais avec un temps de génération de maillage très court (~3s) et une structure de maillage structurée.
3. Plaque avec fissures émanant d'un trou : Analyse de singularités multiples. La méthode capture les singularités sans raffinement local autour des pointes de fissure, validant l'efficacité pour les problèmes de rupture complexes.
4. Plaque avec deux fissures croisées : Mise en évidence de la simplicité de la génération de maillage pour des configurations de fissures complexes, comparé aux multiples étapes manuelles requises en FEM classique (ex: commande KSCON dans ANSYS).
5. Réacteur nucléaire fissuré sous pression : Application à une structure industrielle réelle avec géométrie non régulière. La méthode génère un maillage en moins d'une seconde et converge rapidement avec l'ordre des éléments.

Dans tous les cas, les résultats montrent une excellente concordance avec les solutions analytiques et les simulations FEM de référence, avec des taux de convergence élevés grâce à l'utilisation d'éléments d'ordre élevé.

5. Signification et Impact

Cette recherche représente une avancée significative dans l'automatisation de l'analyse par éléments finis. En combinant la flexibilité des maillages quadtree (facilité de génération) avec la puissance analytique de la SBFEM (gestion des singularités et des polygones), la méthode proposée :

• Réduit considérablement le temps de pré-traitement (génération de maillage).
• Élimine les erreurs humaines liées à la définition manuelle des zones de raffinement.
• Permet des analyses de fissuration et de propagation de fissures plus robustes et moins coûteuses en ressources computationnelles.
• Offre une alternative efficace aux méthodes d'enrichissement (XFEM) pour les problèmes de rupture, tout en conservant la simplicité d'un maillage structuré.

En conclusion, cette technique ouvre la voie à des outils de simulation CAE plus automatisés, capables de gérer des géométries complexes et des phénomènes physiques singuliers avec une efficacité et une précision accrues.