6. D modeling#

6.1. Characteristics of modeling#

This modeling is strictly identical to the previous one. In order to test that the HM-XFEM model is functional with PENTA15, we choose to subdivide the 5 HEXA20 from the previous modeling, each into 2 PENTA15. This makes a total of 10 quadratic pentahedral elements. As before, the pentahedral elements contained in extreme hexahedra are not enriched. On the contrario those contained in the 3 central hexahedra are enriched by the Heaviside function.

6.2. Characteristics of the mesh#

The mesh consists of 10 quadratic pentahedral cells (PENTA15).

6.3. Tested sizes and results#

The results (resolution with STAT_NON_LINE) are summarised in the tables below. To test all the nodes in the column at the same time, we calculate MIN and MAX.

Quantities tested

Reference type

Reference value

Tolerance (%)

DZ (below) MIN

“ANALYTIQUE”

-4.323558729E-03

0.0001

DZ (below) MAX

“ANALYTIQUE”

-4.323558729E-03

0.0001

DZ (above) MIN

“ANALYTIQUE”

4.297130927-03

0.0001

DZ (above) MAX

“ANALYTIQUE”

4.297130927-03

0.0001

Quantities tested

Reference type

Reference value

Tolerance (%)

DX (below) MIN

“ANALYTIQUE”

0

0.0001

DX (below) MAX

“ANALYTIQUE”

0

0.0001

DX (above) MIN

“ANALYTIQUE”

0

0.0001

DX (above) MAX

“ANALYTIQUE”

0

0.0001

Quantities tested

Reference type

Reference value

Tolerance (%)

DY (below) MIN

“ANALYTIQUE”

0

0.0001

DY (below) MAX

“ANALYTIQUE”

0

0.0001

DY (above) MIN

“ANALYTIQUE”

0

0.0001

DY (above) MAX

“ANALYTIQUE”

0

0.0001

_images/10000000000004B2000002F35083EC4F3879FAC7.png

Figure 6.3-a : Field of movement by direction (Oz)

A second modeling is then carried out which uses the same parameters as the previous one, but this time by conducting the calculation XFEM (which is then a pure mechanical calculation) with the operator MECA_STATIQUE instead of the operator STAT_NON_LINE in the case HM-XFEM. The results obtained are strictly identical to the previous ones.