]> SALOME platform Git repositories - tools/eficas.git/blob - Cuve2dg/Cuve2dg_Cata_V3.py
Salome HOME
Modif V6_4_°
[tools/eficas.git] / Cuve2dg / Cuve2dg_Cata_V3.py
1 # -*- coding: utf-8 -*-
2
3 # --------------------------------------------------
4 # Definition de variables sous forme de tuple
5 # --------------------------------------------------
6
7 import types
8 class Tuple:
9   def __init__(self,ntuple):
10     self.ntuple=ntuple
11
12   def __convert__(self,valeur):
13     if type(valeur) == types.StringType:
14       return None
15     if len(valeur) != self.ntuple:
16       return None
17     return valeur
18
19   def info(self):
20     return "Tuple de %s elements" % self.ntuple
21
22   __repr__=info
23   __str__=info
24
25 # --------------------------------------------------
26 # debut entete
27 # --------------------------------------------------
28
29 import Accas
30 from Accas import *
31
32 class loi ( ASSD ) : pass
33 class variable ( ASSD ) : pass
34
35
36 #CONTEXT.debug = 1
37 JdC = JDC_CATA ( code = 'CUVE1D-DEFAILLGLOB',
38                  execmodul = None,
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') )
47                           )
48                  ) # Fin JDC_CATA
49
50 # --------------------------------------------------
51 # fin entete
52 # --------------------------------------------------
53
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='**',
59
60   # A1
61   A1 = SIMP ( statut="o", typ="R", defaut=21.263, 
62               fr="coef du coef a(T) d'une Weibull générale", ),
63   # A2
64   A2 = SIMP ( statut="o", typ="R", defaut=9.159, 
65               fr="coef du coef a(T) d'une Weibull générale", ),
66   # A3
67   A3 = SIMP ( statut="o", typ="R", defaut=0.04057, 
68               fr="coef du coef a(T) d'une Weibull générale", ),
69   # B1
70   B1 = SIMP ( statut="o", typ="R", defaut=17.153, 
71               fr="coef du coef b(T) d'une Weibull générale", ),
72   # B2
73   B2 = SIMP ( statut="o", typ="R", defaut=55.089, 
74               fr="coef du coef b(T) d'une Weibull générale", ),
75   # B3
76   B3 = SIMP ( statut="o", typ="R", defaut=0.0144, 
77               fr="coef du coef b(T) d'une Weibull générale", ),
78   # C1
79   C1 = SIMP ( statut="o", typ="R", defaut=4., 
80               fr="coef du coef c(T) d'une Weibull générale", ),
81   # C2
82   C2 = SIMP ( statut="o", typ="R", defaut=0., 
83               fr="coef du coef c(T) d'une Weibull générale", ),
84   # C3
85   C3 = SIMP ( statut="o", typ="R", defaut=0., 
86               fr="coef du coef c(T) d'une Weibull générale", ),
87
88 ); # FIN def Coef_WeibGen
89
90
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='**',
95
96   # COEFFLUENCE1
97   Azimut_0deg  = SIMP ( statut="o", typ="R", defaut=5.8, 
98                         fr="Fluence à l'azimut 0 (10^19 n/cm2)", ),
99   # COEFFLUENCE2
100   Azimut_5deg  = SIMP ( statut="o", typ="R", defaut=5.48, 
101                         fr="Fluence à l'azimut 5 (10^19 n/cm2)", ),
102   # COEFFLUENCE3
103   Azimut_10deg = SIMP ( statut="o", typ="R", defaut=4.46, 
104                         fr="Fluence à l'azimut 10 (10^19 n/cm2)", ),
105   # COEFFLUENCE4
106   Azimut_15deg = SIMP ( statut="o", typ="R", defaut=3.41, 
107                         fr="Fluence à l'azimut 15 (10^19 n/cm2)", ),
108   # COEFFLUENCE5
109   Azimut_20deg = SIMP ( statut="o", typ="R", defaut=3.37, 
110                         fr="Fluence à l'azimut 20 (10^19 n/cm2)", ),
111   # COEFFLUENCE6
112   Azimut_25deg = SIMP ( statut="o", typ="R", defaut=3.16, 
113                         fr="Fluence à l'azimut 25 (10^19 n/cm2)", ),
114   # COEFFLUENCE7
115   Azimut_30deg = SIMP ( statut="o", typ="R", defaut=2.74, 
116                         fr="Fluence à l'azimut 30 (10^19 n/cm2)", ),
117   # COEFFLUENCE8
118   Azimut_35deg = SIMP ( statut="o", typ="R", defaut=2.25, 
119                         fr="Fluence à l'azimut 35 (10^19 n/cm2)", ),
120   # COEFFLUENCE9
121   Azimut_40deg = SIMP ( statut="o", typ="R", defaut=1.89, 
122                         fr="Fluence à l'azimut 40 (10^19 n/cm2)", ),
123   # COEFFLUENCE10
124   Azimut_45deg = SIMP ( statut="o", typ="R", defaut=1.78, 
125                         fr="Fluence à l'azimut 45 (10^19 n/cm2)", ),
126
127 ); # FIN def Coef_Fluence
128
129 #================================
130 # 1. Definition des OPTIONS
131 #================================
132
133 OPTIONS = PROC ( nom = "OPTIONS",
134                  op = 68,
135                  repetable = 'n',
136                  fr = "Définitions des options", 
137
138 #===
139 # Liste des paramètres
140 #===
141
142 # INCRTPS
143   IncrementTemporel          = SIMP ( statut="o", typ="I", max=1, defaut=1, 
144                                       fr="Incrément temporel (=1 pour calcul déterministe)", ),
145
146 # DTPREC
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)", ),
149
150 # DTARCH
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)", ),
153
154 # NBO
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", ),
158
159
160   ListeInstants              = SIMP ( statut="o", typ="R", max="**",
161                                       fr = "Liste des instants pour lesquels la température et les contraintes seront archivés", ),
162
163 ) # Fin PROC OPTIONS
164
165 #================================
166 # 2. Caracteristiques du DEFAUT
167 #================================
168
169 DEFAUT = PROC ( nom = "DEFAUT",
170                 op = 68,
171                 repetable = 'n',
172                 fr = "Caractéristiques du défaut", 
173
174 #===
175 # Liste des paramètres
176 #===
177
178 # TYPEDEF
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
183                        ),
184
185 #====
186 # Definition des parametres selon le type du defaut
187 #====
188
189   SousRevetement = BLOC ( condition = "TypeInitial=='Sous Revetement'",
190
191     # ORIEDEF into LONGITUD, CIRCONF
192     Orientation              = SIMP ( statut="o", typ="TXM", 
193                                       fr="Orientation du défaut",
194                                       into=( "Longitudinale", 
195                                              "Circonferentielle" ), ),
196
197     # PROFDEF
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)", ),
200
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",
204                                       into = ( "Valeur", 
205                                                "Relation lineaire avec la longueur" ), ),
206
207     Option_Valeur            = BLOC ( condition = "OptionCalculLongueur=='Valeur'",
208       # LONGDEF
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
212
213     Option_Rapport           = BLOC ( condition = "OptionCalculLongueur=='Relation lineaire avec la longueur'",
214       # LONGSURPROF
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
218
219     # DECADEF
220     DecalageRadial           = SIMP ( statut="o", typ="R", defaut=-0.00001, 
221                                       fr="Décalage radial du defaut sous revêtement (en m)", ),
222     # ANGLDEF
223     Azimut                   = SIMP ( statut="o", typ="R", defaut=0., 
224                                       fr="Position angulaire du défaut dans le cylindre de cuve (en degrés)", ),
225     # ALTIDEF
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)", ),
228     # POINDEF
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" ), ),
232     # CORRECPLASTIC
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" ), ),
236
237   ), # Fin BLOC SousRevetement
238
239   Debouchant = BLOC ( condition = "TypeInitial=='Debouchant'",
240
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" ), ),
246     # PROFDEF
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)", ),
249     # ANGLDEF
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)", ),
252     # ALTIDEF
253     Altitude                 = SIMP ( statut="o", typ="R", defaut=2., 
254                                       fr="Altitude du défaut sur le cylindre de la cuve (en m)", ),
255     # IRWIN
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" ), ),
259     # CORRECPLASTIC
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" ), ),
263
264   ), # Fin BLOC debouchant
265
266 ) # Fin PROC DEFAUT
267
268
269 #================================
270 # 3. Caracteristiques de la CUVE
271 #================================
272
273 CUVE = PROC (nom = "CUVE",
274              op = 68,
275              repetable = 'n',
276              fr = "Caractéristiques de la cuve", 
277
278 #===
279 # Liste des paramètres
280 #===
281
282   # TYPEGEOM
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
287                                ),
288
289 #====
290 # Definition des parametres selon le type de traitement de la geometrie
291 #====
292
293   Geometrie = BLOC ( condition = "TraitementGeometrie=='Topologie'",
294
295     # RINT
296     RayonInterne        = SIMP ( statut="o", typ="R", val_min=0.,  defaut=1.994, 
297                                  fr="Rayon interne de la cuve (en m)", ),
298
299     # DTV : comment preciser que REXT > RINT ?
300     # REXT
301     RayonExterne        = SIMP ( statut="o", typ="R", defaut=2.2015, 
302                                  fr="Rayon externe de la cuve (en m)", ),
303
304     # DTV : comment preciser que LREV < RINT ?
305     # LREV
306     EpaisseurRevetement = SIMP ( statut="o", typ="R", defaut=0.0075, 
307                                  fr="Epaisseur du revêtement (m)", ),
308
309     # LIGMIN
310     LigamentExterneMin  = SIMP ( statut="o", typ="R", defaut=0.75, 
311                                  fr="Ligament externe minimal avant rupture (% de l'épaisseur de cuve)", ),
312
313   ), # Fin BLOC  Geometrie
314
315   Maillage  = BLOC ( condition = "TraitementGeometrie=='Maillage'",
316
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
322
323 ) # Fin PROC CUVE
324
325 #====================================================
326 # 4. Modeles de fluence, d'irradiation et de tenacite
327 #====================================================
328
329 #=======================
330 # 4.1 Modeles de fluence
331 #=======================
332
333 MODELES = PROC ( nom = "MODELES",
334                  op = 68,
335                  repetable = 'n',
336                  fr = "Modèles de fluence, d'irradiation et de ténacité", 
337
338
339 #===
340 # Liste des paramètres
341 #===
342
343   Fluence = FACT ( statut="o",
344   
345     # MODELFLUENCE
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 
357                          ),
358
359 #====
360 # Definition des parametres selon le modele de fluence
361 #====
362
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 ?
365       # fmax
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
369
370     France        = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=12.7 (France)', ) ",
371       # fmax
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", ),
374       # KPFRANCE
375       KPFrance      = SIMP ( statut="o", typ="R",
376                              fr="Paramètre exponentiel du modèle France",
377                              defaut = 12.7, ),
378     ), # Fin BLOC France
379
380     ValeurImposee = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=0. (ValeurImposee)', ) ",
381       # fmax
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
385
386     SDM           = BLOC ( condition = " Modele in ( 'Donnees francaises du palier CPY (SDM)', ) ",
387       # fmax
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", ),
390     ), # Fin BLOC SDM
391
392     USNRC         = BLOC ( condition = " Modele in ( 'Regulatory Guide 1.99 rev 2 (USNRC)', ) ",
393       # fmax
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",
398                              defaut=9.4488, ),
399     ), # Fin BLOC USNRC
400
401     REV_2         = BLOC ( condition = " Modele in ( 'Dossier 900 MWe AP9701 rev 2 (REV_2)', ) ",
402       # fmax
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", ),
405     ), # Fin BLOC REV_2
406
407     SDM_Lissage   = BLOC ( condition = " Modele in ( 'Lissage du modele ajuste (SDM_Lissage)', ) ",
408       # fmax
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
412
413     GrandeDev     = BLOC ( condition = " Modele in ( 'Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)', ) ",
414       # fmax
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
418
419     GD_Cuve       = BLOC ( condition = " Modele in ( 'Grand developpement (GD_Cuve)', ) ",
420       # fmax
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
425
426   ), # Fin FACT Fluence
427
428 #==========================
429 # 4.2 Modeles d'irradiation
430 #==========================
431
432   Irradiation = FACT ( statut="o",
433
434     # TYPIRR
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
439                        ),
440
441 #====
442 # Definition des parametres selon le type d'irradiation
443 #====
444
445     RTndt = BLOC ( condition = "Option=='RTndt de la cuve a l instant de l analyse'",
446  
447       RTNDT = SIMP ( statut="o", typ="R", defaut=73., 
448                      fr="RTNDT de la cuve à l'instant de l'analyse (°C)", ),
449
450     ), # Fin BLOC RTndt
451
452     ModeleIrradiation = BLOC ( condition = "Option=='Modele d irradiation'",
453  
454       # MODELIRR
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
463                                 ),
464       # CU
465       TeneurCuivre     = SIMP ( statut="o", typ="R", defaut=0., 
466                                 fr="Teneur en cuivre (%)", ),
467       # NI
468       TeneurNickel     = SIMP ( statut="o", typ="R", defaut=0., 
469                                 fr="Teneur en nickel (%)", ),
470       # P
471       TeneurPhosphore  = SIMP ( statut="o", typ="R", defaut=0., 
472                                 fr="Teneur en phosphore (%)", ),
473
474       Parametres_USNRC = BLOC ( condition = " Modele in ( 'Joint Soude : Regulatory Guide 1.00 rev 2', 'Metal de Base : Regulatory Guide 1.00 rev 2' , ) ",
475         # RTimoy
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)", ),
478         # RTicov
479         CoefVariationRTndt   = SIMP ( statut="o", typ="R", defaut=0., 
480                                       fr="Coefficient de variation de la RTNDT initiale", ),
481         # USectDRT
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)", ),
484         # nbectDRTNDT
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 
488
489     ), # Fin BLOC ModeleIrradiation
490
491   ), # Fin FACT Irradiation
492
493 #========================
494 # 4.3 Modeles de tenacite
495 #========================
496
497   Tenacite = FACT ( statut = "o",
498
499     # MODELKIC
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
514                     ),
515
516
517 #====
518 # Definition des parametres selon le modele de tenacité
519 #====
520
521 # Modeles type RCC-M
522
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', ) ",
524
525       # nbectKIc
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 ", ),
528
529       # KICPAL
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)) ", ),
532
533       # KICCDVD
534       CoefficientVariation_KIc = SIMP ( statut="o", typ="R", defaut = 0.15, 
535                                         fr = "Coefficient de variation de la loi normale de K1c ", ),
536
537     ), # Fin BLOC KIc_RCCM
538
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', )",
540
541       # ARRETFISSURE
542       ArretDeFissure           = SIMP ( statut="o", typ="TXM", defaut="NON", 
543                                         fr="Prise en compte de l'arrêt de fissure",
544                                         into=( "OUI", "NON" ), ),
545
546       Parametres_Fissure       = BLOC ( condition = "ArretDeFissure=='OUI'",
547         # INCRDEF
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
551
552       KIa_RCCM = BLOC ( condition = "ArretDeFissure=='OUI'",
553
554         # nbectKIa
555         NbEcartType_MoyKIa       = SIMP ( statut="o", typ="R", defaut=-2., 
556                                           fr="Nombre d'écart-type par rapport à la moyenne de KIa (nb sigma) ", ),
557
558         # KIACDV
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- ", ),
561
562         # KIAPAL
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
566
567     ), # Fin BLOC Fissure
568
569 # Modeles type Weibull
570
571     KIc_MasterCuve = BLOC ( condition = " Modele in ( 'Weibull basee sur la master cuve', ) ",
572  
573       # T0WALLIN
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", ),
576
577       # fractKIc
578       Fractile_KIc       = SIMP ( statut="o", typ="R", defaut=5., 
579                                   fr="Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", ),
580
581     ), # Fin BLOC KIc_MasterCuve
582
583     KIc_Weibull1 = BLOC ( condition = " Modele in ( 'Weibull basee sur la master cuve (REME)', 'Weibull n°1 (etude ORNL)', ) ",
584
585       # fractKIc
586       Fractile_KIc       = SIMP ( statut="o", typ="R", defaut=5., 
587                                   fr="Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", ),
588
589     ), # Fin BLOC KIc_Weibull1
590
591     KIc_Weibull2 = BLOC ( condition = " Modele in ( 'Weibull n°2', 'Weibull n°3', ) ",
592
593       # nbectKIc
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 ", ),
596
597     ), # Fin BLOC KIc_Weibull2
598
599     Weibull_Generalisee = BLOC ( condition = " Modele in ( 'Weibull generalisee',) ",
600  
601       Coefficients       = Coef_WeibGen(),
602       # nbectKIc
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 ", ),
605
606     ), # Fin BLOC Weibull_Generalisee
607
608 # Modeles type exponentiel (Frama, LOGWOLF)
609
610     KIc_Exponentielle = BLOC ( condition = " Modele in ( 'Exponentielle n°1 (Frama)', 'Exponentielle n°2 (LOGWOLF)', ) ",
611
612       # nbectKIc
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 ", ),
615
616     ), # Fin BLOC KIc_Exponentielle
617
618   ), # Fin FACT Tenacite
619
620 ) # Fin PROC MODELES
621
622
623 #==================
624 # 5. Initialisation
625 #==================
626
627 INITIALISATION = PROC ( nom = "INITIALISATION",
628                         op = 68,
629                         repetable = 'n',
630                         fr = "Initialisation : instant initial, profils radiaux de température et contraintes", 
631
632   TemperatureInitiale = FACT ( statut = "o",
633
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' ), ),
642
643   ), # Fin FACT TemperatureInitiale
644
645   ContraintesInitiales = FACT ( statut = "o",
646
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' ), ),
655
656   ), # Fin FACT ContraintesInitiales
657
658   # INSTINIT
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) ", ),
661
662 ) # Fin PROC INITIALISATION
663
664
665 #==================================
666 # 6. CARACTERISTIQUES DU REVETEMENT
667 #==================================
668
669 REVETEMENT = PROC ( nom = "REVETEMENT",
670                     op = 68,
671                     repetable = 'n',
672                     fr = "Caracteristiques du revêtement", 
673
674   # KTHREV
675   ConditionLimiteThermiqueREV = SIMP ( statut="o", typ="TXM",
676                                        fr="Option pour définir les caractéristiques du revêtement ",
677                                        into=( "ENTHALPIE", "CHALEUR",),
678                                        ),
679
680   EnthalpieREV = BLOC ( condition = "ConditionLimiteThermiqueREV=='ENTHALPIE'",
681
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' ), ),
690
691   ), # Fin BLOC EnthalpieREV
692
693
694   ChaleurREV = BLOC ( condition = "ConditionLimiteThermiqueREV=='CHALEUR'",
695
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' ), ),
704
705   ), # Fin BLOC ChaleurREV
706
707   ConductiviteREV = FACT (statut = "o",
708
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' ), ),
717
718   ), # Fin FACT ConductiviteREV
719
720   ModuleYoungREV = FACT (statut = "o",
721
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' ), ),
730
731   ), # Fin FACT ModuleYoungREV
732
733   CoeffDilatThermREV = FACT (statut = "o",
734
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' ), ),
743
744   ), # Fin FACT CoeffDilatThermREV
745
746   LimiteElasticiteREV = FACT (statut = "o",
747
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' ), ),
756
757   ), # Fin FACT LimiteElasticiteREV
758
759   AutresParametresREV = FACT (statut = "o",
760
761     # TREFREV
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) ", ),
764     # TDETREV
765     TemperaturePourCoefDilatThermREV = SIMP ( statut="o", typ="R", defaut=287.,
766                                               fr="Température de définition du coefficient de dilatation thermique (°C) ", ),
767     # NUREV
768     CoefficientPoissonREV            = SIMP ( statut="o", typ="R", defaut=0.3,
769                                               fr="Coefficient de Poisson ", ),
770
771   ), # Fin FACT AutresParametresREV
772
773 ) # Fin PROC REVETEMENT
774
775
776 #=====================================
777 # 7. CARACTERISTIQUES DU METAL DE BASE
778 #=====================================
779
780 METAL_BASE = PROC ( nom = "METAL_BASE",
781                     op = 68,
782                     repetable = 'n',
783                     fr = "Caracteristiques du metal de base", 
784
785   # KTHMDB
786   ConditionLimiteThermiqueMDB = SIMP ( statut="o", typ="TXM",
787                                        fr="Option pour definir les caractéristiques du revêtement ",
788                                        into=( "ENTHALPIE", "CHALEUR",), ),
789
790   EnthalpieMDB = BLOC ( condition = "ConditionLimiteThermiqueMDB=='ENTHALPIE'",
791
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' ), ),
800
801   ), # Fin BLOC EnthalpieMDB
802
803   ChaleurMDB = BLOC ( condition = "ConditionLimiteThermiqueMDB=='CHALEUR'",
804
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' ), ),
813
814   ), # Fin BLOC ChaleurMDB
815
816   ConductiviteMDB = FACT ( statut = "o",
817
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' ), ),
826
827   ), # Fin FACT ConductiviteMDB
828
829   ModuleYoungMDB = FACT ( statut="o",
830
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' ), ),
839
840   ), # Fin FACT ModuleYoungMDB
841
842   CoeffDilatThermMDB = FACT ( statut="o",
843
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' ), ),
852
853   ), # Fin FACT CoeffDilatThermMDB
854
855   AutresParametresMDB = FACT ( statut = "o",
856
857     # TREFMDB
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) ", ),
860     # TDETMDB
861     TemperaturePourCoefDilatThermMDB = SIMP ( statut="o", typ="R", defaut=287.,
862                                               fr="Température de définition du coefficient de dilatation thermique (°C) ", ),
863     # NUMDB
864     CoefficientPoissonMDB            = SIMP ( statut="o", typ="R", defaut=0.3,
865                                               fr="Coefficient de Poisson ", ),
866
867   ), # Fin FACT TemperatureDeformationNulleMDB
868
869 ) # Fin PROC METAL_BASE
870
871
872 #===============================
873 # 8. TRANSITOIRE THERMOMECANIQUE
874 #===============================
875
876 TRANSITOIRE = PROC ( nom = "TRANSITOIRE",
877                      op = 68,
878                      repetable = 'n',
879                      fr = "Description du transitoire thermohydraulique", 
880
881   Pression = FACT ( statut = "o",
882
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' ), ),
891
892   ), # FIN FACT Pression
893
894   # TYPCLTH
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
902                                         ),
903
904   TemperatureImposeeParoi      = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee en paroi', ) ",
905
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' ), ),
914
915   ), # Fin BLOC TemperatureImposeeParoi
916
917   FluxChaleur                  = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Flux de chaleur impose en paroi', ) ",
918
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' ), ),
927
928   ), # Fin BLOC FluxChaleur
929
930   TemperatureImposeeFluide     = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', 'Temperature imposee du fluide et debit d injection de securite', ) ",
931
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' ), ),
940
941   ), # Fin BLOC TemperatureImposeeFluide
942
943   CoefficientEchange          = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', ) ",
944
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' ), ),
953
954   ), # Fin BLOC CoefficientEchange
955
956   DebitMassique               = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ",
957
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' ), ),
966
967   ), # Fin BLOC DebitMassique
968
969   TemperatureInjection        = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ",
970
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
980
981   Creare                      = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ",
982
983     # DH
984     DiametreHydraulique              = SIMP ( statut="o", typ="R", defaut=-2.,
985                                               fr="Diamètre hydraulique (m) ", ),
986     # SECTION
987     SectionEspaceAnnulaire           = SIMP ( statut="o", typ="R", defaut=-2.,
988                                               fr="Section espace annulaire (m2) ", ),
989     # DELTA
990     HauteurCaracConvectionNaturelle  = SIMP ( statut="o", typ="R", defaut=-2.,
991                                               fr="Hauteur caractéristique convection naturelle (m) ", ),
992     # ALPHA_CF
993     CoeffVestale_ConvectionForcee    = SIMP ( statut="o", typ="R", defaut=1.,
994                                               fr="Coefficient Vestale convection forcée (-) ", ),
995     # ALPHA_CN
996     CoeffVestale_ConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=1.,
997                                               fr="Coefficient Vestale convection naturelle (-) ", ),
998     # EPS
999     CritereConvergenceRelative       = SIMP ( statut="o", typ="R", defaut=0.00001,
1000                                               fr = "Critère convergence relative (-) ", ),
1001     # VM
1002     VolumeMelange_CREARE             = SIMP ( statut="o", typ="R", defaut=-2.,
1003                                               fr = "Volume de mélange CREARE (m3) ", ),
1004     # T0
1005     TemperatureInitiale_CREARE       = SIMP ( statut="o", typ="R", defaut=-2.,
1006                                               fr="Température initiale CREARE (°C) ", ),
1007     # SE
1008     SurfaceEchange_FluideStructure   = SIMP ( statut="o", typ="R", defaut=-2.,
1009                                               fr="Surface d'échange fluide/structure (m2) ", ),
1010
1011   ), # Fin BLOC Creare
1012
1013   DebitInjection              = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ",
1014
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' ), ),
1023
1024   ), # Fin BLOC DebitInjection
1025
1026
1027   Vestale                     = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ",
1028
1029     # DH
1030     DiametreHydraulique              = SIMP ( statut="o", typ="R", defaut=-2.,
1031                                               fr="Diamètre hydraulique (m) ", ),
1032     # SECTION
1033     SectionEspaceAnnulaire           = SIMP ( statut="o", typ="R", defaut=-2.,
1034                                               fr="Section espace annulaire (m2) ", ),
1035     # DELTA
1036     HauteurCaracConvectionNaturelle  = SIMP ( statut="o", typ="R", defaut=-2.,
1037                                               fr="Hauteur caractéristique convection naturelle (m) ", ),
1038     # ALPHA_CF
1039     CoeffVestale_ConvectionForcee    = SIMP ( statut="o", typ="R", defaut=1.,
1040                                               fr="Coefficient d'échange Vestale convection forcée (-) ", ),
1041     # ALPHA_CN
1042     CoeffVestale_ConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=1.,
1043                                               fr="Coefficient d'échange Vestale convection naturelle (-) ", ),
1044     # EPS
1045     CritereConvergenceRelative       = SIMP ( statut="o", typ="R", defaut=0.00001,
1046                                               fr="Critère convergence relative (-) ", ),
1047
1048   ), # Fin BLOC Vestale
1049
1050 ) # Fin PROC TRANSITOIRE