Documentation

3D Triaxial Load Secondary Creep (NAFEMS Test 6(a))

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 analytical 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 hex 20 Strain Hardening 3D standard
(F) 2nd order hex 20 Strain Hardening 3D reduced
(G) 2nd order tet 21 Norton 3D standard
(H) 2nd order tet 21 Norton 3D reduced
(I) 2nd order tet 21 Time Hardening 3D standard
(J) 2nd order tet 21 Time Hardening 3D reduced
Mesh used for the SimScale cases (G) – (J)

Simulation Setup

Material:

  • linear elastic, isotropic: E = 200 GPa, ν

    ν

     = 0.3

  • creep formulation:

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

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

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

Important

The material Data of the reference was given in MPamm and h units and were converted to SI units using the following relation:
ASI=136001106nAMPah

ASI=136001106nAMPah

 

Constraints:

  • face ADEH: dx=0.0
    dx=0.0

     

    m

  • face ABFE: dy=0.0
    dy=0.0

     

    m

  • face ABCD: dz=0.0
    dz=0.0

     

    m

Loads:

  • face BCGF: σxx
    σxx

     

    = 300 MPa

  • face CDHG: σyy
    σyy

     

    = 200 MPa

  • face EFGH: σzz
    σzz

     

    = 100 MPa

Advanced Automatic Time Stepping:

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

Reference Solution

ϵcxx=ϵczz=0.0042183600t(1)

(1)ϵxxc=ϵzzc=0.0042183600t

 

ϵceff=0.0048713600t(2)

(2)ϵeffc=0.0048713600t

 

ϵcyy=0.0(3)

(3)ϵyyc=0.0

 

The equations used to solve the problem are derived in [NAFEMS_R27]. As the SimScale solution is calculated using SI units, the reference solution was adopted to a time unit of seconds instead of hours.

Results

Comparison of the average creep strain ϵcxx

ϵxxc

 of the cube after a creep time of 3.6e6 s (1000 h) with the analytical solution presented in [NAFEMS_R27]:

Comparison of the creep strain ϵcxx

ϵxxc

 

Case time [s] [NAFEMS_R27] SimScale Error (%)
(A) 3600000 4.218 4.21875 -0.0178
(B) 3600000 4.218 4.21875 -0.0178
(C) 3600000 4.218 4.21875 -0.0178
(D) 3600000 4.218 4.21875 -0.0178
(E) 3600000 4.218 4.21875 -0.0178
(F) 3600000 4.218 4.21875 -0.0178
(G) 3600000 4.218 4.21875 -0.0178
(H) 3600000 4.218 4.21875 -0.0178
(I) 3600000 4.218 4.21875 -0.0178
(J) 3600000 4.218 4.21875 -0.0178
Contents