Material properties of aluminum alloy and alumina.

## Abstract

In this chapter, extended finite element method (XFEM) has been used to simulate the fatigue crack growth problems in functionally graded material (FGM) in the presence of hole, inclusion and minor crack under elastic and plastic conditions. The fatigue crack growth analysis of alloy/ceramic FGMs, alloy and equivalent composite is done by XFEM in the presence of multiple discontinuities under mode-I mechanical load. The validity of linear elastic fracture mechanics (LEFM) theory is limited to the brittle materials. Therefore, the elastic plastic fracture mechanics (EPFM) theory needs to be utilized to characterize the plastic behavior of the material. A generalized Ramberg-Osgood material model has been used to model the stress-strain behavior of the material. Plasticity has been checked by Von Mises Yield criteria. J-integral has been used to calculate the SIF. Crack growth direction is determined by maximum principal stress criteria.

### Keywords

- FGM
- composite materials XFEM
- elastic-plastic loading
- fatigue fracture
- crack propagation
- discontinuities
- inclusions
- holes
- minor cracks

## 1. Introduction

Development of novel materials improves performance and efficiency of the structures, and also leads to development of advanced and sophisticated structures. This complex process of materials, structures and technology has led to the development of composite materials. Strength and stiffness plays a key role in evaluating the worth of the material. These characteristics provide strength to the structure to retain its desired shape and size under loading or any other external action.

Cracks/flaws are inevitable in all engineering materials. Loading under severe environmental conditions may either initiate new cracks or may cause the propagation of pre-existing cracks in the structures. Theoretically, fracture can be defined as the breaking or rupturing of a material resulting into its separation into two or more pieces.

Composite materials manifested in the middle of the twentieth century. Composites are naturally occurring or engineered materials made from two or more constituents with different chemical or physical properties distinct boundary among constituents. Lightweight composite materials with high strength to weight and stiffness to weight ratios have been used successfully in aircraft industry and other engineering applications. Under high temperature conditions the strength of the metal is deteriorated whereas, ceramics have excellent resistance to heat.

FGMs can be referred as multiphase composite materials in which the composition or microstructure or both are spatially varied which lead to a certain gradation in the local material properties. FGMs can be defined as multi-phase composites. FGMs are synthesized such that they own continuous variations in volume fractions of their components in space to return a pre-established composition. FGMs possess continuously varying properties in one or more than one direction and the form non-homogeneous macrostructure due to these variations. By gradually varying the volume fraction of the constituents, FGMs exhibit a smooth and continuous change from one surface to another, thus reducing interface problems, and minimizing thermal stress concentrations. The ceramic phase of FGMs provides a good resistance to heat, while the metal phase provides a strong mechanical performance and hence reduces the possibility of catastrophic failure.

The major advantages of FGM over conventional materials are firstly, FGM satisfies the working conditions for which it is specifically developed. Secondly, it is economical as it reduces material costs for particular engineering applications. Thirdly, it can reduce the magnitude of residual and thermal stresses generated under working conditions. Finally, FGMs exhibit better fracture toughness and bond strength. This is normally achieved by using a ceramic layer connected with a metallic layer. FGMs have wide area of engineering applications like in the computer circuit and aerospace industries. FGMs have typical applications is in aircraft and automotive industries as thermal barrier coatings (TBCs).

In general, all structural components are subjected to thermo-mechanical cyclic load. The fatigue life of these components is generally predicted without considering the effect of defects/discontinuities present in component. However, FGMs are commonly made by sintering process, which are porous in nature. These discontinuities at the vicinity of a major crack tip lead to increase the effective SIF at the major crack tip due to which the life of the components get depreciated. Hence, the analysis of FGMs in the vicinity of discontinuities becomes very important from the design point of view. To widen the spectrum of applications of FGMs, the fatigue/fracture behavior should be properly evaluated.

Over the years, greater understanding of fracture mechanics has undoubtedly prevented a significant number of structural failures. Fracture mechanics approach for the design of structures includes flaw size as one of the key variables. Fracture toughness replaces strength of material as a relevant material attribute, and its evaluation is mainly done in composites using the

Fatigue and quasi-static fracture are two forms of crack growth phenomenon. Fatigue fracture refers to the slow propagation of cracks under cyclic loading conditions where the stress intensity factors are below the fracture toughness of the material. Quasi-static fracture is observed near the end of the fatigue life when the increased crack length leads to stress intensity factors which are above the fracture toughness [2, 3].

