EDF souhaite conforter la justification de l'intégrité des cuves de réacteurs nucléaires à eau pressurisée, basées sur l'amorçage de fissure, par le concept d'arrêt de fissure. L'objectif de la thèse est de consolider les travaux de prévision de l'arrêt des fissures fragiles dans l'acier de cuve. Des essais d'arrêt de fissure isothermes réalisés entre -150°C et -50°C sur des éprouvettes CT et sur des anneaux en compression ont été analysés. Les analyses fractographiques montrent que la propagation et l'arrêt de fissure se produisent en clivage même si de la déchirure ductile précède la propagation instable de fissure. Un critère local d'arrêt de fissure fragile est appliqué lors de calculs par éléments finis en élastoviscoplasticité et en dynamique : la rupture se produit si la contrainte principale maximale est égale à une contrainte critique. Ce critère permet de bien prédire les vitesses de propagation ainsi que les longueurs de fissure à l'arrêt. La contrainte critique augmente avec la température en relation avec les signes de dissipation observés sur les surfaces de rupture. Une dépendance du critère à la géométrie est observée. Elle peut être liée aux hypothèses de calculs réalisés en bidimensionnel. L'étude de la dynamique de la structure montre que le phénomène d'arrêt de fissure est fortement lié à la dynamique globale de la structure : l'arrêt et la fermeture de la fissure interviennent simultanément.

EDF souhaite conforter la justification de l’intégrité des cuves de réacteurs nucléaires à eau pressurisée, basées sur l’amorçage de fissure, par le concept d’arrêt de fissure. L'objectif de la thèse est de consolider les travaux de prévision de l’arrêt des fissures fragiles dans l’acier de cuve. Des essais d'arrêt de fissure isothermes réalisés entre -150°C et -50°C sur des éprouvettes CT et sur des anneaux en compression ont été analysés. Les analyses fractographiques montrent que la propagation et l’arrêt de fissure se produisent en clivage même si de la déchirure ductile précède la propagation instable de fissure. Un critère local d’arrêt de fissure fragile est appliqué lors de calculs par éléments finis en élastoviscoplasticité et en dynamique : la rupture se produit si la contrainte principale maximale est égale à une contrainte critique. Ce critère permet de bien prédire les vitesses de propagation ainsi que les longueurs de fissure à l’arrêt. La contrainte critique augmente avec la température en relation avec les signes de dissipation observés sur les surfaces de rupture. Une dépendance du critère à la géométrie est observée. Elle peut être liée aux hypothèses de calculs réalisés en bidimensionnel. L’étude de la dynamique de la structure montre que le phénomène d’arrêt de fissure est fortement lié à la dynamique globale de la structure : l’arrêt et la fermeture de la fissure interviennent simultanément. / EDF wants to complete the assessment of reactor pressure vessels usually based on crack initiation concept by crack arrest concept. The work aims at improving the knowledge of cleavage crack arrest in reactor pressure vessel steel. For that purpose, isothermal crack arrest experiments were performed for temperatures ranging from –150°C up to –50°C on compact tensile specimens and on pre-cracked rings submitted to compressive loading. Fractographic observations revealed that the whole crack propagation and arrest occurs by cleavage even if ductile tearing occurs before initiation of the unstable crack propagation. A local cleavage crack arrest criterion is applied in finite element computations carried out in elastoviscoplasticity and in full dynamics: the crack propagates since the largest principal stress reaches a critical stress. The application of this criterion on the experiments leads to a good prediction of the crack speed and of the crack length and shows that the critical stress increases with the temperature in relation with dissipation features observed on the fracture surfaces. Dependence to the geometry is observed; it can be due to the assumption used for the 2D computations. The study of the structural dynamic shows that the crack arrest phenomenon is very linked to the global dynamics of the structure: crack arrest and crack closure occur approximately at the same time.

