Open access peer-reviewed chapter

Numerical Analysis Methods of Structural Fatigue and Fracture Problems

Written By

Qiu Zhiping, Zhang Zesheng and Wang Lei

Submitted: August 18th, 2017 Reviewed: November 7th, 2017 Published: December 20th, 2017

DOI: 10.5772/intechopen.72285

Chapter metrics overview

1,194 Chapter Downloads

View Full Metrics


Fatigue and fracture problems, which lead to 95% of structural failure, have attracted much attention of engineers and researchers all over the world. Compared with experimental method, numerical simulation method based on empirical models shows its remarkable advantages in structure design because of less cost and higher efficiency. However, the application of numerical simulation method in fatigue lifetime prediction is restricted by low accuracy and poor applicability in some circumstances. Most numerical method is based on empirical models. This chapter first reviews various kinds of empirical models of fatigue and fracture problems, including some modifying methods of basic empirical models, which have been widely applied to fatigue lifetime prediction and indicated their advantages and disadvantages. Then, FEM is introduced as an important method to obtain stress intensity factor or crack growth route. At last, this chapter is finished with existing problems and current trends in fatigue lifetime prediction via numerical method.


  • fatigue lifetime prediction
  • crack propagation
  • numerical method
  • empirical model
  • Paris law
  • perturbation approach
  • extended finite element method
  • fractal geometry

1. Introduction

With the development of mechanical engineering and manufacturing technology, engineering structures applied in aircrafts and huge machines become much more complex. These structures usually bear constantly changing loads in tour of duty. Although the max stress in structure caused by these dynamic loads is much lower than yield limit and ultimate strength of material, structure is destroyed after a long time. Internal defects in engineering structures appear in producing, processing, and assembling process. Internal defects lead to stress concentration, crack initiation, and propagation and even fatigue failure under dynamic load. According to statistical data, loss caused by improper structural fatigue lifetime design in America equals 4.4% of gross national product, and 95% of structure failures are related to fatigue break caused by alternating dynamic loads [1]. There are numerous historical examples that result in great loss of human life and economic value. For example, two Comet aircraft crashed in 1954, and the main reason is fatigue of fuselage structure [2]. Mechanical failure caused by fatigue, which concentrates much attention of engineers and researches, has been studied for more than 150 years [3]. However, it is still much difficult to prevent fatigue failure because fatigue of materials is far from being completely comprehended [4].

Metallic materials are widely applied in design of structures and parts in present days; therefore, fatigue of metals is a problem deserving efforts. In fact, the fatigue process is constitutive of crack initiation and crack propagation to total failure, as shown in Figure 1, and fatigue lifetime should conclude crack initiation life and crack propagation life.

Figure 1.

Schematic illustration of crack length versus time/cycles.

On one hand, it is widely accepted that the crack initiation phase costs a majority of fatigue lifetime in a high-cycle fatigue regime [5]. Furthermore, crack initiation behavior has a great influence on crack growth prediction in a unified approach for fatigue lifetime prediction [6]. Therefore, knowledge and technology of crack initiation life prediction are significant for evaluation of fatigue lifetime of structures and deserve our efforts to study deeply. On the other hand, there are frequently small cracks and defects in engineering structures due to manufacturing and environment factors; therefore, fatigue crack propagation prediction plays an important role in estimating the structural safety under dynamic loads.

Therefore, people divide structural life prediction problem into two problems: fatigue problem and fracture problem. People pay attention to crack initiation life in fatigue problem and make efforts to construct the relationship between structure life and stress or strain in structure. It is assumed in fatigue problem that there is a small crack existing in structure, and crack propagation behavior is studied in order to predict the remaining life of structure. These two problems have aroused widespread concern nowadays.

Experimental method and numerical method are two significant ways to analyze fatigue lifetime of structures. Experimental method has been widely applied since a long time ago. However, it is much expensive to predict structural life via experimental method. Furthermore, it is difficult to execute experiments for some complicated structures. Therefore, numerical method based on empirical models becomes much more popular in structural life prediction, and in some cases, those do not need high accuracy because of less cost and higher efficiency.


