Modélisation double-échelle de la rupture des roches : influence du frottement sur les micro-fissures

Double-scale modelling of failure in rocks : influence of micro-cracks friction

• Auteur / Autrice : Aleksandra Wrzesniak
• Direction : Cristian Dascalu
• Type : Thèse de doctorat
• Discipline(s) : Matériaux, mécanique, génie civil, électrochimie
• Date : Soutenance le 14/12/2012
• Etablissement(s) : Grenoble
• Ecole(s) doctorale(s) : École doctorale Ingénierie - matériaux mécanique énergétique environnement procédés production (Grenoble ; 2008-....)
• Partenaire(s) de recherche : Laboratoire : Sols, solides, structures - risques (Grenoble)
• Jury : Président / Présidente : Irini Djeran-Maigre
• Jury : Examinateurs / Examinatrices : Cristian Dascalu, Pierre Bésuelle, Darius Seyedi
• Jury : Rapporteur / Rapporteuse : André Dragon, Jian-Fu Shao

Mots clés

Mots clés contrôlés

Géomatériaux

Frottement

Méthodes d'homogénéisation numérique

Mots clés libres

Micro-fissures

Homogénéisation asymptotique

Endommagement

Contact unilatéral

Frottement

Géomatériaux

Comportement quasi-fragile

Propagation sous-critique

Effet du temps

Micro-cracks

Asymptoptic homogenization

Damage

Unilateral contact

Friction

Geomaterials

Quasi-brittle behaviour

Sub-critic

Time effects

Résumé

Propagation des fissures microscopiques, est représentée par des variables d’endommagement. L’évolution de la variable d’endommagement est généralement formulée sur la base d’observations expérimentales. De nombreux modèles phénoménologiques d’endommagement ont été proposés dans la littérature. L’objet de cette thèse est de développer une nouvelle procédure pour obtenir des lois d’évolution macroscopique d’endommagement,dans lesquelles l’évolution de l’endommagement est entièrement déduite de l’analyse de la microstructure. Nous utilisons une homogénéisation basée sur des développements asymptotiques pour décrire le comportement global à partir de la description explicite d’un volume élémentaire microfissuré.Nous considérons d’une part un critère quasi-fragile (indépendant du temps) puis un critère sous-critique(dépendant du temps) pour décrire la propagation des microfissures. De plus, le frottement entre les lèvres des microfissures est pris en compte. Une analyse énergétique est proposée, conduisant à une loi d’évolution d’endommagement qui intègre une dégradation de la rigidité, un adoucissement du comportement du matériau, des effets de taille et d’unilatéralité, mettant en avant un comportement différent à la rupture en contact avec et sans frottement. L’information sur les micro-fissures est contenue dans les coefficients homogénéisés et dans la loi d’évolution de l’endommagement. Les coefficients homogénéisés décrivent la réponse globale en présence de micro-fissures (éventuellement statiques), tels qu’ils sont calculées avec la(quasi-) solution microscopique statique. La loi d’endommagement contient l’information sur l’évolution des micro-fissures, résultant de l’équilibre énergétique dans le temps pendant la propagation microscopique.La loi homogénéisée est formulée en incrément de contrainte. Les coefficients homogénéisés sont calculées numériquement pour des longueurs de fissures et des orientations différentes. Cela permet la construction complète des lois macroscopiques. Une première analyse concerne le comportement local macroscopique, pour des trajets de chargement complexes, afin de comprendre le comportement prédit par le modèle à deux échelles et l’influence des paramètres micro structuraux, comme par exemple le coefficient de frottement. Ensuite, la mise en œuvre en éléments finis des équations macroscopiques est effectuée et des simulations pour différents essais de compression sont réalisées. Les résultats des simulations numériques sont comparés avec les résultats expérimentaux obtenus en utilisant un nouvel appareil triaxial récemment mis au point au Laboratoire 3SR à Grenoble (France).

In continuum damage models, the degradation of the elastic moduli, as the results of microscopic crackgrowth, is represented through damage variables. The evolution of damage variable is generally postulatedbased on the results of the experimental observations. Many such phenomenological damage modelshave been proposed in the literature. The purpose of this contribution is to develop a new procedurein order to obtain macroscopic damage evolution laws, in which the damage evolution is completelydeduced from micro-structural analysis. We use homogenization based on two-scale asymptotic developmentsto describe the overall behaviour starting from explicit description of elementary volumes withmicro-cracks. We consider quasi-brittle (time independent) and sub-critical (time dependent) criteria formicro-cracks propagation. Additionally, frictional contact is assumed on the crack faces. An appropriatemicro-mechanical energy analysis is proposed, leading to a damage evolution law that incorporates stiffnessdegradation, material softening, size effect, and unilaterality, different fracture behaviour in contactwithout and with friction. The information about micro-cracks is contained in the homogenized coefficientsand in the damage evolution law. The homogenized coefficients describe the overall response inthe presence of (possibly static) micro-cracks, as they are computed with the (quasi-) static microscopicsolution. The damage law contains the information about the evolution of micro-cracks, as a result ofthe energy balance in time during the microscopic propagation. The homogenized law is obtained in therate form. Effective coefficients are numerically computed for different crack lengths and orientations.This allows for the complete construction of the macroscopic laws. A first analysis concerns the localmacroscopic behaviour, for complex loading paths, in order to understand the behaviour predicted bythe two-scale model and the influence of micro structural parameters, like for example friction coefficient.Next, the FEM implementation of the macroscopic equations is performed and simulations for variouscompression tests are conducted. The results of the numerical simulations are compared with the experimentalresults obtained using a new true-triaxial apparatus recently developed at the Laboratory 3SRin Grenoble (France).