FR 3085788 A1 20200313 FR 1858153 A 20180912 Procédé de protection d'un réacteur nucléaire et réacteur nucléaire correspondant L'invention concerne en général la protection des réacteurs nucléaires contre les risques de rupture de gaine de crayons de combustible par interaction pastille-gaine (IPG). Aujourd'hui, sur les réacteurs à eau pressurisée français équipés de quatre boucles primaires, la protection vis-à-vis du risque de rupture de gaine par IPG est assurée par un arrêt automatique du réacteur sur puissance linéique élevée. En d'autres termes, le système de protection du réacteur évalue périodiquement, avec une constante de temps courte, la puissance linéique des crayons du combustible nucléaire et déclenche un arrêt automatique du réacteur si la puissance linéique dans un tronçon de crayon dépasse une valeur maximale prédéterminée. La rupture de la gaine par IPG résulte d'un phénomène thermomécanique. La puissance linéique est un paramètre qui n'est pas entièrement représentatif du phénomène physique conduisant au risque. Pour tenir compte de cet écart, il est nécessaire d'établir une corrélation enveloppe entre le phénomène physique conduisant au risque et la puissance linéique. Ceci est pénalisant vis-à-vis des marges de manoeuvres pour le pilotage du réacteur. Dans ce contexte, l'invention vise à proposer un procédé de protection d'un réacteur nucléaire contre les risques de rupture de gaines de crayons de combustible par IPG qui est moins pénalisant pour le pilotage de réacteur. A cette fin, l'invention porte selon un premier aspect sur un procédé de protection d'un réacteur nucléaire comportant un coeur ayant une pluralité d'assemblages de combustible nucléaire, chaque assemblage comportant une pluralité de crayons, chaque crayon comprenant une gaine et un combustible nucléaire enfermé dans la gaine, le procédé comportant les étapes suivantes : - Etape (S12) : Reconstruction d'une puissance linéique maximale dégagée parmi les crayons des assemblages de combustible nucléaire du coeur ; - Etape (S14) : Calcul de l'état thermomécanique et du taux d'épuisement des crayons ; - Etape (S24) : Calcul d'une contrainte mécanique ou d'une densité d'énergie de déformation dans la gaine d'un des crayons en utilisant ladite puissance linéique maximale reconstruite, les états thermomécaniques calculés et les taux d'épuisement calculés, à l'aide d'un méta-modèle d'un code de thermomécanique ; - Etape (S30) : Comparaison de la contrainte mécanique calculée ou de la densité d'énergie de déformation calculée avec un seuil respectif ; 2 - Etape (S32) : Arrêt du réacteur nucléaire si la contrainte mécanique calculée ou la densité d'énergie de déformation calculée dépasse ledit seuil respectif. Comme expliqué plus haut, la rupture de la gaine par IPG est due à un phénomène thermomécanique, le paramètre représentatif étant la contrainte mécanique 5 dans les gaines en zircaloy, ou la densité d'énergie de déformation (DED) pour les gainages M5. L'invention permet de calculer en ligne, dans le système de protection du réacteur nucléaire, l'évolution de la contrainte ou de la densité d'énergie de déformation, et de la comparer à un seuil. C'est le dépassement du seuil qui déclenche l'arrêt automatique du 10 réacteur. Ce seuil correspond directement à la limite technologique conduisant à la rupture. En d'autres termes, la limite technologique du combustible devient directement le seuil déclenchant l'arrêt automatique du réacteur. Le réacteur est donc protégé de la rupture par IPG à partir de la surveillance directe du paramètre physique thermomécanique conditionnant la rupture. 15 Le calcul en ligne de la contrainte ou de la DED est effectué avec une constante de temps compatible avec les exigences de protection du réacteur grâce à l'utilisation d'un méta-modèle du code de thermomécanique. Un tel méta-modèle permet de calculer la contrainte mécanique ou la densité d'énergie de déformation extrêmement rapidement, beaucoup plus rapidement qu'un code de thermomécanique classique fonctionnant 20 typiquement par éléments finis. Du fait que le réacteur est protégé directement à partir de la surveillance du paramètre physique conditionnant la rupture, il n'est pas nécessaire d'utiliser une corrélation enveloppe du type mentionné ci-dessus, de telle sorte que les marges de manoeuvrabilité du réacteur sont améliorées. 25 Le procédé de protection peut également présenter une ou plusieurs des caractéristiques ci-dessous, considérées individuellement ou selon toutes les combinaisons techniquement possibles : - l'étape de reconstruction de la puissance linéique maximale est effectuée en utilisant des mesures fournies en continu par des sondes de mesure d'un flux neutronique ; 30 - l'étape de reconstruction est effectuée avec une constante de temps inférieure à 1 minutes ; - l'étape de calcul de l'état thermomécanique et du taux d'épuisement des crayons comprend une sous-étape de calcul de la distribution de puissance en trois dimensions dans le coeur, et une sous-étape de calcul du taux d'épuisement des crayons 35 en utilisant la distribution de puissance en trois dimensions calculée ; 3 - l'étape de calcul de l'état thermomécanique et du taux d'épuisement des crayons comprend une sous-étape de calcul de de l'état thermomécanique des crayons à l'aide d'un code de thermomécanique, en utilisant la distribution de puissance en trois dimensions calculée ; 5 - l'étape de calcul de l'état thermomécanique et du taux d'épuisement des crayons comprend une sous-étape de stockage de l'état thermomécanique et du taux d'épuisement calculé pour un tronçon témoin d'un des crayons ; - le tronçon témoin est le tronçon dégageant la puissance linéique maximale ; - le méta-modèle du code de thermomécanique utilise au moins les entrées suivantes : 10 contrainte mécanique ou densité d'énergie de déformation calculée au pas de temps précédent, puissance linéique maximale reconstruite, dérivée par rapport au temps de la puissance linéique maximale reconstruite, taux d'épuisement calculé pour le tronçon témoin ; - le méta-modèle du code de thermomécanique est un réseau de neurones multicouches ; 15 - un biais prédéterminé est appliqué au seuil considéré pour l'étape de comparaison ; - le procédé comporte une étape de détermination du biais, comprenant les sous-étapes suivantes : - Sous-étape (S34): Calcul de l'état thermomécanique initial et du taux d'épuisement initial des crayons ; 20 - Sous-étape (S36) : Calcul de l'évolution dans le temps de la contrainte mécanique ou de la densité d'énergie de déformation dans la gaine des crayons à partir de l'état thermomécanique initial et du taux d'épuisement initial des crayons calculés, pour au moins un transitoire accidentel, à l'aide dudit code de thermomécanique ; - Sous-étape (S38) : Détermination du maximum de la contrainte mécanique ou de 25 la densité d'énergie de déformation calculée à l'aide dudit code de thermomécanique ; - Sous-étape (S42) : Calcul de l'évolution dans le temps de la contrainte mécanique ou de la densité d'énergie de déformation dans la gaine des crayons à partir de l'état thermomécanique initial et du taux d'épuisement initial du combustible à l'endroit de la puissance linéique maximale des crayons calculés, pour ledit au moins un transitoire 30 accidentel, à l'aide dudit méta-modèle du code de thermomécanique, en prenant en compte une simulation de l'évolution de ladite puissance linéique maximale reconstruite à l'étape (S12); - Sous-étape (S46) : Détermination du maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculée à l'aide dudit méta-modèle du code de 35 thermomécanique ; 4 - Sous-étape (S48) : Détermination du biais, sur la base de l'écart entre d'une part la contrainte mécanique ou de la densité d'énergie de déformation calculé à l'aide dudit code de thermomécanique et d'autre part le maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculé à l'aide dudit méta-modèle du code de 5 thermomécanique. Selon un second aspect, l'invention porte sur un réacteur nucléaire comprenant : - un coeur ayant une pluralité d'assemblages de combustible nucléaire, chaque assemblage comportant une pluralité de crayons, chaque crayon comprenant une gaine et un combustible nucléaire enfermé dans la gaine, 10 - un ensemble de surveillance et de protection du réacteur, comprenant un calculateur configuré pour mettre en oeuvre le procédé ayant les caractéristiques ci-dessus. Selon un troisième aspect, l'invention porte sur un programme d'ordinateur comprenant les instructions pour l'exécution des étapes du procédé ayant les caractéristiques ci-dessus. 15 Selon un quatrième aspect, l'invention porte sur un support utilisable dans un ordinateur et sur lequel est enregistré le programme ci-dessus. D'autres caractéristiques et avantages ressortiront de la description détaillée qui en est donné ci-dessous, à titre indicatif également limitatif, en référence aux figures annexées, parmi lesquelles : 20 - la figure 1 est une vue schématique illustrant un réacteur nucléaire à eau pressurisée ; - la figure 2 est une vue schématique latérale d'un assemblage de combustible du coeur de réacteur de la figure 1 ; - la figure 3 est une vue schématique en coupe longitudinale d'un crayon de combustible de l'assemblage de la figure 2 ; 25 - la figure 4 est un diagramme d'étape illustrant le procédé de protection de l'invention ; - la figure 5 est une représentation schématique simplifiée d'un méta-modèle du type utilisé dans le procédé de la figure 4 ; et - la figure 6 est un diagramme d'étape détaillant l'étape préalable de détermination du biais appliqué à l'état thermomécanique considéré pour l'étape de calcul de la contrainte 30 mécanique ou de la DED dans le procédé de la figure 4. La figure 1 illustre schématiquement un réacteur nucléaire 1 qui comprend un coeur 2. Le réacteur nucléaire 1 est typiquement un réacteur à eau pressurisée, et comprend : 35 - au moins un générateur de vapeur 3, - une turbine 4 couplée à un générateur 5 d'énergie électrique, et 5 - un condenseur 6. Le réacteur 1 comprend un circuit primaire 8 équipé d'au moins une pompe 9 et dans lequel circule de l'eau sous pression, selon le trajet matérialisé par les flèches sur la figure 1. Cette eau remonte notamment au travers du coeur 2 pour y être réchauffée en 5 assurant la réfrigération du coeur 2. Le circuit primaire 8 comprend en outre un pressuriseur 10 permettant de mettre sous pression l'eau circulant dans le circuit primaire 8 et de contrôler ladite pression. L'eau du circuit primaire 8 alimente également le ou chaque générateur de vapeur 3 où elle est refroidie en assurant la vaporisation d'eau circulant dans un circuit 10 secondaire 12. La vapeur produite par le ou chaque générateur 3 est canalisée par le circuit secondaire 12 vers la turbine 4 puis vers le condenseur 6 où cette vapeur est condensée par échange de chaleur indirect avec de l'eau de refroidissement circulant dans le condenseur 6. 15 Le circuit secondaire 12 comprend en aval du condenseur 6 une pompe 13 et un réchauffeur 14. De manière classique également, le coeur 2 comprend une pluralité d'assemblages 16 de combustible nucléaire. Ces assemblages 16 sont chargés dans une cuve 18. Un seul assemblage 16 est 20 représenté sur la figure 1, mais le coeur 2 comprend par exemple 157 assemblages 16. Le réacteur 1 comprend des grappes de commande 20 qui sont disposées dans la cuve 18 au-dessus de certains assemblages 16. Une seule grappe 20 est représentée sur la figure 1, mais le coeur 2 peut comprendre par exemple environ 60 grappes 20. Les grappes 20 sont être déplacées par des mécanismes 22 pour venir s'insérer 25 dans les assemblages de combustible 16 qu'elles surplombent. Classiquement, chaque grappe de commande 20 comprend des crayons comportant un ou plusieurs matériaux absorbant les neutrons. Ainsi, le déplacement vertical de chaque grappe 20 permet de régler la réactivité du réacteur 1 et autorise des variations de la puissance globale P fournie par le coeur 2 30 depuis la puissance nulle jusqu'à la puissance nominale PN, en fonction de l'enfoncement des grappes 20 dans les assemblages 16. Certaines de ces grappes 20 sont destinées à assurer la régulation du fonctionnement du coeur 2, par exemple en puissance ou en température, et sont dénommées grappes de contrôle. D'autres sont destinées uniquement à l'arrêt du 35 réacteur 1 et sont dénommées grappes d'arrêt. 6 Les grappes 20 sont réunies en groupes en fonction de leur nature et de leur destination. Par exemple pour la majorité des réacteurs de type 900 MWe, ces groupes sont dénommés groupes G1, G2, N1, N2, R, SA, SB, SC, SD... Le réacteur 1 comprend également un certain nombre de capteurs de mesure de 5 valeurs effectives de paramètres de fonctionnement du réacteur, notamment un thermocouple 21A de mesure de la température moyenne de l'eau du circuit primaire en sortie de la cuve 18 et un thermocouple 21B de mesure de la température moyenne de l'eau du circuit primaire en entrée de la cuve 18. De manière classique également, le réacteur nucléaire 1 comprend des sondes de 10 mesure du flux neutronique dans le coeur 2 du réacteur, fonctionnant de manière continue. Typiquement ces sondes sont des chambres externes 21 C disposées autour de la cuve 18. Le nombre et les positions des chambres 21C, habituellement dénommées « chambres ex-core », varient en fonction du modèle du réacteur 1. 15 De manière classique encore, le réacteur 1 comprend des thermocouples 21D disposés dans le coeur 2 au-dessus d'assemblages 16 pour mesurer la température de l'eau du circuit primaire en sortie des assemblages 16. Une seule chambre 21C et un seul capteur 21 D ont été représentés sur la figure 1. En plus ou à la place des chambres ex-core, certains réacteurs comprennent des 20 chambres internes fixes de mesure en continu du flux neutronique (non représentées), lesquelles chambres sont disposées à l'intérieur de certains assemblages combustibles placés dans le coeur. Le nombre et les positions de ces chambres fixes, habituellement dénommées « in-core », varient en fonction du modèle du réacteur. Ces chambres sont par exemple des collectrons au cobalt. 25 Les chambres ex-core 21C et les thermocouples 21 D ou les chambres « in-core » fixes fournissent en continu des informations relatives à la répartition tridimensionnelle de la puissance locale dans le coeur. Le réacteur nucléaire 1 comprend encore des organes permettant de mesurer ou de reconstruire les grandeurs suivantes : positions des grappes 20 dans le coeur, pression 30 dans le pressuriseur 10, vitesse de rotation de la ou chaque pompe primaire 9. Comme illustré par la figure 2, chaque assemblage 16 comprend de manière classique une pluralité de crayons de combustible 24. Ces crayons sont disposés en réseau et maintenus en position par un squelette 26 de support des crayons 24. Le squelette 26 comprend classiquement un embout inférieur 28, un embout 35 supérieur 30, des tubes-guides 31 reliant les deux embouts 30 et 28 et destinés à recevoir des crayons des grappes de commande 20 et des grilles entretoises 32. 7 Comme illustré par la figure 3, chaque crayon de combustible 24 comprend de manière classique, une gaine 33 et un combustible nucléaire enfermé dans la gaine 33. La gaine 33 est un tube longitudinal fermé à son extrémité longitudinale inférieure par un bouchon inférieur 34 et à son extrémité longitudinale supérieure par un bouchon 5 supérieur 35. Le combustible nucléaire comprend une série de pastilles 36 empilées dans la gaine 33 et prenant appui contre le bouchon inférieur 34. Un ressort de maintien 40 est disposé dans le tronçon supérieur de la gaine 33 pour prendre appui sur le bouchon supérieur 35 et sur la pastille 36 supérieure. 10 Classiquement, les pastilles 36 sont à base d'oxyde d'uranium et la gaine 33 est par exemple en un alliage de zirconium appelé zircaloy. Sur la figure 3, qui correspond à un crayon de combustible 24 issu de fabrication et avant irradiation, il existe un jeu radial J entre les pastilles 36 et la gaine 33. Cela est illustré plus particulièrement par la partie cerclée agrandie de la figure 3. 15 Lorsque le réacteur 1 va fonctionner, par exemple à sa puissance globale nominale PN, le combustible nucléaire des pastilles 36 va être, selon le terme employé dans la technique, conditionné. Le conditionnement se caractérise essentiellement par la fermeture du jeu J entre les pastilles 36 et la gaine 33, due au fluage de la gaine 33 et au gonflement des 20 pastilles 36. Plus spécifiquement, on peut distinguer les étapes suivantes pour chaque pastille 36 : 1) Sous l'effet de la différence de pression entre l'extérieur (eau du circuit primaire 8) et l'intérieur du crayon 24, la gaine 33 se déforme progressivement par fluage 25 radialement vers l'intérieur du crayon 24. Toutes choses égales par ailleurs, la vitesse de fluage de la gaine 33 est une caractéristique du matériau qui la constitue. Par ailleurs, les produits de fission, qui sont en majorité retenus dans la pastille 36, induisent un gonflement de la pastille 36. Durant cette phase, la sollicitation de la gaine 33 du point de vue des contraintes est uniquement le fait de la pression différentielle existant entre 30 l'extérieur et l'intérieur du crayon 24. Les contraintes dans la gaine 33 sont des contraintes de compression (négatives par convention). 2) Le contact entre la pastille 36 et la gaine 33 s'amorce au bout d'un laps de temps qui dépend essentiellement des conditions locales d'irradiation (puissance, flux neutronique, température...) et du matériau de la gaine 33. En réalité, le contact s'établit 35 de façon progressive sur une période qui débute par un contact mou suivi par l'établissement d'un contact fort. La pression de contact de l'oxyde de la pastille 36 sur la 8 face interne de la gaine 33 conduit à une inversion des contraintes dans la gaine 33 qui deviennent positives et tendent à solliciter la gaine en traction. 3) Le gonflement de la pastille 36 se poursuit et impose alors sa déformation à la gaine 33 vers l'extérieur. En régime permanent établi, cette expansion est suffisamment 5 lente pour que la relaxation du matériau de la gaine 33 permette un équilibre des efforts dans la gaine 33. Une analyse montre que dans ces conditions le niveau des contraintes en traction est modéré (quelques dizaines de MPa) et ne présente pas de risque vis-à-vis de l'intégrité de la gaine 33. S'il n'y a pas de risque de rupture de la gaine 33 en régime permanent du fait 10 de l'équilibre thermomécanique dans la gaine 33 à des niveaux de contraintes assez faibles, un risque apparaît par contre dès que la puissance fournie par le crayon 24 varie fortement. En effet, une augmentation de puissance engendre une augmentation de la température dans le crayon 24. Etant donnée la différence des caractéristiques 15 mécaniques (coefficient de dilatation thermique, module d'Young) et la différence de température entre la pastille 36 en oxyde d'uranium et la gaine 33 en alliage de zirconium, la pastille 36 va se dilater davantage que la gaine 33 et imposer sa déformation à cette dernière. Par ailleurs, la présence dans l'espace entre la gaine 33 et la pastille 36 de 20 produits de fission corrosifs, tels que l'iode, crée les conditions de la corrosion sous contrainte. Ainsi, la déformation imposée par la pastille 36 à la gaine 33 lors d'un transitoire, ou variation, de puissance peut provoquer une rupture de la gaine 33. Or, une telle rupture de la gaine 33 n'est pas admissible pour des raisons de sûreté puisqu'elle pourrait résulter en la libération de produits de fission dans le circuit 25 primaire 8. Des transitoires de puissance peuvent se produire lors d'un fonctionnement normal du réacteur 1, c'est-à-dire dans les situations dites de catégorie 1. En effet, des variations de puissance peuvent être nécessaires notamment pour s'adapter aux besoins en énergie électrique du réseau que le générateur 5 alimente. Des transitoires de 30 puissance peuvent également se produire dans les situations accidentelles dites de catégorie 2, comme l'augmentation excessive de charge, le retrait incontrôlé de groupe(s) de grappes 20 en puissance, la dilution d'acide borique ou encore la chute de grappes 20 non détectée. Le réacteur 1 comporte également un ensemble 41 de surveillance et de 35 protection du réacteur, comprenant une unité de traitement d'informations 42. L'ensemble 41 comporte un système de protection du réacteur 43. 9 Le système de protection du réacteur 43 est configuré pour protéger le réacteur nucléaire 1 contre un certain nombre de risques, et en particulier contre les risques de rupture de gaine de crayon de combustible par IPG. Le système de protection du réacteur 43 reçoit les informations relevées par 5 les capteurs équipant le réacteur nucléaire, et déclenche automatiquement l'arrêt du réacteur nucléaire en cas de risque, notamment en cas de risque de rupture de gaine de crayon de combustible par IPG. L'acquisition des informations des capteurs et le déclenchement de l'arrêt automatique du réacteur sont effectués avec une constante de temps courte, compatible avec la cinétique des phénomènes physiques à l'origine des 10 risques à traiter. L'ensemble 41 comporte également un système de surveillance du réacteur 44. Le système de surveillance du réacteur 44 est configuré pour surveiller en permanence le fonctionnement du coeur. 15 Le système de surveillance du réacteur 44 reçoit les informations relevées par les capteurs équipant le réacteur nucléaire, et calcule des grandeurs caractéristiques du fonctionnement du réacteur. Il déclenche des alarmes quand ces grandeurs sortent du domaine de fonctionnement autorisé. L'ensemble de surveillance et de protection du réacteur 41 est renseigné 20 notamment par les thermocouples 21A, 21 B et 21 D. Il est également renseigné, notamment, par les chambres ex-core 21C et/ou les chambres in-core, et par les organes agencés pour mesurer ou reconstruire les positions des grappes 20 dans le coeur, la pression dans le pressuriseur 10, et la vitesse de rotation de la ou des pompes primaires 9. 25 L'unité de traitement d'information 42 est typiquement formée d'un processeur et d'une mémoire associé au processeur. Différents logiciels ou codes de calcul sont stockés dans la mémoire, comme décrit plus bas. En variante, l'unité de traitement d'information est réalisée sous la forme de composants logiques programmables, tels que des FPGA (de l'anglais Field- 30 Programmable Gate Array), ou encore sous forme de circuit intégré dédié, tels que des ASIC (de l'anglais Application-Specific Integrated Circuit). Pour garantir l'intégrité des crayons 24 vis-à-vis de l'interaction pastille 36 gaine 33, l'invention propose un procédé de protection au cours duquel la contrainte 35 mécanique ou la densité d'énergie de déformation dans la gaine d'au moins un des crayons est calculée, le procédé comprenant une étape d'arrêt du réacteur si la contrainte 10 mécanique calculée ou la densité d'énergie de déformation calculée dépasse un seuil prédéterminé respectif. L'unité de traitement d'information 42 est configurée pour mettre en oeuvre le procédé de protection du réacteur nucléaire qui va être décrit ci-dessous. 5 Ce procédé est illustré sur la figure 4. Le procédé de protection comporte les étapes suivantes : - étape S12 : reconstruction d'une puissance linéique maximale P,;, max dégagée parmi les crayons 24 des assemblages de combustible nucléaire 16 du coeur 2 ; - étape S14 : calcul de l'état thermomécanique et du taux d'épuisement des crayons 24 ; 10 - étape S24: calcul d'une contrainte mécanique ou d'une densité d'énergie de déformation (DED) dans la gaine 33 d'un des crayons 24 en utilisant ladite puissance linéique maximale reconstruite, les états thermomécaniques calculés et les taux d'épuisement calculés, à l'aide d'un méta-modèle d'un code de thermomécanique ; - étape S30 : comparaison de la contrainte mécanique calculée ou de la densité d'énergie 15 de déformation calculée avec un seuil respectif ; - étape S32 : arrêt du réacteur 1 si la contrainte mécanique calculée ou la densité d'énergie de déformation calculée dépasse ledit seuil respectif. Les différentes étapes ci-dessus sont effectuées de manière répétée, avec une constante de temps qui sera décrite ci-dessous. 20 Au cours de l'étape S12, le système de protection 43 reconstruit notamment la puissance linéique maximale dégagée par chacun des crayons des assemblages de combustible nucléaire, en utilisant typiquement les informations fournies par les chambres ex-core 21C ou par les chambres in-core. Cette reconstruction est effectuée par un logiciel 46, intégré dans l'unité de 25 traitement d'informations 42. Ce logiciel est connu en soi et ne sera pas décrit en détail ici. On considère ici la puissance dégagée par unité de longueur longitudinale du crayon. En d'autres termes, l'étape S12 comporte une opération permettant de déterminer 30 comment la puissance maximale dégagée par unité de longueur varie au cours du temps. L'étape de reconstruction de la puissance linéique maximale P n max est ainsi effectuée en utilisant des mesures fournies en continu par les sondes de mesure du flux neutronique (chambres ex-core 21C ou chambres in-core). L'étape de reconstruction S12 est effectuée périodiquement, avec une constante 35 de temps inférieure à une minute, de préférence inférieure à 10 secondes, encore de préférence inférieure à 1 sseconde, typiquement de l'ordre de 0.5 s. 11 L'étape S14 de calcul de l'état thermomécanique et du taux d'épuisement des crayons comprend : - une sous-étape S16 de calcul de la distribution de puissance en trois dimensions dans le coeur 2 ; et 5 - une sous-étape S18 de calcul du taux d'épuisement des crayons en utilisant la distribution de puissance en trois dimensions calculée. L'état thermomécanique considéré ici est le la contrainte mécanique a, par exemple la contrainte tangentielle en peau interne aeo ou la différence entre la contrainte tangentielle aeo et la contrainte radiale radiale a,. 10 L'étape S14 est effectuée en utilisant des logiciels intégrés dans le système de surveillance 44. Les sous-étapes S16 et S18 sont effectuées par un logiciel 47 utilisant un code neutronique du coeur. Ce logiciel utilise par exemple le code SMART de la chaîne SCIENCE. Il est connu et ne sera pas décrit en détail ici. 15 Le logiciel 47 calcule la distribution de puissance en trois dimensions dans le coeur en utilisant les mesures fournies en continue par les chambres ex-core 21C et les thermocouples 21 D ou les chambres in-core fixes. La demande FR 2787196 décrit par exemple une méthode pour obtenir la répartition de puissance en trois dimensions dans le coeur en utilisant au moins ces informations. 20 Typiquement, le logiciel 47 utilise la puissance thermique du coeur qui est calculée par bilan enthalpique, le déséquilibre axial de puissance fourni par les chambres ex-core 21C, la position des grappes, les température et débit d'entrée cuve, la pression du pressuriseur et les températures de sortie du coeur fournies par les thermocouples 21 D. A la sous-étape S16, la puissance est calculée en tout point du coeur 2. 25 A la sous-étape 518, le taux d'épuisement TE est calculé en tout point du coeur 2. Il est calculé plus précisément selon un maillage axial fixé par le maillage axial du code de neutronique utilisé par le système de surveillance 44, typiquement par le code SMART. Le taux d'épuisement local résulte de l'intégration au cours du temps de la puissance locale. 30 L'étape S14 comporte encore une sous-étape S20 de calcul de l'état thermomécanique des crayons 24 à l'aide d'un code de thermomécanique, en utilisant la distribution de puissance en trois dimensions calculée. Le code de thermomécanique est par exemple le code COPERNIC 48 développé par la société Framatome. Ce code de calcul 48 est connu et ne sera pas détaillé ici. Il est 35 typiquement implémenté dans le système de surveillance 44. 12 A la sous-étape S20, l'état thermomécanique est calculé en tout point du coeur, c'est-à-dire en tout point de chacun des crayons 24 du coeur 2. L'étape S14 comporte également une sous-étape S22 de stockage de l'état thermomécanique et du taux d'épuisement calculé pour un tronçon témoin d'un des 5 crayons 24. Le tronçon témoin est de préférence le tronçon dégageant la puissance linéique maximale. Ce tronçon témoin est déterminé en utilisant les données fournies par le logiciel 47 intégré au système de surveillance 44. Ces informations sont stockées dans une mémoire 49 de l'unité de traitement 10 d'informations 42. L'étape S14 est effectuée de manière répétée, avec une constante de temps de l'ordre de l'heure. En effet, l'état thermomécanique des crayons 24 est essentiellement conditionné par la durée de fonctionnement du réacteur nucléaire à puissance intermédiaire, qui est normalement d'un ou plusieurs jours. Le calcul de l'état 15 thermomécanique des crayons à l'aide du code de thermomécanique COPERNIC 48 demande environ une heure, ce qui est largement satisfaisant. Du fait que l'état thermomécanique et le taux d'épuisement calculés pour le tronçon témoin sont utilisés par le système de protection 43 du réacteur nucléaire, il est nécessaire que cette information soit disponible de manière pérenne, même en cas de 20 perte du système de surveillance 44, qui n'a pas la même classification que le système de protection 43. Ainsi, il est prévu que le stockage de ces informations dans la mémoire 49 ait la même classification que le système de protection 43 du point de vue de la sureté. En variante, il est possible de prévoir que le système de surveillance 44 ait la même classification que le système de protection 43 du point de vue de la sureté. 25 La contrainte mécanique a calculée à la sous-étape S24 est par exemple la contrainte circonférentielle et normale a9e dans la gaine 33. En variante, il s'agit de la différence entre la contrainte circonférentielle et normale a00 et la contrainte radiale et normale ar, appelée contrainte de TRESCA. Typiquement, la contrainte mécanique a est calculée à l'étape S24 quand la gaine 30 33 est en zircaloy. Pour des combustibles à gainage M5, il est prévu à l'étape S24 de calculer la densité d'énergie de déformation dans la gaine d'un des crayons 33. La densité d'énergie de déformation (DED) est définie de la manière suivante : 13 DED = I croodsoo où 098 est la contrainte tangentielle en peau interne de la gaine 33 en MPa, et e68 est la déformation tangentielle totale en peau interne de la gaine 33. L'intégration débute lorsque la vitesse de déformation en peau interne dépasse 5 une valeur de seuil fixée à 10-5 s-1, et s'arrête lorsque la vitesse passe en dessous de ce seuil. Du fait que le calcul de la contrainte ou de la DED doit être utilisé pour la protection du réacteur 1 lors d'accidents pouvant avoir une cinétique rapide, celui-ci ne peut pas être effectué par un code de thermomécanique tel que COPERNIC. 10 Pour répondre à l'exigence de temps de réponse de la protection, il est prévu d'utiliser un méta-modèle 50 du code de thermomécanique à l'étape S24. Le méta-modèle 50 est par exemple un méta-modèle du code de thermomécanique COPERNIC. Le méta-modèle 50 est typiquement un réseau de neurones multicouches. Le méta-modèle 50 utilise une fonction non-linéaire comme fonction de transfert, 15 par exemple sigmoïde. Le méta-modèle 50 utilise à chaque pas de temps, c'est-à-dire à chaque itération du calcul, au moins les entrées suivantes : - contrainte mécanique ou densité d'énergie de déformation calculée au pas de temps précédent: a(t-dt) ; 20 - puissance linéique maximale reconstruite : Plinm"(t) ; - dérivée par rapport au temps de la puissance linéique maximum reconstruite : d131,,'"(t)/dt, - taux d'épuisement calculé pour le tronçon témoin : TE(t). La contrainte mécanique ou la densité d'énergie de déformation au pas de temps 25 courant a(t) constitue la seule sortie du méta-modèle. La puissance linéique maximale reconstruite est celle calculée à l'étape S12. Le taux d'épuisement calculé pour le tronçon témoin est celui mémorisé à la sous-étape S22. Le calcul de la contrainte mécanique ou de la densité d'énergie de déformation est effectué par le méta-modèle 50 selon l'algorithme général : 30 a(t) = 0-(t - dt) + A (CiPhlaxdt (t)) lit B Prez (t) où A et B sont des constantes. Plus précisément, comme illustré sur la figure 5, le réseau de neurones comporte une couche d'entrée 52 située à gauche sur la figure 5, une couche de sortie 54 située à 14 droite de la figure 5, et une ou plusieurs couches intermédiaires 56 appelées encore couches cachées. Sur la figure 5, deux couches intermédiaires 56 sont représentées. Toutefois, le réseau de neurones pourrait comporter un nombre quelconque de couches intermédiaires. 5 La couche d'entrée comporte un neurone pour chacune des entrées scalaires mentionnées ci-dessus, à savoir quatre neurones. La couche de sortie 54 comporte un unique neurone, correspondant à la grandeur d'intérêt Y reconstruite en ligne à savoir la contrainte ou la DED selon la formule factorisée suivante : me) rn i=nH+lo-i. Ail. X0 1-1 10 La ou chaque couche cachée 56 comporte un nombre déterminé de neurones, par exemple au moins cinq neurones. Le nombre de couches cachées et le nombre de neurones dans la ou chaque couche cachée est déterminé pendant la phase de paramétrage ou d'apprentissage du réseau de neurones. 15 Les entrées constituent un vecteur, appelé vecteur d'entrée et noté Xo. Dans le cas dpnlax de l'invention : X0 6(t I"' (t),11;:ax (1),TE dt Chaque couche cachée i du réseau de neurones élabore un vecteur de valeur Xi, par traitement du vecteur de valeur X,_, issu de la couche immédiatement précédente i-1. Ainsi, pour chaque couche cachée i, Xi = fi. Ai. Xi_l 20 où i est compris entre 1 et nH le nombre de couches cachées, A est une matrice prédéterminée de dimension Ni*Nki, N; étant le nombre de neurones de la couche i et N,_, étant le nombre de neurones de la couche i-1. f, est un opérateur non linéaire prédéterminé. L'opérateur f est par exemple une fonction de transfert sigmoïde agissant sur chaque composante du vecteur. 25 Les opérateurs non linéaires f, agissent chacun comme une fonction non linéaire par exemple sigmoïde sur toutes les composantes du vecteur A,X0 . Une fonction sigmoïde est une fonction s'exprimant de la manière suivante : f(x)=1/(1+e-lx) où I est une constante prédéterminée 30 Chaque matrice A agit comme une transformation affine sur les composantes du vecteur XHI. En d'autres termes : AtXt_i = WiXt_i + bt 15 où Wi est la matrice des poids synaptiques qui relient les N, neurones de la couche i aux NINI neurones i-1, et bi est le vecteur des biais de dimension Ni de la couche i. Pour la dernière couche on notera que la matrice AnH+1 est de dimension 1 x Nriry et que anH+1= 1 : il n'y a pas d'application de fonction de transfert non linéaire pour établir le 5 signal Y de la dernière couche. A l'étape S30, la contrainte mécanique calculée à l'étape S24 est comparée à un seuil prédéterminé, ou la densité d'énergie de déformation calculée à l'étape S24 est comparée à un seuil prédéterminé. 10 Ledit seuil est égal à la limite technologique du gainage considéré : contrainte de TRESCA conduisant à la rupture pour le gainage zircaloy, DED conduisant à la rupture pour le gainage M5. A l'étape S32, une procédure d'arrêt automatique du réacteur est enclenchée si la contrainte mécanique calculée ou la densité d'énergie de déformation calculée dépasse 15 ledit seuil respectif. Cet arrêt automatique est réalisé en descendant à grande vitesse les grappes d'arrêt 20 dans le coeur 2 du réacteur nucléaire. Cette opération est entièrement automatique. Le risque de rupture de gaine ne survient pas nécessairement au point chaud du 20 coeur, c'est-à-dire au niveau du tronçon de gaine voyant la puissance linéique maximale. Le risque de rupture est plus grand au contraire au niveau d'un tronçon de crayon un peu moins puissant, mais d'avantage fragilisé par un taux d'épuisement plus élevé. Le calcul de la contrainte ou de la DED au point chaud du coeur, c'est-à-dire dans le tronçon de gaine voyant la puissance linéique maximale, n'est donc pas entièrement 25 conservatif. Avantageusement, dans le processus de l'invention, un biais prédéterminé est appliqué à la limite technologique du gainage considéré, constituant le seuil de déclenchement de l'arrêt automatique du réacteur à l'étape S32. Le biais consiste en une diminution de la limite technologique. 30 Dans ce cas, le procédé comporte avantageusement une étape S10 de détermination du biais. L'étape S10 est effectuée au cours d'une phase d'étude, appelée phase d'étude IPG, qui intervient avant la mise en oeuvre des étapes S12, S14, S24, S30 et S32. Cette étape S10 est effectuée une seule fois et n'est pas répétée. 35 Le biais est prédéterminé au sens où le même biais est utilisé à chacune des itérations du procédé de l'invention. 16 Comme illustré sur la figure 6, l'étape S10 de détermination du biais comprend les sous-étapes suivantes : - sous-étape S34 : calcul de l'état thermomécanique initial et du taux d'épuisement initial des crayons 24 ; 5 - sous-étape S36 : calcul de l'évolution dans le temps de la contrainte mécanique ou de la densité d'énergie de déformation dans la gaine 33 des crayons 24 à partir de l'état thermomécanique initial et du taux d'épuisement initial des crayons 24 calculés, pour au moins un transitoire accidentel, à l'aide d'un code de thermomécanique ; - sous-étape S38: détermination du maximum de la contrainte mécanique ou de la 10 densité d'énergie de déformation calculée à l'aide dudit code de thermomécanique ; - sous-étape S42 : calcul de l'évolution dans le temps de la contrainte mécanique ou de la densité d'énergie de déformation dans la gaine 33 des crayons 24 à partir de l'état thermomécanique initial et du taux d'épuisement initial du combustible à l'endroit de la puissance linéique maximale des crayons 24 calculés, pour ledit au moins un transitoire 15 accidentel, à l'aide du méta-modèle 50 du code de thermomécanique, en prenant en compte une simulation de l'évolution de ladite puissance linéique maximale reconstruite à l'étape S12; - sous-étape S44: détermination du maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculée à l'aide du méta-modèle 50 du code de 20 thermomécanique ; - sous-étape S46 ; détermination du biais, sur la base de l'écart entre d'une part le maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculée à l'aide dudit code de thermomécanique et d'autre part le maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculée à l'aide du méta-modèle 50 25 du code de thermomécanique. Ces sous-étapes sont toutes effectuées en bureau d'étude, lors des études de sûreté IPG. La sous-étape S34 est effectuée en simulant l'évolution des crayons 36 au cours de l'enchaînement des cycles de gestion du combustible pour le réacteur nucléaire. On 30 considère ici les cycles standards d'exploitation du réacteur nucléaire. Cette simulation est effectuée en utilisant par exemple un logiciel de modélisation du comportement neutronique des crayons 24. L'état thermomécanique initial des crayons est déterminé en utilisant un code de thermomécanique, par exemple COPERNIC. 17 Le ou les transitoires accidentels simulés à la sous-étape 36 sont des transitoires accidentels dits de catégorie 2, induisant les variations de puissance les plus fortes et les plus rapides dans le coeur 2. Ces transitoires sont par exemple les suivants : - augmentation excessive de charge ; 5 - retrait incontrôlé de groupes de grappes de contrôle 20, alors que le réacteur est en puissance ; - chute de grappes 20. Ces transitoires sont décrit plus en détail dans FR2846139. A la sous-étape S36, l'évolution de la contrainte mécanique ou de la densité 10 d'énergie de déformation est calculée en utilisant un code de thermomécanique tel que COPERNIC. Cette évolution est calculée typiquement pour la gaine 33 de tous les crayons 24. Elle est calculée en tous points de la gaine de chaque crayon. Le maximum déterminé à la sous-étape S38 correspond à la valeur la plus élevée 15 de contrainte mécanique ou de DED rencontrée pour un tronçon de gaine au cours du ou de chaque transitoire. La longueur du tronçon considéré est fixée par le maillage axial du code de thermomécanique. L'étape S10 comporte encore une étape S40 de calcul de l'évolution dans le temps de la puissance linéique maximale reconstruite par le système de protection 43 20 pour ledit au moins un transitoire accidentel. A la sous-étape S40, on simule donc la puissance linéique maximale qui serait reconstruite à l'étape S12 du procédé. L'étape S10 comporte également une sous-étape S41 de stockage de l'état thermomécanique et du taux d'épuisement d'un tronçon témoin. Ce tronçon témoin 25 correspond au tronçon de crayon 24 dégageant la puissance linéique maximale à l'état initial des crayons 24. Cet état initial est celui calculé à la sous-étape S34. La sous-étape S42 est effectuée à l'aide du méta-modèle 50 utilisé à l'étape S24. L'évolution de la puissance linéique maximale reconstruite, déterminée à la sous-étape 40, constitue une des entrées de ce méta-modèle, pour la sous-étape S42. L'état 30 thermomécanique et le taux d'épuisement du tronçon témoin obtenu à la sous-étape S41 constitue une autre entrée du méta-modèle 50 pour la sous-étape S42. A la sous-étape S44, le maximum déterminé correspond à la valeur la plus élevée de contrainte mécanique ou de DED rencontrée au point chaud du coeur tel que l'estime le système de protection 43. 35 A la sous-étape S46, le biais est déterminé en comparant le maximum obtenu à la sous-étape S38 et le maximum obtenu à la sous-étape S44, pour le ou chaque transitoire. 18 Le biais à appliquer à l'état thermomécanique initial du tronçon témoin est choisi pour que le maximum calculé par le méta-modèle 50 à l'étape S44 soit toujours supérieur ou égal au maximum calculé par le code de thermomécanique à l'étape S38, pour tous les transitoires considérés. 5 Le biais est par exemple un écart de contrainte ou de DED défalqué de la limite technologique correspondante. Le procédé de protection décrit ci-dessus a été décrit en mentionnant l'utilisation de codes de calcul tels que SMART et COPERNIC. En variante, il est mis en oeuvre en utilisant d'autres codes de calcul équivalents. 10 Le procédé de l'invention peut être mis en oeuvre en utilisant tout type de système de surveillance, sous réserve que celui-ci permette de déterminer en ligne la distribution de puissance en trois dimensions, et qu'il intègre un code de thermomécanique configuré pour déterminer le conditionnement du combustible. De même, il a été décrit ci-dessus que le procédé de protection prévoyait le calcul 15 de la contrainte mécanique ou de la densité d'énergie de déformation, en fonction du type de gaine des crayons de combustible. Il serait possible de calculer à la fois la contrainte et la densité d'énergie de déformation. Le procédé de protection a été décrit dans une application à réacteur nucléaire de type à eau pressurisée (REP ou PWR). En variante, REVENDICATIONS1.- Procédé de protection d'un réacteur nucléaire (1 ), le réacteur nucléaire (1) comportant un coeur (2) ayant une pluralité d'assemblages de combustible nucléaire (16), chaque assemblage (16) comportant une pluralité de crayons (24), chaque crayon (24) comprenant une gaine (33) et un combustible nucléaire enfermé dans la gaine (33), le procédé comportant les étapes suivantes : - Etape (S12) : Reconstruction d'une puissance linéique maximale dégagée parmi les crayons (24) des assemblages de combustible nucléaire (16) du coeur (2) ; - Etape (S14) : Calcul de l'état thermomécanique et du taux d'épuisement des crayons (24); - Etape (S24) : Calcul d'une contrainte mécanique ou d'une densité d'énergie de déformation dans la gaine (33) d'un des crayons (24) en utilisant ladite puissance linéique maximale reconstruite, les états thermomécaniques calculés et les taux d'épuisement calculés, à l'aide d'un méta-modèle (50) d'un code de thermomécanique ; - Etape (S30) : Comparaison de la contrainte mécanique calculée ou de la densité d'énergie de déformation calculée avec un seuil respectif ; - Etape (S32) : Arrêt du réacteur nucléaire (1) si la contrainte mécanique calculée ou la densité d'énergie de déformation calculée dépasse ledit seuil respectif. 2.- Procédé selon la revendication 1, dans lequel l'étape (S12) de reconstruction de la puissance linéique maximale est effectuée en utilisant des mesures fournies en continu par des sondes de mesure d'un flux neutronique (21C). 3.- Procédé selon la revendication 1 ou 2, dans lequel l'étape de reconstruction (S12) est effectuée avec une constante de temps inférieure à 1 minutes. 4.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape (S14) de calcul de l'état thermomécanique et du taux d'épuisement des crayons (24) comprend une sous-étape (S16) de calcul de la distribution de puissance en trois dimensions dans le coeur (2), et une sous-étape (S18) de calcul du taux d'épuisement des crayons (24) en utilisant la distribution de puissance en trois dimensions calculée. 5.- Procédé selon la revendication 4, dans lequel l'étape (S14) de calcul de l'état thermomécanique et du taux d'épuisement des crayons (24) comprend une sous-étape (S20) de calcul de de l'état thermomécanique des crayons (24) à l'aide d'un code de thermomécanique, en utilisant la distribution de puissance en trois dimensions calculée. 6.- Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape (S14) de calcul de l'état thermomécanique et du taux d'épuisement des crayons (24) 20 comprend une sous-étape (S22) de stockage de l'état thermomécanique et du taux d'épuisement calculé pour un tronçon témoin d'un des crayons (24). 7.- Procédé selon la revendication 6, dans lequel le tronçon témoin est le tronçon dégageant la puissance linéique maximale. 5 8.- Procédé selon la revendication 6 ou 7, dans lequel le méta-modèle du code de thermomécanique utilise au moins les entrées suivantes : contrainte mécanique ou densité d'énergie de déformation calculée au pas de temps précédent, puissance linéique maximale reconstruite, dérivée par rapport au temps de la puissance linéique maximale reconstruite, taux d'épuisement calculé pour le tronçon témoin. 10 9.- Procédé selon l'une quelconque des revendications précédentes, dans lequel le méta- modèle du code de thermomécanique est un réseau de neurones multicouches. 10.- Procédé selon l'une quelconque des revendications précédentes, dans lequel un biais prédéterminé est appliqué au seuil considéré pour l'étape de comparaison (S30). 11.- Procédé selon la revendication 10, dans lequel le procédé comporte une étape (510) 15 de détermination du biais, comprenant les sous-étapes suivantes : - Sous-étape (S34) : Calcul de l'état thermomécanique initial et du taux d'épuisement initial des crayons (24) ; - Sous-étape (S36) : Calcul de l'évolution dans le temps de la contrainte mécanique ou de la densité d'énergie de déformation dans la gaine (33) des crayons (24) à partir de l'état 20 thermomécanique initial et du taux d'épuisement initial des crayons (24) calculés, pour au moins un transitoire accidentel, à l'aide dudit code de thermomécanique ; - Sous-étape (S38) : Détermination du maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculée à l'aide dudit code de thermomécanique ; - Sous-étape (S42) : Calcul de l'évolution dans le temps de la contrainte mécanique ou de 25 la densité d'énergie de déformation dans la gaine (33) des crayons (24) à partir de l'état thermomécanique initial et du taux d'épuisement initial du combustible à l'endroit de la puissance linéique maximale des crayons (24) calculés, pour ledit au moins un transitoire accidentel, à l'aide dudit méta-modèle du code de thermomécanique, en prenant en compte une simulation de l'évolution de ladite puissance linéique maximale reconstruite à 30 l'étape (S12); - Sous-étape (S44) : Détermination du maximum de la contrainte mécanique ou de la densité d'énergie de déformation calculée à l'aide dudit méta-modèle du code de thermomécanique ; - Sous-étape (S46) : Détermination du biais, sur la base de l'écart entre d'une part la 35 contrainte mécanique ou de la densité d'énergie de déformation calculé à l'aide dudit code de thermomécanique et d'autre part le maximum de la contrainte mécanique ou de 21 la densité d'énergie de déformation calculé à l'aide dudit méta-modèle du code de thermomécanique. 12.- Réacteur nucléaire (1) comprenant : - un coeur (2) ayant une pluralité d'assemblages de combustible nucléaire (16), chaque 5 assemblage (16) comportant une pluralité de crayons (24), chaque crayon (24) comprenant une gaine (33) et un combustible nucléaire enfermé dans la gaine (33) ; - un ensemble (41) de surveillance et de protection du réacteur, comprenant une unité de traitement d'informations (42) configurée pour mettre en oeuvre le procédé de l'une quelconque des revendications 1 à 11. 10