4. Exemples d’utilisation

On trouvera des exemples dans les cas test zzzz163 et sslv101. Des conseils d’utilisation de ce modèle sont donnés dans la documentation [U2.05.08].

4.1. Exemple d’utilisation pour le modèle de Beremin

Dans le cas ci-dessous, la contrainte de Weibull est calculée sur le groupe de mailles « CT », en appliquant les options suivantes :

• Contrainte moyenne par élément

• Filtre en plasticité active, sur la base de la variable interne V2

• Maximum des contraintes sur l’historique (HIST_MAXI = « OUI » par défaut)

BERE = POST_BEREMIN(
    RESULTAT=U,
    COEF_MULT=1,
    DEFORMATION="PETIT",
    FILTRE_SIGM="SIGM_ELMOY",
    NUME_VARI=2,
    WEIBULL=_F(GROUP_MA="CT", M=8.0, VOLU_REFE=1.25e-4, SIGM_REFE=2630.0e0),
)

Dans l’exemple ci-dessous en grandes déformations logarithmiques, la contrainte de Weibull est calculée à partir des composantes du tenseur T stockées dans les variables internes 4 à 9 (à l’utilisateur d’adapter les numéros de variables internes en fonction de la loi de comportement utilisée par exemple)

BERE = POST_BEREMIN(
    RESULTAT=U,
    COEF_MULT=1,
    DEFORMATION="GDEF_LOG",
    LIST_NUME_SIEF=[4,5,6,7,8,9],
    FILTRE_SIGM="SIGM_ELMOY",
    NUME_VARI=2,
    WEIBULL=_F(GROUP_MA="CT", M=8.0, VOLU_REFE=1.25e-4, SIGM_REFE=2630.0e0),
)

4.2. Exemple d’utilisation pour la déclinaison industrielle de Beremin

Dans le cas de plans définis par des groupes de noeuds :

fiss = DEFI_FOND_FISS(
    MAILLAGE=mail_meca,
    FOND_FISS=_F(
        GROUP_MA="FRONT0",
        TYPE_FOND="OUVERT",
    ),
    LEVRE_SUP=_F(GROUP_MA="LEVRE"),
    SYME="OUI",
)


liste_plans = ["Plan_" + str(i) for i in range(0, 51)]

beremin_2D = POST_BEREMIN(
    COEF_MULT=2.0,
    DEFORMATION="PETIT",
    FILTRE_SIGM="SIGM_ELMOY",
    HIST_MAXI="NON",
    RESULTAT=evol_meca,
    WEIBULL=_F(
        GROUP_MA="TORE_COMPLET",
        M=3,
        SIGM_SEUIL=1300.0,
        VOLU_REFE=0.000125,
    ),
    METHODE_2D=_F(GROUP_NO_PLAN=liste_plans , FISSURE = fiss),
)

Dans l’exemple ci-dessus, on définit un front de fissure par le biais de la commande DEFI_FOND_FISS puis on applique le post-traitement de Beremin sur 51 plans définis par des groupes de noeuds nommés « Plan_i », avec i allant de 0 à 50.

Dans le cas d’un plan défini par un maillage 2D, pour lequel on demande par ailleurs d’imprimer le résultat de la projection dans un fichier .med :

mesh_2D_tore_c0 = LIRE_MAILLAGE(FORMAT="MED", UNITE=22)

beremin_2d_c0_tore = POST_BEREMIN(
    RESULTAT=resu_elas,
    COEF_MULT=2,
    DEFORMATION="PETIT",
    FILTRE_SIGM="SIGM_ELMOY",
    WEIBULL=_F(GROUP_MA="CT", M=18.0, VOLU_REFE=1.25e-4, SIGM_REFE=2630.0e0),
    METHODE_2D=_F(  MAILLAGE_PLAN=mesh_2D_tore_c0,
                    NOM_MAIL_MED="Couche_0_TORE",
                    UNITE_RESU = 80),
)