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:
ϵ˙cr=f(σ,ϵ,T,t)(1)
$$\begin{array}{}\text{(1)}& {\dot{\u03f5}}_{cr}=f(\sigma ,\u03f5,T,t)\end{array}$$
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]
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.
Creep can be divided in three different stages: primary creep, secondary creep and tertiary creep
Primary creep (0<m<1
$0<m<1$) starts rapidly with an infinite creep rate at the initialization. Here is m
$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
$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.
The creep strain equations are based on an additive strain decomposition:
ϵtot=ϵel+ϵpl+ϵcr(2)
$$\begin{array}{}\text{(2)}& {\u03f5}_{tot}={\u03f5}_{el}+{\u03f5}_{pl}+{\u03f5}_{cr}\end{array}$$
σ=D:ϵel(3)
$$\begin{array}{}\text{(3)}& \sigma =\mathbf{D}:{\u03f5}_{el}\end{array}$$
ϵ˙cr=32ϵc˙σ~σeq(4)
$$\begin{array}{}\text{(4)}& {\dot{\u03f5}}_{cr}=\frac{3}{2}\dot{{\u03f5}^{c}}\frac{\stackrel{~}{\sigma}}{{\sigma}_{eq}}\end{array}$$
where we have
: total strain tensor
: elastic strain tensor
: plastic strain tensor
: creep strain tensor
: equivalent creep strain
: deviatoric stress tensor
: equivalent stress
: elasticity tensor
The elastic strain ϵel
${\u03f5}_{el}$occurs right after apllying the load and is connected with the stresses through the elasticity tensor D
$\mathbf{D}$in equation (3)
$\text{(3)}$. The plastic strain ϵpl
${\u03f5}_{pl}$results from excessive loading and is irreversible. The creep strain tensor can be calculated with the the equivalent creep strain ϵc
${\u03f5}^{\mathbf{c}}$, the deviatoric stress σ~
$\stackrel{~}{\sigma}$and the equivalent stress σeq
${\sigma}_{eq}$. The deviatoric stress σ~
$\stackrel{~}{\sigma}$causes distortion of the body and results after splitting the stress matrix in an isotropic and an deviatoric component.
Two fundamental creep law types defined as creep strain rate equations are available in SimScale:
ϵc˙=A⋅σn⋅tm⋅ϵk
$$\dot{{\u03f5}^{c}}=A\cdot {\sigma}^{n}\cdot {t}^{m}\cdot {\u03f5}^{k}$$
ϵc˙=ϵc0⋅sinh(σσ0)n
$$\dot{{\u03f5}^{c}}={\u03f5}_{0}^{c}\cdot \mathrm{sinh}(\frac{\sigma}{{\sigma}_{0}}{)}^{n}$$
where we have
: a constant, depends on problem
: creep stress index, 1<=n<=10
$1<=n<=10$
: time index, 0<m<=1
$0<m<=1$
: 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.
[NAFEMS_HT21]In this formulation the creep strain rate only depends on the stresses. Therefore m
$m$and k
$k$equal zero.
ϵc˙=A⋅σn
$$\dot{{\u03f5}^{c}}=A\cdot {\sigma}^{n}$$
In this formulation the creep strain rate depends also on time.
ϵc˙=A⋅σn⋅tm
$$\dot{{\u03f5}^{c}}=A\cdot {\sigma}^{n}\cdot {t}^{m}$$
In this formulation the creep strain rate depends on the stresses and the creep strains. (For convenience k
$k$is referred to as m
$m$on the platform.)
ϵc˙=A⋅σn⋅(ϵc)k
$$\dot{{\u03f5}^{c}}=A\cdot {\sigma}^{n}\cdot ({\u03f5}^{c}{)}^{k}$$
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.
ϵc˙=ϵc0⋅sinh(σσ0)n
$$\dot{{\u03f5}^{c}}={\u03f5}_{0}^{c}\cdot \mathrm{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.
The creep formulation parameters (A,m,n)
$(A,m,n)$and (ϵc0,σ0)
$({\u03f5}_{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 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: ASI=13600⋅13600m⋅1106n⋅AMPah
${A}_{SI}=\frac{1}{3600}\cdot \frac{1}{{3600}^{m}}\cdot \frac{1}{{10}^{6n}}\cdot {A}_{MPa\phantom{\rule{thickmathspace}{0ex}}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.
[NAFEMS_HT21] | (1, 2) (2001) “How to Undertake Creep Analysis with FInite Elements”, NAFEMS Publication Ref: HT21 |