L’objectif de cette thèse est de comprendre les micro-mécanismes physiques de propagation et d’arrêt de fissure de clivage dans l’acier de cuve 16MND5 et de proposer un modèle de prédiction robuste et physiquement fondé, en s’appuyant sur une campagne d’essais de rupture fragile sur éprouvettes de laboratoire finement instrumentées, associée à la modélisation numérique de ces essais. Dans un premier temps, des expériences ont été menées sur des éprouvettes CT25 de différentes épaisseurs à cinq températures (-150°C, -125°C, -100°C, -75°C, -50°C). Des trajets de fissures rectilignes et branchées (deux fissures se développant de manière quasi-symétrique) ont été observés. Pour estimer la vitesse de propagation, une caméra ultra-rapide a été utilisée, associée à la mise au point d’un protocole expérimental permettant d’observer la face de l’éprouvette dans l’enceinte thermique, sans givrage. Des observations à 500 000 images.s-1 ont permis de caractériser finement la vitesse instantanée de la fissure sur le ligament complet de la CT (~25 mm). En parallèle, pour pouvoir analyser les essais et l’impact de la viscosité sur la réponse mécanique autour de la fissure, le comportement élasto-viscoplastique du matériau a été étudié jusqu’à une vitesse de déformation de 104 s-1 pour les températures étudiées. La méthode des éléments finis étendus (X-FEM) a été utilisée dans le code de calcul CAST3M pour modéliser la propagation de fissure. Les simulations numériques associent l’approche locale de la rupture en dynamique non linéaire et un critère de propagation en contrainte critique de type RKR à une distance caractéristique. Les travaux réalisés ont permis de confirmer la forme du critère proposé par Prabel à -125°C, et d’identifier les dépendances de ce critère à la température et à la vitesse de déformation. A partir d’analyses numériques en 2D et 3D, un critère multi-température fonction croissante de la vitesse de déformation est proposé. Des modélisations prédictives ont permis de valider le critère sur deux géométries d’éprouvettes (CT et anneau) en mode I à différentes températures. Des observations MEB et des analyses 3D au microscope optique montrent que le mécanisme de rupture est le clivage associé à des zones de cisaillement ductile entre les différents plans de fissuration. L’étude de la fraction surfacique des marches de cisaillement et des contraintes de fermeture associées tend à justifier le critère mis en place. Un modèle analytique est proposé permettant de justifier le critère déduit des modélisations numériques. Ce modèle considère que les ligaments retiennent la lèvre de la fissure et induisent donc des contraintes de fermeture au niveau de la pointe de fissure qu’il faut compenser pour atteindre la contrainte de clivage effective en pointe de fissure. Cette résistance des ligaments est directement reliée à la loi de comportement du matériau et justifie la dépendance du critère de rupture identifié à la vitesse de déformation. Enfin, les branchements de fissure ont été analysés via le dépouillement des vidéos obtenues avec la caméra rapide qui mettent en évidence un amorçage initial rectiligne, puis un amorçage de fissures multiples de part et d’autre du plan de fissure qui conduisent à l’arrêt de la fissure initiale, l’une de ces fissures ‘secondaires’ conduisant ensuite à la rupture de l’éprouvette. Les rôles essentiels de l’épaisseur et du chargement dans ce mécanisme de branchement sont soulignés. L’augmentation de l’épaisseur réduit la fréquence d’apparition de ce mécanisme et finit même par l’annuler. Logiquement l’intensité du chargement doit être suffisamment importante pour créer cette zone plastique étendue : les essais qui présentent une propagation rectiligne sont les essais pour lesquels les chargements à l’amorçage sont les plus faibles. / The purpose of this PhD thesis is to understand physical mechanisms of cleavage crack propagation and arrest in the 16MND5 PWR vessel steel and to propose a robust predicting model based on a brittle fracture experimental campaign of finely instrumented laboratory specimens associated with numerical computations. First, experiments were carried out on thin CT25 specimens at five temperatures (-150°C, -125°C, - 100°C, -75°C, -50°C). Two kinds of crack path, straight or branching path, have been observed. To characterize crack propagation and to measure crack speed, a high-speed framing camera system was used, combined with the development of an experimental protocol which allowed to observe CT surface without icing inside the thermal chamber and on the specimen. The framing camera (520 000 fps) has allowed to have a very accurate estimation of crack speed on the complete ligament of CT (~ 25 mm). Besides, to analyse experiments and to study the impact of viscosity on the mechanical response around the crack tip, the elastic-viscoplastic behavior of the ferritic steel has been studied up to a strain rate of 104 s-1 for the tested temperatures.The eXtended Finite Element Method (X-FEM) was used in CAST3M FE software to model crack propagation. Numerical computations combine a local non linear dynamic approach with a RKR type fracture stress criterion to a characteristic distance. The work carried out has confirmed the form of the criterion proposed by Prabel at -125°C, and has identified the dependencies of the criterion on temperature and strain rate. From numerical analyzes in 2D and 3D, a multi-temperature fracture stress criterion, increasing function of the strain rate, was proposed. Predictive modeling were used to confirm the identified criterion on two specimen geometries (CT and compressive ring) in mode I at different temperatures. SEM observations and 3D analyzes made with optical microscope showed that the fracture mechanism was the cleavage associated with ductile shear zones between the different planes of cracking. The study of the surface fraction of ductile shear zones and associated closing stress tends to justify the established criterion. An analytical model is proposed to justify the criterion deduced from numerical modeling. This model assumes that the ligaments hold the crack lips and therefore induce closing stress at the crack tip which has to be compensated to achieve the effective cleavage stress at the crack tip. This resistance of ligaments is directly related to the mechanical behavior of the material and justifies the dependence of fracture stress criterion with strain rate. Lastly, the crack branching was analyzed thanks to movies obtained with the high-speed camera. The cleavage crack propagates in a straight way over few millimeters. Then, new cracks appear on both sides of the initial crack lips which leads to the arrest of the initial crack. One of the new cracks leads to the failure of the CT specimen. The essential role of the thickness and loading on this branching mechanism is emphasized. The increased thickness reduces the frequency of occurrence of this mechanism and eventually even cancel. Low thicknesses lead to more extensive plasticity at the crack tip and generates the germs for the triggering of secondary cracks that appear. Logically, the intensity of loading must be large enough to create this extended plastic zone : the experiments with a straight path are the tests for which the initial loadings are the lowest.

