]> SALOME platform Git repositories - tools/eficas.git/commitdiff
Salome HOME
*** empty log message ***
authorPascale Noyret <pascale.noyret@edf.fr>
Tue, 27 Jan 2009 09:29:57 +0000 (09:29 +0000)
committerPascale Noyret <pascale.noyret@edf.fr>
Tue, 27 Jan 2009 09:29:57 +0000 (09:29 +0000)
Cuve2dg/Cuve2dg_Cata_V1.py
Cuve2dg/Cuve2dg_Cata_V2.py [new file with mode: 0644]
Cuve2dg/catalogues_cuve2dg.ini

index 61d9b721a4e76e828d7a121edaf64e060b7c38ea..47e23641bbcaeb92eb80e5c1c441f8540549826f 100644 (file)
@@ -12,9 +12,12 @@ class variable ( ASSD ) : pass
 
 
 #CONTEXT.debug = 1
-JdC = JDC_CATA ( code = 'OPENTURNS_STUDY',
+JdC = JDC_CATA ( code = 'CUVE1D-DEFAILLGLOB',
                  execmodul = None,
-                 regles = ( AU_MOINS_UN ( 'CRITERIA' ), ),
+                 regles = ( AU_MOINS_UN ('OPTIONS','DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'), 
+                            AU_MOINS_UN ( 'FIN' ), 
+                            A_CLASSER ( ('OPTIONS', 'DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'),'FIN')
+                          )
                  ) # Fin JDC_CATA
 
 # --------------------------------------------------
@@ -27,628 +30,1461 @@ JdC = JDC_CATA ( code = 'OPENTURNS_STUDY',
 
 
 #================================
-# 1. Definition des LOIS
+# 1. Definition des OPTIONS
 #================================
 
 # Nota : les variables de type OPER doivent etre en majuscules !
 # Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
-DISTRIBUTION = OPER ( nom = "DISTRIBUTION",
-                      sd_prod = loi,
-                      op = 68,
-                      fr = "Definitions des lois marginales utilisees par les variables d'entree", 
-                      
-                      
-#====
-# Type de distribution
-#====
 
-  Kind = SIMP ( statut = "o", typ = "TXM",
-                into = ( "Beta",
-                         "Exponential",
-                         "Gamma",
-                         "Geometric",
-                         "Gumbel",
-                         "Histogram",
-                         "Logistic",
-                         "LogNormal",
-                         "MultiNomial",
-                         "Normal",
-                         "TruncatedNormal",
-                         "Poisson",
-                         "Student",
-                         "Triangular",
-                         "Uniform",
-                         "UserDefined",
-                         "Weibull",
-                         ),
-                fr = "Choix du type de la loi marginale",
-                ang = "1D marginal distribution",
+OPTIONS = OPER ( nom = "OPTIONS",
+                 sd_prod = loi,
+                 op = 68,
+                 fr = "Definitions des options", 
+
+#===
+# Liste des paramètres
+#===
+
+  INCRTPS = SIMP ( statut = "o",
+                   typ = "I",
+                   defaut = "1",
+                  max = 1,
+                   val_max = 100,
+                   fr = "Increment temporel (=1 pour calcul deterministe)",
+                  ),
+
+  DTPREC = SIMP ( statut = "o",
+                  typ = "R",
+                  defaut = "0.1",
+                  max = 1,
+                  val_max = 1.,
+                  fr = "Increment maximum d'evolution de la temperature par noeud et par instant (°C)",
                 ),
 
+  DTARCH = SIMP ( statut = "o",
+                  typ = "R",
+                  defaut = "1000.",
+                 max = 1,
+                  val_max = 1000.,
+                  fr = "Increment maximum de temps pour l'affichage (s)",
+                ),
+
+  NBO = SIMP ( statut = "o",
+               typ = "R",
+              max=1,
+               val_max = 1000.,
+               fr = "Nombre de noeuds a considerer dans le maillage interne",
+              ),
+
+  Liste_instants = SIMP ( statut = "o",
+                          typ = "R",
+                          max = "**",
+                          fr = "Liste des instants ",
+                        ),
+
+) # Fin OPER OPTIONS
+
+#================================
+# 2. Caracteristiques du DEFAUT
+#================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+DEFAUT = OPER ( nom = "DEFAUT",
+                sd_prod = loi,
+                op = 68,
+                fr = "Caracteristiques du defaut", 
+
+#===
+# Liste des paramètres
+#===
+
+  TYPEDEF = SIMP ( statut = "o", typ = "TXM",
+                   into = ( "DSR",
+                            "DD",
+                          ),
+                   #defaut = "DSR",
+                   fr = "Type de defaut : sous revetement ou debouchant",
+                  ),
+
 #====
-# Definition des parametres selon le type de la loi
+# Definition des parametres selon le type du defaut
 #====
 
-  BETA = BLOC ( condition = " Kind in ( 'Beta', ) ",
+  Parametres_DSR = BLOC ( condition = " TYPEDEF in ( 'DSR', ) ",
 
-                  Settings = SIMP ( statut = "o",
-                                       typ = "TXM",
-                                       max = 1,
-                                       into = ( "RT", "MuSigma" ),
-                                       defaut = "RT",
-                                       fr = "Parametrage de la loi beta",
-                                       ang = "Beta distribution parameter set",
-                                       ),
+                  ORIEDEF = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( "LONGITUD", "CIRCONF" ),
+                                   #defaut = "LONGITUD",
+                                   fr = "Orientation du defaut : longitudinale ou circonferentielle",
+                                 ),
 
-                  RT_Parameters = BLOC ( condition = " Settings in ( 'RT', ) ",
+                  PROFDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   #defaut = "0.006",
+                                   max = 1,
+                                   val_max = 1.,
+                                   fr = "Profondeur radiale du defaut (m)",
+                                 ),
 
-                                      R = SIMP ( statut = "o",
-                                                 typ = "R",
-                                                 max = 1,
-                                                 val_min = 0.,
-                                                 fr = "Parametre R de la loi",
-                                                 ang = "R parameter",
-                                                 ),
+                  OPTLONG = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( "VALEUR", "RAPPORT" ),
+                                   #defaut = "VALEUR",
+                                   fr = "Option pour caracteriser la longueur du defaut : soit par valeur, soit par un rapport LONG/PROF",
+                                 ),
+
+                    Option_Valeur = BLOC ( condition = "OPTLONG in ( 'VALEUR', ) ",
+
+                            LONGDEF = SIMP ( statut = "o",
+                                             typ = "R",
+                                             #defaut = "0.060",
+                                             max = 1,
+                                             val_max = 1.,
+                                             fr = "Longueur du defaut sous revetement (m)",
+                                            ),
+
+                                         ), # Fin BLOC Option_Valeur
+
+                    Option_Rapport = BLOC ( condition = "OPTLONG in ( 'RAPPORT', ) ",
 
-                                      # T > R
-                                      T = SIMP ( statut = "o",
+                            LONGSURPROF = SIMP ( statut = "o",
                                                  typ = "R",
+                                                 #defaut = "6.",
                                                  max = 1,
-                                                 val_min = 0.,
-                                                 fr = "Parametre T de la loi | T > R",
-                                                 ang = "T parameter | T > R",
-                                                 ),
+                                                 val_max = 100.,
+                                                 fr = "Rapport longueur/profondeur du defaut sous revetement",
+                                                ),
 
-                                      ), # Fin BLOC RT_Parameters
+                                          ), # Fin BLOC Option_Rapport
 
+                  DECADEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   #defaut = "-0.00001",
+                                   fr = "Decalage radial du defaut sous revetement (m)",
+                                  ),
 
-                  MuSigma_Parameters = BLOC ( condition = " Settings in ( 'MuSigma', ) ",
+                  ANGLDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.",
+                                   fr = "Coordonnee angulaire du defaut (degres)",
+                                  ),
 
-                                      Mu = SIMP ( statut = "o",
-                                                  typ = "R",
-                                                  max = 1,
-                                                  fr = "Parametre Mu de la loi",
-                                                  ang = "Mu parameter",
-                                                  ),
+                  ALTIDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "2.",
+                                   fr = "Altitude du defaut (m)",
+                                  ),
 
-                                      Sigma = SIMP ( statut = "o",
-                                                     typ = "R",
-                                                     max = 1,
-                                                     val_min = 0.,
-                                                     fr = "Parametre Sigma de la loi | Sigma > 0",
-                                                     ang = "Sigma parameter | Sigma > 0",
-                                                     ),
+                  POINDEF = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( "A", "B" ),
+                                   defaut = "A",
+                                   fr = "Choix du point considere du defaut sous revetement",
+                                  ),
 
-                                      ), # Fin BLOC MuSigma_Parameters
+                  ARRETFISSURE = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( "OUI", "NON" ),
+                                        defaut = "NON",
+                                        fr = "Prise en compte de l arret de fissure",
+                                       ),
 
+                  INCRDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.005",
+                                   fr = "Increment de la taille de fissure (m)",
+                                  ),
 
-                  A = SIMP ( statut = "o",
-                             typ = "R",
-                             max = 1,
-                             fr = "Parametre A de la loi",
-                             ang = "A parameter",
-                             ),
+                  CORRECPLASTIC = SIMP ( statut = "o",
+                                         typ = "TXM",
+                                         into = ( "OUI", "NON" ),
+                                         defaut = "NON",
+                                         fr = "Prise en compte de la correction plastique BETA ",
+                                        ),
 
-                  # B > A
-                  B = SIMP ( statut = "o",
-                             typ = "R",
-                             max = 1,
-                             fr = "Parametre B de la loi | B > A",
-                             ang = "B parameter | B > A",
-                             ),
+                         ), # Fin BLOC Parametres_DSR
 
-  ), # Fin BLOC BETA
+  Parametres_DD = BLOC ( condition = " TYPEDEF in ( 'DD', ) ",
 
+                  ORIEDEF = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( "LONGITUD", "CIRCONF" ),
+                                   #defaut = "LONGITUD",
+                                   fr = "Orientation du defaut : longitudinale ou circonferentielle",
+                                  ),
 
+                  PROFDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   #defaut = "0.006",
+                                   max = 1,
+                                   val_max = 1.,
+                                   fr = "Profondeur radiale du defaut (m)",
+                                  ),
 
-  EXPONENTIAL = BLOC ( condition = " Kind in ( 'Exponential', ) ",
+                  ANGLDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.",
+                                   fr = "Coordonnee angulaire du defaut (degres)",
+                                  ),
 
-                         Lambda = SIMP ( statut = "o",
-                                         typ = "R",
-                                         max = 1,
-                                         val_min = 0.,
-                                         fr = "Parametre Lambda | Lambda > 0",
-                                         ang = "Lambda parameter | Lambda > 0",
-                                         ),
+                  ALTIDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "2.",
+                                   fr = "Altitude du defaut (m)",
+                                  ),
 
-                         Gamma = SIMP ( statut = "o",
-                                        typ = "R",
-                                        max = 1,
-                                        fr = "Parametre Gamma",
-                                        ang = "Gamma parameter",
-                                        ),
+                  ARRETFISSURE = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( "OUI", "NON" ),
+                                        defaut = "NON",
+                                        fr = "Prise en compte de l arret de fissure",
+                                       ),
 
-  ), # Fin BLOC EXPONENTIAL
+                  INCRDEF = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.005",
+                                   fr = "Increment de la taille de fissure (m)",
+                                  ),
 
+                  IRWIN = SIMP ( statut = "o",
+                                 typ = "TXM",
+                                 into = ( "OUI", "NON" ),
+                                 defaut = "NON",
+                                 fr = "Prise en compte de la correction plastique d'Irwin ",
+                                ),
 
+                  CORRECPLASTIC = SIMP ( statut = "o",
+                                         typ = "TXM",
+                                         into = ( "OUI", "NON" ),
+                                         defaut = "NON",
+                                         fr = "Prise en compte de la correction plastique BETA ",
+                                        ),
 
-  GAMMA = BLOC ( condition = " Kind in ( 'Gamma', ) ",
+                        ), # Fin BLOC Parametres_DD 
 
-                   Settings = SIMP ( statut = "o",
-                                        typ = "TXM",
-                                        max = 1,
-                                        into = ( "KLambda", "MuSigma" ),
-                                        defaut = "KLambda",
-                                        fr = "Parametrage de la loi gamma",
-                                        ang = "Gamma distribution parameter set",
-                                        ),
+) # Fin OPER DEFAUT
 
-                   KLambda_Parameters = BLOC ( condition = " Settings in ( 'KLambda', ) ",
 
-                                       K = SIMP ( statut = "o",
-                                                  typ = "R",
-                                                  max = 1,
-                                                  val_min = 0.,
-                                                  fr = "Parametre K de la loi | K > 0",
-                                                  ang = "K parameter | K > 0",
-                                                  ),
+#================================
+# 3. Caracteristiques de la CUVE
+#================================
 
-                                       Lambda = SIMP ( statut = "o",
-                                                       typ = "R",
-                                                       max = 1,
-                                                       val_min = 0.,
-                                                       fr = "Parametre Lambda de la loi | Lambda > 0",
-                                                       ang = "Lambda parameter | Lambda > 0",
-                                                       ),
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+CUVE = OPER (nom = "CUVE",
+             sd_prod = loi,
+             op = 68,
+             fr = "Caracteristiques de la cuve", 
 
-                                       ), # Fin BLOC KLambda_Parameters
+#===
+# Liste des paramètres
+#===
 
+  TYPEGEOM = SIMP ( statut = "o",
+                    typ = "TXM",
+                   into = ( "GEOMETRIE", "MAILLAGE"), 
+                    #defaut = "GEOMETRIE",
+                    fr = "Traitement de la geometrie d'une cuve",
+                   ),
 
-                   MuSigma_Parameters = BLOC ( condition = " Settings in ( 'MuSigma', ) ",
 
-                                       Mu = SIMP ( statut = "o",
-                                                   typ = "R",
-                                                   max = 1,
-                                                   defaut = 0.0,
-                                                   fr = "Parametre Mu de la loi",
-                                                   ang = "Mu parameter",
-                                                   ),
+#====
+# Definition des parametres selon le type de traitement de la geometrie
+#====
 
-                                       Sigma = SIMP ( statut = "o",
-                                                      typ = "R",
-                                                      max = 1,
-                                                      defaut = 1.0,
-                                                      val_min = 0.,
-                                                      fr = "Parametre Sigma de la loi | Sigma > 0",
-                                                      ang = "Sigma parameter | Sigma > 0",
-                                                      ),
+  Geometrie = BLOC ( condition = " TYPEGEOM in ( 'GEOMETRIE', ) ",
+
+                  RINT = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "1.994",
+                                fr = "Rayon interne de la cuve (m)",
+                               ),
+
+                  REXT = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "2,2015",
+                                fr = "Rayon externe de la cuve (m)",
+                               ),
 
-                                       ), # Fin BLOC MuSigma_Parameters
+                  LREV = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "0.0075",
+                                fr = "Epaisseur du revetement (m)",
+                               ),
 
-                   Gamma = SIMP ( statut = "o",
+                  LIGMIN = SIMP ( statut = "o",
                                   typ = "R",
-                                  max = 1,
-                                  fr = "Parametre Gamma",
-                                  ang = "Gamma parameter",
-                                  ),
+                                  defaut = "0.75",
+                                  fr = "Ligament externe minimal avant rupture (% de l'epaisseur de cuve)",
+                                 ),
 
+                     ), # Fin BLOC  Geometrie
 
