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.
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]
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]:
where,
are the creep rates of GBS, IDG and IDC, respectively;
are the proportional rate constants,
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
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]
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]
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]
with
where ε
0 is the initial elastic-plastic strain corresponding to offsetting the experimental compliance,
and
are the primary creep strain and transient time, respectively. Actually, the primary creep time
is measured when the primary strain
is attained in 99% such that
= 4.6
[
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]
It can be inferred from Eq. (16) that the creep rate would first drop from an initial value of k +
/
to some minimum value, and as the first exponential term diminishes to zero quickly, the creep rate increases exponentially with time when t >
. 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 ε
0 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]
where
is critical GBS ductility, and
is critical grain ductility,
is failure strain, and
is rupture time or creep life. Note that the primary strain
arises from grain boundary hardening and it is anelastic in nature and therefore not damaging, while
represents the ability of grain boundaries to accommodate voids/cracks.
In the case of GBS dominance, intragranular deformation is frozen, i.e.,
,
, and M*= p from Eq. (13c), then the creep life can be obtained from Eq. (17) as [
47]
In the case of intragranular deformation dominance, the ID-dominated creep life can be obtained from Eq. (17) as [
47]
Then, creep rupture occurs whichever of or is reached first, that is, when < , intergranular fracture occurs first, = ; and when ≥ , fracture is caused by exhaustion of intragranular ductility, = . 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
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
(= 0.3% for Haynes
® 282
TM) corresponds to the low-level ductility ε
cu,L, and the intragranular ductility
(= 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
2 = 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
3 (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
23C
6. In heat resistant steels such as modified 9Cr1Mo steels, it involves M
23C
6 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
23C
6 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) (
) and GBS causes intergranular fracture (
).
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
2(Mo, W) was found mostly around Cr-rich carbide M
23C
6, 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
23C
6/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
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]
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]
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]
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]
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
5 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.