3. B modeling

3.1. Characteristics of the mesh

The geometry is identical to modeling A. The mesh represents the blank piece of any crack (it will be represented by the X- FEM method); the mesh is linear.

3.2. Calculations made

In order to test the various possibilities offered by the CALC_G_XFEM command, several calculation configurations are implemented.

The first is the reference result obtained by thermal loading without an initial state.

The others are derived from the calculation with initial state without thermal loading; they differ in the values of the radius of the crowns supplied to order CALC_G_XFEM. As for modeling A, stress intensity factors are also tested (the reference is obtained by POST_K1_K2_K3). Note that only the initial stress field at the gauss points was used, the definition of a stress field at the nodes not being available with X- FEM.

3.3. Tested sizes and results

3.3.1. Tested sizes

The value of the energy return rate resulting from the operator is tested. They should all be equal within a low tolerance.

3.3.2. Results

Testing energy return rates and stress intensity factors:

Identification	Reference	Reference type	Tolerance
CALC_G_XFEM
\(G\) from thermal calculation (CALC_G)		“NON_REGRESSION”	“”
\(G\) with automatic spokes (CALC_G)	55.36	“AUTRE_ASTER”	2.0%
\(G\) with R_ INF = 2.4, R_ SUP = 4.8 (CALC_G)	55.36	“AUTRE_ASTER”	2.0%
\(G\) with R_ INF =4,228, R_ SUP = 9 (CALC_G)	55.36	“AUTRE_ASTER”	3,5%
\({K}_{I}\) with automatic spokes (CALC_K_G)	3574	“AUTRE_ASTER”	1.0%
\({K}_{I}\) with R_ INF = 2.4 R_ SUP = 4.8 (CALC_K_G)	3574	“AUTRE_ASTER”	1.0%
\({K}_{I}\) with R_ INF =4,228, R_ SUP = 9 (CALC_K_G)	3574	“AUTRE_ASTER”	1,5%