3. B modeling#

3.1. Characteristics of modeling#

_images/100000000000045F0000030353DEDBF695FE6177.png

The field \(\theta\) is calculated then the energy release rate \(G\), the stress coefficients \({K}_{1}\) and \({K}_{2}\), the energy recovery rate obtained by the formula d’IRWIN, the direction of propagation of the crack, are calculated.

3.2. Characteristics of the mesh#

Number of knots: 673

Number of meshes and types: 112 QUAD8 meshes and 142 TRIA6 meshes

3.3. Tested sizes and results#

The values tested are the energy recovery rate calculated by the theta method and the energy recovery rate calculated by the formula IRWIN using the stress intensity coefficients \({K}_{1}\) and \({K}_{2}\).

Identification

Reference

Tolerance

Crowns 1 to 6 \(G\)

\(2.3093{10}^{\mathrm{-}3}\)

< 1%

Crowns 1 to 6 \(G\) (IRWIN)

\(2.3093{10}^{\mathrm{-}3}\)

< 1%

Crowns 1 to 6 \({K}_{1}\)

\(24.491\)

< 1%

Crowns 1 to 6 \({K}_{2}\)

\(0.\)

absolute

3.4. note#

The calculation of \(G\), \({K}_{1}\), \({K}_{2}\), \(G\) (IRWIN) \(\mathrm{=}\frac{1}{E}({K}_{1}^{2}+{K}_{2}^{2})\) was carried out from 6 different \(\theta\) fields, each corresponding to a circular crown centered in \(C\).