4. Modeling A#

4.1. Characteristics of modeling#

3D modeling

Item MECA_TETRA4.

4.2. Characteristics of the mesh#

Number of knots: 4

Number of meshes and types: 1 TETRA4

4.3. Charging path#

Charging is divided into two phases:

      • Phase 1: Traction in imposed displacement

\(\mathit{DX}\mathrm{=}{F}^{\mathit{trac}}\) imposed on \(\mathit{N1}\)

      • Phase 2: Traction/Shear in imposed displacement

\(\mathit{DX}\mathrm{=}\mathit{DY}\mathrm{=}0.5\mathrm{\ast }{F}^{\mathit{cisa}}\) imposed on \(\mathit{N1}\)

\(\mathit{DY}\mathrm{=}0.75\mathrm{\ast }{F}^{\mathit{cisa}}\) imposed on \(\mathit{N2}\)

where \({F}^{\mathit{trac}}\) and \({F}^{\mathit{cisa}}\) are time-increasing affine functions

4.4. Tested sizes and results#

The non-regression test is carried out on the value of the rotation angles of the deformation reference frames and the damage tensor.

To do this, we extract the deformation (EPSI_ELGA) and damage (VARI_ELGA) fields at time 2, and we create the matrices (in Python) corresponding to the deformation and damage tensors. Next, we use Python’s LinearAlgebra library to calculate the eigenvectors of matrices associated with deformation and damage tensors. Finally, the angle of rotation of these eigenvectors with respect to the initial coordinate system is calculated.

Instant

Field Name

Component

Location

Aster

2

EPSI_ELGA

Angle of rotation of the proper coordinate system

\(\mathit{VOLUME}\), point 1

45.7160

2

VARI_ELGA

Angle of rotation of the proper coordinate system

\(\mathit{VOLUME}\), point 1

23.0825