1. Introduction
It is widely acknowledged that thermo-mechanical fatigue (TMF) is one of the most significant challenges facing the gas turbine sector as designers strive to meet the demanding emissions and fuel burn ACARE (Advisory Council for Aeronautics Research and Innovation in Europe) “Flightpath 2050” targets. Increasing operating temperatures to improve cycle efficiency, in conjunction with weight reduction strategies that include thinner disc rims, have led not only to more aggressive TMF cycles, but also to the requirement for assessment of the behaviour of components for which TMF was not previously considered significant.
TMF testing however is difficult, expensive and time consuming and there is a wide range of variables that require consideration. Such variables include peak cycle temperature, temperature cycle range, applied strain range, R ratio, phase angle (between thermal and mechanical strains), waveform and applied strain rate. Accordingly, a data generation programme that encompasses all of these factors is expensive and unrealistic. Therefore there is a clear requirement for a robust predictive model which can provide accurate extrapolation from a known set of results. Furthermore, it is desirable that the model, if possible, utilise isothermal data in order to further simplify the input data requirements. However, this is obviously a significant extension of models which predict only isothermal data and it is necessary to investigate whether such models are compromised by the damage mechanisms involved in TMF.
This paper focuses on modelling the effects of phase angle on TMF life in the polycrystalline nickel alloy RR1000. RR1000 is commonly produced through a powder metallurgy process, which is widely anticipated to play a significant role in upcoming engines as a turbine disc material. However studies into the TMF behaviour of the alloy have been extremely limited and TMF lifing of many alloys tends to follow the approach of attempting extrapolation from isothermal data and particularly assuming isothermal data at the peak cycle temperature to be a worst case scenario. Within the gas turbine, however, this approach may be extremely limited due to the wide range of phase angles possible and the marked effect of a TMF cycle.
The phase angle as defined for a TMF test is the shift which occurs between the mechanical strain and the temperature applied to the material. If the mechanical strains and temperature increase/decrease in proportion with each other, the cycle is considered to be In-Phase (IP). If one always increases as the other decreases, the cycle is Out-Of-Phase (OOP) or more accurately, −180° Out-Of-Phase. Clearly there are also an infinite number of potential phase angles which could be applied between these two extremes and furthermore an added complication is that the cycle can take one of two directions around the strain-temperature loop, either clockwise (CW) or anticlockwise (ACW), as shown in
Figure 1.
Laboratory simulations of TMF testing are often, by necessity, idealised versions of in service cycles. It is impossible to undertake a TMF testing programme which will encompass all of the possible variations which may occur in service. However, the range of phase angles that may occur in component TMF cycles will be covered by tests performed for a number of standard phase angles together with modelling work that allows the results to be interpolated between phase angles and, ideally, to different maximum and minimum cycle temperatures. This will give designers the ability to predict lives for all component loading cycles. As such it is critical that any model which is sought to describe TMF behaviour is capable of interpolation across phase angles, and if at all possible different maximum/minimum temperatures, in order to allow designers the capacity for accurate life prediction of critical parts.
Figure 1.
Phase angles (ϕ) and directions (clockwise (CW) or anticlockwise (ACW)).
Figure 1.
Phase angles (ϕ) and directions (clockwise (CW) or anticlockwise (ACW)).
Previous attempts to model TMF behaviour in the literature [
1,
2,
3,
4,
5], especially in nickel based superalloys for disc applications, are quite limited. They do however, fall broadly into two categories; the first assumes that under strain control loading a significant fraction of the TMF life is spent in crack initiation and as such traditional approaches such as strain–energy partitioning, Manson-Coffin [
6]
etc., can be applied [
7,
8,
9,
10]. Conversely, the second approach assumes that cracks form very early in the test, possibly even on the first cycle [
11,
12,
13]. On this basis the fatigue life of the specimen will be dominated by crack propagation. Each type of approach has been shown to have its own benefits and drawbacks and it is also clear that the most appropriate method will be material and temperature dependent.
4. Discussion
Based on the wide range of variables which occur during a TMF test it is clear that a lifing approach capable of holistic data fitting may not be possible. However, an objective of the current work is to assess possible lifing strategies in order to provide a basis for a more widespread approach. As such, the discussion will focus on the advantages and drawbacks of methods based around mean stress, stabilised stress range and crack propagation life. Based on the results produced within this work, no difference was seen between specimens loaded in a clockwise (CW) or anticlockwise (ACW) direction around the strain–temperature loop, for −90° OOP tests. As such the approaches here are simplified by not considering directionality in these tests.
One relatively straightforward approach to ordering TMF results can be made by simply considering the stabilised mean stress achieved during the test, with the expectation being that an increased mean stress for a constant stress range will lead to a reduction in fatigue life. Due to the non-isothermal nature of the TMF test, stress–strain hysteresis loops can vary dramatically in shape for different phase angles, leading to significant changes in mean stress. The most common cause of this is the alternating contributions between creep strain and plasticity, typical in the TMF cycle. By considering an IP cycle, in which the strain and temperature increase in proportion, it can be appreciated that yield will occur within the material and dislocations will be generated within the grains, resulting in plastic deformation. However, at higher temperatures creep deformation becomes more significant in the alloy, which acts to lower the applied stress through stress relaxation, since the test is under strain control. Under an IP TMF cycle, with creep therefore dominating at the peak strain, this stress relaxation acts to decrease the mean stress throughout the test. Conversely in a −180° OOP test, creep is most prevalent at the minimum strain, and therefore acts to increase the mean stress, as can be observed in
Figure 4.
Clearly for a series of tests performed at different strain ranges, this approach is overly simplistic in the presence of varying peak stresses and stress ranges. For example a −180° OOP test with a high strain range would not be expected to show a longer fatigue life than an IP test with a low strain range, simply because it has a lower value of mean stress. However, for the current evaluation of the effect of phase angle on fatigue life, the method can be instructive for a series of experiments where a constant strain range is utilised.
Following initial plasticity on loading, as demonstrated in
Figure 4, similar stress ranges occur across the range of phase angles, enabling the effect of mean stress to be more clearly evaluated.
Figure 6 shows a comparison of mean stress against normalised cycles to failure for the range of tests performed at Δε
mech = 1%. The data is normalised by the number of fatigue cycles to failure for an isothermal test, performed at the peak temperature (700 °C) of the TMF cycle and clearly demonstrates the effects of mean stress with values <1 showing a detriment to life, and values >1 showing an extension of life. With an
R2 value of 0.81 it is clear that a relationship exists and that the TMF life is strongly dependent on mean stress. However, it should again be emphasised that this relationship exists only for specimens with a comparable stabilised stress range, brought about in these experiments by the application of a constant applied strain range of 1% throughout the tests. In the case of extrapolation to a varying applied strain range (and hence stabilised stress range) the relationship holds little value because of the tendency of stabilised stress range to become the dominant factor.
Figure 6.
Mean stress vs. cycles to failure for varying phase angle tests at Δε of 1%.
Figure 6.
Mean stress vs. cycles to failure for varying phase angle tests at Δε of 1%.
Typically, more holistic fatigue lifing approaches will be either strain range [
6,
7,
16] or stress range [
17] based. In the current work a strong correlation between stabilised stress range and TMF life in the IP and −180° OOP experiments indicates that a Basquin [
17] type approach could potentially form the basis of a model which is capable of TMF life prediction and extrapolation across a range of phase angles.
As mentioned previously however, TMF experiments are costly, difficult and time consuming and it is preferable that predictions can be made based on isothermal data. However, as indicated in the literature [
18], simple “worst case scenario” and generalised fatigue life reduction factors fail to capture the complexities of the effect of phase angle, potentially leading to overly conservative life predictions or even premature component failures. In the present work, efforts have been made to effect a creep and strain controlled fatigue prediction based on TMF hysteresis loops which have been derived from a database of isothermal RR1000 test results. SC03 is a non-linear finite element analysis program which utilises the Mroz multilayer plasticity model to predict cyclic deformation based on stress–strain data obtained from isothermal testing [
19]. The wide range of input data used to produce the deformation model for RR1000 provides confidence in its ability to predict conditions outside of the envelope of this original data. It also allows for an iterative approach over both temperature and strain conditions by which the replication of TMF stress–strain hysteresis loops can be achieved.
Examples of loops produced by SC03 are shown in
Figure 7 and compared with experimental data generated under TMF loading conditions. It can be seen that peak stress, minimum stress and loop shape are all captured well by the model, providing confidence in its ability to extrapolate to alternative phase angles and produce a model which can provide TMF predictions over a wide range of phase angles when used in association with predicted stresses from SC03 to estimate TMF test failure lives. As described earlier a strong relationship exists between the stabilised stress range and the fatigue life, and as such Basquin’s law [
17] has been applied to derive a predictive equation for the experiments.
Figure 7.
Modelling of TMF loops (a) IP Δε of 0.8%; (b) −180° OOP Δε of 1.2%.
Figure 7.
Modelling of TMF loops (a) IP Δε of 0.8%; (b) −180° OOP Δε of 1.2%.
In order to prove the effectiveness of the technique, the predictive equation has been derived based purely on IP experiments. Stress–strain hysteresis loops were then produced for each of the remaining tested conditions (4 × −180° OOP, −45° CW, −90° CW, −90° ACW, −135° ACW) and predictions made of TMF life using the Basquin approach [
17]. The equation takes the form
Where ∆σ is the stress range and Nf is the number of TMF cycles to failure again normalised by the fatigue life of an isothermal test at ∆ε = 1% at peak cycle temperature (700 °C), for consistency. The constants were optimised at σ′f = 989.5, b = −0.0818.
Figure 8 shows the predictions of phase angles between 0° and −180° using the above equation. It is clear that reasonable predictions are obtained for all of the data. If anything the method has the tendency to slightly underpredict the results. This may be a result of differences in the crack propagation life of the IP and −180° OOP tests brought about by the difference in temperature at peak stress. Furthermore it could be argued that when considering the range of phase angle tests performed at ∆ε
mech = 1% the method fails to capture differences between the phase angles, and certainly does not separate the results in the manner in which the mean stress approach did. It should be noted however, that only a factor of 2 difference is observed in the failure life of these experiments, and for this reason it is difficult to draw strong conclusions. Furthermore, it is also necessary to consider that due to a combination of limited creep strain accumulation at the peak cycle temperature, the accumulation of significant mean stresses does not occur in the current tests. Future approaches for higher peak cycle temperatures may be well advised to consider either an adaptation of the Basquin approach to also consider mean stress, or the application of an approach based on both stress range and mean stress.
Figure 8.
Basquin fit predictions of TMF behavior.
Figure 8.
Basquin fit predictions of TMF behavior.
Both of the approaches detailed so far consider only a total life philosophy. Previous work [
20] has indicated that the ordering of IP, OOP and IF (Isothermal Fatigue) tests may prove difficult for short TMF life experiments when using only methods such as Basquin [
17] or Manson-Coffin [
6,
7]. Indeed the ordering of these test results may prove counterintuitive to the ordering seen for longer life experiments. As such it is useful to further investigate the length of the crack propagation phase within the current series of tests, so that the ratio of crack initiation/propagation can be considered.
Although a facility for investigating TMF crack propagation rates is under development at Swansea University, information regarding these rates is extremely limited in the published literature [
21,
22,
23,
24]. A fractographic study was undertaken of the tested specimens and observations have revealed extensive striations in the −180° OOP specimens,
Figure 9. Past observations in the literature have indicated that striation counting in nickel based superalloys is a reliable means of estimation of crack propagation rates [
25,
26]. As such it was possible to identify a crack initiation site and count striations until failure occurred in the specimen by ductile overload.
The number of striations were recorded for two −180° OOP specimens (∆ε
mech = 1.2% and 1.4%) and one −90° OOP ACW specimen (∆ε
mech = 1%). The fractured IP testpieces were examined for striations, but determination was found to be difficult due to the clear interactions of other high temperature damage mechanisms (
i.e., creep and environmental damage) at peak stress in these experiments. It should be noted that a comprehensive review of damage mechanisms occurring under TMF loading in nickel based superalloys is available in the open literature [
27]. However, in order to simplify the lifing approaches applied here, time dependent damage mechanisms (
i.e., creep and environmental damage) are considered as a single entity.
Observations show that significant fractions of these experiments are spent in crack propagation, since the striation count approaches the total specimen life. Indeed all three specimens showed striations that accounted for >80% of the total fatigue life. Furthermore if the striation counts are believed to be accurate, the initiation phase diminishes with increasing strain until, allowing for minor counting errors, the material cracks on the first cycle in the ∆εmech = 1.4%, −180° OOP test, and spends the remainder of its life in propagation.
Figure 9.
Fracture surface of RR1000 from an OOP180° anti-clockwise test (TMF013), with ∆εmech = 1.4% resulting in a life of Nf = 312.
Figure 9.
Fracture surface of RR1000 from an OOP180° anti-clockwise test (TMF013), with ∆εmech = 1.4% resulting in a life of Nf = 312.
Whilst the previous approach based on Basquin’s method shows a reasonable fit, it is unlikely that a single fit would prove adequate if more results with TMF lives >1000 cycles were generated. With crack initiation likely to dominate life, the crack propagation phase would become negligible in the determination of a predictive equation, and hence the prediction would most likely underpredict some of the short term test data observed in the current work. As such, in order to produce a total life prediction methodology, it is clear that there is a requirement for an approach which can also account for crack propagation, rather than assuming total life to be a result of a long crack initiation period followed by rapid crack growth to failure. An existing database of isothermal data in RR1000 can again be accessed for this purpose.
In order to derive a prediction for the propagation rates of TMF experiments a model has been produced based on the aforementioned isothermal data. Provided the strain controlled TMF cycle is accurately described, in terms of time, temperature, peak stabilised stress, strain range and phase angle, these simple variables can be used as inputs for the model. Furthermore the aspect ratio of the crack is required in order to calculate the geometry factor,
Y, such that
(Δ
K is stress intensity range, Δσ is stabilised stress range,
a is crack length). Based on typical crack shapes observed in the current specimens, the variation of Δ
K was estimated with increasing crack depth, a. In order to achieve this, a simplified approach was utilised where the hollow cylindrical test piece was approximated by a rectangular geometry. Analysis of the relationship between crack size and test piece geometry meant that this approximation was deemed reasonable. Due to the relatively long gauge length of the specimen, no significant effects of boundary conditions influenced the calculations. A stress intensity profile for a crack growing in the current hollow cylindrical specimens can be seen in
Figure 10.
Time dependent and time independent crack growth rates are then calculated, so that the contribution from the two damage mechanisms at high temperature can be understood. The time dependent crack growth term, based on an Arrhenius type expression, is of the form:
Where da/dt is the time dependent rate of crack growth, A*and n are material fitting constants, Q is the activation energy for time dependent crack growth (265 kJ/mol−1), R is the gas constant, T is temperature and K is the stress intensity factor.
The time dependent crack growth rate is integrated through the whole loading cycle and therefore implicitly includes any minor cycle effects. Time dependent crack growth will occur when both high temperatures and high tensile stresses are present, such as in an IP test.
Figure 10.
Stress intensity (
K)
vs. crack length (
a) for a semi-elliptical crack as shown in
Figure 9.
Figure 10.
Stress intensity (
K)
vs. crack length (
a) for a semi-elliptical crack as shown in
Figure 9.
The model identifies minor cycles based on stress turning points and accounts for them in the time independent crack growth rate calculation through an exchange rate term. The exchange rate varies slightly at different values of major cycle ΔK, but in practice the exchange rate at a Kmax of 20 MPa
is used throughout the life calculation.
If the peak stress in the loading cycle occurs at a temperature which is lower than the cycle maximum, the time independent crack growth is calculated by averaging the rates from the temperature associated with the peak cycle stress and the peak cycle temperature.
In the −180° OOP tests considered in the current fractography study it is clear that time independent contributions dominate fatigue crack propagation life, since the maximum temperature occurs at the minimum stress. Subsequently a relatively slow crack growth rate is obtained. The OOP experiments performed at −90° CW and ACW give very small time dependent contributions since the peak stress in the cycle occurs at the relatively low temperature of 500 °C and the stresses are low close to the peak cycle temperature. The IP tests performed at 0° give large contributions from both time dependent and time independent crack growth, since the maximum stress occurs at the maximum temperature. As a consequence the crack propagation life is decreased due to interactions between fatigue and other damage mechanisms such as environmental degradation and creep.
The results of the model can be compared with the experimental data,
Figure 11.
Figure 11a shows a reasonable comparison between predictions and experimental data for a −180° OOP test performed at ∆ε = 1.2%, although the crack growth rates in the specimen do exceed the model predictions slightly. However, within the bounds of typical experimental scatter in crack propagation testing, this type of variation is not untypical. Further confidence from the model is gained from the predictions made of a −90° OOP test at ∆ε = 1%,
Figure 11b, since a larger element of time dependence of the crack growth rate must be considered, with maximum stress occurring at 500 °C. Clearly, further evaluation of the model is necessary, and since striation counting may be limited in scope in the current specimens, particularly IP specimens where creep/environment interactions make striation recognition difficult, this comparison will be required to be made with data from the TMF crack propagation facility currently under development.
Figure 11.
TMF crack propagation predictions for (a) −90° OOP test at Δε = 1%; (b) −180° OOP test at Δε = 1.2%.
Figure 11.
TMF crack propagation predictions for (a) −90° OOP test at Δε = 1%; (b) −180° OOP test at Δε = 1.2%.
Each of the methods detailed above clearly have their own advantages and disadvantages. Fatigue life predictions based purely on mean stress provide reasonable ordering of results for a single strain range, but this approach is clearly limited for different strain ranges as the applied stress range varies. However, the trends indicate the importance of the interactions between creep and plasticity within the TMF tests, in particular the shape of the hysteresis loops produced and the maximum and minimum applied stresses. Subsequently the second approach involving the Basquin equation shows the dependence of TMF life on the stabilised stress range where a strong relationship is evident.
The Basquin equation derived for the current material clearly gives a reasonable methodology for prediction of alternative phase angles, based solely on the IP data set. However, it should be acknowledged that this relationship exists only over a limited range of TMF lives, with only two tests showing lives in excess of 3000 cycles. As shown by the striation counting exercise, these experiments are heavily influenced by crack propagation, with even the −90° OOP test at ∆ε = 1% spending greater than 80% of its life in propagation. As peak strains and/or strain ranges are reduced however, this trend is unlikely to continue, with significant fractions of the total life spent initiating a crack, and as such the relationship is unlikely to remain as strong.
In order to overcome this issue it is necessary to consider a combination of methods. The generation of longer life TMF data would undoubtedly result in a shallower fatigue curve at extended lives, which when extrapolated to shorter lives would tend to underestimate the fatigue life response. However, proportionally the propagation life which would be insignificant at longer lives would become increasingly more important at shorter lives. As such a Basquin type approach with an added crack propagation phase is likely to produce the most appropriate methodology for lifing a wide variety of TMF experiments, particularly in the case of notched specimens where crack initiation is usually followed by a significant period of crack propagation.