Shock and Vibration
Volume 2016 (2016), Article ID 3287645, 10 pages
http://dx.doi.org/10.1155/2016/3287645

Research Article

Fundamental Frequency Analysis of Sandwich Beams with Functionally Graded Face and Metallic Foam Core

Lin Mu and Guiping Zhao

State Key Laboratory of Mechanical Structure Strength and Vibration, School of Aerospace, Xi’an Jiaotong University, Xi’an 710049, China

Received 5 August 2015; Revised 7 November 2015; Accepted 11 November 2015

Academic Editor: Jussi Sopanen

Copyright © 2016 Lin Mu and Guiping Zhao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This study is interested in assessing a way to analyze fundamental frequency of sandwich beams with functionally graded face sheet and homogeneous core. The face sheet, which is an exponentially graded material (EGM) varying smoothly in the thickness direction only, is composed of a mixture of metal and ceramic. The core which is made of foam metal is homogeneous. The classical plate theory (CPT) is used to analyze the face sheet and a higher-order theory (HOT) is used to analyze the core of sandwich beams, in which both the transverse normal and shear strains of the core are considered. The extended Galerkin method is used to solve the governing equations to obtain the vibration equations of the sandwich beams suitable for numerical analysis. The fundamental frequency obtained by the theoretical model is validated by using the finite element code ABAQUS and comparison with earlier works. The influences of material and geometric properties on the fundamental frequency of the sandwich beams are analyzed.

1. Introduction

Functionally graded materials (FGM), consisting of two or more constituent phases with continuous and smoothly varying composition, are one of advanced materials in the family of engineering composites [1]. FGM refers to a heterogeneous composite material with gradient variation from one surface of the material to the other which results in continuously varying material properties [2]. Compared with fiber-reinforced composite laminates, FGM can eliminate the stress concentration at the interface of the layers and overcome the delaminating in composite laminates [3]. Functionally gradient ceramic materials are usually made from a mixture of metal and ceramics. In the mixture, the ceramic can resist high temperature in thermal environments, while the metal can decrease the tensile stress on the ceramic surface at the earlier state of cooling. These FGM are widely used in great application potential in many engineering sectors such as aerospace vehicles and power generators.

Many researchers have investigated the static and the dynamic characteristics of functionally graded sandwich structures. Zenkour [4, 5] presented a two-dimensional solution for bending analysis of simply supported functionally graded ceramic-metal sandwich plate and studied the bucking and free vibration of the simply supported functionally graded sandwich plate. In his study, the faces of sandwich plate are made of functionally graded material, and the core layer is made of the isotropic ceramic material. Woodward and Kashtalyan [6] presented a 3D elasticity solution for bending response of sandwich plate, in which the stiffness of the core varies gradually in the thickness direction. The use of graded core to improve performance of sandwich structures, especially under localized loading, is examined and discussed. Zenkour and Alghamdi [7] studied the bending issue of sandwich plate with functionally graded faces and homogeneous ceramic core under mechanical and thermal loads. Cinefra and Soave [8] presented the closed form solutions for free vibration of simply supported sandwich plate with isotropic faces and functionally graded core, which are composed of two different materials. Kashtalyan and Menshykova [9] presented a three-dimensional elasticity solution for sandwich plate with a functionally graded core subjected to transverse loading. Their studies indicated that the use of graded core instead of a conventional homogeneous one can eliminate discontinuity of the in-plane normal and shear stresses across the adjoining layers (between the face sheets and the core). Based on the three-dimensional elastic theory, Li et al. [10] studied the free vibration of functionally graded sandwich plates with simply supported and clamped edges. Neves et al. [11, 12] analyzed the free vibration of sandwich plate/shell with gradient core using the Carrera Unified Formulation (CUF) [13, 14] and meshless method [15]. Thai et al. [16] presented a new first-order shear deformation theory for functionally graded sandwich plates composed of functionally graded face sheets and an isotropic homogeneous core. And analytical solutions for bending, buckling, and free vibration analysis of rectangular plates under various boundary conditions are presented. Based on a refined shear deformation theory, Vo et al. [17] presented a finite element model for vibration and buckling of functionally graded sandwich beams, which is made of metallic or ceramic core and functional graded face sheets.

Although many researchers have investigated the mechanical properties of the functionally graded sandwich structures, the vibration characteristics of the sandwich structures with graded face sheets and metallic foam cores have been of little concern up to now. And the sandwich structures with graded face sheets and metallic foam core were widely used for the leading edges of the swing on airplane or submarine hull, which require high specific stiffness and specific strength.