In many cases, multiple cracks may exist in the components. There interaction resulting in the variation of stress intensity factor, stress distribution and propagation direction of the major crack. In the past, the failure of structures was analyzed in the presence of multiple cracks [4]. Some efforts have been made using analytical, experimental and simulation techniques to analyze the effect of interaction among multiple cracks [5, 6, 7].

Although, many analytical [8, 9] and experimental methods [10, 11] have been explored for the calculation of fracture parameters even then the drawbacks associated with experimental investigation and scarcity of analytical solution have impelled the analysts towards alternative techniques. Numerical methods hold the promise in this regard.

Many numerical methods are available to simulate the problems of fatigue failure in materials. These include finite element method (FEM), boundary element method, hybrid boundary node method [12, 13], meshfree methods [14, 15, 16, 17, 18] and extended finite element method [19, 20]. Out of these methods, FEM has been widely used for solving a wide variety of engineering and industrial problems [21, 22, 23, 24, 25, 26]. It has achieved a remarkable success in solving various linear and non-linear problems [27, 28, 29, 30, 31, 32, 33, 34]. Despite its numerous advantages and unparalleled success, it is not well-suited for solving the problems involving crack propagation. In crack growth problems, element edges provide natural lines along which a crack can grow. This is advantageous if the crack path is known a priori, but in most of the fracture phenomenon, the crack path is unknown. Thus, FEM requires a conformal mesh and re-meshing to ensure that the element boundaries coincide with the moving discontinuities (crack). Moreover, crack tip singularity cannot be accurately modeled by standard finite element approximation. Therefore, the modeling of crack growth becomes quite tedious and time consuming due to the modification in mesh topology at each stage of crack propagation. To overcome this difficulty, a new method known as extended finite element method (XFEM) has been developed to model arbitrary discontinuities without a need of conformal mesh or re-meshing. Level set method (LSM) is used in conjunction with XFEM for defining as well as tracking the geometry of cracks and other discontinuities like holes and inclusions. To cope up with these problems, XFEM has been adopted as a tool for the analysis of fatigue crack propagation in FGM.

## 2. Calculation of SIF for FGM

A domain based interaction integral approach can be used for calculating the stress intensity factors for homogeneous, bi-layer and functionally graded materials under thermal as well as mechanical loading. In this chapter, interaction integral approach will be extended to calculate the SIFs for FGM and bi-layered FGM under mechanical loads. The interaction integral is calculated based on *J*-integral. The *J*-integral for an elastic body subjected to thermo-mechanical load is given as,

For the interaction integral calculation of an elastic body, consider two equilibrium states i.e. state 1, the actual state with given boundary conditions and state 2, an auxiliary state of the cracked body. The parameters for auxiliary state are represented with superscript

where, the auxiliary field for the FGM may be taken from [36] as.

The SIFs are calculated from the interaction integral as [36]:

Mode-I SIF is given as,

Mode-II SIF is given as,

where,

## 3. Fatigue crack growth

Here we use Paris law for stable crack propagation, the generalized Paris’s law is given as:

Where,

The solution of Eq. (6) gives

According to this criterion, the equivalent mode-I SIF is obtained as

For stable crack propagation, the generalized Paris’ law for FGM is given as

where,

In the numerical example, the crack growth value Δ*a* is assumed and the corresponding number of cycles Δ*N* is computed from Eq. (9). When multiple crack tips are present, the crack growth value Δ*a* is assumed for the most dominant crack tip, corresponding Δ*N* is computed and then at the other crack tips the crack growth is computed corresponding to Δ*N*. Eventually, when the maximum value of

## 4. Modeling of the properties of FGM

In this chapter, the results have been presented for a FGM plate as shown in Figure 1. The FGM plate is manufactured by reinforcing an alloy with ceramic. The volume fraction of ceramic is varied in the

Material properties | Aluminum alloy | Alumina |
---|---|---|

Elastic modulus | 70 | 300 |

Poisson’s ratio, | 0.33 | 0.21 |

Coefficient of thermal expansion ^{°}C) | 25 × 10^{−6} | 8.2 × 10^{−6} |

Fracture toughness | 29 | 3.5 |

