Search tips
Search criteria 


Logo of materialsLink to Publisher's site
Materials (Basel). 2013 March; 6(3): 1118–1137.
Published online 2013 March 20. doi:  10.3390/ma6031118
PMCID: PMC5512967

Recent Advances in Creep Modelling of the Nickel Base Superalloy, Alloy 720Li


Recent work in the creep field has indicated that the traditional methodologies involving power law equations are not sufficient to describe wide ranging creep behaviour. More recent approaches such as the Wilshire equations however, have shown promise in a wide range of materials, particularly in extrapolation of short term results to long term predictions. In the aerospace industry however, long term creep behaviour is not critical and more focus is required on the prediction of times to specific creep strains. The current paper illustrates the capability of the Wilshire equations to recreate full creep curves in a modern nickel superalloy. Furthermore, a finite-element model based on this method has been shown to accurately predict stress relaxation behaviour allowing more accurate component lifing.

Keywords: creep, stress rupture, stress relaxation, superalloy

1. Introduction

Nickel-base superalloys perform a key role in gas turbine aero engines due to their superior mechanical properties at elevated temperatures and good corrosion resistance. This high temperature strength is usually provided by a distribution of γ’(Ni3(Al,Ti)) precipitates which hinder dislocation movement. The polycrystalline superalloy, Alloy 720Li, is an example of a workable alloy used in the manufacture of turbine discs. In this case a heat treatment procedure has been devised to provide a microstructure consisting of primary γ’(1–10 μm), secondary γ’(70–120 nm) and tertiary γ’(5–50 nm depending on whether the material is quenched or aged).

Turbine components are subjected to high temperatures and high centripetal force which limit their life. However, because they are safety critical, it is imperative that accurate predictions can be made for the safe operating lives of these components. Whilst the operating life of the disc will be fatigue dominated, in order to affect an accurate fatigue life prediction it is critical that the redistribution of stresses at high temperatures be well understood. This stress redistribution will be as a result of creep deformation at the high temperatures experienced in the component, particularly at stress raising features such as blade loading slots, which experience the highest temperatures due to their proximity to the disc rim, and hence the hot gas stream.

Until recent years, creep models have been based on the traditional assumptions of creep deformation, in particular the accuracy of power law type approaches and the application of fracture mechanism maps defining regimes where specific deformation types (i.e., dislocation or diffusional creep) occur. However, the limitations of these methods have been exposed, with concerns raised particularly about issues such as power law breakdown at high stresses and the lack of accuracy of power law techniques such as the Larson-Miller approach [1].

More recent techniques such as the Theta projection method [2,3], Hyperbolic tangent [4] and Wilshire equations [5,6,7,8,9] have offered alternative approaches, with the hyperbolic tangent and Wilshire equations methods showing particular aptitude in fitting stress-rupture data.

In particular the Wilshire equations require further investigation due to their wide range of early successes. These include accurate predictive methodologies for copper [5], aluminium alloys [6] and a wide range of steels [7,8,9]. Also, more recently a process of full creep curve derivation based on numerical analysis of selected times to strain [10,11] has shown promise as an improvement on the Theta projection curve fitting technique. The equations also benefit from close correspondence of observed trends with physical phenomena, i.e., breaks seen in the equations regularly occurring at the yield stress of the material. Such correlations add confidence to the technique and aid in the application of the equations to a wide range of materials.

In the aerospace sector, it is the focus on accurate modelling of the creep curve and in particular times to critical strains which is considered most necessary, since extrapolation of data to long lives is irrelevant given the dominance of fatigue processes over this time period. The current submission offers an interesting perspective of the applicability of the Wilshire equations, since their application to precipitate hardened superalloys has been limited, and certainly creep curve shape has not been previously considered.

2. Experimental Data

Alloy 720Li (Low interstitial) was developed for gas turbine disc components by modifying the chemical composition of Alloy 720 [12]. The more recent Alloy 720Li has reduced Cr content to suppress sigma phase formation and C and B contents are reduced to improve forgeability. The chemical composition is of the alloy was (wt%) 16.3Cr, 14.7Co, 3.00Mo, 1.31W, 5.02Ti, 2.57Al, 0.011C, 0.026Zr, 0.015B with balance Ni.