L'objectif de cette thèse est de comprendre les micro-mécanismes physiques de propagation et d'arrêt de fissure de clivage dans l'acier de cuve 16MND5 et de proposer un modèle de prédiction robuste et physiquement fondé, en s'appuyant sur une campagne d'essais de rupture fragile sur éprouvettes de laboratoire finement instrumentées, associée à la modélisation numérique de ces essais. Dans un premier temps, des expériences ont été menées sur des éprouvettes CT25 de différentes épaisseurs à cinq températures (-150°C, -125°C, -100°C, -75°C, -50°C). Des trajets de fissures rectilignes et branchées (deux fissures se développant de manière quasi-symétrique) ont été observés. Pour estimer la vitesse de propagation, une caméra ultra-rapide a été utilisée, associée à la mise au point d'un protocole expérimental permettant d'observer la face de l'éprouvette dans l'enceinte thermique, sans givrage. Des observations à 500 000 images.s-1 ont permis de caractériser finement la vitesse instantanée de la fissure sur le ligament complet de la CT (~25 mm). En parallèle, pour pouvoir analyser les essais et l'impact de la viscosité sur la réponse mécanique autour de la fissure, le comportement élasto-viscoplastique du matériau a été étudié jusqu'à une vitesse de déformation de 104 s-1 pour les températures étudiées. La méthode des éléments finis étendus (X-FEM) a été utilisée dans le code de calcul CAST3M pour modéliser la propagation de fissure. Les simulations numériques associent l'approche locale de la rupture en dynamique non linéaire et un critère de propagation en contrainte critique de type RKR à une distance caractéristique. Les travaux réalisés ont permis de confirmer la forme du critère proposé par Prabel à -125°C, et d'identifier les dépendances de ce critère à la température et à la vitesse de déformation. A partir d'analyses numériques en 2D et 3D, un critère multi-température fonction croissante de la vitesse de déformation est proposé. Des modélisations prédictives ont permis de valider le critère sur deux géométries d'éprouvettes (CT et anneau) en mode I à différentes températures. Des observations MEB et des analyses 3D au microscope optique montrent que le mécanisme de rupture est le clivage associé à des zones de cisaillement ductile entre les différents plans de fissuration. L'étude de la fraction surfacique des marches de cisaillement et des contraintes de fermeture associées tend à justifier le critère mis en place. Un modèle analytique est proposé permettant de justifier le critère déduit des modélisations numériques. Ce modèle considère que les ligaments retiennent la lèvre de la fissure et induisent donc des contraintes de fermeture au niveau de la pointe de fissure qu'il faut compenser pour atteindre la contrainte de clivage effective en pointe de fissure. Cette résistance des ligaments est directement reliée à la loi de comportement du matériau et justifie la dépendance du critère de rupture identifié à la vitesse de déformation. Enfin, les branchements de fissure ont été analysés via le dépouillement des vidéos obtenues avec la caméra rapide qui mettent en évidence un amorçage initial rectiligne, puis un amorçage de fissures multiples de part et d'autre du plan de fissure qui conduisent à l'arrêt de la fissure initiale, l'une de ces fissures 'secondaires' conduisant ensuite à la rupture de l'éprouvette. Les rôles essentiels de l'épaisseur et du chargement dans ce mécanisme de branchement sont soulignés. L'augmentation de l'épaisseur réduit la fréquence d'apparition de ce mécanisme et finit même par l'annuler. Logiquement l'intensité du chargement doit être suffisamment importante pour créer cette zone plastique étendue : les essais qui présentent une propagation rectiligne sont les essais pour lesquels les chargements à l'amorçage sont les plus faibles.

Ce mémoire de thèse présente l'étude de la propagation et de l'arrêt d'une fissure de clivage dans un acier de cuve REP. Une bonne compréhension des phénomènes en jeu nécessite de solides données expérimentales, ainsi qu'un outil de modélisation performant. La méthode X-FEM est implantée dans Cast3m, permettant de simuler efficacement la propagation de fissure par éléments finis. Deux techniques concernant l'actualisation des fonctions de niveau, ainsi qu'une intégration non conforme sont proposées. La campagne d'essais de propagation de fissure concerne trois géométries : des éprouvettes CT, et des anneaux comprimés en mode I et mixte. On relève les vitesses de propagation. Des fractographies désignent le clivage comme responsable de la ruine. Un modèle de propagation basé sur la contrainte principale en pointe de fissure est identifié. La contrainte critique de clivage dépend de la vitesse de sollicitation. Ce modèle permet de prédire précisément par la simulation numérique, le comportement de la fissure observé expérimentalement.