Propagation robuste de défauts en 3D

Robust 3D crack propagation

• Auteur / Autrice : Matthieu Le Cren
• Direction : Nicolas Moës, Patrick Massin, Alexandre Martin
• Type : Thèse de doctorat
• Discipline(s) : Mécanique des solides, des matériaux, des structures et des surfaces
• Date : Soutenance le 18/10/2018
• Etablissement(s) : Ecole centrale de Nantes
• Ecole(s) doctorale(s) : Sciences de l'ingénierie et des systèmes (Nantes Université)
• Partenaire(s) de recherche : Laboratoire : Institut de Recherche en Génie Civil et Mécanique (Nantes)
• Jury : Président / Présidente : Anthony Gravouil
• Jury : Examinateurs / Examinatrices : Nicolas Moës, Patrick Massin, Alexandre Martin, Anthony Gravouil, Sylvie Pommier, Stéphane Pierre Alain Bordas, Claude Stolz
• Jury : Rapporteur / Rapporteuse : Sylvie Pommier, Stéphane Pierre Alain Bordas

Mots clés

Mots clés contrôlés

Modélisation tridimensionnelle

Matériaux -- Fissuration

Rupture, Mécanique de la

Mots clés libres

X-FEM

Rupture fragile

Taux de restitution d'énergie

Facteurs d'intensité de contrainte

X-FEM

Brittle fracture

Energy release rate

Stress intensity factors

Résumé

Afin d'assurer le contrôle de son parc de production d'électricité, EDF doit maîtriser le vieillissement de ses installations pour en garantir le bon fonctionnement dans la durée. Dans ce but, il est nécessaire de disposer d’outils performants pour le modéliser et simuler la propagation des défauts dans les structures.Dans ces travaux de thèse, on s’intéresse à la propagation de fissures avec la méthode X-FEM et notamment à l’étape de localisation de la fissure par une technique de courbes de niveau. Nous avons proposé une approche fondée sur une méthode de propagation d’information de distance dite fast marching method pour rendre cette étape plus robuste. Elle est applicable à tous types de mailles,linéaires ou quadratiques.De plus, le calcul du taux de restitution d’énergie et des facteurs d’intensité de contrainte en pointe de fissure doit être suffisamment précis pour permettre de calculer la direction et l’avancée de la fissure. Dans ce but, nous avons proposé d’étudier une méthode d’intégrale de domaine pour laquelle on soulève plusieurs difficultés liées à la représentation de la fissure dans un espace tridimensionnel. Plusieurs améliorations sont proposées pour rendre les calculs plus précis et plus robustes.Dans le cas des fissures à front courbe, nous avons identifié les limites de l'utilisation des champs asymptotiques obtenus en pointe de fissure sous l'hypothèse des déformations planes comme champs auxiliaires d’une méthode d’intégrale d’interaction et nous avons proposé de nouveaux champs de déplacements auxiliaires qui prennent en compte la courbure du front de fissure. Toutes ces approches sont développées et validées dans le logiciel code_aster.

In order to ensure the control of its nuclear power plants, EDF must guarantee that they function effectively over the long term. For this purpose, it is necessary to have efficient tools tomodel and simulate crack propagation in structures. In this PhD work, we are interested in the propagation of cracks with the X-FEM method which allows using the same mesh as for a structure without default. We target especially the reconstruction of thelevel sets that characterize the position of the crack after propagation. We have proposed a fast marching method approach based on the propagation of distance information from the crack surface to the whole structure to make this step more robust in the X-FEM propagation process. It is applicable to all types of meshes, linear or quadratic. The calculation of information characteristic of thecrack status such as the energy release rate and the stress intensity factors must be accurate enough to obtain the direction and advance of the crack front ateach propagation step. For this purpose, we proposed to study a domain integral method, for which several difficulties related to the representation of the crackin a three-dimensional space are identified. Several improvements are proposed to make the calculations more accurate and more robust. In the case of curved cracks front, we have identified the limitations of using asymptotic fields obtained under the plane deformation hypothesis as auxiliary fields of an interaction integral method and we have proposed new auxiliary displacement fields that take into account the curvature of the crack front. All these methods are developed and validated with EDF software code_aster.