The semi-explicit nonsmooth Newmark time integrator for robust unilateral contact in dynamic fragmentation simulations

Cet article introduit et valide un schéma d'intégration temporelle de type Newmark semi-explicite non lisse (NSN) qui gère de manière robuste le contact unilatéral dans les simulations de fragmentation dynamique en imposant strictement les contraintes, atteignant ainsi une stabilité et une précision supérieures aux méthodes basées sur la pénalité tout en révélant que la dissipation de contact peut paradoxalement augmenter le nombre de fragments en améliorant la localisation des dommages.

L'essentiel

Imaginez que vous essayez de simuler ce qui se passe lorsqu'un objet solide, comme une assiette en céramique ou un rocher, est frappé si fort qu'il vole en éclats en des milliers de petits morceaux. Il ne s'agit pas d'une simple cassure ; c'est une explosion chaotique où les morceaux s'envolent, s'entrechoquent, rebondissent sur les murs et se broient les uns les autres.

Le document présente un nouveau programme informatique (un « intégrateur temporel ») conçu pour simuler ce chaos sans que l'ordinateur ne plante ou ne donne des résultats absurdes. Voici comment cela fonctionne, décomposé en concepts simples :

1. Le problème : Le piège du « ressort »

Pour simuler la rupture, les scientifiques utilisent généralement une méthode où ils font semblant que le matériau est composé de minuscules ressorts. Lorsque le matériau se brise, les ressorts se rompent. Lorsque les morceaux s'entrechoquent, ils utilisent des « ressorts de pénalité » pour les repousser afin qu'ils ne se traversent pas les uns les autres.

L'analogie : Imaginez essayer d'arrêter une boule de bowling avec un élastique.

• Si l'élastique est trop lâche (faible rigidité), la boule passe directement à travers (non physique).
• Si l'élastique est super tendu (haute rigidité) pour arrêter la boule parfaitement, il agit comme un mur rigide. Mais si vous le rendez trop tendu, l'ordinateur doit effectuer des étapes de calcul minuscules, minuscules pour calculer le rebond, ce qui rend la simulation interminable.
• La prétention du papier : L'ancienne méthode (utilisant ces ressorts serrés) est instable. Elle provoque une dérive de l'ordinateur, une perte d'énergie ou un plantage, surtout lorsqu'il y a des millions de collisions simultanées.

2. La solution : Le « agent de circulation » (Newmark non lisse)

Les auteurs ont créé une nouvelle méthode appelée Newmark Non Lisse (NSN - Nonsmooth Newmark). Au lieu d'utiliser des élastiques pour repousser les pièces, cette méthode agit comme un agent de circulation strict à une intersection très fréquentée.

L'analogie :

