## Abstract

In this study, the stochastic finite element solution is given to obtain the variability in the natural frequency of functionally graded material (FGM) beam. The elastic modulus is assumed to vary in the thickness direction and the width of the beam to vary as well in the longitudinal direction following the exponential law. The random material properties of elastic modulus and mass density of the FGM beam are assumed to be one-dimensional homogeneous stochastic processes. The stochastic finite element analysis of FGM beam is performed in conjunction with Monte Carlo simulation (MCS) employing the spectral representation method for 16, the description of random processes of the random material properties under consideration. The response variability of the natural frequency due to random elastic modulus is evaluated for various states of randomness. Furthermore, the investigation on the effect of the correlation between random elastic modulus and random mass density on the response variability is addressed in detail as well.

### Keywords

- functionally graded materials
- finite element method
- FGM beam
- Monte Carlo simulation

## 1. Introduction

Functionally graded materials (FGM) have received considerable attention in many engineering applications, since the theory of FGM was firstly introduced in 1984. In general, FGM is made from the volume fractions of two or more material components that have continuous variation of material properties from one surface to another [1]. Nowadays, FGM suits the specific demand in different engineering applications, especially for high temperature environment applications of heat exchanger tubes, thermal barrier coating for turbine blades, thermoelectric generators, furnace linings, electrically insulated metal ceramic joints, space/aerospace industries, automotive applications, biomedical area, etc.

The manufacturing of FGM with fully specified profile of material gradation, however, is very difficult causing significant variability in their mechanical and structural properties. Therefore, proper handling of the randomness in the material properties is required for accurate prediction of structural response for safe and reliable design. The stochastic analysis is a useful analytical tool to predict the response of structures with random material properties. In this direction, there is a reasonable body of recent research on the effect of uncertainties in material properties on the mechanical behavior of FGM. Investigators used stochastic simulation to study the effect of microstructural randomness on stress in FGM [2]. Ferrante et al. studied the effect of non-Gaussian porosity randomness on the response of functionally graded plate [3]. Yang et al. dealt with the stochastic bending response of moderately thick FGM plates [4, 5]. The effect of random material properties on post buckling response of FGM plate are presented in Lal et al. [6].

However, the above mentioned literatures are for the static analysis. To the best of author’s knowledge, few limited works have been done on the eigen analysis of FGM structures involving randomness in system parameters. Certain efforts have been made in the past to predict the dynamic behavior of structures with randomness. In most of the studies conducted, investigators dealt with the free vibration of functionally graded laminates with random material properties using first-order perturbation technique (FOPT) incorporating mixed type and semi-analytical approach to derive the standard eigenvalue problem [7]. Some of these papers presented the stochastic finite element method (SFEM) to investigate the natural frequency of functionally graded plates based on the higher-order shear deformation theory (HSDT) utilizing first-order reliability method and second-order reliability method [8]. In most cases, Jagtap et al. [9] examined the stochastic nonlinear free vibration response of FGM plate using HSDT with von-Karman kinematic nonlinearity via direct iterative stochastic finite element method. Shegokar et al. investigated the stochastic finite element nonlinear free vibration analysis of FGM beam with random material properties due to thermo-piezoelectric loadings [10]. The above mentioned literatures investigated the free vibration and nonlinear behavior of FGM beam and plate. The material properties, such as Young’s modulus, shear modulus, and Poisson’s ratio of FGM, are modeled as independent random variables.

In this chapter, the stochastic finite element solution is suggested to obtain the variability in the natural frequency of functionally graded material (FGM) beam. The elastic modulus and width of the FGM beam are assumed to vary in thickness and longitudinal directions following the exponential law. The uncertain material properties, such as modulus of elasticity and mass density of the FGM beam, are considered to be a one-dimensional homogeneous stochastic process. The stochastic finite element analysis of FGM beam is performed using the spectral representation method for the description of randomness in conjunction with Monte Carlo simulation (MCS). The response variability of natural frequency due to random elastic modulus and mass density in FGM beam is given. Furthermore, the effect of correlation between the two random parameters is observed as well.

## 2. Theory formulation of FGM

An FGM is defined to be a material which has a continuous gradation through-the-thickness (h). One side of the material is typically ceramic and the other side is typically metal. A mixture of the two materials composes the through-the-thickness characteristics. Let us consider a functionally graded beam shown in Figure 1. The parameters of the model FGM beam are as follows: *L* is the length of the beam,

The elastic material properties vary through-the-beam thickness according to the volume fractions of the constituents using power law distribution (as shown in Figure 2).

with

where *E* denotes the effective Young’s modulus of elasticity, and *Em* and *Ec* represent the Young’s modulus of metal and ceramic, respectively. The parameters *g*(*z*) and *p* represent the volume fraction of the metal and ceramic exponent, respectively.

In exponential law, for the material parameter of Young’s modulus *E*, mass density *b*, with absolute values for *z* coordinate, which endows the symmetric characteristic to the beam with respect to mid-plane.

In Eq. (3), *z =* 0) of the beam. The parameter

## 3. Finite element formulation for FGM beam and frequency analysis

### 3.1 Finite elements

In case of four degree of freedom beam element, as shown in Figure 4, the transverse displacement function may be assumed as a cubic polynomial in

The width and mass per unit length of the element are

where

The nodal displacement vector of the element is

then, the displacement field is

where

In this case, the stiffness of the beam

where

### 3.2 Application of Hamilton’s principle

Hamilton’s principle may be a theoretical base for dynamical systems by its nature of integral form in time with Lagrangian density to account for continuous space. In this paper, the analysis of natural frequency of FGM beam is performed using Hamilton’s principle. The strain energy expression