2. Empirical models in fatigue problem

Approaches to predict fatigue initiation life in literature can be classified into several types. These approaches study the fatigue problem from different perspectives, involving the average or local values of stresses and strains, the initiation of crack and defects, and macro- and microanalysis [7]. Nevertheless, people prefer to use phenomenological models, which reflect general material response at macroscopic scale under cyclic loads, rather than complex micro- or mesoscopic model of material fatigue behavior in structure design [8].

2.1. Empirical models of high-cycle fatigue

Wöhler is the pioneer in this field, who established the traditional stress-based approach in the nineteenth century [9]. He carried out a few fatigue experiments on metallic materials and indicated the relationship between fatigue crack initiation life and cyclic stress. He proposed to apply SNcurves in description of fatigue behavior of metals in his paper. Effectiveness of this method in high-cycle fatigue analysis is demonstrated afterward by many researchers. There are several kinds of expression of SNcurve, mainly including exponential function expression and power function expression. Basquin was the first person who suggested using exponential function to construct the expression of SNcurve in the twentieth century. The typical exponential function expression is written as follows:


where mand Care constants, which can be determined based on experiment data, Nstands for the number of loading cycles, and Smaxis the maximum value of stress at specific stress ratio. The power function expression with two parameters is usually expressed in the following form:


where Sais the stress amplitude at specific ratio. The power function expression with three parameters is expressed as




where D, A, and αare constants. The parameter Cin Eqs. (3) and (4) nearly equals fatigue limit.

2.2. Empirical models of low-cycle fatigue

Stress level is usually high in low-cycle fatigue, and the maximum value of stress is nearly close to the ultimate strength of material. The number of loading cycles in low-cycle fatigue, which is not more than 103 times, is much less than that in high-cycle fatigue. Plastic deformation plays an important role in low-cycle fatigue, in which the accumulation of plastic deformation results in structural failure. Because low-cycle fatigue lifetime is much sensitive to the change of stress level, SNcurve is unable to reflect the low-cycle fatigue performance of material. Therefore, εNcurve is applied to low-cycle fatigue analysis. The most widely accepted low-cycle fatigue lifetime model based on εNcurve is proposed by Basquin [10], which is expressed as follows:


where εeis the amplitude of elastic strain, Eis the elasticity modulus of material, σfis the fatigue strength coefficient of material, and bis the fatigue strength exponent. Because the relationship between plastic strain and fatigue lifetime is not taken into consideration in Basquin formula, Coffin [11] and Manson [12] proposed an empirical model when studying the relationship between fatigue lifetime and plastic strain amplitude. The expression of Coffin-Manson model is.


in which εastands for the amplitude of total strain and εfand cstand for the fatigue ductility coefficient and fatigue ductility exponent separately. The relationships between plastic strain, elastic strain, total strain, and fatigue lifetime are shown in Figure 2.

Figure 2.

Elastic strain-life curve and plastic strain-life curve.

2.3. Improved models considering mean stress or stress ratio

There are many factors, such as residual stress, temperature, multiaxial stress, and geometrical feature, that influence structural fatigue lifetime, in which mean stress or stress ratio concentrates the most attention.

2.3.1. Walker formula

Mean stress and stress ratio are of great significance for structural fatigue lifetime. Walker formula considers sensitivity of different materials to mean stress; therefore, it shows well effect for all materials [13, 14]. An equivalent local strain parameter is defined in Walker formula; its expression is


ris the material parameter. In order to construct the relationship between Walker formula and fatigue lifetime, Jaske et al. [15] carried out many experiments on different kinds of materials and proposed following expression based on experimental data:


where A0and A1are regression coefficients and εuand εeare the upper and lower limits of this reverse hyperbolic tangent function, respectively. The strain-life curve is shown in Figure 3. There are too many parameters to be fitted in this method, which need plenty of experimental data. That disadvantage constricts badly the application of Walker formula in engineering.

Figure 3.

Strain-life curve of Walker formula.

2.3.2. Morrow’s modifying method

Morrow’s modifying method and SWT modifying method are two commonly used methods. Morrow mean stress modifying formula is shown as follows [16]:


Considering the greater influence made by mean stress in long life period, further modifying method is given:


where εais strain amplitude and σmis mean stress. Morrow’s modifying method aims at elastic strain; therefore, it is only suitable when stress amplitude is constant or mean stress is compression stress.

2.3.3. SWT modifying method

Expression of Smith-Watson-Topper (SWT) parameter modifying method is [17]




SWT mean stress modifying method is not valid for compression mean stress, and it will obtain too conservative result when the stretching mean stress is large.

2.3.4. Goodman’s modifying method

We can acquire the fatigue limit points of material at different stress ratio r=σmin/σmaxunder infinite life requirement with the support of large amount of experimental data. Draw these points in rectangular coordinate system whose X-axis is mean stress σm=σmin+σmax/2and Y-axis is stress amplitude σa=σmaxσmin/2; thus, the fatigue limit curve is fitted based on these points. It is unpractical to carry out many experiments on all materials and structures in engineering, so we usually use a simplified straight line to replace the fatigue limit curve. Goodman simplified straight line, which is one of these straight lines, is widely accepted due to its simplicity and conservative estimation [18], as shown in Figure 4. Goodman simplified straight line can be expressed in the following relationship:


where Sestands for the fatigue strength of material and Sustands for the ultimate tensile strength of material. However, it has been proved that Goodman modifying method is only appropriate for low-ductility material, such as high-strength steel and cast iron.

Figure 4.

Fatigue limit curve and Goodman simplified straight line.


3. Empirical models in fracture problem

3.1. Paris law

Paris et al. [19] made great contribution in this field who was pioneer suggesting that crack growth rate, da/dN, was a function of the maximum stress intensity factor Kmaxin 1961. Then, Liu [20] related the crack growth to the stress intensity factor range ΔKsubsequently. Paris and Erdogan [21] proposed the well-known Paris law, which can be presented as follows:


where Cand mcan be obtained from experiment data, and they are usually considered as constants for a particular metal and environment [22]. Since then researchers have made efforts to study on Paris law and its deviation; however, we are still far from a complete comprehension [23].

It is believed that the relationship between crack propagation and ΔKcan be divided into three distinct regions, as shown in Figure 5. The crack propagation is slow in region A, and concept of a fatigue threshold stress intensity factor range ΔKthis proposed by Mcclintock [24], beneath which cracks are regarded not to grow. In region B, the “mid growth” range, crack propagation is stable, and Paris law is supposed to be held. Region C is associated with fast crack propagation leading to final failure. Therefore, calculation of number of loading cycles in region B, which could be gained from Paris law, is significant for prediction of fatigue crack growth life.

Figure 5.

Schematic diagram of the relationship between crack growth andΔK.

3.2. Improved models

3.2.1. Models considering mean stress or stress ratio

Since Paris law is proposed, much related work is done, and many modifying methods are put forward [22, 25, 26, 27]. It is commonly accepted that crack growth rate of material is related to mean stress or stress ratio. Several models, in which Forman formula [28] and Walker formula [29] are most famous, take this factor into consideration. Forman formula also considers the fracture toughness as an important factor; its expression is


Forman formula is valid for dealing with experimental data of many kinds of materials, especially high-hardness alloy, but it is hard to obtain the fracture toughness Kcfor high-ductility material. According to following relationship:


Forman formula can be transformed as follows:


Forman formula explains the reason why crack growth enlarges sharply when stress intensity factor is close to fracture toughness.

