2. Modeling A#

2.1. But#

The aim of this test is to validate the sending and the transformation (for law VMIS_JOHN_COOK) of a field of internal variables into the initial state of a EPX calculation via CALC_EUROPLEXUS. It also validates the transformation of internal variables in the EPX direction to Code_Aster for the VMIS_JOHN_COOK law and the use of LIRE_EUROPLEXUS outside of CALC_EUROPLEXUS.

2.2. Description#

This test is the equivalent of « bench » EPX bm_ini_med_ecro_vmjc.

It is a cube formed by a single element. The lower face is embedded and a loading (displacement imposed along X) is applied to the upper face. This calculation is done with the operator STAT_NON_LINE in order to produce an initial state for the calculation EPX.

We then launch CALC_EUROPLEXUS with this initial state (displacements + constraints + internal variables) and without additional loading than those that allowed us to obtain the initial state.

We then check that the initial internal variables and in final state are indeed the same as in bm_ini_med_ecro_vmjc using the keyword COURBE (and TEST_TABLE) and that the internal variables retrieved by Code_Aster at the end of the calculation are indeed those expected.

The internal variables should normally not change since no additional load has been added. However, some components that do not depend on the calculated values are completely recalculated by EPX at each iteration. This is the case with ECR5.

2.3. Principle of validation#

Comparison with test case EPX or analytical references.

2.4. Tested values#

In this first table the numbers of the internal variables are those of EPX.

Mesh

Instant

Component

Component

Reference Value

Tolerance

M1

\(\mathit{V1}=\mathit{ECR1}\)

1

4.078819E+2

1E-4

M1

\(\mathit{V2}=\mathit{ECR2}\)

2

1.420337E+2

1E-4

M1

\(\mathit{V3}=\mathit{ECR3}\)

3

5.773488E-3

1E-4

M1

\(\mathit{V4}=\mathit{ECR4}\)

4

1.561430E+2

1E-4

M1

0.004

\(\mathit{V1}=\mathit{ECR1}\)

1

1

4.078819E+2

1E-4

M1

0.004

\(\mathit{V2}=\mathit{ECR2}\)

2

1.420337E+2

1E-4

M1

0.004

\(\mathit{V3}=\mathit{ECR3}\)

3

3

5.773488E-3

1E-4

M1

0.004

\(\mathit{V4}=\mathit{ECR4}\)

4

4

1.561430E+2

1E-4

M1

0.004

\(\mathit{V5}=\mathit{ECR5}\)

5

5

4019.1847623425019

1E-4

We then test the value of the internal variables in the Aster result at the output of CALC_EUROPLEXUS.

Mesh

Instant

Component

Component

Reference Value

Tolerance

M1

\(\mathit{V1}\)

3

5.773488E-3

1E-4

M1

0.004

\(\mathit{V1}\)

3

3

5.773488E-3

1E-4

Finally, we test the value of the internal variables in the Aster result at the output of LIRE_EUROPLEXUS.

Mesh

Instant

Component

Component

Reference Value

Tolerance

M1

\(\mathit{V1}\)

3

5.773488E-3

1E-4

M1

0.004

\(\mathit{V1}\)

3

3

5.773488E-3

1E-4