4. C modeling#
4.1. Characteristics of modeling#
Modeling \(C\) is*tridimensional* and*static nonlinear*.
We first perform an elastic preconsolidation (\(\mathrm{ELAS}\)) of the sample up to \({p}_{c}=200\mathrm{kPa}\) (1st phase of the calculation). This preconsolidation takes place in \(1\), a time step between \(t=-2\) and \(t=0\).
The vertical displacement imposed on the upper facet varies between \(0\) and \(-0.2\) (2nd phase of the calculation) in \(100\) time steps between \(t=0\) and \(t=10\). During this second phase, the automatic subdivision of the time step is activated to manage the situations of non-convergence of local integration. In the integration of the equilibrium equations, an update of the tangent matrix is requested, which is provided by the routines of Hujeux’s law and significantly accelerates the convergence. We also require the subdivision of the time step (command DEFI_LIST_INST) to deal with situations where local integration fails due to excessively large load increments. This feature is highly recommended.
4.2. Tested sizes and results#
4.2.1. Tested values#
The solutions are calculated at point \(C\) and compared to references GEFDYN. They are given in terms of equivalent stress \(Q\), total volume deformation \({\epsilon }_{v}\) and isotropic work-hardening coefficients \(\left({r}_{\text{ela}}^{\text{iso},m}+{r}_{\text{iso}}^{m}\right)\) and deviatory \(\left({r}_{\text{ela}}^{d,m}+{r}_{\text{dev}}^{m}\right)\), and summarized in the following tables:
\(Q=\sqrt{\frac{1}{2}\underline{\underline{s}}\mathrm{:}\underline{\underline{s}}}[\mathit{Pa}]\)
\({\epsilon }_{\mathit{zz}}\) |
Code_Aster |
GEFDYN |
relative error |
-1% |
|
||
-2% |
|
||
-5% |
|
||
-10% |
|
||
-20% |
|
\({\epsilon }_{V}=\text{trace}\left(\epsilon \right)\)
\({\epsilon }_{\mathit{zz}}\) |
Code_Aster |
GEFDYN |
relative error |
-1% |
-7.389E-3 |
-7.47E-3 |
|
-2% |
-1.001E-2 |
-1.005E-2 |
|
-5% |
-1.229E-2 |
-1.227E-2 |
|
-10% |
-1.096E-2 |
-1.092E-2 |
|
-20% |
-4.88E-3 |
-4.88E-3 |
|
\(\left({r}_{\text{ela}}^{d,m}+{r}_{\text{dev}}^{m}\right)\)
\({\epsilon }_{\mathit{zz}}\) |
Code_Aster |
GEFDYN |
relative error |
-1% |
0.642 |
0.648 |
|
-2% |
0.761 |
0.765 |
|
-5% |
0.877 |
0.878 |
|
-10% |
0.931 |
0.932 |
|
-20% |
0.964 |
0.964 |
|
\(\left({r}_{\text{ela}}^{\text{iso},m}+{r}_{\text{iso}}^{m}\right)\)
\({\epsilon }_{\mathit{zz}}\) |
Code_Aster |
GEFDYN |
relative error |
-1% |
0.102 |
0.102 |
|
-2% |
0.107 |
0.108 |
|
-5% |
0.115 |
0.115 |
|
-10% |
0.125 |
0.126 |
|
-20% |
0.147 |
0.147 |
|
4.2.2. Comments#
The relative error is higher when the values tested are lower. In short, the difference between the two codes is very reasonable.