-  ), # Fin BLOC GAMMA
+  Maillage = BLOC ( condition = " TYPEGEOM in ( 'MAILLAGE', ) ",
 
+                  Liste_abscisses = SIMP ( statut = "o",
+                                           typ = "R",
+                                           max = "**",
+                                           fr = "Liste des abscisses (m) A FAIRE",
+                                          ),
+                   ), # Fin BLOC Maillage
 
+) # Fin OPER CUVE
 
-  GEOMETRIC = BLOC ( condition = " Kind in ( 'Geometric', ) ",
+#====================================================
+# 4. Modeles de fluence, d'irradiation et de tenacite
+#====================================================
 
-                       P = SIMP ( statut = "o",
-                                  typ = "R",
-                                  max = 1,
-                                  val_min = 0.,
-                                  val_max = 1.,
-                                  fr = "Parametre P | 0 < P < 1",
-                                  ang = "P parameter | 0 < P < 1",
-                                  ),
+#=======================
+# 4.1 Modeles de fluence
+#=======================
 
-  ), # Fin BLOC GEOMETRIC
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+MODELES = OPER ( nom = "MODELES",
+                 sd_prod = loi,
+                 op = 68,
+                 fr = "Modeles de fluence, d'irradiation et de tenacite", 
 
 
+#===
+# Liste des paramètres
+#===
 
-  GUMBEL = BLOC ( condition = " Kind in ( 'Gumbel', ) ",
+  MODELFLUENCE = SIMP ( statut = "o",
+                        typ = "TXM",
+                       into = ( "Reglementaire", "France", "ValeurImposee", "SDM", "USNRC", "REV_2", "SDM_Lissage", "GrandeDev", "GD_Cuve"), 
+                        #defaut = "Reglementaire",
+                        fr = "Modele d'attenuation de la fluence dans l'epaisseur de la cuve",
+                       ),
 
-                    Settings = SIMP ( statut = "o",
-                                         typ = "TXM",
-                                         max = 1,
-                                         into = ( "AlphaBeta", "MuSigma" ),
-                                         defaut = "AlphaBeta",
-                                         fr = "Parametrage de la loi gumbel",
-                                         ang = "Gumbel distribution parameter set",
-                                         ),
 
-                    AlphaBeta_Parameters = BLOC ( condition = " Settings in ( 'AlphaBeta', ) ",
+#====
+# Definition des parametres selon le modele de fluence
+#====
 
-                                        Alpha = SIMP ( statut = "o",
-                                                       typ = "R",
-                                                       max = 1,
-                                                       val_min = 0.,
-                                                       fr = "Parametre Alpha de la loi | Alpha > 0",
-                                                       ang = "Alpha parameter | Alpha > 0",
-                                                       ),
+  Reglementaire = BLOC ( condition = " MODELFLUENCE in ( 'Reglementaire', ) ",
 
-                                        Beta = SIMP ( statut = "o",
-                                                      typ = "R",
-                                                      max = 1,
-                                                      fr = "Parametre Beta de la loi",
-                                                      ang = "Beta parameter",
-                                                      ),
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
 
-                                        ), # Fin BLOC AlphaBeta_Parameters
+                        ), # Fin BLOC Reglementaire
 
+  France = BLOC ( condition = " MODELFLUENCE in ( 'France', ) ",
 
-                    MuSigma_Parameters = BLOC ( condition = " Settings in ( 'MuSigma', ) ",
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
 
-                                        Mu = SIMP ( statut = "o",
-                                                    typ = "R",
-                                                    max = 1,
-                                                    fr = "Parametre Mu de la loi",
-                                                    ang = "Mu parameter",
-                                                    ),
+                 KPFRANCE = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "12.7",
+                                    fr = "Parametre exponentiel du modele France",
+                                   ),
 
-                                        Sigma = SIMP ( statut = "o",
-                                                       typ = "R",
-                                                       max = 1,
-                                                       val_min = 0.,
-                                                       fr = "Parametre Sigma de la loi | Sigma > 0",
-                                                       ang = "Sigma parameter | Sigma > 0",
-                                                       ),
+                   ), # Fin BLOC France
 
-                                        ), # Fin BLOC MuSigma_Parameters
+  ValeurImposee = BLOC ( condition = " MODELFLUENCE in ( 'ValeurImposee', ) ",
 
-  ), # Fin BLOC GUMBEL
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
 
+                         ), # Fin BLOC ValeurImposee
 
+  SDM = BLOC ( condition = " MODELFLUENCE in ( 'SDM', ) ",
 
-  HISTOGRAM = BLOC ( condition = " Kind in ( 'Histogram', ) ",
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
 
-                       Sup = SIMP ( statut = "o",
-                                    typ = "R",
-                                    max = 1,
-                                    fr = "Borne superieure de la distribution",
-                                    ang = "Upper bound",
-                                    ),
-
-                       # Il faut definir une collection de couples ( x,p ) 
-                       Values = SIMP ( statut = 'o',
-                                       typ = 'R',
-                                       max = '**',
+              ), # Fin BLOC SDM
+
+  USNRC = BLOC ( condition = " MODELFLUENCE in ( 'USNRC', ) ",
+
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+                 KPUS = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "9.4488",
+                                fr = "Parametre exponentiel du modele US",
+                               ),
+
+                ), # Fin BLOC USNRC
+
+  REV_2 = BLOC ( condition = " MODELFLUENCE in ( 'REV_2', ) ",
+
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+                ), # Fin BLOC REV_2
+
+  SDM_Lissage = BLOC ( condition = " MODELFLUENCE in ( 'SDM_Lissage', ) ",
+
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+                      ), # Fin BLOC SDM_Lissage
+
+  GrandeDev = BLOC ( condition = " MODELFLUENCE in ( 'GrandeDev', ) ",
+
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+                     ), # Fin BLOC GrandeDev
+
+  GD_Cuve = BLOC ( condition = " MODELFLUENCE in ( 'GD_Cuve', ) ",
+
+                  fmax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+                  COEFFLUENCE1 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "5.8",
+                                        fr = "Fluence a l'azimut 0 (10^19 n/cm)",
                                        ),
 
-  ), # Fin BLOC HISTOGRAM
+                  COEFFLUENCE2 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "5.48",
+                                        fr = "Fluence a l'azimut 5 (10^19 n/cm)",
+                                       ),
 
+                  COEFFLUENCE3 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "4.46",
+                                        fr = "Fluence a l'azimut 10 (10^19 n/cm)",
+                                       ),
 
+                  COEFFLUENCE4 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "3.41",
+                                        fr = "Fluence a l'azimut 15 (10^19 n/cm)",
+                                       ),
 
-  LOGNORMAL = BLOC ( condition = " Kind in ( 'LogNormal', ) ",
+                  COEFFLUENCE5 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "3.37",
+                                        fr = "Fluence a l'azimut 20 (10^19 n/cm)",
+                                       ),
 
-                     Settings = SIMP ( statut = "o",
-                                       typ = "TXM",
-                                       max = 1,
-                                       into = ( "MuSigmaLog", "MuSigma", "MuSigmaOverMu" ),
-                                       defaut = "MuSigmaLog",
-                                       fr = "Parametrage de la loi lognormale",
-                                       ang = "Lognormal distribution parameter set",
+                  COEFFLUENCE6 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "3.16",
+                                        fr = "Fluence a l'azimut 25 (10^19 n/cm)",
                                        ),
 
-                     MuSigma_Parameters = BLOC ( condition = " Settings in ( 'MuSigma', ) ",
+                  COEFFLUENCE7 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "2.74",
+                                        fr = "Fluence a l'azimut 30 (10^19 n/cm)",
+                                       ),
 
-                                                 Mu = SIMP ( statut = "o",
-                                                             typ = "R",
-                                                             max = 1,
-                                                             fr = "Parametre Mu de la loi | Mu > Gamma",
-                                                             ang = "Mu parameter | Mu > Gamma",
-                                                             ),
+                  COEFFLUENCE8 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "2.25",
+                                        fr = "Fluence a l'azimut 35 (10^19 n/cm)",
+                                       ),
 
-                                                 Sigma = SIMP ( statut = "o",
-                                                                typ = "R",
-                                                                max = 1,
-                                                                val_min = 0.,
-                                                                fr = "Parametre Sigma de la loi | Sigma > 0",
-                                                                ang = "Sigma parameter | Sigma > 0",
-                                                                ),
+                  COEFFLUENCE9 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "1.89",
+                                        fr = "Fluence a l'azimut 40 (10^19 n/cm)",
+                                       ),
 
-                                                 ), # Fin BLOC MuSigma_Parameters
+                  COEFFLUENCE10 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "1.78",
+                                        fr = "Fluence a l'azimut 45 (10^19 n/cm)",
+                                       ),
 
-                     MuSigmaOverMu_Parameters = BLOC ( condition = " Settings in ( 'MuSigmaOverMu', ) ",
+                  ), # Fin BLOC GD_Cuve
 
-                                                 Mu = SIMP ( statut = "o",
-                                                             typ = "R",
-                                                             max = 1,
-                                                             fr = "Parametre Mu de la loi | Mu > Gamma",
-                                                             ang = "Mu parameter | Mu > Gamma",
-                                                             ),
+#==========================
+# 4.2 Modeles d'irradiation
+#==========================
 
-                                                 SigmaOverMu = SIMP ( statut = "o",
-                                                                typ = "R",
-                                                                max = 1,
-                                                                val_min = 0.,
-                                                                fr = "Parametre SigmaOverMu de la loi | SigmaOverMu > 0",
-                                                                ang = "SigmaOverMu parameter | SigmaOverMu > 0",
-                                                                ),
+  TYPIRR = SIMP ( statut = "o",
+                  typ = "TXM",
+                 into = ( "RTNDT", "FLUENCE" ),
+                  #defaut = "RTNDT",
+                  fr = "Type d'irradiation ",
+                 ),
 
-                                                 ), # Fin BLOC MuSigmaOverMu_Parameters
+#====
+# Definition des parametres selon le type d'irradiation
+#====
 
-                     MuSigmaLog_Parameters = BLOC ( condition = " Settings in ( 'MuSigmaLog', ) ",
+  Parametres_RTNDT = BLOC ( condition = " TYPIRR in ( 'RTNDT', ) ",
+                  RTNDT = SIMP ( statut = "o",
+                                 typ = "R",
+                                 defaut = "73.",
+                                 fr = "RTNDT finale (°C)",
+                                ),
 
-                                                    MuLog = SIMP ( statut = "o",
-                                                                   typ = "R",
-                                                                   max = 1,
-                                                                   fr = "Parametre Mu log de la loi",
-                                                                   ang = "Mu log parameter",
-                                                                   ),
+                            ), # Fin BLOC Parametres_RTNDT
 
-                                                    SigmaLog = SIMP ( statut = "o",
-                                                                      typ = "R",
-                                                                      max = 1,
-                                                                      val_min = 0.,
-                                                                      fr = "Parametre Sigma log de la loi | SigmaLog > 0",
-                                                                      ang = "Sigma log parameter | SigmaLog > 0",
-                                                                      ),
-                                            
-                                                    ), # Fin BLOC MuSigmaLog_Parameters
-
-                     Gamma = SIMP ( statut = "o",
-                                    typ = "R",
-                                    max = 1,
-                                    fr = "Parametre Gamma",
-                                    ang = "Gamma parameter",
-                                    ),
+  Parametres_FLUENCE = BLOC ( condition = " TYPIRR in ( 'FLUENCE', ) ",
+                  MODELIRR = SIMP ( statut = "o",
+                                    typ = "TXM",
+                                   into = ( "HOUSSIN", "PERSOZ", "LEFEBVRE", "USNRCmdb", "BRILLAUD", "USNRCsoud" ),
+                                    #defaut = "HOUSSIN",
+                                    fr = "Modele d'irradiation pour virole ou joint soude",
+                                   ),
+
+                  CU = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "Teneur en cuivre (%)",
+                             ),
 
-   ), # Fin BLOC LOGNORMAL
+                  Ni = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "Teneur en nickel (%)",
+                             ),
 
+                  P = SIMP ( statut = "o",
+                             typ = "R",
+                             defaut = "0.",
+                             fr = "Teneur en phosphore (%)",
+                            ),
 
+                  RTimoy = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.",
+                                  fr = "Moyenne de RTNDT : virole C1 de cuve Chinon : mdb=>-17.°C et js=>42.°C (HT-56/05/038 : p.52)",
+                                 ),
+
+                  RTicov = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.",
+                                  fr = "Coefficient de variation de la RTNDT initiale",
+                                 ),
+
+  Parametres_USNRC = BLOC ( condition = " MODELIRR in ( 'USNRCsoud', 'USNRCmdb' , ) ",
+                  USectDRT = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "28.",
+                                    fr = "pour modeles USNRCsoud ou USNRCmdb, ecart-type du decalage de RTNDT (°F) (28. pour js et 17. pour mdb)",
+                                   ),
 
-   LOGISTIC = BLOC ( condition = " Kind in ( 'Logistic', ) ",
+                           ), # Fin BLOC Parametres_USNRC 
 
-                       Alpha = SIMP ( statut = "o",
-                                      typ = "R",
-                                      max = 1,
-                                      fr = "Parametre Alpha de la loi",
-                                      ang = "Alpha parameter",
+                  nbectDRTNDT = SIMP ( statut = "o",
+                                       typ = "R",
+                                       defaut = "2.",
+                                       fr = "Nombre d ecart-type par rapport a la moyenne de DRTNDT",
                                       ),
 
-                       Beta = SIMP ( statut = "o",
-                                     typ = "R",
-                                     max = 1,
-                                     val_min = 0.,
-                                     fr = "Parametre Beta de la loi | Beta > = 0",
-                                     ang = "Beta parameter | Beta > = 0",
-                                     ),
+                           ), # Fin BLOC Parametres_FLUENCE
 
-   ), # Fin BLOC LOGISTIC
+#========================
+# 4.3 Modeles de tenacite
+#========================
 
+  MODELKIC = SIMP ( statut = "o",
+                    typ = "TXM",
+                   into = ( "RCC-M", "RCC-M_pal", "RCC-M_exp", "Houssin_RC", "Wallin", "REME", "ORNL", "Frama", "WEIB3", "WEIB2", "LOGWOLF", "WEIB_GEN" ),
+                    #defaut = "RCC-M",
+                    fr = "Modele de tenacite ",
+                   ),
 
+#====
+# Definition des parametres selon le modele de tenacité
+#====
 
-   MULTINOMIAL = BLOC ( condition = " Kind in ( 'MultiNomial', ) ",
-                         
-                         N = SIMP ( statut = "o",
-                                    typ = "I",
-                                    max = 1,
-                                    fr = "Dimension de la loi",
-                                    ang = "Distribution dimension",
-                                    ),
+  Parametres_RCCM = BLOC ( condition = " MODELKIC in ( 'RCC-M', ) ",
 
-                         # Il faut un vecteur P de taille N
-                         Values = SIMP ( statut = 'o',
-                                         typ = 'R',
-                                         max = '**',
-                                         ),
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
 
-   ), # Fin BLOC MULTINOMIAL
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
 
+                  KICPAL = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "195.",
+                                  fr = "palier deterministe de K1c quand modele RCC-M  (MPa(m^0.5)) ",
+                                 ),
 
+                  KICCDV = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.15",
+                                  fr = "coef de variation de la loi normale de K1c quand modele RCC-M ",
+                                 ),
 