Walker formula is another wide-applied crack propagation model in engineering, which expresses the influence made by stress ratio on crack growth rate. Furthermore, it takes maximum of stress intensity factor into consideration:


Three parameters C, m, and ncan be acquired based on experimental data of crack propagation experiments with different stress ratios. Walker formula is valid when R>0and R<0. According to the relationship between stress ratio and amplitude of stress intensity factor, another commonly used form of Walker formula is obtained:


3.2.2. Model based on crack closure theory

In 1971, Elber [30] found that crack opened completely only when the stress was larger than a certain value, and he developed a modified Paris law based on this theory. The stress when crack is completely open is defined as crack opening stress σop, and the stress when crack begins to close is defined as crack closing stress σcl. It has been demonstrated that crack opening stress is nearly equal to crack closing stress. The modified formula is written as follows:




Uis the crack closure parameter, and its expression is


where efficient stress amplitude Δσeffis the difference between maximum stress σmaxand crack opening stress σop.

3.2.3. Model considering crack retardation caused by high load

In Weeler’s opinion [31], when structure bears cyclic load with constant amplitude; an occasional overload enlarges the size of plastic zone on crack tip, which would prevent crack from growing to some degree. On the basis of Weeler’s research, Willenberg [32] assumed that crack retardation is due to residual compression stress σres, which is related to plastic deformation caused by high load. Combining the expression of Forman formula, crack growth rate in retardation period is acquired:


The effective stress intensity factor range is


and the effective stress ratio is


The maximum and minimum values of effective cyclic stress are


Then, crack growth rate in retardation period can be estimated as the residual stress σresis known. However, the residual stress σrescan only be obtained via experimental method.

3.2.4. Model considering crack propagation threshold

In 1972, Donahue [33] took threshold of stress intensity factor range ΔKthinto consideration and proposed a generalized Paris law. The modified expression is


The following expression was proposed by McEvily and Greoeger [34] in their research about fatigue crack propagation threshold in 1977:


in which material constant mequals 2.

Furthermore, if considering stress ratio at the same time, Paris law can be modified into the following expression:


It can be figured out that the above equation is further modified on the basis of Forman formula.

In 1999, McEvily found it out that the following modification is suitable for many alloys’ fatigue crack propagation:


where ΔKeffthstands for the effective stress intensity factor range near crack propagation threshold. This modifying method considers the influences created by crack closure and small crack’s elastic-plastic behavior, and it is useful to predict the long crack propagation under cyclic positive stress.

3.2.5. Model based on perturbation series expansion method

Perturbation series expansion method, which is a common method to deal with nonlinear problems, has been widely used in fluid mechanics, structure dynamics, and damage identification. In this method, the parameter in ideal model is regarded to have a small perturbation in order to study the properties of system. This parameter can be expanded into series form:


where εis a positive small constant.

Qiu and Zheng [35] proposed a novel numerical calculation method to investigate the fatigue crack growth evolution in aluminum alloy sheets accounting for the measurement error. The initial crack length is considered as a modified parameter with a small correction term due to the measurement error; the solution to the crack growth equation is expressed in the form of a perturbation series, and a series of modified equations for predicting the crack length history are derived. The proposed method is verified to be indeed feasible and effective for predicting fatigue crack growth evolution by comparing numerical results with experimental data, as shown in Figure 6.

Figure 6.

Comparison of the measured and predicted crack length history in Ref. [35].


4. Finite element method

There are many kinds of numerical method to obtain stress intensity factor or crack growth route after continuous study of many researchers. Finite difference method (FDM), boundary element method (BEM), mesh-less method, and finite element method (FEM) are four common methods. Many studies have been carried out based on these numerical methods: Christen applied FEM to two-dimensional crack problem and obtained the displacement field and stress field; Nayroles [36] combined the moving least square method (MLSM) with mesh-less method to solve boundary problem. FEM is the most widely used method in above four methods at present [37, 38]. Considering singularity on crack tip, element’s density is increased in order to obtain the precious results. Therefore, FEM’s rate of convergence is low, and precision is unsatisfactory. People developed precious numerical solution methods based on several kinds of theories, in which semi-analytic numerical solution and new type elements are hot issues.

