u4.81.41 Opérateur POST_BORDET#

Syntaxe#

Détail de la syntaxe
table = POST_BORDET(
    ◆ / TOUT = "OUI" (ou non renseigné),
      / GROUP_MA = grma,
    ◆ / INST = float,
      / NUME_ORDRE = int,
    ◇ PRECISION = float (défaut: 1e-06),
    ◇ CRITERE = / "RELATIF",
                / "ABSOLU" (par défaut),
    ◇ PROBA_NUCL = / "NON" (par défaut),
                   / "OUI",
    # Si: equal_to("PROBA_NUCL", 'OUI')
        ◆ PARAM = _F(
             ◆ M = float,
             ◆ SIGM_REFE = fonction,
             ◆ VOLU_REFE = float,
             ◆ SIG_CRIT = float,
             ◆ SEUIL_REFE = float,
             ◆ SEUIL_CALC = nappe / fonction,
             ◆ DEF_PLAS_REFE = float,
          ),
    # Si: equal_to("PROBA_NUCL", 'NON')
        ◆ PARAM = _F(
             ◆ M = float,
             ◆ SIGM_REFE = fonction,
             ◆ VOLU_REFE = float,
             ◆ SIG_CRIT = float,
             ◆ SEUIL_REFE = float,
             ◆ SEUIL_CALC = nappe / fonction,
          ),
    ◆ RESULTAT = resultat,
    ◆ TEMP = fonction,
    ◇ COEF_MULT = float (défaut: 1.0),
)


◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles

Opérandes#

Opérande RESULTAT#

Désigne le résultat du calcul thermo-mécanique pour lequel on calcule les grandeurs de Bordet. La structure de données RESULTAT fournie doit comporter un et un seul modèle et un et un seul champ matériau.

Opérande TOUT/GROUP_MA#

♦ /TOUT = “OUI”,

/GROUP_MA = group_ma, [l_group_ma]

Désigne le domaine du modèle sur lequel le calcul des grandeurs de Bordet sera effectué. Par défaut, il est réalisé sur l’ensemble du modèle. Notons que comme le modèle fait intervenir la variation de la déformation plastique, le groupe de mailles sur lequel le calcul est effectué doit contenir a minima tout le domaine sur lequel la déformation plastique a évolué durant le chargement, et que prendre un domaine plus grand ne changera pas le résultat (mais la durée du calcul).

Opérande INST/NUME_ORDRE#

♦ /INST = instant, [R]

/NUME_ORDRE = ordre [I]

Le calcul sera réalisé pour tous les instants ou numéro d’ordre jusqu’à l’instant ou numéro ordre désigné par cet opérande.

Le résultat sera une table contenant les différents instants de calcul, la contrainte et la probabilité de Bordet.

Opérande PRECISION#

◊ /PRECISION = 1E-06 [DEFAUT]

prec, [R]

Permet de définir, dans le cas de l’utilisation du mot clé INST, une précision dans la recherche du dernier instant de calcul des grandeurs de Bordet. Si la précision n’est pas définie par l’utilisateur, une précision de 1E-06 est appliquée.

Si l’instant demandé par l’utilisateur ne correspond à aucun intervalle [instant-prec ; instant+prec]sauvegardé dans le résultat, le calcul s’arrête en erreur fatale.

Opérande CRITERE#

◊ /CRITERE = “ABSOLU” [DEFAUT]

“RELATIF”,

Permet de définir, dans le cas de l’utilisation des mots clés INST et PRECISION ,si la précision est absolue ou relative. Par défaut, elle est absolue.

Opérande PROBA_NUCL#

◊ PROBA_NUCL = /”NON [DEFAUT]

/”OUI”

Indique si l’utilisateur souhaite prendre en compte le terme exponentiel dans son calcul. Ce terme, comme précisé plus haut, est d’influence négligeable si \(\frac{{\sigma}_{\mathrm{ys},0}{\varepsilon}_{p,0}}{{\sigma}_{\mathrm{ys}}(T,{\dot{\varepsilon}}_{p})}\gg {\varepsilon}_{p}\) . Si l’utilisateur souhaite prendre en compte ce terme, il doit alors renseigner un paramètre supplémentaire qui est la déformation plastique équivalente de référence. Sinon, ce paramètre est inutile.

Opérande PARAM#

♦ PARAM = _F(

♦ M = m, [R]

♦ SIGM_REFE = sigma_ref, [fonction]

♦ VOLU_REFE = r_sup, [R]

♦ SIG_CRIT = sigma_crit, [R]

♦ SEUIL_REFE = lim_elas_ref, [R]

♦ SEUIL_CALC = lim_elas, [fonction,nappe]

♦ DEF_PLAS_REFE = lim_elas, [R]

),

Désigne l’ensemble des paramètres matériau nécessaires au calcul des grandeurs du modèle de Bordet. Les trois premiers ont des équivalents dans le modèle de Beremin (avec des valeurs qui peuvent varier d’un modèle à l’autre), les quatre suivants étant spécifiques au modèle de Bordet.

Mot clé M#

♦ M = m, [R]

Désigne l’exposant \(m\) de la loi de type Weibull (parfois appelé facteur de forme).

