u4.84.21 Opérateur CALC_FORC_NONL#
Syntaxe#
Détail de la syntaxe
dyna_trans = CALC_FORC_NONL(
◆ RESULTAT = resultat,
◇ / TOUT_ORDRE = "OUI" (ou non renseigné),
/ NUME_ORDRE = int,
/ INST = float,
/ FREQ = float,
/ NUME_MODE = int,
/ NOEUD_CMP = text,
/ LIST_INST = listr8,
/ LIST_FREQ = listr8,
/ LIST_ORDRE = listis,
/ NOM_CAS = text,
◇ CRITERE = / "RELATIF" (par défaut),
/ "ABSOLU",
# Si: (equal_to("CRITERE", 'RELATIF'))
◇ PRECISION = float (défaut: 1e-06),
# Si: (equal_to("CRITERE", 'ABSOLU'))
◆ PRECISION = float,
◆ MODELE = modele,
◇ CHAM_MATER = cham_mater,
◇ CARA_ELEM = cara_elem,
⟐ 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 = / "ELAS" (par défaut),
/ "ELAS_VMIS_LINE",
/ "ELAS_VMIS_TRAC",
/ "ELAS_VMIS_PUIS",
/ "ELAS_HYPER",
/ "ELAS_POUTRE_GR",
/ "CABLE",
/ "ARME",
/ "ASSE_CORN",
/ "BARCELONE",
/ "BETON_DOUBLE_DP",
/ "BETON_RAG",
/ "BETON_REGLE_PR",
/ "BETON_UMLV",
/ "CABLE_GAINE_FROT",
/ "CAM_CLAY",
/ "CJS",
/ "CORR_ACIER",
/ "CZM_EXP_REG",
/ "CZM_EXP_MIX",
/ "CZM_FAT_MIX",
/ "CZM_LIN_REG",
/ "CZM_OUV_MIX",
/ "CZM_TAC_MIX",
/ "CZM_ELAS_MIX",
/ "CZM_LAB_MIX",
/ "CZM_TRA_MIX",
/ "CZM_TURON",
/ "DASHPOT",
/ "DIS_BILI_ELAS",
/ "DIS_CHOC",
/ "DIS_CONTACT",
/ "DIS_ECRO_CINE",
/ "DIS_GOUJ2E_ELAS",
/ "DIS_GOUJ2E_PLAS",
/ "DIS_GRICRA",
/ "CHOC_ENDO",
/ "CHOC_ELAS_TRAC",
/ "CHOC_ENDO_PENA",
/ "JONC_ENDO_PLAS",
/ "DIS_VISC",
/ "DIS_ECRO_TRAC",
/ "DRUCK_PRAGER",
/ "DRUCK_PRAG_N_A",
/ "ELAS_MEMBRANE_NH",
/ "ELAS_MEMBRANE_SV",
/ "ENDO_PORO_BETON",
/ "ENDO_CARRE",
/ "ENDO_FISS_EXP",
/ "ENDO_HETEROGENE",
/ "ENDO_ISOT_BETON",
/ "ENDO_LOCA_EXP",
/ "ENDO_LOCA_TC",
/ "ENDO_ORTH_BETON",
/ "ENDO_SCALAIRE",
/ "FLUA_PORO_BETON",
/ "FLUA_ENDO_PORO",
/ "FONDATION",
/ "GLRC_DAMAGE",
/ "GLRC_DM",
/ "GTN",
/ "DHRC",
/ "BETON_GRANGER",
/ "BETON_GRANGER_V",
/ "GRAN_IRRA_LOG",
/ "GRILLE_CINE_LINE",
/ "GRILLE_ISOT_LINE",
/ "HAYHURST",
/ "HOEK_BROWN",
/ "HOEK_BROWN_EFF",
/ "HOEK_BROWN_TOT",
/ "HUJEUX",
/ "IRRAD3M",
/ "JOINT_BA",
/ "JOINT_BANDIS",
/ "JOINT_MECA_RUPT",
/ "JOINT_MECA_FROT",
/ "JOINT_MECA_ENDO",
/ "KICHENIN_NL",
/ "KIT_CG",
/ "KIT_DDI",
/ "KIT_HH",
/ "KIT_H",
/ "KIT_HHM",
/ "KIT_HM",
/ "KIT_THH",
/ "KIT_THHM",
/ "KIT_THM",
/ "KIT_THV",
/ "KIT_THH2M",
/ "KIT_HH2M",
/ "KIT_HH2",
/ "KIT_THH2",
/ "LAIGLE",
/ "LEMAITRE",
/ "LEMAITRE_IRRA",
/ "LEMA_SEUIL",
/ "LETK",
/ "LKR",
/ "VISC_MAXWELL",
/ "VISC_MAXWELL_MT",
/ "MAZARS",
/ "MAZARS_UNIL",
/ "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",
/ "MOHR_COULOMB",
/ "RANKINE",
/ "RGI_BETON",
/ "RGI_BETON_BA",
/ "MONOCRISTAL",
/ "MULTIFIBRE",
/ "NORTON",
/ "NORTON_HOFF",
/ "POLYCRISTAL",
/ "RELAX_ACIER",
/ "ROUSSELIER",
/ "ROUSS_PR",
/ "ROUSS_VISC",
/ "RUPT_FRAG",
/ "SANS",
/ "VENDOCHAB",
/ "VISC_ENDO_LEMA",
/ "VISCOCHAB",
/ "VISC_CIN1_CHAB",
/ "VISC_CIN2_CHAB",
/ "VISC_CIN2_MEMO",
/ "VISC_CIN2_NRAD",
/ "VISC_MEMO_NRAD",
/ "VISC_DRUC_PRAG",
/ "VISC_GTN",
/ "VISC_IRRA_LOG",
/ "VISC_ISOT_LINE",
/ "VISC_ISOT_NL",
/ "VISC_ISOT_TRAC",
/ "VISC_TAHERI",
/ "VMIS_ASYM_LINE",
/ "VMIS_CIN1_CHAB",
/ "VMIS_CIN2_CHAB",
/ "VMIS_CINE_GC",
/ "VMIS_CINE_LINE",
/ "VMIS_CIN2_MEMO",
/ "VMIS_CIN2_NRAD",
/ "VMIS_MEMO_NRAD",
/ "VMIS_CINE_LINE",
/ "VMIS_ECMI_LINE",
/ "VMIS_ECMI_TRAC",
/ "VMIS_ISOT_LINE",
/ "VMIS_ISOT_NL",
/ "VMIS_ISOT_PUIS",
/ "VMIS_ISOT_TRAC",
/ "VMIS_JOHN_COOK",
/ "UMAT",
/ "MFRONT",
/ "BETON_BURGER",
/ "Iwan",
/ "META_LEMA_ANI",
/ "MetaAcierEPIL_PT",
/ "MohrCoulombAS",
/ "GonfElas",
/ "NLH_CSRM",
/ "MCC",
/ "CSSM",
◇ 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),
/ "SANS",
/ "MESSAGE",
◇ ALGO_CPLAN = / "DEBORST" (par défaut),
/ "ANALYTIQUE",
◇ SYME_MATR_TANG = / "OUI" (par défaut),
/ "NON",
# Si: not equal_to('RELATION', 'MFRONT')
# Si: equal_to("RELATION", 'KIT_DDI')
◆ RELATION_KIT = / "VMIS_CINE_LINE",
/ "VMIS_ISOT_TRAC",
/ "VMIS_ISOT_LINE",
/ "VMIS_ISOT_PUIS",
/ "GLRC_DM",
/ "BETON_GRANGER",
/ "BETON_GRANGER_V",
/ "BETON_UMLV",
/ "ROUSS_PR",
/ "BETON_DOUBLE_DP",
/ "ENDO_ISOT_BETON",
/ "MAZARS",
# Si: equal_to("RELATION", 'KIT_CG')
◆ RELATION_KIT = / "CABLE_GAINE_FROT",
/ "VMIS_ISOT_LINE",
/ "VMIS_ISOT_TRAC",
/ "VMIS_CINE_LINE",
/ "ELAS",
/ "SANS",
# 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",
/ "CZM_EXP_REG",
/ "CZM_LIN_REG",
/ "DRUCK_PRAGER",
/ "DRUCK_PRAG_N_A",
/ "ELAS",
/ "ENDO_ISOT_BETON",
/ "GonfElas",
/ "HOEK_BROWN_EFF",
/ "HOEK_BROWN_TOT",
/ "HUJEUX",
/ "Iwan",
/ "JOINT_BANDIS",
/ "NLH_CSRM",
/ "LAIGLE",
/ "LETK",
/ "LKR",
/ "MAZARS",
/ "MFRONT",
/ "MOHR_COULOMB",
/ "RANKINE",
/ "VISC_DRUC_PRAG",
/ "VISC_MAXWELL",
/ "VISC_MAXWELL_MT",
/ "GAZ",
/ "LIQU_SATU",
/ "LIQU_GAZ_ATM",
/ "LIQU_VAPE_GAZ",
/ "LIQU_AD_GAZ_VAPE",
/ "LIQU_AD_GAZ",
/ "LIQU_VAPE",
/ "LIQU_GAZ",
/ "HYDR_UTIL",
/ "HYDR_TABBAL",
/ "HYDR_VGM",
/ "HYDR_VGC",
/ "HYDR_ENDO",
# Si: 'MFRONT' in value("RELATION_KIT")
◆ COMPOR_MFRONT = compor_mgis,
◇ VERI_BORNE = / "ARRET" (par défaut),
/ "SANS",
/ "MESSAGE",
◇ ALGO_CPLAN = / "DEBORST" (par défaut),
/ "ANALYTIQUE",
◇ RESI_INTE = float (défaut: 1e-08),
◇ ITER_INTE_MAXI = int,
◇ SYME_MATR_TANG = / "OUI" (par défaut),
/ "NON",
# 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 = / "PETIT" (par défaut),
/ "PETIT_REAC",
/ "GROT_GDEP",
/ "SIMO_MIEHE",
/ "GDEF_LOG",
/ "GREEN_LAGRANGE",
◇ 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,
◇ ALGO_INTE = / "ANALYTIQUE",
/ "SECANTE",
/ "DEKKER",
/ "NEWTON_1D",
/ "BRENT",
/ "NEWTON",
/ "NEWTON_RELI",
/ "NEWTON_PERT",
/ "RUNGE_KUTTA",
/ "SPECIFIQUE",
/ "SEMI_EXPLICITE",
/ "BASCULE_EXPLICITE",
/ "SANS_OBJET",
# 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,
),
)
◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles
Opérandes#
La commande permet le calcul du complément de forces internes dû aux non linéarités de comportement. Ce calcul concerne les lois de comportement où les forces internes dépendent essentiellement des champs cinématiques (déplacements, vitesses, accélérations); elle concerne donc surtout les non linéarités localisées comme par exemple le contact pénalisé entre éléments discrets. Les termes calculés aux nœuds sont obtenus par la différence entre le résidu intégré de forces internes (dépendant à la fois des champs cinématiques et des paramètres de loi de comportement) et la force interne qui serait obtenue à partir des mêmes champs cinématiques en considérant un comportement linéaire pour la structure.
Remarque 1 :
Lenom du concept résultat créé de type dyna_transest obligatoirement différent du nom de concept résultat de type dyna_transutilisé en entrée sous l’opérande RESULTATet constitue une évolution de champs de type ‘ DEPL’.Ce nouveau résultat de forces nodales complémentaires s’utilise sous le mot-clé EXCIT_RESU dans la commande DYNA_VIBRA avec TYPE_RESU=”HARM”(dans ce cas après transformation par l’opérateur REST_SPEC_TEMP [U4.63.34]) et DYNA_LINE_TRAN ou DYNA_VIBRA et TYPE_RESU=”TRAN”. Un exemple est fourni dans le cas test SDLS119A.
Remarque 2 :
Dans le cas précis du calcul avec la relation ‘DIS_CHOC’affectée à des éléments discrets, Il est alors nécessaire de toujours orienter l’élément discret de contact dans la direction structure vers butée. Que ce soit avec des éléments POI1, comme toujours par le mot clé ORIENTATIONde AFFE_CARA_ELEM, mais aussi avec des éléments SEG2où on doit bien définir l’élément dans le sens structure candidate au contact, associée au nœud 1, vers la butée, associée au nœud 2.
Opérandes TOUT_ORDRE / NUME_ORDRE / LIST_ORDRE / NUME_MODE/ NOEUD_CMP / NOM_CAS / INST / LIST_INST / FREQ / LIST_FREQ / PRECISION / CRITERE#
Voir [U4.71.00] pour la description de ces opérandes.
Opérande MODELE#
♦ MODELE = mo,
Nom du modèle, nécessaire à entrer
Opérande CHAM_MATER#
◊ CHAM_MATER = chmater,
Nom du champ de matériau où sont définies les caractéristiques de matériau des éléments. Cet argument est nécessaire pour le calcul car les lois de comportement définies dans les mots clés COMPORTEMENT nécessitent toujours un champ de matériau.
Opérande CARA_ELEM#
◊ CARA_ELEM = carac,
Le concept des caractéristiques élémentaires de type cara_elem est nécessaire pour le calcul s’il existe dans le modèle des éléments de structure.
Opérande COMPORTEMENT#
◊ COMPORTEMENT = _F
Mot clé facteur permettant d’affecter des lois de comportement à des mailles ou des groupes de mailles du maillage.
On définit une affectation par occurrence du mot clé COMPORTEMENT.
Opérande RELATION#
◊ RELATION,
Nom d’une relation de loi de comportement sous format texte. Les relations admises sont celles où les forces internes dépendent essentiellement des champs cinématiques (déplacements, vitesses, accélérations); outre la relation ‘ELAS’, la première relation prévue est la relation ‘DIS_CHOC’ affectée à des éléments discrets.
Opérande DEFORMATION#
◊ DEFORMATION,
Nom du type de déformation sous format texte compatible avec le type des éléments affectés par la loi de comportement.
Opérandes TOUT/GROUP_MA/MAILLE#
◊ TOUT = “OUI” ,
La relation est affectée sur tout le maillage.
◊ GROUP_MA = lgma ,
La relation est affectée sur les groupes de mailles contenus dans la liste lgma.
◊ MAILLE = lma ,
La relation est affectée sur les mailles contenues dans la liste lma.