In this study, we develop a way to analyze composite sandwich beams consisting of graded face sheets and metallic foam core, which is very different from metallic or ceramic core that has been analyzed in literature. The face sheets of the sandwich beams are a nonhomogeneous composite which feature a gradual change in properties with the continuous change of component (ceramic-metal) from one surface to the other. Young’s modulus and density are characterized by a power law distribution. The core of the sandwich beam is isotropic aluminum foams. The fundamental frequency of the sandwich beams is calculated and influences of material and geometric properties on the fundamental frequency of the sandwich beams are analyzed.

2. Mechanical Model

The mechanical model is described in Figure 1. In this model, is the thickness of face sheet and measures the thickness of the core. is the length of the sandwich beam and is the thickness of the sandwich beam. In the local Cartesian coordinate system , - and -axis are set to the centerline of global mid-surface of the beam and thickness direction, respectively. So the global mid-surface is the beam reference surface.

Figure 1: Functionally graded sandwich beam.

The basic assumptions of the model are as follows:

(1)The face sheets are composed of functionally graded materials through thickness direction, Young’s modulus and material densities vary continuously in thickness direction, and Poisson’s ratio is assumed constant.(2)The face sheets are very thin so that they can carry the tangential and bending loads and follow Kirchhoff hypothesis.(3)The core is made of isotropic soft material and capable of carrying transverse normal and shear stresses.(4)The deflection in transverse direction is large, while the tangential deformation of beam is small as the tangential stiffness of the sandwich beam is large. In the process of deformation, the face sheets are considered to be incompressible in the transverse direction while the core, due to its soft material, is considered to be compressible.

2.1. Displacements of the Model

The displacements of the sandwich beam can be properly described by and in and directions of the top face sheet, respectively:

where and are the mid-surface displacements in and directions of the top face sheet, respectively. By analogy, and are the displacements in and directions of the bottom face sheet, respectively:

where and are the mid-surface displacements in and directions of the bottom face sheet, respectively.

For the core, and represent the displacements in and directions of the core, respectively. The tangential displacement is expanded in a second-order power series to the thickness direction; only a first-order term is used for the transverse displacement [18]. Consider

where is an additional warping function describing the deformation of the sandwich core. and are defined as

2.2. Functionally Graded Materials

The face sheets of the sandwich beam are gradient composite made of ceramic (AlO₂) and metal (Al). Young’s modulus and mass density vary continuously in the thickness direction (-axis direction).

The effective material properties can be expressed as

where , are Young’s modulus for the ceramic (AlO₂) and metal (Al), and

where , are the mass density for the ceramic (AlO₂) and metal (Al) and is the gradient exponent and .

2.3. Lagrangian of the Deformed System

The governing equations of sandwich beams are derived by Hamilton’s principle of the minimization of the Lagrangian of the deformed system:

where , , and denote the variations of the kinetic energy, strain energy, and external work performed during arbitrary instants of time , respectively.

The variation of the kinetic energy and strain energy can be written as

where is the mass density function of the face sheet, is the mass density of the core, and , and denote the speeds of the top face sheet, bottom face sheet, and the core, respectively. In addition, it is very direct to express the variation of the work performed by external loads as

where and are the uniformly distributed pressure on the top and bottom face sheets, respectively.

Substituting ,  , and into Hamilton’s principle in (7) and considering stress integrals through the thickness of the beam, the governing equations are obtained by setting the coefficients of , , , , and over the field of integration to zero separately:

where masses are given in Appendix A. The force resultants are given by

From (11a), (11b), (11c), (11d), (11e), and (11f) we determined

where , , and are given in Appendix B.

After substituting (12a), (12b), (12c), and (12d) into (10a), (10b), (10c), (10d), and (10e) we arrived at

2.4. Extended Galerkin Method

Now , , , and are cast into state-space form:

At the same time, the boundary conditions of the sandwich beams are considered to be clamped at their two ends and ; then

with the initial conditions

We could make the boundary conditions satisfied by taking

where represents the th normalized eigenmode of a clamped-clamped Euler-Bernoulli beam, defined as [19]

where − and is the positive th root of the equation .

Solving (13a), (13b), (13c), and (13d) is tackled by using the extended Galerkin method and then the mass matrix and stiffness matrix for the numerical analysis are determined as a matrix equation:

The entities in mass matrix and stiffness matrix are listed in Appendix C.

3. Validation of Gradient Sandwich Beam Model

To assess the validity of the model of the sandwich beam with gradient face sheet, we calculate the natural frequencies and compare the results with those obtained by using ABAQUS code. The face sheet is divided into ten layers and the thickness of each layer is the same (see Figure 2(a)). Each layer of the face sheet model has 100 elements and the core model has 1000 elements. The sizes of meshes are 10 mm × 1 mm for the face sheet and 10 mm × 3 mm for the core. The CPE8R quadratic plane-strain elements are used for both core and face sheets. The nodes of the core and face sheets at their common interfaces are tied to ensure continuity of the displacements. Ends of the face sheets and core are taken to be clamped in the mode. The mesh of the sandwich beam is shown in Figure 2(b).

(a) The face sheets with ten layers

(b) FE model and mesh

Figure 2: Model of the sandwich beam.

The material properties of the composite face sheets and the core are as follows.

The face sheets [3]:

The core [20]:

Geometric parameters of the sandwich beam:

The first four frequencies of sandwich beam calculated by the present analytical model and numerical simulation are presented in Table 1. means the th natural frequency. For example, “” means the second-order natural frequency. It can be found that the natural frequencies of the gradient sandwich beam predicted by the analytical model are very close to those simulated by ABAQUS.

Table 1: The first four frequencies (rad/s) with increasing Ns (comparison between the present result and that of ABAQUS).

Ns
th natural frequency (rad/s)

2
2416.0 (3.46%)
6050.6 (20.12%)
—
—

3
2389.5 (2.33%)
5248.6 (4.20%)
9200.2 (12.20%)
—

4
2371.8 (1.57%)
5185.0 (2.94%)
8471.7 (3.32%)
12544.5 (8.64%)

5
2358.4 (0.99%)
5152.6 (2.29%)
8403.1 (2.48%)
11832.7 (2.47%)

6
2354.7 (0.84%)
5119.5 (1.64%)
8365.6 (2.02%)
11773.4 (1.96%)

7
2348.5 (0.57%)
5110.7 (1.46%)
8328.5 (1.57%)
11732.7 (1.61%)

8
2347.0 (0.51%)
5094.2 (1.13%)
8316.4 (1.42%)
11700.2 (1.33%)

9
2343.9 (0.37%)
5090.0 (1.05%)
8297.1 (1.19%)
11685.6 (1.20%)

By ABAQUS
2335.2
5037.1
8199.7
11547.0

Note. Relative error, ([prediction by the present model] − [result by ABAQUS])/[result by ABAQUS] × 100%.

In order to further verify the present model, the fundamental frequencies of sandwich beams with different values of thickness ratio in comparison with earlier works [17] are given in Tables 2 and 3. means the gradient exponent and means the nondimensional fundamental natural frequencies. It can be seen that the frequencies predicted by the present model do have an excellent agreement with these from Vo and Thai’s model. Therefore, it can be concluded that the present analytical model is reliable and applicable.

Table 2: Comparison of the nondimensional fundamental natural frequencies () with various thickness ratios of the sandwich beam ().

1-1-1
1-2-1
1-8-1

Ref. []
Present
Ref. []
Present
Ref. []
Present

0.5
7.3914
7.5552 (2.22%)
7.0539
7.0856 (0.45%)
6.3333
5.8824 (−7.12%)

1
7.9726
8.0061 (0.42%)
7.5187
7.5140 (−0.06%)
6.6705
6.3057 (−5.47%)

2
8.4655
8.3244 (−1.67%)
7.8960
7.8126 (−1.06%)
6.9233
6.6750 (−3.59%)

5
8.9201
8.5768 (−3.85%)
8.2498
8.0356 (−2.60%)
7.1155
6.9648 (−2.12%)

10
9.1193
8.6898 (−4.71%)
8.4162
8.1250 (−3.46%)
7.1884
7.0768 (−1.55%)

Note. Ref. [] ≡ Vo et al. [17], and relative error: [(prediction by the present model) − (Ref. [])]/(Ref. []) × 100%. .

Table 3: Comparison of the nondimensional fundamental natural frequencies () with various thickness ratios of the sandwich beam ().

1-1-1
1-2-1
1-8-1

Ref. []
Present
Ref. []
Present
Ref. []
Present

