A simple fatigue life prediction algorithm using the modified NASGRO equation is proposed in this paper. The NASGRO equation is modified by introducing the concept of intrinsic effective threshold stress intensity factor (SIF) range ΔKeffth. One advantage of the proposed method is that the complex growth behavior analysis of small cracks can be avoided, and then the fatigue life can be calculated by directly integrating the crack growth model from the initial defect size to the critical crack size. The fatigue limit and the intrinsic effective threshold SIF range ΔKeffth are used to calculate the initial defect size or initial flaw size. The value of ΔKeffth is determined by extrapolating the crack propagation rate curves. Instead of using the fatigue limit determined by the fatigue strength at the specific fatigue life, the fatigue limit is selected based on the horizontal tendency of the S-N curve. The calculated fatigue lives are compared to the experimental data of two different alloys. The predicted S-N curves agree with the test data well. Besides, the prediction results are compared with that calculated using the FASTRAN code. Results indicate that the proposed life prediction algorithm is simple and efficient.
1. Introduction
Despite more than one hundred years of research, the fatigue life prediction is still a challenging problem. The traditional cumulative fatigue damage (CFD) theories are often used for predicting the fatigue life of structural materials. Different CFD formulas are proposed to predict the S-N curve, such as the famous Basquin form [1] in the stress-based method and the Coffin-Manson relationship [2, 3] in the strain-based method. It should be noted that many factors can affect the fatigue life, such as the material geometry [4, 5], the stress state [6–9], and the effect of environment [10–12]. But these factors are not taken into account in the CFD-based life prediction method. Besides, there are no clear definitions for the fatigue failure in above methods. In a sense, the CFD-based life prediction models are more experience-based rather than physical.
With the advent and technical developments of Linear Elastic Fracture Mechanics (LEFM) theories and computer science, large efforts have been devoted to predict or to understand the crack growth behavior. The Irwin stress intensity factor (SIF) range was first used by Paris and Erdogan [13] to characterize the crack growth behavior under external cycle loading, and the crack propagation rate is characterized using the SIF range concept. Since then, several crack propagation models were proposed, such as the Priddle equation [14] and the Forman model [15], and the crack growth behavior was divided into three stage: the near threshold region, the Paris region, and the unstable crack propagation region. The discovery of crack closure by Elber [16] was another major step in the crack growth analysis. The crack closure is a physical phenomenon that the crack is already closed before the external tensile loading is completely unloaded, and also the crack will still remain closed until the external tensile loading reaches a certain level. It is reported that many mechanism models can explain the crack closure phenomenon, such as the plasticity-induced crack closure [16], the oxide-induced crack closure [17], and the roughness-induced crack closure [18]. Among them, Newman’s finite element model is the best and deserves future developments [19]. Newman’s strip yield model [20, 21] was based on the Dugdale model [22] and was then further developed and implemented into the NASGRO software [23, 24] and the FASTRAN code [25]. Even when these advanced crack growth models which are based on the crack closure are proposed, it is still a tough task today to calculate the fatigue life by directly integrating the crack propagation model from the initial defect size to the critical crack length.
Newman Jr. et al. [26–28] calculated the fatigue lives of various materials and specimen configurations using the FASTRAN code and small-crack theory. The life prediction model imbedded in the FASTRAN code is Newman’s crack closure model. It is known that the growth behavior of small cracks is far different from that of large cracks [29, 30], and the small-crack growth mechanism has not been fully understood yet. Besides, the traditional SIF range ΔK solution may not be applicable to a small crack if the crack initiates from a subsurface inclusion [31]. The other problem is that the value of ΔKeffth and initial defect size which are the key parameters in the LEMF-based life prediction method are selected via trial and error. This paper is aimed at simplifying the life prediction algorithm by using the modified NASGRO equation and the equivalent initial flaw size (EIFS) concept [32]. The fatigue lives can be calculated by directly integrating the crack growth model from the initial defect size to the critical crack size. The complex growth mechanism of a small crack can be ignored in the proposed method. In addition, the value of ΔKeffth is determined by extrapolating the da/dN~ΔKeff relation and the EIFS concept is used to calculate the initial defect size.
2. Fatigue Life Prediction Algorithm
The plasticity-induced crack closure was discovered by Elber [16], and then many numerical models were proposed to simulate the crack closure phenomenon. Newman Jr. firstly developed a ligament or strip yield model [20, 21] based on the Dugdale model framework, and the model was further developed and implemented into the NASGRO code [23, 24] which is widely used for the crack growth analysis; the equation is expressed as(1)dadN=C1ΔKeffC21-ΔKth/ΔKp1-Kmax/Kcq,where C1, C2, p, and q are materials constants. ΔKth is the threshold SIF range, ΔK is the SIF range, Kmax is the maximum SIF, Kc is the critical SIF, and ΔKeff is the effective SIF range and can be expressed as a function of ΔK and U:(2)ΔKeff=UΔK,where U is a function of the stress ratio R and Newman’s crack opening function f and is expressed as(3)U=1-f1-R.
For constant loading, f is expressed as [21](4)f=maxR,A0+A1R+A2R2+A3R3R≥0A0+A1R-2≤R<0,where (5)A0=0.825-0.34α+0.05α2cosπ/2Smaxσ01/α,A1=0.415-0.071αSmaxσ0,A2=1-A0-A1-A3,A3=2A0+A1-1.Smax/σ0 and α are used to calculate the crack open function f in NASGRO software. According to the NASGRO user manual, Smax/σ0 is selected as 0.3 in general [23, 24]. α is the constraint factor which is used to consider the stress state effect. Theoretically, α is 1 in the plane stress condition and 3 in the plane strain condition. For the aluminum alloy, α were selected from 1.7 to 2 in general [26]. In this paper, α is assumed to be 2. p controls the shape in the threshold region and should be calibrated using the experimental data in the threshold region. For the reason that the crack growth data in the threshold region are not sufficient or even not available, an initial value of p is guessed. In the NASGRO database [23, 24], some typical values of p can be used here. In this paper, p is selected as 1 for Al2024-T3 and 1.5 for Al7075-T6. C1 is the material constant and C2 is the Paris exponent, and they can be calibrated using the crack propagation data in the Paris region.
The value of ΔKth is the key parameter for life prediction algorithm based on the LEFM and is often used in the fatigue crack growth analysis. It defines a value below which the fatigue crack will cease to grow or the growth rate is slow enough to be neglected. It is known that the crack nucleation and the small-crack growth consume most of the fatigue life [33]. Hence, the fatigue life is very sensitive to the threshold SIF range. The accuracy of the value of ΔKth will have a large impact on model predictions based on the LEMF. Experiments are the main methods to obtain the values of ΔKth in the current research, such as the load-reduction test method in ASTM E647 and the compression precracking method [34]. In other words, the value of ΔKth highly depends on the experimental techniques which may lead to quite different measurements of ΔKth. In the current work, the concept of “intrinsic effective threshold stress intensity factor” ΔKeffth is introduced considering the crack closure effect and is considered to be the material property which is independent of the stress ratios or experimental technique. Similar concept is also used in Newman’s fatigue life prediction model [26, 27, 35], but the value of ΔKeffth in this paper is not determined via trial and error and is obtained by extrapolating the crack growth rata curve, that is, da/dN~ΔKeff in the logarithm linear region to 10−10 m/cycle. The NASGRO equation is hence modified as(6)dadN=C1ΔKeffC21-ΔKeffth/ΔKeffp1-Kmax/Kcq.
Since only a very small part of fatigue life is consumed in the fast crack propagation stage, (6) is hence simplified as(7)dadN=C1ΔKeffC21-ΔKeffthΔKeffp.
Fatigue life N is expressed as(8)N=∫aiac1C1ΔKeffC21-ΔKeffth/ΔKeffpda,where ac is the critical crack size when the fatigue failure occurs. ac can be set as a constant in the numerical simulation, since the fatigue life consumed in the fast crack propagation stage is negligible. ai is the initial defect size and is often determined by solving the El Haddad equation [36]; the equation is expressed as(9)ΔKth=ΔσfπcF,where Δσf is the fatigue limit, c represents the surface crack length, and F represents the boundary correction factor. Multiply both sides of (9) by U, and (9) can be expressed as(10)ΔKeffth=ΔσfπcFU,where ΔKeffth=UΔKth and U=(1-f)/(1-R). ΔKeffth is the intrinsic effective threshold SIF range. ai in this paper is determined by solving (10).
In the current study, smooth and circular-hole specimens are used to validate the proposed method, as shown in Figure 1. The surface crack for a smooth specimen is modeled as a semicircular surface crack; that is, the crack depth is equal to the surface crack length. The boundary correction factor F for a surface crack is expressed as [37] (11)F=M0g1fφfwfx,where(12)M0=m1+m2at2+m3at4,m1=1.13-0.09ac,m2=-0.54+0.890.2+a/c,m3=0.5-10.65+a/c+141-ac24,g1=1+0.1+0.35at21-sinφ2,fφ=sinφ+ac2cos2φ0.25,fw=secπcWat,fx=1+1.464ac1.65-0.5,where a is the surface crack depth, 2c is the surface crack length, W is the plate width, and t is the plate thickness. In this paper, only crack growth in the surface crack tip is considered and its SIF range (φ=0) is used in (10). The surface crack for a circular-hole specimen is modeled as a quarter-circular corner crack, as shown in Figure 1, and the corresponding boundary correction factor is available in the literature [38].
Schematic representation of a semicircular surface crack and a corner crack from a hole.
For low-cycle fatigue condition, the plastic zone effect cannot be neglected. The SIF range considering the plastic zone effect is expressed as(13)ΔK′=ΔσfπdF′,where d=c+ρ, F′ is the boundary correction factor considering the plasticity correction, and the plasticity correction factor ρ is given by [32] (14)ρ=csecπσmax1-R4σ0-1,where σmax is the maximum stress and σ0 is the flow stress and is equal to the average of the ultimate strength and yield strength.
3. Model Validation for Smooth and Circular-Hole Specimens3.1. Experimental Data
The crack growth rate data of Al2024-T3 and Al7075-T6 under different stress ratios are collected from the literature [39]. They are used to calibrate parameters C1 and C2 in (8). Figure 2 shows the da/dN~ΔKeff relation of Al2024-T3 and Al7075-T6. For Al2024-T3, the calibrated parameters C1=3.5535E-11, C2=4.059, and the corresponding ΔKeffth is 1.2903 MPa. For Al7075-T6, C1=1.5779E-10, C2=3.477, and the corresponding ΔKeffth is 0.8771 MPa.
Crack growth data of Al2024-T3 and Al7075-T6.
Al2024-T3
Al7075-T6
The fatigue life test data of Al2024-T3 and Al7075-T6 are collected from the open literature [26, 40, 41], and the corresponding geometric sizes are listed in Tables 1 and 2. The fatigue properties for two different aluminum alloys used here are listed in Table 3. Δσs and Δσh are the fatigue limits of smooth and circular-hole specimens, respectively. EIFSs and EIFSh calculated by using (10) are the corresponding EIFS.
Specimen geometry of Al2024-T3.
Specimen types
R
Width (mm)
Thickness (mm)
Hole diameter (mm)
Smooth
0
25.4
2.3
/
−1
25.4
2.3
/
0.1
10.0
5.0
/
Circular-hole
0/−1
50.8
2.3
3.2
Specimen geometry of Al7075-T6.
Specimen types
R
Width (mm)
Thickness (mm)
Hole diameter (mm)
Smooth
0
25.4
2.3
/
−1
50.8
4.8
/
Circular-hole
0/−1
50.8
2.3
1.6
Mechanical properties and fatigue properties for Al2024-T3 and Al7075-T6.
Materials
R
σy (MPa)
σu (MP)
Δσs (MPa)
Δσh (MPa)
EIFSs (μm)
EIFSh (μm)
Al2024-T3
0
360
490
200
110
50.68
24.77
0.1
360
490
170
/
61.04
/
−1
360
490
285
110
73.19
86.57
Al7075-T6
0
520
575
205
120
22.95
9.75
−1
501
569
415
130
15.14
27.96
3.2. Life Prediction for Smooth Specimens
The predicted S-N curves using the proposed method are compared with experimental data under different stress ratios, as shown in Figure 3. The solid lines represent predicted S-N curves and the solid marks are experimental data. The predicted fatigue lives are consistent with experimental data.
Experimental data and calculated fatigue lives for smooth specimens.
Al2024-T3
Al7075-T6
3.3. Life Prediction for Circular-Hole Specimens
Figure 4 shows the calculated fatigue lives and test data. The predicted fatigue lives agree well with the experimental data.
Experimental data and calculated fatigue lives for circular-hole specimens.
Al-2024-T3
Al7075-T6
4. Discussion and Conclusions
Newman Jr. et al. [26] calculated the fatigue lives of two aluminum alloys using the FASTRAN code, and the life prediction model imbedded is Newman’s crack closure model. In Newman’s life prediction approach, the values of ΔKeffth and initial defect size, which are the key parameters in the LEMF-based life prediction method, are empirically selected via trial and error to fit the experimental data. In a sense, it is more experience-based rather than physical. In this study, the fatigue life is calculated by integrating the modified NASGRO equation from the initial defect size to the final crack size. The extrapolation method is used to determine the value of ΔKeffth, and the initial defect size is estimated using the EIFS concept. In Figure 5, the life predictions of the proposed model are compared with the experimental data and Newman Jr. et al’s results [26]. Solid lines are the numerical calculation results using the proposed algorithm, and the dashed lines are Newman’s prediction results. The calculated fatigue lives are a little conservative compared with Newman’s model predictions in the low-cycle fatigue region. In the middle and high cycle fatigue region, the prediction results of both methods are nearly the same. However, Newman’s life prediction algorithm is very complicated and experience-based. On the whole, the life prediction algorithm proposed in this paper is more simple and efficient.
The calculated fatigue lives and test data for smooth Al2024-T3 specimens.
Fatigue limit is the main parameter used to calculate the EIFS. However, studies show that the fatigue limit may not exist for metallic materials in the very high cycle fatigue (fatigue life longer than 107 cycles in general) [42, 43]. In this paper, the fatigue limit is determined based on the test data in the region where the fatigue life is larger than 106 cycles. If the experimental data show obvious horizontal tendency in this region and have little dispersion, the conservative data points are used in fitting the fatigue limit, such as the smooth Al7075-T6 specimens under R=-1. If the experimental data show larger discrepancies, the conservative data point is considered in determining the fatigue limit, and the weight should also be taken into account, such as the smooth Al7075-T6 specimens under R=0. The little change in the value of fatigue limit will not have a significant effect on model predictions. Figure 6 shows the calculated fatigue lives of smooth Al2024-T3 specimens under R=0 using different values of fatigue limit. It can be seen from Figure 6 that the life prediction results are almost the same.
Life predictions using different values of fatigue limit.
There are some discrepancies existing between the calculated fatigue lives and experimental data, especially under high stress range. One possible reason is that the coalescence of multiple cracks is predominant in the early crack propagation stage [44, 45], and the crack size may increase abruptly. The abrupt change in the crack size is not considered in the proposed method. Besides, the crack depth is assumed to be equal to the crack length which may also cause discrepancies.
In this study, a simple and efficient fatigue life prediction algorithm based on the LEFM is proposed. The crack propagation model is based on the modified NASGRO equation. The concept of intrinsic effective threshold SIF range ΔKeffth is introduced in the modified NASGRO equation. The value of ΔKeffth is determined by back-extrapolation. The fatigue limit and intrinsic effective threshold SIF range ΔKeffth are employed to estimate the initial defect size. The following conclusions can be drawn:
Only large crack propagation data are required in the proposed life prediction algorithm, and the algorithm can avoid the complicated analysis of small-crack propagation. A number of test data verified the effectiveness of the proposed life prediction algorithm.
ΔKeffth is considered as the material property which is independent of stress ratios, and it is obtained from large crack growth data by back-extrapolation.
The proposed algorithm is validated by two aluminum alloys (Al2024-T3 and Al7075-T6) under constant amplitude loading, and it can calculate fatigue lives of any structure materials if the corresponding accurate ΔK solutions are available.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The research is financially supported by the National Natural Science Foundation of China (no. 51405009) and the Fundamental Research Funds for the Central Universities.
BasquinO. H.The exponential law of endurance tests1910102625630CoffinL.TavernelliJ.The cyclic straining and fatigue of metals1959215794807MansonS. S.Fatigue: a complex subject—some simple approximations19655719322610.1007/bf023210562-s2.0-34250605357NinhN. T.WahabM. A.The effect of residual stresses and weld geometry on the improvement of fatigue life1995481–458158810.1016/0924-0136(94)01697-y2-s2.0-0029196419ShiraniM.HärkegårdG.Fatigue life distribution and size effect in ductile cast iron for wind turbine components2011181122410.1016/j.engfailanal.2010.07.0012-s2.0-78649358089TorresM. A. S.VoorwaldH. J. C.An evaluation of shot peening, residual stress and stress relaxation on the fatigue life of AISI 4340 steel200224887788610.1016/S0142-1123(01)00205-52-s2.0-0036680749WalkerK.The effect of stress ratio during crack propagation and fatigue for 2024-T3 and 7075-T6 aluminum1970West Conshohocken, Pa, USAASTM InternationalNorooziA. H.GlinkaG.LambertS.A study of the stress ratio effects on fatigue crack growth using the unified two-parameter fatigue crack growth driving force2007299–111616163310.1016/j.ijfatigue.2006.12.0082-s2.0-34548042098JiangS.ZhangW.LiX.SunF.An analytical model for fatigue crack propagation prediction with overload effect20142014971367810.1155/2014/7136782-s2.0-84934929641RokhlinS. I.KimJ.-Y.NagyH.ZoofanB.Effect of pitting corrosion on fatigue crack initiation and fatigue life1999624-542544410.1016/S0013-7944(98)00101-52-s2.0-0032987242CarterD. R.HayesW. C.Fatigue life of compact bone—I effects of stress amplitude, temperature and density197691273410.1016/0021-9290(76)90136-62-s2.0-0016891897StewartA. T.The influence of environment and stress ratio on fatigue crack growth at near threshold stress intensities in low-alloy steels198013346347810.1016/0013-7944(80)90078-82-s2.0-0019141739ParisP.ErdoganF.A critical analysis of crack propagation laws196385452853310.1115/1.3656900PriddleE. K.High cycle fatigue crack propagation under random and constant amplitude loadings1976428911710.1016/0308-0161(76)90014-42-s2.0-0016940344FormanR. G.KearneyV. E.EngleR. M.Numerical analysis of crack propagation in cyclic-loaded structures196789345946310.1115/1.3609637ElberW.The significance of fatigue crack closure1971486230243SureshS.ZamiskiG. F.RitchieD. R. O.Oxide-induced crack closure: an explanation for near-threshold corrosion fatigue crack growth behavior19811281435144310.1007/bf026436882-s2.0-0019599985GrayG. T.WilliamsJ. C.ThompsonA. W.Roughness-induced crack closure: an explanation for microstructurally sensitive fatigue crack growth198314242143310.1007/bf026442202-s2.0-0001294381ParisP. C.TadaH.DonaldJ. K.Service load fatigue damage—a historical perspective1999211S35S462-s2.0-0033227521NewmanJ. C.Jr.A crack-closure model for predicting fatigue crack growth under aircraft spectrum loading19817485384NewmanJ. C.Jr.A crack opening stress equation for fatigue crack growth1984244R131R13510.1007/BF000207512-s2.0-0021405273DugdaleD. S.Yielding of steel sheets containing slits19608210010410.1016/0022-5096(60)90013-22-s2.0-50549180512ConsortiumN.Fatigue crack growth computer program NASGRO Ⓡ version 3.0. user manual2001JSC-22267BFormanR.ShivakumarV.CardinalJ.WilliamsL.McKeighanP.Fatigue crack growth database for damage tolerance analysis2005DOT/FAA/AR-05/15Washington, DC, USAUS Department of Transportation Federal Aviation Administration (FAA), Office of Aviation ResearchNewmanJ. C.Jr.FASTRAN-2: a fatigue crack growth structural analysis program199292:30964NewmanJ. C.Jr.PhillipsE. P.SwainM. H.Fatigue-life prediction methodology using small-crack theory199921210911910.1016/S0142-1123(98)00058-92-s2.0-0033078998NewmanJ. C.Jr.Fatigue-life prediction methodology using a crack-closure model1995117443343910.1115/1.28047362-s2.0-0029386973NewmanJ. C.Jr.AnagnostouE. L.RuskD.Fatigue and crack-growth analyses on 7075-T651 aluminum alloy coupons under constant- and variable-amplitude loading20146213314310.1016/j.ijfatigue.2013.04.0202-s2.0-84897973919MillerK. J.De Los RiosE.1986Suffolk, UKMechanical Engineering PublicationsKitagawaH.TakahashiS.Applicability of fracture mechanics to very small cracks or the cracks in the early stageProceedings of the 2nd International Conference on Mechanical Behavior of Materials ASM1976Metals Park, Ohio, USA627631SchijveJ.1984Delft University of TechnologyLiuY.MahadevanS.Probabilistic fatigue life prediction using an equivalent initial flaw size distribution200931347648710.1016/j.ijfatigue.2008.06.005ZBL1273.745092-s2.0-56249130742SchijveJ.2001Berlin, GermanySpringerNewmanJ.Jr.SchneiderJ.DanielA.McKnightD.Compression pre-cracking to generate near threshold fatigue-crack-growth rates in two aluminum alloys20052710–121432144010.1016/j.ijfatigue.2005.07.0062-s2.0-27144434031NewmanJ. C.Jr.AnnigeriB. S.Fatigue-life prediction method based on small-crack theory in an engine material2012134303250110.1115/1.40042612-s2.0-84855406580El HaddadM. H.TopperT. H.SmithK. N.Prediction of non propagating cracks197911357358410.1016/0013-7944(79)90081-X2-s2.0-0018295787NewmanJ.Jr.RajuI.Stress-intensity factor equations for cracks in three-dimensional finite bodies subjected to tension and bending loads19862311334NewmanJ.Jr.RajuI.Stress-intensity factor equations for cracks in three-dimensional finite bodies1983791238265HuangX.TorgeirM.CuiW.An engineering model of fatigue crack growth under variable amplitude loading200830121010.1016/j.ijfatigue.2007.03.0042-s2.0-35548989760RodopoulosC. A.CurtisS. A.de Los RiosE. R.SolisRomeroJ.Optimisation of the fatigue resistance of 2024-T351 aluminium alloys by controlled shot peening—methodology, results and analysis200426884985610.1016/j.ijfatigue.2004.01.0032-s2.0-2342481721ZhaoT.JiangY.Fatigue of 7075-T651 aluminum alloy200830583484910.1016/j.ijfatigue.2007.07.0052-s2.0-38949160739BathiasC.There is no infinite fatigue life in metallic materials199922755956510.1046/j.1460-2695.1999.00183.x2-s2.0-0032626508PyttelB.SchwerdtD.BergerC.Very high cycle fatigue-is there a fatigue limit?2011331495810.1016/j.ijfatigue.2010.05.0092-s2.0-77956612866FriseP. R.BellR.Modelling fatigue crack growth and coalescence in notches199251110712610.1016/0308-0161(92)90008-42-s2.0-0026617102ParkinsR. N.The application of stress corrosion crack growth kinetics to predicting lifetimes of structures19892981019103810.1016/0010-938X(89)90091-72-s2.0-0024888951