2. Benchmark solution#
For each modeling, this test carries out an inter-comparison between the reference solution (obtained with a very fine time step), the solution with a moderately coarse discretization, the solution with the effect of temperature (or another control variable), the solution by changing the system of units (\(\mathrm{Pa}\) in \(\mathrm{MPa}\)), and the solution obtained after rotation or symmetry.
2.1. Definition of robustness test cases#
We propose 3 analysis angles to test the robustness of the integration of laws of behavior:
studies of equivalent problems
tangent matrix check
study of the discretization of the time step
For each of them, we study the evolution of the relative differences between several calculations using the same law but presenting different parameters or calculation options. The exploitation focuses on the invariants of the stress tensor: tensor trace, Von-Mises constraint and scalar internal variables: generally it is the cumulative plasticity.
The global convergence criteria are the values provided by default by ASTER. (RESI_GLOB_RELA =10-6, ITER_GLOB_MAXI =10). We adopted a usual Newton scheme for updating the tangent matrix: at each converged increment (REAC_INC =1) and every 1 iterations (REAC_ITER =1).
2.2. Equivalent problem studies#
For a rough discretization of the paths: 1 time step for each segment of the journey, the solution obtained for each law is compared to 3 strictly equivalent problems for the state of the material point:
\(\mathit{Tpa}\), same path with a change of unit, we substitute \(\mathrm{Pa}\) for \(\mathrm{MPa}\) in the material data and the possible parameters of the law,
\(\mathit{Trot}\), path by imposing the same \(\stackrel{ˉ}{\varepsilon }\) tensor after a rotation: \({\text{}}^{t}R\cdot \stackrel{ˉ}{\varepsilon }\cdot R\) where \(R\) is a rotation matrix. For the 2D case, the angle of rotation will be \(\mathrm{\alpha }=0.9\mathit{rad}\), for the 3D configuration, we chose the Euler angles with the arbitrary values {\(\mathrm{\Psi }=0.9\mathit{rad}\), \(\mathrm{\theta }=0.7\mathit{rad}\), and \(\mathrm{\phi }=0.4\mathit{rad}\)},
\(\mathit{Tsym}\), path by imposing the tensor \(\stackrel{ˉ}{\varepsilon }\) after symmetry: permutation of the axes \(x\) and \(y\) in 2D, permutation of \(x\) in \(y\), \(y\) in \(z\) and \(z\) in and in \(x\) in in 3D.
For each of these problems, the solution (stress invariants, cumulative equivalent plastic deformation) must be identical to the basic solution, obtained with the same time discretization. The reference value of the difference is therefore 0. In practice, this means that the difference found must be in the order of machine precision, i.e. around \(1.E-15\).
2.3. Tangent matrix test#
For each behavior, the tangent matrix is also tested, by difference with the matrix obtained by disturbance. Again, the reference value is 0.
2.4. Study of the discretization of the time step (A2)#
We study the behavior of law integration as a function of discretization. For the same modeling, and therefore a given behavior, several discretizations are studied here at different times, by multiplying the number of steps in the loading path by 5. In the reference [1], discretization is extended up to 3125 increments per segment on the same principle. Here, in order to limit the duration of the tests, we limit ourselves to 3 successive refinements. This leads to the following discretization:
Number of intervals per load segment |
1 |
5 |
25 |
Total number of steps over the entire trip |
8 |
40 |
200 |
Calculation |
\(\mathit{T0}\) |
|
|
The reference solution, \(\mathit{Tréf}\), is the one obtained for \(N=25\), i.e. 200 steps for the entire journey. These different solutions make it possible to judge the sensitivity to large time steps and the robustness of the integration. To show the convergence speed as a function of the time step, we report here the solutions presented in [1], up to 3125 time steps for each of the 8 segments of the loading path.
2.4.1. VMIS_LINE#
Gaps |
\(\mathrm{N1}\) |
|
|
|
|
|
||
\({\mathrm{V1}}_{N}\) |
3.70e-02 |
1.38e-02 |
1.38e-02 |
3.37e-03 |
6.82e-04 |
1.14e-04 |
0.00e+00 |
|
\(\mathrm{VMIS}\) |
4.34e-03 |
1.86e-03 |
1.86e-03 |
4.72e-04 |
9.72e-05 |
1.64e-05 |
0.00e+00 |
|
\(\mathrm{TRAC}\) |
1.19e-01 |
6.89e-02 |
1.70e-02 |
1.70e-02 |
3.45e-03 |
5.80e-04 |
0.00e+00 |
2.4.2. VMIS_ISOT_TRAC#
Gaps |
\(\mathrm{N1}\) |
|
|
|
|
|
||
\({\mathrm{V1}}_{N}\) |
3.58e-02 |
1.34e-02 |
1.34e-02 |
3.26e-03 |
6.60e-04 |
1.11e-04 |
0.00e+00 |
|
\(\mathrm{VMIS}\) |
6.38e-03 |
2.38e-03 |
2.38e-03 |
5.81e-04 |
1.18e-04 |
2.00e-05 |
0.00e+00 |
|
\(\mathrm{TRAC}\) |
1.20e-01 |
7.69e-02 |
1.67e-02 |
1.67e-02 |
3.30e-03 |
5.53e-04 |
0.00e+00 |
2.4.3. VMIS_CINE_LINE#
Gaps |
\(\mathrm{N1}\) |
|
|
|
|
|
||
\(\mathrm{VMIS}\) |
3.91e-03 |
1.05e-03 |
1.05e-03 |
2.16e-04 |
4.15e-05 |
6.91e-06 |
0.00e+00 |
|
\(\mathrm{TRAC}\) |
7.48e-14 |
7.44e-14 |
7.44e-14 |
7.44e-14 |
7.69e-14 |
5.87e-14 |
0.00e+00 |
2.4.4. VMIS_ECMI_LINE#
Gaps |
\(\mathrm{N1}\) |
|
|
|
|
|
||
\({\mathrm{V1}}_{N}\) |
3.71e-02 |
1.39e-02 |
1.39e-02 |
3.40e-03 |
6.88e-04 |
1.15e-04 |
0.00e+00 |
|
\(\mathrm{VMIS}\) |
3.63e-03 |
2.00e-03 |
2.00e-03 |
4.79e-04 |
9.53e-05 |
1.61e-05 |
0.00e+00 |
|
\(\mathrm{TRAC}\) |
1.64e-01 |
9.16e-02 |
2.18e-02 |
2.18e-02 |
4.33e-03 |
7.30e-04 |
0.00e+00 |
2.4.5. VMIS_ECMI_TRAC#
Gaps |
\(\mathrm{N1}\) |
|
|
|
|
|
||
\({\mathrm{V1}}_{N}\) |
3.70e-02 |
1.38e-02 |
1.38e-02 |
3.38e-03 |
6.84e-04 |
1.15e-04 |
0.00e+00 |
|
\(\mathrm{VMIS}\) |
2.36e-03 |
1.18e-03 |
1.18e-03 |
2.98e-04 |
6.00e-05 |
1.01e-05 |
0.00e+00 |
|
\(\mathrm{TRAC}\) |
1.46e-01 |
8.51e-02 |
2.13e-02 |
2.13e-02 |
4.31e-03 |
7.22e-04 |
0.00e+00 |
2.4.6. VMIS_CIN1_CHAB#
Gaps |
\(\mathrm{N1}\) |
|
|
|
|
|
||
\({\mathrm{V1}}_{N}\) |
3.32e-02 |
1.12e-02 |
1.12e-02 |
2.57e-03 |
5.10e-04 |
8.52e-05 |
0.00e+00 |
|
\(\mathrm{VMIS}\) |
9.04e-02 |
3.24e-02 |
3.24e-02 |
7.45e-03 |
1.48e-03 |
2.49e-04 |
0.00e+00 |
|
\(\mathrm{TRAC}\) |
3.34e-14 |
3.31e-14 |
3.31e-14 |
3.27e-14 |
3.48e-14 |
3.86e-14 |
0.00e+00 |
2.4.7. VMIS_CIN2_CHAB#
Gaps |
\(\mathit{N1}\) |
|
|
|
|
|
||
\({\mathit{V1}}_{N}\) |
3.32e-02 |
1.12e-02 |
1.12e-02 |
2.57e-03 |
5.10e-04 |
8.52e-05 |
0.00e+00 |
|
\(\mathit{VMIS}\) |
9.04e-02 |
3.24e-02 |
3.24e-02 |
7.45e-03 |
1.48e-03 |
2.49e-04 |
0.00e+00 |
|
\(\mathit{TRAC}\) |
3.72e-14 |
3.69e-14 |
3.69e-14 |
3.69e-14 |
3.83e-14 |
4.75e-14 |
0.00e+00 |
2.4.8. MOHR_COULOMB#
The internal variables have the following meanings:
\(\mathit{V1}\): plastic volume deformation \({\epsilon }_{v}^{p}=\frac{1}{3}\mathit{trace}\left({\epsilon }^{p}\right)\)
\(\mathit{V2}\): plastic volume deformation \(\mid {\epsilon }_{d}^{p}\mid =\sqrt{\frac{2}{3}\left({\epsilon }^{p}-{\epsilon }_{v}^{p}I\right)\mathrm{:}\left({\epsilon }^{p}-{\epsilon }_{v}^{p}I\right)}\)
\(\mathit{V3}\): plasticity indicator
Gap \(\left(\text{\%}\right)\) |
|
|
|
|
|
|
||
\(V1\) |
1.34 |
0.57 |
0.57 |
0.57 |
0.57 |
0.57 |
0.6 |
0 |
\(V2\) |
1.34 |
0.57 |
0.57 |
0.57 |
0.57 |
0.57 |
0.6 |
0 |
\(V3\) |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
|
\(\mathit{VMIS}\) |
6.41 |
2.34 |
2.34 |
0.34 |
0.57 |
0.1 |
0.02 |
0 |
\(\mathit{TRAC}\) |
2.06 |
0.87 |
0.87 |
0.24 |
0.05 |
0.009 |
0 |
2.4.9. Mohr Coulombas#
The internal variables tested are as follows:
\(V7\): plastic volume deformation \(\mid {\epsilon }_{d}^{p}\mid =\sqrt{\frac{2}{3}\left({\epsilon }^{p}-{\epsilon }_{v}^{p}I\right)\mathrm{:}\left({\epsilon }^{p}-{\epsilon }_{v}^{p}I\right)}\)
\(V8\): plastic volume deformation \({\epsilon }_{v}^{p}=\frac{1}{3}\mathit{trace}\left({\epsilon }^{p}\right)\)
\(V9\): dissipation.
Gap \(\left(\text{\%}\right)\) |
|
|
|
|
|
|
||
\(V7\) |
1.54 |
0.62 |
0.62 |
0.17 |
0.037 |
0.0026 |
0 |
|
\(V8\) |
1.54 |
0.62 |
0.62 |
0.17 |
0.037 |
0.0026 |
0 |
|
\(V9\) |
0.23 |
0.20 |
0.20 |
0.088 |
0.021 |
0.0053 |
0 |
|
\(\mathrm{VMIS}\) |
2.95 |
0.31 |
0.31 |
0.076 |
0.013 |
0.0015 |
0 |
|
\(\mathrm{TRAC}\) |
1.54 |
0.62 |
0.62 |
0.17 |
0.037 |
0.0026 |
0 |
2.5. REFERENCES#
LEVASSEUR: « Third-party Application Maintenance of the _Aster code » Verification of the robustness and reliability of the integration of behavioral laws in ASTER. Report PRINCIPIA RET .693.127.01 December 2006.