Creep is the inelastic, irreversible deformation of structures during time. It is a life limiting factor and depends on stress, strain, temperature and time. This dependency can be modeled as followed:

\[\begin{equation}\label{creep_dependency_eq} \dot{\mathbf{\epsilon}}_{cr} = f(\sigma,\epsilon,T,t) \end{equation}\]

Creep can occur in all crystalline materials, such as metal or glass, has various impacts on the behavior of the material and can ultimately lead to following problems: [NAFEMS_HT21]

  1. Disproportionate deformation leads to buckling or failure. In machines with moving parts this can reduce the size of gaps and lead to frictional abrasion.
  2. The connection of bolted parts can loosen over time, because of relaxations of the pre-stressing. This relaxation can also occur in cables, gaskets and flanges.
  3. Excessive heat and load cycles can favor the propagation of cracks, which can lead to damage and ultimately to failure of the part.

Because of these extremely different factors the effects of creep in certain situations can be highly complex. By undertaking creep analysis the influence of these effects can be evaluated and the life time of parts can be estimated. Especially in the hot section of turbines of commercial and military plains this can be a crucial factor for the design process and save a lot of money.

Three stages of creep

Creep can be divided in three different stages: primary creep, secondary creep and tertiary creep

Primary creep (\(0 < m < 1\)) starts rapidly with an infinite creep rate at the initialization. Here is \(m\) the time index. It occurs after a certain amount of time and slows down constantly. It occurs in the first hour after applying the load and is essential in calculating the relaxation over time.

Secondary creep (\(m = 1\)) follows right after the primary creep stage. The strain rate is now constant over a long period of time.

The strain rate in the tertiary creep stage is growing rapidly until failure. This happens in a short period of time and is not of great interest. Therefore only primary and secondary creep are modeled on the SimScale platform.


Creep strain and the three main creep stages

Calculating the creep strain

The creep strain equations are based on an additive strain decomposition:

\[\begin{equation}\label{strain_decomposition} \mathbf{\epsilon}_{tot} = \mathbf{\epsilon}_{el} + \mathbf{\epsilon}_{pl} + \mathbf{\epsilon}_{cr} \end{equation}\]\[\begin{equation}\label{small_strain_hypo} \mathbf{\sigma} = \mathbf{D} : \mathbf{\epsilon}_{el} \end{equation}\]\[\begin{equation}\label{creep_strain_eq} \dot{\mathbf{\epsilon}}_{cr} = \frac{3}{2} \dot{\epsilon^c} \frac {\tilde{\mathbf{\sigma}}}{\sigma_{eq}} \end{equation}\]

where we have

  • \(\mathbf{\epsilon}_{tot}\) : total strain tensor
  • \(\mathbf{\epsilon}_{el}\) : elastic strain tensor
  • \(\mathbf{\epsilon}_{pl}\) : plastic strain tensor
  • \(\mathbf{\epsilon}_{cr}\) : creep strain tensor
  • \(\mathbf{\epsilon^c}\) : equivalent creep strain
  • \(\tilde{\mathbf{\sigma}}\) : deviatoric stress tensor
  • \(\sigma_{eq}\) : equivalent stress
  • \(\mathbf{D}\) : elasticity tensor

The elastic strain \(\mathbf{\epsilon}_{el}\) occurs right after apllying the load and is connected with the stresses through the elasticity tensor \(\mathbf{D}\) in equation \(\eqref{small_strain_hypo}\). The plastic strain \(\mathbf{\epsilon}_{pl}\) results from excessive loading and is irreversible. The creep strain tensor can be calculated with the the equivalent creep strain \(\mathbf{\epsilon^c}\), the deviatoric stress \(\tilde{\mathbf{\sigma}}\) and the equivalent stress \(\sigma_{eq}\). The deviatoric stress \(\tilde{\mathbf{\sigma}}\) causes distortion of the body and results after splitting the stress matrix in an isotropic and an deviatoric component.

Calculating the equivalent creep strain

Two fundamental creep law types defined as creep strain rate equations are available in SimScale:

  • Power or Bailey-Norton Law:
\[\dot{\epsilon^c} = A \cdot \sigma^{n} \cdot t^m \cdot \epsilon^{k}\]
  • Hyperbolic Sine Law:
\[\dot{\epsilon^c} = \epsilon_{0}^c \cdot \operatorname{sinh} ( \frac{\sigma}{\sigma_0} )^{n}\]

where we have

  • \(A\) : a constant, depends on problem
  • \(n\) : creep stress index, \(1 <= n <= 10\)
  • \(m\) : time index, \(0 < m <= 1\)
  • \(k\) : creep strain index

The creep stress index depends on the temperature and the stress level and can be determined with help of an Ashby deformation mechanism map.


Example of an Ashby deformation mechanism map [NAFEMS_HT21]

On basis of these laws four different creep law formulations are currently available on the platform:

  • Norton (Power Law):
In this formulation the creep strain rate only depends on the stresses. Therefore \(m\) and \(k\) equal zero.
\[\dot{\epsilon^c} = A \cdot \sigma^{n}\]
  • Time Hardening (Power Law):
In this formulation the creep strain rate depends also on time.
\[\dot{\epsilon^c} = A \cdot \sigma^{n} \cdot t^m\]
  • Strain Hardening (Power Law):
In this formulation the creep strain rate depends on the stresses and the creep strains. (For convenience \(k\) is referred to as \(m\) on the platform.)
\[\dot{\epsilon^c} = A \cdot \sigma^{n} \cdot (\epsilon^c)^k\]
  • Garofalo (Hyperbolic Sine Law)
The hyperbolic sine law accounts the stresses \(\sigma\) differently. For small values of \(\sigma\) the stress dependency behaves similar and when \(\sigma\) gets zero the stress dependency becomes zero as well.
\[\dot{\epsilon^c} = \epsilon_{0}^c \cdot \operatorname{sinh} ( \frac{\sigma}{\sigma_0} )^{n}\]

The choice of the formulation to be used varies from problem to problem and depends on the data of the temperature and stress field of the problem.

Important remarks

  1. The creep formulation parameters \((A, m, n)\) and \((\epsilon_0^c, \sigma_0)\) are defined in the sense of the above equations with respect to the rate form of the creep laws.

  2. It is important to have a consistent unit system when converting parameters found in literature or material supplier or test data. Often the creep parameters are are given relate to length measure in mm, stresses in MPa and time in hours.

    An example of how to convert the parameter A of the Time Hardening or Strain Hardening formulation would be: \(A_{SI} = \frac{1}{3600} \cdot \frac{1}{3600^m} \cdot \frac{1}{10^{6n}} \cdot A_{MPa\;h}\)

  3. Especially when primary creep is observed it is very important to refine the time stepping sufficiently to capture the creep behavior correctly. Therefore it is advised to use the Field Change criteria of the Automatic Time Stepping available for Code_Aster.