u7.03.41 Procédure MACR_ECREVISSE#
Résumé:
La macro-commande MACR_ECREVISSE a pour but de calculer le débit de fluide (air/vapeur/liquide) à travers une (ou plusieurs) fissure(s) traversante(s) dans une structure modélisée en 2 dimensions. Elle réalise pour cela, à chaque pas de temps, le chaînage de deux codes :
Code_Aster qui permet de connaître l’état thermo-mécanique de la structure
Ecrevisse qui réalise le calcul thermo-hydraulique d’écoulement à travers la fissure.
En pratique, la macro-commande se charge d’effectuer successivement pour tous les pas de temps :
un calcul thermique linéaire
un calcul mécanique quasi-statique
d’appeler la macro-commande CALC_ECREVISSE, qui elle va être en charge :
de récupérer le profil de la fissure
d’appeler une troisième macro-commande MACR_ECRE_CALC qui va générer le fichier de commande et lancer Ecrevisse
d’extraire et copier les résultats
de vérifier les critères imposés sont vérifiés et si besoin redécouper le pas de temps.
Le concept de sortie principal est la structure de données résultat du calcul mécanique (evol_noli). Il est également possible d’obtenir la structure de données résultats du calcul thermique, ainsi que 2 tables, l’une contenant les débits à chaque instant, l’autre récapitulant les données concernant les fissures. En créant un répertoire de sortie dans le profil d’étude, on pourra y récupérer toutes les sorties concernant Ecrevisse.
Il est possible d’effectuer des poursuites du calcul.
Remarque
Il est nécessaire de préciser le mot-clé DEBUG=_F(HIST_ETAPE=”OUI”) dans DEBUT ou POURSUITE pour utiliser MACR_ECREVISSE .
Syntaxe#
Détail de la syntaxe
/ evol_noli
/ evol_ther
/ table = MACR_ECREVISSE(
◇ reuse = <objet de ETAT_INIT/EVOL_NOLI>,
◇ TABLE = co,
◇ DEBIT = co,
/ TEMPER = co,
◇ / ETAT_INIT = _F(
◆ EVOL_NOLI = evol_noli,
◆ EVOL_THER = evol_ther,
◆ NUME_ORDRE = int,
),
◆ MODELE_MECA = modele,
◆ MODELE_THER = modele,
◆ FISSURE = _F(
◆ PREFIXE_FICHIER = text,
◆ GROUP_MA = grma,
◆ GROUP_NO_ORIG = grno,
◆ GROUP_NO_EXTR = grno,
◆ ZETA = float,
◆ RUGOSITE = float,
◇ TORTUOSITE = float (défaut: 1.0),
◆ OUVERT_REMANENTE = float,
◆ SECTION = / "ELLIPSE",
/ "RECTANGLE",
# Si: equal_to("SECTION", 'ELLIPSE')
◇ LISTE_COTES_BL = float,
◆ LISTE_VAL_BL = float,
# Si: equal_to("SECTION", 'RECTANGLE')
◇ LISTE_COTES_BL = float,
◆ LISTE_VAL_BL = float,
),
◆ ECOULEMENT = _F(
◆ / PRES_ENTREE = float,
/ PRES_ENTREE_FO = fonction / formule / nappe,
◆ / PRES_SORTIE = float,
/ PRES_SORTIE_FO = fonction / formule / nappe,
◆ FLUIDE_ENTREE = / 1,
/ 2,
/ 3,
/ 4,
/ 5,
/ 6,
# Si: equal_to("FLUIDE_ENTREE", 1)
◆ / TEMP_ENTREE = float,
/ TEMP_ENTREE_FO = fonction / formule / nappe,
# Si: equal_to("FLUIDE_ENTREE", 2)
◆ / TITR_MASS = float,
/ TITR_MASS_FO = fonction / formule / nappe,
# Si: equal_to("FLUIDE_ENTREE", 3)
◆ / TEMP_ENTREE = float,
/ TEMP_ENTREE_FO = fonction / formule / nappe,
# Si: equal_to("FLUIDE_ENTREE", 4)
◆ / TEMP_ENTREE = float,
/ TEMP_ENTREE_FO = fonction / formule / nappe,
◆ / PRES_PART = float,
/ PRES_PART_FO = fonction / formule / nappe,
# Si: equal_to("FLUIDE_ENTREE", 5)
◆ / TITR_MASS = float,
/ TITR_MASS_FO = fonction / formule / nappe,
◆ / PRES_PART = float,
/ PRES_PART_FO = fonction / formule / nappe,
# Si: equal_to("FLUIDE_ENTREE", 6)
◆ / TEMP_ENTREE = float,
/ TEMP_ENTREE_FO = fonction / formule / nappe,
),
◇ LIST_INST = listr8,
◆ MODELE_ECRE = _F(
◇ IVENAC = / 0 (par défaut),
/ 1,
◆ ECOULEMENT = / "GELE",
/ "SATURATION",
# Si: equal_to("ECOULEMENT", 'GELE')
◆ PRESS_EBULLITION = float,
◆ FROTTEMENT = / -1,
/ -2,
/ -3,
/ -4,
/ 0,
/ 1,
/ 11,
/ 12,
/ 13,
/ 14,
/ 2,
/ 21,
/ 22,
/ 23,
/ 24,
/ 3,
/ 4,
# Si: less_than('FROTTEMENT', 0)
◆ REYNOLDS_LIM = float,
◆ FROTTEMENT_LIM = float,
◆ TRANSFERT_CHAL = / -1,
/ -11,
/ -12,
/ -2,
/ 0,
/ 1,
/ 11,
/ 12,
/ 2,
# Si: less_than('TRANSFERT_CHAL', 0)
◆ XMINCH = float,
◆ XMAXCH = float,
),
◆ CONV_CRITERE = _F(
◆ TEMP_REF = float,
◆ PRES_REF = float,
◇ CRITERE = / "EXPLICITE",
/ "PRESS",
/ "TEMP",
/ "TEMP_PRESS" (par défaut),
# Si: equal_to("CRITERE", 'TEMP_PRESS') or equal_to("CRITERE", 'TEMP') or equal_to("CRITERE", 'PRESS')
◇ SUBD_NIVEAU = int (défaut: 3),
◇ SUBD_PAS_MINI = float,
◇ NUME_ORDRE_MIN = int (défaut: -1),
◇ PREC_CRIT = float (défaut: 1.0),
),
◇ CONVERGENCE_ECREVISSE = _F(
◇ KGTEST = float (défaut: 0.5),
◇ ITER_GLOB_MAXI = int (défaut: 400),
◇ CRIT_CONV_DEBI = float (défaut: 1e-05),
),
◇ COURBES = / "AUCUNE" (par défaut),
/ "INTERACTIF",
/ "POSTSCRIPT",
◆ / LOGICIEL = text,
/ VERSION = "3.2.2" (ou non renseigné),
◇ ENTETE = text (défaut: "Titre du calcul Ecrevisse"),
◇ IMPRESSION = / "NON" (par défaut),
/ "OUI",
◇ INFO = / 1 (par défaut),
/ 2,
◆ CHAM_MATER = cham_mater,
◆ TEMP_INIT = float,
◇ CARA_ELEM = cara_elem,
◆ EXCIT_MECA = _F(
◆ CHARGE = char_cine_meca / char_meca,
◇ FONC_MULT = fonction / formule / nappe,
◇ TYPE_CHARGE = / "DIDI",
/ "FIXE_CSTE" (par défaut),
/ "SUIV",
),
◆ CONTACT = char_contact,
⟐ COMPORTEMENT = _F(
# Si: is_in("RELATION", ('VMIS_ISOT_LINE','VMIS_ISOT_TRAC','VISCOCHAB','VISC_ISOT_LINE','VISC_ISOT_TRAC',))
◇ POST_ITER = "CRIT_RUPT" (ou non renseigné),
# Si: is_in("RELATION", ('VMIS_ISOT_LINE','VMIS_CINE_LINE','VMIS_ECMI_LINE','VMIS_ISOT_TRAC','VMIS_CIN1_CHAB','VMIS_CIN2_CHAB'))
◇ POST_INCR = / "REST_ECRO",
/ "SANS",
◆ / TOUT = "OUI" (ou non renseigné),
/ GROUP_MA = grma,
◇ RELATION = / "ARME",
/ "ASSE_CORN",
/ "BETON_BURGER",
/ "BETON_DOUBLE_DP",
/ "BETON_GRANGER",
/ "BETON_GRANGER_V",
/ "BETON_RAG",
/ "BETON_REGLE_PR",
/ "BETON_UMLV",
/ "Barcelone",
/ "CABLE",
/ "CABLE_GAINE_FROT",
/ "CAM_CLAY",
/ "CHOC_ELAS_TRAC",
/ "CHOC_ENDO",
/ "CHOC_ENDO_PENA",
/ "CJS",
/ "CORR_ACIER",
/ "CSSM",
/ "CZM_ELAS_MIX",
/ "CZM_EXP_MIX",
/ "CZM_EXP_REG",
/ "CZM_FAT_MIX",
/ "CZM_LAB_MIX",
/ "CZM_LIN_REG",
/ "CZM_OUV_MIX",
/ "CZM_TAC_MIX",
/ "CZM_TRA_MIX",
/ "CZM_TURON",
/ "DASHPOT",
/ "DHRC",
/ "DIS_BILI_ELAS",
/ "DIS_CHOC",
/ "DIS_CONTACT",
/ "DIS_ECRO_CINE",
/ "DIS_ECRO_TRAC",
/ "DIS_GOUJ2E_ELAS",
/ "DIS_GOUJ2E_PLAS",
/ "DIS_GRICRA",
/ "DIS_VISC",
/ "DRUCK_PRAGER",
/ "DRUCK_PRAG_N_A",
/ "ELAS" (par défaut),
/ "ELAS_HYPER",
/ "ELAS_HYPER_VISC",
/ "ELAS_MEMBRANE_NH",
/ "ELAS_MEMBRANE_SV",
/ "ELAS_POUTRE_GR",
/ "ELAS_VMIS_LINE",
/ "ELAS_VMIS_PUIS",
/ "ELAS_VMIS_TRAC",
/ "ENDO_CARRE",
/ "ENDO_FISS_EXP",
/ "ENDO_HETEROGENE",
/ "ENDO_ISOT_BETON",
/ "ENDO_LOCA_EXP",
/ "ENDO_LOCA_TC",
/ "ENDO_ORTH_BETON",
/ "ENDO_PORO_BETON",
/ "ENDO_SCALAIRE",
/ "FLUA_ENDO_PORO",
/ "FLUA_PORO_BETON",
/ "FONDATION",
/ "GLRC_DAMAGE",
/ "GLRC_DM",
/ "GRAN_IRRA_LOG",
/ "GRILLE_CINE_LINE",
/ "GRILLE_ISOT_LINE",
/ "GTN",
/ "GonfElas",
/ "HAYHURST",
/ "HOEK_BROWN",
/ "HOEK_BROWN_EFF",
/ "HOEK_BROWN_TOT",
/ "HUJEUX",
/ "IRRAD3M",
/ "Iwan",
/ "JOINT_BA",
/ "JOINT_BANDIS",
/ "JOINT_MECA_ENDO",
/ "JOINT_MECA_FROT",
/ "JOINT_MECA_RUPT",
/ "JONC_ENDO_PLAS",
/ "KICHENIN_NL",
/ "KIT_CG",
/ "KIT_DDI",
/ "KIT_H",
/ "KIT_HH",
/ "KIT_HH2",
/ "KIT_HH2M",
/ "KIT_HHM",
/ "KIT_HM",
/ "KIT_THH",
/ "KIT_THH2",
/ "KIT_THH2M",
/ "KIT_THHM",
/ "KIT_THM",
/ "KIT_THV",
/ "LAIGLE",
/ "LEMAITRE",
/ "LEMAITRE_IRRA",
/ "LEMA_SEUIL",
/ "LETK",
/ "LKR",
/ "MAZARS",
/ "MAZARS_UNIL",
/ "MCC",
/ "META_LEMA_ANI",
/ "META_LEMA_ANI",
/ "META_P_CL",
/ "META_P_CL_PT",
/ "META_P_CL_PT_RE",
/ "META_P_CL_RE",
/ "META_P_IL",
/ "META_P_IL_PT",
/ "META_P_IL_PT_RE",
/ "META_P_IL_RE",
/ "META_P_INL",
/ "META_P_INL_PT",
/ "META_P_INL_PT_RE",
/ "META_P_INL_RE",
/ "META_V_CL",
/ "META_V_CL_PT",
/ "META_V_CL_PT_RE",
/ "META_V_CL_RE",
/ "META_V_IL",
/ "META_V_IL_PT",
/ "META_V_IL_PT_RE",
/ "META_V_IL_RE",
/ "META_V_INL",
/ "META_V_INL_PT",
/ "META_V_INL_PT_RE",
/ "META_V_INL_RE",
/ "MFRONT",
/ "MOHR_COULOMB",
/ "MONOCRISTAL",
/ "MULTIFIBRE",
/ "MetaAcierEPIL_PT",
/ "MohrCoulombAS",
/ "NLH_CSRM",
/ "NORTON",
/ "NORTON_HOFF",
/ "POLYCRISTAL",
/ "RANKINE",
/ "RELAX_ACIER",
/ "RGI_BETON",
/ "RGI_BETON_BA",
/ "ROUSSELIER",
/ "ROUSS_PR",
/ "ROUSS_VISC",
/ "RUPT_FRAG",
/ "SANS",
/ "UMAT",
/ "VENDOCHAB",
/ "VISCOCHAB",
/ "VISC_CIN1_CHAB",
/ "VISC_CIN2_CHAB",
/ "VISC_CIN2_MEMO",
/ "VISC_CIN2_NRAD",
/ "VISC_DRUC_PRAG",
/ "VISC_ENDO_LEMA",
/ "VISC_GTN",
/ "VISC_IRRA_LOG",
/ "VISC_ISOT_LINE",
/ "VISC_ISOT_NL",
/ "VISC_ISOT_TRAC",
/ "VISC_MAXWELL",
/ "VISC_MAXWELL_MT",
/ "VISC_MEMO_NRAD",
/ "VISC_TAHERI",
/ "VMIS_ASYM_LINE",
/ "VMIS_CIN1_CHAB",
/ "VMIS_CIN2_CHAB",
/ "VMIS_CIN2_MEMO",
/ "VMIS_CIN2_NRAD",
/ "VMIS_CINE_GC",
/ "VMIS_CINE_LINE",
/ "VMIS_CINE_LINE",
/ "VMIS_ECMI_LINE",
/ "VMIS_ECMI_TRAC",
/ "VMIS_ISOT_LINE",
/ "VMIS_ISOT_NL",
/ "VMIS_ISOT_PUIS",
/ "VMIS_ISOT_TRAC",
/ "VMIS_JOHN_COOK",
/ "VMIS_MEMO_NRAD",
◇ REGU_VISC = / "NON" (par défaut),
/ "OUI",
# Si: equal_to("RELATION", 'MONOCRISTAL')
◆ COMPOR = compor,
# Si: equal_to("RELATION", 'MULTIFIBRE')
◇ RIGI_GEOM = / "DEFAUT" (par défaut),
/ "OUI",
# Si: equal_to("RELATION", 'POLYCRISTAL')
◆ COMPOR = compor,
# Si: equal_to("RELATION", 'UMAT')
◆ NB_VARI = int,
◆ LIBRAIRIE = text,
◆ NOM_ROUTINE = text,
# Si: equal_to("RELATION", 'MFRONT')
◆ COMPOR_MFRONT = compor_mgis,
◇ VERI_BORNE = / "ARRET" (par défaut),
/ "MESSAGE",
/ "SANS",
◇ ALGO_CPLAN = / "ANALYTIQUE",
/ "DEBORST" (par défaut),
◇ SYME_MATR_TANG = / "NON",
/ "OUI" (par défaut),
# Si: not equal_to('RELATION', 'MFRONT')
# Si: equal_to("RELATION", 'KIT_DDI')
◆ RELATION_KIT = / "BETON_DOUBLE_DP",
/ "BETON_GRANGER",
/ "BETON_GRANGER_V",
/ "BETON_UMLV",
/ "ENDO_ISOT_BETON",
/ "GLRC_DM",
/ "MAZARS",
/ "ROUSS_PR",
/ "VMIS_CINE_LINE",
/ "VMIS_ISOT_LINE",
/ "VMIS_ISOT_PUIS",
/ "VMIS_ISOT_TRAC",
# Si: equal_to("RELATION", 'KIT_CG')
◆ RELATION_KIT = / "CABLE_GAINE_FROT",
/ "ELAS",
/ "SANS",
/ "VMIS_CINE_LINE",
/ "VMIS_ISOT_LINE",
/ "VMIS_ISOT_TRAC",
# Si: is_in("RELATION", ['KIT_HHM','KIT_HH', 'KIT_H','KIT_HM','KIT_THHM', 'KIT_THH','KIT_THM','KIT_THV','KIT_THH2M','KIT_HH2M','KIT_HH2','KIT_THH2'])
◆ RELATION_KIT = / "Barcelone",
/ "CAM_CLAY",
/ "CJS",
/ "CSSM",
/ "CZM_EXP_REG",
/ "CZM_LIN_REG",
/ "DRUCK_PRAGER",
/ "DRUCK_PRAG_N_A",
/ "ELAS",
/ "ENDO_ISOT_BETON",
/ "GAZ",
/ "GonfElas",
/ "HOEK_BROWN_EFF",
/ "HOEK_BROWN_TOT",
/ "HUJEUX",
/ "HYDR_ENDO",
/ "HYDR_TABBAL",
/ "HYDR_UTIL",
/ "HYDR_VGC",
/ "HYDR_VGM",
/ "Iwan",
/ "JOINT_BANDIS",
/ "LAIGLE",
/ "LETK",
/ "LIQU_AD_GAZ",
/ "LIQU_AD_GAZ_VAPE",
/ "LIQU_GAZ",
/ "LIQU_GAZ_ATM",
/ "LIQU_SATU",
/ "LIQU_VAPE",
/ "LIQU_VAPE_GAZ",
/ "LKR",
/ "MAZARS",
/ "MCC",
/ "MFRONT",
/ "MOHR_COULOMB",
/ "NLH_CSRM",
/ "RANKINE",
/ "VISC_DRUC_PRAG",
/ "VISC_MAXWELL",
/ "VISC_MAXWELL_MT",
# Si: 'MFRONT' in value("RELATION_KIT")
◆ COMPOR_MFRONT = compor_mgis,
◇ VERI_BORNE = / "ARRET" (par défaut),
/ "MESSAGE",
/ "SANS",
◇ ALGO_CPLAN = / "ANALYTIQUE",
/ "DEBORST" (par défaut),
◇ RESI_INTE = float (défaut: 1e-08),
◇ ITER_INTE_MAXI = int,
◇ SYME_MATR_TANG = / "NON",
/ "OUI" (par défaut),
# Si: 'MFRONT' not in value("RELATION_KIT")
# Si: value("RELATION").startswith('META_') and not value("RELATION").startswith('META_LEMA_ANI')
◆ RELATION_KIT = / "ACIER",
/ "ZIRC",
◇ DEFORMATION = / "GDEF_LOG",
/ "GREEN_LAGRANGE",
/ "GROT_GDEP",
/ "PETIT" (par défaut),
/ "PETIT_REAC",
/ "SIMO_MIEHE",
◇ RESI_CPLAN_MAXI = float,
# Si: not exists("RESI_CPLAN_MAXI")
◇ RESI_CPLAN_RELA = float (défaut: 1e-06),
◇ ITER_CPLAN_MAXI = int (défaut: 10),
# Si: (equal_to("RELATION", 'MFRONT'))
◇ RESI_INTE = float,
◇ ITER_INTE_MAXI = int,
# Si: is_in("RELATION", ('RGI_BETON','FLUA_PORO_BETON','FLUA_ENDO_PORO', 'RGI_BETON_BA',))
◇ RESI_INTE = float (défaut: 1e-06),
◇ ITER_INTE_MAXI = int (défaut: -1),
# Si: not is_in("RELATION", ('MFRONT','RGI_BETON','FLUA_PORO_BETON','FLUA_ENDO_PORO','RGI_BETON_BA',))
◇ RESI_INTE = float (défaut: 1e-06),
◇ ITER_INTE_MAXI = int (défaut: 20),
# Si: is_in("DEFORMATION", ('PETIT','PETIT_REAC','GROT_GDEP'))
◇ ITER_INTE_PAS = int (défaut: 0),
◇ ALGO_INTE = / "ANALYTIQUE",
/ "BASCULE_EXPLICITE",
/ "BRENT",
/ "DEKKER",
/ "NEWTON",
/ "NEWTON_1D",
/ "NEWTON_PERT",
/ "NEWTON_RELI",
/ "RUNGE_KUTTA",
/ "SANS_OBJET",
/ "SECANTE",
/ "SEMI_EXPLICITE",
/ "SPECIFIQUE",
# Si: not is_in("RELATION", ('RGI_BETON','FLUA_PORO_BETON','FLUA_ENDO_PORO', 'RGI_BETON_BA', 'ENDO_PORO_BETON'))
◇ TYPE_MATR_TANG = / "PERTURBATION",
/ "VERIFICATION",
# Si: (exists("TYPE_MATR_TANG"))
◇ VALE_PERT_RELA = float (défaut: 1e-05),
# Si: equal_to("TYPE_MATR_TANG", 'TANGENTE_SECANTE')
◇ SEUIL = float (défaut: 3.0),
◇ AMPLITUDE = float (défaut: 1.5),
◇ TAUX_RETOUR = float (défaut: 0.05),
◇ PARM_THETA = float (défaut: 1.0),
# Si: not exists("TYPE_MATR_TANG")
◇ RESI_RADI_RELA = float,
# Si: is_in("RELATION", ('RGI_BETON','FLUA_PORO_BETON','FLUA_ENDO_PORO', 'RGI_BETON_BA', 'ENDO_PORO_BETON'))
◇ TYPE_MATR_TANG = / "MATR_ELAS" (par défaut),
/ "MATR_ENDO",
◇ PARM_THETA = float (défaut: 1.0),
# Si: not exists("TYPE_MATR_TANG")
◇ RESI_RADI_RELA = float,
),
⟐ NEWTON = _F(
◇ REAC_INCR = int (défaut: 1),
◇ PREDICTION = / "DEPL_CALCULE",
/ "ELASTIQUE",
/ "EXTRAPOL",
/ "TANGENTE",
◇ MATRICE = / "ELASTIQUE",
/ "TANGENTE" (par défaut),
◇ PAS_MINI_ELAS = float (défaut: 0.0),
◇ REAC_ITER = int (défaut: 0),
◇ REAC_ITER_ELAS = int (défaut: 0),
◇ EVOL_NOLI = evol_noli,
),
⟐ CONVERGENCE = _F(
# Si: exists("RESI_REFE_RELA")
◆ | SIGM_REFE = float,
| EPSI_REFE = float,
| FLUX_THER_REFE = float,
| FLUX_HYD1_REFE = float,
| FLUX_HYD2_REFE = float,
| VARI_REFE = float,
◇ / RESI_REFE_RELA = float,
/ | RESI_GLOB_MAXI = float,
| RESI_GLOB_RELA = float,
◇ ITER_GLOB_MAXI = int (défaut: 10),
◇ ITER_GLOB_ELAS = int (défaut: 25),
◇ ARRET = / "NON",
/ "OUI" (par défaut),
),
◇ ENERGIE = _F(
◇ CALCUL = "OUI",
),
◆ EXCIT_THER = _F(
◆ CHARGE = char_cine_ther / char_ther,
◇ FONC_MULT = fonction / formule / nappe,
),
)
◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles
Opérandes#
Opérandes TABLE / TEMPER / DEBIT#
◊ TEMPER = CO(“temp”)
◊ TABLE = CO(“table”)
◊ DEBIT = CO(“debit”)
Ces mots-clés permettent de donner le nom des concepts sortant de la macro-commande, à savoir
TABLE : une table de 7 colonnes donnant pour chaque pas de temps : la cote \(z\) , le flux, la pression totale (\(\mathrm{Pa}\) ), la température du fluide (\(°C\) ), le coefficient d’échange par convection (\({\mathit{W.m}}^{-2.}°{C}^{-1}\) ), le nom de la fissure.
TEMPER : la structure de données résultat du calcul thermique
DEBIT : la valeur du débit (\({\mathit{kg.s}}^{-1}\) ) à chaque pas de temps et pour chaque fissure
Remarques:
|
Mot cléETAT_INIT#
◊ ETAT_INIT = _F ( ♦ EVOL_NOLI = evol_noli ♦ EVOL_THER = evol_ther ♦ NUME_ORDRE = nume_ordre
Mot-clé permettant de définir l’état initial dans le cadre d’une poursuite de la macro-commande.
Opérande EVOL_NOLI#
Nom du concept de type evol_noli d’où sera extrait l’état mécanique.
Opérande EVOL_THER#
Nom du concept de type evol_ther d’où sera extrait l’état thermique. Ce concept sera enrichi au cours du calcul.
Opérande NUME_ORDRE#
L’état initial pour poursuivre le calcul sera défini à partir du numéro d’archivage NUME_ORDRE pour evol_noli et pour evol_ther.
Opérandes MODELE_MECA / MODELE_THER#
♦ MODELE_MECA = m_meca
♦ MODELE_THER = m_ther
Ces mots-clés permettent de renseigner le nom du modèle (m_meca) dont les éléments font l’objet du calcul mécanique et le nom du modèle (m_ther) dont les éléments font l’objet du calcul thermique.
Opérande LIST_INST#
♦LIST_INST = list_inst
Liste des pas de temps qui correspondent aux instants de calcul imposés. En cas de redécoupage, les instants sont insérés dans cette liste.
Opérande CHAM_MATER#
♦ CHAM_MATER = chmat
Nom du champ de matériau affecté au maillage. Attention, ce champ doit comprendre les données associées au comportement mécanique et au comportement thermique (mot-clé THER de DEFI_MATERIAU).
Opérande TEMP_INIT#
♦ TEMP_INIT = temp_init
Valeur de la température initiale sur tout le domaine.
Opérande CARA_ELEM#
◊ CARA_ELEM = carac
Ce mot-clé permet de renseigner, les caractéristiques des éléments de poutre, barre, coque, tuyau, élément discret, lorsqu’ils sont présents dans le modèle. Attention, ce mot-clé n’est utilisé que pour le calcul mécanique mais n’est pas transmis pour l’opérateur de thermique (pour éviter tout problème avec les barres).
Mot clé EXCIT_MECA#
♦ EXCIT_MECA = _F ( ♦ CHARGE = chi ◊ FONC_MULT = fi ◊ TYPE_CHARGE = / « FIXE_CSTE » / « SUIV » / « DIDI » )
Il s’agit ici de renseigner les conditions aux limites du problème mécanique, ainsi que les éventuels chargements mécaniques qui s’appliquent à la structure. La syntaxe est identique à celle du mot clé EXCIT de STAT_NON_LINE mais il n’est pas possible de piloter un chargement. Voir le document [U4.51.03].
Mot clé CONTACT#
♦ CONTACT = char_contact
On renseigne ici les conditions de contact entre les lèvres de la fissure, qui peuvent se refermer pendant le calcul. Le but de ce mot-clé est celui d’éviter l’interpénétration des lèvres qui peut se produire sous l’effet de l’écoulement du fluide.
Mot clé EXCIT_THER#
♦ EXCIT_THER = _F ( ♦ CHARGE = chi ◊ FONC_MULT = fi )
Il s’agit ici de renseigner les conditions aux limites du problème thermique ainsi que les éventuels chargements thermiques. La syntaxe est identique à celle du mot clé EXCIT de THER_LINEAIRE.
Voir le document [U4.54.01].
Mot clé COMPORTEMENT#
Il s’agit ici de renseigner la loi de comportement qui sera utilisé pour résoudre le problème mécanique. Voir le document [U4.51.03].
Mot clé NEWTON#
Il s’agit ici de renseigner les paramètres de l’algorithme de Newton pour résoudre le problème mécanique. Ce mot clé est identique à celui de STAT_NON_LINE, il convient donc de se reporter au document [U4.51.03].
Mot clé CONVERGENCE#
Il s’agit ici de renseigner les critères de convergence pour le calcul mécanique. La syntaxe correspond à celle de STAT_NON_LINE (voir le document [U4.51.03]).
Remarque:
Il peut être nécessaire d’utiliser RESI_GLOB_MAXI plutôt que RESI_GLOB_RELA , lorsqu’il n’y a pas de chargement mécanique (sinon, effort nul, donc RESI_GLOB_RELA indéfini). |
Mot clé ENERGIE#
Ce mot-clé permet d’activer le calcul du bilan d’énergie et son affichage au cours du calcul mécanique (voir le document [R4.09.01]). Ce bilan est stocké dans la table de nom PARA_CALC d’où il peut être extrait à l’aide de la commande RECU_TABLE [U4.71.02].
Mot cléFISSURE#
♦ FISSURE
Ce mot clé permet de définir tous les paramètres ayant trait à une fissure. Plusieurs occurrences de ce mot-clé sont possibles s’il y a plusieurs fissures.
Opérande GROUP_MA#
♦ GROUP_MA = gma
Groupes des mailles définissant la ou les lèvres de la fissure. On donne le groupe de la lèvre inférieure et supérieure.
Opérandes GROUP_NO_ORIG et GROUP_NO_EXTR#
♦ GROUP_NO_ORIG = ogno
♦ GROUP_NO_EXTR = egno
Permet de définir les deux groupes de noeuds qui définissent les deux extrémités de la fissures afin de l’orienter.
Opérande SECTION#
♦ SECTION = /"ELLIPSE"
/ »RECTANGLE »
Il s’agit ici de définir la section de passage de l’écoulement (plan perpendiculaire au plan de modélisation). Celle-ci peut être elliptique ou rectangulaire. L’une des dimensions de cette section correspond à la taille de la fissure dans le plan de la modélisation. Celle-ci est déterminée à chaque instant par la macro-commande. L’autre dimension correspond à la taille de la fissure hors du plan de modélisation. Celle-ci est fixe pendant tout le calcul et doit être renseigné grâce aux opérandes LISTE_COTES_BL et LIST_VAL_BL.
Opérande RUGOSITE#
♦ RUGOSITE = eps
Rugosité absolue de la paroi (en mètres). Cette valeur n’influe le résultat que si le frottement est calculé par Ecrevisse (opérande FROTTEMENT >0) et si l’écoulement n’est plus laminaire.
Pour les options de frottement \(21,22,23,24\) il s’agit du diamètre du plus gros granulat (écoulement dans une fissure dans du béton).
Opérandes LISTE_COTES_BL / LISTE_VAL_BL#
◊ LISTE_COTES_BL = lcbl
♦ LISTE_VAL_BL = lvbl
Figure 3.2.5-1: Représentation de la fissure
Pour caractériser la taille de la fissure dans le plan non modélisé, il faut fournir la valeur du«petit» axe de la section elliptique ou la largeur de la section rectangulaire au minimum en 2 points.
LISTE_COTES_BL permet de renseigner les abscisses curvilignes où l’on donne la dimension de la fissure et LIST_VAL_BL permet de donner la dimension. Si LISTE_COTES_BL n’est pas renseigné, LISTE_VAL_BL doit contenir la largeur de la fissure à l’entrée et à la sortie de la fissure.
Opérande ZETA#
♦ ZETA = zeta
Coefficient de perte de charge singulière à l’entrée (sans dimension).
Opérande OUVERT_REMANENTE#
♦ OUVERT_REMANENTE = ouv_rem
Valeur de l’ouverture rémanente (en mètres) pour le calcul Ecrevisse correspondant à l’ouverture hydraulique réelle de la fissure lorsque les deux lèvres sont au contact. Ceci traduit le fait que, même lorsque la fissure est refermée, une petite quantité de fluide transitera, en particulier à cause de la rugosité. Ce paramètre dépend du matériau : par exemple, c’est de l’ordre de \(10\mu m\) pour le béton.
Ce mot-clé est donc pris en compte donc lors du calcul hydraulique (Ecrevisse), mais non par Code_Aster : il peut donc bien arriver que le calcul hydraulique soit effectué avec l’ouverture rémanente, est que le contact soit activé dans le calcul mécanique.
Opérande TORTUOSITE#
◊ TORTUOSITE = tort
Si la fissure est tortueuse, on peut admettre que, vue du fluide, la longueur de la fissure est supérieure à l’épaisseur de la paroi traversée. Ce coefficient, qui par défaut vaut 1, permet de prendre en compte ce phénomène. La longueur de la fissure sera donc égale à \(\frac{L}{\mathit{tort}}\) avec \(\mathit{tort}\le 1\)
Opérande PREFIXE_FICHIER#
◊ PREFIXE_FICHIER = prefix
Permet de définir le nom des fichiers de sortie Ecrevisse dans le dossier REPE_OUT (courbes et fichiers textes), et ainsi de distinguer les fichiers associés à chaque fissure.
Mot cléECOULEMENT#
♦ ECOULEMENT
Il s’agit dans ce mot clé de définir la nature du fluide et de ses caractéristiques.
Opérandes PRES_ENTREE / PRES_ENTREE_FO#
♦ / PRES_ENTREE = pe
/ PRES_ENTREE_FO = fpe
Cet opérande permet de donner la pression totale (\(\mathit{Pa}\) ) ou son évolution au cours du temps, en amont de la fissure. Elle doit être comprise entre \(215.{10}^{5}\mathit{Pa}\) et \({10}^{5}\mathit{Pa}\) .
Opérandes PRES_SORTIE / PRES_SORTIE_FO#
♦ / PRES_SORTIE = pe
/ PRES_SORTIE_FO = fpe
Cet opérande permet de donner la pression totale (\(\mathit{Pa}\) ) ou son évolution au cours du temps, en sortie de la fissure. Elle doit être comprise entre \(215.{10}^{5}\mathit{Pa}\) et \(800\mathit{Pa}\) et être inférieure à la pression d’entrée (sinon, il n’y a pas d’écoulement et le calcul s’arrête).
Opérande FLUIDE_ENTREE#
♦ FLUIDE_ENTREE = 1,2,3,4,5 ou 6
Cet opérande permet de définir l’état du fluide à l’entrée de la fissure :
1 : Eau sous-refroidie ou saturée;
2 : Fluide diphasique;
3 : Vapeur saturée ou surchauffée;
4 : Air + vapeur surchauffée;
5 : Air + vapeur saturée;
6 : Air seul.
En fonction de la situation du fluide, seules certaines des caractéristiques suivantes sont à renseigner.
Opérandes TEMP_ENTREE / TEMP_ENTREE_FO#
♦ / TEMP_ENTREE = te
/ TEMP_ENTREE_FO = fte
Cet opérande permet de renseigner la température (\(°C\) ) ou son évolution au cours du temps, en amont de la fissure. Il doit être renseigné si le fluide à l’entrée est «Eau sous-refroidie ou saturée», «Vapeur saturée ou surchauffée», «Air + vapeur surchauffée» ou «Air seul».
(opérande FLUIDE_ENTREE = 1,3,4,ou 6).
Opérandes TITR_MASS / TITR_MASS_FO#
♦ / TITR_MASS = xe
/ TITR_MASS_FO = fxe
Cet opérande permet de renseigner le titre massique de vapeur, ou son évolution au cours du temps, en amont de la fissure. Il correspond au rapport entre la masse de vapeur et la masse d’eau liquide et vapeur. Il doit être renseigné si le fluide à l’entrée est «Fluide diphasique» ou «Air + vapeur saturée» (opérande FLUIDE_ENTREE = 2 ou 5).
Opérandes PRES_PART / PRES_PART_FO#
♦ / PRES_PART = pae
/ PRES_PART_FO = fpae
Cet opérande permet de renseigner la pression partielle d’air (\(\mathit{Pa}\) ) ou son évolution au cours du temps, en amont de la fissure. Il doit être renseigné si le fluide à l’entrée est «Air + vapeur surchauffée» ou «Air + vapeur saturée» (opérande FLUIDE_ENTREE=4 ou 5).
Mot cléMODELE_ECRE#
♦ MODELE_ECRE
Ce mot clé permet de renseigner les modèles qui seront utilisés par Ecrevisse pour calculer l’écoulement du fluide.
Opérande IVENAC#
◊ IVENAC = 0 ou 1,
Lorsque cet opérande est activé (IVENAC=1), le calcul ECREVISSE s’effectue avec prise en compte de la vena contracta (perte de charge en entrée). Sinon (cas par défaut), le calcul s’effectuera sans cette option de modélisation.
Opérande ECOULEMENT#
♦ ECOULEMENT = / 'SATURATION',
/ “GELE”
En l’absence d’air et en présence d’eau (opérande FLUIDE_ENTREE=1 ou 2), cet opérande permet de choisir entre le modèle homogène à l’équilibre et le modèle d’écoulement «GELE» caractérisé par une fraction non nulle de liquide métastable.
Opérande PRESS_EBULLITION#
Pression d’ébullition (\(\mathit{Pa}\) ) à fournir uniquement pour le modèle d’écoulement gelé.
Opérande FROTTEMENT#
Il s’agit ici de définir comment est calculé le frottement. Les valeurs valides sont: \(-14,-12,-11,-4,-3,-2,-1,0,1,2,3,4,11,12,14,21,22,23,24\) .
La valeur \(0\) correspond à un calcul sans frottement.
Pour les valeurs négatives, l’utilisateur fixe la valeur du coefficient de frottement pour les écoulements turbulents.
Pour les valeurs positives, le frottement en écoulement turbulent est calculé à partir du coefficient de rugosité. En particulier, les options \(21,22,23,24\) correspondent à la relation de frottement de Greiner et Ramm [Bib.6] pour l’écoulement dans une fissure dans le béton; la rugosité est dans ce cas le diamètre du plus gros granulat. De plus, pour les valeurs supérieures à \(10\) , la loi utilisée pour l’écoulement laminaire (frottement proportionnel à l’inverse du nombre de Reynolds) est raccordée avec celle en écoulement turbulent (valeur constante du coefficient de frottement). Autrement, une discontinuité est présente.
Le deuxième chiffre (1,2,3,4) de chaque option (unités) détermine l’option de calcul de la viscosité dynamique pour les écoulements diphasiques, cf. documentation Ecrevisse [Bib.1] [Bib.2] [Bib.3] [Bib.4] [Bib.5] .
Opérandes REYNOLDS_LIM / FROTTEMENT_LIM#
Coefficient de Reynolds limite et coefficient de frottement imposé pour un Reynolds supérieur au Reynolds limite.
A fournir seulement si FROTTEMENT < 0.
Opérande TRANSFERT_CHAL#
Cet opérande permet de définir si on veut ou non prendre en compte le transfert de chaleur par convection entre le fluide et la paroi. Les valeurs valides sont: \(-12,-11,-2,-1,0,1,2,11,12\) .
La valeur 0 correspond à un calcul sans transfert de chaleur (adiabatique).
Les autres valeurs correspondent à différentes options pour le calcul du coefficient de convection (différences uniquement pour le régime laminaire). Se reporter à la documentation Ecrevisse.
Opérandes XMINCH / XMAXCH#
Cet opérateur n’a de sens que dans les cas diphasiques, si TRANSFERT_CHAL < 0. Il permet de renseigner la valeur des titres massiques en gaz délimitant la zone de transition entre les modèles de calcul du coefficient de convection, cf. documentation Ecrevisse.
Mot clé CONV_CRITERE#
La macro-commande calcule à chaque instant trois critères permettant d’estimer l’importance des changements intervenus entre deux pas de temps et ainsi la validité du chaînage effectué. Pour calculer ces critères, l’utilisateur doit définir des valeurs de référence pour la pression et la température, permettant de quantifier l’écart acceptable.
\(\begin{array}{c}{e}_{T}=\frac{\mathit{Max}({T}_{t}-{T}_{\text{t-1}})}{{T}_{\text{REF}}}\\ {e}_{P}=\frac{\mathit{Max}({P}_{t}-{P}_{\text{t-1}})}{{P}_{\text{ref}}}\\ {e}_{G}=\sqrt{{e}_{T}^{2}+{e}_{P}^{2}}\end{array}\)
Ensuite, il peut activer s’il le souhaite le redécoupage du pas de temps si la valeur du critère est inférieure à une valeur donnée.
L’expérience acquise montre toutefois qu’il vaut mieux essayer d’optimiser la liste d’instant que de compter sur le redécoupage du pas de temps, car les calculs Ecrevisse sont relativement longs.
.
Opérande CRITERE#
Permet de définir la nature du critère utilisé pour gérer les pas de temps dans la macro-commande. Les valeurs possibles sont :
EXPLICITE: pas de redécoupage quelque soit la valeur des indicateurs
TEMP: il y aura redécoupage en fonction de la valeur de \({e}_{T}\)
PRES: il y aura redécoupage en fonction de la valeur de \({e}_{P}\)
TEMP_PRES: il y aura redécoupage en fonction de la valeur de \({e}_{G}\)
OpérandeTEMP_REF#
Température de référence pour le calcul du critère en température de la macro-commande.
OpérandePRES_REF#
Pression de référence pour le calcul du critère en pression de la macro-commande.
Opérande PREC_CRIT#
Valeur à laquelle on compare l’erreur obtenue pour activer ou non le redécoupage du pas de temps (1 par défaut). Par exemple pour la température: \({e}_{T}<\text{PREC\_CRIT}\)
Opérande NUME_ORDRE_MIN#
Numéro d’ordre à partir duquel le critère d’erreur est pris en compte. La convergence est forcée sur les numéros d’ordre inférieurs. On utilise souvent cette opérande pour démarrer l’activation du critère d’erreur à partir du deuxième pas de temps et forcer la convergence sur le premier pas qui est souvent délicat à traiter, puisqu’on impose des chargements brutaux à l’instant initial.
Opérande SUBD_NIVEAU#
Nombre de redécoupage du pas de temps autorisé. Au delà de cette valeur, on ne redécoupe plus et on sort de la macro-commande.
OpérandeSUBD_PAS_MINI#
Pas de temps en dessous duquel on ne découpe plus.
Mot cléCONVERGENCE_ECREVISSE#
Opérande KGTEST#
Paramètre de l’algorithme itératif d’Ecrevisse.
Doit être compris entre \(0\) et \(1.\) . Fixé par défaut à \(0,5\) , on peut être amené à mettre une valeur légèrement inférieure si des difficultés de convergence sur le débit sont observées.
Opérande ITER_GLOB_MAXI#
Nombre maximum d’itérations autorisé pour le calcul du débit (400 par défaut).
Opérande CRIT_CONV_DEBI#
Précision utilisée pour la convergence du calcul du débit. C’est la valeur par rapport à laquelle Ecrevisse teste le débit max et min dans son procédé de calcul.
\(\frac{{G}_{max}-{G}_{min}}{G}⩽\text{CRIT\_CONV\_DEBI}\)
Par défaut à \({10}^{\text{-5}}\) , des valeurs jusqu’à \({10}^{\text{-2}}\) améliore le temps de calcul sans perdre beaucoup de précision.
Opérande COURBES#
◊ COURBES = / « AUCUNE » [DEFAUT] / « POSTSCRIPT » / « INTERACTIF »
Permet de générer ou non les courbes de sortie d” Ecrevisse. Si les courbes sont générées, elles sont stockées dans le dossier REPE_OUT au format postscript (tous les instants sont dans le même fichier postscript). En mode interactif, les courbes sont affichées pour chaque pas de temps.
Opérande LOGICIEL#
Permet de spécifier «en dur» le chemin pour trouver l’exécutable d’Ecrevisse. A présent, seules les versions 3.2, 3.2.1 et 3.2.2 sont compatibles avec la MACR_ECREVISSE. Il faut renseigner obligatoirement cet opérande ou alternativement l’opérande VERSION.
Opérande VERSION#
Indique la version d’Ecrevisse utilisée dans le couplage. Il faut renseigner obligatoirement cet opérande ou alternativement l’opérande LOGICIEL. A présent, la version 3.2.2 est utilisée. Elle est la seule qu’on peut renseigner sous l’opérande VERSION. Si on souhaite utiliser une autre version (3.2 ou 3.2.1), celle-ci doit être renseignée sous l’opérande LOGICIEL.
Opérande ENTETE#
Sert à donner un titre au calcul Ecrevisse.
Opérande IMPRESSION#
Lorsque l’opérande vaut “OUI”, la macro imprime une table récapitulative des résultats d’Ecrevisse le long de l’abscisse curviligne (température, pression, flux, vitesse, coefficient de convection, …).
Opérande INFO#
INFO=1: presque aucun affichage.
INFO=2: affiche les informations concernant le critère de convergence.
Fonctionnement du couplage#
Principe#
Cette macro-commande a été développée dans un cadre précis : estimer les débits d’une fissure dans un massif en béton modélisé en 2D, contenant éventuellement des aciers modélisés sous forme de barre (non pris en compte dans le calcul thermique). L’utilisation dans tout autre cadre doit donc se faire avec la plus grande précaution. Toutefois, il y a en fait 3 macro-commande imbriquées, ce qui peut permettre de n’utiliser qu’une partie des fonctionnalités.
L’algorithme est le suivant :
Début de la macro Boucle sur la liste d’instants MACR_ECREVISSE.py Récupération chargement thermique et mécanique Exécution de THER_LINEAIRE Projection du chargement thermique et définition des chargements venant d’Ecrevisse Exécutionde STAT_NON_LINE Entrée dans CALC_ECREVISSE.py Récupération de l’ouverture de la fissure (POST_RELEVE_T) Récupération de la température en bord de fissure Entrée dans MACR_ECRE_CALC.py Exécution d’ECREVISSE Récupération du débit calculé et des tables de résultats Extraction et copie des données sortantes (CL,thermique, résultat Ecrevisse…) Évaluation du critère de convergence (erreur en température et/ou pression) Redécoupage éventuel du pas de temps voire arrêt suivant les conditions Récupération des infos et concepts (chargement,…) |
Le détail des macro-commandes intermédiaires, CALC_ECREVISSE et MACR_ECRE_CALC, est donné en annexe.
Unités#
On peut habituellement se passer des unités dans Code_Aster en restant cohérent. Cependant, pour les besoins d’Ecrevisse, il est important ici de respecter le système d’unités international (SI), en particulier pour les longueurs (\(m\) ), les pressions (\(\mathit{Pa}\) ), et les températures (\(°C\) ).
Quelques conseils d’utilisation#
Le maillage doit être en 2 dimensions, et la fissure maillée explicitement. Pour que le calcul soit de bonne qualité, il est nécessaire d’optimiser le maillage pour qu’il soit compatible avec les phénomènes que l’on veut représenter : les chargements mécaniques, les chargements thermiques mais aussi l’écoulement du fluide. Il est ainsi important de raffiner suffisamment le maillage aux bords de la fissure si on veut voir s’échauffer le matériau. En entrée de fissures, les phénomènes thermiques et hydrauliques sont rapides, il est donc nécessaire là aussi de raffiner suffisamment le maillage.
La discrétisation en temps est également très importante et doit être suffisamment fine, sous peine d’avoir des écarts importants dans les réponses. Mieux vaut essayer d’optimiser cette liste d’instants plutôt que de se reposer uniquement sur les critères pour redécouper car cette stratégie est beaucoup plus coûteuse en temps.
Pour des calculs avec condensation de la vapeur sur le chemin de fissuration, le coefficient de convection peut augmenter jusqu’à deux ordres de grandeurs au passage à l’état liquide. Comme la condensation a lieu sur une zone très petite en espace, la variation du coefficient de convection et donc du flux de chaleur transféré à la paroi peut être très rapide (voir la ).
Cela peut engendrer des oscillations spatiales dans la température du massif.
Figure 4.3-1: Coefficients de convection à plusieurs instants d’un calcul avec condensation de la vapeur.
La macro-commande peut être réentrante . Les résultats thermiques et mécaniques seront bien enrichis, par contre de nouvelles tables de sortie seront créés pour la nouvelle liste de pas de temps.
Remarque :
Les chargements sont détruits à chaque pas de temps, ainsi que le champ de matériaux, les post-traitements possibles sont de fait limités (déplacement, contrainte et variables internes aux points de Gauss). |
Exemple d’utilisation#
Voir tests zzzz218a, b et c.
Voir également les tests zzzz354(a-h) et zzzz355(a-h).
Référence Ecrevisse#
Hervouet (2004). «Simulation d’écoulement air/eau/vapeur à travers une fissure : logiciel ECREVISSE version 3.0. Note de principe». Note interne EDF no. H-I86-2004-00741-FR.
Hervouet (2008). «Modifications apportées au logiciel ECREVISSE version 3.1 pour obtenir ECREVISSE version 3.1.1.». Note interne EDF no. H-I81-2008-03647-FR.
Hervouet (2009). «Simulation d’écoulement air/eau/vapeur à travers une fissure : logiciel ECREVISSE version 3.1.1 Manuel utilisateur». Note interne EDF no. H-I81-2009-00243-FR.
Hervouet (2016). «Modifications apportées au logiciel ECREVISSE version 3.2.1 pour obtenir ECREVISSE version 3.2.2.». Note interne EDF no. H-I8B-2016-00390-FR.
Hervouet (2016). «Simulation en 1D d’écoulement à travers une fissure : logiciel ECREVISSE version 3.2.2. Manuel utilisateur». Note interne EDF no. H-I8B-2016-00399-FR.
Greiner, W. Ramm (1995). «Air leakage characteristics in cracked concrete», Nuclear Engineering and Design, Volume 156, Issue 1, Pages 167-172.
Annexes#
Syntaxe de la procédure CALC_ECREVISSE#
CALC_ECREVISSE =
(
# CONCEPT SORTANT
♦ CHARGE_MECA = CO(“char_meca”) [CO] ♦ CHARGE_THER1 = CO(“char_ther”) [CO] ♦ CHARGE_THER2 = CO(“char_ther”) [CO] ♦ TABLE = CO(“table”) [CO] ♦ DEBIT = CO(“debit”) [CO]
# MODELES MECANIQUES
♦ MODELE_MECA = m_meca [modele] ♦ MODELE_THER = m_ther [modele]
# DONNEES GEOMETRIQUES RELATIVES AUX RESULTATS
♦ RESULTAT = _F ( ♦ MECANIQUE = rmeca [resultat] ♦ THERMIQUE = rther [resultat] / ◊ NUME_ORDRE = nume_ordre [I] ◊ INST = instant [R] ),
# DONNEES GEOMETRIQUES RELATIVES A LA FISSURE
♦ FISSURE = _F ( ♦ GROUP_MA = gma [grma] ♦ GROUP_NO_ORIG = ogno [grno] ♦ GROUP_NO_EXTR = egno [grno] ♦ ZETA = zeta [R] ♦ RUGOSITE = rug [R] ♦ OUVERT_REMANENTE = ouv_rem [R] ◊ PREFIXE_FICHIER = / “FISSURE1”, [DEFAUT] /prefix ◊ TORTUOSITE = tort [R] ♦ SECTION = / « ELLIPSE » / « RECTANGLE » ◊ LISTE_COTES_BL = / (0, max(abs_curv)) [DEFAUT] / lcbl listr8 ♦ LISTE_VAL_BL = lvbl listr8 ),
# DONNEES RELATIVES A L’ECOULEMENT
♦ ECOULEMENT = … idem MACR_ECREVISSE/ECOULEMENT
♦ MODELE_ECRE = … idem MACR_ECREVISSE/MODELE_ECRE
# DONNEES RELATIVES A LA CONVERGENCE NUMERIQUE
◊ CONVERGENCE_ECREVISSE = … idem MACR_ECREVISSE/CONVERGENCE_ECREVISSE
# GENERAL
COURBES, LOGICIEL, VERSION… = … idem MACR_ECREVISSE/GENERAL
◊ INFO = / 1 [DEFAUT] / 2 )
Opérandes spécifiques à CALC_ECREVISSE#
Opérande CHARGE_MECA#
Concept sortant contenant le chargement mécanique fourni par Ecrevisse.
Opérandes CHARG_THER1 et CHARG_THER2#
Concepts sortants contenant les chargements thermiques fournis par Ecrevisse.
Mot cléRESULTAT#
Opérande MECANIQUE#
Contient le résultat de STAT_NON_LINE
Opérande THERMIQUE#
Contient le résultat de THER_LINEAIRE
Opérande NUME_ORDRE#
Numéro d’ordre du pas de temps de calcul en cours.
Opérande INST#
Valeur de l’instant (en secondes) du pas de temps de calcul en cours.
Syntaxe de la procédure MACR_ECRE_CALC#
MACR_ECRE_CALC =
(
# CONCEPT SORTANT
♦ TABLE = CO(“table”) [CO] ♦ DEBIT = CO(“debit”) [CO]
# DONNEES GEOMETRIQUES RELATIVES A LA FISSURE
♦ FISSURE = _F ( ♦ LONGUEUR = long [R] ♦ RUGOSITE = rug [R] ♦ ANGLE = alpha [R] ♦ ZETA = zeta [R] ♦ SECTION = / « ELLIPSE » / « RECTANGLE » ♦ LISTE_COTES_BL = lcbl listr8 ♦ LISTE_VAL_BL = lvbl listr8 ♦ LISTE_COTES_AH = lcah listr8 ♦ LISTE_VAL_AH = lvah listr8 ),
# DONNEES RELATIVES A L’ECOULEMENT
♦ ECOULEMENT = idem MACR_ECREVISSE/ECOULEMENT
# DONNEES RELATIVES AU PROFIL DE TEMPERATURE A TRAVERS LA PAROI
◊ TEMPERATURE = _F ( / ♦ GRADIENT = « FOURNI » ♦ LISTE_COTES_TEMP = lct [R] ♦ LISTE_VAL_TEMP = lvt [R] / ♦ GRADIENT = « IMPOSE » ♦ TEMP1 = tm1 [R] ♦ TEMP2 = tm2 [R] / ♦ GRADIENT = « CALCULE » ♦ EPAISSEUR_PAROI = epp [R] ♦ CONVECTION_AMONT = alphe [R] ♦ CONVECTION_AVAL = alphs [R] ♦ LAMBDA = lambd [R] ♦ TEMP_FLUIDE_AVAL = ts [R] ),
# DONNEES RELATIVES A L’ECOULEMENT
♦ ECOULEMENT = … idem MACR_ECREVISSE/ECOULEMENT
♦ MODELE_ECRE = … idem MACR_ECREVISSE/MODELE_ECRE
# DONNEES RELATIVES A LA CONVERGENCE NUMERIQUE
◊ CONVERGENCE_ECREVISSE = … idem MACR_ECREVISSE/CONVERGENCE_ECREVISSE
# GENERAL
COURBES, LOGICIEL, VERSION… = … idem MACR_ECREVISSE/GENERAL )
Opérandes spécifiques à MACR_ECRE_CALC#
Mot cléFISSURE#
Opérande ANGLE#
Angle de l’écoulement par rapport à l’ascendante verticale (degrés).
Opérande LONGUEUR#
Définit la longueur réelle de la fissure (tenant compte de sa tortuosité)
Opérande LISTE_COTES_AH#
LISTE_COTES_AH permet de renseigner les abscisses curvilignes où l’on donne la dimension du grand axe de la fissure.
Opérande LISTE_VAL_AH#
LISTE_VAL_AH permet de renseigner les valeurs du grand axe de section de la fissure / du conduit. Selon la valeur de l’opérande SECTION:
ELLIPSE : Liste des valeurs des points définissant le grand axe de la section;
RECTANGLE : Liste des valeurs des points définissant la longueur de la section.
Mot cléTEMPERATURE#
Opérande GRADIENT#
Modèle du gradient de température :
Fourni : Distribution de température fournie;
Imposé : Distribution imposée de température;
Calculé : Profil de température calculé;
Opérande LISTE_COTES_TEMP#
Liste des cotes pour les températures.
Opérande LISTE_VAL_TEMP#
Liste des valeurs de température.
Opérande TEMP1#
Gradient de température de la paroi le long de l’écoulement.
Opérande TEMP2#
Température de la paroi à l’entrée.
Opérande EPAISSEUR_PAROI#
Épaisseur de la paroi.
Opérande CONVECTION_AMONT#
Coefficient de convection à la surface de la paroi coté amont.
Opérande CONVECTION_AVAL#
Coefficient de convection à la surface de la paroi coté aval.
Opérande LAMBDA#
Conduction thermique de la paroi.
Opérande TEMP_FLUIDE_AVAL#
Température du fluide coté aval.