The kinematic energy

where

Substituting Eq. (7) into Eqs. (9) and (10), the following can be obtained:

The governing differential equations of motion and the related governing equation can be derived using Hamilton’s principle

where

Substituting Eqs. (13) and (14) into Eq. (15), the following can be obtained:

here,

For simple harmonic vibration, we assume the displacements to be

where

## 4. Modeling of randomness

### 4.1 Mathematical expression

In order to model the randomness in the material properties, the modulus of elasticity and mass density along the mid-plane are assumed to vary along its length of FGM beam in a random manner. We can model these variations as one-dimensional univariate (1D-1V) homogeneous stochastic processes. The simple mathematical expressions for the randomly varying modulus of elasticity and mass density can be written as

where

The numerical generation of sample functions of Gaussian zero-mean homogeneous stochastic processes, which describe the randomness in parameters of the structure, is accomplished using the spectral representation method. For a one-dimensional univariate (1D-1V) stochastic process, we have [11]

In Eq. (21),

where

The uniform random phase angle

here,

In all examples, the coefficient of variation

here,

and

### 4.2 Monte Carlo analysis

In order to obtain the response variability in the natural frequency of the FGM beam, we employed the scheme of Monte Carlo simulation (MCS). As a matter of fact, the MCS corresponds to the deterministic analyses on a set of heterogeneous models of the given structure, in which the material properties have different values depending on the position in the domain of the structure.

The generation of heterogeneous random samples is accomplished by the aforementioned spectral representation scheme, and we use 10,000 samples for respective analyses. In particular, we adopt the local average scheme other than the mid-point rule in applying the MCS, with which better results can be obtained especially for the processes with small correlation distance. Figure 5 shows an example plot of the processes employed to model the randomness in the system parameters.

## 5. Numerical example

The geometric dimensions of the example FGM beam are: ^{3}, and

### 5.1 Deterministic analysis results

The results in Figure 7 correspond to the prismatic homogeneous beam since the parameters in exponents,

The discrepancies between exact and finite element solutions for the frequencies for the first three modes are shown in Figure 8. The differences given in percentile tends to zero as the number of finite elements is increased, meaning the results are converging to exact solutions.

Figure 9 shows the first three normalized natural frequencies of uniform FGM beams for three cases of modulus ratio (

### 5.2 Variability of natural frequency due to randomness in elastic modulus

Figure 10 shows the COV of response versus the correlation distance (

The relationship between COV of natural frequency and the COV of stochastic process is shown in Figure 11. The standard deviation of stochastic process is changed from 0.0 to 0.25. As seen in Figure 11, the COV of response shows a slightly nonlinear pattern in all the cases of

The effect of mesh refinement on the COV of natural frequency is shown in Figure 12. The correlation distance log (

Figure 13 shows the variation of the coefficient of variation (COV) depending on the non-uniformity parameter (

The overall features of the effect of Young’s modulus ratio on COV of natural frequency are shown in Figure 14. The COV of natural frequency is not affected by the parameter

### 5.3 FGM beam having correlation multiple randomness

It is natural to have preposition that not only the elastic modulus, but also the mass density of the material can have randomness. Therefore, we need to consider the effect of correlation between two random parameters of elastic modulus and mass density. To this aim, we consider three correlation cases: +1.0, 0.0, and −1.0. When the random processes for elastic modulus,

As shown in Figure 15, the maximum COV exceeds the input standard deviation of the stochastic process in the case of negative perfect correlation, while it is about 75% in the case of no correlation. However, in case of positive perfect correlation (Figure 15c), the response COV of natural frequency of FGM beam is small enough to be ignored. These results can easily be understood because the stiffness and mass matrix of each elements in Eq. (18) increase or decrease with the same rate in the case of positive perfect correlation. In case of negative perfect correlation, the ratio of random parts in the stiffness and mass are relatively large since the random parts have opposite sign, which makes the response variability large.

In the case, when we take the correlated multiple random material properties into account, we also obtained the slight nonlinear pattern of variation of COV in terms of COV of stochastic fields, as shown in Figure 16.

## 6. Conclusions

To evaluate response variability due to a single parameter of the random Young’s modulus and multiple uncertain material properties, a formulation in the context of stochastic finite element solution is suggested for the natural frequency of FGM beam. In deriving the formula for the covariance of the response, modified power spectral density and correlation function are defined by using the general formula of random processes. The Monte Carlo simulation is performed employing the statistical preconditioning scheme as a random process generation technique. The local average method is employed instead of mid-point rule in Monte Carlo simulation.

In FGM beam natural frequency, the response COV for correlation between a random of single parameter and two uncertain parameters is observed. The coefficient of variation of natural frequency can only reach about 50% of the input standard deviation of the stochastic process in a single parameter of the random elastic modulus. However, the number of values is increased over 100% of the input standard deviation of the stochastic process in multiple uncertain material properties, when the correlation distance tends to infinity. The results showed that the COV of natural frequency of FGM beam in a single parameter of the random Young’s modulus and multiple uncertain material properties achieved maximum variability for *d* about 1.0.

There is a very small difference between deterministic natural frequency and probabilistic natural frequency of FGM beam for the case of positive perfect correlation. Also, the COV of natural frequency does not depend on the number of elements, Young’s modulus ratio, and the ratio of non-uniformity parameter of FGM beam. The importance of these parameters needs to be studied as a further work.

Hermite shape functions of beam finite element:

Stiffness matrix:

Mass matrix: