Soutenance de thèse de Vincent GAUTHIER
Ecole Doctorale
SCIENCES POUR L'INGENIEUR : Mécanique, Physique, Micro et Nanoélectronique
Spécialité
Sciences pour l'ingénieur : spécialité Mécanique des Solides
établissement
Aix-Marseille Université
Mots Clés
endommagement,hétérogène,MOX,,
Keywords
damage,heterogeneous,MOX,,
Titre de thèse
Modélisation des effets de l'endommagement dans les milieux hétérogènes visco-élastiques; simulation du comportement des combustibles MOX
Modelling the effects of mechanical damage in heterogeneous viscoelastic media; MOX fuel behavior simulation.
Date
Wednesday 15 December 2021
à 14:00
Adresse
Château de Cadarache,
Route de Vinon Sur Verdon,
13 115 SAINT PAUL LEZ DURANCE
Amphithéâtre
Jury
Directeur de these | M. Mihail GARAJEU | Aix Marseille Université |
Rapporteur | M. Eric LORENTZ | IMSIA |
Rapporteur | M. Martin I IDIART | Université de La Plata |
Examinateur | Mme Carole NADOT | Ecole Nationale Supérieure de Mécanique et d'Aérotechnique (ENSMA) |
Examinateur | M. Djimédo KONDO | Sorbonne Université |
Examinateur | M. Bruno MICHEL | CEA Cadarache |
Examinateur | M. Thomas HELFER | CEA Cadarache |
CoDirecteur de these | M. Renaud MASSON | CEA Cadarache |
Résumé de la thèse
La thèse sintéresse au comportement des combustibles oxydes mixtes uranium plutonium (MOX), utilisés dans les réacteurs nucléaires REP du parc français. Ce combustible ressemble à un composite constitué dune phase matrice et de deux phases dinclusion; les phases se distinguant par leur teneur en plutonium. Lors de transitoires thermiques menés sur ces combustibles irradiés, transitoires représentatifs de situations accidentelles de type perte de réfrigérant primaire, le MOX peut être le siège de dilatations différentielles pouvant potentiellement conduire à son endommagement. Lobjectif de ce travail est donc détudier les conditions dapparition et de propagation de lendommagement dans un composite particulaire représentatif du combustible MOX.
Tout dabord, on aborde la question du critère damorçage de lendommagement. On génère, à cet effet, des microstructures périodiques biphasées (matrice et une phase dinclusions) simplifiant la microstructure du composite représentatif du combustible. Les inclusions sont tirées de façon aléatoire dans lélément de volume avec lalgorithme RSA. Après avoir passé en revue les différents critères damorçage possibles pour les structures homogènes, on évalue la possibilité de les étendre aux milieux hétérogènes considérés. Le critère permettant de décrire au mieux lamorçage est un critère basé sur un bilan énergétique dans une microstructure simple : la sphère composite. Le critère damorçage obtenu est fonction uniquement de paramètres matériaux (modules élastiques, taux de restitution dénergie critique) et microstructuraux (taille des inclusions, fraction volumique).
Ensuite, on met en uvre des calculs à champs complets sur les microstructures simplifiées générées, pour une situation de chargement simplifiée, mais représentative du problème initial : composite simulé libre de contraintes sur son bord et déformation libre imposée dans les inclusions. Les calculs à champs complets sont menés avec la méthode Eléments Finis. Les deux phases ont un comportement élastique et la matrice a un comportement fragile. Deux modèles dendommagement fragile sont comparés : un modèle dendommagement local anisotrope (modèle de référence pour la simulation du comportement du combustible) et un modèle dendommagement non local isotrope (modèle champs de phase). Lors de cette comparaison, deux géométries distinctes de surface de fissure ont été observées (localisée pour DDIF-2, diffuse pour le modèle champ de phase). On retient, à lissue de cette étude, le modèle isotrope (modèle champ de phase) qui sert de référence à laquelle se rapporteront les modèles développés dans la suite de la thèse, qui sont basés sur des approches en champs moyens. Ce choix est motivé par le fait quon ne dispose pas de cadre variationnel évident pour le modèle DDIF-2.
Enfin, dans la dernière partie de la thèse, sont présentées plusieurs approches à champ moyen basées sur des principes variationnels. Tout dabord, en conservant le principe variationnel utilisé pour les calculs à champs complets, une borne supérieure du comportement effectif peut être établie. Ensuite, nous investiguons des formulations compatibles avec des principes variationnels formulés de manière incrémentale . Les potentiels correspondant à ces formulations sont ensuite linéarisés suivant deux approches distinctes : lapproche sécante modifiée, et lapproche du second ordre . Les avantages et inconvénients de chacun sont discutés pour un cas de chargement particulier. Ces considérations nous permettent de définir le modèle dendommagement qui sera à terme implémenté dans les codes de calcul du CEA afin de simuler le comportement du crayon combustible dans les situations les plus variées possibles.
Thesis resume
The thesis focuses on the behavior of mixed oxide uranium-plutonium (MOX) fuel, used in the French PWR nuclear reactors. This fuel resembles a composite consisting of a matrix phase and two inclusion phases; the phases are distinguished by their plutonium content. During thermal transients conducted on these irradiated fuels, transients being representative of accidental situations such as loss of primary coolant, the MOX fuel can be the seat of differential expansions that can potentially lead to its damage. The objective of this work is therefore to study the conditions of damage occurrence and propagation in a particulate composite representative of MOX fuel.
First of all, we address the question of the damage initiation criterion. For this purpose, periodic two-phase microstructures (matrix and a phase of inclusions) simplifying the microstructure of the composite representative of the fuel are generated. The inclusions are randomly drawn in the volume element with the RSA algorithm. After reviewing the different possible onset of damage criteria for homogeneous structures, the possibility of extending them to the heterogeneous media considered is evaluated. The criterion that best describes the initiation is a criterion based on an energy balance in a simple microstructure: the composite sphere. The criterion obtained depends only on material parameters (elastic moduli, critical energy restitution rate) and microstructural parameters (inclusion size, volume fraction).
Then, full field calculations are performed on the simplified microstructures generated, for a simplified loading situation, but representative of the initial problem: simulated composite free of stresses on its edge and free strain imposed in the inclusions. The full field calculations are carried out with the Finite Element Method. The two phases have an elastic behavior and the matrix has a brittle behavior. Two brittle damage models are compared: a local anisotropic damage model (reference model for the simulation of the fuel behavior) and a non local isotropic damage model (phase field model). During this comparison, two distinct crack surface geometries were observed (localized for DDIF-2, diffuse for the phase field model). At the end of this study, the isotropic model (phase field model) is retained as a reference to which the models developed in the rest of the thesis, which are based on mean field approaches, will refer. This choice is motivated by the fact that there is no obvious variational framework for the DDIF-2 model.
Finally, in the last part of the thesis, several mean-field approaches based on variational principles are presented. First, by retaining the variational principle used for full-field calculations, an upper bound on the effective behavior can be established. Then, we investigate formulations compatible with variational principles formulated in an incremental way . The potentials corresponding to these formulations are then linearized following two distinct approaches: the modified secant approach , and the second order approach. The advantages and disadvantages of each are discussed for a particular loading case. These considerations allow us to define the damage model that will eventually be implemented in the CEA calculation codes in order to simulate the behavior of the fuel rod in the most varied situations possible.