5. C modeling#
5.1. Characteristics of modeling#
Method MAILLAGE is used by PROPA_FISS. The CALC_G operator is used to calculate stress intensity factors.
5.2. Characteristics of the mesh#
We use the same mesh as that of modeling \(A\).
5.3. Tested sizes and results#
In this modeling, we do not use the COMP_LINE keyword from PROPA_FISS, but we give as input to PROPA_FISS a table containing several times. We therefore test by computer that the calculation of the cycle is correct.
The values of \({K}_{I}\) and \({K}_{\mathrm{II}}\) are tested for both cracks after each propagation. To check if these values are correct, a relative tolerance equal to \(\text{1\%}\) is used for the values of \({K}_{I}\). On the other hand, to check if the value of \({K}_{\mathrm{II}}\) is zero, we use an absolute tolerance (threshold value) linked to the value of \({K}_{I}\): it is considered that \({K}_{\mathrm{II}}\) is zero if its value is less than \(\text{0,1\%}\) of the value of \({K}_{I}\). Indeed, in this case we can overlook the value of \({K}_{\mathrm{II}}\).
Propagation |
Crack |
\({K}_{I}\) reference [\(\mathrm{Pa}\sqrt{\mathrm{mm}}\)] |
Tolerance |
1 |
left |
3.80E+007 |
< 1% |
righthand |
3.80E+007 |
< 1% |
|
2 |
left |
4.18E+007 |
< 1% |
righthand |
4.18E+007 |
< 1% |
|
3 |
left |
4.63E+007 |
< 1% |
righthand |
4.63E+007 |
< 1% |
Propagation |
Crack |
\({K}_{\mathrm{II}}\) reference [\(\mathrm{Pa}\sqrt{\mathrm{mm}}\)] |
Tolerance [\(\mathrm{Pa}\sqrt{\mathrm{mm}}\)] |
1 |
left |
0 |
< \({K}_{I\mathit{Réf}}\mathrm{/}1000\) |
righthand |
0 |
< \({K}_{I\mathit{Réf}}\mathrm{/}1000\) |
|
2 |
left |
0 |
< \({K}_{I\mathit{Réf}}\mathrm{/}1000\) |
righthand |
0 |
< \({K}_{I\mathit{Réf}}\mathrm{/}1000\) |
|
3 |
left |
0 |
< \({K}_{I\mathit{Réf}}\mathrm{/}1000\) |
righthand |
0 |
< \({K}_{I\mathit{Réf}}\mathrm{/}1000\) |
5.4. notes#
All values tested are within the tolerances used. This means that method MAILLAGE correctly calculates both the position of the two cracks and the level sets.
The error obtained on the values of \({K}_{I}\) is less than \(\text{1\%}\) and the values of \({K}_{\mathrm{II}}\) are always less than \({K}_{I}\mathrm{/}1000\). The results obtained are therefore very satisfactory.