Ce travail est dédié au traitement numérique de la fissuration au sein d'une structure dont la rupture résulte d'un endommagement ductile. Trois points importants sont ainsi à considérer : la modélisation du comportement du matériau, la modélisation de la fissure au sein de la structure, la méthode de propagation couplant les deux modélisations précédentes. L'étude est conduite en 2D en déformations planes (DP), dans le cadre des petites perturbations. Une procédure dynamique a été choisie avec un schéma de résolution implicite, en utilisant le code de calculs industriel Abaqus. Le comportement du matériau est modélisé via un modèle de type GTN modifié (GTNm), voir Longère et al, Int. J. Dam. Mech. (2012), prenant en compte les effets combinés de l'écrouissage, de l'adoucissement thermique, de la viscoplasticité et de l'endommagement ductile. Le modèle GTNm a été implanté en 2D DP dans Abaqus via une routine utilisateur (umat). Pour décrire les conséquences cinématiques de la présence d'une fissure au sein de la structure ainsi que de sa propagation, nous avons retenu la méthode des éléments finis étendus (X-FEM), voir Moes et al, Int. J. Numer Meth. Engng (1999). La méthode X-FEM a été implantée en 2D DP dans Abaqus via une routine utilisateur (uel). Afin de permettre le couplage entre la X-FEM et le comportement fortement non linéaire du matériau, certains choix ont dû être faits, notamment en ce qui concerne la forme du champ de déplacement enrichi ainsi que l'intégration numérique. Comme dans le travail de Haboussa et co-auteurs, voir Haboussa et al, Int. J. Numer Meth. Engng (2011), certaines grandeurs sont moyennées dans une pastille située en pointe de fissure dans le but d'atténuer la dépendance des résultats numériques au maillage. L'originalité de ce travail réside sans doute dans la méthode de propagation, dont le but est de reproduire une rupture induite par endommagement ductile. En premier lieu, la direction de propagation potentielle de la fissure est calculée via une analyse de bifurcation dans le cas d'un état matériau figé. Le principe est en effet de considérer que la direction de propagation de la fissure est la direction dans laquelle la déformation s'est localisée, comme c'est le cas dans la phase de coalescence Dans un second temps, un critère d'amorçage basé sur une énergie stockée en pointe de fissure a été défini. Si cette énergie dépasse une valeur critique, nous considérons que la fissure se propage. Nous calculons alors la longueur du nouvel incrément de fissure via une des deux méthodes retenues dans notre étude : une méthode dite par épuisement (MEPUI) et une méthode par contrôle de la vitesse de fissuration (MCVF). La MEPUI est basée sur le principe que la fissure se propage tant que le critère d'amorçage de la fissure est vérifié et est valable quelle que soit la vitesse de chargement de la structure. A contrario, la MCVF est basée sur la loi de Kanninen qui permet de contrôler la vitesse de fissuration et n'est valable que pour des vitesses de chargement élevées. Ces deux méthodes de propagation (MEPUI et MCVF) ont été implantées dans Abaqus via une routine utilisateur (uel). Afin de tester les deux méthodes mises en place dans notre étude (MEPUI et MCVF), des simulations numériques sur des plaques fissurées et entaillées ont été conduites. Considérant quelques simplifications, ce travail permet de reproduire la propagation d'une fissure dans une structure 2D dont la rupture résulte d'un endommagement par germination et croissance de cavités.
Identifer | oai:union.ndltd.org:CCSD/oai:tel.archives-ouvertes.fr:tel-01070096 |
Date | 10 December 2013 |
Creators | Crété, Jean-Philippe |
Publisher | Université de Bretagne Sud |
Source Sets | CCSD theses-EN-ligne, France |
Language | fra |
Detected Language | French |
Type | PhD thesis |
Page generated in 0.0015 seconds