Creep Phenomena, Mechanisms and Modeling

1. Introduction

Creep refers to time-dependent deformation under sustained loading, which is an important subject in engineering design for high-temperature applications. The phenomenon was first observed by Andrade (1910) [1,2], more than 110 years ago. To date, it remains to be and becomes even more critical to evaluate high-temperature components especially those operating in gas turbines, steam turbines, and pressure vessels, etc., as they tend to be exposed to higher and higher temperatures to meet the CO₂ emission reduction target.

Depending on the application, engineering creep failure criteria could be given as [3]:

I. creep strain to a specified level, e. g. 1%.

II. creep life for intended hours of service, e. g. 100,000 hours.

Type I criteria apply to cases where change of shape/dimension exceeds the tolerance limit to cause unwanted interference with adjacent components, e.g., elongation of rotating turbine blades rubbing the casing. Type II criteria generally apply to static high-temperature components for which economic safe-operation is of paramount concern. Then, in design data sheets, material creep strength can be defined either as the stress to cause a certain amount of creep strain for a prescribed time, e.g., 100 h, 1,000h, and 10,000h, etc.; or as the stress to cause creep rupture for a prescribed life, e.g., 100,000h, at given temperatures. These data are often obtained through creep coupon testing. Sometimes, creep strength is also defined in terms of the temperature capability for a given performance (time at stress). For example, the creep strength of a new generation single crystal Ni-base superalloys is marked by an increase of 25-30^oC for 1,000h to reach 1% strain at 137 MPa [4].

In addition, modern machineries generally operate under complex loading conditions consisting of both mechanical and thermal load cycles, which would result in creep-fatigue interaction, or thermomechanical fatigue. The description of such complicated problems relies on a thorough understanding of deformation and damage mechanisms in both fatigue and creep. However, due to the limited space, the topic of creep-fatigue interaction will be left out of the scope.

In general, creep design often concerns the following questions:

1. What is the creep deformation behavior, i.e., creep strain/rate as a function of time, stress and temperature, to the point of rupture?

2. How does the creep life vary with stress and temperature?

3. What are the effects of chemical composition and microstructure (via manufacturing process and heat treatment, etc.) on the material’s creep resistance?

4. How does the material microstructure change during long-time service at high temperatures, and if changed, what is its effect on the material’s creep performance?

5. Does the operating environment have an effect on the material’s creep performance?

Numerous studies have been dedicated to address the above questions in the past, but mostly on one or two aspects in each case study (article). Some fundamentals of creep have been covered in a recent book [5], yet unresolved issues still need further discussion. This article summarizes the observed creep phenomena and reviews the state-of-the-art of creep modeling. Particularly, the emphasis is given to the mechanism-based modeling approach describing the observed phenomena of creep deformation, fracture mode and life with constitutive equations derived from the underlying mechanisms for complex engineering alloys over a wide range of stress and temperature conditions. This way, it is hoped that design engineers can have a full physical picture of creep in order to achieve optimal design for various industrial applications.

2. Creep Phenomena

In this section, creep phenomena are depicted by empirical equations, as when they were first observed. Interpretations of these phenomena will be made using mechanism-based creep models with experimental data discussed in later sections.

2.1. Creep Deformation

Creep deformation occurs in the homologous temperature range of 0.3 – 0.9 T/T_m, where T is the application temperature and T_m is the material’s melting temperature in Kelvin degree. It generally proceeds in three stages: (1) primary; (2) secondary, or the steady-state; and (3) tertiary stage, as schematically shown in Figure 1. In the primary stage, creep rate is decelerating; in the secondary stage it exhibits a nearly constant creep rate, which is often known as the minimum creep rate; and in the tertiary stage creep rate is accelerating to the point of rupture at last. Creep deformation can occur under stresses well below the material’s yield strength. Therefore, it is distinguished from rate-independent plasticity, even though both lead to permeant shape change. Basically, all creep phenomena such as creep strain, creep rate, creep life and creep ductility are attributes of the creep curve (process). Metallurgical examinations using scanning electron microscopy (SEM) and transmission electron microscopy (TEM) on the microstructure and fracture surfaces of crept specimens are often performed to obtain information for understanding the mechanisms of creep deformation and fracture of the materials.

Many empirical equations have been proposed to describe creep curves. Norton (1929) and Bailey (1935) proposed a generalized power-law [6,7]:

$$\epsilon =A{\sigma }^n{t}^m$$

(1)

where t is time, A, n, and m are temperature dependent material constants (in Andrade’s equation, m = 1/3).

Graham-Walles (1955) expanded the Norton-Bailey equation in a series form to describe the three-stage creep deformation (m = 1/3 for the primary stage, m = 1 for the secondary stage, and m > 1 for the tertiary stage) [8].

Evans and Wilshire (1985) proposed a θ-projection method which consists of two exponential terms as [9]

$$\epsilon ={\theta }_1\left(1-e^{{\theta }_{2}t}\right)+{\theta }_3\left(e^{{\theta }_{4}t}-1\right)$$

(2)

where the θ_i parameters (i = 1, 2, 3, 4) are polynomial functions of stress (σ) and temperature (T), fitting to the experimental data.

These empirical equations were proposed only to depict creep deformation, but bore no meaning on the underlying mechanisms, the rupture (intergranular/transgranular) mode and time to failure of the creep process. So, they are not able to predict creep rupture life, unless Type I failure criteria are prescribed.

2.2. Creep Life

