FEMxDEM double scale approach with second gradient regularization applied to granular materials modelization

Approche double échelle de type FEMxDEM avec régularisation second gradient pour la modélisation des géomatériaux L'approche multi-échelle FEMxDEM est une méthode numérique innovante pour les problèmes géotechniques impliquant des matériaux granulaires. La méthode des éléments finis (FEM) et la méthode des éléments discrets (DEM) sont simultanément appliquées à résoudre, respectivement, le problème structurel à la macro-échelle et la microstructure du matériau à la micro-échelle. L'avantage d'utiliser une telle configuration à double échelle est de permettre d'étudier un problème d'ingénierie sans la nécessité de lois de comportement standard, capturant ainsi l'essence des propriétés des matériaux. Le lien entre les échelles est obtenu par homogénéisation numérique, de sorte que la loi de comportement continu numérique et la matrice tangente correspondante sont obtenues directement à partir de la réponse discrète de la microstructure.En règle générale, l'approche FEMxDEM présente quelques inconvénients; la vitesse de convergence et la robustesse de la méthode ne sont pas aussi efficaces que dans les modèles FEM classiques. De plus, le coût de calcul de l'intégration de la micro-échelle et la dépendance du maillage typique de la macro-échelle, rendent l'approche multi-échelle FEMxDEM discutable pour des utilisations pratiques. Le but de ce travail est de se concentrer sur ces questions théoriques et numériques avec l'objectif de rendre l'approche multi-échelle FEMxDEM robuste et applicable à des configurations à l'échelle réelle. Une variété d'opérateurs est proposée afin d'améliorer la convergence et la solidité de la méthode dans un cadre quasi-Newton. L'indépendance de l’intégration des différents points de Gauss et les caractéristiques d’intensivité sur les d'éléments sont exploités par l'utilisation d’une parallélisation en utilisant un paradigme OpenMP. Au niveau macro, une relation constitutive second gradient est mise en œuvre afin d'enrichir la relation de Cauchy de premier gradient apportant indépendance du maillage au modèle.Les améliorations susmentionnées rendent l'approche FEMxDEM compétitive avec les modèles FEM classiques en termes de coût de calcul permettant ainsi d'effectuer des simulations multi-échelle FEMxDEM robustes et indépendantes du maillage, depuis l'échelle du laboratoire (par exemple essaie biaxiale test) jusqu’à celle du problème à l'échelle de l'ingénierie (par exemple, excavation d’une galerie).Mots clés:Double échelle, homogénéisation numérique, loi constitutive numérique, élasto-plasticité, second gradient, matériaux microstructurés, grande déformation, éléments finis, éléments discrets, méthode de Newton, parallélisation, unicité.