1. Problème de référence#

1.1. Géométrie#

Selon la modélisation \(\mathrm{2D}\) ou \(\mathrm{3D}\) , on considère respectivement un carré ou un cube de côté 2 \(\text{mm}\) .

1.2. Propriétés du matériau#

Le matériau obéit à la loi de comportement élastique fragile ENDO_FISS_EXP à gradient d’endommagement (D_PLAN_GRAD_VARI et 3D_GRAD_VARI). Les données macroscopiques correspondent à:

\(E=30000\text{MPa}\)

Module d’Young

\(\nu =0.2\)

Coefficient de Poisson

\({G}_{f}=0.1{\text{N/mm}}^{2}\)

Énergie de fissuration

\(p=5\)

Paramètre de forme

\({f}_{t}=2.986\text{MPa}\)

Limite en traction

\({f}_{c}=29.86\text{MPa}\)

Limite en compression

\(D=50\text{mm}\)

Demi-largeur de la bande d’endommagement

Le choix de ces valeurs de ft et fc en lien avec le critère de François dans ENDO_FISS_EXP correspond à une limite en traction confinée de s c = 3 MPa. Par ailleurs, les paramètres internes du modèle qui interviennent dans la solution analytique sont égaux à:

\(k=1.5\times {10}^{-3}\text{MPa}\) ; \(m=11.111\) ; \(\gamma =9.534\times {10}^{3}\)

Dans le test, ils sont calculés par la macro-commande DEFI_MATER_GC.

1.3. Conditions aux limites et chargements#

Les déplacements sont imposés en tous les nœuds de la structure, de sorte à correspondre à la déformation homogène souhaitée. Plus précisément, le déplacement en un nœud de coordonnées \(X\) vaut : \(u(x)=\epsilon \cdot x\) . On commence par imposer une déformation uniaxiale de traction puis on décharge jusqu’à imposer une déformation de compression pour vérifier l’effet de la restauration de rigidité en compression.

1.4. Conditions initiales#

Aucune.