-   NORMAL = BLOC ( condition = " Kind in ( 'Normal', ) ",
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
 
-                    Mu = SIMP ( statut = "o",
-                                typ = "R",
-                                max = 1,
-                                fr = "Parametre Mu de la loi",
-                                ang = "Mu parameter",
-                                ),
 
-                   Sigma = SIMP ( statut = "o",
+                  KIAPAL = SIMP ( statut = "o",
                                   typ = "R",
-                                  max = 1,
-                                  val_min = 0.,
-                                  fr = "Parametre Sigma de la loi | Sigma > 0",
-                                  ang = "Sigma parameter | Sigma > 0",
-                                  ),
+                                  defaut = "195.",
+                                  fr = "palier deterministe de K1a -tenacite a l'arret- quand modele RCC-M  (MPa(m^0.5)) ",
+                                 ),
 
-   ), # Fin BLOC NORMAL
+                  KIACDV = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.10",
+                                  fr = "coef de variation de la loi normale de K1a -tenacite a l'arret- quand modele RCC-M ",
+                                 ),
 
+                          ), # Fin BLOC Parametres_RCC-M
 
+  Parametres_RCCMpal = BLOC ( condition = " MODELKIC in ( 'RCC-M_pal', ) ",
 
-   POISSON = BLOC ( condition = " Kind in ( 'Poisson', ) ",
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
 
-                     Lambda = SIMP ( statut = "o",
-                                     typ = "R",
-                                     max = 1,
-                                     val_min = 0.,
-                                     fr = "Parametre Lambda de la loi | Lambda > 0",
-                                     ang = "Lambda parameter | Lambda > 0",
-                                     ),
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
 
-   ), # Fin BLOC POISSON
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
 
+                             ), # Fin BLOC Parametres_RCCMpal
 
+  Parametres_RCCMexp = BLOC ( condition = " MODELKIC in ( 'RCC-M_exp', ) ",
 
-   STUDENT = BLOC ( condition = " Kind in ( 'Student', ) ",
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
 
-                     Mu = SIMP ( statut = "o",
-                                 typ = "R",
-                                 max = 1,
-                                 fr = "Parametre Mu de la loi",
-                                 ang = "Mu parameter",
-                                 ),
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
 
-                     Nu = SIMP ( statut = "o",
-                                 typ = "R",
-                                 max = 1,
-                                 val_min = 2.,
-                                 fr = "Parametre Nu de la loi | V > = 2",
-                                 ang = "Nu parameter | V > = 2",
-                                 ),
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
 
-   ), # Fin BLOC STUDENT
+                             ), # Fin BLOC Parametres_RCCMexp
 
+  Parametres_Houssin_RC = BLOC ( condition = " MODELKIC in ( 'Houssin_RC', ) ",
 
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
 
-   TRIANGULAR = BLOC ( condition = " Kind in ( 'Triangular', ) ",
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
 
-                         A = SIMP ( statut = "o",
+                  nbectKIa = SIMP ( statut = "o",
                                     typ = "R",
-                                    max = 1,
-                                    fr = "Borne inferieure de la loi | A < = M < = B",
-                                    ang = "Lower bound | A < = M < = B",
-                                    ),
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
 
-                         M = SIMP ( statut = "o",
+                                ), # Fin BLOC Parametres_Houssin_RC
+      
+  Parametres_Wallin = BLOC ( condition = " MODELKIC in ( 'Wallin', ) ",
+                  T0WALLIN = SIMP ( statut = "o",
+                                    typ = "I",
+                                    defaut = "-27",
+                                    fr = "parametre T0 du modele Wallin (°C)",
+                                   ),
+
+                  nbectKIc = SIMP ( statut = "o",
                                     typ = "R",
-                                    max = 1,
-                                    fr = "Mode de la loi | A < = M < = B",
-                                    ang = "Mode | A < = M < = B",
-                                    ),
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
 
-                         B = SIMP ( statut = "o",
+                  fractKIc = SIMP ( statut = "o",
                                     typ = "R",
-                                    max = 1,
-                                    fr = "Borne superieure de la loi | A < = M < = B",
-                                    ang = "Upper bound | A < = M < = B",
-                                    ),
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
 
-   ), # Fin BLOC TRIANGULAR
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
 
+                            ), # Fin BLOC Parametres_Wallin
 
+  Parametres_REME = BLOC ( condition = " MODELKIC in ( 'REME', ) ",
 
-   TRUNCATEDNORMAL = BLOC ( condition = " Kind in ( 'TruncatedNormal', ) ",
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
 
-                             MuN = SIMP ( statut = "o",
-                                          typ = "R",
-                                          max = 1,
-                                          fr = "Parametre Mu de la loi",
-                                          ang = "Mu parameter",
-                                          ),
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
 
-                             SigmaN = SIMP ( statut = "o",
+                          ), # Fin BLOC Parametres_REME
+
+  Parametres_ORNL = BLOC ( condition = " MODELKIC in ( 'ORNL', ) ",
+
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+                          ), # Fin BLOC Parametres_ORNL
+
+  Parametres_Frama = BLOC ( condition = " MODELKIC in ( 'Frama', ) ",
+
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+                           ), # Fin BLOC Parametres_Frama
+
+  Parametres_WEIB3 = BLOC ( condition = " MODELKIC in ( 'WEIB3', ) ",
+
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+                           ), # Fin BLOC Parametres_WEIB3
+
+  Parametres_WEIB2 = BLOC ( condition = " MODELKIC in ( 'WEIB2', ) ",
+
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+                           ), # Fin BLOC Parametres_WEIB2
+
+  Parametres_LOGWOLF = BLOC ( condition = " MODELKIC in ( 'LOGWOLF', ) ",
+
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+                             ), # Fin BLOC Parametres_LOGWOLF
+
+  Parametres_WEIB_GEN = BLOC ( condition = " MODELKIC in ( 'WEIB_GEN',) ",
+                  A1 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "21.263",
+                              fr = "coef du coef a(T) d'une Weibull generale",
+                             ),
+                  A2 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "9.159",
+                              fr = "coef du coef a(T) d'une Weibull generale",
+                             ),
+                  A3 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.04057",
+                              fr = "coef du coef a(T) d'une Weibull generale",
+                             ),
+                  B1 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "17.153",
+                              fr = "coef du coef b(T) d'une Weibull generale",
+                             ),
+                  B2 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "55.089",
+                              fr = "coef du coef b(T) d'une Weibull generale",
+                             ),
+                  B3 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.0144",
+                              fr = "coef du coef b(T) d'une Weibull generale",
+                             ),
+                  C1 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "4.",
+                              fr = "coef du coef c(T) d'une Weibull generale",
+                             ),
+                  C2 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "coef du coef c(T) d'une Weibull generale",
+                             ),
+                  C3 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "coef du coef c(T) d'une Weibull generale",
+                             ),
+
+                  nbectKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+                  fractKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+                  nbectKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+                              ), # Fin BLOC Parametres_WEIB_GEN
+
+) # Fin OPER MODELES
+
+
+#==================
+# 5. Initialisation
+#==================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+INITIALISATION = OPER ( nom = "INITIALISATION",
+                        sd_prod = loi,
+                        op = 68,
+                        fr = "Initialisation", 
+
+  Liste_tempinit = SIMP ( statut = "o",
+                          typ = "R",
+                          max = "**",
+                          fr = "Profil radial de la temperature initiale dans la cuve (m) (°C) ",
+                         ),
+
+  Prolongation_tempsinit = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                 ),
+
+  Liste_sigmainit = SIMP ( statut = "o",
+                           typ = "R",
+                           max = "**",
+                           fr = "Profil radial des contraintes residuelles dans la cuve (m) (xx) (xx) (xx) ",
+                          ),
+
+  Prolongation_sigmainit = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                 ),
+
+  INSTINIT = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "-1.",
+                    fr = "Instant initial (s) ",
+                   ),
+
+) # Fin OPER INITIALISATION
+
+
+#==================================
+# 6. CARACTERISTIQUES DU REVETEMENT
+#==================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+REVETEMENT = OPER ( nom = "REVETEMENT",
+                    sd_prod = loi,
+                    op = 68,
+                    fr = "Caracteristiques du revetement", 
+
+  KTHREV = SIMP ( statut = "o",
+                  typ = "TXM",
+                  into = ( "ENTHALPIE", "CHALEUR",),
+                  #defaut = "CHALEUR",
+                  fr = "Options pour definir les caracteristiques du revetement ",
+                 ),
+
+  Parametres_ENTHALPIErev = BLOC ( condition = " KTHREV in ( 'ENTHALPIE', ) ",
+
+                  Liste_enthalpie_rev = SIMP ( statut = "o",
+                                               typ = "R",
+                                               max = "**",
+                                               fr = "Temperature (°C) / enthalpie (J/kg) ",
+                                              ),
+
+                  Prolongation_enthalpie_rev = SIMP ( statut = "o",
+                                                      typ = "TXM",
+                                                      into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                      fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                     ),
+
+                                  ), # Fin BLOC Parametres_ENTHALPIErev
+
+
+  Parametres_CHALEURrev = BLOC ( condition = " KTHREV in ( 'CHALEUR', ) ",
+
+                  Liste_chaleur_rev = SIMP ( statut = "o",
                                              typ = "R",
-                                             max = 1,
-                                             val_min = 0.,
-                                             fr = "Parametre SigmaN de la loi | SigmaN > 0",
-                                             ang = "SigmaN parameter | SigmaN> 0",
-                                             ),
+                                             max = "**",
+                                             fr = "Temperature (°C) / chaleur volumique (J/kg/K) ",
+                                            ),
 
-                             A = SIMP ( statut = "o",
-                                        typ = "R",
-                                        max = 1,
-                                        fr = "Borne inferieure de la loi | A < = B",
-                                        ang = "Lower bound | A < = B",
-                                        ),
+                  Prolongation_chaleur_rev = SIMP ( statut = "o",
+                                                    typ = "TXM",
+                                                    into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                    fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                   ),
 
-                             B = SIMP ( statut = "o",
-                                        typ = "R",
-                                        max = 1,
-                                        fr = "Borne superieure de la loi | A < = B",
-                                        ang = "Upper bound | A < = B",
-                                        ),
+                                ), # Fin BLOC Parametres_CHALEURrev
 
-   ), # Fin BLOC TRUNCATEDNORMAL
+   Liste_conductivite_rev = SIMP ( statut = "o",
+                                   typ = "R",
+                                   max = "**",
+                                   fr = "Temperature (°C) / conductivite thermique (W/m/°C) ",
+                                  ),
 
+   Prolongation_conductivite_rev = SIMP ( statut = "o",
+                                          typ = "TXM",
+                                          into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                          fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                         ),
 
+   Liste_young_rev = SIMP ( statut = "o",
+                            typ = "R",
+                            max = "**",
+                            fr = "Temperature (°C) / module d'Young (MPa) ",
+                           ),
 
-   UNIFORM = BLOC ( condition = " Kind in ( 'Uniform', ) ",
+   Prolongation_young_rev = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                  ),
 
-                     A = SIMP ( statut = "o",
-                                typ = "R",
-                                max = 1,
-                                fr = "Borne inferieure de la loi | A < = B",
-                                ang = "Lower bound | A < = B",
+   Liste_dilatation_rev = SIMP ( statut = "o",
+                                 typ = "R",
+                                 max = "**",
+                                 fr = "Temperature (°C) / coefficient de dilatation thermique (°C-1) ",
                                 ),
 
-                     B = SIMP ( statut = "o",
-                                typ = "R",
-                                max = 1,
-                                fr = "Borne superieure de la loi | A < = B",
-                                ang = "Upper bound | A < = B",
+   Prolongation_dilatation_rev = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                       ),
+
+   Liste_elasticite_rev = SIMP ( statut = "o",
+                                 typ = "R",
+                                 max = "**",
+                                 fr = "Temperature (°C) / limite d'elasticite (MPa) ",
                                 ),
 
-   ), # Fin BLOC UNIFORM
+   Prolongation_elasticite_rev = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                       ),
 
+   TREFREV = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "20.",
+                    fr = "Temperature de deformation nulle (°C) ",
+                   ),
 
+   TDETREV = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "287.",
+                    fr = "Temperature de definition du coefficient de dilatation thermique (°C) ",
+                   ),
 
-   USERDEFINED = BLOC ( condition = " Kind in ( 'UserDefined', ) ",
+   NUREV = SIMP ( statut = "o",
+                  typ = "R",
+                  defaut = "0.3",
+                  fr = "Coefficient de Poisson ",
+                 ),
 
-                           # Il faut definir une collection de couples ( x,p ) 
-                         Values = SIMP ( statut = 'o',
-                                         typ = 'R',
-                                         max = '**',
-                                         ),
+) # Fin OPER REVETEMENT
+
+
+#=====================================
+# 7. CARACTERISTIQUES DU METAL DE BASE
+#=====================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+METAL_BASE = OPER ( nom = "METAL_BASE",
+                    sd_prod = loi,
+                    op = 68,
+                    fr = "Caracteristiques du metal de base", 
+
+  KTHMDB = SIMP ( statut = "o",
+                  typ = "TXM",
+                  into = ( "ENTHALPIE", "CHALEUR",),
+                  #defaut = "CHALEUR",
+                  fr = "Options pour definir les caracteristiques du revetement ",
+                 ),
+
+  Parametres_ENTHALPIEmdb = BLOC ( condition = " KTHMDB in ( 'ENTHALPIE', ) ",
+
+                  Liste_enthalpie_mdb = SIMP ( statut = "o",
+                                               typ = "R",
+                                               max = "**",
+                                               fr = "Temperature (°C) / enthalpie (J/kg) ",
+                                              ),
 
-   ), # Fin BLOC USERDEFINED
+                  Prolongation_enthalpie_mdb = SIMP ( statut = "o",
+                                                      typ = "TXM",
+                                                      into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                      fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                     ),
+
+                                  ), # Fin BLOC Parametres_ENTHALPIEmdb
+
+
+  Parametres_CHALEURmdb = BLOC ( condition = " KTHMDB in ( 'CHALEUR', ) ",
 
+                  Liste_chaleur_mdb = SIMP ( statut = "o",
+                                             typ = "R",
+                                             max = "**",
+                                             fr = "Temperature (°C) / chaleur volumique (J/kg/K) ",
+                                            ),
+
+                  Prolongation_chaleur_mdb = SIMP ( statut = "o",
+                                                    typ = "TXM",
+                                                    into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                    fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                   ),
 
+                                ), # Fin BLOC Parametres_CHALEURmdb
 