Paris law parameter | 10^{−12} | 2.8 × 10^{−10} |

Paris law parameter, | 3 | 10 |

The variation of the elastic modulus for FGM is modeled as

A plot of *L* = 100 mm is shown in Figure 2. The fatigue life of FGM has been compared with the same of the aluminum alloy and an equivalent composite of aluminum alloy/alumina. The equivalent composite considered in this example has the same overall volume fractions of aluminum alloy and ceramic as the FGM. The volume fractions of ceramic and aluminum alloy in the FGM are obtained as

In this example, the equivalent composite is assumed to have the same amount of metal and ceramic. The volume fraction of alumina in the equivalent composite is calculated as

where,

we get

The Poisson’s ratio is shown in Figure 4.

The coefficient of thermal expansion for the FGM by the rule of mixtures is calculated as

The coefficient of thermal expansion for the equivalent composite is given by

The value of the coefficient of thermal expansion for the equivalent composite is calculated using Eq. (14b) and is found to be ^{−6}/^{°}C. A variation of coefficient of thermal expansion for the FGM is shown in Figure 5.

The fracture toughness of the FGM as well as the equivalent composite may be expressed as a function of the volume fraction of the ceramic by the following formula given by [41]

The variation of fracture toughness is shown in Figure 6.

The Paris law parameters are assumed to have exponential variation in a manner similar to the elastic modulus. Thus, the variation in the parameters of Paris equation is taken as

For the equivalent composite, we find the location at which the volume fraction of ceramic in the FGM is same as that of the equivalent composite. This location

The values of *C* and

Plastic behaviour for FGM can be modeled using Ramberg Osgood equation [42]

Here,

The crack tip extension at the principal crack tip is

where

The constitutive relation for the elastic-plastic material is given as

where x is the vector of *x* and *y***-**coordinates,

and for plane strain condition as

The incremental theory of plasticity [43] has been used to model the elastic-plastic constitutive relation for a material. An incremental stress vector

Total strain increment is the sum of elastic and plastic strains

Elastic incremental strain and stress is determined

where

Plastic modulus

For a given strain energy

According to the Von Mises criteria,

After taking the derivatives from both sides of failure criteria equation

For simplicity we take

## 5. XFEM: Introduction and formulation for cracks and discontinuities in FGM

XFEM or the extended finite element method is a numerical technique which allows crack modeling irrespective of the mesh, and eliminates the cumbersome process of remeshing in problems involving change in the crack geometry like crack growth. XFEM models a crack by enriching the standard finite element approximation with some functions, which are obtained from the theoretical background of the problem. Moving discontinuities are tracked by the level set method. XFEM is a numerical method, based on the finite element method (FEM) that is especially designed for treating discontinuities. The formulation is done as discussed in [35, 44]. The solution of FGM differs from homogeneous materials only in the spatial gradation in the material properties. After calculating the values of stress and strain, the SIF is determined.

## 6. Numerical examples and discussion

The FGM plate considered in all the numerical simulations has 100% aluminum alloy on one side and 100% alumina on the other side. The volume fraction of alumina changes from 0% on one side to 100% on the other side so as to produce an FGM. The equivalent composite is equivalent to the FGM in the sense that both the FGM and the composite plate contain the same amount of aluminum alloy and alumina. The fatigue crack growth analysis of alloy/ceramic FGMs, aluminum alloy and equivalent composite is done by XFEM in the presence of multiple cracks, holes and inclusions under mode-I mechanical load and their fatigue life are compared. The constituents of the FGM plate are aluminum alloy and alumina. A major crack of large initial length is assumed to exist at the edge of the plate. The major crack is assumed to be in the direction of material gradation. The fatigue crack growth analyses of the FGM, the equivalent composite and the aluminum alloy plates have been carried out in the presence of minor cracks, holes and inclusions till the final failure of the plate under mode-I mechanical load. The effect of these small defects on the fatigue life as well as on the crack path has been investigated in detail.

### 6.1. Example 1

