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:

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]

Disproportionate deformation leads to buckling or failure. In machines with moving parts this can reduce the size of gaps and lead to frictional abrasion.

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.

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.

Calculating the creep strain

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

The elastic strain \(\mathbf{\epsilon}_{el}\) occurs right after applying the load and is connected with the stresses through the elasticity tensor \(\mathbf{D}\) in equation 1. 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 a 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}$$

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.

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

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.

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 given related 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}\)

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.

This website uses cookies so that we can provide you with the best user experience possible. Cookie information is stored in your browser and performs functions such as recognising you when you return to our website and helping our team to understand which sections of the website you find most interesting and useful.

Strictly Necessary Cookies

Strictly Necessary Cookie should be enabled at all times so that we can save your preferences for cookie settings.

If you disable this cookie, we will not be able to save your preferences. This means that every time you visit this website you will need to enable or disable cookies again.