-   WEIBULL = BLOC ( condition = " Kind in ( 'Weibull', ) ",
+   Liste_conductivite_mdb = SIMP ( statut = "o",
+                                   typ = "R",
+                                   max = "**",
+                                   fr = "Temperature (°C) / conductivite thermique (W/m/°C) ",
+                                  ),
 
-                     Settings = SIMP ( statut = "o",
+   Prolongation_conductivite_mdb = SIMP ( statut = "o",
                                           typ = "TXM",
-                                          max = 1,
-                                          into = ( "AlphaBeta", "MuSigma" ),
-                                          defaut = "AlphaBeta",
-                                          fr = "Parametrage de la loi weibull",
-                                          ang = "Weibull distribution parameter set",
+                                          into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                          fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                         ),
+
+   Liste_young_mdb = SIMP ( statut = "o",
+                            typ = "R",
+                            max = "**",
+                            fr = "Temperature (°C) / module d'Young (MPa) ",
+                           ),
+
+   Prolongation_young_mdb = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                  ),
+
+   Liste_dilatation_mdb = SIMP ( statut = "o",
+                                 typ = "R",
+                                 max = "**",
+                                 fr = "Temperature (°C) / coefficient de dilatation thermique (°C-1) ",
+                                ),
+
+   Prolongation_dilatation_mdb = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                       ),
+
+   TREFMDB = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "20.",
+                    fr = "Temperature de deformation nulle (°C) ",
+                   ),
+
+   TDETMDB = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "287.",
+                    fr = "Temperature de definition du coefficient de dilatation thermique (°C) ",
+                   ),
+
+   NUMDB = SIMP ( statut = "o",
+                  typ = "R",
+                  defaut = "0.3",
+                  fr = "Coefficient de Poisson ",
+                 ),
+
+) # Fin OPER METAL_BASE
+
+
+#================
+# 8. TRANSITOIRES
+#================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+TRANSITOIRE = OPER ( nom = "TRANSITOIRE",
+                     sd_prod = loi,
+                     op = 68,
+                     fr = "Description du transitoire", 
+
+  Liste_pression = SIMP ( statut = "o",
+                          typ = "R",
+                          max = "**",
+                          fr = "instant (s) / pression (MPa) ",
+                         ),
+
+  Prolongation_pression = SIMP ( statut = "o",
+                                 typ = "TXM",
+                                 into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                 fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                ),
+
+  TYPCLTH = SIMP ( statut = "o",
+                   typ = "TXM",
+                   into = ( "TEMP_IMPO", "FLUX_REP", "ECHANGE", "DEBIT", "TEMP_FLU"),
+                   #defaut = "ECHANGE",
+                   fr = "Type de condition thermique en paroi interne ",
+                  ),
+
+  Parametres_TEMP_IMPO = BLOC ( condition = " TYPCLTH in ( 'TEMP_IMPO', ) ",
+
+                  Liste_temp_impo = SIMP ( statut = "o",
+                                           typ = "R",
+                                           max = "**",
+                                           fr = "Instant (s) / Temperature imposee (°C) ",
                                           ),
 
-                     AlphaBeta_Parameters = BLOC ( condition = " Settings in ( 'AlphaBeta', ) ",
+                  Prolongation_temp_impo = SIMP ( statut = "o",
+                                                  typ = "TXM",
+                                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                 ),
 
-                                         Alpha = SIMP ( statut = "o",
-                                                        typ = "R",
-                                                        max = 1,
-                                                        val_min = 0.,
-                                                        fr = "Parametre Alpha de la loi | Alpha > 0",
-                                                        ang = "Alpha parameter | Alpha > 0",
-                                                        ),
+                               ), # Fin BLOC Parametres_TEMP_IMPO
 
-                                         Beta = SIMP ( statut = "o",
-                                                       typ = "R",
-                                                       max = 1,
-                                                       val_min = 0.,
-                                                       fr = "Parametre Beta de la loi | Beta > 0",
-                                                       ang = "Beta parameter | Beta > 0",
-                                                       ),
+  Parametres_FLUX_REP = BLOC ( condition = " TYPCLTH in ( 'FLUX_REP', ) ",
 
-                                         ), # Fin BLOC AlphaBeta_Parameters
+                  Liste_flux_rep = SIMP ( statut = "o",
+                                          typ = "R",
+                                          max = "**",
+                                          fr = "Instant (s) / Flux de chaleur impose (W/m2) ",
+                                         ),
 
+                  Prolongation_flux_rep = SIMP ( statut = "o",
+                                                 typ = "TXM",
+                                                 into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                 fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                ),
 
-                     MuSigma_Parameters = BLOC ( condition = " Settings in ( 'MuSigma', ) ",
+                              ), # Fin BLOC Parametres_FLUX_REP
 
-                                         Mu = SIMP ( statut = "o",
-                                                     typ = "R",
-                                                     max = 1,
-                                                     fr = "Parametre Mu de la loi",
-                                                     ang = "Mu parameter",
-                                                     ),
+  Parametres_ECHANGE = BLOC ( condition = " TYPCLTH in ( 'ECHANGE', ) ",
 
-                                         Sigma = SIMP ( statut = "o",
-                                                        typ = "R",
-                                                        max = 1,
-                                                        val_min = 0.,
-                                                        fr = "Parametre Sigma de la loi | Sigma > 0",
-                                                        ang = "Sigma parameter | Sigma > 0",
-                                                        ),
+                  Liste_echange_temp = SIMP ( statut = "o",
+                                              typ = "R",
+                                              max = "**",
+                                              fr = "Instant (s) / Temperature impose (°C) ",
+                                             ),
+
+                  Prolongation_echange_temp = SIMP ( statut = "o",
+                                                     typ = "TXM",
+                                                     into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                     fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                    ),
+
+                  Liste_echange_coef = SIMP ( statut = "o",
+                                              typ = "R",
+                                              max = "**",
+                                              fr = "Instant (s) / Coefficient d echange (W/m2/K) ",
+                                             ),
 
-                                         ), # Fin BLOC MuSigma_Parameters
+                  Prolongation_echange_coef = SIMP ( statut = "o",
+                                                     typ = "TXM",
+                                                     into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                     fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                    ),
+
+                             ), # Fin BLOC Parametres_ECHANGE
+
+  Parametres_DEBIT = BLOC ( condition = " TYPCLTH in ( 'DEBIT', ) ",
+
+                  Liste_debit = SIMP ( statut = "o",
+                                       typ = "R",
+                                       max = "**",
+                                       fr = "Instant (s) / Debit massique (kg/s) ",
+                                      ),
+
+                  Prolongation_debit = SIMP ( statut = "o",
+                                              typ = "TXM",
+                                              into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                              fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                             ),
+
+                  Liste_temp_injection = SIMP ( statut = "o",
+                                                typ = "R",
+                                                max = "**",
+                                                fr = "Instant (s) / Temperature d injection de securite  (°C) ",
+                                               ),
+
+                  Prolongation_temp_injection = SIMP ( statut = "o",
+                                                       typ = "TXM",
+                                                       into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                       fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                      ),
+
+                  DH = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Diametre hydraulique (m) ",
+                             ),
+
+                  SECTION = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "-2.",
+                                   fr = "Section espace annulaire (m2) ",
+                                  ),
+
+                  DELTA = SIMP ( statut = "o",
+                                 typ = "R",
+                                 defaut = "-2.",
+                                 fr = "Hauteur caracteristique convection naturelle (m) ",
+                                ),
+
+                  ALPHA_CF = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection forcee (-) ",
+                                   ),
 
-                     Gamma = SIMP ( statut = "o",
+                  ALPHA_CN = SIMP ( statut = "o",
                                     typ = "R",
-                                    max = 1,
-                                    fr = "Parametre Gamma",
-                                    ang = "Gamma parameter",
-                                    ),
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection naturelle (-) ",
+                                   ),
+
+                  EPS = SIMP ( statut = "o",
+                               typ = "R",
+                               defaut = "0.00001",
+                               fr = "Critere convergence relative (-) ",
+                              ),
+
+                  VM = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Volume de melange CREARE (m3) ",
+                             ),
+
+                  T0 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Temperature initiale CREARE (degC) ",
+                             ),
+
+                  SE = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Surface d'echange fluide/structure (m2) ",
+                             ),
+
+                           ), # Fin BLOC Parametres_DEBIT
 
-    ), # Fin BLOC WEIBULL
+  Parametres_TEMP_FLU = BLOC ( condition = " TYPCLTH in ( 'TEMP_FLU', ) ",
+
+                  Liste_temp_flu = SIMP ( statut = "o",
+                                          typ = "R",
+                                          max = "**",
+                                          fr = "Instant (s) / Debit massique (kg/s) ",
+                                         ),
 
-) # Fin OPER DISTRIBUTION
+                  Prolongation_temp_flu = SIMP ( statut = "o",
+                                                 typ = "TXM",
+                                                 into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                 fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                ),
 
+                  Liste_debit_injection = SIMP ( statut = "o",
+                                                 typ = "R",
+                                                 max = "**",
+                                                 fr = "Instant (s) / Debit d injection de securite (kg/s) ",
+                                                ),
+
+                  Prolongation_debit_injection = SIMP ( statut = "o",
+                                                        typ = "TXM",
+                                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclus, L=lineaire ",
+                                                       ),
+
+                  DH = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Diametre hydraulique (m) ",
+                             ),
+
+                  SECTION = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "-2.",
+                                   fr = "Section espace annulaire (m2) ",
+                                  ),
+
+                  DELTA = SIMP ( statut = "o",
+                                 typ = "R",
+                                 defaut = "-2.",
+                                 fr = "Hauteur caracteristique convection naturelle (m) ",
+                                ),
+
+                  ALPHA_CF = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection forcee (-) ",
+                                   ),
+
+                  ALPHA_CN = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection naturelle (-) ",
+                                   ),
 
+                  EPS = SIMP ( statut = "o",
+                               typ = "R",
+                               defaut = "0.00001",
+                               fr = "Critere convergence relative (-) ",
+                              ),
 
+                              ), # Fin BLOC Parametres_TEMP_FLU
 
+) # Fin OPER TRANSITOIRE
 
 
 #================================