In order to produce a homogenous microstructure of γ’, standard thermomechanical treatments were used. Post heat treatment, average primary, secondary and tertiary γ’ particles sizes of 5 μm, 100 nm and <10 nm respectively were present. The alloy had a total volume fraction of γ’ particles of ~30%. The resultant microstructure of the material, shown in Figure 1 had a grain size of approximately 10 μm.

Figure 1
Microstructure of Alloy 720Li.

Creep tests were performed at constant stress at a range of temperatures between 550 °C and 750 °C with applied stresses chosen to give rupture lives from 5 h to 5,000 h. Specimen strain was recorded for the duration of each test. A stress relaxation test was completed at 650 °C in air [13]. For this test, specimen strain was held at 0.5% for 20 h, followed by an increase in strain to 0.7% for a further 20 h. Specimen load was monitored throughout the test.

Tensile strength, σTS, values were obtained from tensile tests at each creep temperature (550–750 °C). These tests were performed under strain control using a constant strain rate of 5 × 10−5 s−1 increasing upon yielding to a faster rate of ~0.002 s−1 until failure.

3. Results and Discussion

For all creep tests, normal creep curves were recorded. On loading, an initial strain (ε0) is observed, followed by primary creep where creep rate (ε˙=dε/dt) decreases until a minimum value is obtained (ε˙m). After ε˙m, the creep rate accelerates during tertiary creep until failure. Rupture time (tf) and strain to failure (εf) were recorded and utilised in the plotting of traditional creep data analysis. Rupture time decreases with increasing ε˙m (Figure 2) and the data can be represented using the Monkman-Grant relation [14]


where the product, M, is the Monkman-Grant constant. For Alloy 720Li a value of M0.272 was obtained.

Figure 2
Relationship of minimum creep rate (ε˙m) with rupture time (tf) for Alloy 720Li (550 °C to 750 °C).

In line with common practice, and using the data shown in Figure 2, the stress (σ) and temperature (T) dependence of ε˙m and tf can be calculated using the power law


where R = 8.314 Jmol−1K−1. A and n are material parameters dependent on stress and temperature. The stress exponent n varies from n14 for the tests with conditions giving short lives to n5 for the long duration tests. Qc is the activation energy for creep with a value of >400 kJ mol−1. However, this value also varies with test conditions. The stress and temperature dependence of ε˙m and tf can be displayed on single curve by normalising σ by the ultimate tensile stress σTS of the material. Therefore Equation (2) becomes


where A*A and Qc*Qc; Qc* is lower than Qc since σTS partially allows for the dependence of creep properties on temperature. For Alloy 720Li, a value of Qc* 330 kJ mol−1 was established. Values for σTS obtained from tensile test data at each creep temperature along with values of 0.2% proof stress, σY [15] are shown in Table 1.

Table 1
Tensile properties of Alloy 720Li (550 °C to 750 °C).

The dependence of temperature compensated ε˙m with σ/σTS is shown in Figure 3. This plot shows that n in Equation (3) can not be considered constant and varies with applied σ. The variability in n values has been attributed to varying creep mechanisms. This inconsistency in n makes predicting creep properties over a wide range of test conditions difficult. For example the ability to predict long term creep properties from short term test data is very important in reducing the duration of component design, and is essentially impossible since no extrapolation can be robustly predicted due to the unknown curvature of the graph. Other methods have attempted to address this issue, such as the Wilshire equations [5,6,7,8,9], which have proved effective in extrapolation of creep properties over a range of creep conditions.

Figure 3
Power Law representation of Alloy 720Li (550 °C to 750 °C).

3.1. Extrapolation of Creep Data

For any material, applying a stress equal to σTS at any given temperature will result in near instantaneous rupture. Furthermore, creep rupture cannot be expected if the applied stress is zero. Therefore creep properties must be considered over the range 0σσTS. Creep prediction methods must allow for the fact that ε˙m and tf0 as (σ/σTS)1. Similarly, the following must be true, as (σ/σTS)0, ε˙m0 and tf. The Wilshire equations make use of the sigmoidal shape of the function f(σ)=ln(ln(σ/σTS)) and have been shown to accurately predict creep properties over a range of creep conditions. These equations relate temperature compensated tf and ε˙m to σ/σTS using