Attention : cette grandeur n’est pas forcément égale à son équivalent dans le modèle de Beremin.

Mot clé SIGM_REFE#

♦ SIGM_REFE = sigma_ref, [fonction]

Désigne la contrainte de référence \({\sigma}_{u}(T)\) de la loi de type de Weibull (parfois appelé facteur d’échelle). C’est la contrainte pour laquelle la probabilité de rupture cumulée des sites potentiels de clivage vaut 1.

Cette contrainte dépend de la température ; on attend ici une fonction de la température.

Attention : cette grandeur n’est pas forcément égale à son équivalent dans le modèle de Beremin.

Mot clé VOLU_REFE#

♦ VOLU_REFE = V0 [R]

Désigne le volume élémentaire de référence \({V}_{0}\) de la zone plastique.

Attention : cette grandeur n’est pas forcément égale à son équivalent dans le modèle de Beremin.

Mot clé SIG_CRIT#

♦ SIG_CRIT = sigma_crit, [R]

Désigne la contrainte critique \({\sigma}_{\mathrm{th}}\) en dessous de laquelle la propagation des microfissures ferritiques ne peut être significative ; si \({\sigma}_{1}<{\sigma}_{\mathrm{th}}\) en tout point, alors nécessairement \({P}_{\mathrm{Bordet}}=0\) .

Mot clé SEUIL_REFE#

♦ SEUIL_REFE = lim_elas_ref, [R]

Désigne la limite élastique \({\sigma}_{\mathrm{ys},0}\) à une température de référence à utiliser dans le modèle.

Mot clé SEUIL_CALC#

♦ SEUIL_CALC = lim_elas, [fonction,nappe]

Désigne la limite élastique du matériau \({\sigma}_{\mathrm{ys}}(T,{\dot{\varepsilon}}_{p})\) , qui en toute rigueur dépend de la température et de la vitesse de déformation plastique.

Si on ne connaît pas la dépendance à la vitesse de déformation plastique, on peut utiliser pour SEUIL_CALC une simple fonction de la température.

Si on connaît la dépendance à la fois à la température et à la vitesse de déformation plastique, on peut définir une nappe ; le paramètre de la nappe doit être la vitesse de déformation plastique, et la variable pour chaque fonction la température (cf. document utilisateur de DEFI_NAPPE)

Exemple :

SIGY1=DEFI_FONCTION(NOM_PARA=”TEMP”,

VALE=(0.,200.,100.,200.,),

PROL_DROITE=”CONSTANT”,PROL_GAUCHE=”CONSTANT”,);

SIGY2=DEFI_FONCTION(NOM_PARA=”TEMP”,

VALE=(0.,300.,100.,300.,),

PROL_DROITE=”CONSTANT”,PROL_GAUCHE=”CONSTANT”,);

SIGY=DEFI_NAPPE(NOM_PARA=”EPSI”,

PROL_DROITE=”CONSTANT”,PROL_GAUCHE=”CONSTANT”,

PARA=(0.0005,0.001),

FONCTION=(SIGY1,SIGY2),);

CALC_BORDET(…..SEUIL_CALC=SIGY,…);

Mot clé DEF_PLAS_REFE#

♦ DEF_PLAS_REFE = lim_elas, [R]

Ce mot clé est OBLIGATOIRE si PROBA_NUCL=”OUI” et INTERDIT si PROBA_NUCL=”NON”.

Il désigne la déformation plastique équivalente de référence \({\varepsilon}_{p,0}\) qui intervient uniquement dans le terme exponentiel.

Opérande TEMP#

♦ TEMP = temperature, [fonction]

Comme précisé dans les paragraphes ci-avants, certains paramètres matériau dépendent de la température. TEMP désigne la température, considérée uniforme pour l’instant sur la zone de calcul de Bordet. L’utilisateur peut renseigner une fonction constante, auquel cas la température est considérée uniforme en espace et constante en temps, ou une fonction du temps (paramètre INST), auquel cas la température est considérée uniforme en espace mais évolutive en temps.

Opérande COEF_MULT#

◊ COEF_MULT = 1, [DEFAUT]

coefficient, [R]

La valeur par défaut de ce coefficient est 1.0.

Le tableau suivant, dans lequel l’épaisseur est notée \(e\) , indique des valeurs typiques du coefficient \(C\) en fonction du type de symétrie:

  • symétrie simple : le plan de symétrie du maillage passe par le plan du défaut et le défaut est entièrement maillé,

  • symétrie double : le plan de symétrie du maillage passe également par le plan du défaut mais une seule moitié du défaut est maillé.

3D et 3D_SI

AXIS et AXIS_SI

D_PLAN et D_PLAN_SI

C_PLAN

SIMPLE

2

\(4\pi\)

2e

2e

DOUBLE

4

sans objet

sans objet

sans objet

NON

1

\(2\pi\)

e

e

Tableau 3.9-1 : Valeurs du coefficient multiplicateur de symétrie-épaisseur

Exemple d’utilisation#

On trouvera des exemples dans le cas test élémentaire zzzz268 et dans le cas test ssna108a.

Des conseils d’utilisation de ce modèle sont donnés dans la documentation [U2.05.08].