Creep rupture life is an important design consideration for high temperature components such as gas turbine components, pressure vessels and nuclear reactors, etc. Larson and Miller (1952) proposed a data collation method with the Larson-Miller parameter (LMP) defined as [10]

$$LMP=T\left(C_{LM}+\mathrm{l}\mathrm{o}\mathrm{g}{t}_r\right)=f\left(\sigma \right)$$

(3)

where C_LM is called Larson-Miller constant which often takes a value of ~ 20, f(σ) is a function of applied stress.

Using LMP, the creep strength can be determined from σ-LMP correlation (usually, a polynomial fit). The LMP offers a method to predict creep life assuming that the σ-LMP correlation is a master curve for the material. It implies that short-term data obtained at higher temperatures can be used to predict long-term creep life at lower temperatures, if the LMP of the two conditions is the same. However, due to the scatter in experimental lives, the uniqueness of Eq. (3) has never been proven. In addition, the LM method has inherent ambiguities relating to physical mechanisms. Thus, the reliability of LMP method is still in question. Nevertheless, this method is popularly adopted in creep design, perhaps owing to its straightforwardness, as long as enough data are generated to cover the LMP range of application.

Another popular parametric representation of creep strength is the Orr-Sherby-Dorn (OSD) parameter method defined as [11]

$$P_{\mathrm{OSD}}=\log t_{\mathrm{r}}-\frac{Q_{\mathrm{c}}}{2.3RT}=f\left(\sigma \right)$$

(4)

where Q_c is the activation energy for creep (J·mol⁻¹), R = 8.314 J·K^-1mol^-1 is the universal gas constant. P_OSD is usually a negative number. Similar to the LMP method, the OSD approach also suffers from its empirical nature.

Monkman and Grant (1956) proposed a direct correlation of creep life to the steady-state creep rate, as [12]

$${\dot{\epsilon }}_{min}^m{t}_r=C$$

(5)

where ${\dot{\epsilon }}_{min}$ is the minimum creep rate, as

$${\dot{\epsilon }}_m=A_0\mathrm{e}\mathrm{x}\mathrm{p}\left(\frac{Q_c}{RT}\right){\left(\frac{\sigma }{{\sigma }_0}\right)}^n$$

(6)

where A₀ is a proportional constant, Q_c is the activation energy, n is the stress exponent, σ₀ is the normalization stress.

A more sophisticated creep life equation was proposed by Wilshire [13], as

$$\sigma /{\sigma }_{TS}=\mathrm{e}\mathrm{x}\mathrm{p}\left\{-k_1{\left[t_r\mathrm{e}\mathrm{x}\mathrm{p}\left(\frac{Q_c}{RT}\right)\right]}^u\right\}$$

(7)

and

$$\sigma /{\sigma }_{TS}=\mathrm{e}\mathrm{x}\mathrm{p}\left\{-k_2{\left[{\dot{\epsilon }}_m\mathrm{e}\mathrm{x}\mathrm{p}\left(\frac{Q_c}{RT}\right)\right]}^v\right\}$$

(8)

where σ_TS is the material’s tensile strength, k₁, k₂, u, and v are material constants, Q_c is the activation energy determined from the temperature dependence of έ_m and/or t_r at constant (σ/σ_TS). Eq. (7) describes the creep rupture life, and Eq. (8) describes the minimum creep rate. It infers that t_r → 0 as σ/σ_TS → 1, whereas έ_m → 0 and t_r → ∞ when σ/σ_TS → 0. These equations have been applied to many materials including 9Cr steel [13], 1Cr-1Mo-0.25V steel [14], and Waspaloy [15]. In this method, a linear regression of log (t_r·exp(Q_c/RT)) vs. log(−log(σ/σ_TS)) data is employed to obtain the slope u, and intercept logk₁.

Abdallah et al. performed a critical analysis of the above and a few more empirical creep lifing equations against Grade 22 steel creep data from National Institute of Materials Science (NIMS), Japan [16]. They pointed out that in the Wilshire approach when plotting to obtain k₁ and u, multiple linear regimes may exist, which is thought to due to changes in creep mechanism. To ascertain the stress threshold value of regime-change and determine k₁ and u values in different regimes, more creep rupture tests have to be conducted.

For single-mechanism-controlled processes, which occur in a narrow stress range, both LMP/OSD and Wilshire approaches are practically equivalent. But for multi-mechanism creep over a wide range of stress and temperature, the stress-dependence of creep life turns into either a polynomial function (LM) or multi-regime exponential function (Wilshire), which require a lot of creep data to establish the correlation and therefore diminish the predictive capability of the model. Most of the time, the above relationships are used to correlate short-term creep test data, and extrapolation of such correlation for long-term creep life prediction is desirable, as a design life of 100,000 h can hardly be experimentally-validated before the product is put into the market for commercial competition. The reliability of extrapolation of those short-term empirical relationships is often in question, because many factors such as microstructural evolution and environment cannot be fully taken into account in short-term creep tests, which can significantly affect the material’s long-term creep performance. Furthermore, if short-term creep tests are conducted in different mechanism regimes from the long-term service conditions, extrapolations from such tests are not physically warranted. Hence, prediction of long-term creep properties from short-term experiments is still rated as one of the most important challenges in assessing structural durability, as regarded by the energy sector of the United Kingdom [17]. To fully resolve the issue(s) requires complete understanding of the potential service failure as well as the deformation and failure mechanism(s) in the design creep tests.

2.3. Creep Ductility and Fracture Mode

Creep ductility refers to the strain level at creep rupture, i.e., the last point on the creep curve, as shown in Figure 1. It may also be defined as the cross-sectional area reduction at fracture [18]. To keep in consistency with deformation modeling, the strain definition is used hereafter. One of the notorious effects of creep is to cause an otherwise ductile material become very “brittle” (with very low ductility) after creep for a period of time under certain stress/temperature conditions. Spindler proposed a creep ductility equation for steels, as [19]

$${\epsilon }_{cu}=min\left\{{\epsilon }_{cu,U},max\left[{\epsilon }_{cu,L},C_1{e}^{Q^{*}/T}{\left({\dot{\epsilon }}_{min}\right)}^{n_1}{\sigma }^{m_1}\right]\right\}$$

(9)

where ε_cu,L/U are the lower/upper shelf creep rupture strains, C₁, n₁, m₁ are material constants, Q* is related to diffusion activation energy. Holdsworth used Eq. (9) to draw a schematic trending curve [20], as shown in Figure 2, and divided the curve into three regimes: Regime I was plastic hole growth leading to ductile dimpled failure; Regime II was diffusion-controlled cavity growth and creep ductility decreases with reducing strain rate, whereas in regime III cavity growth was constrained and creep ductility is independent of creep strain rate.

It has been commonly observed that low creep ductility is often associated with intergranular fracture and high creep ductility is associated with transgranular fracture. Eq. (9) describes the variation of creep ductility with stress and temperature for steels. All the creep equations in this section were not derived from deformation mechanisms, but rather mathematical correlations of experimental observations, so they are empirical in nature. As one can see, phenomenological equations from Eq. (1) to Eq. (9) do not relate to each other. If one chose Eq. (1-2) to depict creep strain, he/she would not know when the life would end. If one chose Eq. (3-7) to estimate the creep life, he/she would have no idea how creep strain accumulated to the point of fracture. And, Eq. (9) certainly would not help in understanding the strain accumulation and life. To overcome the conundrum, one ought to resort to understanding the creep deformation and damage mechanisms and derive physics-based equations to link the above phenomena as discussed below.

3. Creep Mechanisms

Creep deformation mechanisms have been studied extensively using microstructural characterization means including SEM, TEM, etc. Generally, it is found that intragranular dislocation movements such as intragranular dislocation glide (IDG) and intragranular dislocation climb (IDC), diffusion, and grain boundary sliding (GBS) could all operate, leading to time-dependent creep deformation. The dominance of these deformation mechanisms depends on the stress and temperature of application. This is often reflected in the apparent power-law exponent and activation energy.

Based on the thermodynamics and kinetics of slip, an IDG model was proposed by Kocks, Argon and Ashby [21]. Mechanism-based models were also developed to show that IDC proceeds by the power law, as Eq. (6), with an exponent of 3 for dilute solution alloys [22,23] and 4.5 for metals with formation of sub-grains [24]. Creep may also occur as a result of matter flow by either lattice diffusion (Nabarro-Herring creep) or grain boundary diffusion (Coble creep) with a stress exponent of 1 [23]. In addition, GBS may proceed by dislocation glide plus climb along grain boundaries and it is shown to have a stress-exponent of 2 for materials with clean and planar grain boundaries (no grain boundary precipitates) [25]. In complex engineering alloys containing solute atoms, precipitates in the matrix as well as at grain boundaries interacting with moving dislocations, IDC and GBS tend to have higher stress exponents [26,27,28]. GBS usually replaces diffusion creep, as grain boundary precipitate denuded zones are rarely observed, which means no grain boundary migration by diffusion due to particle-strengthening in complex engineering alloys [29,30]. IDG in complex engineering alloys may also be affected by the presence of strengthening particles through Orowan-looping or precipitate-cutting [31].

Frost and Ashby (1982) constructed deformation-mechanism maps for many metals [23], which graphically shows the regions of dominance of particular mechanisms in the stress–temperature field with the boundaries being defined by the line of equal contributions from both sides. An example of deformation mechanism map for pure nickel is shown in Figure 3. Usually, IDC (power-law) creep is dominant in the high-temperature mid-stress region. As the stress increases, exponential stress-dependence of creep rate is manifested, which is called power-law breakdown, meaning that deformation enters into the IDG region. Full discussions of deformation mechanism maps can be referred to Frost and Ashby (1982) for pure metals [21]and Wu (2019) for some engineering alloys [5].

It has long been noticed that metals with good tensile ductility could fail with a relatively low elongation after creep for a prolonged period in the homologous temperature range of 0.3 to 0.9 [32]. The reduction in creep ductility is thought to be caused by several grain boundary weakening mechanisms: (1) cavity nucleation and growth by vacancy condensation under stress, either by free diffusion of vacancies [33] or constrained diffusion accommodated by power-law creep [34,35]; (2) grain boundary sliding (GBS) leading to cavity growth from ledges, grain boundary precipitates and/or triple junction (wedge) cracks [36,37,38,39]; and (3) cavity nucleation from Zener–Stroh dislocation pileups against high-angle grain boundaries [40].

4. Creep Modeling

Since the mid-twentieth century, a tremendous amount of efforts have been made on creep modeling, because of the importance of the problem in engineering design. Usually, empirical relations were developed first to describe one aspect of creep phenomena. When deviation occurred, as more materials were being examined, either the existing model was amended or new empirical models were proposed to describe the “new” creep data. This is how the section-2 equations were developed, and those are only a small fraction of what have been proposed in the literature. In the end, one has to manage a set of unrelated empirical relations with many parameters. This requires a significant effort of testing and analysis for a complex engineering alloy.

This review paper tries to link various creep attributes through mechanism-based modeling, in particular, deformation-mechanism-based true-stress model (DMTS) [5]. This model simulates creep strain accumulation by IDG, IDC and GBS under true stress conditions with a failure criterion that either critical intragranular deformation (IDG+IDC) or GBS levels are reached first to cause fracture. Hence, creep curve, ductility, life and fracture mode are evaluated from the partition of these physical mechanisms. By the kinetics, the dominance of each mechanism manifests naturally at each (primary, secondary, tertiary) stage of creep before unstable deformation (necking) and fracture occurs. Using the mechanism-based approach, the effects of microstructure and environment can also be identified in association with particular deformation mechanisms. Therefore, mechanism-informed insights can be gained through the DMTS model analysis. The model has been applied to many engineering alloys including Waspaloy [41], Mar-M 509 [42], modified Grade 91 steel [43], G115 [44], SA508Gr3 steel [45], Alloy 800 [46], and Haynes^® 282^TM [47].

4.1. Deformation-Mechanism-Based True-Stress Model

In general, IDG/IDC and GBS proceed independently in a deformation process, leading to the total strain rate given by [5,41,42,43,47]

$$\dot{\epsilon }={\dot{\epsilon }}_g+{\dot{\epsilon }}_c+{\dot{\epsilon }}_{gbs}$$

(10)

For consistency and simplicity in mathematics, each mechanism strain rate is assumed to take the power-law form of Eq. (6) [5,41,42,43,47]:

$$\begin{array}{c}{\dot{\epsilon }}_{gbs}=A_{gbs}\mathrm{e}\mathrm{x}\mathrm{p}\left(-\frac{Q_{gbs}}{RT}\right){\left(\frac{\sigma }{{\sigma }_y}\right)}^p\\ {\dot{\epsilon }}_g=A_g\mathrm{e}\mathrm{x}\mathrm{p}\left(-\frac{Q_g}{RT}\right){\left(\frac{\sigma }{{\sigma }_y}\right)}^n\\ {\dot{\epsilon }}_c={(1+M\epsilon )A}_c\mathrm{e}\mathrm{x}\mathrm{p}\left(-\frac{Q_c}{RT}\right){\left(\frac{\sigma }{{\sigma }_y}\right)}^m\end{array}$$

(11)

where, ${\dot{\epsilon }}_{gbs},{\dot{\epsilon }}_g,\dot{{\dot{\epsilon }}_c}$ are the creep rates of GBS, IDG and IDC, respectively; $A_{gbs},A_g,A_c$are the proportional rate constants, $Q_{gbs},Q_g,Q_c$ are the activation energies, p, n, m are the stress exponents for the respective mechanism, and M is the dislocation multiplication factor [5,41,42,43,47].

As creep deformation is an incompressible flow, the uniaxial true stress-strain evolves with the creep elongation as

$$\begin{array}{l}\epsilon =\ln (1+e)\\ \sigma ={\sigma }_a\exp (\epsilon )\end{array}$$

(12)

where e is the engineering strain, and σ_a is the initial engineering stress.

Integration of Eqs. (10-12) leads to [5,41,42,43,47]

$$\epsilon =\frac{1}{M^{*}}[\exp \left(M^{*}kt\right)-1]$$

(13a)

$$k={\dot{\epsilon }}_{s0}+{\dot{\epsilon }}_{g0}+{\dot{\epsilon }}_{c0}$$

(13b)

$$M^{*}=[p{\dot{\epsilon }}_{s0}+n{\dot{\epsilon }}_{g0}+\left(m+M\right){\dot{\epsilon }}_{c0}]/k$$

(13c)

where k is the absolute minimum creep rate, and M* is the tertiary shape factor, the subscript 0 signifies values at t = 0, and σ = σ_a.

The transient creep equation for GBS is given by [28]

$${\epsilon }_{gbs}^p=\frac{\sigma }{{\beta }^{2}H}\left[1-\exp \left(-\frac{{\beta }^{2}H{\dot{\epsilon }}_{s}t}{\sigma \left(\beta -1\right)}\right)\right]$$

(14)

where β is the transient shape factor and H is the work hardening coefficient of GBS.

Adding up Eqs. (13) and (14) results in a creep equation that depicts the entire creep curve (Stages I, II and III) as [5,41,42,43,47]

$${\epsilon }_{total}={\epsilon }_0+{\epsilon }_{tr}^P\left[1-\exp \left(-\frac{t}{t_T}\right)\right]+\frac{1}{M^{*}}[\exp \left(M^{*}kt\right)-1]$$

(15a)

with

$$\left\{\begin{array}{c}{\epsilon }_{tr}^P=\frac{\sigma }{{\beta }^{2}H}\\ t_T=\frac{\sigma (\beta -1)}{{\beta }^{2}H{\dot{\epsilon }}_s}\end{array}\right.$$

(15b)

where ε₀ is the initial elastic-plastic strain corresponding to offsetting the experimental compliance, ${\epsilon }_{tr}^P$ and $t_T$ are the primary creep strain and transient time, respectively. Actually, the primary creep time $t_{tr}^P$is measured when the primary strain ${\epsilon }_{tr}^P$ is attained in 99% such that $t_{tr}^P$ = 4.6$t_T$ [28]. Note that Eq. (15) is similar to the θ-projection equation, Eq. (2), but each component is physically- defined in terms of relevant deformation mechanisms.

From Eq. (15), the instantaneous total creep rate is given by [5,41]

$$\dot{\epsilon }=\frac{{\epsilon }_{tr}^P}{t_T}\exp \left(-\frac{t}{t_T}\right)+k\mathrm{e}\mathrm{x}\mathrm{p}\left(M^{*}kt\right)$$

(16)

It can be inferred from Eq. (16) that the creep rate would first drop from an initial value of k + ${\epsilon }_{tr}^P$/$t_T$ to some minimum value, and as the first exponential term diminishes to zero quickly, the creep rate increases exponentially with time when t > $t_{tr}^P$. In the tertiary stage, this relationship appears to be log-linear, so that both k and M* can be obtained by linear regression of the tertiary creep data. The absolute minimum creep rate k, which is actually very close to the conventionally-measured minimum creep rate, is plotted against σ_a in a log-log plot, from which the power-law proportional constant and exponent can be obtained for each mechanism. Detailed analysis steps of creep curves can be found in the book [5]. The absolute minimum creep rate vs. normalized stress relationship is shown in Figure 4 and the creep curves predicted by Eq. (15) are shown in Figure 5, respectively, in comparison with experimental data for Haynes^® 282^TM.

The DMTS model considers the concurrent operation of intragranular deformation (ID = IDG + IDC) and GBS. Naturally, the fracture mode and life will be determined based on which mechanism reaches its ductility limit first. For prediction of long-time creep performance, the initial elastic strain ε₀ and the primary strain can be negligible, because t_r >> t_T, thus, at the time of creep rupture, Eq. (15) reduces to [5,47]

$${\epsilon }_r-{\epsilon }_{tr}^P\approx \frac{1}{M^{*}}\left[\exp \left(M^{*}k{t}_r\right)-1\right]=\left\{\begin{array}{c}{\epsilon }_{gbs,cr}\\ {\epsilon }_{g,cr}\end{array}\right.$$

(17)

where ${\epsilon }_{gbs,cr}$ is critical GBS ductility, and ${\epsilon }_{g,cr}$ is critical grain ductility, ${\epsilon }_r$ is failure strain, and $t_r$ is rupture time or creep life. Note that the primary strain ${\epsilon }_{tr}^P$ arises from grain boundary hardening and it is anelastic in nature and therefore not damaging, while ${\epsilon }_{gbs,cr}$ represents the ability of grain boundaries to accommodate voids/cracks.

In the case of GBS dominance, intragranular deformation is frozen, i.e., ${\dot{\epsilon }}_{c0}=0$, ${\dot{\epsilon }}_{g0}=0$, and M*= p from Eq. (13c), then the creep life can be obtained from Eq. (17) as [47]

$$t_{r,gbs}=\frac{1}{p{\dot{\epsilon }}_{s0}}\mathrm{l}\mathrm{n}\left[1+p{\epsilon }_{gbs,cr}\right]$$

(18a)

In the case of intragranular deformation dominance, the ID-dominated creep life can be obtained from Eq. (17) as [47]

$$t_{r,i}=\frac{1}{M^{*}k}\mathrm{l}\mathrm{n}\left[1+M^{*}{\epsilon }_{g,cr}\right]$$

(18b)

Then, creep rupture occurs whichever of $t_{r,gbs}$ or $t_{r,i}$ is reached first, that is, when $t_{r,gbs}$< $t_{r,i}$, intergranular fracture occurs first, $t_r$ = $t_{r,gbs}$; and when $t_{r,gbs}$ ≥ $t_{r,i}$, fracture is caused by exhaustion of intragranular ductility, $t_r$ = $t_{r,i}$. In reality, the fracture may appear to be in a mix-mode, because both ID and GBS occur concurrently, thus the fracture surface may be composed of both grain boundary facets and regions of ductile tearing.

The creep ductility can be calculated by substituting $t_r$ into Eq. (15). The variation of creep ductility as function of stress and temperature is shown in Figure 6 for Haynes^® 282^TM. Comparing the trend in Figure 2, the GBS ductility ${\epsilon }_{gbs,cr}$ (= 0.3% for Haynes^® 282^TM) corresponds to the low-level ductility ε_cu,L, and the intragranular ductility ${\epsilon }_{g,cr}$ (= 30% for Hayne^s® 282^TM) corresponds to the high-level ductility ε_cu,U in Figure 2. At high stress, the experimental fracture strain is much higher than the low-level ductility, even though it varies a lot because of the deformation instability under high stresses. The transition between the two levels of creep ductility is a region of mix-mode fracture, with nearly equal partition of IDC and GBS. In Figure 2 the transition is attributed to diffusion-controlled cavity growth. Note that diffusional cavitation has a stress-dependence of power 1 [33], which was not the case for Haynes^® 282^TM. Fractographic examinations of crept specimens corresponding to the data points in Figure 6 showed that indeed the fracture mode transitioned from intragranular fracture to intergranular fracture as the stress went from high to low [47].

The DMTS model-predicted creep lives of Haynes^® 282^TM are shown in Figure 7 in comparison with the creep rupture data from Oak Ridge National Laboratory (ORNL) [48] and the National Research Council Canada-Carleton University-Velan collaboration (NRC-CU-Velan) [49]. The DMTS model predictions were made corresponding to each test condition, using the mechanism parameter values given in [47]. First of all, the independent tests proved the observed behavior to be true. Second, the DMTS model prediction agree very well with the experimental data in terms of LMP with C_LM = 20. A very high coefficient of determination (R² = 0.98) is obtained between the DMTS model and the experimental LMP. Unlike the irrelevance of LMP to the underlying mechanisms, the DMTS model carries the mechanism-partitioning information, which is insightful for further discussion.

Engineering design by predicting long-term creep life from short-term creep test data using the LMP - σ relation is based on the presumption that this relationship represents a one-to-one correspondence for creep rupture, but it has never been proven by theory and experiments. To investigate this matter, isothermal LMP values at each stress level are computed for different temperatures using the DMTS model, as illustrated in Figure 8. The symbols represent the indicated conditions with a rupture life < 10,000 h, which are typical of short-term tests. The solid line represents the LMP correlation (polynomial fit) for the ORNL data. Assuming that the DMTS model represents the “ideal material” as represented by the set of mechanism parameters, which has no material variability. As one can see, the isothermal short-life LMPs still fall close to the ORNL-LMP - σ line (solid line), but the mechanism-based isothermal LMP deviate from the short-term LMP best-fit line towards high LMP values, which means that the LMP - σ relation is not a one-to-one correspondence relationship. The deviated lines are governed by GBS. This is the regime where most service conditions fall in, but “predictions” are often attempted with short-term creep tests. Short-term tests at higher temperatures tend to generate a larger LMP value at a given stress level, or achieve the target LMP value at higher stresses. For example, the LMP for a service life of 100,000 h at 760 ^oC is equal to 25.8 × 10³ (with C_LM = 20). The same LMP value was produced from the short-term test with a test life of 5058.3 h under 132 MPa at 816 ^oC. If the LMP - σ relation were one-to-one, σ = 132 MPa would be taken as the design stress. However, at σ = 132 MPa and T = 1033 K (760 ^oC), the rupture life of Haynes 282 is 33,144 h, which falls quite short of the target service life of 100,000h, and hence the design is potentially dangerous. Although the short-term creep tests (< 10,000 h) would result in a rather close correlation, it should be very cautious using the empirical short-term LMP method for long-term creep life prediction where GBS plays a dominant role.

4.2. Effects of Composition and Microstructure

In general, composition and microstructure have strong effects on materials creep performance. Usually, the base metal imparts the basic physical properties such as elasticity, slip and diffusion properties, while alloying elements contribute to additional strengthening with solute atoms and precipitates/phase changes. For example, in Ni-base superalloys, this involves γ and γ’ phases and grain boundary M₂₃C₆. In heat resistant steels such as modified 9Cr1Mo steels, it involves M₂₃C₆ phase, MX phase, and Cu-rich phase. In analytical models, the effect of particle-strengthening is often considered as back-stress through either Orowan-looping or precipitate-cutting mechanisms [26,31]. Effects of grain size, grain boundary precipitates and grain boundary serrations have also been considered in creep behavior modeling [28,50]. More recent work uses molecular dynamics to simulate dislocation-particle interaction [51,52]. Molecular dynamics simulation of GBS has been done for pure aluminum without grain boundary precipitates [53]. The dynamic behavior was studied in picoseconds (ps), which is way below the engineering timeframe. As far as of engineering concerns, this article reviews a few examples of how DMTS can be used to delineate microstructure effects.

An example of microstructural effect on creep can be found for Alloy 800. A series of Alloy 800 were made with the high/low carbon content and increasing Ti+Al content in ascending order from HB647, 75193, HN823 to HH3283 [54]. The mechanism parameters for the creep rate in different batches of Alloy 800 at 550^oC are given in Table 1 [46], which are determined to match the minimum creep rate as a function of stress, as shown in Figure 9. The trend of increasing the power-law exponent for the IDG mechanism is consistent with the trend of increasing Ti+Al content that enables the intragranular γ’ precipitate strengthening mechanism. On the other hand, GBS is assumed to be the same for the high C content alloys: 75193, HN823 and HH3283, with abundant grain boundary Cr₂₃C₆ precipitates, while HB647 has a lower resistance to GBS due to its much less C content than the other batches. In these alloys, Ti and Al would have little contribution to GBS, because they do not participate in grain boundary precipitation. The impact of microstructural variation on the rupture life for Alloy 800 is shown in Figure 10. IDG is responsible for intragranular deformation instability (necking) (${\epsilon }_{g,cr}=1$) and GBS causes intergranular fracture (${\epsilon }_{gbs,cr}=0.005$).

Microstructural evolution in F91 was studied by Zhang et al. through aging heat treatment at 550^oC and 600^oC for up to 5,000 hours [55]. Notably, a considerable amount of Mo-rich Laves phase formed in F91 during the thermal exposure. The Laves phase Fe₂(Mo, W) was found mostly around Cr-rich carbide M₂₃C₆, which pinned on the prior-austenite grain boundaries and martensitic lath boundaries, and it coarsened rapidly, leading to premature creep rupture than the unaged F91. By the nature of the Laves phase formation, it would affect GBS much more than IDG. Therefore, assuming that both aged and unaged F91 have the same IDG component unaffected, subtracting the IDG rate from the total creep rate, the GBS+IDC rates are presented in Figure 11. It is seen that the Laves phase formation proportionally increased GBS, keeping the power-exponent fairly constant, since it weakens the pinning effect of M₂₃C₆/MX on lath/grain boundaries.

Weldments also represent another form of microstructural change, as compared to the base metal. The creep behavior of Inconel 740H weldment at the temperature of 760 ^oC was investigated experimentally and analytically using DMTS model [56]. Inconel 740H is the only material approved by ASME for advanced ultra-supercritical (A-USC) steam turbine applications [57]. Its weldment is as important as the base metal in such applications. The Inconel 740H weldment specimens were prepared with the gas tungsten arc welding (GTAW) technique. In comparison with the base metal, the 740H weld material (after a post weld heat treatment) had apparently wider grain boundaries with coarsened/elongated grain boundary precipitates and γ’ denuded zone. This would potentially weaken the resistance to GBS. In contrast, weldment of Haynes^® 282^TM (after post-weld heat treatment) has no γ’ denuded zone along the grain boundaries, so its creep properties are very similar to the base metal [49]. The creep data generated on Inconel 740H weldment at 760 ^oC were analyzed using DMTS model to delineate intragranular deformation and GBS, and the activation energies for the respective mechanism were assumed to be the same as Haynes 282, a similar alloy, then the DMTS model is used to predict the creep rupture lives within a temperature range of 700 – 800 ^oC, as shown in Figure 12. The creep life predictions of the DMTS model for Inconel 740H weldment agree very well with creep rupture test data from both literature [58] and NRC-CU-Velan in-house investigation [56]. Although only a few short-term creep data at 760^oC are used to calibrate the mechanism parameters of the DMTS model, the agreement between the DMTS model and test data is remarkably good. The turning points on the DMTS model curves correspond to the transition from intragranular fracture to intergranular fracture, as the stress goes from high to low. The fractured surfaces and longitudinal sections of creep-tested Inconel 740H weldment specimens were examined using SEM, which corroborated the DMTS model inference that the creep failure of Inconel 740H weldment was dominated by GBS with an intergranular fracture mode [56].

The above examples demonstrate that using the mechanisms-delineation approach such as the DMTS approach, the effects of composition and microstructure of the materials can be reconciled with the physical deformation mechanisms over a wide range of stress and homologous temperature. This approach removes the ambiguities associated with empirical equations and links the material behavior closely to the physical metallurgy of materials, which should be the focus of future creep studies.

4.3. Environmental Effects

Most high-temperature components operate in open air, which would inevitably suffer from oxidization. The question is to what extent oxidation may affect the material’s long-term creep performance? Bueno conducted short-term creep tests on 2¼Cr-1Mo steel in vacuum and air, which showed that oxidation had an effect of increasing the creep rate in air as compared to that in vacuum [59]. Such oxidation-free creep data are rarely available for existing materials, because long-term vacuum creep tests are cost-prohibitive. In conventional creep tests, creep data are obtained with coupon-borne influence of oxidation, but oxidation is a time-dependent process, therefore, empirical extrapolation of short-term creep data for long-term creep life prediction without considering oxidation is questionable. For example, Wheaton (1965) found that MarM 509 alloy exhibited a higher creep rate ~1×10^-6 h^-1 but a longer creep life of 983 h at 172.4 MPa /871^oC, compared to the creep rate of 5×10^-7 h^-1 and rupture life of ~250 h at 48.3MPa/1093^oC [60]. The Monkman-Grant relation would not work for this phenomenon. When the critical GBS ductility ${\epsilon }_{gbs,cr}$ was modified with consideration of oxidation, the failure time was correctly predicted [42]. Dyson and Osgerby (1995) considered creep-environment interaction using the continuum damage-mechanics (CDM) approach [61]. However, CDM does not attribute the effect of oxidation to specific deformation mechanisms, that is, CDM does not distinguish intragranular deformation or GBS, so the treatment is still ambiguous.

Oxidation damage can be assessed through post-mortem examination of creep-fracture coupons. In general, the growth of oxide scale, δ, follows the parabolic law as expressed by [62]

$$\delta =\sqrt{2k_{ox}t}$$

(19)

where k_ox is the oxidation coefficient.

When an oxide scale forms, if it is intact, the stress in the oxide scale is given by the composite rule as [5]

$${\sigma }_{ox}=\frac{E_{ox}}{f{E}_{ox}+(1-f)E_s}[\sigma +(1-f)E_s({\alpha }_s-{\alpha }_{ox})(T-T_0)]$$

(20)

where α_ox and α_s are the coefficients of thermal expansion (CTE), E_ox and E_s are the elastic moduli of the oxide and substrate, respectively; f is the volume fraction of the oxide scale. For a round bar coupon, f = 2δ/πr, r is the cross-sectional radius. It is noticed that, even when a very thin film forms with f ≈ 0, the stress in the oxide film is magnified by E_ox/E_s, while the stress in the substrate remains nearly the same. Thus, the oxide scale tends to break first.

If the oxide scale breaks, the load bearing area of the coupon (component) is effectively reduced, and the true stress becomes [62]

$$\sigma =\frac{{\sigma }_a}{1-f_{ox}}\mathrm{e}\mathrm{x}\mathrm{p}\left(\epsilon \right)$$

(21)

Then, following the same derivation procedure of the DMTS model, the same creep strain equation, Eq. (13) can be obtained but with oxidation-modified parameters [62]

$$k=\frac{1}{{(1-f_{ox})}^p}{\dot{\epsilon }}_{s0}+\frac{1}{{(1-f_{ox})}^n}{\dot{\epsilon }}_{g0}+\frac{1}{{\left(1-f_{ox}\right)}^m}{\dot{\epsilon }}_{c0}$$

(22)

$$M^{*}=[\frac{p}{{(1-f_{ox})}^p}{\dot{\epsilon }}_{s0}+\frac{n}{{(1-f_{ox})}^n}{\dot{\epsilon }}_{g0}+\frac{m+M}{{\left(1-f_{ox}\right)}^m}{\dot{\epsilon }}_{c0}]/k$$

(23)

The above equations reduce to the basic DMTS model when oxidation is absent, and hence it is called the oxidation-modified DMTS (O-DMTS) model. Eq. (22) is shown to predict the creep rate increasing with formation and breaking of oxide scale, which agrees with the experimental observation for F91 [62].

For the purpose of comparison, the predicted creep lives of NIMS Plate MgB by the O-DMTS model are plotted in the LMP diagram in comparison with the experimental data (the rupture lives are taken to be 90% time to rupture (TTR)), because necking occurred at last before rupture. The Larson-Miller constant, C__LM, is 33 [63]. Figure 13 shows that the predicted lives agree very well with the experimental data from the NIMS database which contains creep rupture life data up to 10⁵ h. While the LMP method collapse the creep rupture data, the relationship between LMP and stress is not as simple as discussed in section 2.2, one would have to use polynomials to draw a line (σ vs. LMP relationship) between the data. The problem with polynomial fitting is that it cannot be established without a large amount of data, and extrapolation of polynomials is usually invalid beyond the fitted range. Therefore, long-term creep test data are still needed, when using the LMP method. Interestingly, the oxidation-modified DMTS (O-DMTS) creep model predicted lives also fall well with LMP, but it is drawn from fundamental deformation mechanisms. Furthermore, the O-DMTS model provides a quantitative partition of GBS/IDG/IDC mechanisms, and it also predicts the rupture mode. Figure 14 shows two microstructures at creep rupture under 600℃/160 MPa and 650℃/100 MPa, along with the mechanism pie-chart made up by the contributions from the three mechanisms to the total strain rate. In the 600℃/160 MPa case, GBS was predominant (~94%) and the microstructure retained its original lath structure; whereas in the 650℃/100 MPa case, as the model predicted that a significant portion (~38%) of creep deformation would occur by IDC, the lath structure indeed became more elongated (note that the O-DMTS prediction of mechanism partitioning ratio was for the initial condition, and IDC would take more dominance as creep proceeded to the last point of fracture as the true stress increases with elongation). The model predication was thus supported by the metallurgical evidence. Using the O-DMTS model, such pie-charts can be drawn for every creep condition. By this way, engineers are able to understand the physics of failure and identify the failure mode in quantitative details, whereas the LMP plot does not provide mechanism information, from life prediction point of view.

5. Concluding Remarks

This review summarizes creep phenomena including creep strain, creep rate, creep ductility, creep life and fracture mode, and uses a mechanism-delineation approach via the deformation-mechanism-based true-stress (DMTS) creep model and the oxidation modified DMTS (O-DMTS) creep model to describe/interpret the experimental observations. Over the past 110 years, these phenomena were treated with various empirical equations, but those equations do not connect to each other so that they do not provide a complete picture of creep. Especially, the empirical long-term creep life prediction was not warranted for the conditions where the mechanism(s) deviate from the tests from which the data for extrapolation were generated. The main points drawn from this review are summarized below.

As a general creep model for metallic materials, the DMTS model takes the contributions of dislocation mechanisms such as GBS, IDC and IDG into consideration, which enables describing all salient creep phenomena including creep strain, creep rate, creep ductility, creep life and fracture mode all together, as demonstrated for complex engineering alloys investigated in this paper and in the references. By the general validity of the physical mechanisms involved, the DMTS model is applicable to all metallic materials.

The DMTS model interprets the microstructural effect with respect to the mechanism contribution that may be affected by the change in microstructure (grain size, intragranular/grain boundary precipitates, etc.) either through compositional changes, or thermomechanical processing and welding, or long-term aging. This point of view is not provided by other models. Microstructural evolution would especially influence long-term creep performance of materials. However, it lacks systematic and quantitative information to cover all service conditions for engineering alloys. More studies should be carried out for high-temperature structural materials over a wider range of stress and temperature, incorporating quantitative description of phase transformation and growth. Linking microstructural evolution with deformation mechanisms should be the focus of future creep studies.

The oxidation-modified DMTS (O-DMTS) creep model further takes into account the effect of oxide scale formation on the true stress during creep process. This model gives a quantitative description of the environmental effect separated from the material-intrinsic deformation mechanisms: IDG, IDC and GBS. It then predicts long-term creep performance with true oxidation contribution, as opposed to other empirical methods where short-term data are extrapolated without oxidation in consideration. As a demonstration, the O-DMTS model is validated with the NIMS creep data for modified 9Cr-1Mo steels. The model shows promising for long-term creep life prediction and failure mode using short-term creep test data with the identified deformation mechanisms, which is also unique compared to other creep models. Again, in principle, this model can be applied to all metallic materials which experience IDG, IDC and GBS as common deformation mechanisms and suffer oxidation at high temperatures.

Using the DMTS/O-DMTS creep model, the Larson-Miller parameter (LMP) method has been analyzed and furbished with mechanism partition information, providing insights into the controlling deformation mechanisms and potential failure mode. The DMTS model is applied to construct the isothermal LMP - σ curves of 700 ℃, 760 ℃, 788 ℃, and 927 ℃ for Haynes 282. While the LMP approach may be used to construct an experimental best-fit line for short-term creep rupture (< 10⁴ h), the LMP results of long-term creep rupture deviate from the short-term correlation, and the disparity becomes larger towards large LMP. This implies that the LMP method for long-term creep life prediction based on the short-term data has limitations. The DMTS model certainly amends the traditional LMP for extrapolation to long-term creep life. This is a significant improvement to the state-of-the art creep design methodology. By providing mechanism-informed and reliable creep life prediction, it saves a significant testing effort to generate massive experimental data as the Larson-Miller parameter method would require.