Télécharger le fichier pdf d’un mémoire de fin d’études
Le circuit primaire
Le coeur du r´eacteur est plong´e dans un fluide caloporteur dont le rˆole est double : d’une part, ´evacuer la chaleur produite, et d’autre part, entretenir les r´eactions nucl´eaires en ralentissant les neutrons. Dans les r´eacteurs a` eau pressuris´ee, ce rˆole est jou´e par de l’eau a` 15.5 MPa. Le fluide atteint la temperature de 320˚C en sortie de coeur, avant de c´eder de la chaleur au circuit secondaire par le biais d’un ´echangeur thermique et recommence chaque cycle a` la temp´erature de 285˚C.
Le circuit secondaire
L’eau du circuit secondaire est vaporis´ee en recevant la chaleur du circuit primaire au coeur des g´en´erateurs de vapeur : la vapeur produite active les turbines, coupl´ees a` un alternateur et dont le mouvement permet de produire l’´electricit´. La vapeur est ensuite condens´ee et renvoy´ee vers l’´echangeur pour un nouveau cycle.
La condensation de l’eau du circuit secondaire est r´ealis´ee par le biais d’un circuit tertiaire, dont l’eau est soit refroidie a` son tour dans des tours a´eror´efrig´erantes, soit directement d´evers´ee dans un fleuve ou dans la mer.
Le crayon combustible
Le crayon combustible est l’´el´ement central du coeur du r´eacteur. Il s’agit d’un ensemble cylin-drique form´e d’un tube de gainage dans lequel sont empil´ees des pastilles de combustible en dioxyde d’uranium. Les dimensions caract´eristiques d’un crayon combustible sont donn´ees dans la figure 1.3. Initialement, les pastilles et la gaine ne sont pas en contact mais s´epar´ees par un espace appel´ jeu pastille-gaine.
Les ressorts de maintien sont n´ecessaires pour accommoder le gonflement du combustible sous l’effet de l’augmentation de la temp´erature. Le crayon est sous une pression interne d’helium d’environ 25 bar et clos herm´etiquement par deux bouchons en inconel.
Le rendement local du crayon est caract´eris´ par sa puissance lin´eique (Plin en W.cm−2) d´efinie comme la puissance g´en´er´ee par le crayon par unit´e de longueur fissile (ie. la longueur accumul´ee de l’empilement de pastilles combustibles). Les diff´erents r´egimes de fonctionnement du crayon combustible seront plus longuement abord´es par la suite, en particulier au paragraphe 1.1.3.
La gaine
La gaine constitue la premi`ere barri`ere de protection de l’environnement vis `a vis des produits de fission et son int´egrit´ doit ˆetre assur´ee quelles que soient les conditions de fonctionnement du r´eacteur. Le mat´eriau retenu est le Zircaloy, un alliage de zirconium qui poss`ede des propri´et´es ad´equates en termes de neutronique, de tenue m´ecanique et de r´esistance `a la corrosion par l’eau pr´essuris´ee `a haute temp´erature du circuit primaire.
La compr´ehension des ph´enom`enes `a l’origine de la ruine de la gaine est l’enjeu de nombreux efforts de recherche et d´eveloppement dans lesquels s’inscrit le travail pr´esent´ ici.
La pastille combustible
Dans les R´eacteurs a` Eau Pressuris´ee, le combustible est pr´esent sous forme de pastilles cylindriques de dioxyde d’uranium (voir figure 1.4) empil´ees au sein du crayon.
Les pastilles combustibles sont chanfrein´ees de fa¸con a` faciliter leur introduction dans la gaine. L’´evidement central a pour vocation d’accomoder les d´eformations du combustible pendant l’irradi-ation. Le comportement de la pastille combustible d’uranium sous irradiation sera plus longuement d´ecrit au paragraphe 1.2.
|
Table des matières
I Contexte – Etat de l’art
1 Probl´ematique industrielle
1.1 Le r´eacteur `a eau pressuris´ee
1.1.1 Principe de fonctionnement du r´eacteur
1.1.2 Le crayon combustible
1.1.3 R´egimes de puissance nominal et incidentel
1.2 Comportement du combustible UO2 sous irradiation
1.2.1 R´eactions de fission et ´evolution de la microstructure
1.2.2 Comportement thermo-m´ecanique de la pastille et Interaction Pastille-Gaine
1.3 Vers la mod´elisation du combustible : le projet PLEIADES
2 Etat de l’art pour la mod´elisation du comportement de l’UO2
2.1 Cadre de la mod´elisation
2.1.1 Sollicitations thermo-m´ecaniques recontr´ees au cours de l’irradiation
2.1.2 Essais de caract´erisation m´ecanique du mat´eriau vierge
2.1.3 Prise en compte des effets de l’irradiation sur le combustible
2.2 Mod`eles de fluage ´el´ementaires
2.2.1 Mod`eles fluage-diffusion
2.2.2 Mod`eles de fluage dislocations
2.2.3 Cartes de m´ecanismes de d´eformation
2.2.4 Fluage et irradiation
2.3 Loi Gatt-Monerie
2.3.1 Mod´elisation du fluage thermique stationnaire
2.3.2 Traitement de la porosit´e
2.4 Mod´elisation du combustible UO2 en r´eacteur
2.4.1 Mod´elisation m´ecanique de l’UO2 sous irradiation
2.4.2 Prise en compte de l’´evolution de la microstructure
2.4.3 Conclusions : les limites de l’approche multi-´echelles
II Mod`ele d’agr´egat polycristallin pour l’UO2
3 Repr´esentation de l’agr´egat
3.1 Un premier mod`ele de polycristal : l’approche `a champ moyen
3.1.1 Formulation auto-coh´erente
3.1.2 Mod´elisation de l’agr´egat d’UO2 : vers l’approche `a champ complet
3.2 Mosa¨ıque de Vorono¨ı
3.2.1 D´efinition
3.2.2 Mod´elisation du polycristal
3.2.3 Triangulation de Delaunay
3.2.4 G´en´eration et maillage de l’agr´egat
3.3 P´eriodicit´e
3.3.1 Homog´en´eisation p´eriodique
3.3.2 G´en´eration d’une mosa¨ıque de Vorono¨ı p´eriodique
3.4 Microstructures g´en´er´ees
3.4.1 Propri´et´es morphologiques du mod`ele polycristallin
3.4.2 Propri´et´es du maillage par ´el´ements finis
3.5 Conclusion sur la repr´esentation g´eom´etrique du polycristal
4 M´ecanismes de d´eformation du combustible
4.1 Elasticit´e de l’agr´egat d’UO2
4.1.1 Elasticit´e du cristal d’UO2
4.1.2 Homog´en´eisation p´eriodique
4.1.3 Analyse de l’anisotropie du polycristal d’UO2
4.1.4 Bilan des r´esultats obtenus dans le domaine ´elastique
4.2 Premi`ere approche pour la plasticit´e cristalline
4.2.1 Mouvements de dislocations dans l’UO2
4.2.2 Pic de compression et avalanche de dislocations
4.2.3 Identification des param`etres du mod`ele PCU
4.2.4 Analyse des r´esultats obtenus avec le mod`ele PCU
4.3 Mod´elisation du fluage intragranulaire
4.3.1 Prise en compte du fluage-diffusion
4.3.2 Modification du fluage-dislocations
4.3.3 Identification des param`etres du mod`ele ld2
4.3.4 Mod`eles de fluage intragranulaire : bilan
4.4 Glissement et d´ecoh´esion intergranulaire
4.4.1 Les joints de grains dans l’UO2
4.4.2 Mod`eles de zones coh´esives
4.4.3 Impl´ementation dans le mod`ele d’agr´egat
4.5 Conclusions sur les m´ecanismes de d´eformation du polycristal
III Sensibilit´e num´erique et repr´esentativit´e du VER 117
5 Premi`ere analyse du comportement local du VER 119
5.1 Pr´esentation de la simulation
5.1.1 Conditions du calcul
5.1.2 Origine des ph´enom`enes de localisation
5.1.3 R´eponse effective de l’agr´egat
5.2 Analyse des distributions de contraintes locales
5.2.1 Localisation de la contrainte intragranulaire
5.2.2 Contraintes aux joints de grains
5.3 Conclusions de l’analyse du comportement local du VER
5.3.1 Revue des r´esultats obtenus
5.3.2 Probl`emes de sensibilit´e et choix du VER
10TABLE DES MATIERES `
6 Etude de sensibilit´e au choix du maillage 131
6.1 Plasticit´e incompressible et ´el´ements finis
6.1.1 Blocage volum´etrique des ´el´ements
6.1.2 Test du poin¸con de Prandtl
6.1.3 Incompressibilit´e et effet “checkerboard”
6.2 Sensibilit´e au choix du maillage ´el´ements finis
6.2.1 Description de l’´etude
6.2.2 Sensibilit´e du comportement macroscopique
6.2.3 Comportement intragranulaire
6.2.4 Comportement aux joints de grains
6.3 Conclusions de l’´etude de sensibilit´e
7 Calcul des contraintes intergranulaires 149
7.1 Positionnement du probl`eme
7.1.1 Introduction
7.1.2 Pr´esentation de l’´etude
7.2 M´ethodes de calcul des contraintes intergranulaires
7.2.1 Extrapolation des contraintes aux noeuds
7.2.2 Une m´ethode alternative de calcul des contraintes interfaciales
7.2.3 Utilisation d’´el´ements coh´esifs
7.2.4 Analyse de la contrainte aux joints de grains par les diff´erentes m´ethodes . 157
7.2.5 Une quatri`eme m´ethode de calcul des contraintes interfaciales
7.3 Bilan de l’´etude : comparaison entre les diff´erentes m´ethodes
IV Validation du mod`ele d’agr´egat 165
8 D´ecoh´esion intergranulaire au cours d’un essai de fluage 167
8.1 Mesures exp´erimentales de la d´ecoh´esion
8.1.1 Caract´erisation des essais
8.1.2 Mesures de d´ecoh´esion intergranulaire
8.2 Simulation de l’essai
8.2.1 Un premier calcul en l’absence de d´ecoh´esion
8.2.2 Vers la prise en compte de la d´ecoh´esion intergranulaire
8.2.3 R´esultats de la simulation
8.3 Conclusions
8.3.1 Bilan des simulations de d´ecoh´esion
8.3.2 Perspectives pour la mod´elisation du comportement intergranulaire
9 Mod´elisation d’une rampe de puissance 179
9.1 Positionnement du probl`eme
9.1.1 D´efinition d’un chargement `a partir des r´esultats d’ALCYONE
9.1.2 Cas du comportement ´elastique monocristallin
9.2 Comportement effectif du polycristal
9.3 Analyse des contraintes locales
9.3.1 Pression hydrostatique intragranulaire
9.3.2 Comportement intergranulaire
9.3.3 Approche microm´ecanique et couplage multi-physique : bilan
Bibliographie 199
11TABLE DES MATIERES `
V Annexes 205
A Algorithme de tirage des germes 207
B Orientations cristallines 209
B.1 R´epartition al´eatoire
B.2 Orientations ´equir´eparties
C Conditions de p´eriodicit´e pour le calcul EF 213
C.1 Ecriture des conditions de p´eriodicit´e dans le cas g´en´eral
C.1.1 Position du probl`eme
C.1.2 Utilisation de noeuds fictifs
C.2 Elements de p´eriodicit´e pour la g´eom´etrie irr´eguli`ere
C.2.1 Cas 3D
C.3 Comparaison des deux types de g´eom´etries du VER
D Validation de la mise en oeuvre d’´el´ements coh´esifs ´elastiques 217
Table des figures 224
Liste des tableaux 22
Télécharger le rapport complet