4.1. Extended finite element method

Collapsed singular isoparametric elements, which can reflect the singularity on crack tip correctly, were introduced by Barsoum [39]. This method is popular because of its high precision and executing simplicity. In this method, planar eight-node isoparametric element is degenerated into singular isoparametric element, as shown in Figure 7. Stress intensity factor is calculated based on the displacements of nodes A and B; the expression is.


Figure 7.

Eight-node singular isoparametric element.

In plane stress problem, E=E; in plane strain problem, E=E1μ2. E, μ, and vare, respectively, elasticity modulus, Poisson ratio, and displacement perpendicular to crack surface. Chen and Kuang [40] use interpolation method to acquire the displacements of nodes A and B on the basis of Barsoum’s research and obtain the following expression of stress intensity factor:


Lin [41] proposed the 1/4 node displacement method, as shown in Figure 8; the corresponding calculation equation of stress intensity factor is


Figure 8.

Mesh of 1/4 node element displacement method.

Belytschko [42] applied extended finite element method (XFEM) to calculating stress intensity factor and neglected the high-order terms of asymptotic displacement function. The calculation results were not satisfying enough. Karihaloo and Xiao [43] took high-order terms of asymptotic displacement function and outer elements of crack tip into consideration, thus obtaining results of high accuracy. However, calculation efficiency of this method is relatively low. Although researchers have obtained precious results with the help of new type elements, there are still many factors that influence calculation results that need to be studied.

4.2. Fractal finite element method

In the aspect of semi-analytic numerical method, weighted function method and boundary collocation method develop fast. These methods are able to acquire results of high accuracy when dealing with particular models; however, calculation accuracy cannot be guaranteed when dealing with general models.

Fractal finite element method is also a semi-analytic method. Fractal geometry is introduced into ordinary FEM, which not only improves calculation accuracy but also shortens calculation time and saves storage capacity of a computer. In fractal finite element method, an artificial boundary Γ0is introduced to divide the structure with crack into two parts: singular field Dnear crack tip and normal field Ωfar away from crack tip, as shown in Figure 9. Ordinary finite element mesh is constructed in normal field; self-similar mesh needs to be constructed based on fractal theory in singular field.

Figure 9.

Illustration of division of structure with crack.

Self-similar mesh is shown in Figure 10. In singular field, infinite curves Γ1Γ2Γ3similar to Γ0are generated based on the proportionally coefficient ξ0<ξ<1regarding crack tip as centre. The density of fractal mesh is controlled by ξ. Based on appropriate global interpolation function and fractal transforming technique, plenty of unknown degrees on slave nodes are transformed into a series of generalized coordinates. Stress intensity factor on crack tip can be calculated via solving generalized coordinates, thus saving calculation time and storage capacity obviously.

Figure 10.

Self-similar mesh in singular field.


5. Conclusion

This chapter reviews the most common empirical models and numerical methods of structural fatigue lifetime prediction. The main advantages and disadvantages of these methods are discussed.

Numerical method based on empirical models, as one of significant ways to analyze structural fatigue life, becomes popular in structural life prediction nowadays because of less cost and higher efficiency.

SNcurve and εNare applied to high-cycle and low-cycle fatigue problems, respectively. And there are many modified models considering mean stress or stress ratio. However, this chapter further shows that part of these models are too complicated to apply to engineering, and other models are only valid in some specific cases.

Paris law is the most significant model of crack propagation problem. But it only considers the stress intensity factor as the factors make influences on crack propagation. Many improved models considering stress ratio, crack closure, crack retardation, and crack propagation threshold have been put forward.

