Open access peer-reviewed chapter

# Numerical Analysis Methods of Structural Fatigue and Fracture Problems

By Qiu Zhiping, Zhang Zesheng and Wang Lei

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

DOI: 10.5772/intechopen.72285

## Abstract

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.

### Keywords

• 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 . 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 . Mechanical failure caused by fatigue, which concentrates much attention of engineers and researches, has been studied for more than 150 years . However, it is still much difficult to prevent fatigue failure because fatigue of materials is far from being completely comprehended .

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.

On one hand, it is widely accepted that the crack initiation phase costs a majority of fatigue lifetime in a high-cycle fatigue regime . Furthermore, crack initiation behavior has a great influence on crack growth prediction in a unified approach for fatigue lifetime prediction . 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 . 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 .

### 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 . 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:

emSmaxN=CE1

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:

SamN=CE2

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

SmaxCmN=DE3

or

Smax=C1+ANαE4

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 , which is expressed as follows:

εe=σaE=σfE2NfbE5

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  and Manson  proposed an empirical model when studying the relationship between fatigue lifetime and plastic strain amplitude. The expression of Coffin-Manson model is.

εa=σfE2Nfb+εf2NcE6

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.

### 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

εeq=2εarσmaxE1rE7

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

logNf=A0+A1tanh1logεuεeεeq2logεuεeE8

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.

#### 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 :

εa=σfE1σmσf2Nfb+εf1σmσf2NfcE9

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

εa=σfE1σmσf2Nfb+εf2NfcE10

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 

σmaxεa=σf′2E2Nf2b+σf2Nfb+cE11

where

σmax=σm+σaE12

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 , as shown in Figure 4. Goodman simplified straight line can be expressed in the following relationship:

σaSe+σmSu=1E13

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.

## 3. Empirical models in fracture problem

### 3.1. Paris law

Paris et al.  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  related the crack growth to the stress intensity factor range ΔKsubsequently. Paris and Erdogan  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 . Since then researchers have made efforts to study on Paris law and its deviation; however, we are still far from a complete comprehension .

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 , 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  and Walker formula  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:

R=KminKmaxE16
ΔK=KmaxKminE17

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  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:

and

Uis the crack closure parameter, and its expression is

U=ΔKeffΔK=ΔσeffΔσ=σmaxσopΔσ<1E23

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 , 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  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

ΔKeff=fσmaxeffσmineffπaE25

and the effective stress ratio is

Reff=σmineff/σmaxeffE26

The maximum and minimum values of effective cyclic stress are

σmaxeff=σmaxσresE27
σmineff=σminσresE28

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  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  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:

a=i=0aiεiE33

where εis a positive small constant.

Qiu and Zheng  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. .

## 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  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 . 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.

KI=E42πL4vAvBE34

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  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:

KI=E122πL8vAvBE35

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

KI=E22πLvAE36

Belytschko  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  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.

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.

## 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.

## More

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Qiu Zhiping, Zhang Zesheng and Wang Lei (December 20th 2017). Numerical Analysis Methods of Structural Fatigue and Fracture Problems, Contact and Fracture Mechanics, Pranav H. Darji and Veera P. Darji, IntechOpen, DOI: 10.5772/intechopen.72285. Available from:

### Related Content

#### Contact and Fracture Mechanics

Edited by Pranav H. Darji

Next chapter

#### Accelerated Fatigue Test in Mechanical Components

By Moises Jimenez

#### Related Book

Edited by Pranav H. Darji

First chapter

#### Manipulation of Tribological Properties of Metals by Ultrashort Pulsed Laser Micro-/Nanostructuring

By Quan-Zhong Zhao and Zhuo Wang

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books