v1.03.127 MFRON02 - Test de l’interface Code_Aster-MFront : pour des lois de comportement de béton et de sols#

Résumé:

Ce test valide certains comportements de béton ou de métaux avec endommagement définis à l’aide de MFront par comparaison avec des comportement similaires de Code_Aster .

Modélisation C: cette modélisation permet de valider le modèle de fluage de béton de Burger avec vieillissement.

Modélisation G: cette modélisation permet de valider le modèle de Drucker-Prager.

Modélisation H: cette modélisation permet de valider l’appel à mtest sur un test similaire à la modélisation c.

Modélisation J: cette modélisation permet de valider l’utilisation d’Iwan (loi MFront) dans un KIT_THM

Chaque modélisation (sauf J) valide une loi de comportement, généralement par comparaison avec les résultats de la loi équivalente de Code_Aster .

Modélisation G#

Caractéristiques de la modélisation#

  • Comportement testé: DruckerPrager.mfront. Loi élasto-plastique associée de Drucker-Prager [cf. R7.01.16], le critère étant défini par: \(F(\sigma )={\sigma}_{\mathit{eq}}+a\mathit{tr}(\sigma )/3–k\le 0\)

Les coefficients a et k sont définis à partir des propriétés matériau de la façon suivante:

\(a=-3\frac{\tau –1}{\tau +1}\) et \(k=2\frac{{\sigma}_{c}\tau }{\tau +1}\)

  • La modélisation est équivalente à celle d’un fichier mtest (point matériel) de Mfront:

un point matériel est soumis à une déformation imposée \({\epsilon}_{xx}=5{10}^{-3}t\)

  • les propriétés matériau sont: \(E=200\mathit{GPa},\nu =0,3,\tau =0.6,{\sigma}_{c}=150\mathit{MPa}\) ,

Grandeurs testées et résultats#

La solution de référence est celle des tests de référence Mfront:

Identification (t=1)

Référence

Tolérance

\({\sigma}_{xx}(\mathit{MPa})\)

90

0,1%

Modélisation H#

Caractéristiques de la modélisation#

Cette modélisation permet de valider l’appel à mtest. Les données de mtest (fichier mfron01h.22) correspondent à celle de la modélisation C.

Grandeurs testées et résultats#

On teste ici simplement le code retour de mtest: il vaut 0 si l’exécution s’est bien passée.

Modélisation i#

Caractéristiques de la modélisation#

Cette modélisation permet de valider RELATION_KIT avec une loi Mfront. La mise en données a été pris dans le cas test wtnp114a où trois lois (ELAS, LIQU_SATU, HYDR_UTIL) sont combinées par RELATION_KIT. La loi élastique a été remplacée par une loi plastique Mfront.

Grandeurs testées et résultats#

On teste ici simplement des valeurs non régression.

Modélisation J#

Caractéristiques de la modélisation#

Comportement testé: Iwan

Les propriétés élastiques isotropes du matériau sont :

  • \(E=5.8\mathit{GPa}\)

  • \(\nu =0.3\)

  • \(\rho =2764\mathit{kg}.{m}^{-3}\)

Les paramètres du modèle d’Iwan (modélisation B) sont:

  • \({\gamma}_{\mathrm{ref}}=2.00e-4\)

  • \(\mathrm{n}=0.1\)

Grandeurs testées et résultats#

Les valeurs testées sont des valeurs de non-régression car ce test ne fait que vérifeir la possiibltié de brancher une loi MFront dans KIT_THM.

Localisation

Grandeur

\(A\)

\({\varepsilon}_{xx}\)

\({\varepsilon}_{yy}\)

\(B\)

\({\varepsilon}_{xx}\)

\({\varepsilon}_{yy}\)

\(C\)

\({u}_{x}\)

\({u}_{y}\)

\({\varepsilon}_{xx}\)

\({\varepsilon}_{yy}\)

D

\({\varepsilon}_{xx}\)

\({\varepsilon}_{yy}\)

Synthèse des résultats#

Les résultats sont satisfaisants et valident l’interface entre Code_Aster et MFRONT en 3D, axisymétrique et contraintes planes, pour des comportements de béton ou de sols.