1 # -*- coding: utf-8 -*-
3 # --------------------------------------------------
4 # Definition de variables sous forme de tuple
5 # --------------------------------------------------
9 def __init__(self,ntuple):
12 def __convert__(self,valeur):
13 if type(valeur) == types.StringType:
15 if len(valeur) != self.ntuple:
20 return "Tuple de %s elements" % self.ntuple
25 # --------------------------------------------------
27 # --------------------------------------------------
32 class loi ( ASSD ) : pass
33 class variable ( ASSD ) : pass
37 JdC = JDC_CATA ( code = 'CUVE1D-DEFAILLGLOB',
39 regles = ( AU_MOINS_UN ('OPTIONS','DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'),
40 A_CLASSER ( ('OPTIONS'), ('DEFAUT') ),
41 A_CLASSER ( ('DEFAUT'), ('CUVE') ),
42 A_CLASSER ( ('CUVE'), ('MODELES') ),
43 A_CLASSER ( ('MODELES'), ('INITIALISATION') ),
44 A_CLASSER ( ('INITIALISATION'), ('REVETEMENT') ),
45 A_CLASSER ( ('REVETEMENT'), ('METAL_BASE') ),
46 A_CLASSER ( ('METAL_BASE'), ('TRANSITOIRE') )
50 # --------------------------------------------------
52 # --------------------------------------------------
54 # --------------------------------------------------
55 # RESPONSABLE D. THAI VAN
56 # Ce fichier contient la liste des coefficients pour un
57 # modele de Weibull generalise
58 def Coef_WeibGen() : return FACT(statut='o',min=1,max='**',
61 A1 = SIMP ( statut="o", typ="R", defaut=21.263,
62 fr="coef du coef a(T) d'une Weibull générale", ),
64 A2 = SIMP ( statut="o", typ="R", defaut=9.159,
65 fr="coef du coef a(T) d'une Weibull générale", ),
67 A3 = SIMP ( statut="o", typ="R", defaut=0.04057,
68 fr="coef du coef a(T) d'une Weibull générale", ),
70 B1 = SIMP ( statut="o", typ="R", defaut=17.153,
71 fr="coef du coef b(T) d'une Weibull générale", ),
73 B2 = SIMP ( statut="o", typ="R", defaut=55.089,
74 fr="coef du coef b(T) d'une Weibull générale", ),
76 B3 = SIMP ( statut="o", typ="R", defaut=0.0144,
77 fr="coef du coef b(T) d'une Weibull générale", ),
79 C1 = SIMP ( statut="o", typ="R", defaut=4.,
80 fr="coef du coef c(T) d'une Weibull générale", ),
82 C2 = SIMP ( statut="o", typ="R", defaut=0.,
83 fr="coef du coef c(T) d'une Weibull générale", ),
85 C3 = SIMP ( statut="o", typ="R", defaut=0.,
86 fr="coef du coef c(T) d'une Weibull générale", ),
88 ); # FIN def Coef_WeibGen
91 # --------------------------------------------------
92 # RESPONSABLE D. THAI VAN
93 # Ce fichier contient la liste des coefficients
94 def Coef_Fluence() : return FACT(statut='o',min=1,max='**',
97 Azimut_0deg = SIMP ( statut="o", typ="R", defaut=5.8,
98 fr="Fluence à l'azimut 0 (10^19 n/cm2)", ),
100 Azimut_5deg = SIMP ( statut="o", typ="R", defaut=5.48,
101 fr="Fluence à l'azimut 5 (10^19 n/cm2)", ),
103 Azimut_10deg = SIMP ( statut="o", typ="R", defaut=4.46,
104 fr="Fluence à l'azimut 10 (10^19 n/cm2)", ),
106 Azimut_15deg = SIMP ( statut="o", typ="R", defaut=3.41,
107 fr="Fluence à l'azimut 15 (10^19 n/cm2)", ),
109 Azimut_20deg = SIMP ( statut="o", typ="R", defaut=3.37,
110 fr="Fluence à l'azimut 20 (10^19 n/cm2)", ),
112 Azimut_25deg = SIMP ( statut="o", typ="R", defaut=3.16,
113 fr="Fluence à l'azimut 25 (10^19 n/cm2)", ),
115 Azimut_30deg = SIMP ( statut="o", typ="R", defaut=2.74,
116 fr="Fluence à l'azimut 30 (10^19 n/cm2)", ),
118 Azimut_35deg = SIMP ( statut="o", typ="R", defaut=2.25,
119 fr="Fluence à l'azimut 35 (10^19 n/cm2)", ),
121 Azimut_40deg = SIMP ( statut="o", typ="R", defaut=1.89,
122 fr="Fluence à l'azimut 40 (10^19 n/cm2)", ),
124 Azimut_45deg = SIMP ( statut="o", typ="R", defaut=1.78,
125 fr="Fluence à l'azimut 45 (10^19 n/cm2)", ),
127 ); # FIN def Coef_Fluence
129 #================================
130 # 1. Definition des OPTIONS
131 #================================
133 OPTIONS = PROC ( nom = "OPTIONS",
136 fr = "Définitions des options",
139 # Liste des paramètres
143 IncrementTemporel = SIMP ( statut="o", typ="I", max=1, defaut=1,
144 fr="Incrément temporel (=1 pour calcul déterministe)", ),
147 IncrementMaxTemperature = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1., defaut=0.1,
148 fr="Incrément maximum d'évolution de la température par noeud et par instant (°C)", ),
151 IncrementMaxTempsAffichage = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1000., defaut="1000.",
152 fr="Incrément maximum de temps pour l'affichage (s)", ),
155 # Question : NBO depend-il de TYPGEOM ??
156 NombreNoeudsMaillage = SIMP ( statut="o", typ="R", max=1, val_min=1., val_max=1000.,
157 fr = "Nombre de noeuds à considérer dans le maillage interne", ),
160 ListeInstants = SIMP ( statut="o", typ="R", max="**",
161 fr = "Liste des instants pour lesquels la température et les contraintes seront archivés", ),
165 #================================
166 # 2. Caracteristiques du DEFAUT
167 #================================
169 DEFAUT = PROC ( nom = "DEFAUT",
172 fr = "Caractéristiques du défaut",
175 # Liste des paramètres
179 TypeInitial = SIMP ( statut="o", typ="TXM",
180 fr="Type initial du défaut : sous revêtement ou débouchant",
181 into=( "Sous Revetement", # DSR
182 "Debouchant", ), # DD
186 # Definition des parametres selon le type du defaut
189 SousRevetement = BLOC ( condition = "TypeInitial=='Sous Revetement'",
191 # ORIEDEF into LONGITUD, CIRCONF
192 Orientation = SIMP ( statut="o", typ="TXM",
193 fr="Orientation du défaut",
194 into=( "Longitudinale",
195 "Circonferentielle" ), ),
198 ProfondeurRadiale = SIMP ( statut="o", typ="R", max=1, val_min=1., val_max=1., defaut=0.006,
199 fr="Profondeur radiale du défaut (m)", ),
201 # OPTLONG into VALEUR, RAPPORT
202 OptionCalculLongueur = SIMP ( statut="o", typ="TXM",
203 fr="Option pour caractériser la longueur du défaut : soit par valeur, soit par un rapport LONG/PROF",
205 "Relation lineaire avec la longueur" ), ),
207 Option_Valeur = BLOC ( condition = "OptionCalculLongueur=='Valeur'",
209 Longueur = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1., defaut = 0.060,
210 fr = "Longueur du défaut sous revêtement (m)", ),
211 ), # Fin BLOC Option_Valeur
213 Option_Rapport = BLOC ( condition = "OptionCalculLongueur=='Relation lineaire avec la longueur'",
215 LongueurSurProfondeur = SIMP ( statut="o", typ="R", max=1, val_max=100.,
216 fr="Rapport longueur/profondeur du défaut sous revêtement", ),
217 ), # Fin BLOC Option_Rapport
220 DecalageRadial = SIMP ( statut="o", typ="R", defaut=-0.00001,
221 fr="Décalage radial du defaut sous revêtement (en m)", ),
223 Azimut = SIMP ( statut="o", typ="R", defaut=0.,
224 fr="Position angulaire du défaut dans le cylindre de cuve (en degrés)", ),
226 Altitude = SIMP ( statut="o", typ="R", val_min=0., val_max=4., defaut=2.,
227 fr="Altitude du défaut sur le cylindre de cuve (en m)", ),
229 Pointe = SIMP ( statut="o", typ="TXM", defaut="A",
230 fr="Choix du point considéré du défaut sous revêtement",
231 into=( "A", "B" ), ),
233 CorrectionPlastiqueBeta = SIMP ( statut="o", typ="TXM", defaut="NON",
234 fr="Prise en compte de la correction plastique BETA dans le calcul du facteur d'intensité de contraintes",
235 into=( "OUI", "NON" ), ),
237 ), # Fin BLOC SousRevetement
239 Debouchant = BLOC ( condition = "TypeInitial=='Debouchant'",
241 # ORIEDEF into LONGITUD, CIRCONF
242 Orientation = SIMP ( statut="o", typ="TXM",
243 fr="Orientation du défaut : longitudinale ou circonférentielle",
244 into=( "Longitudinale",
245 "Circonferentielle" ), ),
247 ProfondeurRadiale = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1., defaut=0.006,
248 fr="Profondeur radiale du défaut (en m)", ),
250 Azimut = SIMP ( statut="o", typ="R", defaut=0.,
251 fr="Position angulaire du défaut dans le cylindre de la cuve (en degrés)", ),
253 Altitude = SIMP ( statut="o", typ="R", defaut=2.,
254 fr="Altitude du défaut sur le cylindre de la cuve (en m)", ),
256 CorrectionPlastiqueIrwin = SIMP ( statut="o", typ="TXM", defaut="NON",
257 fr="Prise en compte de la correction plastique d'Irwin dans le calcul du facteur d'intensité de contraintes",
258 into=( "OUI", "NON" ), ),
260 CorrectionPlastiqueBeta = SIMP ( statut="o", typ="TXM", defaut="NON",
261 fr="Prise en compte de la correction plastique BETA dans le calcul du facteur d'intensité de contraintes",
262 into=( "OUI", "NON" ), ),
264 ), # Fin BLOC debouchant
269 #================================
270 # 3. Caracteristiques de la CUVE
271 #================================
273 CUVE = PROC (nom = "CUVE",
276 fr = "Caractéristiques de la cuve",
279 # Liste des paramètres
283 TraitementGeometrie = SIMP ( statut="o", typ="TXM",
284 fr="Choix de la définition de la geométrie d'une cuve",
285 into=( "Topologie", # GEOMETRIE
286 "Maillage"), # MAILLAGE
290 # Definition des parametres selon le type de traitement de la geometrie
293 Geometrie = BLOC ( condition = "TraitementGeometrie=='Topologie'",
296 RayonInterne = SIMP ( statut="o", typ="R", val_min=0., defaut=1.994,
297 fr="Rayon interne de la cuve (en m)", ),
299 # DTV : comment preciser que REXT > RINT ?
301 RayonExterne = SIMP ( statut="o", typ="R", defaut=2.2015,
302 fr="Rayon externe de la cuve (en m)", ),
304 # DTV : comment preciser que LREV < RINT ?
306 EpaisseurRevetement = SIMP ( statut="o", typ="R", defaut=0.0075,
307 fr="Epaisseur du revêtement (m)", ),
310 LigamentExterneMin = SIMP ( statut="o", typ="R", defaut=0.75,
311 fr="Ligament externe minimal avant rupture (% de l'épaisseur de cuve)", ),
313 ), # Fin BLOC Geometrie
315 Maillage = BLOC ( condition = "TraitementGeometrie=='Maillage'",
317 # DTV : comment preciser que c'est une suite de liste de nombres croissants ?
318 # Question : quel rapport avec NBO ??
319 Liste_abscisses = SIMP ( statut="o", typ="R", max="**",
320 fr = "Liste des abscisses (m)", ),
321 ), # Fin BLOC Maillage
325 #====================================================
326 # 4. Modeles de fluence, d'irradiation et de tenacite
327 #====================================================
329 #=======================
330 # 4.1 Modeles de fluence
331 #=======================
333 MODELES = PROC ( nom = "MODELES",
336 fr = "Modèles de fluence, d'irradiation et de ténacité",
340 # Liste des paramètres
343 Fluence = FACT ( statut="o",
346 Modele = SIMP ( statut="o", typ="TXM",
347 fr="Modèle d'atténuation de la fluence dans l'épaisseur de la cuve",
348 into=( "Exponentiel sans revetement k=9.7 (Reglementaire)", # Reglementaire
349 "Exponentiel sans revetement k=12.7 (France)", # France
350 "Exponentiel sans revetement k=0. (ValeurImposee)", # ValeurImposee
351 "Donnees francaises du palier CPY (SDM)", # SDM
352 "Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)", # GrandeDev
353 "Regulatory Guide 1.99 rev 2 (USNRC)", # USNRC
354 "Dossier 900 MWe AP9701 rev 2 (REV_2)", # REV_2
355 "Lissage du modele ajuste (SDM_Lissage)", # SDM_Lissage
356 "Grand developpement (GD_Cuve)"), # GD_Cuve
360 # Definition des parametres selon le modele de fluence
363 Reglementaire = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=9.7 (Reglementaire)', ) ",
364 # DTV : comment proposer une liste de valeurs, tout en proposant de fournir d'autres valeurs ?
366 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
367 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
368 ), # Fin BLOC Reglementaire
370 France = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=12.7 (France)', ) ",
372 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
373 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
375 KPFrance = SIMP ( statut="o", typ="R",
376 fr="Paramètre exponentiel du modèle France",
380 ValeurImposee = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=0. (ValeurImposee)', ) ",
382 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
383 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
384 ), # Fin BLOC ValeurImposee
386 SDM = BLOC ( condition = " Modele in ( 'Donnees francaises du palier CPY (SDM)', ) ",
388 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
389 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
392 USNRC = BLOC ( condition = " Modele in ( 'Regulatory Guide 1.99 rev 2 (USNRC)', ) ",
394 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
395 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
396 KPUS = SIMP ( statut="o", typ="R",
397 fr="Paramètre exponentiel du modèle US",
401 REV_2 = BLOC ( condition = " Modele in ( 'Dossier 900 MWe AP9701 rev 2 (REV_2)', ) ",
403 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
404 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
407 SDM_Lissage = BLOC ( condition = " Modele in ( 'Lissage du modele ajuste (SDM_Lissage)', ) ",
409 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
410 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
411 ), # Fin BLOC SDM_Lissage
413 GrandeDev = BLOC ( condition = " Modele in ( 'Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)', ) ",
415 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
416 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
417 ), # Fin BLOC GrandeDev
419 GD_Cuve = BLOC ( condition = " Modele in ( 'Grand developpement (GD_Cuve)', ) ",
421 FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5,
422 fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ),
423 Coefficients = Coef_Fluence(),
424 ), # Fin BLOC GD_Cuve
426 ), # Fin FACT Fluence
428 #==========================
429 # 4.2 Modeles d'irradiation
430 #==========================
432 Irradiation = FACT ( statut="o",
435 Option = SIMP ( statut = "o", typ = "TXM",
436 fr = "Type d'irradiation",
437 into = ( "RTndt de la cuve a l instant de l analyse", # RTNDT
438 "Modele d irradiation" ), # FLUENCE
442 # Definition des parametres selon le type d'irradiation
445 RTndt = BLOC ( condition = "Option=='RTndt de la cuve a l instant de l analyse'",
447 RTNDT = SIMP ( statut="o", typ="R", defaut=73.,
448 fr="RTNDT de la cuve à l'instant de l'analyse (°C)", ),
452 ModeleIrradiation = BLOC ( condition = "Option=='Modele d irradiation'",
455 Modele = SIMP ( statut="o", typ="TXM",
456 fr="Modèle d'irradiation pour virole ou joint soudé",
457 into=( "Metal de Base : formule de FIM/FIS Houssin", # HOUSSIN
458 "Metal de Base : formule de FIM/FIS Persoz", # PERSOZ
459 "Metal de Base : formule de FIM/FIS Lefebvre", # LEFEBVRE
460 "Metal de Base : Regulatory Guide 1.00 rev 2", # USNRCmdb
461 "Joint Soude : formulation de FIM/FIS Brillaud", # BRILLAUD
462 "Joint Soude : Regulatory Guide 1.00 rev 2" ), # USNRCsoud
465 TeneurCuivre = SIMP ( statut="o", typ="R", defaut=0.,
466 fr="Teneur en cuivre (%)", ),
468 TeneurNickel = SIMP ( statut="o", typ="R", defaut=0.,
469 fr="Teneur en nickel (%)", ),
471 TeneurPhosphore = SIMP ( statut="o", typ="R", defaut=0.,
472 fr="Teneur en phosphore (%)", ),
474 Parametres_USNRC = BLOC ( condition = " Modele in ( 'Joint Soude : Regulatory Guide 1.00 rev 2', 'Metal de Base : Regulatory Guide 1.00 rev 2' , ) ",
476 MoyenneRTndt = SIMP ( statut="o", typ="R", defaut=0.,
477 fr="Moyenne de RTNDT : virole C1 de cuve Chinon : mdb=>-17.°C et js=>42.°C (HT-56/05/038 : p.52)", ),
479 CoefVariationRTndt = SIMP ( statut="o", typ="R", defaut=0.,
480 fr="Coefficient de variation de la RTNDT initiale", ),
482 EcartTypeRTndt = SIMP ( statut="o", typ="R", defaut=28.,
483 fr="Ecart-type du décalage de RTNDT (°F) (28. pour js et 17. pour mdb)", ),
485 NombreEcartTypeRTndt = SIMP ( statut="o", typ="R", defaut=2.,
486 fr="Nombre d'écart-type par rapport à la moyenne de DRTNDT", ),
487 ), # Fin BLOC Parametres_USNRC
489 ), # Fin BLOC ModeleIrradiation
491 ), # Fin FACT Irradiation
493 #========================
494 # 4.3 Modeles de tenacite
495 #========================
497 Tenacite = FACT ( statut = "o",
500 Modele = SIMP ( statut="o", typ="TXM",
501 fr="Modèle de calcul de la ténacité à l'amorçage KIc",
502 into=( "RCC-M/ASME coefficient=2", # RCC-M
503 "RCC-M/ASME coefficient=2.33 (Houssin)", # Houssin_RC
504 "RCC-M/ASME avec KI=KIpalier", # RCC-M_pal
505 "RCC-M/ASME avec KI~exponentiel", # RCC-M_exp
506 "Weibull basee sur la master cuve", # Wallin
507 "Weibull basee sur la master cuve (REME)", # REME
508 "Weibull n°1 (etude ORNL)", # ORNL
509 "Weibull n°2", # WEIB2
510 "Weibull n°3", # WEIB3
511 "Weibull generalisee", # WEIB_GEN
512 "Exponentielle n°1 (Frama)", # Frama
513 "Exponentielle n°2 (LOGWOLF)" ), # LOGWOLF
518 # Definition des parametres selon le modele de tenacité
523 KIc_RCCM = BLOC ( condition = " Modele in ( 'RCC-M/ASME coefficient=2', 'RCC-M/ASME coefficient=2.33 (Houssin)', 'RCC-M/ASME avec KI=KIpalier', 'RCC-M/ASME avec KI~exponentiel', ) ",
526 NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2.,
527 fr = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ),
530 PalierDuctile_KIc = SIMP ( statut="o", typ="R", defaut=195. ,
531 fr="Palier déterministe de K1c ou valeur du palier ductile plafonnant la courbe (en MPa(m^0.5)) ", ),
534 CoefficientVariation_KIc = SIMP ( statut="o", typ="R", defaut = 0.15,
535 fr = "Coefficient de variation de la loi normale de K1c ", ),
537 ), # Fin BLOC KIc_RCCM
539 Fissure = BLOC ( condition = " Modele in ( 'RCC-M/ASME coefficient=2', 'RCC-M/ASME coefficient=2.33 (Houssin)', 'RCC-M/ASME avec KI=KIpalier', 'RCC-M/ASME avec KI~exponentiel', )",
542 ArretDeFissure = SIMP ( statut="o", typ="TXM", defaut="NON",
543 fr="Prise en compte de l'arrêt de fissure",
544 into=( "OUI", "NON" ), ),
546 Parametres_Fissure = BLOC ( condition = "ArretDeFissure=='OUI'",
548 IncrementTailleFissure = SIMP ( statut="o", typ="R", defaut=0.005,
549 fr="Incrément de la taille de fissure pour la propagation (en m)", ),
550 ), # Fin BLOC Parametres_Fissure
552 KIa_RCCM = BLOC ( condition = "ArretDeFissure=='OUI'",
555 NbEcartType_MoyKIa = SIMP ( statut="o", typ="R", defaut=-2.,
556 fr="Nombre d'écart-type par rapport à la moyenne de KIa (nb sigma) ", ),
559 CoefficientVariation_KIa = SIMP ( statut="o", typ="R", defaut=0.10,
560 fr="Coefficient de variation de la loi normale de K1a -ténacite à l'arrêt- ", ),
563 PalierDuctile_KIa = SIMP ( statut="o", typ="R", defaut=195.,
564 fr="Palier déterministe de K1a -ténacite à l'arrêt- (en MPa(m^0.5)) ", ),
565 ), # Fin BLOC KIa_RCCM
567 ), # Fin BLOC Fissure
569 # Modeles type Weibull
571 KIc_MasterCuve = BLOC ( condition = " Modele in ( 'Weibull basee sur la master cuve', ) ",
574 Temperature_KIc100 = SIMP ( statut="o", typ="I", defaut=-27,
575 fr="Paramètre T0 du modèle Wallin (°C) : température pour laquelle la téncité du matériau vaut en moyenne 100MPa.m^5", ),
578 Fractile_KIc = SIMP ( statut="o", typ="R", defaut=5.,
579 fr="Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", ),
581 ), # Fin BLOC KIc_MasterCuve
583 KIc_Weibull1 = BLOC ( condition = " Modele in ( 'Weibull basee sur la master cuve (REME)', 'Weibull n°1 (etude ORNL)', ) ",
586 Fractile_KIc = SIMP ( statut="o", typ="R", defaut=5.,
587 fr="Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", ),
589 ), # Fin BLOC KIc_Weibull1
591 KIc_Weibull2 = BLOC ( condition = " Modele in ( 'Weibull n°2', 'Weibull n°3', ) ",
594 NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2.,
595 fr="Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ),
597 ), # Fin BLOC KIc_Weibull2
599 Weibull_Generalisee = BLOC ( condition = " Modele in ( 'Weibull generalisee',) ",
601 Coefficients = Coef_WeibGen(),
603 NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2.,
604 fr="Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ),
606 ), # Fin BLOC Weibull_Generalisee
608 # Modeles type exponentiel (Frama, LOGWOLF)
610 KIc_Exponentielle = BLOC ( condition = " Modele in ( 'Exponentielle n°1 (Frama)', 'Exponentielle n°2 (LOGWOLF)', ) ",
613 NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2.,
614 fr="Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ),
616 ), # Fin BLOC KIc_Exponentielle
618 ), # Fin FACT Tenacite
627 INITIALISATION = PROC ( nom = "INITIALISATION",
630 fr = "Initialisation : instant initial, profils radiaux de température et contraintes",
632 TemperatureInitiale = FACT ( statut = "o",
634 ProfilRadial_TemperatureInitiale = SIMP ( statut="o", typ=Tuple(2), max="**",
635 fr="Profil radial de la température initiale dans la cuve (en m : °C) ", ),
636 Amont_TemperatureInitiale = SIMP ( statut="o", typ="TXM",
637 fr="Prolongation à la frontière amont",
638 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
639 Aval_TemperatureInitiale = SIMP ( statut="o", typ="TXM",
640 fr="Prolongation à la frontière aval",
641 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
643 ), # Fin FACT TemperatureInitiale
645 ContraintesInitiales = FACT ( statut = "o",
647 ProfilRadial_ContraintesInitiales = SIMP ( statut="o", typ=Tuple(4), max="**",
648 fr="Profil radial des contraintes radiale, circonférentielle et longitudinale dans la cuve (en m : xx : xx : xx) ", ),
649 Amont_ContraintesInitiales = SIMP ( statut="o", typ="TXM",
650 fr="Prolongation à la frontière amont",
651 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
652 Aval_ContraintesInitiales = SIMP ( statut="o", typ="TXM",
653 fr="Prolongation à la frontière aval",
654 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
656 ), # Fin FACT ContraintesInitiales
659 InstantInitialisation = SIMP ( statut="o", typ="R", defaut = -1.,
660 fr="Instant initial auquel sont définies la température, ainsi que les contraintes initiales (en s) ", ),
662 ) # Fin PROC INITIALISATION
665 #==================================
666 # 6. CARACTERISTIQUES DU REVETEMENT
667 #==================================
669 REVETEMENT = PROC ( nom = "REVETEMENT",
672 fr = "Caracteristiques du revêtement",
675 ConditionLimiteThermiqueREV = SIMP ( statut="o", typ="TXM",
676 fr="Option pour définir les caractéristiques du revêtement ",
677 into=( "ENTHALPIE", "CHALEUR",),
680 EnthalpieREV = BLOC ( condition = "ConditionLimiteThermiqueREV=='ENTHALPIE'",
682 EnthalpieREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
683 fr="Température (°C) / enthalpie massique (J/kg) ", ),
684 Amont_EnthalpieREV = SIMP ( statut="o", typ="TXM",
685 fr="Prolongation à la frontière amont",
686 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
687 Aval_EnthalpieREV = SIMP ( statut="o", typ="TXM",
688 fr="Prolongation à la frontière aval",
689 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
691 ), # Fin BLOC EnthalpieREV
694 ChaleurREV = BLOC ( condition = "ConditionLimiteThermiqueREV=='CHALEUR'",
696 ChaleurREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
697 fr="Température (°C) / chaleur volumique (J/kg/K) ", ),
698 Amont_ChaleurREV = SIMP ( statut="o", typ="TXM",
699 fr="Prolongation à la frontière amont",
700 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
701 Aval_ChaleurREV = SIMP ( statut="o", typ="TXM",
702 fr="Prolongation à la frontière aval",
703 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
705 ), # Fin BLOC ChaleurREV
707 ConductiviteREV = FACT (statut = "o",
709 ConductiviteREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
710 fr="Température (°C) / conductivité thermique (W/m/°C) ", ),
711 Amont_ConductiviteREV = SIMP ( statut="o", typ="TXM",
712 fr="Prolongation à la frontière amont",
713 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
714 Aval_ConductiviteREV = SIMP ( statut="o", typ="TXM",
715 fr="Prolongation à la frontière aval",
716 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
718 ), # Fin FACT ConductiviteREV
720 ModuleYoungREV = FACT (statut = "o",
722 ModuleYoungREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
723 fr="Température (°C) / module d'Young (MPa) ", ),
724 Amont_ModuleYoungREV = SIMP ( statut="o", typ="TXM",
725 fr="Prolongation à la frontière amont",
726 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
727 Aval_ModuleYoungREV = SIMP ( statut="o", typ="TXM",
728 fr="Prolongation à la frontière aval",
729 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
731 ), # Fin FACT ModuleYoungREV
733 CoeffDilatThermREV = FACT (statut = "o",
735 CoeffDilatThermREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
736 fr="Température (°C) / coefficient de dilatation thermique (°C-1) ", ),
737 Amont_CoeffDilatThermREV = SIMP ( statut="o", typ="TXM",
738 fr="Prolongation à la frontière amont",
739 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
740 Aval_CoeffDilatThermREV = SIMP ( statut="o", typ="TXM",
741 fr="Prolongation à la frontière aval",
742 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
744 ), # Fin FACT CoeffDilatThermREV
746 LimiteElasticiteREV = FACT (statut = "o",
748 LimiteElasticiteREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
749 fr="Température (°C) / limite d'élasticite (MPa) ", ),
750 Amont_LimiteElasticiteREV = SIMP ( statut="o", typ="TXM",
751 fr="Prolongation à la frontière amont",
752 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
753 Aval_LimiteElasticiteREV = SIMP ( statut="o", typ="TXM",
754 fr="Prolongation à la frontière aval",
755 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
757 ), # Fin FACT LimiteElasticiteREV
759 AutresParametresREV = FACT (statut = "o",
762 TemperatureDeformationNulleREV = SIMP ( statut="o", typ="R", defaut=20.,
763 fr="Température de référence pour laquelle les déformations thermiques sont nulles (°C) ", ),
765 TemperaturePourCoefDilatThermREV = SIMP ( statut="o", typ="R", defaut=287.,
766 fr="Température de définition du coefficient de dilatation thermique (°C) ", ),
768 CoefficientPoissonREV = SIMP ( statut="o", typ="R", defaut=0.3,
769 fr="Coefficient de Poisson ", ),
771 ), # Fin FACT AutresParametresREV
773 ) # Fin PROC REVETEMENT
776 #=====================================
777 # 7. CARACTERISTIQUES DU METAL DE BASE
778 #=====================================
780 METAL_BASE = PROC ( nom = "METAL_BASE",
783 fr = "Caracteristiques du metal de base",
786 ConditionLimiteThermiqueMDB = SIMP ( statut="o", typ="TXM",
787 fr="Option pour definir les caractéristiques du revêtement ",
788 into=( "ENTHALPIE", "CHALEUR",), ),
790 EnthalpieMDB = BLOC ( condition = "ConditionLimiteThermiqueMDB=='ENTHALPIE'",
792 EnthalpieMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
793 fr="Température (°C) / enthalpie massique (J/kg) ", ),
794 Amont_EnthalpieMDB = SIMP ( statut="o", typ="TXM",
795 fr="Prolongation à la frontière amont",
796 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
797 Aval_EnthalpieMDB = SIMP ( statut="o", typ="TXM",
798 fr="Prolongation à la frontière aval",
799 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
801 ), # Fin BLOC EnthalpieMDB
803 ChaleurMDB = BLOC ( condition = "ConditionLimiteThermiqueMDB=='CHALEUR'",
805 ChaleurMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
806 fr="Température (°C) / chaleur volumique (J/kg/K) ", ),
807 Amont_ChaleurMDB = SIMP ( statut="o", typ="TXM",
808 fr="Prolongation à la frontière amont",
809 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
810 Aval_ChaleurMDB = SIMP ( statut="o", typ="TXM",
811 fr="Prolongation à la frontière aval",
812 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
814 ), # Fin BLOC ChaleurMDB
816 ConductiviteMDB = FACT ( statut = "o",
818 ConductiviteMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
819 fr="Température (°C) / conductivité thermique (W/m/°C) ", ),
820 Amont_ConductiviteMDB = SIMP ( statut="o", typ="TXM",
821 fr="Prolongation à la frontière amont",
822 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
823 Aval_ConductiviteMDB = SIMP ( statut="o", typ="TXM",
824 fr="Prolongation à la frontière aval",
825 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
827 ), # Fin FACT ConductiviteMDB
829 ModuleYoungMDB = FACT ( statut="o",
831 ModuleYoungMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
832 fr="Température (°C) / module d'Young (MPa) ", ),
833 Amont_ModuleYoungMDB = SIMP ( statut="o", typ="TXM",
834 fr="Prolongation à la frontière amont",
835 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
836 Aval_ModuleYoungMDB = SIMP ( statut="o", typ="TXM",
837 fr="Prolongation à la frontière aval",
838 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
840 ), # Fin FACT ModuleYoungMDB
842 CoeffDilatThermMDB = FACT ( statut="o",
844 CoeffDilatThermMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**",
845 fr="Température (°C) / coefficient de dilatation thermique (°C-1) ", ),
846 Amont_CoeffDilatThermMDB = SIMP ( statut="o", typ="TXM",
847 fr="Prolongation à la frontière amont",
848 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
849 Aval_CoeffDilatThermMDB = SIMP ( statut="o", typ="TXM",
850 fr="Prolongation à la frontière aval",
851 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
853 ), # Fin FACT CoeffDilatThermMDB
855 AutresParametresMDB = FACT ( statut = "o",
858 TemperatureDeformationNulleMDB = SIMP ( statut="o", typ="R", defaut=20.,
859 fr="Température de référence pour laquelle les déformations thermiques sont nulles (°C) ", ),
861 TemperaturePourCoefDilatThermMDB = SIMP ( statut="o", typ="R", defaut=287.,
862 fr="Température de définition du coefficient de dilatation thermique (°C) ", ),
864 CoefficientPoissonMDB = SIMP ( statut="o", typ="R", defaut=0.3,
865 fr="Coefficient de Poisson ", ),
867 ), # Fin FACT TemperatureDeformationNulleMDB
869 ) # Fin PROC METAL_BASE
872 #===============================
873 # 8. TRANSITOIRE THERMOMECANIQUE
874 #===============================
876 TRANSITOIRE = PROC ( nom = "TRANSITOIRE",
879 fr = "Description du transitoire thermohydraulique",
881 Pression = FACT ( statut = "o",
883 ProfilTemporel_Pression = SIMP ( statut="o", typ=Tuple(2), max="**",
884 fr = "Instant (s) / pression (MPa) ", ),
885 Amont_Pression = SIMP ( statut="o", typ="TXM",
886 fr="Prolongation à la frontière amont",
887 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
888 Aval_Pression = SIMP ( statut="o", typ="TXM",
889 fr="Prolongation à la frontière aval",
890 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
892 ), # FIN FACT Pression
895 TypeConditionLimiteThermique = SIMP ( statut="o", typ="TXM",
896 fr="Type de condition thermique en paroi interne ",
897 into=( "Temperature imposee en paroi", # TEMP_IMPO
898 "Flux de chaleur impose en paroi", # FLUX_REP
899 "Temperature imposee du fluide et coefficient echange", # ECHANGE
900 "Debit massique et temperature d injection de securite", # DEBIT
901 "Temperature imposee du fluide et debit d injection de securite"), # TEMP_FLU
904 TemperatureImposeeParoi = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee en paroi', ) ",
906 ProfilTemporel_TemperatureImposeeParoi = SIMP ( statut="o", typ=Tuple(2), max="**",
907 fr = "Instant (s) / Température imposée (°C) ", ),
908 Amont_TemperatureImposeeParoi = SIMP ( statut="o", typ="TXM",
909 fr="Prolongation à la frontière amont",
910 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
911 Aval_TemperatureImposeeParoi = SIMP ( statut="o", typ="TXM",
912 fr="Prolongation à la frontière aval",
913 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
915 ), # Fin BLOC TemperatureImposeeParoi
917 FluxChaleur = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Flux de chaleur impose en paroi', ) ",
919 ProfilTemporel_FluxChaleuri = SIMP ( statut="o", typ=Tuple(2), max="**",
920 fr="Instant (s) / Flux de chaleur impose (W/m2) ", ),
921 Amont_FluxChaleur = SIMP ( statut="o", typ="TXM",
922 fr="Prolongation à la frontière amont",
923 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
924 Aval_FluxChaleur = SIMP ( statut="o", typ="TXM",
925 fr="Prolongation à la frontière aval",
926 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
928 ), # Fin BLOC FluxChaleur
930 TemperatureImposeeFluide = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', 'Temperature imposee du fluide et debit d injection de securite', ) ",
932 ProfilTemporel_TemperatureImposeeFluide = SIMP ( statut="o", typ=Tuple(2), max="**",
933 fr = "Instant (s) / Température imposée (°C) ", ),
934 Amont_TemperatureImposeeFluide = SIMP ( statut="o", typ="TXM",
935 fr="Prolongation à la frontière amont",
936 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
937 Aval_TemperatureImposeeFluide = SIMP ( statut="o", typ="TXM",
938 fr="Prolongation à la frontière aval",
939 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
941 ), # Fin BLOC TemperatureImposeeFluide
943 CoefficientEchange = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', ) ",
945 ProfilTemporel_CoefficientEchange = SIMP ( statut="o", typ=Tuple(2), max="**",
946 fr="Instant (s) / Coefficient d'échange (W/m2/K) ", ),
947 Amont_CoefficientEchange = SIMP ( statut="o", typ="TXM",
948 fr="Prolongation à la frontière amont",
949 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
950 Aval_CoefficientEchange = SIMP ( statut="o", typ="TXM",
951 fr="Prolongation à la frontière aval",
952 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
954 ), # Fin BLOC CoefficientEchange
956 DebitMassique = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ",
958 ProfilTemporel_DebitMassique = SIMP ( statut="o", typ=Tuple(2), max="**",
959 fr="Instant (s) / Débit massique (kg/s) ", ),
960 Amont_DebitMassique = SIMP ( statut="o", typ="TXM",
961 fr="Prolongation à la frontière amont",
962 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
963 Aval_DebitMassique = SIMP ( statut="o", typ="TXM",
964 fr="Prolongation à la frontière aval",
965 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
967 ), # Fin BLOC DebitMassique
969 TemperatureInjection = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ",
971 ProfilTemporel_TemperatureInjection = SIMP ( statut="o", typ=Tuple(2), max="**",
972 fr="Instant (s) / Température d'injection de sécurité (°C) ", ),
973 Amont_TemperatureInjection = SIMP ( statut="o", typ="TXM",
974 fr="Prolongation à la frontière amont",
975 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
976 Aval_TemperatureInjection = SIMP ( statut="o", typ = "TXM",
977 fr="Prolongation à la frontière aval",
978 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
979 ), # Fin BLOC TemperatureInjection
981 Creare = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ",
984 DiametreHydraulique = SIMP ( statut="o", typ="R", defaut=-2.,
985 fr="Diamètre hydraulique (m) ", ),
987 SectionEspaceAnnulaire = SIMP ( statut="o", typ="R", defaut=-2.,
988 fr="Section espace annulaire (m2) ", ),
990 HauteurCaracConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=-2.,
991 fr="Hauteur caractéristique convection naturelle (m) ", ),
993 CoeffVestale_ConvectionForcee = SIMP ( statut="o", typ="R", defaut=1.,
994 fr="Coefficient Vestale convection forcée (-) ", ),
996 CoeffVestale_ConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=1.,
997 fr="Coefficient Vestale convection naturelle (-) ", ),
999 CritereConvergenceRelative = SIMP ( statut="o", typ="R", defaut=0.00001,
1000 fr = "Critère convergence relative (-) ", ),
1002 VolumeMelange_CREARE = SIMP ( statut="o", typ="R", defaut=-2.,
1003 fr = "Volume de mélange CREARE (m3) ", ),
1005 TemperatureInitiale_CREARE = SIMP ( statut="o", typ="R", defaut=-2.,
1006 fr="Température initiale CREARE (°C) ", ),
1008 SurfaceEchange_FluideStructure = SIMP ( statut="o", typ="R", defaut=-2.,
1009 fr="Surface d'échange fluide/structure (m2) ", ),
1011 ), # Fin BLOC Creare
1013 DebitInjection = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ",
1015 ProfilTemporel_DebitInjection = SIMP ( statut="o", typ=Tuple(2), max="**",
1016 fr="Instant (s) / Débit d'injection de sécurité (kg/s) ", ),
1017 Amont_DebitInjection = SIMP ( statut="o", typ="TXM",
1018 fr="Prolongation à la frontière amont",
1019 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
1020 Aval_DebitInjection = SIMP ( statut="o", typ="TXM",
1021 fr="Prolongation à la frontière aval",
1022 into=( 'Continu', 'Exclu', 'Lineaire' ), ),
1024 ), # Fin BLOC DebitInjection
1027 Vestale = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ",
1030 DiametreHydraulique = SIMP ( statut="o", typ="R", defaut=-2.,
1031 fr="Diamètre hydraulique (m) ", ),
1033 SectionEspaceAnnulaire = SIMP ( statut="o", typ="R", defaut=-2.,
1034 fr="Section espace annulaire (m2) ", ),
1036 HauteurCaracConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=-2.,
1037 fr="Hauteur caractéristique convection naturelle (m) ", ),
1039 CoeffVestale_ConvectionForcee = SIMP ( statut="o", typ="R", defaut=1.,
1040 fr="Coefficient d'échange Vestale convection forcée (-) ", ),
1042 CoeffVestale_ConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=1.,
1043 fr="Coefficient d'échange Vestale convection naturelle (-) ", ),
1045 CritereConvergenceRelative = SIMP ( statut="o", typ="R", defaut=0.00001,
1046 fr="Critère convergence relative (-) ", ),
1048 ), # Fin BLOC Vestale
1050 ) # Fin PROC TRANSITOIRE