FR 2846139 A1 20040423 FR 0213093 A 20021021 i La présente invention concerne un procédé de détermination d'au moins une valeur limite d'au moins un premier paramètre de fonctionnement d'un réacteur nucléaire comprenant un coeur dans lequel des assemblages combustibles sont chargés, les assemblages combustibles comprenant des 5 crayons combustibles comportant chacun des pastilles de combustible nucléaire et une gaine entourant les pastilles. L'invention s'applique, par exemple, aux réacteurs à eau pressurisée. Un grand nombre de ces réacteurs sont actuellement exploités dans le monde. Il peut être utile, notamment dans des pays comme la France o 80% de l'électricité est produite par des réacteurs nucléaires, que la puissance globale fournie par les réacteurs varie afin de s'adapter aux besoins du réseau électrique qu'ils alimentent. En particulier, il est souhaitable de pouvoir faire fonctionner les 15 réacteurs à puissance réduite pendant une longue période o la demande du réseau est faible, avant de revenir si nécessaire à puissance nominale. Pour autant, une telle exploitation de chaque réacteur, qui permettrait de mieux exploiter ses capacités, ne devrait pas induire de problème de sreté. Un but de l'invention est de résoudre ce problème en fournissant un 20 procédé permettant de déterminer au moins une valeur limite d'un paramètre de fonctionnement d'un réacteur nucléaire, permettant de mieux exploiter les capacités du réacteur, tout en conservant un fonctionnement sr du réacteur. A cet effet, l'invention a pour objet un procédé de détermination d'au moins une valeur limite d'au moins un premier paramètre de fonctionnement d'un 25 réacteur nucléaire comprenant un coeur dans lequel des assemblages combustibles sont chargés, les assemblages combustibles comprenant des crayons combustibles comportant chacun des pastilles de combustible nucléaire et une gaine entourant les pastilles, caractérisé en ce qu'il comprend les étapes de: b) simuler au moins un transitoire de fonctionnement du réacteur nucléaire, c) calculer la valeur atteinte par une grandeur physique au cours du transitoire de fonctionnement dans au moins une gaine d'un crayon combustible, et d) déterminer, en tant que valeur limite, la valeur du premier paramètre de 5 fonctionnement lorsque la valeur calculée à l'étape c) correspond à une valeur de la grandeur physique caractérisant une rupture de la gaine. Selon des modes particuliers de réalisation, le procédé peut comprendre l'une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou selon toutes les combinaisons techniquement possibles 10 - il comprend avant l'étape b) une étape de: a) détermination d'une valeur de rupture de la grandeur physique caractérisant une rupture de la gaine, - l'étape a) comprend les sous-étapes de: ai) soumettre des crayons combustibles à des rampes de puissance 15 nucléaire, a2) calculer les valeurs atteintes par la grandeur physique dans au moins une gaine rompue lors d'une rampe de puissance, a3) sélectionner la valeur minimale parmi les valeurs atteintes calculées à l'étape a2), - la valeur de rupture utilisée à l'étape d) est égale à la valeur minimale sélectionnée à l'étape a3), - la valeur de rupture utilisée à l'étape d) est égale à la valeur minimale sélectionnée à l'étape a3) et corrigée par un facteur représentatif d'un mode de fonctionnement du réacteur, - le procédé comprend avant l'étape c) une étape de: b') détermination d'au moins un crayon combustible dont la gaine est la plus sollicitée lors du transitoire simulé à l'étape b), et l'étape c) est effectuée pour le ou chaque crayon déterminé à l'étape b'), l'étape b') comprend les sous-étapes de: b'1) évaluation de la valeur atteinte par la grandeur physique dans les gaines de plusieurs crayons combustibles, et b'2) sélection, en tant que crayon dont la gaine est la plus sollicitée, du crayon dont la valeur évaluée à l'étape b'1) est la plus élevée, - le premier paramètre de fonctionnement est la puissance linéique fournie par un crayon combustible, - le premier paramètre de fonctionnement est une durée de fonctionnement du réacteur à une puissance intermédiaire inférieure à sa puissance nominale, - la valeur limite est une valeur limite de déclenchement d'un arrêt d'urgence du réacteur, - le procédé comprend en outre au moins une étape de: e) détermination d'une valeur limite de déclenchement d'une alarme à 10 partir de la valeur limite d'arrêt d'urgence déterminée à l'étape d), - l'étape b) est mise en oeuvre pour au moins un mode de fonctionnement du réacteur choisi dans le groupe constitué par: * un mode de fonctionnement à puissance globale du réacteur égale à sa puissance nominale, 15 * un mode de fonctionnement prolongé à puissance intermédiaire o la puissance globale du réacteur est inférieure à sa puissance nominale sur une durée d'au moins 8 h pour une plage de 24 h, * un mode de fonctionnement en suivi de réseau o la 20 puissance globale varie alternativement autour d'une puissance élevée et autour d'une puissance réduite,. un mode de fonctionnement en réglage primaire o la puissance globale du réacteur varie de 0 à 5% autour d'une valeur de référence comprise entre 95 et 100% de 25 la puissance globale nominale du réacteur, - pour au moins un mode de fonctionnement, on utilise les étapes b) à d) d'un autre mode de fonctionnement avec, comme valeur de rupture de la grandeur physique, la valeur de rupture de cet autre mode de fonctionnement corrigé d'une valeur corrective, - le transitoire simulé à l'étape b) est un transitoire choisi dans le groupe consistant en: * une augmentation excessive de charge, * un retrait incontrôlé d'au moins un groupe de grappes de contrôle, ò une chute de l'une des grappes de contrôle, - la grandeur physique est une contrainte ou une fonction de contrainte(s) dans la gaine, et - la grandeur physique est une densité d'énergie de déformation dans la gaine. L'invention a en outre pour objet un système de détermination d'au moins un paramètre de fonctionnement d'un réacteur nucléaire, caractérisé en ce qu'il 10 comprend des moyens de mise en oeuvre des étapes d'un procédé tel que défini ci-dessus. Selon une variante, le procédé comprend au moins un ordinateur et des moyens de stockage dans lesquels sont stockés au moins un programme pour l'exécution d'étapes du procédé de détermination mis en oeuvre par le système. 15 L'invention a en outre pour objet un programme d'ordinateur comprenant des instructions pour l'exécution des étapes d'un procédé tel que défini cidessus. L'invention a en outre pour objet un support utilisable dans un ordinateur et sur lequel est enregistré un programme tel que défini ci-dessus. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple, et faite en se référant aux dessins annexés, sur lesquels: - 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 combustible du coeur du réacteur de la figure 1, - la figure 3 est une vue schématique en coupe longitudinale d'un crayon combustible de l'assemblage de la figure 2, - la figure 4 est un schéma bloc d'un système de détermination de 30 valeurs limites de paramètres de fonctionnement du réacteur de la figure 1, - la figure 5 est un organigramme illustrant des étapes successives du procédé mis en oeuvre par le système de la figure 4, - la figure 6 est une courbe illustrant des essais de rampe de puissance effectués sur des crayons combustibles, - la figure 7 est un graphique illustrant les contraintes maximales atteintes dans les gaines de crayons combustibles lors d'essais de rampe de puissance, - la figure 8 est un diagramme représentant des domaines de fonctionnement du réacteur de la figure 1, - la figure 9 est une courbe illustrant la simulation d'un transitoire de puissance lorsque le réacteur de la figure 1 est en fonctionnement prolongé à 10 puissance intermédiaire, - la figure 10 est une vue analogue à la figure 9, le transitoire étant simulé après un fonctionnement prolongé à puissance intermédiaire puis retour à un fonctionnement à puissance nominale, et - la figure 11 est une courbe illustrant le fonctionnement du réacteur en 15 suivi de réseau. La figure 1 illustre schématiquement un réacteur nucléaire 1 à eau pressurisée qui comprend classiquement: - un coeur 2, - un générateur de vapeur 3, - une turbine 4 couplée à un générateur 5 d'énergie électrique, et - un condenseur 6. Le réacteur 1 comprend un circuit primaire 8 équipé d'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 25 y être réchauffée en 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. L'eau du circuit primaire 8 alimente également le générateur de vapeur 3 o elle est refroidie en assurant la vaporisation d'eau circulant dans un circuit 30 secondaire 12. La vapeur produite par le 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. Le circuit secondaire 12 comprend en aval du condenseur 6 une pompe 13 et un réchauffeur 14. s De manière classique également, le coeur 2 comprend des assemblages combustibles 16 qui sont chargés dans une cuve 18. Un seul assemblage 16 est représenté sur la figure 1, mais le coeur 2 comprend par exemple 157 assemblages 16. Le réacteur 2 comprend des grappes de contrôle 20 qui sont 10 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 peuvent être déplacées par des mécanismes 22 pour venir s'insérer dans les assemblages combustibles 16 qu'elles surplombent. Classiquement, chaque grappe de contrôle 20 comprend des crayons de contrôle en un matériau 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 depuis la puissance nulle jusqu'à la puissance nominale 20 PN, en fonction de l'enfoncement des grappes 20 dans les assemblages combustibles 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 réglage. D'autres sont destinées à l'arrêt du réacteur I 25 et sont dénommées grappes d'arrêt. Les grappes 20 sont réunies en groupes en fonction de leur nature et de leur destination. Par exemple pour les réacteurs de type 900 MWe CPY, ces groupes sont dénommés groupes GI, G2, NI, N2, R, SA, SB, SC, SD... Comme illustré par la figure 2, chaque assemblage combustible 16 30 comprend de manière classique un réseau de crayons combustibles 24 et un squelette 26 de support des crayons 24. Le squelette 26 comprenant classiquement un embout inférieur 28, un embout supérieur 30, des tubes guides reliant les deux embouts 26 et 28 et destinés à recevoir des crayons de contrôle des grappes 20, et des grilles entretoises 32. Comme illustré par la figure 3, chaque crayon combustible 24 comprend de manière classique, une gaine 33 sous forme d'un tube fermé à son 5 extrémité inférieure par un bouchon inférieur 34 et à son extrémité supérieure par un bouchon supérieur 35. Le crayon 24 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. Classiquement, les pastilles 36 sont en oxyde d'uranium et la gaine 33 en alliage de zirconium. Sur la figure 3, qui correspond à un crayon combustible 24 issu de fabrication et avant d'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 15 la figure 3. 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 20 jeu J entre les pastilles 36 et la gaine 33, due au fluage de la gaine 33 et au gonflement des 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 25 circuit primaire 8) et l'intérieur du crayon 24, la gaine 33 se déforme progressivement par fluage 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 30 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 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 de façon progressive sur une période qui débute par un s contact mou suivi par l'établissement d'un contact fort. La pression de contact de l'oxyde de la pastille 36 sur la 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 10 déformation à la gaine 33 vers l'extérieur. En régime permanent établi, cette expansion est suffisamment 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 15 33. S'il n'y a pas de risque de rupture de la gaine 33 en régime permanent du fait 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 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 25 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 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 30 rupture de la gaine 33. Or, une telle rupture de la gaine 33 n'est pas admissible pour des raisons de sreté puisqu'elle pourrait résulter en la libération de produits de fission dans le circuit 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 5 générateur 5 alimente. Des transitoires de 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. Pour garantir l'intégrité des crayons 24 vis-à-vis de l'interaction pastilles 36 - gaines 33, l'invention propose de déterminer des valeurs limites de paramètres de fonctionnement du réacteur 1 en prenant en compte cette interaction pastille - gaine. Pour cela, on peut utiliser, par exemple, un système informatique 40 15 telle que celui de la figure 4. Ce système 40 comprend par exemple une unité de traitement d'informations 42 comprenant un ou plusieurs processeur(s), des moyens 44 de stockage des données, des moyens 46 d'entrée/sortie et éventuellement des moyens 48 d'affichage. Le procédé de détermination des valeurs limites des paramètres de 20 fonctionnement est illustré par l'organigramme de la figure 5. Une première étape illustrée par le bloc 50 consiste à déterminer une valeur d'une grandeur physique pour laquelle une rupture d'une gaine 33 d'un crayon 24 intervient lors d'une augmentation de la puissance nucléaire fournie par le crayon 24. La grandeur physique est par exemple la contrainte circonférentielle et normale ao dans la gaine 33. Dans d'autres variantes, il peut s'agir d'une fonction de contrainte(s), par exemple de la différence entre ce et la contrainte radiale et normale ar ou de la densité d'énergie de déformation dans la gaine 33. Pour déterminer la valeur de rupture, on peut soumettre un certain 30 nombre de tronçons de crayons 24, par exemple pré-irradiés dans des réacteurs de puissance et correspondant donc à des niveaux d'épuisement différents, à des rampes ou augmentation brutale de puissance dans des réacteurs d'essais. De telles rampes sont schématisées sur la figure 6 o la puissance linéique PLIN fournie par un tronçon de crayon 24 est portée en ordonnée et le temps est porté en abscisse. La courbe 52 en trait mixte sur la figure 6 illustre une rampe de s puissance subie par un tronçon de crayon 24 au cours de laquelle aucune rupture de la gaine 33 n'intervient. La forte puissance fournie par le tronçon de crayon 24 à l'issue de la rampe est représentée par la partie 54 de la courbe 52. Cette forte puissance est maintenue pendant plusieurs heures. La courbe 56 en trait plein illustre le cas d'un tronçon de crayon 24 pour lequel une rupture de la gaine 33 intervient. Dans ce cas, la rampe qui est schématisée par la partie 58 de la courbe 56, est interrompue immédiatement et la puissance PLIN ramenée rapidement à 0. On calcule, pour chaque tronçon de crayon 24, la valeur maximale 15 0OMAX atteinte par la grandeur physique ao grâce à un logiciel qui est repéré par la référence 60 sur la figure 5 et qui est par exemple stocké dans les moyens 44 de stockage du système 40. Un tel logiciel 60 peut être un logiciel classique faisant intervenir les éléments finis pour modéliser le comportement thermomécanique du tronçon de crayon 24. La figure 7 fait apparaître les valeurs aOMAx ainsi calculées pour l'ensemble des tronçons de crayons 24 qui correspondent chacun à un épuisement différent. Sur la figure 7, la contrainte est portée en ordonnée et l'épuisement en abscisse. La valeur de rupture retenue sera inférieure ou égale au minimum parmi les valeurs agMAx pour lesquelles une rupture de gaine 33 est effectivement intervenue. Dans l'exemple de la figure 7, une seule rupture est intervenue sur un crayon 24 qui est repéré par le carré plein 61. Ainsi, la valeur de rupture aORUP retenue pour la suite du procédé est d'environ 535 MPa. Ensuite, et comme schématisé par le bloc 62 sur la figure 5, on engendre les historiques de puissance pour chaque crayon 24 présent dans le coeur 2. En effet, l'état thermomécanique local d'un crayon 24 dépend de cet 5 historique. Il peut donc être utile de connaître l'historique de la puissance fournie par chaque crayon 24 depuis son introduction dans le coeur 2 jusqu'au moment d'un accident que l'on va simuler. Ces historiques peuvent être établis avec un logiciel de modélisation du comportement neutronique des crayons 24, désigné par la référence 1o numérique 63 sur la figure 5 et stocké dans les moyens 44. Un tel logiciel peut être un logiciel classique fondé sur les éléments finis. Les historiques de fonctionnement relatifs à chaque crayon 24 sont générés pour différents modes de fonctionnement du coeur 2, à savoir: - le fonctionnement en base o la puissance globale P du coeur 2 est 15 égale à sa puissance nominale PN, - le fonctionnement à puissance réduite avec les grappes de contrôle insérées dans les assemblages combustibles 16, - le fonctionnement à puissance réduite avec les grappes de contrôle extraites des assemblages combustibles 16. Les historiques peuvent être générés en prenant en compte différents niveaux de puissance réduite, par exemple 10% PN, 30% PN, 50% PN... La simulation des conditions initiales d'un accident de catégorie 2, illustrée par le bloc 64 sur la figure 5, est effectuée grâce au logiciel 63 et par exemple à l'aide d'un autre logiciel 65 de simulation du comportement de 25 l'ensemble du réacteur 1. Il peut s'agir d'un logiciel classique incluant notamment une prise en compte des régulations d'un système de protection du réacteur 1. Les conditions simulées lors de l'étape 64 correspondent à une situation dite de catégorie 1. Une telle situation est située dans le domaine 66 du diagramme de la 30 figure 8 ou le pourcentage de la puissance globale P du réacteur 1 par rapport à sa puissance nominale PN est portée en ordonnée et l'écart de puissance AI entre la partie haute et la partie basse du coeur 2 est portée en abscisse. Les conditions initiales dans le plan P/Al sont choisies de manière à maximiser l'accroissement de la puissance locale en cas de transitoire de puissance dans le coeur 2. Ces conditions initiales sont obtenues par des oscillations xénon, 5 générées par des variations instantanées de puissance globale P du réacteur 2 combinée à des mouvements de grappes de contrôle 20, conduisant à des distributions de puissance pénalisantes dans le coeur 2. Ainsi, les conditions initiales peuvent être caractérisées par: - le taux d'épuisement dans les crayons 24, - la puissance globale P fournie par le coeur 2, - la distribution axiale de puissance dans le coeur 2, - la distribution axiale de xénon, - la position des grappes de contrôle 20, et - la concentration en bore dans l'eau du circuit primaire 8. Après l'étape schématisée par le bloc 64, on simule, comme schématisé par le bloc 67, des transitoires accidentels de fonctionnement du réacteur 1 qui induisent des variations brusques de puissance. Ces transitoires sont simulés par exemple grâce aux logiciels 63 et 65. Les transitoires accidentels sont simulés à partir des conditions initiales déterminées à l'étape 64, 20 à plusieurs instants dans chaque cycle. Les transitoires simulés sont les 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 les suivants - l'augmentation excessive de charge, le retrait incontrôlé de groupes de grappes de contrôle 20, alors que le réacteur 1 est en puissance, - la chute de grappe(s) 20. L'augmentation excessive de charge correspond à une augmentation 30 rapide du débit de vapeur dans le générateur de vapeur 3. Une telle augmentation entraîne un déséquilibre entre la puissance thermique du coeur 2 et la charge du générateur de vapeur 3. Ce déséquilibre conduit à un refroidissement du circuit primaire 8. Du fait de l'effet modérateur et/ou de la régulation de la température moyenne dans le coeur 2 par les grappes de contrôle 20, la réactivité, et donc le flux nucléaire, augmentent dans le coeur 2. Ainsi, la puissance globale P fournie par le coeur 2 augmente rapidement. Pour simuler ce transitoire, on considère que le débit de vapeur dans 5 le générateur 3 augmente de sa valeur initiale jusqu'à la valeur maximale autorisée par les caractéristiques du circuit secondaire 12. Cette augmentation est par ailleurs suffisamment lente pour les niveaux de puissance étudiés afin d'éviter l'arrêt automatique du réacteur par basse pression du pressuriseur 10. Le retrait incontrôlé de groupes de grappes de contrôle 20 alors que le 10 réacteur fonctionne provoque une augmentation incontrôlée de la réactivité. Il en résulte une augmentation rapide de la puissance nucléaire globale P et du flux de chaleur dans le coeur 2. Jusqu'à l'ouverture d'une vanne de décharge ou d'une soupape de sreté du circuit secondaire 12, l'extraction de chaleur dans le générateur de vapeur 3 augmente moins vite que la puissance dégagée dans le 15 circuit primaire 8. Il en résulte un accroissement de la température et de la pression de l'eau du circuit primaire 8. Pour simuler ce transitoire, on suppose un retrait des groupes de puissance à la vitesse maximale de 72 pas/mn jusqu'à l'extraction totale. Si une ou plusieurs des grappes de contrôle 20 chute dans le coeur, il 20 en résulte une réduction immédiate de la réactivité et de la puissance globale P dans le coeur 2. En l'absence d'action protectrice, le déséquilibre ainsi induit entre le circuit primaire 8 et le circuit secondaire 12 entraîne une diminution de la température d'entrée de l'eau dans le coeur 2, ainsi qu'une augmentation de la puissance neutronique par les contre-réactions et la régulation de température, 25 jusqu'à l'atteinte d'un nouvel équilibre entre le circuit primaire 8 et le circuit secondaire 12. La présence de la ou des grappes de contrôle 20 ayant chutée provoque une déformation de la distribution radiale de puissance, tandis que le retrait du groupe de régulation conduit à une modification axiale de la puissance. On a constaté que d'autres transitoires accidentels de puissance 30 n'avaient pas besoin d'être simulés puisqu'ils s'avéraient moins pénalisants que ceux cités précédemment. C'est notamment le cas de l'accident de dilution. Après détermination des historiques de puissance et simulation des transitoires accidentels, on va calculer à l'étape 68, la valeur de la grandeur physique ao atteinte dans les crayons 24 au cours des transitoires simulés à l'étape 67. Cette étape 68 utilise les résultats des étapes 62 et 67. En fait, cette étape 68 va comprendre, pour chaque transitoire, une première sous-étape o ce calcul va être effectué grossièrement pour chacun des crayons 24. Ce calcul 5 grossier peut être effectué grâce au logiciel 60. En variante, ce calcul peut n'être effectué que pour 1/8ème des crayons 24 pour des raisons de symétrie du coeur 2. On va alors déterminer dans quel crayon 24 la valeur la plus élevée de la grandeur physique a( sera atteinte. Il s'agira donc du crayon 24 dont la gaine 33 est la plus sollicitée, ou crayon limitant. En variante, on peut sélectionner 10 plusieurs crayons limitants. Une fois le crayon 24 limitant ainsi déterminé, on va calculer de manière précise dans une deuxième sous-étape, pour ce seul crayon 24, la valeur maximale de la grandeur physique ao atteinte dans sa gaine 33. Ce calcul précis peut être effectué par le logiciel 60. Ensuite, et comme schématisé par le bloc 70, on va pouvoir déterminer par comparaison de la valeur maximale calculée à l'étape 68 avec la valeur de rupture (ffRUP déterminée à l'étape 50 si un risque de rupture de gaine 33 est encouru lors d'un transitoire accidentel. Ainsi, on va pouvoir déterminer au cours de l'étape 70 quelles sont les 20 valeurs de paramètres de fonctionnement du coeur 2 au moment d'une rupture de gaine 33 lors des transitoires accidentels simulés. Ces valeurs déterminées seront alors les valeurs limites des paramètres de fonctionnement qu'un exploitant du réacteur 1 devra respecter. Un tel paramètre de fonctionnement peut être par exemple la 25 puissance linéique PLIN dans les crayons 24. La valeur limite du paramètre ainsi déterminée peut être utilisée pour piloter le réacteur 1 et notamment pour établir des seuils d'arrêt d'urgence ou d'alarme. En effet, le seuil d'arrêt d'urgence peut être par exemple égal à la valeur de PLIN déterminée lors de l'étape 70 et les seuils d'alarme correspondent 30 à cette valeur limite minorée par exemple d'un certain pourcentage. Dans des modes particuliers de mise en oeuvre du procédé, on prend en compte lors des étapes 64 et 67, les modes de fonctionnement suivants du réacteur 1: - fonctionnement de base, le réacteur 1 fonctionnant à sa puissance globale nominale PN, - fonctionnement prolongé à puissance globale intermédiaire Pl, - le retour au fonctionnement de base après fonctionnement prolongé s à puissance globale intermédiaire Pl, - le fonctionnement en suivi de réseau, - le fonctionnement en réglage primaire, - le fonctionnement en téléréglage. La simulation d'un transitoire accidentel lors du fonctionnement de 10 base correspond, pour chaque crayon 24, à la simulation d'une rampe de puissance telle que celle représentée sur la figure 6. Le fonctionnement prolongé à puissance intermédiaire est défini comme étant le fonctionnement du réacteur 1 en régime permanent, à une puissance globale PI inférieure ou égale à environ 92% de sa puissance 15 nominale PN sur une durée de plus de 8 h par plage de 24 h. Un tel fonctionnement a pour effet de déconditionner les pastilles 36 dans les crayons 24. En effet, si localement la puissance diminue, il en résulte une diminution de température dans les pastilles 36 et dans la gaine 33, qui entraîne 20 une diminution des dilatations thermiques de ces éléments. Chaque pastille 36 possède un coefficient de dilatation thermique plus important que celui de la gaine 33 et rétrocède donc une dilatation absolue plus importante. Ce phénomène est en outre accentué par le fait que pour une baisse de puissance locale donnée, la variation de température dans chaque pastille 36 25 est plus importante que celle dans la gaine 33. Pour les crayons 24 o le contact entre la gaine 33 et les pastilles 36 n'est pas établi, le jeu radial J s'agrandit. En ce qui concerne les crayons 24 o le jeu J était fermé, le jeu J se réouvre. En cas de réouverture du jeu J, il y a fluage en compression de la 30 gaine 33 par effet de pression. Il en résulte une augmentation de la contrainte GO dans la gaine 33 quand survient le transitoire accidentel. La figure 9 illustre, pour un crayon 24, la simulation d'un transitoire accidentel lors d'un tel fonctionnement. La puissance linéique PLIN dans le crayon 24 y est portée en ordonnée et le temps t en abscisse. Le tronçon 72 de la courbe correspond au fonctionnement prolongé à 5 puissance intermédiaire. La puissance linéique PILIN dans le crayon 24 est alors inférieure à la puissance linéique nominale PNLIN correspondant à la puissance globale nominale PN. Le transitoire, qui correspond à une rampe de puissance, est schématisé par le tronçon 74 de cette même courbe. Ce transitoire se termine par une rupture de la gaine 33 au point 76. La simulation du fonctionnement prolongé à puissance intermédiaire peut être effectuée pour plusieurs valeurs de puissance intermédiaire PILIN afin de couvrir diverses situations susceptibles d'être rencontrées lors de l'exploitation du réacteur 1. En outre, cette simulation peut être effectuée pour déterminer la durée 15 maximale autorisée tmax de ce mode de fonctionnement pour que le transitoire 74 suivant le palier 72 n'induise pas de rupture de gaine 33. Pour cela, on simule des paliers 72 de durée croissante jusqu'à atteindre la valeur de rupture de la grandeur physique sélectionnée. La durée tmax est la durée du palier 72 pour lequel la valeur de rupture est atteinte. Le retour au fonctionnement de base après fonctionnement prolongé à puissance intermédiaire mérite d'être étudié vis-à-vis des interactions pastilles/gaines pour la raison suivante. Si la puissance locale augmente, il en résulte une augmentation des températures dans les pastilles 36 et dans la gaine 33 d'un crayon 24. Le 25 coefficient de dilatation thermique de chaque pastille 36 étant plus important que celui de la gaine 33, la pastille 36 impose sa déformation à la gaine 33 en y engendrant des contraintes relativement élevées. Il se produit donc un reconditionnement du combustible. Un tel fonctionnement est illustré par la courbe de la figure 10 o le 30 palier 72 à la puissance linéique PILIN est suivi d'un palier 78 à la valeur nominale PNLIN. Pour la simulation du retour au fonctionnement de base, le palier 72 de la figure 10 dure de préférence la durée tmax. Le fonctionnement en suivi de réseau est un fonctionnement dans lequel le réacteur 1 fonctionne alternativement autour d'une puissance globale réduite PR et au voisinage de sa puissance nominale globale PN. La puissance réduite peut être par exemple d'environ 30 ou 50% de PN. Par exemple, la durée 5 de fonctionnement autour de la puissance réduite peut être de 8 h et celle au voisinage de la puissance maximale PN de 16 h pour une période de 24 h. La puissance globale P peut varier dans chaque phase autour de la valeur de référence correspondante PR ou PN pour s'adapter aux besoins en énergie électrique du réseau alimenté par le générateur 5. Lorsque la valeur de référence 10 est PN, les variations de puissance ne sont que négatives. La simulation du fonctionnement dans ce mode peut être effectuée à partir du fonctionnement de base. Dans ce cas, on utilise comme valeur de rupture la valeur calculée précédemment diminuée d'une valeur corrective. Cette valeur corrective peut être 15 typiquement de l'ordre de 20 MPa lorsque la grandeur physique considérée Ga est une contrainte ou une différence de contraintes. Pour déterminer la valeur corrective on peut utiliser le logiciel 60. On simule alors grossièrement des transitoires en fonctionnement de base et des transitoires en fonctionnement en suivi de réseau. Pour simplifier les calculs, on peut modéliser le fonctionnement en suivi de réseau selon la courbe de la figure 11 o la puissance linéique PLIN est portée en ordonnée et le temps t est porté en abscisse. Les oscillations autour des valeurs haute et basse ne sont pas représentées. Pour les prendre en compte on effectue plusieurs simulations avec des puissances réduites PRLIN et 25 des puissances élevées PELIN dans les plages de variation. La différence entre les valeurs maximales de ae calculées par le logiciel 60 pour les deux modes de fonctionnement constitue alors la valeur corrective. Dans le mode de fonctionnement en réglage primaire, la puissance globale P oscille autour d'une valeur de référence de l'ordre de 98% de PN. Les 30 variations autour de cette valeur de référence peuvent être de plus 2% et de moins 3% de PN. Dans le mode de fonctionnement en téléréglage, la puissance globale P oscille autour d'une valeur de référence de l'ordre de 92% de PN. Les variations autour de cette valeur de référence peuvent être de plus 8% et de moins 7% de PN. De même que pour le fonctionnement en suivi de réseau, le fonctionnement en réglage primaire et le fonctionnement en téléréglage peuvent 5 être simulés à partir du fonctionnement de base, la valeur de la grandeur physique caractérisant la rupture étant à nouveau affectée d'un biais. Ainsi, on n'effectue qu'un seul calcul en fonctionnement de base et l'on utilise une valeur de rupture corrigée pour prendre en compte le fonctionnement de base, le fonctionnement en suivi de réseau, le fonctionnement en réglage 10 primaire et le fonctionnement en téléréglage. Dans un mode de réalisation plus simple, on peut ne tenir compte que du fonctionnement de base. La simulation du mode de fonctionnement de base, ou les différentes simulations des modes de fonctionnement, peuvent être effectuées pour chaque 15 configuration respective des groupes de grappes de contrôle 20, à différents taux d'épuisement des crayons 24. On dispose ainsi par exemple d'une valeur limite de puissance linéique PLIN pour chaque configuration de grappes de contrôle 20 et par taux d'épuisement. Outre ces valeurs limites que l'on va utiliser comme seuils d'arrêt automatique du réacteur 2, l'étape 70 de la figure 5 permet également d'établir des seuils à partir desquels d'autres actions que l'arrêt du réacteur 2 vont être automatiquement déclenchées. Il s'agit par exemple d'un seuil dit C4, qui correspond aux valeurs du seuil d'arrêt d'urgence minorées de 3%, et à partir 25 desquelles la charge de la turbine 4 sera automatiquement limitée. Il peut également s'agir d'un seuil d'alarme qui peut être réglé aux valeurs du seuil C4 minorées de 3%. Le procédé ainsi mis en oeuvre peut conduire à une limitation des domaines d'exploitation représentés sur la figure 8. Sur cette figure, la référence numérique 90 désigne le domaine correspondant aux situations de fonctionnement de catégorie 2. Ce domaine est situé autour du domaine 66. Les références numériques 92 désignent des frontières (en pointillés) de ce domaine telles qu'elles étaient établies avant l'utilisation du procédé selon l'invention. La mise en oeuvre du procédé de l'invention conduit à redéfinir certaines de ces frontières 92. Les nouvelles frontières 94 sont représentées en 5 trait plein sur la figure 8. Comme on peut le constater, le domaine 90 redéfini à partir du procédé de l'invention est donc plus étroit. Les frontières 94, qui correspondent aux seuils d'arrêt d'urgence, sont donc plus pénalisantes que celles utilisées jusqu'à présent. Le fonctionnement du réacteur 1 utilisant les paramètres limites obtenus grâce au procédé est donc plus sr. On notera que chaque tronçon de droite des frontières 94 correspond à des couples de valeurs limites de P et de AI. La durée tmax déterminée constitue également une valeur limite d'un paramètre de fonctionnement déterminée grâce au procédé. Cette valeur permet à un exploitant d'augmenter la durée de fonctionnement à puissance réduite et 1S donc d'exploiter au mieux les capacités du réacteur 2 en réduisant les risques d'endommagement des crayons 24. Ainsi, le procédé décrit permet de vérifier les valeurs limites de paramètres de fonctionnement et les procédures de sreté, de modifier les valeurs limites de paramètres de fonctionnement du réacteur 1 si nécessaire et 20 de convertir d'anciennes limites trop étroites en marge d'exploitation. Ainsi, il est possible d'assurer un fonctionnement sr de ce dernier tout en exploitant mieux ses capacités. De manière générale, l'étape 50 pourra être mise en oeuvre séparément du reste du procédé de détermination de valeurs limites de 25 paramètres du fonctionnement du réacteur 2. De manière générale également, les différentes étapes de ce procédé peuvent être mises en oeuvre par un seul et même logiciel mémorisé dans les moyens 44. Les principes ci-dessus peuvent s'appliquer à d'autres types de 30 réacteur que les réacteurs à eau pressurisée, par exemple aux réacteurs à eau bouillante. REVENDICATIONS 1. Procédé de détermination d'au moins une valeur limite (tmax) d'au moins un premier paramètre (t, PLIN, P, AI) de fonctionnement d'un réacteur nucléaire (1) comprenant un coeur (2) dans lequel des assemblages 5 combustibles (16) sont chargés, les assemblages combustibles (16) comprenant des crayons combustibles (24) comportant chacun des pastilles (36) de combustible nucléaire et une gaine (33) entourant les pastilles (36), caractérisé en ce qu'il comprend les étapes de: b) simuler au moins un transitoire de fonctionnement du réacteur nucléaire 'o (2), c) calculer la valeur atteinte par une grandeur physique (Qo) au cours du transitoire de fonctionnement dans au moins une gaine (33) d'un crayon combustible (24), et d) déterminer, en tant que valeur limite, la valeur du premier paramètre de 15 fonctionnement lorsque la valeur calculée à l'étape c) correspond à une valeur (GORUP) de la grandeur physique caractérisant une rupture de la gaine (33) . 2. Procédé selon la revendication 1, caractérisé en ce qu'il comprend avant l'étape b) une étape de: a) détermination d'une valeur de rupture (aORUp) de la grandeur physique 20 caractérisant une rupture de la gaine (33). 3. Procédé selon la revendication 2, caractérisé en ce que l'étape a) comprend les sous-étapes de: ai) soumettre des crayons combustibles (24) à des rampes de puissance nucléaire, a2) calculer les valeurs atteintes (aOMAx) par la grandeur physique (ce) dans au moins une gaine (33) rompue lors d'une rampe de puissance, a3) sélectionner la valeur minimale parmi les valeurs atteintes calculées à l'étape a2). 4. Procédé selon la revendication 3, caractérisé en ce que la valeur de 30 rupture utilisée à l'étape d) est égale à la valeur minimale sélectionnée à l'étape a3). 5. Procédé selon la revendication 3, caractérisé en ce que la valeur de rupture utilisée à l'étape d) est égale à la valeur minimale sélectionnée à l'étape a3) et corrigée par un facteur représentatif d'un mode de fonctionnement du réacteur. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend avant l'étape c) une étape de: b') détermination d'au moins un crayon combustible (24) dont la gaine (33) est la plus sollicitée lors du transitoire simulé à l'étape b), et en ce que l'étape c) est effectuée pour le ou chaque crayon (24) déterminé à l'étape b'). 7. Procédé selon la revendication 6, caractérisé en ce que l'étape b') 1o comprend les sous-étapes de: b'1) évaluation de la valeur atteinte par la grandeur physique (Co) dans les gaines (33) de plusieurs crayons combustibles (24), et b'2) sélection, en tant que crayon dont la gaine (33) est la plus sollicitée, du crayon (24) dont la valeur évaluée à l'étape b'1) est la plus élevée. 8. Procédé selon l'une des revendications précédentes, caractérisé en ce que le premier paramètre de fonctionnement est la puissance linéique (PLIN) fournie par un crayon combustible (24). 9. Procédé selon l'une des revendications 1 à 7, caractérisé en ce que le premier paramètre de fonctionnement est une durée (t) de fonctionnement du 20 réacteur (1) à une puissance intermédiaire inférieure à sa puissance nominale (PN). 10. Procédé selon l'une des revendications précédentes, caractérisé en ce que la valeur limite est une valeur limite de déclenchement d'un arrêt d'urgence du réacteur (1). 11. Procédé selon la revendication 10, caractérisé en ce qu'il comprend en outre au moins une étape de: e) détermination d'une valeur limite de déclenchement d'une alarme à partir de la valeur limite d'arrêt d'urgence déterminée à l'étape d). 12. Procédé selon l'une des revendications précédentes, caractérisé en ce 30 que l'étape b) est mise en oeuvre pour au moins un mode de fonctionnement du réacteur (1) choisi dans le groupe constitué par: - un mode de fonctionnement à puissance globale du réacteur (1) égale à sa puissance nominale, - un mode de fonctionnement prolongé à puissance intermédiaire o la puissance globale du réacteur (1) est inférieure à sa puissance nominale sur une durée d'au moins 8 h pour une plage de 24 h, - un mode de fonctionnement en suivi de réseau o la puissance globale 5 varie alternativement autour d'une puissance élevée (PN; PE) et autour d'une puissance réduite (PR), - un mode de fonctionnement en réglage primaire o la puissance globale du réacteur (1) varie de 0 à 5% autour d'une valeur de référence comprise entre et 100% de la puissance globale nominale du réacteur (2). 13. Procédé selon la revendication 12, caractérisé en ce que, pour au moins un mode de fonctionnement, on utilise les étapes b) à d) d'un autre mode de fonctionnement avec, comme valeur de rupture de la grandeur physique, la valeur de rupture de cet autre mode de fonctionnement corrigé d'une valeur corrective. 14. Procédé selon l'une des revendications précédentes, caractérisé en ce que le transitoire simulé à l'étape b) est un transitoire choisi dans le groupe consistant en: - une augmentation excessive de charge, - un retrait incontrôlé d'au moins un groupe de grappes de contrôle (20), 20 - une chute de l'une des grappes de contrôle (20). 15. Procédé selon l'une des revendications précédentes, caractérisé en ce que la grandeur physique est une contrainte ou une fonction de contrainte(s) dans la gaine (33). 16. Procédé selon l'une des revendications précédentes, caractérisé en ce 25 que la grandeur physique est une densité d'énergie de déformation dans la gaine (33). 17. Système (40) de détermination d'au moins une valeur limite d'un paramètre de fonctionnement d'un réacteur nucléaire (2), caractérisé en ce qu'il comprend des moyens (42, 44, 46, 48) de mise en oeuvre des étapes d'un 30 procédé selon l'une quelconque des revendications précédentes. 18. Système selon la revendication 17, caractérisé en ce qu'il comprend au moins un ordinateur (42) et des moyens (44) de stockage dans lesquels sont stockés au moins un programme pour l'exécution d'étapes du procédé de détermination mis en oeuvre par le système. 19. Programme d'ordinateur comprenant des instructions pour l'exécution des étapes d'un procédé selon l'une quelconque des revendications l à 16. 20. Support utilisable dans un ordinateur et sur lequel est enregistré un programme selon la revendication 19.