A Modified Walker Model Dealing with Mean Stress Effect in Fatigue Life Prediction for Aeroengine Disks
Mean stress effect plays an important role in fatigue life prediction, and it is discovered that maximum stress has nonnegligible influence on mean stress effect. Therefore, a modified Walker model is proposed to account for mean stress effect on fatigue life of aeroengine disks, which contains the influence of stress ratio and maximum stress on mean stress effect. Eight sets of fatigue data for standard smooth bars from six kinds of materials commonly used in aeroengine disks as well as two sets of experimental data from simulated specimens of turbine disks were employed to investigate the prediction capability of the proposed model against other candidate mean stress relationships. It is found that Goodman model generates most conservative results, while Morrow model overestimates fatigue life for most cases. SWT model yields similar results to Walker model but with less accuracy. The results of the modified Walker model turn out to be superior to those of any other candidate models for all cases examined, especially for large mean stress ones. Thus, the modified Walker model can be an effective method to predict fatigue lives of aeroengine disks influenced by mean stresses.
In aeroengine, most critical regions of disks are always subject to time varying loads with the presence of mean stresses. Many researchers have found that mean stresses have significant influence on fatigue life [1–5]. It is discovered that tensile mean stresses are usually detrimental to component fatigue while compressive mean stresses are beneficial in terms of fatigue strength. Besides, it is worth noting that mean stresses have significant effect on fatigue behavior in high cycle fatigue (HCF), while they have less influence in low cycle fatigue (LCF). This is because that large amount of plastic deformation will be generated in the LCF process, which will significantly reduce any beneficial or detrimental effect of the mean stresses .
In scientific researches, most fatigue tests are conducted particularly under completely reversed conditions with mean stresses being zero. Therefore, fatigue life predicting methods on the base of these fully reversed fatigue test data should be modified to account for mean stress effect for better accuracy. Since the fatigue process of a component consists of the crack initiation phase and the crack propagation phase, thus many researchers have proposed plenty of models to consider the mean stress effect on the crack initiation life and the crack propagation life separately. The commonly used FCG models considering mean stress effect on fatigue crack growth rate are Priddle model , Collipriest model , McEvily model , Forman model , and so on. These models are used in different situations and turn out to be effective . As to the influence of mean stress on crack initiation life, there are many models proposed, like Goodman model , Gerber model, Morrow model , Manson-Halford model, the Smith-Watson-Topper (SWT) parameter , the Walker equation , the generalized energy damage parameter, and so on . In this paper, we mainly focus on the effect of mean stress on crack initiation life, and the fatigue lives mentioned in rest of the paper are all crack initiation lives.
The prediction accuracy of the methods to assess mean stress effect on crack initiation behavior of various materials under different loading conditions has been investigated by Zhu et al. [12, 13], Correia et al. , Dowling et al. [15, 16], and Burger and Lee . Zhu et al. [12, 13] have proposed the mean stress effect correction in fatigue life predictions based on energy parameters. Correia et al. [18, 19] proposed the generalization of the design fatigue life curve for several fatigue damage parameters including the Walker-like strain damage parameter and others with mean stress effects. Susmel  introduced a mean stress sensitivity index into the modified Wöhler curve method to account for the mean stress effect perpendicular to the critical planes under multiaxial loadings. Niesłony and Böhm  proposed a stress-based approach by employing two curves under alternating stress and stress ratio . Kujawski  made use of analogy with Neuber’s rule and proposed a deviatoric version of the SWT model to consider mean stress effect for relatively large compressive mean stress cases. Kamaya and Kawakubo  investigated the effects of mean stress on fatigue properties of type 316 stainless steel and indicated that mean stress correction was not necessary in component design under load control mode and for constrained ratcheting strain region. Ince  proposed a mean stress fatigue model based on the distortional strain energy to account for the mean stress effect on fatigue life for both positive and compressive mean stress conditions. Besides, researchers have proposed many other equations accounting for mean stress effects, some of which are discussed by Nihei et al.  and Kluger and Lagoda .
2. Mean Stress Models in Fatigue
Even though various models have been proposed to consider the mean stress effect, the Goodman model, Morrow model, SWT model, and Walker model are still considered to be most popular methods while predicting the crack initiation life in practical engineering. Therefore, these four mean stress models will be briefly illustrated as follows.
Goodman relationship employs the ultimate tensile strength, , as an important parameter together with mean stress and stress amplitude to obtain the equivalent stress , as is shown in (1). However, the results from the Goodman method is found to be highly inaccurate [15, 27], especially for large mean stress cases.
The Morrow equation shares the same form with Goodman’s, except for employing the true fracture strength , instead of . In some cases, the true fracture strength may be unavailable. Therefore, an alternate form is proposed by substituting the fitting constant, , the stress intercept at cycle of the curve, to estimate the equivalent completely reversed stress amplitude at various stress ratios . Both of the true fracture form and the stress intercept form of Morrow equation are shown in (2a) and (2b). The assumption of is often accurate for steels, as shown by Landgraf ; however, for aluminums alloys the assumption shows less accuracy, as illustrated by Dowling et al. .
The SWT method employs the maximum stress and the stress amplitude as parameters to calculate the equivalent completely reversed stress amplitude , as is shown in (3). As can be seen, the SWT method has the advantage of simplicity and is not depending on any material constants. Besides, it is found that the SWT model provides good life prediction results in the long fatigue life region but is conservative in the low cycle fatigue life region [30, 31].
Unlike the SWT method, the Walker model supposes that mean stress effect is material dependent. Thus, a fitting parameter , which varies from 0 to 1, is employed in the Walker equation to account for the variance of mean stress effect on different materials, as is shown in (4). The advantage of the Walker model is that it provides an opportunity to fit fatigue test data at various mean stresses all together, with being obtained as part of the fitting process. The higher is, the less sensitive of the material reacts to mean stress effect. On the contrary, a lower means that the material is more sensitive to mean stress effect .
Basquin’s equation  describes the relation of fatigue life and the equivalent fully reversed stress amplitude, which can be expressed aswhere is the fatigue strength coefficient, is the fatigue strength exponent, and is the fatigue life. Usually, Basquin’s equation is employed in conjunction with different mean stress models to estimate fatigue lives only with being replaced.
3. The Proposed Model
From previous investigations, it is known that the Walker method gives better predictions than other models when fatigue data are available to fit the adjustable parameter [15–17, 33]. However, while employing the Walker model to estimate fatigue lives of some commonly used materials, such as FGH4095 powder metallurgy superalloy (600°C) and GH4169 wrought superalloy (650°C), it is noticed that the estimation results of the Walker equation distinctly changed from conservative to overestimated with the increase of stresses for large mean stress cases, which are shown in Figure 1. That is, for large mean stress cases, the Walker equation gives relatively conservative results in the low stress region; conversely, the estimated fatigue life results become nonconservative while the stresses are large, which certainly brings rather great errors.
(a) Walker, FGH4095 600°C
(b) Walker, GH4169 650°C
Therefore, further researches were carried out to investigate the influence of maximum stress on mean stress effect. Based on the Walker model and Basquin’s equation, the fitting parameter for each stress level is determined by the fatigue data on that stress level from different stress ratios through multiple linear regression, as is shown in (6), where and are the fitting results of Walker model for all data sets and will not change with different stress levels.
It can be obviously seen from Figure 2 that the fitting parameter varies with the change of logarithmic of maximum stress almost in a linear relationship for both FGH4095 and GH4169 superalloy. Since acts as an important index of mean stress effect in the Walker model, it can be deduced that maximum stress has nonnegligible influence on mean stress effect as well and should be involved while predicting fatigue life with mean stress effect for better accuracy. Thereupon the modified Walker model is proposed.
The equivalent completely reversed stress amplitude based on the proposed modified Walker model is shown in (7), in which the material dependent fitting parameter is replaced by compared with Walker model, with one additional fitting parameter to account for the influence of maximum stress on mean stress effect. Therefore, the modified Walker model not only contains the influence of material on mean stress effect by employing parameter , just as do in the Walker model, but also reflects the influence of maximum stress through parameter . Taking the logarithm with base of 10 on the combination of (5) and (7), the modified Walker model can be expressed as (8), with being the fitting results of multiple linear regression by all data sets. Obviously, the parameters in (7) and (8) have the following relations: , and .
Though the modified Walker equation contains two fitting parameters while the Walker equation merely includes one, the determination of the two parameters in modified Walker equation is not complicated compared to that of the one parameter in Walker equation. Since (7) shares the same form as (4), there is no difference in the multiple linear regressions of fatigue test data at various mean stresses for the Walker model and modified Walker model. The fitting process of the Walker model has been specifically demonstrated by Dowling et al. in ; therefore it will not be repeated here.
4.1. Comparison of Mean Stress Models with Experimental Data from Standard Smooth Bars
Fatigue test data of standard smooth bars for various stress ratios are employed to verify the prediction capability of the candidate mean stress models in (1) to (4) and (7) by plotting the equivalent fully reversed stress amplitude versus fatigue life . In what follows, of each test data point is calculated by the corresponding mean stress equations mentioned previously.
The plots of versus for FGH4095 (600°C) and GH4169 (650°C) superalloy are shown in Figures 3 and 4 separately, which include the results of Goodman model, Morrow model, SWT model, and the modified Walker model, respectively (the results of Walker model have been plotted in Figure 1 already). In each plot, the degree of resulting data points which consolidate with the fitted line gives an indication of the success of the mean stress relationship. It should be noted that the fitted lines in the Goodman model, Morrow model, and SWT model are obtained just from the fatigue test data of (or ), while the fitted lines of the Walker model and modified Walker model result from fitting the data sets at all mean stresses.
(a) Goodman, FGH4095 600°C
(b) Morrow, FGH4095 600°C
(c) SWT, FGH4095 600°C
(d) Modified Walker, FGH4095 600°C
(a) Goodman, GH4169 650°C
(b) Morrow, GH4169 650°C
(c) SWT, GH4169 650°C
(d) Modified Walker, GH4169 650°C
For FGH4095 superalloy (600°C), it can be recognized in Figure 3(a) that the data of stress ratio distribute far above the fitted line of data in the Goodman relation plot, and the data set of stress ratio also distribute above the fitted line, but they lie comparatively closer than the data set, which indicates that the Goodman relationship is applicable to estimate fatigue life for small mean stress cases; however when the mean stress gets larger, Goodman relation becomes highly inaccurate, and the larger the mean stress is, the less the accuracy may become. The results of Morrow equation for and data sets distribute below the zero mean stress line. Similarly, the results of data set lie much closer than that of data set. In Figure 3(c), reasonable results are obtained by the SWT model, which show better accuracy than the Goodman and Morrow relations. Furthermore, it is discovered that the Walker model provides slightly better results than the SWT model through comparison of Figures 1(a) and 3(c), which can be attributed to the adjustable fitting parameter by varying the resulting curve to fit data in cases where the SWT method is not very accurate . Meanwhile, it is observed that either SWT or Walker model provides conservative results in the low stress region for large mean stress cases, while in the high stress region the opposite happens. Compared with other plots in Figures 1 and 3, it can be clearly observed that the modified Walker model with the consideration of mean stress effect on fatigue life influenced by both materials and maximum stresses gives superior results, which is illustrated by the fatigue test data of all stress ratios in Figure 3(d) which basically equally distribute on both sides of the fitted line.
Similar conclusions can be drawn for GH4169 superalloy (650°C) in Figure 4; only the results of Goodman model distributed much far from the fitted line than Morrow model for data sets of and .
To quantitatively summarize the accuracy of the mean stress equations for different materials, the difference between the experimental logarithmic life and predicted logarithmic life for fatigue data sets is taken to calculate the prediction error of each model, which can be expressed aswhere is the experimental fatigue life and is the predicted fatigue life of a given data point. Obviously, a positive error means that the prediction is conservative, while a negative error indicates the prediction is overestimated.
Besides, a probabilistic method is applied to analyze the prediction errors in the form of fitted probability density functions (PDF). The mean and standard deviation values are effective parameters to illustrate the prediction accuracy of different models. The fitted PDFs for prediction errors of each mean stress relationship are shown in Figures 5 and 6 for FGH4095 (600°C) and GH4169 (650°C) superalloy, separately. As clearly seen in Figure 5, the Goodman model tends to provide conservative predictions with positive mean values of prediction error for FGH4095 superalloy (600°C). On the contrary, the Morrow model tends to overestimate fatigue lives with negative errors. Besides, the standard deviations of Morrow model are comparatively large, which indicates that the prediction accuracy of Morrow model shows less stability. The results of SWT model and Walker model share the same tendency, while the Walker model provides slightly better results with larger probability to locate near the zero-error region. It can be seen that the modified Walker model shows best performance with mean values around zero and small standard deviations compared with other mean stress models. Particularly for the large mean stress cases like data set, the modified Walker model provides much better results than other models.
(a) , FGH4095 600°C
(b) , FGH4095 600°C
(a) , GH4169 650°C
(b) , GH4169 650°C
The prediction accuracy of each mean stress model for GH4169 superalloy is similar to that of FGH4095 superalloy. Only the Goodman model shows less accuracy while Morrow model brings better results compared with FGH95 superalloy. Similarly, the modified Walker model provides better results than SWT and Walker models, even though these two models show much better results than the Goodman and Morrow models.
Besides the FGH4095 (600°C) and GH4169 (650°C) superalloy, some other materials which are commonly used in aeroengines such as FGH4095 (500°C), GH4169 (500°C), and FGH4096 superalloy , TC11 titanium alloy and 35Cr2Ni4MoA high-strength steel , and 2014-T6 aluminum  are also employed to demonstrate the accuracy of the mean stress models. The material properties and the fitting constants of Basquin’s, Walker, and modified Walker equations are shown in Table 1. It should be noted that the parameters in Table 1 are determined by the analysis of fatigue data in corresponding references.
The plots of versus by the newly proposed modified Walker model for other materials examined in this study are shown in Figure 7. Besides, the mean and standard deviation values of prediction errors for all the materials examined are shown in Table 2, in which stands for mean value while is standard deviation value of the prediction errors. It should be noted that all the and in Table 2 are obtained by all stress ratio data sets.
(a) FGH4095 500°C
(b) GH4169 500°C
(c) FGH4096 550°C
(d) TC11 300°C
(e) 35Cr2Ni4MoA 160°C
(f) 2014-T6 Al
By analyzing plots in Figure 7 as well as the mean and standard deviation values in Table 2 for all the materials examined, it can be seen that the results of candidate mean stress models for the materials examined share the same tendency with FGH4095 and GH4169 superalloy. That is, the Goodman relation gives the most conservative results with positive mean values of error, while Morrow model overestimates the fatigue life for most cases. Besides, the standard deviation values of both Goodman and Morrow are comparatively large, which indicates that the prediction accuracy of these two models varies strongly. SWT model, as a special case of Walker equation, fails to reflect the influence of different materials on mean stress effect and thus provides less prediction accuracy than the Walker equation. The modified Walker model brings the least standard deviation values as well as mean values of error being zero, showing favorable capability to estimate fatigue lives for different materials and different mean stress cases. As to the zero values for in the Walker and modified Walker models, because both the Walker and modified Walker model determine parameters through multiple linear regression of all stress ratio data sets, which is essentially based on the least squares principle, there is no wander that the mean values of prediction errors for the Walker and modified Walker model become zero.
4.2. Comparison of Mean Stress Models with Experimental Data from Simulated Specimens
The previous work has proved that the modified Walker model provides favorable results for standard smooth bars of all materials examined. However, it remains to be testified whether the modified Walker model is applicable to estimate fatigue lives of real components in aeroengine. Therefore, two sets of experimental data from a type of simulated bolt-hole specimen and a type of simulated run-way hole specimen of turbine disk at 550°C are analyzed to investigate the prediction capability of candidate models. The main geometric parameters of the two specimens are shown in Table 3 and the finite element (FE) models of the two specimens are shown in Figure 8.
(a) Bolt-hole specimen
(b) Run-way hole specimen
The finite element analysis (FEA) results of the two specimens are shown in Figure 9. It should be noted that Figure 9 shows the contour plot of modified Walker equivalent stress for the two simulated specimens. Since both of the critical regions in the two specimens are obviously notched, the theory of critical distances (TCD) should be employed to properly accounting for the detrimental effect of stress concentration, and the line method (LM) of TCD is found to be an effective method :where and are the polar coordinates of the stress point from the critical point, with the stress gradient being the axis, and is the distance from the notch tip to the minimum of relative stress gradient , as is shown in The predicted fatigue life of the two simulated specimens can be obtained by modified Walker model in conjunction with LM, as is shown in Table 4. Besides, the predicted lives by other models and experimental fatigue lives of the two specimens are also provided for comparison. Obviously, the modified Walker model provides superior results than any other models for both of the two types of simulated specimens. Although the Walker model also gives reasonable fatigue life estimation, it still shows less accuracy for failing to reflect the effect of maximum stress on mean stress effect. Morrow model overestimates the fatigue lives while SWT model generates conservative results for both of the two specimens.
To comprehensively account for the influence of stress ratio, material, and maximum stress on mean stress effect, the modified Walker model is proposed. Then, eight sets of experimental data from standard smooth bars of six different materials and two sets of fatigue test data from simulated bolt-hole and runway-hole specimens of a turbine disk have been employed to investigate the prediction capability of the proposed model against other mean stress relationships commonly used. The results of the modified Walker model turn out to be superior to other candidate models for all materials examined. Goodman model is found to be most inaccurate, and Morrow equation shows better results than Goodman model. However, the results of the two models are both far from satisfactory. The SWT model brings reasonable results; what is more, it has the advantage of brevity; therefore, the SWT method is a good choice for general use. The Walker model provides satisfactory performance when mean stresses are relatively small by employing a material dependent parameter ; however, it fails to correlate the fatigue data well for the large mean stress cases.
The proposed modified Walker model is based on the Walker model with an extra item of stress and stress ratio. Thus the modified Walker model inherits the advantage of Walker model and contains the influence of maximum stress on mean stress effect as well, showing favorable capability to estimate fatigue lives for different materials and different mean stress cases, sharing the same data pool and procedure to determine the fitting parameters, without increasing the calculation cost. Even though being a little complicated than the Walker model, the modified Walker model gives more precise results than the Walker model for all data sets examined, especially in large mean stress cases, showing better accuracy and larger scope of application.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The authors gratefully acknowledge the support from the China Gas Turbine Establishment.
J. Goodman, Mechanics Applied to Engineering. Longmans, Longmans, Green and Co., London, UK, 1919.
J. Morrow, “Fatigue properties of metals, section 3.2,” in Fatigue Design Handbook, Pub. No. AE-4, Society of Automotive Engineers, Warrendale, Pa, USA, 1968, Section 3.2 is a summary of a paper presented at Division 4 of the SAE Iron and Steel Technical Committee, Nov. 4, 1968.View at: Google Scholar
K. N. Smith, P. Watson, and T. H. Topper, “Stress- strain function for the fatigue of metals,” Journal of Materials, vol. 5, no. 4, pp. 767–778, 1970.View at: Google Scholar
K. Walker, “The effect of stress ratio during crack propagation and fatigue for 2024-T3 and 7075-T6 aluminum,” in Effects of Environment and Complex Load History on Fatigue Life, ASTM STP 462. American Society for Testing and Materials, pp. 1–14, Philadelphia, Pa, USA, 1970.View at: Google Scholar
N. E. Dowling, “Mean stress effects in stress-life and strain life fatigue,” in Proceedings of Fatigue 2004: Second SAE Brazil International Conference on Fatigue, Sao Paulo, Brazil, June 2004, Also SAE Paper No. 2004-01-2227, SAE International, Warrendale, PA, 2004.View at: Google Scholar
Y. L. Lee, M. E. Barkey, and H. T. Kang, “Metal fatigue analysis handbook: practical problem-solving techniques for computer-aided engineering,” pp. 151–157, Butterworth-Heinemann Press, UK, 2011.View at: Google Scholar
J. E. Collipriest, R. M. Ehret, and C. Thatcher, “Fracture mechanics equations for cyclic growth,” Tech. Rep. MFS-24447, NASA Technology Utilization Report, 1973.View at: Google Scholar
J. A. F. O. Correia, S. Blasón, A. M. P. De Jesus, A. F. Canteli, P. M. G. P. Moreira, and P. J. Tavares, “Fatigue life prediction based on an equivalent initial flaw size approach and a new normalized fatigue crack growth model,” Engineering Failure Analysis, vol. 69, pp. 15–28, 2016.View at: Publisher Site | Google Scholar
R. W. Landgraf, “The resistance of metals to cyclic deformation,” in Achievement of High Fatigue Resistance in Metals and Alloys, ASTM STP 467. American Society for Testing and Materials, pp. 3–36, Philadelphia, Pa, USA, 1970.View at: Google Scholar
O. H Basquin, “The exponential law of endurance tests,” American Society Testing Mater Proceedings, vol. 10, pp. 625–630, 1910.View at: Google Scholar
The High Temperature Material Branch of Chinese Society for Metals. China Superalloys Handbook, China Standard Press, Beijing, China, 2012.
Editorial Committee for the Material Data Handbook for Areo-Engine Design, Material Data Handbook for Aero-engine Design, Aviation Industry Press, Beijing, China, 2010.
B. J. Lazan and A. A. Blatherwick, Fatigue properties of aluminum alloys at various stress ratios,” Part 1, rolled alloys. WADC TR 52-307, Part 1, Wright Air Development Center. Wright-Patterson AFB, Dayton, Ohio, USA, 1952.