A rectangular FGM plate of length (*a* = 20 mm is analyzed under plane strain condition in the presence of multiple discontinuities. In all simulations, the plate dimensions, initial crack length and material properties are taken to be same. The properties of FGM, composites and aluminum alloy are already described in Table 1. The material properties of the inclusions are taken as *x*-direction and 235 nodes in *y*-direction. The fatigue crack growth analysis is performed by taking a crack increment of

### 6.2. Plate with a major edge crack under linear elastic condition

Figures 9 and 10 show a plate with a major edge crack of length *a* = 20 mm at the left and right edge respectively. These plates have been analyzed under plane strain condition using a uniform mesh of 117 by 235 nodes. The plots of the fatigue life for different materials are shown in Figure 11. From these figures, it is seen that the equivalent composite withstands 7885 cycles before it fails while the FGM with crack on alloy side undergoes 15,561 cycles and pure aluminum alloy undergoes 19,145 cycles before failure. It is also observed that when a major crack initiates from the ceramic (alumina) rich side then it fails much earlier (4872 cycles) as compared to when the crack initiates from the aluminum alloy side.

These plots show that when a crack is present on the ceramic rich side, the life diminishes by a considerable extent as compared to when a crack is present on the alloy rich side. The equivalent composite shows the minimum life except in case when a crack is present on the ceramic side. It is also observed that the crack follows nearly a straight path in all the materials.

### 6.3. Plate with a major edge crack, minor cracks, holes and inclusions under linear elastic condition

In this case, a major crack of length a = 20 mm. is taken at the left and the right edge of the plate (^{°} randomly. The holes and inclusions have variations in their radii from 3 to 4.5 mm randomly. A cyclic mode-I mechanical load is applied at the top edge of the plate. The plots for crack extension with number of cycles are shown in Figure 14.

It is also observed that the crack deflects in all the materials. Moreover, it is seen that the number of cycles to failure in case of aluminum alloy is about 18,111 cycles whereas in case of FGM with crack on the alloy and ceramic rich sides is 14,622 cycles and 3111 cycles respectively. The fatigue life of the composite plate is found to 6956 cycles. Thus, it can be stated that due to the presence of minor cracks, holes and inclusions, the life of the aluminum alloy is reduced by about 5.42%, whereas the fatigue life of the FGM with crack on the alloy and ceramic rich sides goes down by 6.03 and 36.15% respectively. The fatigue life of the equivalent composite is reduced by 11.78%.

### 6.4. Example 2

A rectangular FGM plate of length (

Material properties | Values |
---|---|

Elastic modulus of copper nickel alloy | 160 |

Elastic modulus of alumina (ceramic) | 386 |

Elastic modulus of soft inclusion | 100 |

Elastic modulus of Hard inclusion | 400 |

Poisson’s ratio of copper nickel alloy | 0.35 |

Poisson’s ratio of alumina (ceramic) | 0.21 |

Poisson’s ratio of inclusion | 0.3 |

Poisson’s ratio of inclusion | 0.23 |

Fracture toughness of copper nickel alloy | 79 |

Fracture toughness of alumina (ceramic) | 5 |

Paris constant | 3 × 10^{−11} |

Paris exponent | 3 |

### 6.5. A major crack in FGM plate under elastic: Plastic loading condition

In this case, a major crack of length *×* 200 mm) as shown in Figure 15. Cyclic loading is applied at the top edge of the FGM plate, and a crack propagates due to this loading. The plots of SIF with crack length for an crack configuration is shown in Figure 16. The failure crack length obtained for edge crack is 0.0402 m.

### 6.6. A major edge crack in FGM plate with holes, inclusions and minor cracks under elastic-Plastic loading condition

In this case, a major crack of length

## 7. Conclusions

In this chapter we have discussed the simulation of cracks in a FGM plate has been carried out in the presence of multiple inhomogeneities by XFEM using both linear elastic as well as elastic-plastic formulations. SIF has been calculated at the tip of the major crack using interaction integral approach. The variation in the SIF at the tip of the major crack has been studied when multiple inhomogeneities are present in the domain. From this study it is observed that minor cracks have least effect in the FGM plate’s failure crack length, whereas soft inclusions have moderate effect and holes have the most severe effect. It is found that the FGM plate’s life increases in each case when soft inclusions are replaced by hard inclusions. Hence the presence of the hard inclusions in the plate increases the failure crack length of the plate i.e. plate survives more.

## Nomenclatures

W˜: | Strain energy |

σ: | Stress |

ε: | Strain |

u: | Deformation |

C: | Compliance matrix |

E: | Modulus of elasticity |

γ: | Coefficient of thermal expansion |

KIC: | Critical stress intensity factor (Fracture toughness) |