FEM is the most popular numerical method to obtain stress intensity factor or crack growth route. Extended finite element method and fractal finite element method are two mainly developing trends of FEM. However, it is still difficult to achieve high efficiency and accuracy of numerical method at the same time.


  1. 1. Jiang Y, Ding F, Feng M. An approach for fatigue lifetime prediction. Journal of Engineering Materials and Technology. 2007;129:182-189
  2. 2. Schijve J. Fatigue of aircraft materials and structures. International Journal of Fatigue. 1994;16(1):21-32
  3. 3. Starke P, Eifler D, Boller C. Fatigue assessment of metallic materials beyond strain measurement. International Journal of Fatigue. 2016;82:274-279
  4. 4. Gugno N, Ciavarella M, Cornett P, Carpinteri A. A generalized Paris' law for fatigue growth. Journal of the Mechanics and Physics of Solids. 2006;54:1333-1349
  5. 5. Zhao T, Jiang Y. Fatigue of 7075-T651 aluminum alloy. International Journal of Fatigue. 2008;30:834-849
  6. 6. Jiang Y, Feng M. Modeling of fatigue crack propagation. Journal of Engineering Materials and Technology. 2004;126(1):77-86
  7. 7. Vergani L, Guagliano M. From microstructure to design: Advances in fatigue of metals. International Journal of Fatigue. 2016;83:1
  8. 8. Karolczuk A, Kluger K, Lagoda T. A correction in the algorithm of fatigue lifetime calculation based on the critical plane approach. International Journal of Fatigue. 2016;83:174-183
  9. 9. Schijve J. Fatigue of structures and materials in the 20th century and the state of air. International Journal of Fatigue. 2003;25:679-702
  10. 10. Basquin OH. The experimental law of endurance tests. Proceedings of ASTM. 1910;10(II):625-630
  11. 11. Coffin Jr LF. A study of the effects of cyclic thermal stresses on a ductile metal. Transactions of ASME. 1954;76:931-950
  12. 12. Manson SS. Fatigue: A complex subject-some simple approximations. Experimental Mechanics. 1965;5(4):193-226
  13. 13. Dowling NE, Calhoun CA, Arcari A. Mean stress effects in stress-life fatigue and the Walker formula. Fatigue & Fracture of Engineering Materials & Structures. 2009;32(3):163-179
  14. 14. Duran JAR, Hernandez CT. Evaluation of three current methods for including the mean stress effect in fatigue crack growth rate prediction. Fatigue & Fracture of Engineering Materials & Structures. 2015;38(4):410-419
  15. 15. Jaske CE, Feddersen CE, Davis KB, Rice RC. Analysis of fatigue, fatigue crack propagation and fracture data. In: NASA CR-132332. United States. 1973. pp. 49-54
  16. 16. Morrow J. Fatigue design handbook. Advances in Engineering. 1968;4(3,2):21-29
  17. 17. Smith KN, Watson P, Topper TH. A stress-strain function for the fatigue of materials. Journal of Materials. 1970;5:767-778
  18. 18. Sutherland HJ, Mandell JF. Optimized Goodman diagram for the analysis of fiberglass composites used in wind turbine blades. Journal of Solar Energy Engineering. 2005;127(4):563-569
  19. 19. Paris PC, Gomez RE, Anderson WE. A rational analytic theory of fatigue. The Trend in Engineering. 1961;13(1):9-14
  20. 20. Liu HW. Crack propagation in thin metal sheet under repeated loading. Journal of Fluids Engineering. 1961;83(1):23-31
  21. 21. Paris PC, Erdogan F. Critical analysis of crack growth propagation laws. ASNW Transactions, Journal of Basic Engineering. 1963;85D(4):528-534
  22. 22. Jones R, Molent L, Pitt S. Similitude and the Paris crack growth law. International Journal of Fatigue. 2008;30(10):1873-1880
  23. 23. Pugno N, Ciavarella M, Cometti P, Carpinteri A. A generalized Paris’ law for fatigue crack growth. Journal of the Mechanics and Physics of Solids. 2006;54:1333-1349
  24. 24. Mcclintock FA. On the plasticity of the fatigue cracks. Fracture of Solids. 1963;20(01):65-102
  25. 25. Branco R, Antunes FV, Costa JD, Feng PY, Zhen BK. Determination of the Paris law constants in round bars from beach marks on fracture surfaces. Engineering Fracture Mechanics. 2012;96:96-106
  26. 26. Rajabipour A, Melchers RE. Application of Paris' law for estimation of hydrogen-assisted fatigue crack growth. International Journal of Fatigue. 2015;80:357-363
  27. 27. Carrascal I, Casado JA, Diego S, Lacalle R, Cicero S, Alvarez JA. Determination of the Paris' law constants by means of infrared thermographic techniques. Polymer Testing. 2014;40:39-45
  28. 28. Forman RG. Numerical analysis of crack propagation in cyclic loaded structures. Journal of Basic Engineering, Transaction ASTM (Series D). 1967;89:459-465
  29. 29. Walker EK. 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. ASTM STP462, West Conshohocken; 1970. pp. 1-14
  30. 30. Elber W. The significance of fatigue crack closure. Damage Tolerance in Aircraft Structure. 1971:230-242
  31. 31. Weeler OE. Spectrum loading and crack growth. Journal of Fluids Engineering. 1972;94(1):181-186
  32. 32. Willenberg J, Engle RM, Wood HA. A crack growth retardation model using an effctive stress concept. AFFDL-TM-71-1-FBR:A715659
  33. 33. Donahue RJ et al. Crack opening displacement and the rate of fatigue crack growth. International Journal of Fracture Mechanics. 1972;8:209-219
  34. 34. McEvily AJ, Greoeger J. On the threshold for fatigue crack growth. Advances in Research on the Strength & Fracture of Materials. 1978: 1293-1298
  35. 35. Qiu Z, Zheng Y. Predicting fatigue crack growth evolution via perturbation series expansion method based on the generalized multinomial theorem. Theoretical and Applied Fracture Mechanics. 2016;86:361-369
  36. 36. Nayroles B, Touzoel G, Villon P. Generalizing the finite element method: Diffuse approximation and diffuse elements. Computational Mechanics. 1992;10(5):307-318
  37. 37. Zheng X, Cui H, Su X, Engler-Pinto Jr CC, Wen W. Numerical modeling of fatigue crack propagation based on the theory of critical distances. Engineering Fracture Mechanics. 2013;114:151-165
  38. 38. Proudhon H, Li J, Wang F, Roos A, Chiaruttine V, Forest S. 3D simulation of short fatigue crack propagation by finite element crystal plasticity and remeshing. International Journal of Fatigue. 2016;82:238-246
  39. 39. Roshdy SB. Triangular quarter-point elements as elastic and perfectly-plastic crack tip elements. International Journal for Numerical Methods in Engineering. 1977;11(11):85-98
  40. 40. Chen LS, Kuang JH. A modified linear extrapolation formula for determination of stress intensity factors. International Journal of Fracture. 1992;54(1):3-8
  41. 41. Lin XB, Smith RA. An improved numerical technique for simulating the growth of planar fatigue cracks. Fatigue and Fracture of Engineering Materials and Structures. 1997;20(10):1363-1373
  42. 42. He WJ, Lin Y, Ding HJ. A three-dimensional formula for determining stress intensity factors in finite element analysis of cracked bodies. Engineering Fracture Mechanics. 1997;57(4):409-415
  43. 43. Xiao QZ, Karihaloo BL. Direct evaluation of accurate coefficients of the linear elastic crack tip asymptotic field. Fatigue & Fracture of Engineering Materials & Structures. 2003;26(8):719-729

Written By

Qiu Zhiping, Zhang Zesheng and Wang Lei

Submitted: August 18th, 2017 Reviewed: November 7th, 2017 Published: December 20th, 2017