0.5
10.2298
10.1771 (−0.52%)
9.7587
9.7449 (−0.14%)
8.2818
7.9373 (−4.16%)

1
11.2767
11.0023 (−2.43%)
10.7706
10.6093 (−1.50%)
9.0064
8.6446 (−4.02%)

2
11.9911
11.7095 (−2.35%)
11.5168
11.3559 (−1.40%)
9.6175
9.2733 (−3.58%)

5
12.4199
12.2989 (−0.97%)
12.0423
11.9846 (−0.48%)
10.1297
9.8627 (−2.64%)

10
12.5230
12.5740 (0.41%)
12.2122
12.2204 (0.07%)
10.3315
10.0985 (−2.26%)

Note. Ref. [] ≡ Vo et al. [17], and relative error: [(prediction by the present model) − (Ref. [])]/(Ref. []) × 100%. .

4. Result and Discussion

In order to investigate the influences of material and geometric properties on the fundamental frequency of the sandwich beams, a sandwich beam of length and thickness is considered. The sandwich beam is composed of functionally graded face sheets and aluminum foam core. And boundary conditions are clamped. According to the researches by Ashby [21], the relationship of Young’s modulus and density of closed-cell aluminum foam is as follows:

where and are Young’s modulus and mass density of closed-cell aluminum foam and and are Young’s modulus and mass density of matrix material of the closed-cell aluminum foam and is a coefficient chosen by experience. Simply for ease of accounting, is taken to be 1.0.

Figure 3 shows fundamental frequency of the sandwich beams with the same thickness and different thickness ratio between face sheets and core. In Figure 3, legends 1-8-1, 1-7-1, 1-6-1, 1-5-1, 1-4-1, and 1-3-1 mean relative thickness of the sandwich beam following by top sheet, core, and bottom sheet. For example, 1-8-1 means the thickness of the core is eight times the thickness of face sheet.

(a)

(b)

(c)

Figure 3: Fundamental frequency of the sandwich beams with different thickness ratios () versus density of the core for different gradient exponents. (a) ; (b) ; (c) .

When the gradient exponent is a constant, the fundamental frequency of the sandwich beam increases with the increasing density of the core firstly. And then the fundamental frequency arrives at a peak value. After that the fundamental frequency decreases slowly with the increasing density of the core. An interesting phenomenon is that the fundamental frequency is almost the same value (named the special frequency, for short) when the density of the core is a specific value ( for short) for all different thickness ratios of sandwich beams from 1-8-1 to 1-3-1. The reason is that the mass and stiffness of sandwich beams change and influence each other at the same time with changing thickness ratio. Figure 3 also shows that the special frequency () presents an increasing trend when the gradient exponent changes from 0.1 to 10.

Figure 4 shows that the fundamental frequency increases rapidly with increasing of the gradient exponent as is smaller than 1 for different densities of foam cores. The growth of the fundamental frequency slows down and it gradually reaches a constant value after the gradient exponent is larger than 1. The reason is that the mass and the stiffness of the gradient face sheets do not change obviously after the gradient exponent is larger than 1.

(a)

(b)

(c)

Figure 4: Fundamental frequency of the sandwich beams with different thickness ratios versus gradient exponent for different densities of foam core. (a) ; (b) ; (c) .

5. Conclusions

(1)An analytical model to analyze fundamental frequency of sandwich beams with functionally graded face sheets and metallic foam core is developed.(2)When the thickness of the sandwich beam is a constant and the thickness ratio of each part from bottom to top of the sandwich beam is varying, the fundamental frequency is very close to a fixed value ( for short) when the relative density of the core is taken to a specific value ( for short). This result can be used to optimize the design of the sandwich beams.(3)The influences of the gradient exponent on the fundamental frequency of sandwich beams are very obvious. When the gradient exponent is larger than 1, the influence is slight.

Appendices

A. Stresses in Top and Bottom Face Sheets and Core

Consider

B. Stiffness Quantities

Consider

C. Elements in Mass Matrix and Stiffness Matrix

Consider

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (11372237).

References

M. Koizumi, “The concept of FGM,” Ceramic Transactions of Functional Graded Materials, vol. 34, pp. 3–10, 1993.

D. K. Jha, T. Kant, and R. K. Singh, “Free vibration response of functionally graded thick plates with shear and normal deformations effects,” Composite Structures, vol. 96, pp. 799–823, 2013.

