Abstract

The problem of the static stress intensity factor around antiplane crack in an infinite strip functionally graded material was studied by using the method of integral transform-dual integral equations in this paper. The shear modulus in the two principal directions of the functionally graded material was assumed to vary proportionately as gradient model of double parameters index function. The partial differential equation was first reduced to Euler equation with Fourier cosine transform. By solving dual integral equations that were derived by applying the solution of Euler equation with the method of Copson, stress intensity factor around the crack tip was derived. And the variation curves of the dimensionless stress intensity factor with the strip height, crack length, gradient parameter, and inhomogeneous coefficient are obtained by using the numerical calculation.

1. Introduction

Functionally gradient material is a new type of materials. The advent of functionally gradient material is in order to meet performance requirements of parts in all kinds of special environment. The sudden changes of the material composition and properties in a member often lead to obvious local stress concentration. If a material transition to another kind of material is the step by step, the stress concentration will greatly reduce. A new type of functional gradient materials (FGM) was developed in order to reduce the stress concentration of the material and improve the performance requirements of the material as discussed elsewhere [15]. Most previous studies about FGM in the past were focused on isotropic material as discussed elsewhere [6, 7]. However, a few FGMs were isotropic material. At present, the application of FGM has been expanded to space, energy, traffic, optics, chemistry, biomedical engineering, and so forth. However, due to the influence of the process conditions, FGM will produce some defects in manufacturing process. These defects will inevitably form crack due to stress action and crack expansion will affect the whole member. So, the research about crack and the crack expansion of FGM are very important as discussed elsewhere [811].

In the fracture mechanics study of FGM, most scholars use the power functions model as a material property parameter model. The model of as discussed by Smith [1], the model of as discussed by Smith [12], and the model of as discussed by Smith [13]. However, the models of double parameters index function were adopted in this paper.

Usually the method of studies about stress intensity factor of antiplane crack in FGM is integral transforms-dual-integral equations. At present, most scholars put the partial differential equation into the Bessel equation with Fourier cosine transform in a designed material property parameter model. Then dual integral equations were obtained by using the solution of Bessel equation. But the difference between this paper and above references is that the partial differential equation is first reduced to Euler equation with Fourier cosine transform, and then dual integral equations were obtained by using the solution of Euler equation. Finally the local stress field and stress intensity factor around the crack tip were derived. This method is very few in the study of FGM.

2. Formulation of the Problem

Consider the antiplane problem of an infinite strip functionally graded material containing a crack of length , and the thickness is , as shown in Figure 1. The -axis and -axis are the elastic principal axis.

Shear moduli , are the functions of , and , vary according to the following functions: where parameter , , the shear modulus and , can be determined as .

Suppose that , , and are displacement components of the direction , and . For the antiplane shear deformation, , are 0 everywhere.

Assume that and is the displacement function. The relationship between the displacement function and stress components and    is written as

The stress equilibrium equation can be expressed as

We can obtain the governing equation [1416] of antiplane problem in infinite strip functionally graded materials by substituting stress components (2) into (3):

The corresponding boundary conditions of the crack surface under antiplane shear loading are as follows: Therefore, an analysis of antiplane problem in infinite strip functionally graded materials is to solve the boundary value problem of partial differential equation (4)–(7).

3. Dual Integral Equation

Due to the symmetry, we only need to study the part of , . Introducing the following Fourier cosine transform: where is called Fourier cosine transform of , and is called Fourier cosine inverse transform of .

Define Among them, is the Fourier integral cosine transformation function of displacement function .

By using the Fourier integral transformation, (4) can be converted into as follows: where , and inhomogenous coefficient .

Choosing proper variable substitutions as follows:

So

By substituting (12) into (10), partial differential equation (10) can be translated into Euler equation as follows: where

Introducing transformations of the independent variable

So, (13) can be transformed into the following formula:

The general solution of Euler equation (16) can be written as where and are the coefficients of solution.

Substituting (17) into (11), we can obtain the formula of solution as follows: where

From (9) and (18), the integral expression of displacement function was obtained as follow:

The shear stresses and can be obtained from (20) and the relationship of stress-displacement (2). The stress formulas are listed as follows:

We can obtain the relation between and by substituting (22) into the boundary conditions (7) as follows:

Let

From (20), (21), (24), and boundary conditions (5) and (6), we can obtain the dual integral equations:

in which,

Solving the dual integral equations (25a) and (25b) by using Copson-Sih method [17], the solution can be expressed as where is the so-called zero Bessel function of the first kind. The function can be governed by the Fredholm integral equation of the second kind: where kernel function is By using numerical integral method, the numerical solution of can be given from (28) and (29).

4. The Stress Intensity Factor

By using of the formula of integration by parts, formula (27) can be written as

We can obtain the stress principle by substituting (30) into (21)–(24) as follows: In formula (31),

It is noticed that

Considering the stress singularity around crack tip as , the stress principle can be written as Introducing the complex transformations as follows: For The formula (35) can be simplified to Defining the stress intensity factor as follows: So, we can obtain Formula (41) can be written as where is called dimensionless stress intensity factor.

The variation curves of the dimensionless stress intensity factor with the strip height , crack length , gradient parameter are and , and inhomogenous coefficient are obtained by using the numerical calculation as shown in Figures 24.

Figure 2 shows the variation curves of the dimensionless stress intensity factor with strip height and crack length under the same condition of inhomogeneous coefficient and gradient parameter. It is noticed that the longer crack and the thinner strip lead to the greater the stress intensity factor.

Figure 3 shows the variation curves of the dimensionless stress intensity factor with inhomogeneous coefficient and gradient parameter under the same condition of strip height, crack length, and gradient parameter . It is noticed that the greater gradient parameter and the smaller inhomogenous coefficient lead to the greater stress intensity factor.

Figure 4 shows the variation curves of the dimensionless stress intensity factor with crack length and gradient parameter under the same condition of strip height, inhomogenous coefficient, and gradient parameter . It is noticed that the greater gradient parameter and the longer crack lead to the greater of the stress intensity factor.

5. Conclusions

(1)There is singularity in crack tip of the orthotropic functionally graded material. It is the same as the orthotropic homogeneous material. Thus, we can study the complex fracture problem of the orthotropic functionally graded material in using the conventional method fracture mechanics.(2)Figure 2 shows the stress intensity factor increase with the increasing of crack length in the condition of fixed strip height . The stress intensity factor decreases with the increasing of strip height in the condition of fixed crack length . (3)Figure 3 shows the stress intensity factor decrease with the increasing of inhomogenous coefficient in the condition of fixed gradient parameter , The stress intensity factor increases with the increasing of gradient parameter in the condition of fixed inhomogenous coefficient . (4)Figure 4 shows the stress intensity factor increase with the increasing of gradient parameter in the condition of fixed crack length . The stress intensity factor increases with the increasing of crack length in the condition of fixed gradient parameter .

Acknowledgments

Project supported by the Science and Technology Major Project of the Ministry of Education of China (208022) and the Graduate Scientific and Technological Innovation Project of the Taiyuan University of Science and Technology (20122005 and 20125027).