respectively. The activation energy for creep (Qc*) for the Wilshire equations is equivalent to that obtained for the power law Equation (3) where Qc* 330 kJ mol−1; u, k1, v and k2 are material constants obtained from experimental data. Plotting ln(ln(σ/σTS)) against the natural logarithm of temperature compensated rupture times (ln(tf exp (−Qc/RT)) gives a straight line of gradient u and y-axis intercept of lnk1 (Figure 4). For Alloy 720Li a different set of u and k1 values were obtained above and below σ/σTS0.65 indicated by the lines in Figure 4. This value is similar to the yield stress of the alloy and is consistent with behaviour observed in other creep resistant alloys [4,5,6,7,8]. The values of u and k1 obtained from stresses corresponding to σ/σTS < 0.65 were 0.184 and 105.53 respectively whereas at stresses greater than the break point values of 0.0993 and 8.471 respectively were obtained (Table 1). These values produce stress rupture master curves which represent experimental behaviour over a range of conditions (Figure 5). The quality of fit can be displayed by plotting creep lives predicted using the above method to those observed experimentally (Figure 6). This figure shows a good quality fit over a wide range of test lives.

Figure 4
Dependence of ln(ln(σ/σTS)) on ln(tfexp(Qc/RT)) for Alloy 720Li, with Qc* [congruent with] 330 kJ mol−1.
Figure 5
Stress rupture plots for Alloy 720Li with stress normalised by σTS predicted using Equation (4) (550 °C to 750 °C).
Figure 6
A comparison between experimental creep rupture lives for Alloy 720Li and those obtained using Equation (4) (550 °C to 750 °C). Dashed lines represent a deviation of a factor +/−2.

A similar approach is used for predicting minimum creep rates. v and k2 are obtained from the gradient and intercept of a best fit line through a plot of ln(ln(σ/σTS)) against ln(ε˙mexp(Qc/RT)) (Figure 7). Again, a change in the value of the coefficients is observed above and below σ/σTS0.65. At stresses corresponding to σ/σTS<0.65, values of v = −0.154 and k2 = 23.58 are obtained. At stresses corresponding to σ/σTS>0.65, values of u = −0.064 and k2 = 2.223 are obtained (Table 1). Figure 8 shows minimum creep rates predicted using this method and Figure 9 illustrates that once more, this method gives a good quality representation of the test data.

Figure 7
Dependence of ln(ln(σ/σTS)) on ln(ε˙mexp(Qc/RT)) for Alloy 720Li, with Qc* [congruent with] 330 kJ mol−1.
Figure 8
Minimum creep rate plots for Alloy 720Li with stress normalised by σTS. Solid lines represent the fit obtained using Equation (5) (550 °C to 750 °C).
Figure 9
A comparison between experimental minimum creep rates for Alloy 720Li and those obtained using Equation (5) (550 °C to 750 °C). Dashed lines represent a deviation of a factor +/−2.

3.2. Creep Deformation

The equations above represent a current method for extrapolating tf and ε˙m. However an understanding of these values alone is not sufficient for component design since they do not quantify the full shape of a creep curve. This is illustrated in Figure 10 which shows 3 distinct creep curves with identical values of ε˙m, εf and tf. For design purposes it is often necessary to calculate time to a given strain or to predict stress relaxation, both of which require a method capable of predicting creep curve shapes over a wide range of stress and temperature.

Figure 10
Schematic representation of different creep curves with equal values of ε˙m, εf and tf.

The creep curves obtained for Alloy 720Li all display a period of primary creep where creep rate decreases subsequent to loading, followed by tertiary creep where creep rate accelerates to failure, however the proportion of primary to tertiary creep varies with applied test conditions (Figure 11).

Figure 11
Normalised creep curves showing the variation in the proportion of primary to tertiary creep with test conditions.

The creep curves demonstrate a smooth transition from primary to tertiary creep consistent with the model proposed by Evans [2,3] as opposed to “steady state” creep models with associated secondary creep. Although these models can be used to analyse creep curve shapes, there is a need to quantify creep curve shape. Each creep curve, after removing the strain on loading, may be broken down into primary creep strain, εP, tertiary creep strain, εT, and strain at failure, εF, where


Strain to failure is dependent on applied stress and decreases as (σ/σTS)1 (Figure 12). εP can be approximated by extrapolating ε˙m back to t = 0 (Figure 13). Comparing εP to applied stress helps to quantify the difference in curve shape observed in Figure 11. Below about 700 MPa, the primary component of creep strain is negligible. By plotting εP against stress normalised by tensile stress some of the effects of temperature can be offset (Figure 14). This plot shows that below 0.4σ/σTS, εP is negligible and above 0.6σ/σTS this value increases rapidly, however for all tests εP is small compared to εT and hence εF ≈ εT (Figure 15). These observations show that creep rupture is dominated by tertiary creep and hence rupture can be characterised by micromechanical deformation and damage processes that occur during tertiary creep in Alloy 720Li.

Figure 12
Dependence of strain to failure, εF against normalised stress for Alloy 720Li at 550 °C to 750 °C.
Figure 13
Description of a creep curve using εP, ε˙m, tf and εF.
Figure 14
Dependence of primary creep strain, εP against normalised stress for Alloy 720Li at 550 °C to 750 °C.
Figure 15
Dependence of total creep strain to failure, εF on tertiary creep strain for Alloy 720Li at 550 °C to 750 °C.

Many methods have been proposed to predict creep curve shapes from early experiments on fuse wire by Andrade [16] and methods describing creep deformation of pipe flanges proposed by Bailey [17]. Many methods use a function to describe the shape of the creep curve and relate the parameters of this function to applied conditions. The θ-projection method [2,3] is one such technique which extends the work of McVetty [4] and Garofalo [18]. This method uses the function


to describe the accumulation of creep strain, ε, with time, t, where θ1–4 are material constants dependent on stress and temperature. The parameters θ1 and θ2 describe primary creep and θ3 and θ4 represent tertiary creep. The function


is then used to evaluate θ1–4 with stress and temperature where n = 1 −4. However, this function is not valid for the full range of creep conditions 0σσTS. A constitutive method based on the θ- projection method has been formulated to predict creep rate (ε˙c) over varying stress and temperature based on evolution of internal material state variables [19]. An alternative method for predicting the evolution of creep rate with time is the continuum damage method proposed by Kachanov [20] and later modified by Rabotnov [21]. Othman and Hayhurst [22] extended this method to account for the effects of primary creep.

3.3. Extrapolation of Creep Curves

A more recent approach of predicting creep curve shapes has been proposed based on the Wilshire equations [10,11]. This method consists of relating times to given strains (tε) to applied conditions using an equation similar to Equation (4)


w and k3 are material constants dependent on ε and are obtained by plotting ln (ln(σ/σTS)) against ln (tε exp (−Qc/RT)) for times to finite strains using a similar method as used to calculate u and k1 (Figure 16).

Figure 16
Dependence of ln(ln(σ/σTS)) on ln(tε=1%exp(Qc/RT)) for Alloy 720Li, with Qc* [congruent with] 330 kJ mol−1.

The activation energy used is equal to that found for Equations (4) and (5) (Qc* 330 kJ mol−1). From these values a single value for w is chosen to represent all strain levels, and values of k3 are optimised for each ε using this constant w. Relating k3 to ε allows creep curves to be generated by constructing t/ε pairs. Currently, the proposed method relates k3 to ε using


where k3,0, k3,1 and k3,2 are material constants. Once the constants are known, Equation (9) can be used to calculate “life to critical strain” and stress/time plots for individual values of creep strain can be calculated. A stress/time plot for 1% creep strain is shown in Figure 17 and the quality of fit is displayed in Figure 18.

Figure 17
Time to 1% creep strain, tε = 1%, for Alloy720Li predicted using Equations (15) and (16) (550 °C to 750 °C).
Figure 18
A comparison between experimental times to 1% creep strain for Alloy 720Li and those obtained using Equation (9) (550 °C to 750 °C). Dashed lines represent a deviation of a factor +/−2.

Figure 19 shows an isothermal surface displaying the relationship between stress, strain and time predicted by the model. This method can be used to predict whole creep curves over a wide range of conditions. Figure 20 shows a comparison between predicted and experimentally obtained creep curves.

Figure 19
Dependence of creep strain, εc, on stress, σ, and time for Alloy 720Li at 650 °C as predicted by the Wilshire creep curve method.
Figure 20
Predicted and experimental creep curves for Alloy 720Li at 650 °C and 700 °C.

3.4. Stress Relaxation

For design purposes the ability to predict creep behaviour in engineering components is essential. One of the ultimate goals of the work was develop a method to accurately predict stress relaxation behaviour, in order to improve the accuracy of fatigue lifing of components based on initial stress conditions. To achieve this, a finite element model was derived based on the aforementioned creep curve method. An important aspect of modelling creep behaviour is accounting for prior creep deformation. Using the method described above, the shape of a creep curve can be accurately predicted for a given applied stress and temperature. The creep rate depends on the position on the creep curve with respect to time (or strain).

Methods used to predict creep behaviour may be based on time hardening, strain hardening or life-fraction hardening. Time fraction hardening is the simplest concept where creep rate is calculated based on the creep rate at time equal to tTOT, the total analysis time. This method is only applicable to cases where the stress state remains relatively constant. The strain hardening model is more suitable in instances where the stress varies considerably during an analysis, however, this method may produce inaccuracies if the conditions vary from those displaying a primary dominated creep curve to those displaying a tertiary dominated creep curve or vice versa. Life-fraction hardening attempts to address this issue by calculating creep rate based on effective time, tEFF, where tEFF = w.tF. w is life fraction ratio (t/tF) calculated at the end of each increment. This method has the advantage of predicting creep rupture when w exceeds 1.

The creep rate for any given increment can be calculated by the difference in effective creep strain at the beginning and end of the time increment, and the duration of the time increment, Δt


For strain hardening, the magnitude of creep strain at the beginning of the time increment, εi, is known and so a value of effective time for the beginning of the increment, ti, can be calculated using Equations (9) and (10). The value of effective time at the end of the increment, ti+1, is simply ti+Δt. The value of creep strain at the end of the time increment, εi+1, cannot be calculated directly from time so an iterative method must be used. For this an initial estimate of εi+1 is required. For the first increment, this can be set to an arbitrary value (ε0 = 0.001). For subsequent increments the value of strain from the previous increment can be used. From this value, an initial estimate k3 can be made using Equation (10) and then using Equation (9), an initial value of time, t0, can be calculated. A Gauss-Newton based iterative method is used to find the value for εi+1 that minimises the squared error between effective time at the end of the increment, ti+1, and predicted time, tn.




Forward difference is used to calculate f(ε). A value of εi+1 is obtained when |titn| is less than a time error tolerance, TETOL or if the number of iterations, n, required to calculate εi+1 exceeds a set value. Using this value for εi+1, creep rate, ε˙, can be calculated using Equation (11).

The calculation of creep rate using time hardening is more complex since both εi and εi+1 in Equation (11) cannot be calculated directly. Therefore an iterative method is used to determine εi using Equation (12) with


where ti is total analysis time at the beginning of the increment. The same iterative method is used to calculate εi+1 determined using Equations (9) and (10) akin to strain hardening where ti+1=ti+Δt. The formulation for life fraction hardening is similar to that of time fraction hardening however, ti is substituted with


where tf is creep rupture time calculated using Equation (4) and w is the life-fraction. w is simply t/tf calculated at the end of each increment based on creep conditions for that increment.

The method outlined above has been included in user defined material subroutine (*CREEP) for use with Abaqus finite-element analysis software and applied to the prediction of stress relaxation behaviour for the same alloy. An axisymmetric simulation of the full specimen was performed using linear elastic behaviour and isotropic plasticity data obtained from tensile testing at 650 °C. A comparison of predicted and experimentally observed stress relaxation data for Alloy 720Li at 650 °C is shown in Figure 21. It can be seen that the stress relaxation behaviour predicted using the Wilshire based creep model closely resembles test data, providing confidence in this method for describing stress conditions at component features. Clearly, accurate descriptions such as this provide confidence and would impact favourably on the total life prediction of the component.

Figure 21
Comparison of predicted and experimentally obtained stress relaxation data.

4. Discussion

The stress and temperature dependence of the creep properties of Alloy 720Li predicted using a traditional power law approach are expressed by the stress exponent n and the activation Qc in Equation (2). For medium duration tests (250–5000 h) a stress exponent of ~5 was obtained consistent with values obtained for a number of metals and alloys [23]. However, at higher stresses (for shorter duration tests the stress exponent n14 introducing doubt that Equation (2) can be used to predict creep properties over the full range of test conditions. Furthermore, the activation energy for creep, Qc, was found to be in excess of 400 kJ mol−1 and also varied with test conditions. A lower activation energy for creep, Qc* = 330 kJ mol−1 was found when temperature compensated minimum creep rate, ε˙m was compared to σ/σTS since σTS varies with test temperature. This value is greater than the activation energy for lattice diffusion in nickel (~280 kJ mol−1) and significantly greater than that of grain boundary self diffusion (~115 kJ mol−1) [24].

The Wilshire approach extends this concept by relating creep rupture life and minimum creep rate to ln(σ/σTS) using Equations (4) and (5) respectively. This method limits predictions of creep properties over the range 0>σ>σTS. Using this approach, creep properties for the whole range of test conditions can be predicted over two distinct regimes separated by a break at 0.65σTS. At lower stresses values of u = 0.184 in Equation (4) and v = 0.154 were obtained. Inverting these values gives 1/u5.5 and 1/v6.5 respectively, closely corresponding to values of n for power law creep. At stresses greater than 0.65σTS values of u and v of 0.0993 and 0.064 respectively were obtained. These values are consistent as σ → σTS, stress levels where power law breakdown is observed.

Previous studies on polycrystalline copper [5], steels [25,26,27] and titanium alloys [11] have shown that a change in slope in Figure 4 and Figure 7 and hence a corresponding changing in the parameters u, v, k1 and k2 can be attributed to differences in creep behaviour above and below a level of stress related to yield stress of the alloy. For each creep temperature the 0.2% proof stress 0.72σTS, so it can be assumed that yield occurs at slight lower stresses, consistent with the observed changes in u and v. Also, analysis of the creep curve shapes at different test conditions gives insight into the mechanisms of creep. The relationship of failure strain, εF, on σ/σTS (Figure 12) shows that εF decreases as stress increases, with creep ductilities < 10% at stresses greater than 0.65σTS. Furthermore, Figure 14 shows that εp is negligible below 0.4σTS but inceases rapidly above 0.6σTS. Since primary creep is characterised by intergranular dislocation work hardening processes it can be assumed that dislocation multiplication occurs within the grains above 0.6σTS. Therefore, subsequent creep must occur in grains with higher dislocation densities resulting in a change in material response consistent with previous studies [5]. At lower stresses, below 0.6σTS, it is assumed that primary creep is limited to slip in suitable orientated grains where local stresses exceed the critical resolved shear stress [28] combined with the movement of pre-existing dislocations within the material under the applied stress.

At all applied test conditions Alloy 720Li behaves as a creep-ductile material, since εF/εP 1 however this ratio decreases as σσTS. Goodall et al. [29] proposed that creep-ductile materials are capable of accommodating significant tertiary creep strain prior to failure when εF/εP10. The ratio εF/εP for Alloy 720Li decreases as σ/σTS increases and significantly εF/εP10 when σ/σTS>0.65 (Figure 22). Therefore, during creep in the high stress regime, the higher dislocation density results in a higher probability of stress rupture occurring before substantial tertiary creep.

Figure 22
Dependence of εF/εP against normalised stress for Alloy 720Li at 550 °C to 750 °C.

5. Conclusions

Minimum creep rate and creep life properties for Alloy 720Li obtained from uniaxial creep testing can be represented using a power law approach (Equation (2)). However, due to the variability of the parameters A and n, this method is not effective at extrapolating to longer test lives. An Arrhenius term is used to compensate for temperature and if the applied stress is normalised by σTS, the test data can be superimposed onto a single curve.

The Wilshire equations method (Equations (4) and (5)) offers an alternative approach and has been shown to extrapolate to long term data favourably compared to other methods for other alloys. This approach highlights a predictable change in behaviour above and below 0.65σ/σTS which is consistent across a range of materials, indicating a change in creep behaviour related to σY. In applying the equations to the current material it is clear that the creep properties of the material (time to rupture and minimum creep rate behaviour) are well represented by the equations.

Interestingly, values of primary creep strain, εp, were shown to decrease to essentially zero, at σ0.4σTS. Above this stress, but below σY, it is probable that a limited amount of slip in suitably orientated grains will result in small amounts of primary creep strain. Above 0.65σ/σTS it is assumed that dislocation multiplication within the grains during primary creep will cause a change in the bulk creep properties of the material, demonstrated by a change in parameters u, v, k1, and k2 in Equations (4) and (5).

However, the value of the equations in the current application lies not in life extrapolation, but in their ability to represent the holistic creep curve. The Wilshire method has recently been extended to allow predictions of ‘times to strain’ to be made, thus recreating creep curve shapes using Equations (9) and (10). A variation of constants above and below 0.65σ/σTS was observed for Equation (9) synonymous with Equations (4) and (5).

The current method has the advantage of extrapolating favourably with respect to stress and temperature when compared to other reviewed creep methods. A finite element model has been formulated based on this creep curve method, the accuracy of which has been verified against stress relaxation data, the results of which indicate good correlation.


The authors gratefully acknowledge the financial and technical support from Rolls-Royce plc and also the financial support from EPSRC.


1. Larson F.R., Miller E.J. Time-Temperature relationship for rupture and creep stresses. Trans. ASME. 1952;74:765–775.
2. Evans R.W., Parker J.D., Wilshire B. An extrapolation procedure for long-term creep strain and creep life prediction with special reference to 0.5Cr0.5Mo0.25V ferritic steels; Proceedings of Recent Advances in Creep and Fracture of Engineering Materials and Structures; Swansea, UK. 24–27 March, 1981; Swansea, UK: Pineridge Press; 1982. pp. 135–184.
3. Evans R.W., Beden I., Wilshire B. Creep life prediction for 0.5Cr0.5Mo0.25V ferritic steel; Proceedings of the 2nd International Conference on Creep and Fracture of Engineering Materials and Structures; Swansea, UK. 1-6 April 1984; Swansea, UK: Pineridge Press; 1984. pp. 1277–1296.
4. McVetty P.G. Creep of metals at elevated temperatures—The hyperbolic sine relation between stress and creep rate. Trans. ASME. 1943;65:761–769.
5. Wilshire B., Battenbough A.J. Creep and creep fracture of polycrystalline copper. Mater. Sci. Eng. A. 2007;443:156–166. doi: 10.1016/j.msea.2006.08.094. [Cross Ref]
6. Wilshire B., Whittaker M.T. The role of grain boundaries in creep strain accumulation. Acta Mater. 2009;57:4115–4124. doi: 10.1016/j.actamat.2009.05.009. [Cross Ref]
7. Wilshire B., Scharning P.J. Long-Term creep life prediction for a high chromium steel. Scr. Mater. 2007;56:701–704. doi: 10.1016/j.scriptamat.2006.12.033. [Cross Ref]
8. Wilshire B., Scharning P.J. Extrapolation of creep life data for 1Cr–0.5Mo steel. Int. J. Press. Vessel. Pip. 2008;85:739–743. doi: 10.1016/j.ijpvp.2008.04.002. [Cross Ref]
9. Evans M. Prediction of long-term creep rupture data for 18Cr–12Ni–Mo steel. Int. J. Press. Vessel. Pip. 2011;88:449–451. doi: 10.1016/j.ijpvp.2011.08.001. [Cross Ref]
10. Wilshire B., Scharning P.J. Prediction of long term creep data for forged 1Cr–1Mo–0.25V steel. Mat. Sci. Tech. 2008;24:1–9. doi: 10.1179/174328407X245779. [Cross Ref]
11. Abdallah Z., Perkins K., Williams S. Advances in the Wilshire extrapolation technique—Full creep curve representation for the aerospace alloy Titanium 834. Mater. Sci. Eng. A. 2012;550:176–182. doi: 10.1016/j.msea.2012.04.054. [Cross Ref]
12. Bryant D.J., McIntosh G. The Manufacture and evaluation of a large turbine disc in cast and wrought Alloy 720Li; Proceedings of Superalloy 1996; Seven Springs, PA, USA. 22–26 September 1996; Warrendale, PA, USA: The Minerals, Metals & Materials Society; 1996. pp. 713–722.
13. BSI (British Standards Institution) Tensile Stress Relaxation Testing, Procedure for Testing Machines. BSI; London, UK: 2003. BS EN 10319-1:2003.
14. Monkman F.C., Grant N.J. An empirical relationship between rupture life and minimum creep rate in creep-rupture tests. Proc. ASTM. 1956;56:593–620.
15. Gopinath K., Gogia A.K., Kamat S.V., Balamuralikrishnan R., Ramamurty U. Tensile properties of Ni-Based Superalloy 720Li: Temperature and strain rate effects. Metal. Mater. Trans. A. 2008;39:2340–2350. doi: 10.1007/s11661-008-9585-3. [Cross Ref]
16. Andrade E.N. On the viscous flow in metals, and allied phenomena. Proc. Roy. Soc. Lond. A. 1910;84:1–12. doi: 10.1098/rspa.1910.0050. [Cross Ref]
17. Bailey R.W. The utilization of creep test data in engineering design. Proc. Inst. Mech. Eng. 1935;131:209–284. doi: 10.1243/PIME_PROC_1935_131_012_02. [Cross Ref]
18. Garofalo F. Fundamentals of Creep Rupture in Metals. Macmillan; London, UK: 1965.
19. Evans R.W. A constitutive model for the high-temperature creep of particle-hardened alloys based on the θ projection method. Proc. Roy. Soc. Lond. A. 2000;456:835–868. doi: 10.1098/rspa.2000.0539. [Cross Ref]
20. Kachanov L.M. Time of the rupture process under creep conditions. Izv. Akad. Nauk. SSR Otd. Tech. Nauk. 1958;8:26–31.
21. Rabotnov Y.N. Creep Problems in Structural Members. North-Holland Publishing Company; North-Holland, the Netherlands: 1969.
22. Othman A.M., Hayhurst D.R. Multi-Axial creep rupture of a model structure using a two parameter material model. Int. J. Mech. Sci. 1990;32:35–48. doi: 10.1016/0020-7403(90)90150-H. [Cross Ref]
23. Sherby O.D. Factors affecting the high temperature strength of polycrystalline solids. Acta Metal. 1962;10:135–147. doi: 10.1016/0001-6160(62)90058-5. [Cross Ref]
24. Wazzan A.R. Lattice and grain boundary self-diffusion in nickel. J. Appl. Phys. 1965;36:3596–3599. doi: 10.1063/1.1703047. [Cross Ref]
25. Whittaker M.T., Wilshire B. Creep and creep fracture of 2.25Cr–1.6W steels (Grade 23) Mater. Sci. Eng. A. 2010;527:4932–4938. doi: 10.1016/j.msea.2010.04.033. [Cross Ref]
26. Whittaker M.T., Wilshire B. Long term creep life prediction for Grade 22 (225Cr–1Mo) steels. Mater. Sci. Tech. 2011;27:642–647. doi: 10.1179/026708310X520529. [Cross Ref]
27. Whittaker M.T., Wilshire B. Advanced procedures for long-term creep data prediction for 2.25 chromium steels. Metal. Mat. Trans. A. 2013;44:136–153. doi: 10.1007/s11661-012-1160-2. [Cross Ref]
28. Manonukul A., Dunne F.P.E. High- and low-cycle fatigue crack initiation using polycrystal plasticity. Proc. Roy. Soc. Lond. A. 2004;460:1881–1903. doi: 10.1098/rspa.2003.1258. [Cross Ref]
29. Goodall I.W., Leckie F.A., Ponter A.R.S., Townley C. The development of high temperature design methods based on reference stresses and bounding theorems. J. Eng. Mater. Tech. Trans. ASME. 1979;101:349–355. doi: 10.1115/1.3443701. [Cross Ref]

Articles from Materials are provided here courtesy of Multidisciplinary Digital Publishing Institute (MDPI)