From 3e19deb985129014027a5f126b276f1d21c0ecc7 Mon Sep 17 00:00:00 2001 From: Pascale Noyret Date: Tue, 26 Jan 2010 10:02:42 +0000 Subject: [PATCH] *** empty log message *** --- Cuve2dg/Cuve2dg_Cata_V1.py | 2264 +++++++++++++++++++++++++ Cuve2dg/Cuve2dg_Cata_V2.py | 1851 ++++++++++++++++++++ Cuve2dg/Cuve2dg_Cata_V3.py | 1050 ++++++++++++ Cuve2dg/Cuve2dg_Cata_VPN.py | 111 ++ Cuve2dg/generator_cuve2dg.py | 1537 +++++++++++++++++ Cuve2dg/sdistCuveqt.py | 102 ++ InterfaceQT4/monPlusieursBasePanel.py | 2 + InterfaceQT4/monUniqueBasePanel.py | 1 + 8 files changed, 6918 insertions(+) create mode 100644 Cuve2dg/Cuve2dg_Cata_V1.py create mode 100644 Cuve2dg/Cuve2dg_Cata_V2.py create mode 100644 Cuve2dg/Cuve2dg_Cata_V3.py create mode 100644 Cuve2dg/Cuve2dg_Cata_VPN.py create mode 100644 Cuve2dg/generator_cuve2dg.py create mode 100644 Cuve2dg/sdistCuveqt.py diff --git a/Cuve2dg/Cuve2dg_Cata_V1.py b/Cuve2dg/Cuve2dg_Cata_V1.py new file mode 100644 index 00000000..47e23641 --- /dev/null +++ b/Cuve2dg/Cuve2dg_Cata_V1.py @@ -0,0 +1,2264 @@ +# -*- 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 = 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 du defaut +#==== + + Parametres_DSR = BLOC ( condition = " TYPEDEF in ( 'DSR', ) ", + + 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)", + ), + + 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', ) ", + + LONGSURPROF = 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 = SIMP ( statut = "o", + typ = "R", + #defaut = "-0.00001", + fr = "Decalage radial du defaut sous revetement (m)", + ), + + ANGLDEF = SIMP ( statut = "o", + typ = "R", + defaut = "0.", + fr = "Coordonnee angulaire du defaut (degres)", + ), + + ALTIDEF = SIMP ( statut = "o", + typ = "R", + defaut = "2.", + fr = "Altitude du defaut (m)", + ), + + POINDEF = SIMP ( statut = "o", + typ = "TXM", + into = ( "A", "B" ), + defaut = "A", + fr = "Choix du point considere du defaut sous revetement", + ), + + 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)", + ), + + CORRECPLASTIC = SIMP ( statut = "o", + typ = "TXM", + into = ( "OUI", "NON" ), + defaut = "NON", + fr = "Prise en compte de la correction plastique BETA ", + ), + + ), # Fin BLOC Parametres_DSR + + 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)", + ), + + ANGLDEF = SIMP ( statut = "o", + typ = "R", + defaut = "0.", + fr = "Coordonnee angulaire du defaut (degres)", + ), + + ALTIDEF = SIMP ( statut = "o", + typ = "R", + defaut = "2.", + fr = "Altitude du defaut (m)", + ), + + 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)", + ), + + 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 ", + ), + + ), # Fin BLOC Parametres_DD + +) # 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 = SIMP ( statut = "o", + typ = "TXM", + into = ( "GEOMETRIE", "MAILLAGE"), + #defaut = "GEOMETRIE", + fr = "Traitement de la geometrie d'une cuve", + ), + + +#==== +# Definition des parametres selon le type de traitement de la geometrie +#==== + + 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)", + ), + + LREV = SIMP ( statut = "o", + typ = "R", + defaut = "0.0075", + fr = "Epaisseur du revetement (m)", + ), + + LIGMIN = 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 = " TYPEGEOM in ( 'MAILLAGE', ) ", + + 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 = 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", + ), + + +#==== +# Definition des parametres selon le modele de fluence +#==== + + Reglementaire = BLOC ( condition = " MODELFLUENCE in ( 'Reglementaire', ) ", + + 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 Reglementaire + + France = BLOC ( condition = " MODELFLUENCE in ( 'France', ) ", + + 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", + ), + + KPFRANCE = SIMP ( statut = "o", + typ = "R", + defaut = "12.7", + fr = "Parametre exponentiel du modele France", + ), + + ), # Fin BLOC France + + ValeurImposee = BLOC ( condition = " MODELFLUENCE in ( 'ValeurImposee', ) ", + + 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', ) ", + + 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 + + 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)", + ), + + 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)", + ), + + COEFFLUENCE5 = SIMP ( statut = "o", + typ = "R", + defaut = "3.37", + fr = "Fluence a l'azimut 20 (10^19 n/cm)", + ), + + COEFFLUENCE6 = SIMP ( statut = "o", + typ = "R", + defaut = "3.16", + fr = "Fluence a l'azimut 25 (10^19 n/cm)", + ), + + COEFFLUENCE7 = SIMP ( statut = "o", + typ = "R", + defaut = "2.74", + fr = "Fluence a l'azimut 30 (10^19 n/cm)", + ), + + COEFFLUENCE8 = SIMP ( statut = "o", + typ = "R", + defaut = "2.25", + fr = "Fluence a l'azimut 35 (10^19 n/cm)", + ), + + COEFFLUENCE9 = SIMP ( statut = "o", + typ = "R", + defaut = "1.89", + fr = "Fluence a l'azimut 40 (10^19 n/cm)", + ), + + COEFFLUENCE10 = 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 = SIMP ( statut = "o", + typ = "TXM", + into = ( "RTNDT", "FLUENCE" ), + #defaut = "RTNDT", + fr = "Type d'irradiation ", + ), + +#==== +# Definition des parametres selon le type d'irradiation +#==== + + Parametres_RTNDT = BLOC ( condition = " TYPIRR in ( 'RTNDT', ) ", + + RTNDT = SIMP ( statut = "o", + typ = "R", + defaut = "73.", + fr = "RTNDT finale (°C)", + ), + + ), # Fin BLOC Parametres_RTNDT + + 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 (%)", + ), + + 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)", + ), + + ), # Fin BLOC Parametres_USNRC + + nbectDRTNDT = SIMP ( statut = "o", + typ = "R", + defaut = "2.", + fr = "Nombre d ecart-type par rapport a la moyenne de DRTNDT", + ), + + ), # Fin BLOC Parametres_FLUENCE + +#======================== +# 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é +#==== + + Parametres_RCCM = BLOC ( condition = " MODELKIC in ( 'RCC-M', ) ", + + 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 (%) ", + ), + + 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 ", + ), + + nbectKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma) ", + ), + + + KIAPAL = 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 = 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', ) ", + + 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_RCCMpal + + Parametres_RCCMexp = BLOC ( condition = " MODELKIC in ( 'RCC-M_exp', ) ", + + 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_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 ", + ), + + 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_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", + 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_Wallin + + Parametres_REME = BLOC ( condition = " MODELKIC in ( 'REME', ) ", + + 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_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 = "**", + 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=exclus, 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=exclus, 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=exclus, 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=exclus, 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=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) ", + ), + + NUREV = 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 = 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) ", + ), + + 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 + + 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=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) ", + ), + + 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 ", + ), + + ), # Fin BLOC Parametres_TEMP_IMPO + + Parametres_FLUX_REP = BLOC ( condition = " TYPCLTH 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=exclus, L=lineaire ", + ), + + ), # Fin BLOC Parametres_FLUX_REP + + Parametres_ECHANGE = BLOC ( condition = " TYPCLTH 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=exclus, 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=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 (-) ", + ), + + 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 (-) ", + ), + + 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 + + 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) ", + ), + + 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 + + +#================================ +# 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/Cuve2dg_Cata_V2.py b/Cuve2dg/Cuve2dg_Cata_V2.py new file mode 100644 index 00000000..ab400041 --- /dev/null +++ b/Cuve2dg/Cuve2dg_Cata_V2.py @@ -0,0 +1,1851 @@ +# -*- coding: utf-8 -*- + +# -------------------------------------------------- +# Definition de variables sous forme de tuple +# -------------------------------------------------- + +import types +class Tuple: + def __init__(self,ntuple): + self.ntuple=ntuple + + def __convert__(self,valeur): + if type(valeur) == types.StringType: + return None + if len(valeur) != self.ntuple: + return None + return valeur + + def info(self): + return "Tuple de %s elements" % self.ntuple + + __repr__=info + __str__=info + +# -------------------------------------------------- +# 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'), + A_CLASSER ( ('OPTIONS'), ('DEFAUT') ), + A_CLASSER ( ('DEFAUT'), ('CUVE') ), + A_CLASSER ( ('CUVE'), ('MODELES') ), + A_CLASSER ( ('MODELES'), ('INITIALISATION') ), + A_CLASSER ( ('INITIALISATION'), ('REVETEMENT') ), + A_CLASSER ( ('REVETEMENT'), ('METAL_BASE') ), + A_CLASSER ( ('METAL_BASE'), ('TRANSITOIRE') ) + ) + ) # Fin JDC_CATA + +# -------------------------------------------------- +# fin entete +# -------------------------------------------------- + + + + + + +#================================ +# 1. Definition des OPTIONS +#================================ + +OPTIONS = PROC ( nom = "OPTIONS", + op = 68, + repetable = 'n', + fr = "Définitions des options", + +#=== +# Liste des paramètres +#=== + +# INCRTPS + IncrementTemporel = SIMP ( statut = "o", + typ = "I", + defaut = "1", + max = 1, + #val_max = 100, + fr = "Incrément temporel (=1 pour calcul déterministe)", + ), + +# DTPREC + IncrementMaxTemperature = SIMP ( statut = "o", + typ = "R", + defaut = "0.1", + max = 1, + val_min = 0., + val_max = 1., + fr = "Incrément maximum d'évolution de la température par noeud et par instant (°C)", + ), + +# DTARCH + IncrementMaxTempsAffichage = SIMP ( statut = "o", + typ = "R", + defaut = "1000.", + max = 1, + val_min = 0., + val_max = 1000., + fr = "Incrément maximum de temps pour l'affichage (s)", + ), + +# NBO +# Question : NBO depend-il de TYPGEOM ?? + NombreNoeudsMaillage = SIMP ( statut = "o", + typ = "R", + max=1, + val_min = 1., + val_max = 100., + fr = "Nombre de noeuds à considérer dans le maillage interne", + ), + +# + ListeInstants = SIMP ( statut = "o", + typ = "R", + max = "**", + fr = "Liste des instants pour lesquels la température et les contraintes seront archivés", + ), + +) # Fin PROC OPTIONS + +#================================ +# 2. Caracteristiques du DEFAUT +#================================ + +DEFAUT = PROC ( nom = "DEFAUT", + op = 68, + repetable = 'n', + fr = "Caractéristiques du défaut", + +#=== +# Liste des paramètres +#=== + +# TYPEDEF + TypeInitial = SIMP ( statut = "o", typ = "TXM", + into = ( "Sous Revetement", # DSR + "Debouchant", # DD + ), + fr = "Type initial du défaut : sous revêtement ou débouchant", + ), + +#==== +# 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 défaut", + ), + + # PROFDEF + ProfondeurRadiale = SIMP ( statut = "o", + typ = "R", + #defaut = "0.006", + max = 1, + val_max = 1., + fr = "Profondeur radiale du défaut (m)", + ), + + # OPTLONG + OptionCalculLongueur = SIMP ( statut = "o", + typ = "TXM", + into = ( "Valeur", # VALEUR + "Relation lineaire avec la longueur" ), # RAPPORT + #defaut = "VALEUR", + fr = "Option pour caractériser la longueur du défaut : 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 défaut sous revêtement (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 défaut sous revêtement", + ), + + ), # Fin BLOC Option_Rapport + + # DECADEF + DecalageRadial = SIMP ( statut = "o", + typ = "R", + #defaut = "-0.00001", + fr = "Décalage radial du defaut sous revêtement (m)", + ), + + # ANGLDEF + Azimut = SIMP ( statut = "o", + typ = "R", + defaut = "0.", + fr = "Position angulaire du défaut dans le cylindre de cuve (degrés)", + ), + + # ALTIDEF + Altitude = SIMP ( statut = "o", + typ = "R", + defaut = "2.", + val_min = 0., + val_max = 4., + fr = "Altitude du défaut sur le cylindre de cuve (m)", + ), + + # POINDEF + Pointe = SIMP ( statut = "o", + typ = "TXM", + into = ( "A", "B" ), + defaut = "A", + fr = "Choix du point considéré du défaut sous revêtement", + ), + + # ARRETFISSURE + ArretDeFissure = SIMP ( statut = "o", + typ = "TXM", + into = ( "OUI", "NON" ), + defaut = "NON", + fr = "Prise en compte de l'arrêt de fissure", + ), + + # INCRDEF + IncrementTailleFissure = SIMP ( statut = "o", + typ = "R", + defaut = "0.005", + fr = "Incrément 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'intensité 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 défaut : longitudinale ou circonférentielle", + ), + + # PROFDEF + ProfondeurRadiale = SIMP ( statut = "o", + typ = "R", + #defaut = "0.006", + max = 1, + val_max = 1., + fr = "Profondeur radiale du défaut (m)", + ), + + # ANGLDEF + Azimut = SIMP ( statut = "o", + typ = "R", + defaut = "0.", + fr = "Position angulaire du défaut dans le cylindre de la cuve (en degrés)", + ), + + # ALTIDEF + Altitude = SIMP ( statut = "o", + typ = "R", + defaut = "2.", + fr = "Altitude du défaut 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'arrêt de fissure", + ), + + # INCRDEF + IncrementTailleFissure = SIMP ( statut = "o", + typ = "R", + defaut = "0.005", + fr = "Incrément 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'intensité 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'intensité de contraintes", + ), + + ), # Fin BLOC debouchant + +) # Fin PROC DEFAUT + + +#================================ +# 3. Caracteristiques de la CUVE +#================================ + +CUVE = PROC (nom = "CUVE", + op = 68, + repetable = 'n', + fr = "Caractéristiques 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 définition de la geométrie 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 revêtement (m)", + ), + + # LIGMIN + LigamentExterneMin = SIMP ( statut = "o", + typ = "R", + defaut = "0.75", + fr = "Ligament externe minimal avant rupture (% de l'épaisseur 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)", + ), + ), # Fin BLOC Maillage + +) # Fin PROC CUVE + +#==================================================== +# 4. Modeles de fluence, d'irradiation et de tenacite +#==================================================== + +#======================= +# 4.1 Modeles de fluence +#======================= + +MODELES = PROC ( nom = "MODELES", + op = 68, + repetable = 'n', + fr = "Modèles de fluence, d'irradiation et de ténacité", + + +#=== +# 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 = "Modèle d'atténuation de la fluence dans l'épaisseur 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 assimilée 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 assimilée 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 = "Paramètre exponentiel du modèle 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 assimilée 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 assimilée 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 assimilée 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 = "Paramètre exponentiel du modèle 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 assimilée 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 assimilée 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 assimilée 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 assimilée 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 à l'azimut 0 (10^19 n/cm2)", + ), + + # COEFFLUENCE2 + FluenceAzimut5 = SIMP ( statut = "o", + typ = "R", + defaut = "5.48", + fr = "Fluence à l'azimut 5 (10^19 n/cm2)", + ), + + # COEFFLUENCE3 + FluenceAzimut10 = SIMP ( statut = "o", + typ = "R", + defaut = "4.46", + fr = "Fluence à l'azimut 10 (10^19 n/cm2)", + ), + + # COEFFLUENCE4 + FluenceAzimut15 = SIMP ( statut = "o", + typ = "R", + defaut = "3.41", + fr = "Fluence à l'azimut 15 (10^19 n/cm2)", + ), + + # COEFFLUENCE5 + FluenceAzimut20 = SIMP ( statut = "o", + typ = "R", + defaut = "3.37", + fr = "Fluence à l'azimut 20 (10^19 n/cm2)", + ), + + # COEFFLUENCE6 + FluenceAzimut25 = SIMP ( statut = "o", + typ = "R", + defaut = "3.16", + fr = "Fluence à l'azimut 25 (10^19 n/cm2)", + ), + + # COEFFLUENCE7 + FluenceAzimut30 = SIMP ( statut = "o", + typ = "R", + defaut = "2.74", + fr = "Fluence à l'azimut 30 (10^19 n/cm2)", + ), + + # COEFFLUENCE8 + FluenceAzimut35 = SIMP ( statut = "o", + typ = "R", + defaut = "2.25", + fr = "Fluence à l'azimut 35 (10^19 n/cm2)", + ), + + # COEFFLUENCE9 + FluenceAzimut40 = SIMP ( statut = "o", + typ = "R", + defaut = "1.89", + fr = "Fluence à l'azimut 40 (10^19 n/cm2)", + ), + + # COEFFLUENCE10 + FluenceAzimut45 = SIMP ( statut = "o", + typ = "R", + defaut = "1.78", + fr = "Fluence à l'azimut 45 (10^19 n/cm2)", + ), + + ), # 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 de la cuve à l'instant de l'analyse (°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 = "Modèle d'irradiation pour virole ou joint soudé", + ), + + # 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 décalage de RTNDT (°F) (28. pour js et 17. pour mdb)", + ), + + # nbectDRTNDT + NombreEcartTypeRTndt = SIMP ( statut = "o", + typ = "R", + defaut = "2.", + fr = "Nombre d'écart-type par rapport à 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 = "Modèle de ténacité", + ), + +#==== +# 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # KICPAL + PalierDuctile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "195.", + fr = "Palier déterministe de K1c quand modèle RCC-M (MPa(m^0.5)) ", + ), + + # KICCDVD = SIMP ( statut = "o", + CoefficientVariation_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "0.15", + fr = "Coefficient de variation de la loi normale de K1c quand modèle RCC-M ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à la moyenne de KIa (nb sigma) ", + ), + + + # KIAPAL + PalierDuctile_KIa = SIMP ( statut = "o", + typ = "R", + defaut = "195.", + fr = "Palier déterministe de K1a -ténacite à l'arrêt- quand modèle RCC-M (MPa(m^0.5)) ", + ), + + # KIACDV + CoefficientVariation_KIa = SIMP ( statut = "o", + typ = "R", + defaut = "0.10", + fr = "Coefficient de variation de la loi normale de K1a -ténacite à l'arrêt- quand modèle 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Paramètre T0 du modèle Wallin (°C) : température pour laquelle la téncité du matériau vaut en moyenne 100MPa.m^5", + ), + + # nbectKIc + NbEcartType_MoyKIc = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caracteristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caracteristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caracteristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à 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 = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caracteristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à la moyenne de KIa (nb sigma) ", + ), + + ), # Fin BLOC Parametres_LOGWOLF + + Weibull_Generalisee = BLOC ( condition = " Tenacite in ( 'Weibull generalisee',) ", + + Coefficients = FACT ( statut = "o", + + # A1 + A1 = SIMP ( statut = "o", + typ = "R", + defaut = "21.263", + fr = "coef du coef a(T) d'une Weibull générale", + ), + + # A2 + A2 = SIMP ( statut = "o", + typ = "R", + defaut = "9.159", + fr = "coef du coef a(T) d'une Weibull générale", + ), + + # A3 + A3 = SIMP ( statut = "o", + typ = "R", + defaut = "0.04057", + fr = "coef du coef a(T) d'une Weibull générale", + ), + + # B1 + B1 = SIMP ( statut = "o", + typ = "R", + defaut = "17.153", + fr = "coef du coef b(T) d'une Weibull générale", + ), + + # B2 + B2 = SIMP ( statut = "o", + typ = "R", + defaut = "55.089", + fr = "coef du coef b(T) d'une Weibull générale", + ), + + # B3 + B3 = SIMP ( statut = "o", + typ = "R", + defaut = "0.0144", + fr = "coef du coef b(T) d'une Weibull générale", + ), + + # C1 + C1 = SIMP ( statut = "o", + typ = "R", + defaut = "4.", + fr = "coef du coef c(T) d'une Weibull générale", + ), + + # C2 + C2 = SIMP ( statut = "o", + typ = "R", + defaut = "0.", + fr = "coef du coef c(T) d'une Weibull générale", + ), + + # C3 + C3 = SIMP ( statut = "o", + typ = "R", + defaut = "0.", + fr = "coef du coef c(T) d'une Weibull générale", + ), + + ), # FIN FACT Coefficients + + # nbectKIc + NbEcartType_MoyKIc = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", + ), + + # fractKIc + Fractile_KIc = SIMP ( statut = "o", + typ = "R", + defaut = "5.", + fr = "Valeur caracteristique de KIc exprimée en ordre de fractile (%) ", + ), + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Nombre d'écart-type par rapport à la moyenne de KIa (nb sigma) ", + ), + + ), # Fin BLOC Parametres_WEIB_GEN + +) # Fin PROC MODELES + + +#================== +# 5. Initialisation +#================== + +INITIALISATION = PROC ( nom = "INITIALISATION", + op = 68, + repetable = 'n', + fr = "Initialisation : instant initial, profils radiaux de température et contraintes", + + TemperatureInitiale = FACT ( statut = "o", + + ProfilRadial_TemperatureInitiale = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Profil radial de la température initiale dans la cuve (m) (°C) ", + ), + + Amont_TemperatureInitiale = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_TemperatureInitiale = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT TemperatureInitiale + + ContraintesInitiales = FACT ( statut = "o", + + ProfilRadial_ContraintesInitiales = SIMP ( statut = "o", + typ = Tuple(4), + max = "**", + fr = "Profil radial des contraintes radiale, circonférentielle et longitudinale dans la cuve (m) (xx) (xx) (xx) ", + ), + + Amont_ContraintesInitiales = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ContraintesInitiales = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT ContraintesInitiales + + # INSTINIT + InstantInitialisation = SIMP ( statut = "o", + typ = "R", + defaut = "-1.", + fr = "Instant initial auquel sont définies la température, ainsi que les contraintes initiales (en s) ", + ), + +) # Fin PROC INITIALISATION + + +#================================== +# 6. CARACTERISTIQUES DU REVETEMENT +#================================== + +REVETEMENT = PROC ( nom = "REVETEMENT", + op = 68, + repetable = 'n', + fr = "Caracteristiques du revêtement", + + # KTHREV + ConditionLimiteThermiqueREV = SIMP ( statut = "o", + typ = "TXM", + into = ( "ENTHALPIE", "CHALEUR",), + #defaut = "CHALEUR", + fr = "Option pour définir les caractéristiques du revêtement ", + ), + + EnthalpieREV = BLOC ( condition = " ConditionLimiteThermiqueREV in ( 'ENTHALPIE', ) ", + + EnthalpieREV_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / enthalpie massique (J/kg) ", + ), + + Amont_EnthalpieREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_EnthalpieREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC EnthalpieREV + + + ChaleurREV = BLOC ( condition = " ConditionLimiteThermiqueREV in ( 'CHALEUR', ) ", + + ChaleurREV_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / chaleur volumique (J/kg/K) ", + ), + + Amont_ChaleurREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ChaleurREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC ChaleurREV + + ConductiviteREV = FACT (statut = "o", + + ConductiviteREV_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / conductivité thermique (W/m/°C) ", + ), + + Amont_ConductiviteREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ConductiviteREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT ConductiviteREV + + ModuleYoungREV = FACT (statut = "o", + + ModuleYoungREV_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / module d'Young (MPa) ", + ), + + Amont_ModuleYoungREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ModuleYoungREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT ModuleYoungREV + + CoeffDilatThermREV = FACT (statut = "o", + + CoeffDilatThermREV_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / coefficient de dilatation thermique (°C-1) ", + ), + + Amont_CoeffDilatThermREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_CoeffDilatThermREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT CoeffDilatThermREV + + LimiteElasticiteREV = FACT (statut = "o", + + LimiteElasticiteREV_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / limite d'élasticite (MPa) ", + ), + + Amont_LimiteElasticiteREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_LimiteElasticiteREV = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT LimiteElasticiteREV + + AutresParametresREV = FACT (statut = "o", + + # TREFREV + TemperatureDeformationNulleREV = SIMP ( statut = "o", + typ = "R", + defaut = "20.", + fr = "Température de référence pour laquelle les déformations thermiques sont nulles (°C) ", + ), + + # TDETREV + TemperaturePourCoefDilatThermREV = SIMP ( statut = "o", + typ = "R", + defaut = "287.", + fr = "Température de définition du coefficient de dilatation thermique (°C) ", + ), + + # NUREV + CoefficientPoissonREV = SIMP ( statut = "o", + typ = "R", + defaut = "0.3", + fr = "Coefficient de Poisson ", + ), + + ), # Fin FACT AutresParametresREV + +) # Fin PROC REVETEMENT + + +#===================================== +# 7. CARACTERISTIQUES DU METAL DE BASE +#===================================== + +METAL_BASE = PROC ( nom = "METAL_BASE", + op = 68, + repetable = 'n', + fr = "Caracteristiques du metal de base", + + # KTHMDB + ConditionLimiteThermiqueMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( "ENTHALPIE", "CHALEUR",), + #defaut = "CHALEUR", + fr = "Option pour definir les caractéristiques du revêtement ", + ), + + EnthalpieMDB = BLOC ( condition = " ConditionLimiteThermiqueMDB in ( 'ENTHALPIE', ) ", + + EnthalpieMDB_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / enthalpie massique (J/kg) ", + ), + + Amont_EnthalpieMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_EnthalpieMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC EnthalpieMDB + + ChaleurMDB = BLOC ( condition = " ConditionLimiteThermiqueMDB in ( 'CHALEUR', ) ", + + ChaleurMDB_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / chaleur volumique (J/kg/K) ", + ), + + Amont_ChaleurMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ChaleurMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC ChaleurMDB + + ConductiviteMDB = FACT ( statut = "o", + + ConductiviteMDB_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / conductivité thermique (W/m/°C) ", + ), + + Amont_ConductiviteMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ConductiviteMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT ConductiviteMDB + + ModuleYoungMDB = FACT ( statut = "o", + + ModuleYoungMDB_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / module d'Young (MPa) ", + ), + + Amont_ModuleYoungMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_ModuleYoungMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT ModuleYoungMDB + + CoeffDilatThermMDB = FACT ( statut = "o", + + CoeffDilatThermMDB_Fct_Temperature = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Température (°C) / coefficient de dilatation thermique (°C-1) ", + ), + + Amont_CoeffDilatThermMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_CoeffDilatThermMDB = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin FACT CoeffDilatThermMDB + + AutresParametresMDB = FACT ( statut = "o", + + # TREFMDB + TemperatureDeformationNulleMDB = SIMP ( statut = "o", + typ = "R", + defaut = "20.", + fr = "Température de référence pour laquelle les déformations thermiques sont nulles (°C) ", + ), + + # TDETMDB + TemperaturePourCoefDilatThermMDB = SIMP ( statut = "o", + typ = "R", + defaut = "287.", + fr = "Température de définition du coefficient de dilatation thermique (°C) ", + ), + + # NUMDB + CoefficientPoissonMDB = SIMP ( statut = "o", + typ = "R", + defaut = "0.3", + fr = "Coefficient de Poisson ", + ), + + ), # Fin FACT TemperatureDeformationNulleMDB + +) # Fin PROC METAL_BASE + + +#=============================== +# 8. TRANSITOIRE THERMOMECANIQUE +#=============================== + +TRANSITOIRE = PROC ( nom = "TRANSITOIRE", + op = 68, + repetable = 'n', + fr = "Description du transitoire thermohydraulique", + + Pression = FACT ( statut = "o", + + ProfilTemporel_Pression = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / pression (MPa) ", + ), + + Amont_Pression = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_Pression = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # FIN FACT Pression + + # TYPCLTH + TypeConditionLimiteThermique = SIMP ( statut = "o", + typ = "TXM", + into = ( "Temperature imposee en paroi", # TEMP_IMPO + "Flux de chaleur impose en paroi", # FLUX_REP + "Temperature imposee du fluide et coefficient echange", # ECHANGE + "Debit massique et temperature d injection de securite", # DEBIT + "Temperature imposee du fluide et debit d injection de securite"), # TEMP_FLU + #defaut = "ECHANGE", + fr = "Type de condition thermique en paroi interne ", + ), + + TemperatureImposeeParoi = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee en paroi', ) ", + + ProfilTemporel_TemperatureImposeeParoi = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Température imposée (°C) ", + ), + + Amont_TemperatureImposeeParoi = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_TemperatureImposeeParoi = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC TemperatureImposeeParoi + + FluxChaleur = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Flux de chaleur impose en paroi', ) ", + + ProfilTemporel_FluxChaleur = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Flux de chaleur impose (W/m2) ", + ), + + Amont_FluxChaleur = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_FluxChaleur = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC FluxChaleur + + TemperatureImposeeFluide = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', 'Temperature imposee du fluide et debit d injection de securite', ) ", + + ProfilTemporel_TemperatureImposeeFluide = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Température imposée (°C) ", + ), + + Amont_TemperatureImposeeFluide = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_TemperatureImposeeFluide = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC TemperatureImposeeFluide + + CoefficientEchange = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', ) ", + + ProfilTemporel_CoefficientEchange = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Coefficient d'échange (W/m2/K) ", + ), + + Amont_CoefficientEchange = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_CoefficientEchange = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC CoefficientEchange + + DebitMassique = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ", + + ProfilTemporel_DebitMassique = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Débit massique (kg/s) ", + ), + + Amont_DebitMassique = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_DebitMassique = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC DebitMassique + + TemperatureInjection = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ", + + ProfilTemporel_TemperatureInjection = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Température d'injection de sécurité (°C) ", + ), + + Amont_TemperatureInjection = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_TemperatureInjection = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + ), # Fin BLOC TemperatureInjection + + Creare = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ", + + # DH + DiametreHydraulique = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Diamètre 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 caractéristique convection naturelle (m) ", + ), + + # ALPHA_CF + CoeffVestale_ConvectionForcee = SIMP ( statut = "o", + typ = "R", + defaut = "1.", + fr = "Coefficient Vestale convection forcée (-) ", + ), + + # ALPHA_CN + CoeffVestale_ConvectionNaturelle = SIMP ( statut = "o", + typ = "R", + defaut = "1.", + fr = "Coefficient Vestale convection naturelle (-) ", + ), + + # EPS + CritereConvergenceRelative = SIMP ( statut = "o", + typ = "R", + defaut = "0.00001", + fr = "Critère convergence relative (-) ", + ), + + # VM + VolumeMelange_CREARE = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Volume de mélange CREARE (m3) ", + ), + + # T0 + TemperatureInitiale_CREARE = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Température initiale CREARE (°C) ", + ), + + # SE + SurfaceEchange_FluideStructure = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Surface d'échange fluide/structure (m2) ", + ), + + ), # Fin BLOC Creare + + + DebitInjection = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ", + + ProfilTemporel_DebitInjection = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Instant (s) / Débit d'injection de sécurité (kg/s) ", + ), + + Amont_DebitInjection = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière amont", + ), + + Aval_DebitInjection = SIMP ( statut = "o", + typ = "TXM", + into = ( 'Continu', 'Exclu', 'Lineaire' ), + fr = "Prolongation à la frontière aval", + ), + + ), # Fin BLOC DebitInjection + + + Vestale = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ", + + # DH + DiametreHydraulique = SIMP ( statut = "o", + typ = "R", + defaut = "-2.", + fr = "Diamètre 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 caractéristique convection naturelle (m) ", + ), + + # ALPHA_CF + CoeffVestale_ConvectionForcee = SIMP ( statut = "o", + typ = "R", + defaut = "1.", + fr = "Coefficient d'échange Vestale convection forcée (-) ", + ), + + # ALPHA_CN + CoeffVestale_ConvectionNaturelle = SIMP ( statut = "o", + typ = "R", + defaut = "1.", + fr = "Coefficient d'échange Vestale convection naturelle (-) ", + ), + + # EPS + CritereConvergenceRelative = SIMP ( statut = "o", + typ = "R", + defaut = "0.00001", + fr = "Critère convergence relative (-) ", + ), + + ), # Fin BLOC Vestale + +) # Fin PROC TRANSITOIRE diff --git a/Cuve2dg/Cuve2dg_Cata_V3.py b/Cuve2dg/Cuve2dg_Cata_V3.py new file mode 100644 index 00000000..cc5fc92a --- /dev/null +++ b/Cuve2dg/Cuve2dg_Cata_V3.py @@ -0,0 +1,1050 @@ +# -*- coding: utf-8 -*- + +# -------------------------------------------------- +# Definition de variables sous forme de tuple +# -------------------------------------------------- + +import types +class Tuple: + def __init__(self,ntuple): + self.ntuple=ntuple + + def __convert__(self,valeur): + if type(valeur) == types.StringType: + return None + if len(valeur) != self.ntuple: + return None + return valeur + + def info(self): + return "Tuple de %s elements" % self.ntuple + + __repr__=info + __str__=info + +# -------------------------------------------------- +# 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'), + A_CLASSER ( ('OPTIONS'), ('DEFAUT') ), + A_CLASSER ( ('DEFAUT'), ('CUVE') ), + A_CLASSER ( ('CUVE'), ('MODELES') ), + A_CLASSER ( ('MODELES'), ('INITIALISATION') ), + A_CLASSER ( ('INITIALISATION'), ('REVETEMENT') ), + A_CLASSER ( ('REVETEMENT'), ('METAL_BASE') ), + A_CLASSER ( ('METAL_BASE'), ('TRANSITOIRE') ) + ) + ) # Fin JDC_CATA + +# -------------------------------------------------- +# fin entete +# -------------------------------------------------- + +# -------------------------------------------------- +# RESPONSABLE D. THAI VAN +# Ce fichier contient la liste des coefficients pour un +# modele de Weibull generalise +def Coef_WeibGen() : return FACT(statut='o',min=1,max='**', + + # A1 + A1 = SIMP ( statut="o", typ="R", defaut=21.263, + fr="coef du coef a(T) d'une Weibull générale", ), + # A2 + A2 = SIMP ( statut="o", typ="R", defaut=9.159, + fr="coef du coef a(T) d'une Weibull générale", ), + # A3 + A3 = SIMP ( statut="o", typ="R", defaut=0.04057, + fr="coef du coef a(T) d'une Weibull générale", ), + # B1 + B1 = SIMP ( statut="o", typ="R", defaut=17.153, + fr="coef du coef b(T) d'une Weibull générale", ), + # B2 + B2 = SIMP ( statut="o", typ="R", defaut=55.089, + fr="coef du coef b(T) d'une Weibull générale", ), + # B3 + B3 = SIMP ( statut="o", typ="R", defaut=0.0144, + fr="coef du coef b(T) d'une Weibull générale", ), + # C1 + C1 = SIMP ( statut="o", typ="R", defaut=4., + fr="coef du coef c(T) d'une Weibull générale", ), + # C2 + C2 = SIMP ( statut="o", typ="R", defaut=0., + fr="coef du coef c(T) d'une Weibull générale", ), + # C3 + C3 = SIMP ( statut="o", typ="R", defaut=0., + fr="coef du coef c(T) d'une Weibull générale", ), + +); # FIN def Coef_WeibGen + + +# -------------------------------------------------- +# RESPONSABLE D. THAI VAN +# Ce fichier contient la liste des coefficients +def Coef_Fluence() : return FACT(statut='o',min=1,max='**', + + # COEFFLUENCE1 + Azimut_0deg = SIMP ( statut="o", typ="R", defaut=5.8, + fr="Fluence à l'azimut 0 (10^19 n/cm2)", ), + # COEFFLUENCE2 + Azimut_5deg = SIMP ( statut="o", typ="R", defaut=5.48, + fr="Fluence à l'azimut 5 (10^19 n/cm2)", ), + # COEFFLUENCE3 + Azimut_10deg = SIMP ( statut="o", typ="R", defaut=4.46, + fr="Fluence à l'azimut 10 (10^19 n/cm2)", ), + # COEFFLUENCE4 + Azimut_15deg = SIMP ( statut="o", typ="R", defaut=3.41, + fr="Fluence à l'azimut 15 (10^19 n/cm2)", ), + # COEFFLUENCE5 + Azimut_20deg = SIMP ( statut="o", typ="R", defaut=3.37, + fr="Fluence à l'azimut 20 (10^19 n/cm2)", ), + # COEFFLUENCE6 + Azimut_25deg = SIMP ( statut="o", typ="R", defaut=3.16, + fr="Fluence à l'azimut 25 (10^19 n/cm2)", ), + # COEFFLUENCE7 + Azimut_30deg = SIMP ( statut="o", typ="R", defaut=2.74, + fr="Fluence à l'azimut 30 (10^19 n/cm2)", ), + # COEFFLUENCE8 + Azimut_35deg = SIMP ( statut="o", typ="R", defaut=2.25, + fr="Fluence à l'azimut 35 (10^19 n/cm2)", ), + # COEFFLUENCE9 + Azimut_40deg = SIMP ( statut="o", typ="R", defaut=1.89, + fr="Fluence à l'azimut 40 (10^19 n/cm2)", ), + # COEFFLUENCE10 + Azimut_45deg = SIMP ( statut="o", typ="R", defaut=1.78, + fr="Fluence à l'azimut 45 (10^19 n/cm2)", ), + +); # FIN def Coef_Fluence + +#================================ +# 1. Definition des OPTIONS +#================================ + +OPTIONS = PROC ( nom = "OPTIONS", + op = 68, + repetable = 'n', + fr = "Définitions des options", + +#=== +# Liste des paramètres +#=== + +# INCRTPS + IncrementTemporel = SIMP ( statut="o", typ="I", max=1, defaut=1, + fr="Incrément temporel (=1 pour calcul déterministe)", ), + +# DTPREC + IncrementMaxTemperature = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1., defaut=0.1, + fr="Incrément maximum d'évolution de la température par noeud et par instant (°C)", ), + +# DTARCH + IncrementMaxTempsAffichage = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1000., defaut="1000.", + fr="Incrément maximum de temps pour l'affichage (s)", ), + +# NBO +# Question : NBO depend-il de TYPGEOM ?? + NombreNoeudsMaillage = SIMP ( statut="o", typ="R", max=1, val_min=1., val_max=1000., + fr = "Nombre de noeuds à considérer dans le maillage interne", ), + +# + ListeInstants = SIMP ( statut="o", typ="R", max="**", + fr = "Liste des instants pour lesquels la température et les contraintes seront archivés", ), + +) # Fin PROC OPTIONS + +#================================ +# 2. Caracteristiques du DEFAUT +#================================ + +DEFAUT = PROC ( nom = "DEFAUT", + op = 68, + repetable = 'n', + fr = "Caractéristiques du défaut", + +#=== +# Liste des paramètres +#=== + +# TYPEDEF + TypeInitial = SIMP ( statut="o", typ="TXM", + fr="Type initial du défaut : sous revêtement ou débouchant", + into=( "Sous Revetement", # DSR + "Debouchant", ), # DD + ), + +#==== +# Definition des parametres selon le type du defaut +#==== + + SousRevetement = BLOC ( condition = "TypeInitial=='Sous Revetement'", + + # ORIEDEF into LONGITUD, CIRCONF + Orientation = SIMP ( statut="o", typ="TXM", + fr="Orientation du défaut", + into=( "Longitudinale", + "Circonferentielle" ), ), + + # PROFDEF + ProfondeurRadiale = SIMP ( statut="o", typ="R", max=1, val_min=1., val_max=1., defaut=0.006, + fr="Profondeur radiale du défaut (m)", ), + + # OPTLONG into VALEUR, RAPPORT + OptionCalculLongueur = SIMP ( statut="o", typ="TXM", + fr="Option pour caractériser la longueur du défaut : soit par valeur, soit par un rapport LONG/PROF", + into = ( "Valeur", + "Relation lineaire avec la longueur" ), ), + + Option_Valeur = BLOC ( condition = "OptionCalculLongueur=='Valeur'", + # LONGDEF + Longueur = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1., defaut = 0.060, + fr = "Longueur du défaut sous revêtement (m)", ), + ), # Fin BLOC Option_Valeur + + Option_Rapport = BLOC ( condition = "OptionCalculLongueur=='Relation lineaire avec la longueur'", + # LONGSURPROF + LongueurSurProfondeur = SIMP ( statut="o", typ="R", max=1, val_max=100., + fr="Rapport longueur/profondeur du défaut sous revêtement", ), + ), # Fin BLOC Option_Rapport + + # DECADEF + DecalageRadial = SIMP ( statut="o", typ="R", defaut=-0.00001, + fr="Décalage radial du defaut sous revêtement (en m)", ), + # ANGLDEF + Azimut = SIMP ( statut="o", typ="R", defaut=0., + fr="Position angulaire du défaut dans le cylindre de cuve (en degrés)", ), + # ALTIDEF + Altitude = SIMP ( statut="o", typ="R", val_min=0., val_max=4., defaut=2., + fr="Altitude du défaut sur le cylindre de cuve (en m)", ), + # POINDEF + Pointe = SIMP ( statut="o", typ="TXM", defaut="A", + fr="Choix du point considéré du défaut sous revêtement", + into=( "A", "B" ), ), + # CORRECPLASTIC + CorrectionPlastiqueBeta = SIMP ( statut="o", typ="TXM", defaut="NON", + fr="Prise en compte de la correction plastique BETA dans le calcul du facteur d'intensité de contraintes", + into=( "OUI", "NON" ), ), + + ), # Fin BLOC SousRevetement + + Debouchant = BLOC ( condition = "TypeInitial=='Debouchant'", + + # ORIEDEF into LONGITUD, CIRCONF + Orientation = SIMP ( statut="o", typ="TXM", + fr="Orientation du défaut : longitudinale ou circonférentielle", + into=( "Longitudinale", + "Circonferentielle" ), ), + # PROFDEF + ProfondeurRadiale = SIMP ( statut="o", typ="R", max=1, val_min=0., val_max=1., defaut=0.006, + fr="Profondeur radiale du défaut (en m)", ), + # ANGLDEF + Azimut = SIMP ( statut="o", typ="R", defaut=0., + fr="Position angulaire du défaut dans le cylindre de la cuve (en degrés)", ), + # ALTIDEF + Altitude = SIMP ( statut="o", typ="R", defaut=2., + fr="Altitude du défaut sur le cylindre de la cuve (en m)", ), + # IRWIN + CorrectionPlastiqueIrwin = SIMP ( statut="o", typ="TXM", defaut="NON", + fr="Prise en compte de la correction plastique d'Irwin dans le calcul du facteur d'intensité de contraintes", + into=( "OUI", "NON" ), ), + # CORRECPLASTIC + CorrectionPlastiqueBeta = SIMP ( statut="o", typ="TXM", defaut="NON", + fr="Prise en compte de la correction plastique BETA dans le calcul du facteur d'intensité de contraintes", + into=( "OUI", "NON" ), ), + + ), # Fin BLOC debouchant + +) # Fin PROC DEFAUT + + +#================================ +# 3. Caracteristiques de la CUVE +#================================ + +CUVE = PROC (nom = "CUVE", + op = 68, + repetable = 'n', + fr = "Caractéristiques de la cuve", + +#=== +# Liste des paramètres +#=== + + # TYPEGEOM + TraitementGeometrie = SIMP ( statut="o", typ="TXM", + fr="Choix de la définition de la geométrie d'une cuve", + into=( "Topologie", # GEOMETRIE + "Maillage"), # MAILLAGE + ), + +#==== +# Definition des parametres selon le type de traitement de la geometrie +#==== + + Geometrie = BLOC ( condition = "TraitementGeometrie=='Topologie'", + + # RINT + RayonInterne = SIMP ( statut="o", typ="R", val_min=0., defaut=1.994, + fr="Rayon interne de la cuve (en m)", ), + + # DTV : comment preciser que REXT > RINT ? + # REXT + RayonExterne = SIMP ( statut="o", typ="R", defaut=2.2015, + fr="Rayon externe de la cuve (en m)", ), + + # DTV : comment preciser que LREV < RINT ? + # LREV + EpaisseurRevetement = SIMP ( statut="o", typ="R", defaut=0.0075, + fr="Epaisseur du revêtement (m)", ), + + # LIGMIN + LigamentExterneMin = SIMP ( statut="o", typ="R", defaut=0.75, + fr="Ligament externe minimal avant rupture (% de l'épaisseur de cuve)", ), + + ), # Fin BLOC Geometrie + + Maillage = BLOC ( condition = "TraitementGeometrie=='Maillage'", + + # DTV : comment preciser que c'est une suite de liste de nombres croissants ? + # Question : quel rapport avec NBO ?? + Liste_abscisses = SIMP ( statut="o", typ="R", max="**", + fr = "Liste des abscisses (m)", ), + ), # Fin BLOC Maillage + +) # Fin PROC CUVE + +#==================================================== +# 4. Modeles de fluence, d'irradiation et de tenacite +#==================================================== + +#======================= +# 4.1 Modeles de fluence +#======================= + +MODELES = PROC ( nom = "MODELES", + op = 68, + repetable = 'n', + fr = "Modèles de fluence, d'irradiation et de ténacité", + + +#=== +# Liste des paramètres +#=== + + Fluence = FACT ( statut="o", + + # MODELFLUENCE + Modele = SIMP ( statut="o", typ="TXM", + fr="Modèle d'atténuation de la fluence dans l'épaisseur de la cuve", + 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 + ), + +#==== +# Definition des parametres selon le modele de fluence +#==== + + Reglementaire = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=9.7 (Reglementaire)', ) ", + # DTV : comment proposer une liste de valeurs, tout en proposant de fournir d'autres valeurs ? + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + ), # Fin BLOC Reglementaire + + France = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=12.7 (France)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + # KPFRANCE + KPFrance = SIMP ( statut="o", typ="R", + fr="Paramètre exponentiel du modèle France", + defaut = 12.7, ), + ), # Fin BLOC France + + ValeurImposee = BLOC ( condition = " Modele in ( 'Exponentiel sans revetement k=0. (ValeurImposee)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + ), # Fin BLOC ValeurImposee + + SDM = BLOC ( condition = " Modele in ( 'Donnees francaises du palier CPY (SDM)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + ), # Fin BLOC SDM + + USNRC = BLOC ( condition = " Modele in ( 'Regulatory Guide 1.99 rev 2 (USNRC)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + KPUS = SIMP ( statut="o", typ="R", + fr="Paramètre exponentiel du modèle US", + defaut=9.4488, ), + ), # Fin BLOC USNRC + + REV_2 = BLOC ( condition = " Modele in ( 'Dossier 900 MWe AP9701 rev 2 (REV_2)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + ), # Fin BLOC REV_2 + + SDM_Lissage = BLOC ( condition = " Modele in ( 'Lissage du modele ajuste (SDM_Lissage)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + ), # Fin BLOC SDM_Lissage + + GrandeDev = BLOC ( condition = " Modele in ( 'Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + ), # Fin BLOC GrandeDev + + GD_Cuve = BLOC ( condition = " Modele in ( 'Grand developpement (GD_Cuve)', ) ", + # fmax + FluenceMax = SIMP ( statut="o", typ="R", defaut=6.5, + fr="Fluence maximale en surface interne assimilée par la cuve (10^19 n/cm2) ; HP-26/99/045 : p.31 : fluence max = 7.3/9.125/10/95", ), + Coefficients = Coef_Fluence(), + ), # Fin BLOC GD_Cuve + + ), # Fin FACT Fluence + +#========================== +# 4.2 Modeles d'irradiation +#========================== + + Irradiation = FACT ( statut="o", + + # TYPIRR + Option = SIMP ( statut = "o", typ = "TXM", + fr = "Type d'irradiation", + into = ( "RTndt de la cuve a l instant de l analyse", # RTNDT + "Modele d irradiation" ), # FLUENCE + ), + +#==== +# Definition des parametres selon le type d'irradiation +#==== + + RTndt = BLOC ( condition = "Option=='RTndt de la cuve a l instant de l analyse'", + + RTNDT = SIMP ( statut="o", typ="R", defaut=73., + fr="RTNDT de la cuve à l'instant de l'analyse (°C)", ), + + ), # Fin BLOC RTndt + + ModeleIrradiation = BLOC ( condition = "Option=='Modele d irradiation'", + + # MODELIRR + Modele = SIMP ( statut="o", typ="TXM", + fr="Modèle d'irradiation pour virole ou joint soudé", + 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 + ), + # 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 = " Modele 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 décalage de RTNDT (°F) (28. pour js et 17. pour mdb)", ), + # nbectDRTNDT + NombreEcartTypeRTndt = SIMP ( statut="o", typ="R", defaut=2., + fr="Nombre d'écart-type par rapport à la moyenne de DRTNDT", ), + ), # Fin BLOC Parametres_USNRC + + ), # Fin BLOC ModeleIrradiation + + ), # Fin FACT Irradiation + +#======================== +# 4.3 Modeles de tenacite +#======================== + + Tenacite = FACT ( statut = "o", + + # MODELKIC + Modele = SIMP ( statut="o", typ="TXM", + fr="Modèle de calcul de la ténacité à l'amorçage KIc", + 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 + ), + + +#==== +# Definition des parametres selon le modele de tenacité +#==== + +# Modeles type RCC-M + + KIc_RCCM = BLOC ( condition = " Modele in ( 'RCC-M/ASME coefficient=2', 'RCC-M/ASME coefficient=2.33 (Houssin)', 'RCC-M/ASME avec KI=KIpalier', 'RCC-M/ASME avec KI~exponentiel', ) ", + + # nbectKIc + NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2., + fr = "Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ), + + # KICPAL + PalierDuctile_KIc = SIMP ( statut="o", typ="R", defaut=195. , + fr="Palier déterministe de K1c ou valeur du palier ductile plafonnant la courbe (en MPa(m^0.5)) ", ), + + # KICCDVD + CoefficientVariation_KIc = SIMP ( statut="o", typ="R", defaut = 0.15, + fr = "Coefficient de variation de la loi normale de K1c ", ), + + ), # Fin BLOC KIc_RCCM + + Fissure = BLOC ( condition = " Modele in ( 'RCC-M/ASME coefficient=2', 'RCC-M/ASME coefficient=2.33 (Houssin)', 'RCC-M/ASME avec KI=KIpalier', 'RCC-M/ASME avec KI~exponentiel', )", + + # ARRETFISSURE + ArretDeFissure = SIMP ( statut="o", typ="TXM", defaut="NON", + fr="Prise en compte de l'arrêt de fissure", + into=( "OUI", "NON" ), ), + + Parametres_Fissure = BLOC ( condition = "ArretDeFissure=='OUI'", + # INCRDEF + IncrementTailleFissure = SIMP ( statut="o", typ="R", defaut=0.005, + fr="Incrément de la taille de fissure pour la propagation (en m)", ), + ), # Fin BLOC Parametres_Fissure + + KIa_RCCM = BLOC ( condition = "ArretDeFissure=='OUI'", + + # nbectKIa + NbEcartType_MoyKIa = SIMP ( statut="o", typ="R", defaut=-2., + fr="Nombre d'écart-type par rapport à la moyenne de KIa (nb sigma) ", ), + + # KIACDV + CoefficientVariation_KIa = SIMP ( statut="o", typ="R", defaut=0.10, + fr="Coefficient de variation de la loi normale de K1a -ténacite à l'arrêt- ", ), + + # KIAPAL + PalierDuctile_KIa = SIMP ( statut="o", typ="R", defaut=195., + fr="Palier déterministe de K1a -ténacite à l'arrêt- (en MPa(m^0.5)) ", ), + ), # Fin BLOC KIa_RCCM + + ), # Fin BLOC Fissure + +# Modeles type Weibull + + KIc_MasterCuve = BLOC ( condition = " Modele in ( 'Weibull basee sur la master cuve', ) ", + + # T0WALLIN + Temperature_KIc100 = SIMP ( statut="o", typ="I", defaut=-27, + fr="Paramètre T0 du modèle Wallin (°C) : température pour laquelle la téncité du matériau vaut en moyenne 100MPa.m^5", ), + + # fractKIc + Fractile_KIc = SIMP ( statut="o", typ="R", defaut=5., + fr="Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", ), + + ), # Fin BLOC KIc_MasterCuve + + KIc_Weibull1 = BLOC ( condition = " Modele in ( 'Weibull basee sur la master cuve (REME)', 'Weibull n°1 (etude ORNL)', ) ", + + # fractKIc + Fractile_KIc = SIMP ( statut="o", typ="R", defaut=5., + fr="Valeur caractéristique de KIc exprimée en ordre de fractile (%) ", ), + + ), # Fin BLOC KIc_Weibull1 + + KIc_Weibull2 = BLOC ( condition = " Modele in ( 'Weibull n°2', 'Weibull n°3', ) ", + + # nbectKIc + NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2., + fr="Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ), + + ), # Fin BLOC KIc_Weibull2 + + Weibull_Generalisee = BLOC ( condition = " Modele in ( 'Weibull generalisee',) ", + + Coefficients = Coef_WeibGen(), + # nbectKIc + NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2., + fr="Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ), + + ), # Fin BLOC Weibull_Generalisee + +# Modeles type exponentiel (Frama, LOGWOLF) + + KIc_Exponentielle = BLOC ( condition = " Modele in ( 'Exponentielle n°1 (Frama)', 'Exponentielle n°2 (LOGWOLF)', ) ", + + # nbectKIc + NbEcartType_MoyKIc = SIMP ( statut="o", typ="R", defaut=-2., + fr="Nombre d'écart-type par rapport à la moyenne de KIc (nb sigma) : det = -2 ", ), + + ), # Fin BLOC KIc_Exponentielle + + ), # Fin FACT Tenacite + +) # Fin PROC MODELES + + +#================== +# 5. Initialisation +#================== + +INITIALISATION = PROC ( nom = "INITIALISATION", + op = 68, + repetable = 'n', + fr = "Initialisation : instant initial, profils radiaux de température et contraintes", + + TemperatureInitiale = FACT ( statut = "o", + + ProfilRadial_TemperatureInitiale = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Profil radial de la température initiale dans la cuve (en m : °C) ", ), + Amont_TemperatureInitiale = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_TemperatureInitiale = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT TemperatureInitiale + + ContraintesInitiales = FACT ( statut = "o", + + ProfilRadial_ContraintesInitiales = SIMP ( statut="o", typ=Tuple(4), max="**", + fr="Profil radial des contraintes radiale, circonférentielle et longitudinale dans la cuve (en m : xx : xx : xx) ", ), + Amont_ContraintesInitiales = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ContraintesInitiales = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT ContraintesInitiales + + # INSTINIT + InstantInitialisation = SIMP ( statut="o", typ="R", defaut = -1., + fr="Instant initial auquel sont définies la température, ainsi que les contraintes initiales (en s) ", ), + +) # Fin PROC INITIALISATION + + +#================================== +# 6. CARACTERISTIQUES DU REVETEMENT +#================================== + +REVETEMENT = PROC ( nom = "REVETEMENT", + op = 68, + repetable = 'n', + fr = "Caracteristiques du revêtement", + + # KTHREV + ConditionLimiteThermiqueREV = SIMP ( statut="o", typ="TXM", + fr="Option pour définir les caractéristiques du revêtement ", + into=( "ENTHALPIE", "CHALEUR",), + ), + + EnthalpieREV = BLOC ( condition = "ConditionLimiteThermiqueREV=='ENTHALPIE'", + + EnthalpieREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / enthalpie massique (J/kg) ", ), + Amont_EnthalpieREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_EnthalpieREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC EnthalpieREV + + + ChaleurREV = BLOC ( condition = "ConditionLimiteThermiqueREV=='CHALEUR'", + + ChaleurREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / chaleur volumique (J/kg/K) ", ), + Amont_ChaleurREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ChaleurREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC ChaleurREV + + ConductiviteREV = FACT (statut = "o", + + ConductiviteREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / conductivité thermique (W/m/°C) ", ), + Amont_ConductiviteREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ConductiviteREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT ConductiviteREV + + ModuleYoungREV = FACT (statut = "o", + + ModuleYoungREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / module d'Young (MPa) ", ), + Amont_ModuleYoungREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ModuleYoungREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT ModuleYoungREV + + CoeffDilatThermREV = FACT (statut = "o", + + CoeffDilatThermREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / coefficient de dilatation thermique (°C-1) ", ), + Amont_CoeffDilatThermREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_CoeffDilatThermREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT CoeffDilatThermREV + + LimiteElasticiteREV = FACT (statut = "o", + + LimiteElasticiteREV_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / limite d'élasticite (MPa) ", ), + Amont_LimiteElasticiteREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_LimiteElasticiteREV = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT LimiteElasticiteREV + + AutresParametresREV = FACT (statut = "o", + + # TREFREV + TemperatureDeformationNulleREV = SIMP ( statut="o", typ="R", defaut=20., + fr="Température de référence pour laquelle les déformations thermiques sont nulles (°C) ", ), + # TDETREV + TemperaturePourCoefDilatThermREV = SIMP ( statut="o", typ="R", defaut=287., + fr="Température de définition du coefficient de dilatation thermique (°C) ", ), + # NUREV + CoefficientPoissonREV = SIMP ( statut="o", typ="R", defaut=0.3, + fr="Coefficient de Poisson ", ), + + ), # Fin FACT AutresParametresREV + +) # Fin PROC REVETEMENT + + +#===================================== +# 7. CARACTERISTIQUES DU METAL DE BASE +#===================================== + +METAL_BASE = PROC ( nom = "METAL_BASE", + op = 68, + repetable = 'n', + fr = "Caracteristiques du metal de base", + + # KTHMDB + ConditionLimiteThermiqueMDB = SIMP ( statut="o", typ="TXM", + fr="Option pour definir les caractéristiques du revêtement ", + into=( "ENTHALPIE", "CHALEUR",), ), + + EnthalpieMDB = BLOC ( condition = "ConditionLimiteThermiqueMDB=='ENTHALPIE'", + + EnthalpieMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / enthalpie massique (J/kg) ", ), + Amont_EnthalpieMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_EnthalpieMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC EnthalpieMDB + + ChaleurMDB = BLOC ( condition = "ConditionLimiteThermiqueMDB=='CHALEUR'", + + ChaleurMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / chaleur volumique (J/kg/K) ", ), + Amont_ChaleurMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ChaleurMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC ChaleurMDB + + ConductiviteMDB = FACT ( statut = "o", + + ConductiviteMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / conductivité thermique (W/m/°C) ", ), + Amont_ConductiviteMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ConductiviteMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT ConductiviteMDB + + ModuleYoungMDB = FACT ( statut="o", + + ModuleYoungMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / module d'Young (MPa) ", ), + Amont_ModuleYoungMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_ModuleYoungMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT ModuleYoungMDB + + CoeffDilatThermMDB = FACT ( statut="o", + + CoeffDilatThermMDB_Fct_Temperature = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Température (°C) / coefficient de dilatation thermique (°C-1) ", ), + Amont_CoeffDilatThermMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_CoeffDilatThermMDB = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin FACT CoeffDilatThermMDB + + AutresParametresMDB = FACT ( statut = "o", + + # TREFMDB + TemperatureDeformationNulleMDB = SIMP ( statut="o", typ="R", defaut=20., + fr="Température de référence pour laquelle les déformations thermiques sont nulles (°C) ", ), + # TDETMDB + TemperaturePourCoefDilatThermMDB = SIMP ( statut="o", typ="R", defaut=287., + fr="Température de définition du coefficient de dilatation thermique (°C) ", ), + # NUMDB + CoefficientPoissonMDB = SIMP ( statut="o", typ="R", defaut=0.3, + fr="Coefficient de Poisson ", ), + + ), # Fin FACT TemperatureDeformationNulleMDB + +) # Fin PROC METAL_BASE + + +#=============================== +# 8. TRANSITOIRE THERMOMECANIQUE +#=============================== + +TRANSITOIRE = PROC ( nom = "TRANSITOIRE", + op = 68, + repetable = 'n', + fr = "Description du transitoire thermohydraulique", + + Pression = FACT ( statut = "o", + + ProfilTemporel_Pression = SIMP ( statut="o", typ=Tuple(2), max="**", + fr = "Instant (s) / pression (MPa) ", ), + Amont_Pression = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_Pression = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # FIN FACT Pression + + # TYPCLTH + TypeConditionLimiteThermique = SIMP ( statut="o", typ="TXM", + fr="Type de condition thermique en paroi interne ", + into=( "Temperature imposee en paroi", # TEMP_IMPO + "Flux de chaleur impose en paroi", # FLUX_REP + "Temperature imposee du fluide et coefficient echange", # ECHANGE + "Debit massique et temperature d injection de securite", # DEBIT + "Temperature imposee du fluide et debit d injection de securite"), # TEMP_FLU + ), + + TemperatureImposeeParoi = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee en paroi', ) ", + + ProfilTemporel_TemperatureImposeeParoi = SIMP ( statut="o", typ=Tuple(2), max="**", + fr = "Instant (s) / Température imposée (°C) ", ), + Amont_TemperatureImposeeParoi = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_TemperatureImposeeParoi = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC TemperatureImposeeParoi + + FluxChaleur = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Flux de chaleur impose en paroi', ) ", + + ProfilTemporel_FluxChaleuri = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Instant (s) / Flux de chaleur impose (W/m2) ", ), + Amont_FluxChaleur = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_FluxChaleur = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC FluxChaleur + + TemperatureImposeeFluide = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', 'Temperature imposee du fluide et debit d injection de securite', ) ", + + ProfilTemporel_TemperatureImposeeFluide = SIMP ( statut="o", typ=Tuple(2), max="**", + fr = "Instant (s) / Température imposée (°C) ", ), + Amont_TemperatureImposeeFluide = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_TemperatureImposeeFluide = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC TemperatureImposeeFluide + + CoefficientEchange = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et coefficient echange', ) ", + + ProfilTemporel_CoefficientEchange = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Instant (s) / Coefficient d'échange (W/m2/K) ", ), + Amont_CoefficientEchange = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_CoefficientEchange = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC CoefficientEchange + + DebitMassique = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ", + + ProfilTemporel_DebitMassique = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Instant (s) / Débit massique (kg/s) ", ), + Amont_DebitMassique = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_DebitMassique = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC DebitMassique + + TemperatureInjection = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ", + + ProfilTemporel_TemperatureInjection = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Instant (s) / Température d'injection de sécurité (°C) ", ), + Amont_TemperatureInjection = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_TemperatureInjection = SIMP ( statut="o", typ = "TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + ), # Fin BLOC TemperatureInjection + + Creare = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Debit massique et temperature d injection de securite', ) ", + + # DH + DiametreHydraulique = SIMP ( statut="o", typ="R", defaut=-2., + fr="Diamètre 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 caractéristique convection naturelle (m) ", ), + # ALPHA_CF + CoeffVestale_ConvectionForcee = SIMP ( statut="o", typ="R", defaut=1., + fr="Coefficient Vestale convection forcée (-) ", ), + # ALPHA_CN + CoeffVestale_ConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=1., + fr="Coefficient Vestale convection naturelle (-) ", ), + # EPS + CritereConvergenceRelative = SIMP ( statut="o", typ="R", defaut=0.00001, + fr = "Critère convergence relative (-) ", ), + # VM + VolumeMelange_CREARE = SIMP ( statut="o", typ="R", defaut=-2., + fr = "Volume de mélange CREARE (m3) ", ), + # T0 + TemperatureInitiale_CREARE = SIMP ( statut="o", typ="R", defaut=-2., + fr="Température initiale CREARE (°C) ", ), + # SE + SurfaceEchange_FluideStructure = SIMP ( statut="o", typ="R", defaut=-2., + fr="Surface d'échange fluide/structure (m2) ", ), + + ), # Fin BLOC Creare + + DebitInjection = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ", + + ProfilTemporel_DebitInjection = SIMP ( statut="o", typ=Tuple(2), max="**", + fr="Instant (s) / Débit d'injection de sécurité (kg/s) ", ), + Amont_DebitInjection = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière amont", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + Aval_DebitInjection = SIMP ( statut="o", typ="TXM", + fr="Prolongation à la frontière aval", + into=( 'Continu', 'Exclu', 'Lineaire' ), ), + + ), # Fin BLOC DebitInjection + + + Vestale = BLOC ( condition = " TypeConditionLimiteThermique in ( 'Temperature imposee du fluide et debit d injection de securite', ) ", + + # DH + DiametreHydraulique = SIMP ( statut="o", typ="R", defaut=-2., + fr="Diamètre 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 caractéristique convection naturelle (m) ", ), + # ALPHA_CF + CoeffVestale_ConvectionForcee = SIMP ( statut="o", typ="R", defaut=1., + fr="Coefficient d'échange Vestale convection forcée (-) ", ), + # ALPHA_CN + CoeffVestale_ConvectionNaturelle = SIMP ( statut="o", typ="R", defaut=1., + fr="Coefficient d'échange Vestale convection naturelle (-) ", ), + # EPS + CritereConvergenceRelative = SIMP ( statut="o", typ="R", defaut=0.00001, + fr="Critère convergence relative (-) ", ), + + ), # Fin BLOC Vestale + +) # Fin PROC TRANSITOIRE diff --git a/Cuve2dg/Cuve2dg_Cata_VPN.py b/Cuve2dg/Cuve2dg_Cata_VPN.py new file mode 100644 index 00000000..852c6efd --- /dev/null +++ b/Cuve2dg/Cuve2dg_Cata_VPN.py @@ -0,0 +1,111 @@ +# -*- coding: utf-8 -*- +import types +class Tuple: + def __init__(self,ntuple): + self.ntuple=ntuple + + def __convert__(self,valeur): + if type(valeur) == types.StringType: + return None + if len(valeur) != self.ntuple: + return None + return valeur + + def info(self): + return "Tuple de %s elements" % self.ntuple + + __repr__=info + __str__=info + + +# -------------------------------------------------- +# debut entete +# -------------------------------------------------- + +import Accas +from Accas import * + + +#CONTEXT.debug = 1 +JdC = JDC_CATA ( code = 'CUVE1D-DEFAILLGLOB', + execmodul = None, + regles = ( AU_MOINS_UN ('OPTIONS','DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE'), + A_CLASSER ( ('OPTIONS',), ('DEFAUT', 'CUVE', 'MODELES', 'INITIALISATION', 'REVETEMENT', 'METAL_BASE', 'TRANSITOIRE')) + ) + ) # 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 = PROC ( nom = "OPTIONS", + 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 = Tuple(2), + max = "**", + fr = "Liste des instants ", + validators=VerifTypeTuple(('R','R')), + ), + + ListeInstants2 = SIMP ( statut = "o", + typ = Tuple(2), + max = "**", + fr = "Liste des instants ", + ), + +) # Fin OPER OPTIONS diff --git a/Cuve2dg/generator_cuve2dg.py b/Cuve2dg/generator_cuve2dg.py new file mode 100644 index 00000000..bfb6b444 --- /dev/null +++ b/Cuve2dg/generator_cuve2dg.py @@ -0,0 +1,1537 @@ +# -*- coding: utf-8 -*- +# CONFIGURATION MANAGEMENT OF EDF VERSION +# ====================================================================== +# COPYRIGHT (C) 1991 - 2002 EDF R&D WWW.CODE-ASTER.ORG +# THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY +# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY +# THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR +# (AT YOUR OPTION) ANY LATER VERSION. +# +# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT +# WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF +# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU +# GENERAL PUBLIC LICENSE FOR MORE DETAILS. +# +# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE +# ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER, +# 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE. +# +# +# ====================================================================== +""" + Ce module contient le plugin generateur de fichier au format + DefaillCUVE pour EFICAS. + +""" +import traceback +import types,string,re + +from Noyau import N_CR +from Accas import MCSIMP +from generator_python import PythonGenerator + +def entryPoint(): + """ + Retourne les informations nécessaires pour le chargeur de plugins + + Ces informations sont retournées dans un dictionnaire + """ + return { + # Le nom du plugin + 'name' : 'cuve2dg', + # La factory pour créer une instance du plugin + 'factory' : Cuve2dgGenerator, + } + + +class Cuve2dgGenerator(PythonGenerator): + """ + Ce generateur parcourt un objet de type JDC et produit + un texte au format eficas et + un texte au format DefaillCUVE + + """ + # Les extensions de fichier préconisées + extensions=('.comm',) + + def __init__(self,cr=None): + # Si l'objet compte-rendu n'est pas fourni, on utilise le compte-rendu standard + if cr : + self.cr=cr + else: + self.cr=N_CR.CR(debut='CR generateur format DefaillCUVE pour DefaillCUVE', + fin='fin CR format DefaillCUVE pour DefaillCUVE') + # Le texte au format DefaillCUVE est stocké dans l'attribut textCuve + self.textCuve='' + + # Ce dictionnaire liste le nom des variables utilisees dans le script + self.variable = { + "NiveauImpression" : "MESSAGE_LEVEL", + "FichierDataIn" : "DATARESUME_FILE", + "FichierTempSigma" : "TEMPSIG_FILE", + "FichierResultats" : "RESU_FILE", + "FichierCSV" : "CSV_FILE", + "FichierCREARE" : "CREARE_FILE", + "GrandeurEvaluee" : "GRANDEUR", + "IncrementTemporel" : "INCRTPS", + "IncrementMaxTemperature" : "DTPREC", + "IncrementMaxTempsAffichage" : "DTARCH", + "TraitementGeometrie" : "TYPEGEOM", + "RayonInterne" : "RINT", + "RayonInterne_mess" : "RINT_MESSAGE", + "RayonExterne" : "REXT", + "RayonExterne_mess" : "REXT_MESSAGE", + "EpaisseurRevetement" : "LREV", + "EpaisseurRevetement_mess" : "LREV_MESSAGE", + "LigamentExterneMin" : "LIGMIN", + "LigamentExterneMin_mess" : "LIGMIN_MESSAGE", + "NombreNoeudsMaillage" : "NBNO", + "TypeInitial" : "TYPEDEF", + "Orientation" : "ORIEDEF", + "ProfondeurRadiale" : "PROFDEF", + "ProfondeurRadiale_mess" : "PROFDEF_MESSAGE", + "ModeCalculLongueur" : "OPTLONG", + "Longueur" : "LONGDEF", + "Longueur_mess" : "LONGDEF_MESSAGE", + "CoefDirecteur" : "PROFSURLONG", + "CoefDirecteur_mess" : "PROFSURLONG_MESSAGE", + "Constante" : "LONGCONST", + "ModeCalculDecalage" : "DECATYP", + "DecalageNormalise" : "DECANOR", + "DecalageNormalise_mess" : "DECANOR_MESSAGE", + "DecalageRadial" : "DECADEF", + "DecalageRadial_mess" : "DECADEF_MESSAGE", + "Azimut" : "ANGLDEF", + "Azimut_mess" : "ANGLDEF_MESSAGE", + "Altitude_mess" : "ANGLDEF_MESSAGE", + "Altitude" : "ALTIDEF", + "Altitude_mess" : "ALTIDEF_MESSAGE", + "Pointe" : "POINDEF", + "ModeleFluence" : "MODELFLUENCE", + "ZoneActiveCoeur_AltitudeSup" : "H1COEUR", + "ZoneActiveCoeur_AltitudeInf" : "H2COEUR", + "FluenceMax" : "fmax", + "KPFrance" : "KPFRANCE", + "KPUS" : "KPUS", + "Azimut_0deg" : "COEFFLUENCE1", + "Azimut_5deg" : "COEFFLUENCE2", + "Azimut_10deg" : "COEFFLUENCE3", + "Azimut_15deg" : "COEFFLUENCE4", + "Azimut_20deg" : "COEFFLUENCE5", + "Azimut_25deg" : "COEFFLUENCE6", + "Azimut_30deg" : "COEFFLUENCE7", + "Azimut_35deg" : "COEFFLUENCE8", + "Azimut_40deg" : "COEFFLUENCE9", + "Azimut_45deg" : "COEFFLUENCE10", + "TypeIrradiation" : "TYPEIRR", + "RTNDT" : "RTNDT", + "ModeleIrradiation" : "MODELIRR", + "TeneurCuivre" : "CU", + "TeneurCuivre_mess" : "CU_MESSAGE", + "TeneurNickel" : "NI", + "TeneurNickel_mess" : "NI_MESSAGE", + "TeneurPhosphore" : "P", + "TeneurPhosphore_mess" : "P_MESSAGE", + "MoyenneRTndt" : "RTimoy", + "MoyenneRTndt_mess" : "RTimoy_MESSAGE", + "CoefVariationRTndt" : "RTicov", + "CoefVariationRTndt_mess" : "RTicov_MESSAGE", + "EcartTypeRTndt" : "USectDRT", + "EcartTypeRTndt_mess" : "USectDRT_MESSAGE", + "NombreEcartTypeRTndt" : "nbectDRTNDT", + "NombreEcartTypeRTndt_mess" : "nbectDRTNDT_MESSAGE", + "ModeleTenacite" : "MODELKIC", + "NBRE_CARACTERISTIQUE" : "NBCARAC", + "NbEcartType_MoyKIc" : "nbectKIc", + "NbEcartType_MoyKIc_mess" : "nbectKIc_MESSAGE", + "PalierDuctile_KIc" : "KICPAL", + "CoefficientVariation_KIc" : "KICCDV", + "Fractile_KIc" : "fractKIc", + "Fractile_KIc_mess" : "fractKIc_MESSAGE", + "Temperature_KIc100" : "T0WALLIN", + "A1" : "A1", + "A2" : "A2", + "A3" : "A3", + "B1" : "B1", + "B2" : "B2", + "B3" : "B3", + "C1" : "C1", + "C2" : "C2", + "C3" : "C3", + "AttnCorrBeta" : "ATTNCORRBETA", + "CorrIrwin" : "CORRIRWIN", + "ArretDeFissure" : "ARRETFISSURE", + "IncrementTailleFissure" : "INCRDEF", + "IncrementTailleFissure_mess" : "INCRDEF_MESSAGE", + "NbEcartType_MoyKIa" : "nbectKIa", + "PalierDuctile_KIa" : "KIAPAL", + "CoefficientVariation_KIa" : "KIACDV", + "InstantInitialisation" : "INSTINIT", + "ConditionLimiteThermiqueREV" : "KTHREV", + "TemperatureDeformationNulleREV" : "TREFREV", + "TemperaturePourCoefDilatThermREV" : "TDETREV", + "CoefficientPoissonREV" : "NUREV", + "ConditionLimiteThermiqueMDB" : "KTHMDB", + "TemperatureDeformationNulleMDB" : "TREFMDB", + "TemperaturePourCoefDilatThermMDB" : "TDETMDB", + "CoefficientPoissonMDB" : "NUMDB", + "TypeConditionLimiteThermique" : "TYPCLTH", + "Instant_1" : "INSTANT1", + "Instant_2" : "INSTANT2", + "DebitAccumule" : "QACCU", + "DebitInjectionSecurite" : "QIS", + "TempInjectionSecurite_mess" : "TIS_MESSAGE", + "DiametreHydraulique" : "DH", + "DiametreHydraulique_mess" : "DH_MESSAGE", + "SectionEspaceAnnulaire" : "SECTION", + "SectionEspaceAnnulaire_mess" : "SECTION_MESSAGE", + "HauteurCaracConvectionNaturelle" : "DELTA", + "HauteurCaracConvectionNaturelle_mess" : "DELTA_MESSAGE", + "CritereConvergenceRelative" : "EPS", + "CoefficientsVestale" : "COEFVESTALE", + "VolumeMelange_CREARE" : "VM", + "VolumeMelange_CREARE_mess" : "VM_MESSAGE", + "TemperatureInitiale_CREARE" : "T0", + "TemperatureInitiale_CREARE_mess" : "T0_MESSAGE", + "SurfaceEchange_FluideStructure" : "SE", + "SurfaceEchange_FluideStructure_mess" : "SE_MESSAGE", + } + + # Ce dictionnaire liste le commentaire des variables utilisees dans le script + self.comment = { + "NiveauImpression" : "Niveau d impression des messages a l ecran (=0 : rien, =1 : temps calcul total, =2 : temps intermediaires)", + "FichierDataIn" : "sortie du fichier recapitulatif des donnees d entree {OUI ; NON}", + "FichierTempSigma" : "sortie des fichiers temperature et contraintes {OUI ; NON}", + "FichierResultats" : "sortie du fichier resultat template_DEFAILLCUVE {OUI ; NON}", + "FichierCSV" : "sortie du fichier resultat template_DEFAILLCUVE.CSV {OUI ; NON}", + "FichierCREARE" : "sortie du fichier Tfluide et Coef Echange {OUI ; NON}", + "GrandeurEvaluee" : "choix de la grandeur sous critere evaluee {FM_KICSURKCP ; MARGE_KI ; MARGE_KCP}", + "IncrementTemporel" : "increment temporel pour l analyse PROBABILISTE (si DETERMINISTE, fixer a 1)", + "IncrementMaxTemperature" : "increment max de temp/noeud/instant (°C)", + "IncrementMaxTempsAffichage" : "increment max de temps pour affichage (s)", + "TraitementGeometrie" : "traitement de la geometrie de la cuve : {GEOMETRIE, MAILLAGE}", + "RayonInterne" : "rayon interne (m)", + "RayonInterne_mess" : "affichage ecran du rayon interne (m)", + "RayonExterne" : "rayon externe (m)", + "RayonExterne_mess" : "affichage ecran du rayon externe (m)", + "EpaisseurRevetement" : "epaisseur revetement (m)", + "EpaisseurRevetement_mess" : "affichage ecran de l epaisseur revetement (m)", + "LigamentExterneMin" : "ligament externe minimal avant rupture (% de l'epaisseur de cuve)", + "LigamentExterneMin_mess" : "affichage ecran du ligament externe minimal avant rupture (% de l'epaisseur de cuve)", + "NombreNoeudsMaillage" : "nbre de noeuds dans l'épaisseur de la cuve", + "TypeInitial" : "type initial du defaut : DEBOUCHANT=Defaut Debouchant, DSR=Defaut Sous Revetement, DECALE=Defaut Decale", + "Orientation" : "orientation (LONGITUD / CIRCONF)", + "ProfondeurRadiale" : "profondeur radiale ou encore hauteur (m)", + "ProfondeurRadiale_mess" : "affichage ecran de la profondeur radiale ou encore hauteur (m)", + "ModeCalculLongueur" : "option pour definir la longueur du defaut (VALEUR pour une valeur fixe, FCTAFFINE pour une fct affine de la profondeur)", + "Longueur" : "longueur (m) pour defaut Sous Revetement", + "Longueur_mess" : "affichage ecran de la longueur (m) pour defaut Sous Revetement", + "CoefDirecteur" : "pente de la fonction affine l = h/profsurlong + a0", + "CoefDirecteur_mess" : "affichage ecran de la pente de la fonction affine l = h/profsurlong + a0", + "Constante" : "constante de la fonction affine a0", + "ModeCalculDecalage" : "type de decalage : normalise (NORMALISE) ou reel (VALEUR)", + "DecalageNormalise" : "decalage radial normalise (valeur comprise entre 0. et 1.) pour defaut Sous Revetement", + "DecalageNormalise_mess" : "affichage ecran du decalage radial normalise (valeur comprise entre 0. et 1.) pour defaut Sous Revetement", + "DecalageRadial" : "decalage radial reel (m) pour defaut decale", + "DecalageRadial_mess" : "affichage ecran du decalage radial reel (m) pour defaut decale", + "Azimut" : "coordonnee angulaire (degre)", + "Azimut_mess" : "affichage ecran de la coordonnee angulaire (degre)", + "Altitude" : "altitude (m) : valeur négative", + "Altitude_mess" : "affichage ecran de l altitude (m) : valeur négative", + "Pointe" : "choix du(des) point(s) du defaut considere {'A','B','BOTH'} pour DSR et DECALE (pour DEBOUCHANT : automatiquement 'B')", + "ModeleFluence" : "modele de fluence : {Reglementaire, France, ValeurImposee, SDM, USNRC, REV_2, SDM_Lissage, GrandeDev, GD_Cuve, Cuve1D}", + "ZoneActiveCoeur_AltitudeSup" : "cote superieure de la zone active de coeur (ici pour cuve palier 900Mw)", + "ZoneActiveCoeur_AltitudeInf" : "cote inferieure de la zone active de coeur (ici pour cuve palier 900Mw)", + "FluenceMax" : "fluence maximale assimilee par la cuve (n/cm2)", + "KPFrance" : "parametre exponentiel du modele France", + "KPUS" : "parametre exponentiel du modele US", + "Azimut_0deg" : "fluence a l'azimut 0 (10^19 n/cm)", + "Azimut_5deg" : "fluence a l'azimut 5 (10^19 n/cm)", + "Azimut_10deg" : "fluence a l'azimut 10 (10^19 n/cm)", + "Azimut_15deg" : "fluence a l'azimut 15 (10^19 n/cm)", + "Azimut_20deg" : "fluence a l'azimut 20 (10^19 n/cm)", + "Azimut_25deg" : "fluence a l'azimut 25 (10^19 n/cm)", + "Azimut_30deg" : "fluence a l'azimut 30 (10^19 n/cm)", + "Azimut_35deg" : "fluence a l'azimut 35 (10^19 n/cm)", + "Azimut_40deg" : "fluence a l'azimut 40 (10^19 n/cm)", + "Azimut_45deg" : "fluence a l'azimut 45 (10^19 n/cm)", + "TypeIrradiation" : "type irradiation : {RTNDT, FLUENCE}", + "RTNDT" : "RTNDT finale (°C)", + "ModeleIrradiation" : "modele d irradiation : {HOUSSIN, PERSOZ, LEFEBVRE, USNRCmdb} pour virole et {BRILLAUD,USNRCsoud} pour jointsoude", + "TeneurCuivre" : "teneur en cuivre (%)", + "TeneurCuivre_mess" : "affichage ecran de la teneur en cuivre (%)", + "TeneurNickel" : "teneur en nickel (%)", + "TeneurNickel_mess" : "affichage ecran de la teneur en nickel (%)", + "TeneurPhosphore" : "teneur en phosphore (%)", + "TeneurPhosphore_mess" : "affichage ecran de la teneur en phosphore (%)", + "MoyenneRTndt" : "moyenne de la RTNDT initiale : virole C1 de cuve Chinon : mdb=>-17.°C et js=>42.°C (HT-56/05/038 : p.52)", + "MoyenneRTndt_mess" : "affichage ecran de la moyenne de la RTNDT initiale", + "CoefVariationRTndt" : "coef de variation de la RTNDT initiale", + "CoefVariationRTndt_mess" : "affichage ecran du coef de variation de la RTNDT initiale", + "EcartTypeRTndt" : "pour modeles USNRCsoud ou USNRCmdb, ecart-type du decalage de RTNDT (°F) (28. pour js et 17. pour mdb)", + "EcartTypeRTndt_mess" : "affichage ecran, pour modeles USNRCsoud ou USNRCmdb, ecart-type du decalage de RTNDT (°F) (28. pour js et 17. pour mdb)", + "NombreEcartTypeRTndt" : "Nbre d ecart-type par rapport a la moyenne de DRTNDT si analyse PROBABILISTE (en DETERMINISTE, fixer a 2.)", + "NombreEcartTypeRTndt_mess" : "affichage ecran du nbre d ecart-type par rapport a la moyenne de DRTNDT si analyse PROBABILISTE", + "ModeleTenacite" : "modele de tenacite : {RCC-M, RCC-M_pal, RCC-M_exp, RCC-M_simpl, Houssin_RC, Wallin, REME, ORNL, Frama, WEIB3, WEIB2, LOGWOLF, WEIB-GEN}", + "NBRE_CARACTERISTIQUE" : "Nb caracteristique : ORDRE ou QUANTILE", + "NbEcartType_MoyKIc" : "Nbre d ecart-type par rapport a la moyenne de KIc si analyse PROBABILISTE (en DETERMINISTE, fixer a -2.)", + "NbEcartType_MoyKIc_mess" : "affichage ecran du nbre d ecart-type par rapport a la moyenne de KIc si analyse PROBABILISTE", + "PalierDuctile_KIc" : "palier deterministe de K1c (MPa(m^0.5))", + "CoefficientVariation_KIc" : "coef de variation de la loi normale de K1c", + "Fractile_KIc" : "valeur caracteristique de KIc exprimee en ordre de fractile (%)", + "Fractile_KIc_mess" : "affichage ecran de la valeur caracteristique de KIc exprimee en ordre de fractile (%)", + "Temperature_KIc100" : "parametre T0 du modele Wallin (°C)", + "A1" : "coef des coefs d une WEIBULL generale", + "A2" : "", + "A3" : "", + "B1" : "", + "B2" : "", + "B3" : "", + "C1" : "", + "C2" : "", + "C3" : "", + "AttnCorrBeta" : "Attenuation de la correction plastique : {OUI, NON} ==> uniquement pour DSR ou DECALE", + "CorrIrwin" : "Correction plastique IRWIN : {OUI, NON} ==> uniquement pour DEBOUCHANT", + "ArretDeFissure" : "prise en compte de l arret de fissure {OUI, NON} (en PROBABILISTE, fixer a NON)", + "IncrementTailleFissure" : "increment de la taille de fissure (m)", + "IncrementTailleFissure_mess" : "affichage ecran de l increment de la taille de fissure (m)", + "NbEcartType_MoyKIa" : "Nbre d ecart-type par rapport a la moyenne de KIa (nb sigma)", + "PalierDuctile_KIa" : "palier deterministe de K1a quand modele RCC-M (MPa(m^0.5))", + "CoefficientVariation_KIa" : "coef de variation de la loi normale de K1a", + "InstantInitialisation" : "instant initial (s)", + "ConditionLimiteThermiqueREV" : "Option 'ENTHALPIE' ou 'CHALEUR'", + "TemperatureDeformationNulleREV" : "temperature de deformation nulle (°C)", + "TemperaturePourCoefDilatThermREV" : "temperature de definition du coefficient de dilatation thermique (°C)", + "CoefficientPoissonREV" : "coefficient de Poisson", + "ConditionLimiteThermiqueMDB" : "Option 'ENTHALPIE' ou 'CHALEUR'", + "TemperatureDeformationNulleMDB" : "temperature de deformation nulle (°C)", + "TemperaturePourCoefDilatThermMDB" : "temperature de definition du coefficient de dilatation thermique (°C)", + "CoefficientPoissonMDB" : "coefficient de Poisson", + "TypeConditionLimiteThermique" : "Type de condition thermique en paroi interne {TEMP_IMPO,FLUX_REP,ECHANGE,DEBIT,TEMP_FLU,APRP}", + "Instant_1" : "Borne inferieure de l intervalle de temps du 2nd palier T1", + "Instant_2" : "Borne superieure de l intervalle de temps du 2nd palier T1", + "DebitAccumule" : "Debit accumule (en m3/h)", + "DebitInjectionSecurite" : "Debit injection de securite (en m3/h)", + "TempInjectionSecurite_mess" : "affichage ecran de la temperature injection de securite", + "DiametreHydraulique" : "Diametre hydraulique (m)", + "DiametreHydraulique_mess" : "affichage ecran du diametre hydraulique (m)", + "SectionEspaceAnnulaire" : "Section espace annulaire (m2)", + "SectionEspaceAnnulaire_mess" : "affichage ecran de la section espace annulaire (m2)", + "HauteurCaracConvectionNaturelle" : "Hauteur caracteristique convection naturelle (m)", + "HauteurCaracConvectionNaturelle_mess" : "affichage ecran de la hauteur caracteristique convection naturelle (m)", + "CritereConvergenceRelative" : "Critere convergence relative (-)", + "CoefficientsVestale" : "Application des coefs de Vestale {OUI;NON}", + "VolumeMelange_CREARE" : "Volume de melange CREARE (m3)", + "VolumeMelange_CREARE_mess" : "affichage ecran du volume de melange CREARE (m3)", + "TemperatureInitiale_CREARE" : "Temperature initiale CREARE (degC)", + "TemperatureInitiale_CREARE_mess" : "affichage ecran de la temperature initiale CREARE (degC)", + "SurfaceEchange_FluideStructure" : "Surface d'echange fluide/structure (m2)", + "SurfaceEchange_FluideStructure_mess" : "affichage ecran de la surface d'echange fluide/structure (m2)", + } + + # Ce dictionnaire liste la valeur par defaut des variables utilisees dans le script + self.default = { + "NiveauImpression" : "1", + "FichierDataIn" : "NON", + "FichierTempSigma" : "NON", + "FichierResultats" : "NON", + "FichierCSV" : "NON", + "FichierCREARE" : "NON", + "GrandeurEvaluee" : "FM_KICSURKCP", + "IncrementTemporel" : "1", + "IncrementMaxTemperature" : "0.1", + "IncrementMaxTempsAffichage" : "1000.", + "TraitementGeometrie" : "GEOMETRIE", + "RayonInterne" : "1.994", + "RayonInterne_mess" : "NON", + "RayonExterne" : "2.2015", + "RayonExterne_mess" : "NON", + "EpaisseurRevetement" : "0.0075", + "EpaisseurRevetement_mess" : "NON", + "LigamentExterneMin" : "0.75", + "LigamentExterneMin_mess" : "NON", + "NombreNoeudsMaillage" : "300", + "TypeInitial" : "DSR", + "Orientation" : "LONGITUD", + "ProfondeurRadiale" : "0.006", + "ProfondeurRadiale_mess" : "NON", + "ModeCalculLongueur" : "VALEUR", + "Longueur" : "0.060", + "Longueur_mess" : "NON", + "CoefDirecteur" : "10.", + "CoefDirecteur_mess" : "NON", + "Constante" : "0.", + "ModeCalculDecalage" : "VALEUR", + "DecalageNormalise" : "0.1", + "DecalageNormalise_mess" : "NON", + "DecalageRadial" : "0.", + "DecalageRadial_mess" : "NON", + "Azimut" : "0.", + "Azimut_mess" : "NON", + "Altitude" : "-4.", + "Altitude_mess" : "NON", + "Pointe" : "B", + "ModeleFluence" : "Reglementaire", + "ZoneActiveCoeur_AltitudeSup" : "-3.536", + "ZoneActiveCoeur_AltitudeInf" : "-7.194", + "FluenceMax" : "6.5", + "KPFrance" : "12.7", + "KPUS" : "9.4488", + "Azimut_0deg" : "5.8", + "Azimut_5deg" : "5.48", + "Azimut_10deg" : "4.46", + "Azimut_15deg" : "3.41", + "Azimut_20deg" : "3.37", + "Azimut_25deg" : "3.16", + "Azimut_30deg" : "2.74", + "Azimut_35deg" : "2.25", + "Azimut_40deg" : "1.89", + "Azimut_45deg" : "1.78", + "TypeIrradiation" : "RTNDT", + "RTNDT" : "64.", + "ModeleIrradiation" : "HOUSSIN", + "TeneurCuivre" : "0.0972", + "TeneurCuivre_mess" : "NON", + "TeneurNickel" : "0.72", + "TeneurNickel_mess" : "NON", + "TeneurPhosphore" : "0.00912", + "TeneurPhosphore_mess" : "NON", + "MoyenneRTndt" : "-12.0", + "MoyenneRTndt_mess" : "NON", + "CoefVariationRTndt" : "0.1", + "CoefVariationRTndt_mess" : "NON", + "EcartTypeRTndt" : "-2.", + "EcartTypeRTndt_mess" : "NON", + "NombreEcartTypeRTndt" : "2.", + "NombreEcartTypeRTndt_mess" : "NON", + "ModeleTenacite" : "RCC-M", + "NBRE_CARACTERISTIQUE" : "QUANTILE", + "NbEcartType_MoyKIc" : "-2.", + "NbEcartType_MoyKIc_mess" : "NON", + "PalierDuctile_KIc" : "195.", + "CoefficientVariation_KIc" : "0.15", + "Fractile_KIc" : "5.", + "Fractile_KIc_mess" : "NON", + "Temperature_KIc100" : "-27.", + "A1" : "21.263", + "A2" : "9.159", + "A3" : "0.04057", + "B1" : "17.153", + "B2" : "55.089", + "B3" : "0.0144", + "C1" : "4.", + "C2" : "0.", + "C3" : "0.", + "AttnCorrBeta" : "NON", + "CorrIrwin" : "NON", + "ArretDeFissure" : "NON", + "IncrementTailleFissure" : "0.", + "IncrementTailleFissure_mess" : "NON", + "NbEcartType_MoyKIa" : "0.", + "PalierDuctile_KIa" : "0.", + "CoefficientVariation_KIa" : "0.", + "InstantInitialisation" : "-1.", + "ConditionLimiteThermiqueREV" : "CHALEUR", + "TemperatureDeformationNulleREV" : "20.", + "TemperaturePourCoefDilatThermREV" : "287.", + "CoefficientPoissonREV" : "0.3", + "ConditionLimiteThermiqueMDB" : "CHALEUR", + "TemperatureDeformationNulleMDB" : "20.", + "TemperaturePourCoefDilatThermMDB" : "287.", + "CoefficientPoissonMDB" : "0.3", + "TypeConditionLimiteThermique" : "TEMP_IMPO", + "Instant_1" : "21.", + "Instant_2" : "45.", + "DebitAccumule" : "2.3", + "DebitInjectionSecurite" : "0.375", + "TempInjectionSecurite_mess" : "NON", + "DiametreHydraulique" : "0.3816", + "DiametreHydraulique_mess" : "NON", + "SectionEspaceAnnulaire" : "0.21712", + "SectionEspaceAnnulaire_mess" : "NON", + "HauteurCaracConvectionNaturelle" : "6.", + "HauteurCaracConvectionNaturelle_mess" : "NON", + "CritereConvergenceRelative" : "0.00001", + "CoefficientsVestale" : "NON", + "VolumeMelange_CREARE" : "14.9", + "VolumeMelange_CREARE_mess" : "NON", + "TemperatureInitiale_CREARE" : "250.", + "TemperatureInitiale_CREARE_mess" : "NON", + "SurfaceEchange_FluideStructure" : "0.", + "SurfaceEchange_FluideStructure_mess" : "NON", + } + + # Ce dictionnaire liste la rubrique d'appartenance des variables utilisees dans le script + self.bloc = { + "NiveauImpression" : "OPTIONS", + "FichierDataIn" : "OPTIONS", + "FichierTempSigma" : "OPTIONS", + "FichierResultats" : "OPTIONS", + "FichierCSV" : "OPTIONS", + "FichierCREARE" : "OPTIONS", + "GrandeurEvaluee" : "OPTIONS", + "IncrementTemporel" : "OPTIONS", + "IncrementMaxTemperature" : "OPTIONS", + "IncrementMaxTempsAffichage" : "OPTIONS", + "TraitementGeometrie" : "DONNEES DE LA CUVE", + "RayonInterne" : "DONNEES DE LA CUVE", + "RayonInterne_mess" : "DONNEES DE LA CUVE", + "RayonExterne" : "DONNEES DE LA CUVE", + "RayonExterne_mess" : "DONNEES DE LA CUVE", + "EpaisseurRevetement" : "DONNEES DE LA CUVE", + "EpaisseurRevetement_mess" : "DONNEES DE LA CUVE", + "LigamentExterneMin" : "DONNEES DE LA CUVE", + "LigamentExterneMin_mess" : "DONNEES DE LA CUVE", + "NombreNoeudsMaillage" : "DONNEES DE LA CUVE", + "TypeInitial" : "CARACTERISTIQUES DU DEFAUT", + "Orientation" : "CARACTERISTIQUES DU DEFAUT", + "ProfondeurRadiale" : "CARACTERISTIQUES DU DEFAUT", + "ProfondeurRadiale_mess" : "CARACTERISTIQUES DU DEFAUT", + "ModeCalculLongueur" : "CARACTERISTIQUES DU DEFAUT", + "Longueur" : "CARACTERISTIQUES DU DEFAUT", + "Longueur_mess" : "CARACTERISTIQUES DU DEFAUT", + "CoefDirecteur" : "CARACTERISTIQUES DU DEFAUT", + "CoefDirecteur_mess" : "CARACTERISTIQUES DU DEFAUT", + "Constante" : "CARACTERISTIQUES DU DEFAUT", + "ModeCalculDecalage" : "CARACTERISTIQUES DU DEFAUT", + "DecalageNormalise" : "CARACTERISTIQUES DU DEFAUT", + "DecalageNormalise_mess" : "CARACTERISTIQUES DU DEFAUT", + "DecalageRadial" : "CARACTERISTIQUES DU DEFAUT", + "DecalageRadial_mess" : "CARACTERISTIQUES DU DEFAUT", + "Azimut" : "CARACTERISTIQUES DU DEFAUT", + "Azimut_mess" : "CARACTERISTIQUES DU DEFAUT", + "Altitude" : "CARACTERISTIQUES DU DEFAUT", + "Altitude_mess" : "CARACTERISTIQUES DU DEFAUT", + "Pointe" : "CARACTERISTIQUES DU DEFAUT", + "ModeleFluence" : "MODELES", + "ZoneActiveCoeur_AltitudeSup" : "MODELES", + "ZoneActiveCoeur_AltitudeInf" : "MODELES", + "FluenceMax" : "MODELES", + "KPFrance" : "MODELES", + "KPUS" : "MODELES", + "Azimut_0deg" : "MODELES", + "Azimut_5deg" : "MODELES", + "Azimut_10deg" : "MODELES", + "Azimut_15deg" : "MODELES", + "Azimut_20deg" : "MODELES", + "Azimut_25deg" : "MODELES", + "Azimut_30deg" : "MODELES", + "Azimut_35deg" : "MODELES", + "Azimut_40deg" : "MODELES", + "Azimut_45deg" : "MODELES", + "TypeIrradiation" : "MODELES", + "RTNDT" : "MODELES", + "ModeleIrradiation" : "MODELES", + "TeneurCuivre" : "MODELES", + "TeneurCuivre_mess" : "MODELES", + "TeneurNickel" : "MODELES", + "TeneurNickel_mess" : "MODELES", + "TeneurPhosphore" : "MODELES", + "TeneurPhosphore_mess" : "MODELES", + "MoyenneRTndt" : "MODELES", + "MoyenneRTndt_mess" : "MODELES", + "CoefVariationRTndt" : "MODELES", + "CoefVariationRTndt_mess" : "MODELES", + "EcartTypeRTndt" : "MODELES", + "EcartTypeRTndt_mess" : "MODELES", + "NombreEcartTypeRTndt" : "MODELES", + "NombreEcartTypeRTndt_mess" : "MODELES", + "ModeleTenacite" : "MODELES", + "NBRE_CARACTERISTIQUE" : "MODELES", + "NbEcartType_MoyKIc" : "MODELES", + "NbEcartType_MoyKIc_mess" : "MODELES", + "PalierDuctile_KIc" : "MODELES", + "CoefficientVariation_KIc" : "MODELES", + "Fractile_KIc" : "MODELES", + "Fractile_KIc_mess" : "MODELES", + "Temperature_KIc100" : "MODELES", + "A1" : "MODELES", + "A2" : "MODELES", + "A3" : "MODELES", + "B1" : "MODELES", + "B2" : "MODELES", + "B3" : "MODELES", + "C1" : "MODELES", + "C2" : "MODELES", + "C3" : "MODELES", + "AttnCorrBeta" : "MODELES", + "CorrIrwin" : "MODELES", + "ArretDeFissure" : "MODELES", + "IncrementTailleFissure" : "MODELES", + "IncrementTailleFissure_mess" : "MODELES", + "NbEcartType_MoyKIa" : "MODELES", + "PalierDuctile_KIa" : "MODELES", + "CoefficientVariation_KIa" : "MODELES", + "InstantInitialisation" : "ETAT INITIAL", + "ConditionLimiteThermiqueREV" : "CARACTERISTIQUES DU REVETEMENT", + "TemperatureDeformationNulleREV" : "CARACTERISTIQUES DU REVETEMENT", + "TemperaturePourCoefDilatThermREV" : "CARACTERISTIQUES DU REVETEMENT", + "CoefficientPoissonREV" : "CARACTERISTIQUES DU REVETEMENT", + "ConditionLimiteThermiqueMDB" : "CARACTERISTIQUES DU MDB", + "TemperatureDeformationNulleMDB" : "CARACTERISTIQUES DU MDB", + "TemperaturePourCoefDilatThermMDB" : "CARACTERISTIQUES DU MDB", + "CoefficientPoissonMDB" : "CARACTERISTIQUES DU MDB", + "TypeConditionLimiteThermique" : "TRANSITOIRE", + "Instant_1" : "TRANSITOIRE", + "Instant_2" : "TRANSITOIRE", + "DebitAccumule" : "TRANSITOIRE", + "DebitInjectionSecurite" : "TRANSITOIRE", + "TempInjectionSecurite_mess" : "TRANSITOIRE", + "DiametreHydraulique" : "TRANSITOIRE", + "DiametreHydraulique_mess" : "TRANSITOIRE", + "SectionEspaceAnnulaire" : "TRANSITOIRE", + "SectionEspaceAnnulaire_mess" : "TRANSITOIRE", + "HauteurCaracConvectionNaturelle" : "TRANSITOIRE", + "HauteurCaracConvectionNaturelle_mess" : "TRANSITOIRE", + "CritereConvergenceRelative" : "TRANSITOIRE", + "CoefficientsVestale" : "TRANSITOIRE", + "VolumeMelange_CREARE" : "TRANSITOIRE", + "VolumeMelange_CREARE_mess" : "TRANSITOIRE", + "TemperatureInitiale_CREARE" : "TRANSITOIRE", + "TemperatureInitiale_CREARE_mess" : "TRANSITOIRE", + "SurfaceEchange_FluideStructure" : "TRANSITOIRE", + "SurfaceEchange_FluideStructure_mess" : "TRANSITOIRE", + } + + def gener(self,obj,format='brut'): + self.text='' + self.textCuve='' + self.dico_mot={} + self.dico_genea={} + self.text=PythonGenerator.gener(self,obj,format) + return self.text + + def generMCSIMP(self,obj) : + self.dico_mot[obj.nom]=obj.valeur + clef="" + for i in obj.get_genealogie() : + clef=clef+"_"+i + self.dico_genea[clef]=obj.valeur + s=PythonGenerator.generMCSIMP(self,obj) + return s + + def writeCuve2DG(self): + print "je passe dans writeCuve2DG" + self.genereTexteCuve() + f = open( "/tmp/data_template", 'wb') + print self.texteCuve + f.write( self.texteCuve ) + f.close() + + def entete(self): + ''' + Ecrit l'entete du fichier data_template + ''' + texte = "############################################################################################"+"\n" + texte += "#"+"\n" + texte += "# OUTIL D'ANALYSE PROBABILISTE DE LA DUREE DE VIE DES CUVES REP"+"\n" + texte += "# ---------------"+"\n" + texte += "# FICHIER DE MISE EN DONNEES"+"\n" + texte += "#"+"\n" + texte += "# SI CALCUL DETERMINISTE :"+"\n" + texte += "# - fixer INCRTPS=1, nbectDRTNDT=2., nbectKIc=-2."+"\n" + texte += "# - les calculs ne sont possibles qu'en une seule pointe du defaut (POINDEF<>BOTH)"+"\n" + texte += "# SI CALCUL PROBABILISTE :"+"\n" + texte += "# - fixer ARRETFISSURE=NON"+"\n" + texte += "#"+"\n" + texte += "############################################################################################"+"\n" + texte += "#"+"\n" + return texte + + def rubrique(self, titre): + ''' + Rubrique + ''' + texte = "#"+"\n" + texte += "############################################################################################"+"\n" + texte += "# " + titre + "\n" + texte += "############################################################################################"+"\n" + texte += "#"+"\n" + return texte + + def sousRubrique(self, soustitre, numtitre): + ''' + Sous-rubrique + ''' + texte = "#"+"\n" + texte += "# " + numtitre + soustitre + "\n" + texte += "#==========================================================================================="+"\n" + texte += "#"+"\n" + return texte + + def ecritLigne(self, variablelue): + ''' + Ecrit l'affectation d'une valeur a sa variable, suivie d'un commentaire + ''' + texte = "%s = %s # %s\n" % (self.variable[variablelue], str(self.dico_mot[variablelue]), self.comment[variablelue]) + return texte + + def affecteValeurDefaut(self, variablelue): + ''' + Affecte une valeur par defaut a une variable, suivie d'un commentaire + ''' + print "Warning ==> Dans la rubrique",self.bloc[variablelue],", valeur par defaut pour ",variablelue," = ",self.default[variablelue] + texte = "%s = %s # %s\n" % (self.variable[variablelue], self.default[variablelue], self.comment[variablelue]) + return texte + + def affecteValeur(self, variablelue, valeuraffectee): + ''' + Affecte une valeur a une variable, suivie d'un commentaire + ''' + texte = "%s = %s # %s\n" % (self.variable[variablelue], valeuraffectee, self.comment[variablelue]) + return texte + + def ecritVariable(self, variablelue): + if self.dico_mot.has_key(variablelue): + texte = self.ecritLigne(variablelue) + else : + texte = self.affecteValeurDefaut(variablelue) + return texte + + def amontAval(self, amont, aval): + if str(self.dico_mot[amont])=='Continu': + if str(self.dico_mot[aval])=='Continu': + texte = 'CC'+"\n" + if str(self.dico_mot[aval])=='Lineaire': + texte = 'CL'+"\n" + if str(self.dico_mot[aval])=='Exclu': + texte = 'CE'+"\n" + if str(self.dico_mot[amont])=='Lineaire': + if str(self.dico_mot[aval])=='Continu': + texte = 'LC'+"\n" + if str(self.dico_mot[aval])=='Lineaire': + texte = 'LL'+"\n" + if str(self.dico_mot[aval])=='Exclu': + texte = 'LE'+"\n" + if str(self.dico_mot[amont])=='Exclu': + if str(self.dico_mot[aval])=='Continu': + texte = 'EC'+"\n" + if str(self.dico_mot[aval])=='Lineaire': + texte = 'EL'+"\n" + if str(self.dico_mot[aval])=='Exclu': + texte = 'EE'+"\n" + return texte + + def genereTexteCuve(self): + self.texteCuve = "" + self.texteCuve += self.entete() + + # Rubrique OPTIONS + self.texteCuve += self.rubrique('OPTIONS') + + self.texteCuve += self.sousRubrique('Impression a l ecran', '') + if self.dico_mot.has_key('NiveauImpression'): + if str(self.dico_mot["NiveauImpression"])=='Aucune impression': + self.texteCuve += self.affecteValeur('NiveauImpression', '0') + if str(self.dico_mot["NiveauImpression"])=='Temps total': + self.texteCuve += self.affecteValeur('NiveauImpression', '1') + if str(self.dico_mot["NiveauImpression"])=='Temps intermediaires': + self.texteCuve += self.affecteValeur('NiveauImpression', '2') + + self.texteCuve += self.sousRubrique('Generation de fichiers', '') + self.texteCuve += self.ecritVariable('FichierDataIn') + self.texteCuve += self.ecritVariable('FichierTempSigma') + self.texteCuve += self.ecritVariable('FichierResultats') + self.texteCuve += self.ecritVariable('FichierCSV') + self.texteCuve += self.ecritVariable('FichierCREARE') + + self.texteCuve += self.sousRubrique('Grandeur evaluee', '') + if self.dico_mot.has_key('GrandeurEvaluee'): + if str(self.dico_mot["GrandeurEvaluee"])=='Facteur de marge KIc/KCP': + self.texteCuve += self.affecteValeur('GrandeurEvaluee', 'FM_KICSURKCP') + if str(self.dico_mot["GrandeurEvaluee"])=='Marge KIc-KI': + self.texteCuve += self.affecteValeur('GrandeurEvaluee', 'MARGE_KI') + if str(self.dico_mot["GrandeurEvaluee"])=='Marge KIc-KCP': + self.texteCuve += self.affecteValeur('GrandeurEvaluee', 'MARGE_KCP') + + self.texteCuve += self.sousRubrique('Divers', '') + self.texteCuve += self.ecritVariable('IncrementTemporel') + self.texteCuve += self.ecritVariable('IncrementMaxTemperature') + self.texteCuve += self.ecritVariable('IncrementMaxTempsAffichage') + if self.dico_mot.has_key('ListeInstants'): + self.texteCuve += "# liste des instants pour ecriture des resultats (s)"+"\n" + self.imprime(1,(self.dico_mot["ListeInstants"])) + else : + print "Warning ==> Dans la rubrique OPTIONS, fournir ListeInstants." + self.texteCuve += "# liste des instants pour ecriture des resultats (s)"+"\n" + self.texteCuve += " 0.\n" + self.texteCuve += " 1.\n" + + + # Rubrique DONNEES DE LA CUVE + self.texteCuve += self.rubrique('DONNEES DE LA CUVE') + if self.dico_mot.has_key('TraitementGeometrie'): + if str(self.dico_mot["TraitementGeometrie"])=='Topologie': + self.texteCuve += self.affecteValeur('TraitementGeometrie', 'GEOMETRIE') + self.texteCuve+="# - si MAILLAGE, fournir NBNO et liste des abscisses (m)"+"\n" + self.texteCuve+="# - si GEOMETRIE, fournir (RINT, RINT_MESSAGE),"+"\n" + self.texteCuve+="# (REXT, REXT_MESSAGE),"+"\n" + self.texteCuve+="# (LREV, LREV_MESSAGE),"+"\n" + self.texteCuve+="# (LIGMIN, LIGMIN_MESSAGE),"+"\n" + self.texteCuve+="# NBNO"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('RayonInterne') + self.texteCuve += self.ecritVariable('RayonInterne_mess') + self.texteCuve += self.ecritVariable('RayonExterne') + self.texteCuve += self.ecritVariable('RayonExterne_mess') + self.texteCuve += self.ecritVariable('EpaisseurRevetement') + self.texteCuve += self.ecritVariable('EpaisseurRevetement_mess') + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('LigamentExterneMin') + self.texteCuve += self.ecritVariable('LigamentExterneMin_mess') + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('NombreNoeudsMaillage') + if str(self.dico_mot["TraitementGeometrie"])=='Maillage': + self.texteCuve += self.affecteValeur('TraitementGeometrie', 'MAILLAGE') + self.texteCuve+="# - si MAILLAGE, fournir NBNO et liste des abscisses (m)"+"\n" + self.texteCuve+="# - si GEOMETRIE, fournir (RINT, RINT_MESSAGE),"+"\n" + self.texteCuve+="# (REXT, REXT_MESSAGE),"+"\n" + self.texteCuve+="# (LREV, LREV_MESSAGE),"+"\n" + self.texteCuve+="# (LIGMIN, LIGMIN_MESSAGE),"+"\n" + self.texteCuve+="# NBNO"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('NombreNoeudsMaillage') + self.imprime(1,(self.dico_mot["ListeAbscisses"])) + else : + self.texteCuve += self.affecteValeurDefaut('TraitementGeometrie') + self.texteCuve+="# - si MAILLAGE, fournir NBNO et liste des abscisses (m)"+"\n" + self.texteCuve+="# - si GEOMETRIE, fournir (RINT, RINT_MESSAGE),"+"\n" + self.texteCuve+="# (REXT, REXT_MESSAGE),"+"\n" + self.texteCuve+="# (LREV, LREV_MESSAGE),"+"\n" + self.texteCuve+="# (LIGMIN, LIGMIN_MESSAGE),"+"\n" + self.texteCuve+="# NBNO"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve += self.affecteValeurDefaut('RayonInterne') + self.texteCuve += self.affecteValeurDefaut('RayonInterne_mess') + self.texteCuve += self.affecteValeurDefaut('RayonExterne') + self.texteCuve += self.affecteValeurDefaut('RayonExterne_mess') + self.texteCuve += self.affecteValeurDefaut('EpaisseurRevetement') + self.texteCuve += self.affecteValeurDefaut('EpaisseurRevetement_mess') + self.texteCuve+="#"+"\n" + self.texteCuve += self.affecteValeurDefaut('LigamentExterneMin') + self.texteCuve += self.affecteValeurDefaut('LigamentExterneMin_mess') + self.texteCuve+="#"+"\n" + self.texteCuve += self.affecteValeurDefaut('NombreNoeudsMaillage') + + + # Rubrique CARACTERISTIQUES DU DEFAUT + self.texteCuve += self.rubrique('CARACTERISTIQUES DU DEFAUT') + + if self.dico_mot.has_key('TypeInitial'): + if str(self.dico_mot["TypeInitial"])=='Defaut Sous Revetement': + self.texteCuve += self.affecteValeur('TypeInitial', 'DSR') + if str(self.dico_mot["TypeInitial"])=='Defaut Decale': + self.texteCuve += self.affecteValeur('TypeInitial', 'DECALE') + if str(self.dico_mot["TypeInitial"])=='Defaut Debouchant': + self.texteCuve += self.affecteValeur('TypeInitial', 'DEBOUCHANT') + else : + self.texteCuve += self.affecteValeurDefaut('TypeInitial') + + self.texteCuve+="# Fournir ORIEDEF, (PROFDEF, PROFDEF_MESSAGE)"+"\n" + self.texteCuve+="# - Si DSR, fournir OPTLONG, (LONGDEF,LONGDEF_MESSAGE) ou (PROFSURLONG,PROFSURLONG_MESSAGE,LONGCONST)"+"\n" + self.texteCuve+="# - Si DECALE, fournir OPTLONG, (LONGDEF,LONGDEF_MESSAGE) ou (PROFSURLONG,PROFSURLONG_MESSAGE,LONGCONST), DECATYP, (DECANOR,DECANOR_MESSAGE) ou (DECADEF,DECADEF_MESSAGE)"+"\n" + self.texteCuve+="# - Si DEBOUCHANT, fournir IRWIN"+"\n" + self.texteCuve+="# Fournir (ANGLDEF, ANGLDEF_MESSAGE), (ALTIDEF, ALTIDEF_MESSAGE)"+"\n" + self.texteCuve+="# - Si DSR ou DECALE, fournir POINDEF"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve+="# Remarque :"+"\n" + self.texteCuve+="# - si DSR ou DECALE, dans la rubrique 'Modele de tenacite', fournir ATTNCORRBETA (ne pas fournir CORRIRWIN)"+"\n" + self.texteCuve+="# - si DEBOUCHANT, dans la rubrique 'Modele de tenacite', fournir CORRIRWIN (ne pas fournir ATTNCORRBETA)"+"\n" + + self.texteCuve+="#"+"\n" + + if self.dico_mot.has_key('Orientation'): + if str(self.dico_mot["Orientation"])=='Longitudinale': + self.texteCuve += self.affecteValeur('Orientation', 'LONGITUD') + if str(self.dico_mot["Orientation"])=='Circonferentielle': + self.texteCuve += self.affecteValeur('Orientation', 'CIRCONF') + else : + self.texteCuve += self.affecteValeurDefaut('Orientation') + + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('ProfondeurRadiale') + self.texteCuve += self.ecritVariable('ProfondeurRadiale_mess') + + self.texteCuve+="#"+"\n" + if self.dico_mot.has_key('ModeCalculLongueur'): + if str(self.dico_mot["ModeCalculLongueur"])=='Valeur': + self.texteCuve += self.affecteValeur('ModeCalculLongueur', 'VALEUR') + self.texteCuve+="# - Si VALEUR, fournir (LONGDEF, LONGDEF_MESSAGE)"+"\n" + self.texteCuve+="# - Si FCTAFFINE, fournir (PROFSURLONG, PROFSURLONG_MESSAGE) et LONGCONST : LONGDEF=PROFDEF/PROFSURLONG + LONGCONST"+"\n" + self.texteCuve += self.ecritVariable('Longueur') + self.texteCuve += self.ecritVariable('Longueur_mess') + if str(self.dico_mot["ModeCalculLongueur"])=='Fonction affine de la profondeur': + self.texteCuve += self.affecteValeur('ModeCalculLongueur', 'FCTAFFINE') + self.texteCuve+="# - Si VALEUR, fournir (LONGDEF, LONGDEF_MESSAGE)"+"\n" + self.texteCuve+="# - Si FCTAFFINE, fournir (PROFSURLONG, PROFSURLONG_MESSAGE) et LONGCONST : LONGDEF=PROFDEF/PROFSURLONG + LONGCONST"+"\n" + self.texteCuve += self.ecritVariable('CoefDirecteur') + self.texteCuve += self.ecritVariable('CoefDirecteur_mess') + self.texteCuve += self.ecritVariable('Constante') + else : + self.texteCuve += self.affecteValeurDefaut('ModeCalculLongueur') + self.texteCuve+="# - Si VALEUR, fournir (LONGDEF, LONGDEF_MESSAGE)"+"\n" + self.texteCuve+="# - Si FCTAFFINE, fournir (PROFSURLONG, PROFSURLONG_MESSAGE) et LONGCONST : LONGDEF=PROFDEF/PROFSURLONG + LONGCONST"+"\n" + self.texteCuve += self.affecteValeurDefaut('Longueur') + self.texteCuve += self.affecteValeurDefaut('Longueur_mess') + + if self.dico_mot.has_key('TypeInitial'): + if str(self.dico_mot["TypeInitial"])!='Defaut Sous Revetement': + self.texteCuve+="#"+"\n" + if self.dico_mot.has_key('ModeCalculDecalage'): + if str(self.dico_mot["ModeCalculDecalage"])=='Valeur normalisee': + self.texteCuve += self.affecteValeur('ModeCalculDecalage', 'NORMALISE') + self.texteCuve+="# - Si NORMALISE, fournir (DECANOR, DECANOR_MESSAGE)"+"\n" + self.texteCuve+="# - Si VALEUR, fournir (DECADEF, DECADEF_MESSAGE)"+"\n" + self.texteCuve += self.ecritVariable('DecalageNormalise') + self.texteCuve += self.ecritVariable('DecalageNormalise_mess') + if str(self.dico_mot["ModeCalculDecalage"])=='Valeur': + self.texteCuve += self.affecteValeur('ModeCalculDecalage', 'VALEUR') + self.texteCuve+="# - Si NORMALISE, fournir (DECANOR, DECANOR_MESSAGE)"+"\n" + self.texteCuve+="# - Si VALEUR, fournir (DECADEF, DECADEF_MESSAGE)"+"\n" + self.texteCuve += self.ecritVariable('DecalageRadial') + self.texteCuve += self.ecritVariable('DecalageRadial_mess') + else : + self.texteCuve += self.affecteValeurDefaut('ModeCalculDecalage') + self.texteCuve+="# - Si NORMALISE, fournir (DECANOR, DECANOR_MESSAGE)"+"\n" + self.texteCuve+="# - Si VALEUR, fournir (DECADEF, DECADEF_MESSAGE)"+"\n" + self.texteCuve += self.affecteValeurDefaut('DecalageRadial') + self.texteCuve += self.affecteValeurDefaut('DecalageRadial_mess') + + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('Azimut') + self.texteCuve += self.ecritVariable('Azimut_mess') + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('Altitude') + self.texteCuve += self.ecritVariable('Altitude_mess') + self.texteCuve+="#"+"\n" + if self.dico_mot.has_key('Pointe'): + if str(self.dico_mot["Pointe"])=='A': + self.texteCuve += self.affecteValeur('Pointe', 'A') + if str(self.dico_mot["Pointe"])=='B': + self.texteCuve += self.affecteValeur('Pointe', 'B') + if str(self.dico_mot["Pointe"])=='A et B': + self.texteCuve += self.affecteValeur('Pointe', 'BOTH') + else : + self.texteCuve += self.affecteValeurDefaut('Pointe') + + # Rubrique MODELES FLUENCE, IRRADIATION, TENACITE + self.texteCuve += self.rubrique('MODELES FLUENCE, IRRADIATION, TENACITE') + self.texteCuve += self.sousRubrique('Modele d attenuation de la fluence dans l epaisseur','A.') + + if self.dico_mot.has_key('ModeleFluence'): + if str(self.dico_mot["ModeleFluence"])=='Exponentiel sans revetement k=9.7 (Reglementaire)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'Reglementaire') + if str(self.dico_mot["ModeleFluence"])=='Exponentiel sans revetement k=12.7 (France)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'France') + if str(self.dico_mot["ModeleFluence"])=='Exponentiel sans revetement k=0. (ValeurImposee)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'ValeurImposee') + if str(self.dico_mot["ModeleFluence"])=='Donnees francaises du palier CPY (SDM)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'SDM') + if str(self.dico_mot["ModeleFluence"])=='Regulatory Guide 1.99 rev 2 (USNRC)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'USNRC') + if str(self.dico_mot["ModeleFluence"])=='Dossier 900 MWe AP9701 rev 2 (REV_2)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'REV_2') + if str(self.dico_mot["ModeleFluence"])=='Lissage du modele ajuste (SDM_Lissage)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'SDM_Lissage') + if str(self.dico_mot["ModeleFluence"])=='Donnees francaises du palier CPY ajustees par secteur angulaire (GrandeDev)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'GrandeDev') + if str(self.dico_mot["ModeleFluence"])=='Grand developpement (GD_Cuve)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'GD_Cuve') + if str(self.dico_mot["ModeleFluence"])=='Exponentiel sans revetement k=9.7 (Reglementaire CUVE1D)': + self.texteCuve += self.affecteValeur('ModeleFluence', 'Cuve1D') + else : + self.texteCuve += self.affecteValeurDefaut('ModeleFluence') + + self.texteCuve+="# - si France, fournir KPFRANCE"+"\n" + self.texteCuve+="# - si USNRC, fournir KPUS"+"\n" + self.texteCuve+="# - si modele GD_Cuve, fournir COEFFLUENCE1, COEFFLUENCE2, ..., COEFFLUENCE9, COEFFLUENCE10"+"\n" + self.texteCuve+="#"+"\n" + + self.texteCuve += self.ecritVariable('ZoneActiveCoeur_AltitudeSup') + self.texteCuve += self.ecritVariable('ZoneActiveCoeur_AltitudeInf') + self.texteCuve += self.ecritVariable('FluenceMax') + if self.dico_mot.has_key('ModeleFluence'): + if str(self.dico_mot["ModeleFluence"])=='Exponentiel sans revetement k=12.7 (France)': + self.texteCuve += self.ecritVariable('KPFrance') + if str(self.dico_mot["ModeleFluence"])=='Regulatory Guide 1.99 rev 2 (USNRC)': + self.texteCuve += self.ecritVariable('KPUS') + if str(self.dico_mot["ModeleFluence"])=='Grand developpement (GD_Cuve)': + self.texteCuve += self.ecritVariable('Azimut_0deg') + self.texteCuve += self.ecritVariable('Azimut_5deg') + self.texteCuve += self.ecritVariable('Azimut_10deg') + self.texteCuve += self.ecritVariable('Azimut_15deg') + self.texteCuve += self.ecritVariable('Azimut_20deg') + self.texteCuve += self.ecritVariable('Azimut_25deg') + self.texteCuve += self.ecritVariable('Azimut_30deg') + self.texteCuve += self.ecritVariable('Azimut_35deg') + self.texteCuve += self.ecritVariable('Azimut_40deg') + self.texteCuve += self.ecritVariable('Azimut_45deg') + + self.texteCuve += self.sousRubrique('Irradiation','B.') + + if self.dico_mot.has_key('TypeIrradiation'): + + if str(self.dico_mot["TypeIrradiation"])=='RTndt de la cuve a l instant de l analyse': + self.texteCuve += self.affecteValeur('TypeIrradiation', 'RTNDT') + self.texteCuve+="# - si RTNDT, fournir RTNDT"+"\n" + self.texteCuve+="# - si FLUENCE, fournir MODELIRR, et autres parametres selon MODELIRR (voir ci-dessous)"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('RTNDT') + + if str(self.dico_mot["TypeIrradiation"])=='Modele d irradiation': + self.texteCuve += self.affecteValeur('TypeIrradiation', 'FLUENCE') + self.texteCuve+="# - si RTNDT, fournir RTNDT"+"\n" + self.texteCuve+="# - si FLUENCE, fournir MODELIRR, et autres parametres selon MODELIRR (voir ci-dessous)"+"\n" + self.texteCuve+="#"+"\n" + if self.dico_mot.has_key('ModeleIrradiation'): + if str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : formule de FIM/FIS Houssin': + self.texteCuve += self.affecteValeur('ModeleIrradiation', 'HOUSSIN') + if str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : formule de FIM/FIS Persoz': + self.texteCuve += self.affecteValeur('ModeleIrradiation', 'PERSOZ') + if str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : formule de FIM/FIS Lefebvre': + self.texteCuve += self.affecteValeur('ModeleIrradiation', 'LEFEBVRE') + if str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : Regulatory Guide 1.00 rev 2': + self.texteCuve += self.affecteValeur('ModeleIrradiation', 'USNRCmdb') + if str(self.dico_mot["ModeleIrradiation"])=='Joint Soude : formulation de FIM/FIS Brillaud': + self.texteCuve += self.affecteValeur('ModeleIrradiation', 'BRILLAUD') + if str(self.dico_mot["ModeleIrradiation"])=='Joint Soude : Regulatory Guide 1.00 rev 2': + self.texteCuve += self.affecteValeur('ModeleIrradiation', 'USNRCsoud') + else : + self.texteCuve += self.affecteValeurDefaut('ModeleIrradiation') + self.texteCuve+="# - pour tout modele, fournir (CU, CU_MESSAGE),"+"\n" + self.texteCuve+="# (NI, NI_MESSAGE),"+"\n" + self.texteCuve+="# - si HOUSSIN, PERSOZ, LEFEBVRE, BRILLAUD, fournir (P, P_MESSAGE)"+"\n" + self.texteCuve+="# - pour tout modele, fournir (RTimoy, RTimoy_MESSAGE),"+"\n" + self.texteCuve+="# - si USNRCsoud ou USNRCmdb, fournir (RTicov, RTicov_MESSAGE)"+"\n" + self.texteCuve+="# (USectDRT, USectDRT_MESSAGE)"+"\n" + self.texteCuve+="# - pour tout modele, fournir (nbectDRTNDT, nbectDRTNDT_MESSAGE)"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve += self.ecritVariable('TeneurCuivre') + self.texteCuve += self.ecritVariable('TeneurCuivre_mess') + self.texteCuve += self.ecritVariable('TeneurNickel') + self.texteCuve += self.ecritVariable('TeneurNickel_mess') + if str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : formule de FIM/FIS Houssin' or str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : formule de FIM/FIS Persoz' or str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : formule de FIM/FIS Lefebvre' or str(self.dico_mot["ModeleIrradiation"])=='Joint Soude : formulation de FIM/FIS Brillaud': + self.texteCuve += self.ecritVariable('TeneurPhosphore') + self.texteCuve += self.ecritVariable('TeneurPhosphore_mess') + self.texteCuve += self.ecritVariable('MoyenneRTndt') + self.texteCuve += self.ecritVariable('MoyenneRTndt_mess') + if str(self.dico_mot["ModeleIrradiation"])=='Metal de Base : Regulatory Guide 1.00 rev 2' or str(self.dico_mot["ModeleIrradiation"])=='Joint Soude : Regulatory Guide 1.00 rev 2': + self.texteCuve += self.ecritVariable('CoefVariationRTndt') + self.texteCuve += self.ecritVariable('CoefVariationRTndt_mess') + self.texteCuve += self.ecritVariable('EcartTypeRTndt') + self.texteCuve += self.ecritVariable('EcartTypeRTndt_mess') + self.texteCuve += self.ecritVariable('NombreEcartTypeRTndt') + self.texteCuve += self.ecritVariable('NombreEcartTypeRTndt_mess') + else : + self.texteCuve += self.affecteValeurDefaut('TypeIrradiation') + self.texteCuve+="# - si RTNDT, fournir RTNDT"+"\n" + self.texteCuve+="# - si FLUENCE, fournir MODELIRR, et autres parametres selon MODELIRR (voir ci-dessous)"+"\n" + self.texteCuve+="#"+"\n" + self.texteCuve += self.affecteValeurDefaut('RTNDT') + + self.texteCuve += self.sousRubrique('Modele de tenacite','C.') + self.texteCuve+="# tenacite d amorcage"+"\n" + + if self.dico_mot.has_key('ModeleTenacite'): + if str(self.dico_mot["ModeleTenacite"])=='RCC-M/ASME coefficient=2': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'RCC-M') + if str(self.dico_mot["ModeleTenacite"])=='RCC-M/ASME coefficient=2 CUVE1D': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'RCC-M_simpl') + if str(self.dico_mot["ModeleTenacite"])=='RCC-M/ASME coefficient=2.33 (Houssin)': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'Houssin_RC') + if str(self.dico_mot["ModeleTenacite"])=='RCC-M/ASME avec KI=KIpalier': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'RCC-M_pal') + if str(self.dico_mot["ModeleTenacite"])=='RCC-M/ASME avec KI~exponentiel': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'RCC-M_exp') + if str(self.dico_mot["ModeleTenacite"])=='Weibull basee sur la master cuve': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'Wallin') + if str(self.dico_mot["ModeleTenacite"])=='Weibull basee sur la master cuve (REME)': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'REME') + if str(self.dico_mot["ModeleTenacite"])=='Weibull n\xb01 (etude ORNL)': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'ORNL') + if str(self.dico_mot["ModeleTenacite"])=='Weibull n\xb02': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'WEIB2') + if str(self.dico_mot["ModeleTenacite"])=='Weibull n\xb03': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'WEIB3') + if str(self.dico_mot["ModeleTenacite"])=='Weibull generalisee': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'WEIB_GEN') + if str(self.dico_mot["ModeleTenacite"])=='Exponentielle n\xb01 (Frama)': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'Frama') + if str(self.dico_mot["ModeleTenacite"])=='Exponentielle n\xb02 (LOGWOLF)': + self.texteCuve += self.affecteValeur('ModeleTenacite', 'LOGWOLF') + else : + self.texteCuve += self.affecteValeurDefaut('ModeleTenacite') + self.texteCuve+="# - si RCC-M, RCC-M_pal, Houssin_RC, fournir (nbectKIc, nbectKIc_MESSAGE), KICPAL, KICCDV"+"\n" + self.texteCuve+="# - si RCC-M_exp, fournir (nbectKIc, nbectKIc_MESSAGE), KICCDV"+"\n" + self.texteCuve+="# - si RCC-M_simpl, ne rien fournir"+"\n" + self.texteCuve+="# - si Frama, LOGWOLF, fournir (nbectKIc, nbectKIc_MESSAGE)"+"\n" + self.texteCuve+="# - si REME, ORNL, WEIB3, WEIB2, fournir NBCARAC, puis (nbectKIc, nbectKIc_MESSAGE) ou (fractKIc, fractKIc_MESSAGE) selon valeur de NBCARAC"+"\n" + self.texteCuve+="# - si Wallin, fournir NBCARAC, puis (nbectKIc, nbectKIc_MESSAGE) ou (fractKIc, fractKIc_MESSAGE) selon valeur de NBCARAC,"+"\n" + self.texteCuve+="# puis T0WALLIN"+"\n" + self.texteCuve+="# - si WEIB-GEN, fournir NBCARAC, puis (nbectKIc, nbectKIc_MESSAGE) ou (fractKIc, fractKIc_MESSAGE) selon valeur de NBCARAC,"+"\n" + self.texteCuve+="# puis A1, A2, A3, B1, B2, B3, C1, C2, C3"+"\n" + self.texteCuve+="# loi de Weibull P(K