# 3D triaxial displacement secondary creep (NAFEMS Test 6(b))¶

## Overview¶

The aim of this test case is to validate the following functionalities:

• creep material behavior
• standard and reduced integration elements
• automatic time stepping

The simulation results of SimScale were compared to the results in [NAFEMS_R27]. The tetrahedral meshes used were created with the parametrized-tetrahedralization-tool on the SimScale platform.

Import validation project into workspace

## Geometry¶

Geometry of the Cube

• Cube edge length l = 0.100 m

## Analysis type and Domain¶

Tool Type : Code_Aster

Analysis Type : Static - nonlinear

Mesh and Element types :

Case Mesh type Number of nodes Creep formulation Element type
(A) 2nd order hex 20 Norton 3D standard
(B) 2nd order hex 20 Norton 3D reduced
(C) 2nd order hex 20 Time Hardening 3D standard
(D) 2nd order hex 20 Time Hardening 3D reduced
(E) 2nd order tet 21 Norton 3D standard
(F) 2nd order tet 21 Norton 3D reduced
(G) 2nd order tet 21 Time Hardening 3D standard
(H) 2nd order tet 21 Time Hardening 3D reduced

Mesh used for the SimScale cases (E) - (H)

## Simulation Setup¶

Material:

• linear elastic, isotropic: E = 200 GPa, $$\nu$$ = 0.3

• creep formulation:

– Norton, A = 8.6805556e-48, n = 5

– Time Hardening, A = 8.6805556e-48, n = 5, m = 0

Important

The material Data of the reference was given in MPa, mm and h units and were converted to SI units using the following relation:
$$A_{SI} = \frac{1}{3600} \cdot \frac{1}{10^{6n}} \cdot A_{MPa\;h}$$

Constraints:

• face ADEH: $$d_{x} = 0.0$$ m
• face ABFE: $$d_{y} = 0.0$$ m
• face ABCD: $$d_{z} = 0.0$$ m
• face BCGF: $$d_{x} = 0.3$$ mm
• face CDHG: $$d_{y} = 0.2$$ mm
• face EFGH: $$d_{z} = 0.1$$ mm

Advanced Automatic Time Stepping:

• Retiming Event: Field Change
• Target Field component: internal variable V1 (accumulated unelastic strain)
• Threshold value: 0.00005
• Timestep Calculation Type: mixed
• Field Change Target Value: 0.00004

## Results¶

Comparison of the average Cauchy stress $$\sigma_{xx}$$ of the cube after a creep time of 3.6e6 s (1000 h) with the solution presented in [NAFEMS_R27]:

Comparison of the Cauchy stress $$\sigma_{xx}$$[ MPa ]
Case time [s] [NAFEMS_R27] SimScale Error (%)
(A) 3600000 1007.7 1008.23 -0.053
(B) 3600000 1007.7 1008.23 -0.053
(C) 3600000 1007.7 1008.23 -0.053
(D) 3600000 1007.7 1008.23 -0.053
(E) 3600000 1007.7 1008.23 -0.053
(F) 3600000 1007.7 1008.23 -0.053
(G) 3600000 1007.7 1008.23 -0.053
(H) 3600000 1007.7 1008.23 -0.053