diff --git a/Cuve2dg/Cuve2dg_Cata_V2.py b/Cuve2dg/Cuve2dg_Cata_V2.py
new file mode 100644 (file)
index 0000000..b9c202d
--- /dev/null
@@ -0,0 +1,2433 @@
+# -*- coding: utf-8 -*-
+
+# --------------------------------------------------
+# debut entete
+# --------------------------------------------------
+
+import Accas
+from Accas import *
+
+class loi ( ASSD ) : pass
+class variable ( ASSD ) : pass
+
+
+#CONTEXT.debug = 1
+JdC = JDC_CATA ( code = 'CUVE1D-DEFAILLGLOB',
+                 execmodul = None,
+                 regles = ( AU_MOINS_UN ('OPTIONS','DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'), 
+                            AU_MOINS_UN ( 'FIN' ), 
+                            A_CLASSER ( ('OPTIONS', 'DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'),'FIN')
+                          )
+                 ) # Fin JDC_CATA
+
+# --------------------------------------------------
+# fin entete
+# --------------------------------------------------
+
+
+
+
+
+
+#================================
+# 1. Definition des OPTIONS
+#================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+
+OPTIONS = OPER ( nom = "OPTIONS",
+                 sd_prod = loi,
+                 op = 68,
+                 fr = "Definitions des options", 
+
+#===
+# Liste des paramètres
+#===
+
+# INCRTPS
+  IncrementTemporel = SIMP ( statut = "o",
+                             typ = "I",
+                             defaut = "1",
+                            max = 1,
+                             #val_max = 100,
+                             fr = "Increment temporel (=1 pour calcul deterministe)",
+                             ),
+
+# DTPREC
+  IncrementMaxTemperature = SIMP ( statut = "o", 
+                                   typ = "R", 
+                                  defaut = "0.1", 
+                                  max = 1, 
+                                  val_max = 1., 
+                                  fr = "Increment maximum d'evolution de la temperature par noeud et par instant (°C)",
+                                  ),
+
+# DTARCH
+  IncrementMaxTempsAffichage = SIMP ( statut = "o", 
+                                      typ = "R", 
+                                     defaut = "1000.", 
+                                     max = 1, 
+                                     val_max = 1000., 
+                                     fr = "Increment maximum de temps pour l'affichage (s)",
+                                     ),
+
+# NBO
+# Question : NBO depend-il de TYPGEOM ??
+  NombreNoeudsMaillage = SIMP ( statut = "o", 
+                                typ = "R", 
+                               max=1, 
+                               val_max = 1000., 
+                               fr = "Nombre de noeuds a considerer dans le maillage interne",
+                               ),
+
+# 
+  ListeInstants = SIMP ( statut = "o",
+                         typ = "R",
+                         max = "**",
+                         fr = "Liste des instants ",
+                         ),
+
+) # Fin OPER OPTIONS
+
+#================================
+# 2. Caracteristiques du DEFAUT
+#================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+DEFAUT = OPER ( nom = "DEFAUT",
+                sd_prod = loi,
+                op = 68,
+                fr = "Caracteristiques du defaut", 
+
+#===
+# Liste des paramètres
+#===
+
+# TYPEDEF
+  TypeInitial = SIMP ( statut = "o", typ = "TXM",
+                       into = ( "Sous Revetement", # DSR
+                                "Debouchant", # DD
+                                 ),
+                       fr = "Type initial du defaut : sous revetement ou debouchant",
+                       ),
+
+#====
+# Definition des parametres selon le type du defaut
+#====
+
+  SousRevetement = BLOC ( condition = " TypeInitial in ( 'Sous Revetement', ) ",
+
+    # ORIEDEF
+    Orientation = SIMP ( statut = "o",
+                         typ = "TXM",
+                         into = ( "longitudinale", # LONGITUD
+                                 "circonferentielle" ), # CIRCONF
+                         fr = "Orientation du defaut",
+                         ),
+
+    # PROFDEF
+    ProfondeurRadiale = SIMP ( statut = "o",
+                               typ = "R",
+                               #defaut = "0.006",
+                               max = 1,
+                               val_max = 1.,
+                               fr = "Profondeur radiale du defaut (m)",
+                               ),
+
+    # OPTLONG
+    OptionCalculLongueur = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( "Valeur", # VALEUR
+                                          "Relation lineaire avec la longueur" ), # RAPPORT
+                                  #defaut = "VALEUR",
+                                  fr = "Option pour caracteriser la longueur du defaut : soit par valeur, soit par un rapport LONG/PROF",
+                                  ),
+
+    Option_Valeur = BLOC ( condition = "OptionCalculLongueur in ( 'Valeur', ) ",
+
+      # LONGDEF
+      Longueur = SIMP ( statut = "o",
+                        typ = "R",
+                                             #defaut = "0.060",
+                                             max = 1,
+                                             val_max = 1.,
+                                             fr = "Longueur du defaut sous revetement (m)",
+                                            ),
+
+    ), # Fin BLOC Option_Valeur
+
+    Option_Rapport = BLOC ( condition = "OptionCalculLongueur in ( 'Relation lineaire avec la longueur', ) ",
+
+        # LONGSURPROF
+        LongueurSurProfondeur = SIMP ( statut = "o",
+                                                 typ = "R",
+                                                 #defaut = "6.",
+                                                 max = 1,
+                                                 val_max = 100.,
+                                                 fr = "Rapport longueur/profondeur du defaut sous revetement",
+                                                ),
+
+    ), # Fin BLOC Option_Rapport
+
+    # DECADEF
+    DecalageRadial = SIMP ( statut = "o",
+                                   typ = "R",
+                                   #defaut = "-0.00001",
+                                   fr = "Decalage radial du defaut sous revetement (m)",
+                                  ),
+
+    # ANGLDEF
+    Azimut = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.",
+                                   fr = "Position angulaire du defaut dans le cylindre de cuve (degres)",
+                                  ),
+
+    # ALTIDEF
+    Altitude = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "2.",
+                                  val_min = 0.,
+                                  val_max = 4.,
+                                   fr = "Altitude du defaut sur le cylindre de cuve (m)",
+                                  ),
+
+    # POINDEF
+    Pointe = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( "A", "B" ),
+                                   defaut = "A",
+                                   fr = "Choix du point considere du defaut sous revetement",
+                                  ),
+
+    # ARRETFISSURE
+    ArretDeFissure = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( "OUI", "NON" ),
+                                        defaut = "NON",
+                                        fr = "Prise en compte de l arret de fissure",
+                                       ),
+
+    # INCRDEF
+    IncrementTailleFissure = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.005",
+                                   fr = "Increment de la taille de fissure pour la propagation (m)",
+                                  ),
+
+    # CORRECPLASTIC
+    CorrectionPlastiqueBeta = SIMP ( statut = "o",
+                                         typ = "TXM",
+                                         into = ( "OUI", "NON" ),
+                                         defaut = "NON",
+                                         fr = "Prise en compte de la correction plastique BETA dans le calcul du facteur d intensite de contraintes",
+                                        ),
+
+  ), # Fin BLOC SousRevetement
+
+  Debouchant = BLOC ( condition = " TypeInitial in ( 'Debouchant', ) ",
+
+    # ORIEDEF
+    Orientation = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( "longitudinale", # LONGITUD
+                                           "circonferentielle" ), # CIRCONF
+                                   #defaut = "LONGITUD",
+                                   fr = "Orientation du defaut : longitudinale ou circonferentielle",
+                                  ),
+
+    # PROFDEF
+    ProfondeurRadiale = SIMP ( statut = "o",
+                                   typ = "R",
+                                   #defaut = "0.006",
+                                   max = 1,
+                                   val_max = 1.,
+                                   fr = "Profondeur radiale du defaut (m)",
+                                  ),
+
+    # ANGLDEF
+    Azimut = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.",
+                                   fr = "Position angulaire du defaut dans le cylindre de la cuve (degres)",
+                                  ),
+
+    # ALTIDEF
+    Altitude = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "2.",
+                                   fr = "Altitude du defaut sur le cylindre de la cuve (m)",
+                                  ),
+
+    # ARRETFISSURE
+    ArretDeFissure = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( "OUI", "NON" ),
+                                        defaut = "NON",
+                                        fr = "Prise en compte de l arret de fissure",
+                                       ),
+
+    # INCRDEF
+    IncrementTailleFissure = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "0.005",
+                                   fr = "Increment de la taille de fissure pour la propagation (m)",
+                                  ),
+
+    # IRWIN
+    CorrectionPlastiqueIrwin = SIMP ( statut = "o",
+                                 typ = "TXM",
+                                 into = ( "OUI", "NON" ),
+                                 defaut = "NON",
+                                 fr = "Prise en compte de la correction plastique d'Irwin dans le calcul du facteur d intensite de contraintes",
+                                ),
+
+    # CORRECPLASTIC
+    CorrectionPlastiqueBeta = SIMP ( statut = "o",
+                                         typ = "TXM",
+                                         into = ( "OUI", "NON" ),
+                                         defaut = "NON",
+                                         fr = "Prise en compte de la correction plastique BETA dans le calcul du facteur d intensite de contraintes",
+                                        ),
+
+  ), # Fin BLOC debouchant
+
+) # Fin OPER DEFAUT
+
+
+#================================
+# 3. Caracteristiques de la CUVE
+#================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+CUVE = OPER (nom = "CUVE",
+             sd_prod = loi,
+             op = 68,
+             fr = "Caracteristiques de la cuve", 
+
+#===
+# Liste des paramètres
+#===
+
+  # TYPEGEOM
+  TraitementGeometrie = SIMP ( statut = "o",
+                    typ = "TXM",
+                   into = ( "Topologie", # GEOMETRIE 
+                            "Maillage"), # MAILLAGE
+                    #defaut = "geometrie",
+                    fr = "Choix de la definition de la geometrie d'une cuve",
+                   ),
+
+
+#====
+# Definition des parametres selon le type de traitement de la geometrie
+#====
+
+  Geometrie = BLOC ( condition = " TraitementGeometrie in ( 'Topologie', ) ",
+
+    # RINT
+    RayonInterne = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "1.994",
+                                fr = "Rayon interne de la cuve (m)",
+                               ),
+
+    # REXT
+    RayonExterne = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "2,2015",
+                                fr = "Rayon externe de la cuve (m)",
+                               ),
+
+    # LREV
+    EpaisseurRevetement = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "0.0075",
+                                fr = "Epaisseur du revetement (m)",
+                               ),
+
+    # LIGMIN
+    LigamentExterneMin = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.75",
+                                  fr = "Ligament externe minimal avant rupture (% de l'epaisseur de cuve)",
+                                 ),
+
+  ), # Fin BLOC  Geometrie
+
+  Maillage = BLOC ( condition = " TraitementGeometrie in ( 'Maillage', ) ",
+
+    # Question : quel rapport avec NBO ??
+    Liste_abscisses = SIMP ( statut = "o",
+                                           typ = "R",
+                                           max = "**",
+                                           fr = "Liste des abscisses (m) A FAIRE",
+                                          ),
+  ), # Fin BLOC Maillage
+
+) # Fin OPER CUVE
+
+#====================================================
+# 4. Modeles de fluence, d'irradiation et de tenacite
+#====================================================
+
+#=======================
+# 4.1 Modeles de fluence
+#=======================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+MODELES = OPER ( nom = "MODELES",
+                 sd_prod = loi,
+                 op = 68,
+                 fr = "Modeles de fluence, d'irradiation et de tenacite", 
+
+
+#===
+# Liste des paramètres
+#===
+
+  # MODELFLUENCE
+  Fluence = SIMP ( statut = "o",
+                        typ = "TXM",
+                       into = ( "Exponentiel sans revetement k=9.7 (Reglementaire)", # Reglementaire
+                                "Exponentiel sans revetement k=12.7 (France)", # France
+                                "Exponentiel sans revetement k=0. (ValeurImposee)", # ValeurImposee
+                                "Donnees francaises du palier CPY (SDM)", # SDM 
+                                "Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)", # GrandeDev
+                                "Regulatory Guide 1.99 rev 2 (USNRC)", # USNRC
+                                "Dossier 900 MWe AP9701 rev 2 (REV_2)", # REV_2
+                                "Lissage du modele ajuste (SDM_Lissage)", # SDM_Lissage 
+                                "Grand developpement (GD_Cuve)"), # GD_Cuve
+                        #defaut = "Reglementaire",
+                        fr = "Modele d'attenuation de la fluence dans l'epaisseur de la cuve",
+                       ),
+
+
+#====
+# Definition des parametres selon le modele de fluence
+#====
+
+  Reglementaire = BLOC ( condition = " Fluence in ( 'Exponentiel sans revetement k=9.7 (Reglementaire)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale en surface interne assimilee par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+  ), # Fin BLOC Reglementaire
+
+  France = BLOC ( condition = " Fluence in ( 'Exponentiel sans revetement k=12.7 (France)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+    # KPFRANCE = SIMP ( statut = "o",
+    KPFrance = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "12.7",
+                                    fr = "Parametre exponentiel du modele France",
+                                   ),
+
+  ), # Fin BLOC France
+
+  ValeurImposee = BLOC ( condition = " Fluence in ( 'Exponentiel sans revetement k=0. (ValeurImposee)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+  ), # Fin BLOC ValeurImposee
+
+  SDM = BLOC ( condition = " Fluence in ( 'Donnees francaises du palier CPY' (SDM), ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+  ), # Fin BLOC SDM
+
+  USNRC = BLOC ( condition = " Fluence in ( 'Regulatory Guide 1.99 rev 2 (USNRC)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+    KPUS = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "9.4488",
+                                fr = "Parametre exponentiel du modele US",
+                               ),
+
+  ), # Fin BLOC USNRC
+
+  REV_2 = BLOC ( condition = " Fluence in ( 'Dossier 900 MWe AP9701 rev 2 (REV_2)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+  ), # Fin BLOC REV_2
+
+  SDM_Lissage = BLOC ( condition = " Fluence in ( 'Lissage du modele ajuste (SDM_Lissage)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                                ),
+
+  ), # Fin BLOC SDM_Lissage
+
+  GrandeDev = BLOC ( condition = " Fluence in ( 'Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+  ), # Fin BLOC GrandeDev
+
+  GD_Cuve = BLOC ( condition = " Fluence in ( 'Grand developpement (GD_Cuve)', ) ",
+
+    # fmax
+    FluenceMax = SIMP ( statut = "o",
+                                typ = "R",
+                                defaut = "6.5",
+                                fr = "Fluence maximale assimilee par la cuve (n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95",
+                               ),
+
+    # COEFFLUENCE1
+    FluenceAzimut0 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "5.8",
+                                        fr = "Fluence a l'azimut 0 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE2
+    FluenceAzimut5 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "5.48",
+                                        fr = "Fluence a l'azimut 5 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE3
+    FluenceAzimut10 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "4.46",
+                                        fr = "Fluence a l'azimut 10 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE4
+    FluenceAzimut15 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "3.41",
+                                        fr = "Fluence a l'azimut 15 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE5
+    FluenceAzimut20 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "3.37",
+                                        fr = "Fluence a l'azimut 20 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE6
+    FluenceAzimut25 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "3.16",
+                                        fr = "Fluence a l'azimut 25 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE7
+    FluenceAzimut30 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "2.74",
+                                        fr = "Fluence a l'azimut 30 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE8
+    FluenceAzimut35 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "2.25",
+                                        fr = "Fluence a l'azimut 35 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE9
+    FluenceAzimut40 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "1.89",
+                                        fr = "Fluence a l'azimut 40 (10^19 n/cm)",
+                                       ),
+
+    # COEFFLUENCE10
+    FluenceAzimut45 = SIMP ( statut = "o",
+                                        typ = "R",
+                                        defaut = "1.78",
+                                        fr = "Fluence a l'azimut 45 (10^19 n/cm)",
+                                       ),
+
+  ), # Fin BLOC GD_Cuve
+
+#==========================
+# 4.2 Modeles d'irradiation
+#==========================
+
+  # TYPIRR
+  Irradiation = SIMP ( statut = "o",
+                  typ = "TXM",
+                 into = ( "RTndt de la cuve a l instant de l analyse", # RTNDT 
+                          "Modele d irradiation" ), # FLUENCE
+                  fr = "Type d'irradiation ",
+                 ),
+
+#====
+# Definition des parametres selon le type d'irradiation
+#====
+
+  RTndt = BLOC ( condition = " Irradiation in ( 'RTndt de la cuve a l instant de l analyse', ) ",
+    RTNDT = SIMP ( statut = "o",
+                                 typ = "R",
+                                 defaut = "73.",
+                                 fr = "RTNDT finale (°C)",
+                                ),
+
+  ), # Fin BLOC RTndt
+
+  Modele = BLOC ( condition = " Irradiation in ( 'Modele d irradiation', ) ",
+    # MODELIRR
+    ModeleIrradiation = SIMP ( statut = "o",
+                                    typ = "TXM",
+                                   into = ( "Metal de Base : formule de FIM/FIS Houssin", # HOUSSIN 
+                                            "Metal de Base : formule de FIM/FIS Persoz", # PERSOZ
+                                            "Metal de Base : formule de FIM/FIS Lefebvre", # LEFEBVRE
+                                            "Metal de Base : Regulatory Guide 1.00 rev 2", # USNRCmdb
+                                            "Joint Soude : formulation de FIM/FIS Brillaud", # BRILLAUD
+                                            "Joint Soude : Regulatory Guide 1.00 rev 2" ), # USNRCsoud
+                                    fr = "Modele d'irradiation pour virole ou joint soude",
+                                   ),
+
+    # CU
+    TeneurCuivre = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "Teneur en cuivre (%)",
+                             ),
+
+    # NI
+    TeneurNickel = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "Teneur en nickel (%)",
+                             ),
+
+    # P
+    TeneurPhosphore = SIMP ( statut = "o",
+                             typ = "R",
+                             defaut = "0.",
+                             fr = "Teneur en phosphore (%)",
+                            ),
+
+    Parametres_USNRC = BLOC ( condition = " ModeleIrradiation in ( 'Joint Soude : Regulatory Guide 1.00 rev 2', 'Metal de Base : Regulatory Guide 1.00 rev 2' , ) ",
+      # RTimoy
+      MoyenneRTndt = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.",
+                                  fr = "Moyenne de RTNDT : virole C1 de cuve Chinon : mdb=>-17.°C et js=>42.°C (HT-56/05/038 : p.52)",
+                                 ),
+
+      # RTicov
+      CoefVariationRTndt = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.",
+                                  fr = "Coefficient de variation de la RTNDT initiale",
+                                 ),
+
+      # USectDRT
+      EcartTypeRTndt = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "28.",
+                                    fr = "Ecart-type du decalage de RTNDT (°F) (28. pour js et 17. pour mdb)",
+                                   ),
+
+      # nbectDRTNDT
+      NombreEcartTypeRTndt = SIMP ( statut = "o",
+                                       typ = "R",
+                                       defaut = "2.",
+                                       fr = "Nombre d ecart-type par rapport a la moyenne de DRTNDT",
+                                      ),
+
+    ), # Fin BLOC Parametres_USNRC 
+
+  ), # Fin BLOC Modele
+
+#========================
+# 4.3 Modeles de tenacite
+#========================
+
+  # MODELKIC
+  Tenacite = SIMP ( statut = "o",
+                    typ = "TXM",
+                   into = ( "RCC-M/ASME coefficient=2", # RCC-M
+                            "RCC-M/ASME coefficient=2.33 (Houssin)", # Houssin_RC
+                            "RCC-M/ASME avec KI=KIpalier", # RCC-M_pal
+                            "RCC-M/ASME avec KI~exponentiel", # RCC-M_exp
+                            "Weibull basee sur la master cuve", # Wallin
+                            "Weibull basee sur la master cuve (REME)", # REME
+                            "Weibull n°1 (etude ORNL)", # ORNL
+                            "Weibull n°2", # WEIB2
+                            "Weibull n°3", # WEIB3
+                            "Weibull generalisee", # WEIB_GEN
+                            "Exponentielle n°1 (Frama)", # Frama
+                            "Exponentielle n°2 (LOGWOLF)" ), # LOGWOLF
+                    fr = "Modele de tenacite ",
+                   ),
+
+#====
+# Definition des parametres selon le modele de tenacité
+#====
+
+  RCCM_delta2 = BLOC ( condition = " Tenacite in ( 'RCC-M/ASME coefficient=2', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # KICPAL
+    PalierDuctile_KIc = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "195.",
+                                  fr = "palier deterministe de K1c quand modele RCC-M  (MPa(m^0.5)) ",
+                                 ),
+
+    # KICCDVD = SIMP ( statut = "o",
+    CoefficientVariation_KIc = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.15",
+                                  fr = "coef de variation de la loi normale de K1c quand modele RCC-M ",
+                                 ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+
+    # KIAPAL
+    PalierDuctile_KIa = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "195.",
+                                  fr = "palier deterministe de K1a -tenacite a l'arret- quand modele RCC-M  (MPa(m^0.5)) ",
+                                 ),
+
+    # KIACDV
+    CoefficientVariation_KIa = SIMP ( statut = "o",
+                                  typ = "R",
+                                  defaut = "0.10",
+                                  fr = "coef de variation de la loi normale de K1a -tenacite a l'arret- quand modele RCC-M ",
+                                 ),
+
+  ), # Fin BLOC Parametres_RCC-M
+
+  RCCM_KIpal = BLOC ( condition = " Tenacite in ( 'RCC-M/ASME avec KI=KIpalier', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_RCCMpal
+
+  RCCM_KIexp = BLOC ( condition = " Tenacite in ( 'RCC-M/ASME avec KI~exponentiel', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_RCCMexp
+
+  RCCM_delta233 = BLOC ( condition = " Tenacite in ( 'RCC-M/ASME coefficient=2.33 (Houssin)', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_Houssin_RC
+      
+  Weibull_MasterCuve = BLOC ( condition = " Tenacite in ( 'Weibull basee sur la master cuve', ) ",
+    # T0WALLIN
+    Temperature_KIc100 = SIMP ( statut = "o",
+                                    typ = "I",
+                                    defaut = "-27",
+                                    fr = "parametre T0 du modele Wallin (°C)",
+                                   ),
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_Wallin
+
+  Weibull_MasterCuveREME = BLOC ( condition = " Tenacite in ( 'Weibull basee sur la master cuve (REME)', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_REME
+
+  Weibull1_ORNL = BLOC ( condition = " Tenacite in ( 'Weibull n°1 (etude ORNL)', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_ORNL
+
+  Exponentielle1_Frama = BLOC ( condition = " Tenacite in ( 'Exponentielle n°1 (Frama)', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_Frama
+
+  Weibull3 = BLOC ( condition = " Tenacite in ( 'Weibull n°3', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_WEIB3
+
+  Weibull2 = BLOC ( condition = " Tenacite in ( 'Weibull n°2', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_WEIB2
+
+  Exponentielle2_LOGWOLF = BLOC ( condition = " Tenacite in ( 'Exponentielle n°2 (LOGWOLF)', ) ",
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_LOGWOLF
+
+  Weibull_Generalisee = BLOC ( condition = " Tenacite in ( 'Weibull generalisee',) ",
+    # A1
+    A1 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "21.263",
+                              fr = "coef du coef a(T) d'une Weibull generale",
+                             ),
+    # A2
+    A2 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "9.159",
+                              fr = "coef du coef a(T) d'une Weibull generale",
+                             ),
+    # A3
+    A3 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.04057",
+                              fr = "coef du coef a(T) d'une Weibull generale",
+                             ),
+    # B1
+    B1 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "17.153",
+                              fr = "coef du coef b(T) d'une Weibull generale",
+                             ),
+    # B2
+    B2 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "55.089",
+                              fr = "coef du coef b(T) d'une Weibull generale",
+                             ),
+    # B3
+    B3 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.0144",
+                              fr = "coef du coef b(T) d'une Weibull generale",
+                             ),
+    # C1
+    C1 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "4.",
+                              fr = "coef du coef c(T) d'une Weibull generale",
+                             ),
+    # C2
+    C2 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "coef du coef c(T) d'une Weibull generale",
+                             ),
+    # C3
+    C3 = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "0.",
+                              fr = "coef du coef c(T) d'une Weibull generale",
+                             ),
+
+    # nbectKIc
+    NbEcartType_MoyKIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIc (nb sigma) : det = -2 ",
+                                   ),
+
+    # fractKIc
+    Fractile_KIc = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "5.",
+                                    fr = "Valeur caracteristique de KIc exprimee en ordre de fractile (%) ",
+                                   ),
+
+    # nbectKIa
+    NbEcartType_MoyKIa = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "-2.",
+                                    fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ",
+                                   ),
+
+  ), # Fin BLOC Parametres_WEIB_GEN
+
+) # Fin OPER MODELES
+
+
+#==================
+# 5. Initialisation
+#==================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+INITIALISATION = OPER ( nom = "INITIALISATION",
+                        sd_prod = loi,
+                        op = 68,
+                        fr = "Initialisation : instant initial, profils radiaux de temperature et contraintes", 
+
+  ProfilRadialTemperatureInitiale = SIMP ( statut = "o",
+                          typ = "R",
+                          max = "**",
+                          fr = "Profil radial de la temperature initiale dans la cuve (m) (°C) ",
+                         ),
+
+  Prolongation_tempinit = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                 ),
+
+  ProfilRadialContraintesInitiales = SIMP ( statut = "o",
+                           typ = "R",
+                           max = "**",
+                           fr = "Profil radial des contraintes residuelles dans la cuve (m) (xx) (xx) (xx) ",
+                          ),
+
+  Prolongation_sigmainit = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                 ),
+
+  # INSTINIT
+  InstantInitialisation = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "-1.",
+                    fr = "Instant initial (s) ",
+                   ),
+
+) # Fin OPER INITIALISATION
+
+
+#==================================
+# 6. CARACTERISTIQUES DU REVETEMENT
+#==================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+REVETEMENT = OPER ( nom = "REVETEMENT",
+                    sd_prod = loi,
+                    op = 68,
+                    fr = "Caracteristiques du revetement", 
+
+  # KTHREV
+  OptionREV = SIMP ( statut = "o",
+                  typ = "TXM",
+                  into = ( "ENTHALPIE", "CHALEUR",),
+                  #defaut = "CHALEUR",
+                  fr = "Option pour definir les caracteristiques du revetement ",
+                 ),
+
+  Parametres_ENTHALPIErev = BLOC ( condition = " OptionREV in ( 'ENTHALPIE', ) ",
+
+    Liste_enthalpie_rev = SIMP ( statut = "o",
+                                               typ = "R",
+                                               max = "**",
+                                               fr = "Temperature (°C) / enthalpie (J/kg) ",
+                                              ),
+
+    Prolongation_enthalpie_rev = SIMP ( statut = "o",
+                                                      typ = "TXM",
+                                                      into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                      fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                     ),
+
+  ), # Fin BLOC Parametres_ENTHALPIErev
+
+
+  Parametres_CHALEURrev = BLOC ( condition = " OptionREV in ( 'CHALEUR', ) ",
+
+    Liste_chaleur_rev = SIMP ( statut = "o",
+                                             typ = "R",
+                                             max = "**",
+                                             fr = "Temperature (°C) / chaleur volumique (J/kg/K) ",
+                                            ),
+
+    Prolongation_chaleur_rev = SIMP ( statut = "o",
+                                                    typ = "TXM",
+                                                    into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                    fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                   ),
+
+  ), # Fin BLOC Parametres_CHALEURrev
+
+  Liste_conductivite_rev = SIMP ( statut = "o",
+                                   typ = "R",
+                                   max = "**",
+                                   fr = "Temperature (°C) / conductivite thermique (W/m/°C) ",
+                                  ),
+
+  Prolongation_conductivite_rev = SIMP ( statut = "o",
+                                          typ = "TXM",
+                                          into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                          fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                         ),
+
+  Liste_young_rev = SIMP ( statut = "o",
+                            typ = "R",
+                            max = "**",
+                            fr = "Temperature (°C) / module d'Young (MPa) ",
+                           ),
+
+  Prolongation_young_rev = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                  ),
+
+  Liste_dilatation_rev = SIMP ( statut = "o",
+                                 typ = "R",
+                                 max = "**",
+                                 fr = "Temperature (°C) / coefficient de dilatation thermique (°C-1) ",
+                                ),
+
+  Prolongation_dilatation_rev = SIMP ( statut = "o",
+                                       typ = "TXM",
+                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                       ),
+
+  Liste_elasticite_rev = SIMP ( statut = "o",
+                                 typ = "R",
+                                 max = "**",
+                                 fr = "Temperature (°C) / limite d'elasticite (MPa) ",
+                                ),
+
+  Prolongation_elasticite_rev = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                       ),
+
+  # TREFREV
+  TemperatureDeformationNulleREV = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "20.",
+                    fr = "Temperature de deformation nulle (°C) ",
+                   ),
+
+  # TDETREV
+  TemperatureDeDefinitionPourCoefThermiqueREV = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "287.",
+                    fr = "Temperature de definition du coefficient de dilatation thermique (°C) ",
+                   ),
+
+  # NUREV
+  CoefficientPoissonREV = SIMP ( statut = "o",
+                  typ = "R",
+                  defaut = "0.3",
+                  fr = "Coefficient de Poisson ",
+                 ),
+
+) # Fin OPER REVETEMENT
+
+
+#=====================================
+# 7. CARACTERISTIQUES DU METAL DE BASE
+#=====================================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+METAL_BASE = OPER ( nom = "METAL_BASE",
+                    sd_prod = loi,
+                    op = 68,
+                    fr = "Caracteristiques du metal de base", 
+
+  # KTHMDB
+  OptionMDB = SIMP ( statut = "o",
+                  typ = "TXM",
+                  into = ( "ENTHALPIE", "CHALEUR",),
+                  #defaut = "CHALEUR",
+                  fr = "Option pour definir les caracteristiques du revetement ",
+                 ),
+
+  Parametres_ENTHALPIEmdb = BLOC ( condition = " OptionMDB in ( 'ENTHALPIE', ) ",
+
+    Liste_enthalpie_mdb = SIMP ( statut = "o",
+                                               typ = "R",
+                                               max = "**",
+                                               fr = "Temperature (°C) / enthalpie (J/kg) ",
+                                              ),
+
+    Prolongation_enthalpie_mdb = SIMP ( statut = "o",
+                                                      typ = "TXM",
+                                                      into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                      fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                     ),
+
+  ), # Fin BLOC Parametres_ENTHALPIEmdb
+
+
+  Parametres_CHALEURmdb = BLOC ( condition = " OptionMDB in ( 'CHALEUR', ) ",
+
+    Liste_chaleur_mdb = SIMP ( statut = "o",
+                                             typ = "R",
+                                             max = "**",
+                                             fr = "Temperature (°C) / chaleur volumique (J/kg/K) ",
+                                            ),
+
+    Prolongation_chaleur_mdb = SIMP ( statut = "o",
+                                                    typ = "TXM",
+                                                    into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                    fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                   ),
+
+  ), # Fin BLOC Parametres_CHALEURmdb
+
+  Liste_conductivite_mdb = SIMP ( statut = "o",
+                                   typ = "R",
+                                   max = "**",
+                                   fr = "Temperature (°C) / conductivite thermique (W/m/°C) ",
+                                  ),
+
+  Prolongation_conductivite_mdb = SIMP ( statut = "o",
+                                          typ = "TXM",
+                                          into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                          fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                         ),
+
+  Liste_young_mdb = SIMP ( statut = "o",
+                            typ = "R",
+                            max = "**",
+                            fr = "Temperature (°C) / module d'Young (MPa) ",
+                           ),
+
+  Prolongation_young_mdb = SIMP ( statut = "o",
+                                   typ = "TXM",
+                                   into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                   fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                  ),
+
+  Liste_dilatation_mdb = SIMP ( statut = "o",
+                                 typ = "R",
+                                 max = "**",
+                                 fr = "Temperature (°C) / coefficient de dilatation thermique (°C-1) ",
+                                ),
+
+  Prolongation_dilatation_mdb = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                       ),
+
+  # TREFMDB
+  TemperatureDeformationNulleMDB = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "20.",
+                    fr = "Temperature de deformation nulle (°C) ",
+                   ),
+
+  # TDETMDB
+  TemperatureDeDefinitionPourCoefDilatationThermiqueMDB = SIMP ( statut = "o",
+                    typ = "R",
+                    defaut = "287.",
+                    fr = "Temperature de definition du coefficient de dilatation thermique (°C) ",
+                   ),
+
+  # NUMDB
+  CoefficientPoissonMDB = SIMP ( statut = "o",
+                  typ = "R",
+                  defaut = "0.3",
+                  fr = "Coefficient de Poisson ",
+                 ),
+
+) # Fin OPER METAL_BASE
+
+
+#================
+# 8. TRANSITOIRES
+#================
+
+# Nota : les variables de type OPER doivent etre en majuscules !
+# Nota : les variables de type OPER doivent etre de premier niveau (pas imbriquees dans un autre type)
+TRANSITOIRE = OPER ( nom = "TRANSITOIRE",
+                     sd_prod = loi,
+                     op = 68,
+                     fr = "Description du transitoire thermohydraulique", 
+
+  Liste_pression = SIMP ( statut = "o",
+                          typ = "R",
+                          max = "**",
+                          fr = "instant (s) / pression (MPa) ",
+                         ),
+
+  Prolongation_pression = SIMP ( statut = "o",
+                                 typ = "TXM",
+                                 into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                 fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                ),
+
+  # TYPCLTH
+  TypeConditionLimiteThermique = SIMP ( statut = "o",
+                                        typ = "TXM",
+                                        into = ( "TEMP_IMPO", # TEMP_IMPO 
+                                                "FLUX_REP", # FLUX_REP 
+                                                "ECHANGE", # ECHANGE 
+                                                "DEBIT", # DEBIT 
+                                                "TEMP_FLU"), # TEMP_FLU
+                                        #defaut = "ECHANGE",
+                                        fr = "Type de condition thermique en paroi interne ",
+                                        ),
+
+  Parametres_TEMP_IMPO = BLOC ( condition = " TypeConditionLimiteThermique in ( 'TEMP_IMPO', ) ",
+
+    Liste_temp_impo = SIMP ( statut = "o",
+                                      typ = "R",
+                                           max = "**",
+                                           fr = "Instant (s) / Temperature imposee (°C) ",
+                                          ),
+
+     Prolongation_temp_impo = SIMP ( statut = "o",
+                                                  typ = "TXM",
+                                                  into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                  fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                 ),
+
+  ), # Fin BLOC Parametres_TEMP_IMPO
+
+  Parametres_FLUX_REP = BLOC ( condition = " TypeConditionLimiteThermique in ( 'FLUX_REP', ) ",
+
+    Liste_flux_rep = SIMP ( statut = "o",
+                                     typ = "R",
+                                          max = "**",
+                                          fr = "Instant (s) / Flux de chaleur impose (W/m2) ",
+                                         ),
+
+     Prolongation_flux_rep = SIMP ( statut = "o",
+                                                 typ = "TXM",
+                                                 into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                 fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                ),
+
+  ), # Fin BLOC Parametres_FLUX_REP
+
+  Parametres_ECHANGE = BLOC ( condition = " TypeConditionLimiteThermique in ( 'ECHANGE', ) ",
+
+    Liste_echange_temp = SIMP ( statut = "o",
+                                              typ = "R",
+                                              max = "**",
+                                              fr = "Instant (s) / Temperature impose (°C) ",
+                                             ),
+
+    Prolongation_echange_temp = SIMP ( statut = "o",
+                                                     typ = "TXM",
+                                                     into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                     fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                    ),
+
+    Liste_echange_coef = SIMP ( statut = "o",
+                                              typ = "R",
+                                              max = "**",
+                                              fr = "Instant (s) / Coefficient d echange (W/m2/K) ",
+                                             ),
+
+    Prolongation_echange_coef = SIMP ( statut = "o",
+                                                     typ = "TXM",
+                                                     into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                     fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                    ),
+
+  ), # Fin BLOC Parametres_ECHANGE
+
+  Parametres_DEBIT = BLOC ( condition = " TypeConditionLimiteThermique in ( 'DEBIT', ) ",
+
+    Liste_debit = SIMP ( statut = "o",
+                                       typ = "R",
+                                       max = "**",
+                                       fr = "Instant (s) / Debit massique (kg/s) ",
+                                      ),
+
+    Prolongation_debit = SIMP ( statut = "o",
+                                              typ = "TXM",
+                                              into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                              fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                             ),
+
+    Liste_temp_injection = SIMP ( statut = "o",
+                                               typ = "R",
+                                                max = "**",
+                                                fr = "Instant (s) / Temperature d injection de securite  (°C) ",
+                                               ),
+
+    Prolongation_temp_injection = SIMP ( statut = "o",
+                                                       typ = "TXM",
+                                                       into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                       fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                      ),
+
+    # DH
+    DiametreHydraulique = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Diametre hydraulique (m) ",
+                             ),
+
+    # SECTION
+    SectionEspaceAnnulaire = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "-2.",
+                                   fr = "Section espace annulaire (m2) ",
+                                  ),
+
+    # DELTA
+    HauteurCaracConvectionNaturelle = SIMP ( statut = "o",
+                                 typ = "R",
+                                 defaut = "-2.",
+                                 fr = "Hauteur caracteristique convection naturelle (m) ",
+                                ),
+
+    # ALPHA_CF
+    CoefficientVestaleConvectionForcee = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection forcee (-) ",
+                                   ),
+
+    # ALPHA_CN
+    CoefficientVestaleConvectionNaturelle = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection naturelle (-) ",
+                                   ),
+
+    # EPS
+    CritereConvergenceRelative = SIMP ( statut = "o",
+                               typ = "R",
+                               defaut = "0.00001",
+                               fr = "Critere convergence relative (-) ",
+                              ),
+
+    # VM
+    VolumeMelangeCREARE = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Volume de melange CREARE (m3) ",
+                             ),
+
+    # T0
+    TemperatureInitialeCREARE = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Temperature initiale CREARE (degC) ",
+                             ),
+
+    # SE
+    SurfaceEchangeFluideStructure = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Surface d'echange fluide/structure (m2) ",
+                             ),
+
+  ), # Fin BLOC Parametres_DEBIT
+
+  Parametres_TEMP_FLU = BLOC ( condition = " TypeConditionLimiteThermique in ( 'TEMP_FLU', ) ",
+
+    Liste_temp_flu = SIMP ( statut = "o",
+                                          typ = "R",
+                                          max = "**",
+                                          fr = "Instant (s) / Debit massique (kg/s) ",
+                                         ),
+
+    Prolongation_temp_flu = SIMP ( statut = "o",
+                                                 typ = "TXM",
+                                                 into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                 fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                ),
+
+    Liste_debit_injection = SIMP ( statut = "o",
+                                                 typ = "R",
+                                                 max = "**",
+                                                 fr = "Instant (s) / Debit d injection de securite (kg/s) ",
+                                                ),
+
+    Prolongation_debit_injection = SIMP ( statut = "o",
+                                                        typ = "TXM",
+                                                        into = ( 'CC', 'CE', 'CL', 'EC', 'EE', 'EL', 'LC', 'LE', 'LL'),
+                                                        fr = "Prolongation aux frontieres amont et aval : C=constant, E=exclu, L=lineaire ",
+                                                       ),
+
+    # DH
+    DiametreHydraulique = SIMP ( statut = "o",
+                              typ = "R",
+                              defaut = "-2.",
+                              fr = "Diametre hydraulique (m) ",
+                             ),
+
+    # SECTION
+    SectionEspaceAnnulaire = SIMP ( statut = "o",
+                                   typ = "R",
+                                   defaut = "-2.",
+                                   fr = "Section espace annulaire (m2) ",
+                                  ),
+
+    # DELTA
+    HauteurCaracConvectionNaturelle = SIMP ( statut = "o",
+                                 typ = "R",
+                                 defaut = "-2.",
+                                 fr = "Hauteur caracteristique convection naturelle (m) ",
+                                ),
+
+    # ALPHA_CF
+    CoefficientVestaleConvectionForcee= SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection forcee (-) ",
+                                   ),
+
+    # ALPHA_CN
+    CoefficientVestaleConvectionNaturelle = SIMP ( statut = "o",
+                                    typ = "R",
+                                    defaut = "1.",
+                                    fr = "Coefficient Vestale convection naturelle (-) ",
+                                   ),
+
+    # EPS
+    CritereConvergenceRelative = SIMP ( statut = "o",
+                               typ = "R",
+                               defaut = "0.00001",
+                               fr = "Critere convergence relative (-) ",
+                              ),
+
+  ), # Fin BLOC Parametres_TEMP_FLU
+
+) # Fin OPER TRANSITOIRE
+
+
+#================================
+# 3. Definition de l'etude
+#================================
+
+# Nota : les variables de type PROC doivent etre en majuscules !
+CRITERIA = PROC ( nom = "CRITERIA",
+                  op = None,
+                  docu = "",
+                  fr = "Mise en donnee pour le fichier de configuration de OPENTURNS.",
+                  ang = "Writes the configuration file for OPENTURNS.",
+
+
+
+  Type = SIMP ( statut = "o",
+                typ = "TXM",
+                into = ( "Min/Max", "Central Uncertainty", "Threshold Exceedence" ),
+                fr = "Type d'Analyse",
+                ang = "Analysis",
+                ),
+
+
+
+
+  MinMax = BLOC ( condition = " Type in ( 'Min/Max', ) ",
+
+                  Method = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( "Experiment Plane", "Random Sampling" ),
+                                  fr = "Methode",
+                                  ang = "Method",
+                                  ),
+                  # UC 3.1.1
+                  ExperimentPlaneSettings = BLOC ( condition = " Method in ( 'Experiment Plane', ) ",
+
+                          ExperimentPlane = SIMP ( statut = "o",
+                                                   typ = "TXM",
+                                                   into = ( "Axial", "Factorial", "Composite", ),
+                                                   fr = "Methode",
+                                                   ang = "Method",
+                                                   ),
+
+                          Levels = SIMP ( statut = "o",
+                                          typ = "R",
+                                          val_min = 0.0,
+                                          max = '**',    
+                                          fr = "Nombre de niveaux dans chaque direction",
+                                          ang = "Levels in each direction",
+                                          ),
+
+                          # Scaled Vector
+                          UnitsPerDimension = SIMP ( statut = "o",
+                                          typ = "R",
+                                          max = '**',    
+                                          fr = "Unite par dimension (autant que de variables declarees)",
+                                          ang = "Units per dimension (as much as declared variables)",
+                                          ),
+
+                          # Translation Vector
+                          Center = SIMP ( statut = "o",
+                                          typ = "R",
+                                          max = '**',    
+                                          fr = "Unite par dimension",
+                                          ang = "Units per dimension",
+                                          ),
+
+                    ), # Fin BLOC ExperimentPlaneSettings
+
+
+
+                  RandomSamplingSettings = BLOC ( condition = " Method in ( 'Random Sampling', ) ",
+
+                          PointsNumber = SIMP ( statut = "o",
+                                                typ = "I",
+                                                val_min = 1,
+                                                fr = "Nombre de points",
+                                                ang = "Points number",
+                                                ),
+
+                    ), # Fin BLOC RandomSamplingSettings
+
+                  Result = SIMP (  statut = "o",
+                                   typ = "TXM",
+                                   into = ( "Min/Max", ),
+                                   defaut = "Min/Max",
+                                   fr = "Le minimum et le maximum",
+                                   ang = "The min and max values",
+                                   ),
+
+
+  ), # Fin BLOC MinMax
+
+
+
+
+  CentralUncertainty = BLOC ( condition = " Type in ( 'Central Uncertainty', ) ",
+
+                  Method = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( "Taylor Variance Decomposition", "Random Sampling" ),
+                                  fr = "Methode",
+                                  ang = "Method",
+                                  ),
+                              
+                  # UC 3.2.
+                  TaylorVarianceDecompositionSettings = BLOC ( condition = " Method in ( 'Taylor Variance Decomposition', ) ",
+
+                      Result = FACT ( statut = "o",
+                                      min = 1,
+                                      max = "**",
+
+                              MeanFirstOrder = SIMP ( statut = "o",
+                                                typ = 'TXM',
+                                                into = ( 'yes', 'no' ),
+                                                defaut = 'yes',
+                                                max = 1,
+                                                fr = "Moyenne au premier ordre",
+                                                ang = "MeanFirstOrder",
+                                                ),
+
+                              StandardDeviationFirstOrder = SIMP ( statut = "o",
+                                                                   typ = 'TXM',
+                                                                   into = ( 'yes', 'no' ),
+                                                                   defaut = 'yes',
+                                                                   max = 1,
+                                                                   fr = "Ecart-type au premier ordre",
+                                                                   ang = "StandardDeviationFirstOrder",
+                                                                   ),
+
+                              MeanSecondOrder = SIMP ( statut = "o",
+                                                       typ = 'TXM',
+                                                       into = ( 'yes', 'no' ),
+                                                       defaut = 'yes',
+                                                       max = 1,
+                                                       fr = "Moyenne au second ordre",
+                                                       ang = "MeanSecondOrder",
+                                                       ),
+
+                              ImportanceFactor = SIMP ( statut = "o",
+                                                        typ = 'TXM',
+                                                        into = ( 'yes', 'no' ),
+                                                        defaut = 'no',
+                                                        max = 1,
+                                                        fr = "Facteur d'importance pour variable de sortie scalaire",
+                                                        ang = "ImportanceFactor",
+                                                        ),
+
+                             ImportanceFactorSettings = BLOC ( condition = " ImportanceFactor in ( 'yes', ) ",
+
+                                    NumericalResults  = SIMP ( statut = "o",
+                                                               typ = 'TXM',
+                                                               into = ( 'yes', 'no' ),
+                                                               defaut = 'yes',
+                                                               max = 1,
+                                                               fr = "Resultats numeriques",
+                                                               ang = "NumericalResults",
+                                                               ),
+
+                                     GraphicalResults  = SIMP ( statut = "o",
+                                                                typ = 'TXM',
+                                                                into = ( 'yes', 'no' ),
+                                                                defaut = 'no',
+                                                                max = 1,
+                                                                fr = "Resultats graphiques",
+                                                                ang = "GraphicalResults",
+                                                                ),
+
+                            ), # Fin BLOC ImportanceFactorSettings
+
+                      ), # Fin FACT Result
+                                                               
+                  ), # Fin BLOC TaylorVarianceDecompositionSettings
+
+
+
+                  RandomSamplingSettings = BLOC ( condition = " Method in ( 'Random Sampling', ) ",
+
+                          PointsNumber = SIMP ( statut = "o",
+                                                typ = "I",
+                                                val_min = 1,
+                                                fr = "Nombre de points",
+                                                ang = "Points number",
+                                                ),
+
+                       Result = FACT ( statut = "o",
+                                      min = 1,
+                                      max = "**",
+
+                              EmpiricalMean = SIMP ( statut = "o",
+                                                     typ = 'TXM',
+                                                     into = ( 'yes', 'no' ),
+                                                     defaut = 'yes',
+                                                     max = 1,
+                                                     fr = "Moyenne empirique",
+                                                     ang = "Empirical mean",
+                                                     ),
+
+                              EmpiricalStandardDeviation = SIMP ( statut = "o",
+                                                                  typ = 'TXM',
+                                                                  into = ( 'yes', 'no' ),
+                                                                  defaut = 'yes',
+                                                                  max = 1,
+                                                                  fr = "Ecart-type empirique",
+                                                                  ang = "Empirical standard deviation",
+                                                                  ),
+
+                              EmpiricalQuantile = SIMP ( statut = "o",
+                                                         typ = 'R',
+                                                         defaut = 0.0,
+                                                         max = 1,
+                                                         val_min = 0.0,
+                                                         val_max = 1.0,
+                                                         fr = "Quantile empirique",
+                                                         ang = "Empirical quantile",
+                                                         ),
+
+                              AnalysedCorrelations = SIMP ( statut = "o",
+                                                            typ = 'TXM',
+                                                            into = ( 'yes', 'no' ),
+                                                            defaut = 'no',
+                                                            max = 1,
+                                                            fr = "Correlations analysees",
+                                                            ang = "Analysed correlations",
+                                                            ),
+
+                              KernelSmoothing = SIMP ( statut = "o",
+                                                       typ = 'TXM',
+                                                       into = ( 'yes', 'no' ),
+                                                       defaut = 'no',
+                                                       max = 1,
+                                                       fr = "Kernel smoothing de l'echantillon",
+                                                       ang = "Kernel smoothing of the sample",
+                                                       ),
+
+                      ), # Fin FACT Result
+                                                               
+                  ), # Fin BLOC RandomSamplingSettings
+
+  ), # Fin BLOC CentralUncertainty
+
+
+
+
+  ThresholdExceedence = BLOC ( condition = " Type in ( 'Threshold Exceedence', ) ",
+
+         Event =  FACT ( statut = "o",
+                         min = 1,
+                         max = 1,
+
+                         Threshold = SIMP ( statut = "o",
+                                            typ = "R",
+                                            max = 1,
+                                            fr = "Le seuil de defaillance",
+                                            ang = "Failure threshold",
+                                            ),
+
+                         ComparisonOperator = SIMP ( statut = "o",
+                                                     typ = "TXM",
+                                                     max = 1,
+                                                     into = ( "Less", "LessOrEqual", "Equal", "GreaterOrEqual", "Greater" ),
+                                                     fr = "Que faut-il ne pas depasser : un maximum ou un minimum",
+                                                     ang = "What is the failure threshold : maximum or minimum",
+                                                     ),
+         ), # Fin FACT Event
+                         
+
+         Method = SIMP ( statut = "o",
+                         typ = "TXM",
+                         into = ( "Simulation", "Analytical" ),
+                         fr = "Methode",
+                         ang = "Method",
+                         ),
+
+         SimulationSettings = BLOC ( condition = " Method in ( 'Simulation', ) ",
+
+               Algorithm = SIMP ( statut = "o",
+                                  typ = "TXM",
+                                  into = ( "MonteCarlo", "LHS", "ImportanceSampling" ),
+                                  fr = "Algorithme de simulation",
+                                  ang = "Simulation algorithm",
+                                  ),
+
+                                 
+               RandomGenerator = FACT ( statut = "o",
+                                        min = 1,
+                                        max = 1,
+
+                           SeedToBeSet = SIMP ( statut = "o",
+                                                typ = 'TXM',
+                                                into = ( 'yes', 'no' ),
+                                                defaut = 'no',
+                                                max = 1,
+                                                fr = "La racine du generateur aleatoire doit-elle etre positionnee ?",
+                                                ang = "Does the random generator seed need to be set ?",
+                                                ),
+
+                           SeedSettings = BLOC ( condition = " SeedToBeSet in ( 'yes', ) ",
+
+                                                 RandomGeneratorSeed = SIMP ( statut = "o",
+                                                                              typ = "I",
+                                                                              max = 1,
+                                                                              fr = "Racine du generateur aleatoire",
+                                                                              ang = "Random generator seed",
+                                                                              ),
+
+                                               ), # Fin BLOC SeedSettings
+
+               ), # Fin FACT RandomGenerator
+
+
+               BlockSize = SIMP ( statut = "f",
+                                  typ = "I",
+                                  max = 1,
+                                  val_min = 1,
+                                  defaut = 1,
+                                  fr = "Nombre de calculs realises en bloc",
+                                  ang = "Number of computations as a block",
+                                  ),
+
+               MaximumOuterSampling = SIMP ( statut = "o",
+                                             typ = "I",
+                                             max = 1,
+                                             val_min = 1,
+                                             fr = "Maximum d'iterations externes",
+                                             ang = "Maximum outer Sampling value",
+                                             ),
+
+               MaximumCoefficientOfVariation = SIMP ( statut = "f",
+                                                      typ = "R",
+                                                      max = 1,
+                                                      defaut = 0.1,
+                                                      val_min = 0.0,
+                                                      fr = " maximum ...",
+                                                      ang = "Absolute maximum ...."
+                                                      ),
+
+               ImportanceSamplingSettings = BLOC ( condition = " Algorithm in ( 'ImportanceSampling', ) ",
+
+                            MeanVector = SIMP ( statut = "o",
+                                                typ = "R",
+                                                max = "**",
+                                                fr = "Moyenne",
+                                                ang = "Mean vector",
+                                                ),
+
+                            Correlation = SIMP ( statut = "o",
+                                                 typ = 'TXM',
+                                                 into = ( 'Independent', 'Linear' ),
+                                                 defaut = 'Linear',
+                                                 max = 1,
+                                                 fr = "Le type de correlation entre les variables",
+                                                 ang = "Correlation between variables",
+                                                 ),
+
+               ), # Fin BLOC ImportanceSamplingSettings
+
+               Result = FACT ( statut = "o",
+                                      min = 1,
+                                      max = "**",
+
+                    Probability = SIMP ( statut = "o",
+                                         typ = 'TXM',
+                                         into = ( 'yes', ),
+                                         defaut = 'yes',
+                                         max = 1,
+                                         fr = "Probabiblite",
+                                         ang = "Probability",
+                                         ),
+
+                    ConfidenceInterval = SIMP ( statut = "o",
+                                                typ = 'TXM',
+                                                into = ( 'yes', 'no' ),
+                                                defaut = 'yes',
+                                                max = 1,
+                                                fr = "Ecart-type empirique",
+                                                ang = "Empirical standard deviation",
+                                                ),
+
+                    ConfidenceIntervalSettings = BLOC ( condition = " ConfidenceInterval in ( 'yes', ) ",
+
+                          Level = SIMP ( statut = "o",
+                                         typ = 'R',
+                                         defaut = 0.0,
+                                         max = 1,
+                                         val_min = 0.0,
+                                         val_max = 1.0,
+                                         fr = "Niveau de confiance",
+                                         ang = "Confidence level",
+                                         ),
+                                                     
+                    ), # Fin BLOC ConfidenceIntervalSettings
+                               
+                    VariationCoefficient = SIMP ( statut = "o",
+                                                  typ = 'TXM',
+                                                  into = ( 'yes', 'no' ),
+                                                  defaut = 'yes',
+                                                  max = 1,
+                                                  fr = "Coefficient de variation",
+                                                  ang = "VariationCoefficient",
+                                                  ),
+
+                    IterationNumber = SIMP ( statut = "o",
+                                             typ = 'TXM',
+                                             into = ( 'yes', 'no' ),
+                                             defaut = 'yes',
+                                             max = 1,
+                                             fr = "Nombre d'iterations",
+                                             ang = "Iteration number",
+                                             ),
+
+                    ConvergenceGraph = SIMP ( statut = "o",
+                                             typ = 'TXM',
+                                             into = ( 'yes', 'no' ),
+                                             defaut = 'yes',
+                                             max = 1,
+                                             fr = "Graphe de convergence",
+                                             ang = "Convergence graph",
+                                             ),
+
+               ), # Fin FACT Result
+                                                               
+
+
+         ), # Fin BLOC SimulationSettings
+
+
+                               
+         AnalyticalSettings = BLOC ( condition = " Method in ( 'Analytical', ) ",
+
+                Approximation = SIMP ( statut = "o",
+                                       typ = "TXM",
+                                       into = ( "FORM", "SORM" ),
+                                       fr = "Approximation",
+                                       ang = "Approximation",
+                                       ),
+
+                OptimizationAlgorithm = SIMP ( statut = "o",
+                                               typ = "TXM",
+                                               into = ( "Cobyla", "AbdoRackwitz" ),
+                                               fr = "Methode d'optimisation",
+                                               ang = "Optimisation method",
+                                               ),
+
+                                     
+                PhysicalStartingPoint = SIMP ( statut = "f",
+                                               typ = "R",
+                                               max = "**",
+                                               fr = "Point de demarrage de l'algorithme iteratif",
+                                               ang = "Initial point for iterative process",
+                                               ),
+
+                MaximumIterationsNumber = SIMP ( statut = "f",
+                                                 typ = "I",
+                                                 max = 1,
+                                                 val_min = 1,
+                                                 fr = "Nombre maximum d iterations",
+                                                 ang = "Maximum number of iterations",
+                                                 ),
+
+                regles = ( EXCLUS ( "MaximumAbsoluteError", "RelativeAbsoluteError" ),  ),
+                                     
+                MaximumAbsoluteError = SIMP ( statut = "f",
+                                              typ = "R",
+                                              max = 1,
+                                              defaut = 1E-6,
+                                              val_min = 0.0,
+                                              fr = "Distance maximum absolue entre 2 iterations successifs",
+                                              ang = "Absolute maximum distance between 2 successive iterates",
+                                              ),
+
+                RelativeAbsoluteError = SIMP ( statut = "f",
+                                               typ = "R",
+                                               max = 1,
+                                               defaut = 1E-6,
+                                               val_min = 0.0,
+                                               fr = "Distance maximum relative entre 2 iterations successives",
+                                               ang = "Relative maximum distance between 2 successive iterates",
+                                               ),
+                                     
+                MaximumConstraintError = SIMP ( statut = "f",
+                                                typ = "R",
+                                                max = 1,
+                                                val_min = 0.0,
+                                                fr = "Valeur maximum absolue de la fonction moins la valeur du niveau",
+                                                ang = "Maximum absolute value of the constraint function minus the level value",
+                                                ),
+
+                ImportanceSampling = SIMP ( statut = "o",
+                                            typ = 'TXM',
+                                            into = ( 'yes', 'no' ),
+                                            defaut = 'no',
+                                            max = 1,
+                                            fr = "Tirage d'importance au point de conception",
+                                            ang = "Importance sampling at design point",
+                                            ),
+
+                FORM = BLOC ( condition = " Approximation in ( 'FORM', ) ",
+
+                    Probability = SIMP ( statut = "o",
+                                         typ = 'TXM',
+                                         into = ( 'yes', ),
+                                         defaut = 'yes',
+                                         max = 1,
+                                         fr = "Probabiblite",
+                                         ang = "Probability",
+                                         ),
+
+                    DesignPoint = SIMP ( statut = "o",
+                                         typ = 'TXM',
+                                         into = ( 'yes', 'no' ),
+                                         defaut = 'yes',
+                                         max = 1,
+                                         fr = "Point de conception",
+                                         ang = "Design point",
+                                         ),
+
+                    HasReliabilityIndex = SIMP ( statut = "o",
+                                                 typ = 'TXM',
+                                                 into = ( 'yes', 'no' ),
+                                                 defaut = 'yes',
+                                                 max = 1,
+                                                 fr = "Indice de fiabilite",
+                                                 ang = "Reliability index",
+                                                 ),
+
+                    ImportanceFactor = SIMP ( statut = "o",
+                                              typ = 'TXM',
+                                              into = ( 'yes', 'no' ),
+                                              defaut = 'no',
+                                              max = 1,
+                                              fr = "Facteur d'importance pour variable de sortie scalaire",
+                                              ang = "ImportanceFactor",
+                                              ),
+
+                    ImportanceFactorSettings = BLOC ( condition = " ImportanceFactor in ( 'yes', ) ",
+
+                            NumericalResults  = SIMP ( statut = "o",
+                                                       typ = 'TXM',
+                                                       into = ( 'yes', 'no' ),
+                                                       defaut = 'yes',
+                                                       max = 1,
+                                                       fr = "Resultats numeriques",
+                                                       ang = "NumericalResults",
+                                                       ),
+
+                             GraphicalResults  = SIMP ( statut = "o",
+                                                        typ = 'TXM',
+                                                        into = ( 'yes', 'no' ),
+                                                        defaut = 'no',
+                                                        max = 1,
+                                                        fr = "Resultats graphiques",
+                                                        ang = "GraphicalResults",
+                                                        ),
+
+                    ), # Fin BLOC ImportanceFactorSettings
+
+
+                    SensitivityAnalysis = SIMP ( statut = "o",
+                                                 typ = 'TXM',
+                                                 into = ( 'yes', 'no' ),
+                                                 defaut = 'no',
+                                                 max = 1,
+                                                 fr = "Analyse de sensibilite",
+                                                 ang = "Sensitivity analysis",
+                                                 ),
+
+                    SensitivityAnalysisSettings = BLOC ( condition = " SensitivityAnalysis in ( 'yes', ) ",
+
+                            HasoferReliabilityIndex = SIMP ( statut = "o",
+                                                             typ = 'TXM',
+                                                             into = ( 'yes', 'no' ),
+                                                             defaut = 'no',
+                                                             max = 1,
+                                                             fr = "Indice de fiabilite de Hasofer",
+                                                             ang = "Hasofer reliability index",
+                                                             ),
+        
+                            HasoferReliabilityIndexSettings = BLOC ( condition = " HasoferReliabilityIndex in ( 'yes', ) ",
+        
+                                    NumericalResults  = SIMP ( statut = "o",
+                                                               typ = 'TXM',
+                                                               into = ( 'yes', 'no' ),
+                                                               defaut = 'yes',
+                                                               max = 1,
+                                                               fr = "Resultats numeriques",
+                                                               ang = "NumericalResults",
+                                                               ),
+        
+                                     GraphicalResults  = SIMP ( statut = "o",
+                                                                typ = 'TXM',
+                                                                into = ( 'yes', 'no' ),
+                                                                defaut = 'no',
+                                                                max = 1,
+                                                                fr = "Resultats graphiques",
+                                                                ang = "GraphicalResults",
+                                                                ),
+
+                            ), # Fin BLOC HasoferReliabilityIndexSettings
+                                                         
+                    ), # Fin BLOC SensitivityAnalysisSettings
+
+                    FunctionCallsNumber = SIMP ( statut = "o",
+                                                 typ = 'TXM',
+                                                 into = ( 'yes', 'no' ),
+                                                 defaut = 'no',
+                                                 max = 1,
+                                                 fr = "Nombre d'appels a la fonction",
+                                                 ang = "Function calls number",
+                                                 ),
+
+
+                ), # Fin BLOC FORM
+
+
+                SORM = BLOC ( condition = " Approximation in ( 'SORM', ) ",
+
+
+                    TvedtApproximation = SIMP ( statut = "o",
+                                                typ = 'TXM',
+                                                into = ( 'yes', 'no' ),
+                                                defaut = 'no',
+                                                max = 1,
+                                                fr = "Approximation de Tvedt",
+                                                ang = "Tvedt approximation",
+                                                ),
+
+                    HohenBichlerApproximation = SIMP ( statut = "o",
+                                                       typ = 'TXM',
+                                                       into = ( 'yes', 'no' ),
+                                                       defaut = 'no',
+                                                       max = 1,
+                                                       fr = "Approximation de HohenBichler",
+                                                       ang = "HohenBichler approximation",
+                                                       ),
+
+                    BreitungApproximation = SIMP ( statut = "o",
+                                                   typ = 'TXM',
+                                                   into = ( 'yes', 'no' ),
+                                                   defaut = 'no',
+                                                   max = 1,
+                                                   fr = "Approximation de Breitung",
+                                                   ang = "Breitung approximation",
+                                                   ),
+
+                    DesignPoint = SIMP ( statut = "o",
+                                         typ = 'TXM',
+                                         into = ( 'yes', 'no' ),
+                                         defaut = 'yes',
+                                         max = 1,
+                                         fr = "Point de conception",
+                                         ang = "Design point",
+                                         ),
+
+                    ImportanceFactor = SIMP ( statut = "o",
+                                              typ = 'TXM',
+                                              into = ( 'yes', 'no' ),
+                                              defaut = 'no',
+                                              max = 1,
+                                              fr = "Facteur d'importance pour variable de sortie scalaire",
+                                              ang = "ImportanceFactor",
+                                              ),
+
+                    ImportanceFactorSettings = BLOC ( condition = " ImportanceFactor in ( 'yes', ) ",
+
+                            NumericalResults  = SIMP ( statut = "o",
+                                                       typ = 'TXM',
+                                                       into = ( 'yes', 'no' ),
+                                                       defaut = 'yes',
+                                                       max = 1,
+                                                       fr = "Resultats numeriques",
+                                                       ang = "NumericalResults",
+                                                       ),
+
+                             GraphicalResults  = SIMP ( statut = "o",
+                                                        typ = 'TXM',
+                                                        into = ( 'yes', 'no' ),
+                                                        defaut = 'no',
+                                                        max = 1,
+                                                        fr = "Resultats graphiques",
+                                                        ang = "GraphicalResults",
+                                                        ),
+
+                    ), # Fin BLOC ImportanceFactorSettings
+
+
+                    SensitivityAnalysis = SIMP ( statut = "o",
+                                                 typ = 'TXM',
+                                                 into = ( 'yes', 'no' ),
+                                                 defaut = 'no',
+                                                 max = 1,
+                                                 fr = "Analyse de sensibilite",
+                                                 ang = "Sensitivity analysis",
+                                                 ),
+
+                    SensitivityAnalysisSettings = BLOC ( condition = " SensitivityAnalysis in ( 'yes', ) ",
+
+                            HasoferReliabilityIndex = SIMP ( statut = "o",
+                                                             typ = 'TXM',
+                                                             into = ( 'yes', 'no' ),
+                                                             defaut = 'no',
+                                                             max = 1,
+                                                             fr = "Indice de fiabilite de Hasofer",
+                                                             ang = "Hasofer reliability index",
+                                                             ),
+        
+                            HasoferReliabilityIndexSettings = BLOC ( condition = " HasoferReliabilityIndex in ( 'yes', ) ",
+        
+                                    NumericalResults  = SIMP ( statut = "o",
+                                                               typ = 'TXM',
+                                                               into = ( 'yes', 'no' ),
+                                                               defaut = 'yes',
+                                                               max = 1,
+                                                               fr = "Resultats numeriques",
+                                                               ang = "NumericalResults",
+                                                               ),
+        
+                                     GraphicalResults  = SIMP ( statut = "o",
+                                                                typ = 'TXM',
+                                                                into = ( 'yes', 'no' ),
+                                                                defaut = 'no',
+                                                                max = 1,
+                                                                fr = "Resultats graphiques",
+                                                                ang = "GraphicalResults",
+                                                                ),
+
+                            ), # Fin BLOC HasoferReliabilityIndexSettings
+                                                         
+                    ), # Fin BLOC SensitivityAnalysisSettings
+
+                    FunctionCallsNumber = SIMP ( statut = "o",
+                                                 typ = 'TXM',
+                                                 into = ( 'yes', 'no' ),
+                                                 defaut = 'no',
+                                                 max = 1,
+                                                 fr = "Nombre d'appels a la fonction",
+                                                 ang = "Function calls number",
+                                                 ),
+
+
+                ), # Fin BLOC SORM
+
+
+                                     
+        ), # Fin BLOC AnalyticalSettings
+
+
+                               
+  ), # Fin BLOC ThresholdExceedence
+
+
+
+) # Fin PROC CRITERIA
+
+
+#===============================
+# 5. Definition des parametres
+#===============================
+VARI = OPER ( nom = "VARI",
+                      sd_prod = variable,
+                      op = None,
+                      fr = "Definitions des lois marginales utilisees par les variables d'entree", 
+                      type=SIMP(statut='f',defaut="IN",into=("IN","OUT"), typ = "TXM",)
+              )
+
+ESSAI=PROC(nom="ESSAI",
+       op=None,
+       fr="Essai",
+       ang = "Test",
+       
+       MALOI       = SIMP(statut='o',typ=(loi,),),
+       MAVARIABLE  = SIMP(statut='o',typ=(variable,),),
+) ;
+
+                     
+
+
index ac75b5ba1e043f6ec1a588b2867cedc153eec647..f66095a4db38eb0f341f98eb1731925b7e9cf070 100644 (file)
@@ -5,7 +5,7 @@ rep_cata=os.getcwd()
 
 catalogues = (
 # (code,version,catalogue,formatIn,formatOut)
-  ('CUVE2DG','V1',os.path.join(rep_cata,'Cuve2dg_Cata_V1.py'),'python','cuve2dg'),
-  #('CUVE2DG','V2',os.path.join(rep_cata,'Cuve2dg_Cata_V2.py'),'python','cuve2dg'),
+  #('CUVE2DG','V1',os.path.join(rep_cata,'Cuve2dg_Cata_V1.py'),'python','cuve2dg'),
+  ('CUVE2DG','V2',os.path.join(rep_cata,'Cuve2dg_Cata_V2.py'),'python','cuve2dg'),
 )