S. Xiang, G.-W. Kang, M.-S. Yang, and Y. Zhao, “Natural frequencies of sandwich plate with functionally graded face and homogeneous core,” Composite Structures, vol. 96, pp. 226–231, 2013.

A. M. Zenkour, “A comprehensive analysis of functionally graded sandwich plates: part 1-deflection and stresses,” International Journal of Solids and Structures, vol. 42, no. 18-19, pp. 5224–5242, 2005.

A. M. Zenkour, “A comprehensive analysis of functionally graded sandwich plates. Part 2. Buckling and free vibration,” International Journal of Solids and Structures, vol. 42, no. 18-19, pp. 5243–5258, 2005.

B. Woodward and M. Kashtalyan, “Bending response of sandwich panels with graded core: 3D elasticity analysis,” Mechanics of Advanced Materials and Structures, vol. 17, no. 8, pp. 586–594, 2010.

A. M. Zenkour and N. A. Alghamdi, “Bending analysis of functionally graded sandwich plates under the effect of mechanical and thermal loads,” Mechanics of Advanced Materials and Structures, vol. 17, no. 6, pp. 419–432, 2010.

M. Cinefra and M. Soave, “Accurate vibration analysis of multilayered plates made of functionally graded materials,” Mechanics of Advanced Materials and Structures, vol. 18, no. 1, pp. 3–13, 2011.

M. Kashtalyan and M. Menshykova, “Three-dimensional elasticity solution for sandwich panels with a functionally graded core,” Composite Structures, vol. 87, no. 1, pp. 36–43, 2009.

Q. Li, V. P. Iu, and K. P. Kou, “Three-dimensional vibration analysis of functionally graded material sandwich plates,” Journal of Sound and Vibration, vol. 311, no. 1-2, pp. 498–515, 2008.

A. M. A. Neves, A. J. M. Ferreira, E. Carrera et al., “Static, free vibration and buckling analysis of isotropic and sandwich functionally graded plates using a quasi-3D higher-order shear deformation theory and a meshless technique,” Composites Part B: Engineering, vol. 44, no. 1, pp. 657–674, 2013.

A. M. A. Neves, A. J. M. Ferreira, E. Carrera et al., “Free vibration analysis of functionally graded shells by a higher-order shear deformation theory and radial basis functions collocation, accounting for through-the-thickness deformations,” European Journal of Mechanics A: Solids, vol. 37, pp. 24–34, 2013.

E. Carrera, “Evaluation of layer-wise mixed theories for laminated plates analysis,” AIAA Journal, vol. 36, no. 5, pp. 830–839, 1998.

E. Carrera, S. Brischetto, M. Cinefra, and M. Soave, “Effects of thickness stretching in functionally graded plates and shells,” Composites Part B: Engineering, vol. 42, no. 2, pp. 123–133, 2011.

A. J. M. Ferreira, C. M. C. Roque, and P. A. L. S. Martins, “Analysis of composite plates using higher-order shear deformation theory and a finite point formulation based on the multiquadric radial basis function method,” Composites Part B: Engineering, vol. 34, no. 7, pp. 627–636, 2003.

H.-T. Thai, T.-K. Nguyen, T. P. Vo, and J. Lee, “Analysis of functionally graded sandwich plates using a new first-order shear deformation theory,” European Journal of Mechanics—A/Solids, vol. 45, pp. 211–225, 2014.

T. P. Vo, H.-T. Thai, T.-K. Nguyen, A. Maheri, and J. Lee, “Finite element model for vibration and buckling of functionally graded sandwich beams based on a refined shear deformation theory,” Engineering Structures, vol. 64, pp. 12–22, 2014.

J. Hohe and L. Librescu, “Recent results on the effect of the transverse core compressibility on the static and dynamic response of sandwich structures,” Composites Part B: Engineering, vol. 39, no. 1, pp. 108–119, 2008.

L. Mu, B. He, D. Xiao, G. Zhao, and C. Cho, “Dynamic response of beam-like composite sandwich plates with transversely flexible cores,” Ocean Engineering, vol. 108, pp. 778–788, 2015.

Z. Qin and R. C. Batra, “Local slamming impact of sandwich composite hulls,” International Journal of Solids and Structures, vol. 46, no. 10, pp. 2011–2035, 2009.

M. F. Ashby, Metal Foams: A Design Guide, Butterworth-Heinemann, Boston, Mass, USA, 2000.