Brevet FR2773612A1 — publié le 19990716

L'invention est relative à un procede de détermination par le calcul, en cas d'effet chaudière, de la pression atteinte dans le fluide contenu dans l'espace interopercules d'une vanne à opercules. Les robinets-vannes à double étanchéité tels que vannes à sièges parallèles, vannes à cônes flexibles, sont des organes sensibles à l'effet chaudière qui est la montée en pression du fluide contenu dans l'espace inter-opercules et le dôme corps-chapeau de la vanne lors d'une montée en température de l'environnement extérieur. Ce phénomène est susceptible de se produire lors d'un accident thermodynamique interne à l'enceinte du bâtiment d'un réacteur au cours duquel la température ambiante atteint l56C et la pression atteint une valeur de 5,6 bar, ou par conduction de la chaleur dans les tuyauteries adjacentes. Les essais d'effet chaudière s'étant révélés coûteux et longs car nécessitant des installations et des phases de préparation, de dépouillement et d'analyse importantes, l'intérêt de la simulation numérique s'est imposée rapidement. Une méthodologie originale s'appuyant sur les résultats issus d'essais sur deux vannes permet de déterminer la pression maximale atteinte dans le fluide pour un transitoire thermique imposé. L'invention a donc pour objet un procédé de détèrmi- nation de la pression maximale atteinte par le fluide contenu dans le volume défini par l'espace inter-opercules et le dôme corps/chapeau d'un robinet-vanne à double étanchéité, lors d'une montée de température déterminée de l'environnement extérieur au robinet-vanne, caractérisé en ce qu'il consiste - à effectuer un calcul thermique ou thermohydraulique pour déterminer le champ de température du fluide contenu dans ledit volume lors de ladite montée en température, - à effectuer un calcul mécanique pour déterminer les déplacements de la surface intérieure du robinet-vanne définissant ledit volume pour des changements internes en pression évoluant dans une plage de pression contenant la pression maximum à déterminer, - à calculer à partir des résultats du calcul mécanique la variation dudit volume interne du robinet-vanne pour chaque valeur de pression dans ladite plage de pression, - à calculer à partir des résultats du calcul thermique ou thermohydraulique l'expansion volumique du fluide contenu dans ledit volume, et - à déterminer, à partir des résultats des calculs de la variation du volume interne d'une part et de l'expansion volumique du fluide contenu dans ledit volume, d'autre part ainsi que du fait que le fluide occupe tout ledit volume, la pression maximale atteinte pour ladite montée en température choisie. Cette méthodologie s'appuie sur des calculs thermiques et mécaniques par éléments finis suivis d'un posttraitement analytique. Suivant des caractéristiques particulières de l'invention - les phases de calcul thermique ou thermohydraulique et de calcul mécanique sont des phases de calcul numérique, - le calcul thermique ou thermohydraulique est un calcul numérique transitoire effectué sur un modèle qui comporte le corps de la vanne et le volume fluide qu'elle renferme. L'invention sera mieux comprise à la lecture de la description qui va suivre faite en référence au dessins annexés, donnés uniquement à tire d'exemple et sur lesquels - la Fig.l représente en perspective un pentaèdre élémentaire utilisé pour le calcul de l'expansion volumique du corps de la vanne soumis à la pression interne; - la Fig.2 est une vue en coupe d'une demi-vanne à opercule montrant les conditions limites pour le calcul thermique; - la Fig.3 est une coupe analogue de la Fig.2 montrant l'implantation des éléments de contact et la pression imposée en vue du calcul mécanique; - la Fig.4 est une vue en perspective d'un élément de goujonnerie; - la Fig.5 est une vue en perspective de l'assemblage corps-chapeau montrant un goujon dans un modèle global; - la Fig.6 est une vue extérieure du modèle global comprenant un corps, un chapeau, un siège et un opercule; - la Fig.7 est un exemple de maillage de la vanne de la figure 6; - la Fig.8 est un exemple de maillage pour le calcul mécanique; - la Fig.9 représente le détail du maillage du corps et de la goujonnerie de la vanne; et - la Fig.lO est un organigramme illustrant le logiciel de post-traitement utilisé pour l'exécution des opérations sur les résultats du calcul thermiques et mécaniques exécutés dans le cadre du procédé de l'invention. Les calculs réalisés jusqu'à présent l'on été avec des codes appartenant à la Demanderesse, à savoir - le code N3S qui est un code d'analyse thermohydraulique par éléments finis, et - le code ASTER qui est un code d'analyse thermomécanique par éléments finis. Ces calculs pourraient également être réalisés par des codes industriels possédant les mêmes spécialités telles que ANSYS, NASTRAN, etc. Les calculs thermiques ou thermohydrauliques ont pour but la détermination du champ de température dans le fluide, considére pour l'instant comme un milieu continu où les termes de convection sont négligés dans un but de simplification. La convection peut-être prise en compte par un calcul thermohydraulique avec le code N3S-Syrthès. Le calcul mecanique est mené avec le code ASTER. Le but est de calculer le déplacement des noeuds internes d'un maillage, représentant le corps du robinet pour en déduire la variation de volume interne de ce corps sous l'effet de la pression. Ce calcul permet accessoirement de déterminer des contraintes dans le corps de la vanne et de faire l'analyse de l'intégrité mécanique de l'ensemble. Le post-traitement consiste en une résolution d'équation différentielle qui identifie la variation du volume d'eau, sous l'effet de la dilatation et la variation du volume de l'intérieur du robinet, sous l'effet de la pression afin de calculer la pression atteinte pour un chargement thermique donné. Ce post-traitement est effectué à l'aide d'un logiciel appelé ALADIN écrit à cet effet et décrit par la suite. Le procédé de détermination de la pression du fluide consiste à déterminer la pression atteinte dans le fluide par un post-traitement analytique de deux calculs en éléments finis. - un calcul thermique ou thermohydraulique transitoire sur un modèle qui comporte le corps de vanne et le volume fluide qu'il renferme, - un calcul mécanique sur le corps uniquement. Le calcul thermique ou thermohydraulique a pour but de calculer, pour chaque pas de temps le champ de température dans le fluide pour un chargement en peau externe de la vanne correspondant à la montée en température de l'environnement. Ce champ de température permet de calculer ensuite l'expansion volumique du fluide sous l'effet de la dilatation. Le calcul mécanique se fait à température constante (température maximale de l'environnement) car l'expérience a démontré que le corps du robinet atteint rapidement cette valeur. Par ailleurs on néglige la dilatation du corps de la vanne car il existe un rapport de 10 à 100 (selon les températures) entre le coefficient de dilatation de l'eau et celui de l'acier. Le chargement mécanique consiste à imposer en peau interne une rampe de pression et à y récupérer les déplacements pour en déduire la variation du volume qu'elle engendre. Cette rampe de pression est censee envelopper la plage dans laquelle la pression varie au cours du transitoire thermique. On considère un volume élémentaire dV du fluide. Sa variation relative de volume dépend de la pression et de la température : d(dV)/dV = a(T)dT - idP où a(T) est le coefficient de dilatation du fluide à la température T et 4 le taux de compressibilité de l'eau (cette valeur est quasiconstante entre 0 et 200in et entre 0 et 1000 bar). Si l'on intègre la relation précitée sur l'ensemble du domaine fluide Vf on obtient sa variation relative de volume Le volume du fluide est égal à chaque instant à celui de la coque noté Vc qui le renferme soit : dVf/Vf = dVc/Vc = f(P)dP, cette fonction f(P)étant déduite par post-traitement du calcul mecanique. Le passage d'un état (T,P) à un état (T+dT, P+dP), T désignant le champ de température dans le fluide, impose la relation suivante Le passage d'un état de référence (ToP0) à un état (T1,Pl) conduit à la formule suivante Le post-traitement du calcul thermique consiste à calculer Le calcul analytique final de détermination de la pression consiste à résoudre l'équation On va maintenant exposer le calcul de l'expansion volumique du corps de la vanne soumis à la pression interne. Pour chaque valeur de pression, on détermine les valeurs des deplacements des noeuds des éléments de la peau interne du corps de vanne. Le calcul consiste à faire la somme des volumes des pentaèdres tel que celui représenté à la Fig.l, constitués par trois noeuds à leur position initiale Nl, N2, N3 et trois noeuds à leur position déplacée N1 + dN1, N2 + dN2 et N3 + dN3. Chacun de ces pentaèdres étant irrégulier, le calcul de son volume consiste à décomposer celui-ci en trois tétraèdres : (N1, N2, N3, N1 + dN1), (N1 + dN1, N2, N3, N2 + dN2) et (N1 + dNl,N2 + dN2, N3, N3 + dN3). Le volume d'un tétraèdre (N1, N2, N3, N4) se calcule aisément en posant n = N1N2 fl2 = N1N3 = N1N4 V = (n1 # n2).n3 6 Ainsi on obtient le volume des pentaèdres correspondant à chaque élément du maillage. A partir de ces informations on deduit dVc en utilisant : le 1 pour les points de la surface interne de la coque et on en soustrait les volumes des pentaèdres. Se pose alors le problème du sens de la variation de Vc correspondant à chaque volume de pentaèdre . Enfin on calcule Vc. On va maintenant décrire le calcul de l'expansion volumique du fluide à partir du champ de température calculé. Le calcul final fait intervenir le terme T désignant un champ de température et T désignant une température. La méthode employée consiste à discrétiser le volume fluide à l'aide des éléments finis qui ont servi au calcul thermique. où Vi est le volume de l'élément i calculé selon une méthode décrite ci-après et Ti, la température moyenne de celle des quatre sommets d'un tétraèdre résultant de la décomposition du pentaèdre de la Fig.l. Ensuite on procède à l'analyse mécanique de la vanne. Le calcul mécanique a pour but la détermination du champ de déplacement des noeuds de la peau interne de la structure. Les calculs mécaniques peuvent se faire indifféremment en élastique ou en élastoplastique. On veille simplement à adapter au mieux la loi de comportement du matériau aux conditions de chargement. Si on risque d'atteindre de hautes pressions, un comportement élastoplastique est mieux adapté. Ce calcul permet également la détermination du champ de contraintes dans le corps du robinet pour un changement de pression interne donné. Il s'agit de vérifier que la pression maximale atteinte au cours d'un transitoire, généralement très supérieure à la pression maximale admissible, PMA, ne provoque pas de plastification du corps de la vanne. Ce type de calcul fait intervenir des eléments de contact unilatéral au niveau de la jonction corps/chapeau, qui comprend un joint d'étanchéité non pris en compte ici compte tenu de sa faible importance mecanique et au niveau du contact opercule/portée d'étanchéité. La raison première est qu'un modèle monobloc est trop rigide par rapport à la réalité et qu'un calcul réaliste doit prendre en compte le glissement relatif des pièces. La seconde raison est qu'à partir des décollements éventuels à la liaison corps/chapeau, on peut établir un pronostic de fuite servant en quelque sorte de "soupape" et limitant la montée en pression. La Figure 2 représente en coupe une demi-vanne à opercule comprenant un corps 1 et un chapeau 2 réunis par une jonction corps/chapeau 3. Le corps 1 est pourvu d'un siège 4 pour un opercule 5. Cet ensemble définit un volume interne 6. On suppose qu'il y a continuité entre les deux milieux, c'est-à-dire entre le fluide circulant dans la vanne et le corps metallique de celle-ci. L'expérience a montré que c'est la représentation la plus réaliste. La condition limite pour le calcul thermique est déterminé par une loi de Fourier 2 (grade T. aH (Tt -T) Le coefficient H est une fonction du temps comme le montre le tableau II ci-après. Les conditions limites pour le calcul mécanique expriment - la liaison corps-chapeau - l'interface entre les deux morceaux du maillage fait intervenir des élements de contact unilatéral. Il est important d'en tenir compte car des essais antérieurs ont montré que la structure se déforme différemment si la liaison y est rigide. Tel est le cas pour la zone de contact Z1 sur la figure 3 - la liaison opercule-portée d'étanchéité qui fait également intervenir du contact unilatéral pour simuler le glissement relatif des composants sous un chargement de pression au niveau de la zone de contact Z2 de la figure 3 - la pression interne imposée - le tensionnement de la goujonnerie. La figure 4 est une vue en perspective schématique représentant en détail les conditions aux limites sur la gouj onnerie. Un goujon 10 est représenté traversant le contact entre le corps 1 et le chapeau 2 représentés par des portions de maillage 11 et 12. Un maillage 13 represente le contact entre le chapeau 2 et un écrou (non représenté) vissé sur le goujon 10 en vue d'assembler le chapeau 2 sur le corps 1. L'effort de tensionnement de la goujonnerie est obtenu en imposant un déplacement entre les noeuds du haut du goujon 10 et leur vis-à-vis sur le chapeau 2. Les noeuds du bas du goujon 10 sont lies aux noeuds coincidents du corps 1. On calcule la raideur équivalente du goujon 10 selon la méthode suivante. Un goujon se compose de trois parties, à savoir deux parties filetées situées de part et d'autre d'un fût lisse. La raideur globale K du goujon est telle que 1/K = l/Kl + 1/K2 + 1/K3 où Kl = E Sf(Lf1/2) est la raideur de la première partie filetée, E étant le module d'Young de l'acier et Sf étant la section efficace de la zone filetée. K2 = E Si/Li est la raideur de la partie lisse centrale, S1 etant la section de la partie lisse. K3 = E Sf/(Lf2/2) est la raideur de la seconde partie filetée. La section équivalente est donc donnée par la relation Se = (Lf1 + Li + Lf2)/(2Sf/Lf1 + 2Sf/Lf2 + S1/L1) On calcule ensuite le déplacement g à imposer pour obtenir le tensionnement requis compte tenu de la raideur du goujon et on effectue un calcul "à froid" c'est-à-dire sans pression et avec les caractéristiques des matériaux à la température ambiante. Généralement l'assemblage global est moins rigide que le goujon seul de sorte que la contrainte calculée est inférieure à celle attendue. Un ajustement du tensionnement g conduit alors à la bonne valeur du serrage. Sur la Fig.5 on a représente schématiquement en perspective et en coupe l'assemblage corps-chapeau montrant le détail d'un goujon dans un modèle global. Généralement un quart de modèle suffit pour représenter une vanne. Dans certains cas on doit représenter une demi-vanne. La liaison de la vanne avec la tuyauterie est également à prendre en considération dans certains cas. Les extrématés des tuyauteries sont alors considérées comme des bords libres. Dans le cas d'un calcul où l'échauffement passe par la tuyauterie, il conviendrait de considérer un modèle intégrant la rigidité et les moments imposés par celle-ci, notamment dans le cas de vannes extérieures au bâtiment du réacteur. Le tableau I ci-après donne les caractéristiques physiques qui interviennent dans tous les calculs Tableau I Donnes physiques relatives à l'eau Température Coefficient de Conductivité Chaleur mas (-C) dilatation a thermique sique (10-41/K) (W/mK) (J/kg.K) 20 1,82 0,599 4183 40 3,87 0,635 4174 60 5,11 0,659 4179 80 6,32 0,674 4195 100 7,52 0,683 4220 120 8,64 0,686 4250 140 9,72 0,685 4287 160 10,7 O, 683 4346 La nature de l'acier du corps des vannes dépend du cas traité. Le tableau Il contient les données standard pour le calcul thermique (transitoire et coefficients d'échange en peau externe). Tableau II Instant (sec.) Temp,(-C) H(W/aK) O O 40 11 000 1 40 11 000 2 160 11 000 60 160 10 000 90 155 9 000 900 155 9 000 1200 9 000 1800 140 7 000 On va maintenant décrire quelques modèles de vannes réalisés en CAO ainsi que les maillages nécessaires pour le calcul. La Fig.6 est une vue extérieure d'un modèle global comprenant comme sur les figures précédentes un corps 1, un chapeau 2, un siège 4 et un opercule 5. Les goujons 10 de liaison du corps et du chapeau sont également représentés. La Fig.7 montre un exemple de maillage de la vanne de la Fig.6 en vue d'un calcul thermique. Le calcul thermique est supposé être effectué en code ASTER de sorte que le maillage constitue un seul bloc en considérant qu'il y a continuité thermique entre le corps de vanne et le fluide. Si le calcul thermique se fait avec le code ASTER N3S Syrthès, on utilise deux maillages distincts. - pour le corps-chapeau, les discontinuités prises en compte en mécanique etant négligées du point de vue thermique, - pour le domaine fluide. La Fig.8 représente un exemple de maillage pour le calcul mécanique. Le maillage est compose des éléments distincts suivants - le corps de vanne 1, - le chapeau 2, la liaison corps-chapeau faisant intervenir des éléments de contact unilatéral, - l'opercule 5, la liaison opercule-corps faisant également intervenir des éléments de contact, - le groupe de goujons 10 : dans les modélisations déjà réalisées les goujons sont representés par les groupes de mailles pentaèdriques dont la section est egale à la section résistante des vrais goujons. Cette representation qui a semble la plus réaliste est susceptible d'être améliorée, notamment en représentant les écrous, ce qui assure une meilleure répartition de l'ef- fort de serrage. Les condition limites de liaisons entre corpschapeau et goujons sont décrites en reference à la Fig.2 La Fig.9 représente le détail du maillage du corps de vanne et de la goujonnerie. On notera la couronne d'élements en grisé de la portee d'étanchéité 16 sur lesquels l'opercule vient s'appuyer. Le calcul mécanique précédemment décrit vise initialement à déterminer les déplacements des noeuds internes de la structure. Ce calcul peut servir accessoirement - à analyser la contrainte dans le corps de la vanne pour le maximum de pression atteint et à vérifier si l'inté- grité mécanique est conservée, - à faire l'analyse de l'étanchéité au niveau de la liaison corps-chapeau. L'analyse des decollements entre le corps et le chapeau peut fournir un critère de fuite à une pression donnée. Le logiciel de post-traitement utilisé pour l'execu- tion des opérations sur les résultats des calculs thermiques et mécaniques exposés plus haut va maintenant être décrit en référence à l'organigramme de la Fig.lO. Le logiciel appelé ALADIN et écrit en FORTRAN permet de traiter la phase de post-traitement du procédé de l'invention. I1 permet 1 - de calculer, à partir des déplacements des noeuds de la peau interne d'une vanne, tels que les noeuds N1 à N3 de la Fig.l, la variation du volume interne de celle-ci sous l'effet d'un chargement en pression selon le principe décrit plus haut 2 - de calculer l'expansion volumique de la zone fluide pour un champ de température donné 3 - de déterminer, à partir des résultats des deux étapes précédentes, la pression atteinte pour un champ thermique donné. On detaille ci-après les fichiers d'entrée et de sortie du code, le format des données, et un exemple de chaque fichier est fourni par la suite. Une première étape consiste à relire les fichiers ther.resu et meca.resu pour un reformatage des données. Les fichiers.resu sont structurés de la façon suivante - noeuds, - mailles + connectivité, - groupes de noeuds, - groupes de mailles, - champs aux noeuds, (TEMP pour ther.resu et DX, DY et DZz pour meca.resu). I1 s'agit de créer des fichiers structurés pour pouvoir lire à chaque pas de temps - la connectivité de chaque maille, - les champs calculés aux noeuds de la maille. (cette structure est nécessaire pour faciliter les calculs ultérieurs), et ce uniquement sur les groupes intéressants. Ainsi - ther.resu est transformé en exemp.xtt qui ne contient que les données relatives au volume fluide, - meca.resu est transformé en exemp.xpp qui ne contient que les données relatives aux noeuds de la peau interne de la vanne. Nota : Les fichiers exemp.xtt et exemp.xpp sont temporaires. Ils sont détruits à l'issue du calcul. Les fichiers exemp.xtt contiennent - la liste de mailles et leur connectivité, suivie des coordonnées de chaque noeud nécessaires pour recalculer le volume de chaque tétraèdre, - la liste des mailles du volume fluide et la température moyenne à savoir la moyenne de la valeur calculée aux 4 noeuds sommets. Les fichiers exemp.xpp contiennent pour chaque numéro d'ordre : la liste des données relatives à chaque maille, soit des blocs constitués de la manière suivante - Maille + connectivité pour chaque noeud (3 en tout) - Nom du noeud, X, Y, Z dX dY dZ La suite du post-traitement consiste à : - Pour le fichier exemp.xtt : faire une boucle sur chaque maille M. Celle-ci subit par dilatation une expansion volumique où Ti est la température moyenne de la maille, a(T) le coefficient de dilatation, Vi son volume. La sommes des AVi est l'expansion volumique totale du domaine fluide. Ces résultats sont inscrits dans le fichier exemp.xtO. - Pour le fichier exemp.xpp : à chaque numéro d'ordre, on effectue la boucle suivante Les champs X, Y, Z, dX, dY, dZ des trois noeuds de chaque maille permettent de constituer le pentaèdre décrit en référence à la Fig.l et d'en calculer le volume. La somme de ces volumes élémentaires donne le AV total du volume interne de la vanne pour une pression donnée. Ces résultats sont inscrits dans le fichier exemp.xpO. Les fichiers xtO et xpO sont ensuite relus pour déterminer la pression atteinte pour chaque numéro d'ordre du changement thermique. Cette étape consiste à résoudre l'équation (expansion volumique du domaine fluide pour un champ T1 calculé à un numéro d'ordre i du calcul thermique). PO est la pression initiale dans le fluide X est la pression cherchée est la compressibilite de l'eau dv f(p) = où1 V est le volume interne du corps de la vanne Vdp est quasi-constant dans la plage # 0-1000 C # 0-100 bar donc On résout cette équation par une méthode itérative classique. V(p) (volume interne de la vanne en fonction de la pression) est calculé à l'aide d'un polynôme d'interpolation d'Hermite determiné à l'aide des valeurs issues du fichier exemp.xpO. Les resultats de chaque étape sont volontairement placés dans des fichiers séparés afin de pouvoir reprendre la procédure à n'importe quel niveau sans perdre les résultats des calculs précédents. Le programme peut être exécuté - soit en entrant les paramètres à l'écran, - soit à l'aide d'un fichier de commande tel qu'indiqué au bas de la Fig.10. (Aladon.dat). Enfin, un fichier appelé Aladin.bla permet une restitution détaillée du déroulement du programme. Le procédé qui vient d'être décrit permet donc, grâce à une modelisation appropriée, d'analyser par le calcul l'effet chaudière dans les vannes placées dans un environnement subissant d'importantes variations de température. REVENDICATIONS 1. Procédé de détermination de la pression maximale atteinte par le fluide contenu dans le volume défini par l'espace inter-opercules et le dôme corps/chapeau (1,2) d'un robinet vanne à double étanchéité, lors d'une montée de température déterminée de l'environnement extérieur au robinetvanne, caractérisé en ce qu'il consiste - à effectuer un calcul thermique ou thermohydraulique pour déterminer le champ de température du fluide contenu dans ledit volume lors de ladite montée en température, - à effectuer un calcul mécanique pour déterminer les déplacements de la surface intérieure du robinet-vanne définissant ledit volume pour des changements internes en pression évoluant dans une plage de pression contenant la pression maximum à déterminer, - à calculer à partir des résultats du calcul mécanique la variation dudit volume interne du robinet-vanne pour chaque valeur de pression dans ladite plage de pression, - à calculer à partir des résultats du calcul thermique ou thermohydraulique l'expansion volumique du fluide contenu dans ledit volume, et - à déterminer, à partir des résultats des calculs de la variation du volume interne d'une part et de l'expansion volumique du fluide contenu dans ledit volume, d'autre part ainsi que du fait que le fluide occupe tout ledit volume, la pression maximale atteinte pour ladite montée en température choisie. 2. Procede suivant la revendication 1, caractérisé en ce que les phases de calcul thermique ou thermohydraulique et de calcul mécanique sont des phases de calcul numérique. 3. Procédé suivant la revendication 2, caractérisé en ce que le calcul thermique ou thermohydraulique est un calcul numérique transitoire effectué sur un modèle qui comporte le corps de la vanne et le volume fluide qu'elle renferme. 4. Procédé suivant l'une des revendications 2 et 3 caractérisé en ce que le calcul mécanique est un calcul en éléments finis effectué uniquement sur le corps de la vanne. 5. Procédé suivant l'une des revendications 1 à 4, caractérisé en ce que par le calcul thermique on calcule pour chaque pas de temps, le champ de température dans le fluide pour un chargement en peau externe de la vanne correspondant à la montée en température de l'environnement et on calcule à partir de ce champ de température l'expansion volumique du fluide sous l'effet de la dilatation. 6. Procédé suivant l'une des revendications 1 à 5, caractérisé en ce que le calcul mécanique est effectué à température constante. 7. Procédé suivant la revendication 6, caractérisé en ce que la température constante à laquelle est effectuée le calcul mécanique est la température maximale de l'environnement. 8. Procédé suivant l'une des revendications 2 à 7, caractérisé en ce que la décomposition en éléments finis est une décomposition en tétraèdres de pentaèdres élémentaires (N1, N2, N3, Nl+dNl, N2+dN2, N3+ dN3) de la peau interne du corps de vanne à leur position initiale et des noeuds correspondants (N1+dNl, N2+dN2, N3+dN3) à leur position déplacée. 9.Procédé suivant l'une quelconque des revendications 1 à 8, caractérisé en ce que les calculs mécaniques et thermiques sont effectués avec le code ASTER. 10. Procédé suivant l'une des revendications 1 à 9, caractérisé en ce que les calculs de la variation du volume interne du robinet-vanne, de l'expansion volumique du fluide continu dans ledit volume la détermination à partir des résultats desdits calculs de la pression maximale atteinte par ladite montées en température, sont effectués par un posttraitement analytique des résultats desdits calculs thermiques et mécaniques. 11. Procédé suivant la revendication 10, caractérisé en ce que le post-traitement analytique est effectué au moyen d'un logiciel permettant - de calculer à partir du déplacement des noeuds (N1, N2, N3) de la peau interne de la vanne, la variation du volume interne de celle-ci sous l'effet d'un chargement en pression - de calculer l'expansion volumique de la zone fluide de la vanne pour un champ de température donné - de déterminer, à partir des résultats des deux étapes précédentes, la pression atteinte pour un champ thermique donné.