• La masse (La voiture) : Le corps principal de l'objet se déplace librement et de manière fluide. L'ordinateur prédit où la voiture irait si elle n'y avait pas d'obstacles. Cette partie est calculée très rapidement (explicitement).
• Le contact (L'intersection) : Si la voiture percute un mur ou une autre voiture, l'agent de circulation intervient. Au lieu de repousser la voiture avec un ressort, l'agent dit instantanément : « Stop ! Vous ne pouvez pas aller là. » Il impose une règle stricte : Pas de traversée.
• La magie : Cette méthode traite la règle du « ne pas passer à travers » comme une loi dure de la physique plutôt que comme un ressort mou. Cela permet à l'ordinateur de prendre des pas de temps beaucoup plus grands car il n'a pas à se soucier de l'élastique qui deviendrait trop tendu.

3. L'approche à la « double personnalité »

Le papier décrit cette méthode comme étant « semi-explicite ». Voyez cela comme une danse en deux étapes :

1. Étape A (La prédiction) : L'ordinateur devine où tout se trouvera au moment suivant, en ignorant les collisions.
2. Étape B (La correction) : Si la prédiction montre que deux pièces se chevauchent, l'ordinateur corrige instantanément leur vitesse et leur position pour réparer le chevauchement, tout comme une bille de billard qui frappe une autre et change instantanément de direction.

Cela permet à la simulation d'être rapide (comme la prédiction) mais précise et stable (comme la correction).

4. Ce qu'ils ont trouvé (Les expériences)

Les auteurs ont testé cette nouvelle méthode d'« agent de circulation » contre les anciennes méthodes de « l'élastique » en utilisant trois scénarios :

• La balle rebondissante : Une balle simple rebondissant sur le sol. La nouvelle méthode était aussi précise que les meilleures méthodes existantes, mais gérait les rebonds sans perdre d'énergie ou devenir instable.
• La barre percutante : Une barre métallique frappant un mur. Les anciennes méthodes avaient du mal avec la vitesse de l'impact, mais la nouvelle méthode a géré le « broyage » parfaitement, en gardant les calculs d'énergie corrects.
• La barre qui vole en éclats : Une barre déjà fissurée qui frappe un mur. Les anciennes méthodes nécessitaient des pas de temps si minuscules pour rester stables qu'elles étaient incroyablement lentes. La nouvelle méthode pouvait prendre des pas énormes, courant 27 fois plus vite tout en étant plus précise.

5. La découverte surprenante : Le fractionnement confiné

La partie la plus intéressante du papier concerne une expérience « confinée ». Imaginez une barre qui vole en éclats à l'intérieur d'une petite boîte plutôt que dans un espace ouvert.

• L'ancienne intuition : Vous pourriez penser que si les morceaux rebondissent sur les murs et perdent de l'énergie (dissipation), il y aura moins d'énergie restante pour briser le matériau, ce qui résulterait en moins de morceaux, mais plus gros.
• La conclusion du papier : C'est l'inverse qui s'est produit. Lorsque les morceaux rebondissaient sur les murs et perdaient un peu d'énergie (dissipation de contact), le matériau se brisait en plus de morceaux, plus petits.
• Pourquoi ? Les auteurs expliquent que le « rebond » agit comme un filtre. Dans un monde parfaitement élastique (rebondissant), les ondes de contrainte rebondissent de manière sauvage, provoquant une « confusion » du matériau qui développe de nombreuses petites fissures fragiles qui ne se séparent pas complètement. Lorsque les murs absorbent une partie de cette énergie, les ondes se calment. Cela permet à la contrainte de se concentrer sur des points spécifiques, poussant les fissures à travers tout le matériau pour créer des fragments nets et séparés.

Résumé

Le papier présente un nouvel outil mathématique pour simuler la rupture d'objets en traitant les collisions comme des règles dures et instantanées plutôt que comme des ressorts mous. Cela rend la simulation informatique :

1. Plus stable : Elle ne plante pas et ne dérive pas.
2. Plus rapide : Elle peut prendre des pas de temps plus grands.
3. Plus précise : Elle prédit correctement en combien de morceaux un objet va se briser.

Les auteurs concluent que cet outil est prêt à être utilisé pour des simulations 3D complexes, telles que comprendre comment les débris spatiaux se fragmentent ou comment les roches se brisent lors d'avalanches, en fournant un moyen robuste de gérer la danse chaotique de millions de fragments en collision.

Résumé technique

Résumé Technique : L'intégrateur temporel de Newmark semi-explicite non lisse pour une gestion robuste du contact unilatéral dans les simulations de fragmentation dynamique

Énoncé du problème
Les simulations numériques de fragmentation dynamique présentent un double défi : capturer l'évolution topologique rapide d'un continuum solide se brisant en fragments discrets et résoudre les contacts denses et simultanés qui se produisent entre ces fragments et les faces des fissures. Bien que la méthode des éléments finis (FEM) avec des modèles de zone cohésive (CZM) extrinsèques traite efficacement la fracture, les traitements de contact conventionnels posent des difficultés significatives. Les méthodes basées sur la pénalité, souvent utilisées pour leur compatibilité avec l'intégration temporelle explicite, souffrent d'un compromis fondamental : une faible rigidité permet une interpénétration non physique, tandis qu'une rigidité élevée restreint le pas de temps stable à des niveaux de calcul prohibitifs. De plus, l'interaction entre la fracture et le contact régularisé introduit une non-linéarité (non-lissage) qui entraîne une dérive d'énergie numérique persistante et une instabilité, masquant ainsi la physique sous-jacente.

Méthodologie
Les auteurs proposent et valident un nouveau schéma d'intégration temporelle semi-explicite : la méthode de Newmark non lisse (NSN). Cette approche adapte le cadre de la Dynamique de Contact Non Lisse (NSCD) à la fragmentation dynamique en découplant le traitement de la dynamique de masse et des interactions de contact :

1. Formulation semi-explicite : Le schéma sépare la dynamique en une partie « lisse » (déformation de masse, chargement externe et fracture cohésive) et une partie « non lisse » (réactions de contact).
   • La dynamique lisse est intégrée explicitement à l'aide d'un schéma de Newmark-$\beta$ de second ordre ($\beta=0, \gamma=1/2$), maintenant une grande précision lors des phases de vol libre.
   • La dynamique non lisse (contact unilatéral) est traitée implicitement au niveau de la vitesse à l'aide d'une approche de Euler implicite de premier ordre. Cela impose rigoureusement les contraintes de vide de Signorini sans compliance artificielle.
2. Mise en œuvre algorithmique : Le problème de contact est formulé comme un Problème de Complémentarité Linéaire (LCP), qui est résolu comme un problème de Programmation Quadratique (QP) convexe restreint à l'ensemble de contact actif. Cela permet des solutions efficaces et parallélisables qui évoluent selon le nombre de contacts actifs plutôt que selon la taille totale du système.
3. Loi de traction-séparation (TSL) modifiée : Pour répondre au défi numérique de la rigidité cohésive infinie à mesure que l'endommagement approche de zéro dans les schémas explicites, les auteurs introduisent un plafonnement de la rigidité ($\tilde{k}$). Cette régularisation empêche le pas de temps de tendre vers zéro tout en préservant la convexité du problème de contact.
4. Validation et application : Le cadre est implémenté dans la bibliothèque FEM open-source Akantu. Il est validé par rapport à des solutions analytiques et des méthodes non lisses établies (Moreau-Jean et CD-Lagrange) en utilisant des tests de référence tels que la balle rebondissante et la barre d'impact. Il est ensuite appliqué à la fragmentation dynamique 1D sous expansion libre et confinée.

Contributions clés et résultats

• Précision et stabilité : Les tests de référence démontrent que le schéma NSN atteint une précision comparable aux méthodes non lisses établies (Moreau-Jean et CD-Lagrange) et surpasse de manière significative les approches par pénalité. Dans le test de la barre d'impact, le NSN maintient une précision de second ordre pour le déplacement et une précision de vitesse supérieure à celle de CD-Lagrange, particulièrement dans les régimes élastiques.
• Efficacité computationnelle : Bien que la méthode NSN incurre un coût de calcul par pas de temps plus élevé (environ 5 fois celui d'un schéma de pénalité purement explicite en raison de la résolution du QP), sa stabilité robuste permet des pas de temps nettement plus grands. Pour les benchmarks 1D, cela se traduit par un temps de résolution total comparable ou meilleur que les approches purement explicites, avec des erreurs de conservation d'énergie réduites de plusieurs ordres de grandeur (proche de la précision machine).
• Perspectives physiques dans la fragmentation confinée :
  • Décalage du budget énergétique : Dans les scénarios d'expansion confinée (où les fragments entrent en collision avec des parois rigides), le budget énergétique pour la fracture passe de l'énergie cinétique locale des fragments individuels à l'énergie cinétique globale de l'ensemble du système. Cela conduit à une réduction drastique de la taille moyenne des fragments et à une augmentation de l'énergie de fracture totale par rapport à l'expansion libre.
  • Rôle contre-intuitif de la dissipation de contact : L'étude révèle que l'introduction de la dissipation de contact (collisions inélastiques, $e < 1$) réduit l'énergie totale dissipée par la fracture, mais augmente le nombre final de fragments. Cela se produit car la dissipation de contact agit comme un filtre passe-bas, supprimant les oscillations d'ondes de contrainte à haute fréquence. Ce filtrage favorise une propagation d'ondes de contrainte plus propre et une meilleure localisation des dommages, poussant les fissures vers une séparation complète plutôt que de laisser le dommage rester distribué dans des zones partiellement endommagées ($d < 1$).

Signification
Cet article établit le schéma NSN comme un outil robuste pour générer des statistiques de fragmentation de haute fidélité dans des scénarios de contacts denses. En traitant rigoureusement le contact comme un phénomène non lisse tout en conservant l'efficacité explicite pour la dynamique de masse, la méthode surmonte les instabilités numériques inhérentes aux formulations par pénalité. Les auteurs affirment que ce cadre est essentiel pour simuler des interactions multi-corps complexes, telles que celles rencontrées dans les scénarios de débris orbitaux, où les interactions fréquentes entre fragments entraînent une fragmentation secondaire. Ce travail constitue une étape critique vers l'extension de l'analyse de fragmentation de haute fidélité à des dimensions supérieures, en tirant parti de la capacité de la méthode à préserver la fidélité physique sans le coût prohibitif des résolutions implicites globales.