From: Pascale Noyret Date: Tue, 27 Jan 2009 09:29:57 +0000 (+0000) Subject: *** empty log message *** X-Git-Tag: V0_a1_SEPQT4~34 X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=4dfb22f16a318b08ff8e7b454fe569f28d8ce2a9;p=tools%2Feficas.git *** empty log message *** --- diff --git a/Cuve2dg/Cuve2dg_Cata_V1.py b/Cuve2dg/Cuve2dg_Cata_V1.py index 61d9b721..47e23641 100644 --- a/Cuve2dg/Cuve2dg_Cata_V1.py +++ b/Cuve2dg/Cuve2dg_Cata_V1.py @@ -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 index 00000000..b9c202df --- /dev/null +++ b/Cuve2dg/Cuve2dg_Cata_V2.py @@ -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,),), +) ; + + + + diff --git a/Cuve2dg/catalogues_cuve2dg.ini b/Cuve2dg/catalogues_cuve2dg.ini index ac75b5ba..f66095a4 100644 --- a/Cuve2dg/catalogues_cuve2dg.ini +++ b/Cuve2dg/catalogues_cuve2dg.ini @@ -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'), )