- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents
Advances in Mechanical Engineering
Volume 2013 (2013), Article ID 634584, 10 pages
Free Vibration Analyses of FGM Thin Plates by Isogeometric Analysis Based on Classical Plate Theory and Physical Neutral Surface
Department of Engineering Mechanics, Hohai University, Nanjing 210098, China
Received 2 January 2013; Revised 23 May 2013; Accepted 6 June 2013
Academic Editor: Indra Vir Singh
Copyright © 2013 Shuohui Yin et al. 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.
The isogeometric analysis with nonuniform rational B-spline (NURBS) based on the classical plate theory (CPT) is developed for free vibration analyses of functionally graded material (FGM) thin plates. The objective of this work is to provide an efficient and accurate numerical simulation approach for the nonhomogeneous thin plates and shells. Higher order basis functions can be easily obtained in IGA, thus the formulation of CPT based on the IGA can be simplified. For the FGM thin plates, material property gradient in the thickness direction is unsymmetrical about the midplane, so effects of midplane displacements cannot be ignored, whereas the CPT neglects midplane displacements. To eliminate the effects of midplane displacements without introducing new unknown variables, the physical neutral surface is introduced into the CPT. The approximation of the deflection field and the geometric description are performed by using the NURBS basis functions. Compared with the first-order shear deformation theory, the present method has lower memory consumption and higher efficiency. Several numerical results show that the present method yields highly accurate solutions.
Although numerous plate theories have been developed since the late 19th century, the classical plate theory (CPT)  and the first-order shear deformation theory (FSDT)  are the most widely accepted and applied theories in engineering.
The CPT is the simplest theory with high computational efficiency that neglects transverse shear strains and midplane displacements. The effects of transverse shear strains and midplane displacements in homogeneous thin plates are small; hence, the CPT is widely used in homogeneous thin plates . However, nonhomogeneous thin plates such as functionally graded material (FGM) plates have material property gradient in the thickness direction and this gradient is not symmetric about the midplane, so the results from the CPT may be incorrect. Our studies show that the natural frequencies of FGM thin plates based on the CPT and FSDT are close to each other for small or big gradient indices, while a significant difference exists for moderate gradient indices. As we all know, the transverse shear deformations have small effects in thin plates, so neglecting midplane displacements in the CPT causes differences in the FSDT solution for FGM thin plates. The material appears relatively homogeneous for small or big gradient indices; in other words, the material property gradient in the thickness direction is almost symmetric about the midplane, thus midplane displacements can be neglected, whereas the material is obviously nonhomogeneous for moderate gradient indices; this means that the material property gradient in the thickness direction is apparently unsymmetrical about the midplane, so neglecting midplane displacements in the CPT causes significant differences in the FSDT solution for FGM thin plates with moderate gradient indices. Abrate  indicated that stretching-bending couplings in constitutive equations of FGM plates do not exist if the proper reference surface is selected. Zhang and Zhou  recently proposed that the physical neutral surface is the proper reference surface and presented a theoretical analysis for FGM thin plates based on the physical neutral surface. To eliminate the effects of midplane displacements without introducing new unknown variables, we integrated the physical neutral surface into the CPT, that is, CPT_neu.
Regarding FGM plates of general geometries and loading conditions, analytical methods are usually infeasible, and numerical approaches are thus required. The finite element method is considered as one of the most versatile and powerful numerical methods. However, element-based polynomial approximation in the standard finite element method inevitably leads to discretization errors in complex geometry; the gap between the finite element analysis (FEA) and computer-aided design (CAD) consumes much time on the mesh generation of FEA. In 2005, Hughes et al.  developed a new numerical method called isogeometric analysis (IGA) to achieve a seamless integration between CAD and FEA. The principle of the IGA is that spline (such as NURBS and T-splines) basis functions are employed as shape functions for geometric description and field approximation. The unique advantages of the IGA are as follows : (1) the same basis function space is adopted for geometric and mesh models, thus avoiding mesh discretization errors even on coarse mesh; (2) the natural division in the spline parametric domain is used instead of the mesh, thereby avoiding the traditional mesh generation procedure; (3) mesh refinement does not necessitate an external description of the geometry, thereby simplifying the adaptive mesh refinement; (4) the splines have higher order continuity, thus simplifying the construction of global higher order continuous conforming elements. In previous years, considerable work has been conducted in improving or applying the original IGA in some fields, including structural vibrations , plates and shells [8–10], fluid mechanics , fluid-structure interaction problems [12, 13], damage and fracture mechanics , and structural shape optimization .
The CPT needs higher order elements, so it is difficult to achieve it in tradition finite element method. In the IGA, higher order NURBS basis functions can be easily obtained, so the formulation of CPT can be simplified [1, 3]. The main objective of this research work is to extend CPT-based IGA to study free vibration of FGM thin plates. To serve this purpose, the CPT_neu is introduced into the NURBS-based IGA thus the method can analyze nonhomogeneous thin plates. Compared with the FSDT-based IGA, the present method has lower memory consumption and higher efficiency.
The paper is organized as follows. The NURBS-based IGA is described briefly in Section 2. In Section 3, the effects of midplane displacements on the natural frequency of FGM thin plates are investigated in detail, the physical neutral surface is then introduced into the CPT, and the CPT_neu-based eigenvalue equations of the IGA are developed. The numerical validation is presented in Section 4. The method is illustrated on several numerical examples in Section 5. Finally, conclusions are drawn in Section 6.
2. NURBS-Based Isogeometric Analysis
2.1. NURBS Basis Functions
NURBS basis functions are linear combinations of B-spline basis functions . B-spline basis functions are used to parameterize curves and are generally defined at the knot vector of the parametric domain . To formulate B-spline basis functions of order , is defined as
The knot vector is called an “open-knot vector” when the end knots are repeated times. Equation (1) reflects a one-to-one relationship between and the set of B-spline basis functions. With the given knot vector , the B-spline basis function is constructed by 
If is the weight associated with the B-spline basis function , then the NURBS basis function , which is a weighted average of the B-spline basis functions, is defined as
The NURBS basis functions satisfy the partition of unity, that is,
Similarly, the two- or three-dimensional NURBS basis functions can be constructed by taking the tensor product of two or three one-dimensional B-spline basis functions: where and denote the two- and three-dimensional geometry related weights, respectively; , , and are the B-spline basis functions of order in the direction, order in the direction, and order in the direction, respectively; both and follow the recursive formula shown in (2) with knot vectors and defined as
By using the NURBS basis functions, a NURBS curve of order can be constructed as where is the control point coordinate.
In a similar method, a two-dimensional NURBS surface and a three-dimensional NURBS solid can be constructed as where and , respectively, represent the coordinates of control point in two and three dimensions.
2.2. Isoparametric Discretization
Similar to the isoparametric finite element method, mapping from the parametric domain to the physical domain in the IGA is expressed as where is the parametric coordinate, that is, in one dimension, in two dimensions, and in three dimensions; and denote the shape function and coordinates at control point ; is the number of control points.
Higher order basis functions can be easily obtained in the IGA, so only the deflection of plate is selected as the independent variable in the CPT based on the IGA, while the other two displacement components and can be obtained through . Similar to the geometry discretization of (9), the deflection field is approximated as where is the deflection at control point .
3. Fundamental Equations of FGM Thin Plate
3.1. FGM Plate
Consider a metal/ceramic FGM plate of thickness . The metal is located at the bottom face of the plate, and the ceramic is located at the top face of the plate. The material is assumed to be isotropic, and the material properties are assumed to vary smoothly and continuously from the ceramic to the metal through the thickness. The -plane is the midplane of the plate with the -axis in a positive upward from the midplane. Young’s modulus and mass density are considered inhomogeneous material properties, whereas Poisson’s ratio is assumed to be constant. The inhomogeneous material properties are expressed in power law form : where is the volume fraction exponent or the gradient index; is the thickness coordinate variable; ; subscripts and represent the ceramic and metal constituents, respectively.
3.2. Effects of Midplane Displacements on FGM Thin Plate
The CPT based on the Kirchhoff hypothesis neglects transverse shear deformations and midplane displacements but can produce reasonably accurate results for homogeneous thin plates. To consider the effect of transverse shear deformation, shear deformation theories were proposed [2, 18]. For homogeneous thin plates, the accuracy of the CPT is almost equal to that of the shear deformation theories. Some studies show that the shear deformation theories can analyse FGM plates . In the following section, we investigate whether the CPT can also effectively analyse FGM thin plates.
A simply supported aluminium/alumina (Al/Al2O3) square plate with length m is considered. The plate is modelled with cubic NURBS elements. Length-to-thickness ratios of and 100 are considered. A mesh with linear parameterization and control points is used as shown in Figure 1. The material properties of Al and Al2O3 are , , and and , , and , respectively. Figure 2 presents a comparison of the first normalized free vibration frequency of the plate obtained by the NURBS-based IGA based on the CPT and FSDT. The solutions based on the CPT and the FSDT are close to each other for small or big gradient indices. A significant difference exists between the two theories for moderate gradient indices. The length-to-thickness ratio has insignificant effects on the difference between the CPT and the FSDT for the FGM thin plate, and this means that transverse shear deformations can be neglected in thin plate. The material appears relatively homogeneous for small or big gradient indices, whereas the material is obviously nonhomogeneous for moderate gradient indices. Differences exist between the CPT and the FSDT because the FSDT considers the effects of transverse shear deformations and midplane displacements. The previous analysis shows that neglecting the midplane displacements caused by inhomogeneous materials in the CPT may cause significant differences in the FSDT solution for moderate gradient indices.
3.3. FGM Thin Plate Equations Based on the Physical Neutral Surface
Considering a FGM rectangular plate under pure bending with the theory of elasticity, the coordinate system and the material parameter variation are the same as those in Section 3.1. One surface () whose normal strains and stresses are zero can be found, and the surface is called as physical neutral surface . The definition of physical neutral surface is expressed as 
From (12), the physical neutral surface and geometric middle surface are the same for homogeneous isotropic plates.
Using the physical neutral surface instead of geometry neutral surface in the CPT, the effects of midplane displacements can be avoided. To eliminate the effects of the midplane displacements in the FGM thin plate without introducing new unknown variables, the physical neutral surface concept is integrated into the CPT, that is, CPT_neu. By using the physical neutral surface and the CPT, the displacement fields can be assumed as the following:
The strains are defines as
According to the generalized Hooke’s law, the stresses can be obtained as where
The pseudostrains and pseudostresses of the plate are denoted as follows: where
For the free vibration analyses, the Galerkin weak form of the elastodynamic undamped equilibrium equation can be expressed as follows :
Substituting (10), (17) into (19) introduces the dynamic discrete system of equations: where denotes the vector of deflections at the control points and and are the global mass matrix and stiffness matrix, respectively.
The element contribution to is and the element contribution to is with
The th normal mode can be obtained by separating the variables: where is the th natural frequency.
4. Numerical Validation
To validate the accuracy of the proposed method, the Al/Al2O3 thin plate considered in the previous section is again employed. All conditions are the same as those in the previous section. Figure 3 shows the comparison of the first normalized natural frequencies of the plate obtained by the proposed method and the NURBS-based IGA based on the FSDT. For the length-to-thickness ratios of and 100, the results obtained by the NURBS-based IGA based on the CPT_neu match well with those from the NURBS-based IGA based on the FSDT. It is concluded that the effects of midplane displacements vanish in the CPT_neu.
5. Numerical Results
In this section, several examples with reference solutions are provided to assess the robustness of the proposed method. All examples reported later are calculated by using the cubic order NURBS basis functions with linear parameterization . A Gauss quadrature scheme is used in each NURBS element. The same Al/Al2O3 material considered in the previous section is employed.
5.1. A Simply Supported Al/Al2O3 Square Thin Plate with Length m
The length-to-thickness ratios of and 100 are considered. The natural frequencies of the plates are investigated, and the results are compared with the NURBS-based IGA based on the FSDT or analytical solutions [17, 21]. Figure 4 shows the comparison of the first normalized natural frequencies from the proposed method and the reference solutions with different control points. The following features are evident: (1) results based on the CPT_neu and FSDT match well with the analytical solutions [17, 21]; (2) as the CPT, natural frequencies obtained by the CPT_neu are larger than the analytical solutions; (3) the accuracy of the CPT_neu is higher than that of the FSDT when the number of control points is small, whereas the accuracy of the CPT_neu is almost the same as that of the FSDT with increasing number of control points.
The computational time versus the number of control points is plotted in Figure 5. The results indicate that the CPT_neu is much more efficient than the FSDT because five degrees of freedom exist in one point for the FSDT, whereas only one degree of freedom exists in one point for the CPT_neu.
5.2. An Al/Al2O3 Rectangular Thin Plate with a Length and Width under Different Boundary Conditions
The effects of the gradient indices, length-to-width ratios (), and boundary conditions of the FGM thin plates on the natural frequencies are investigated to demonstrate the validity of the present method. For this purpose, the normalized natural frequencies of two length-to-width ratio plates with four gradient indices and three boundary conditions are examined; the results are then compared with the analytical solutions  and the NURBS-based IGA based on the CPT and the FSDT. The three boundary conditions are as follows: (1) fully simply supported (SSSS); (2) left and right boundaries are simply supported, whereas bottom and upper boundaries are free (SFSF); (3) left and right boundaries are simply supported, whereas bottom and upper boundaries are clamped (SCSC). The plate length is m, and a mesh with control points is used.
Tables 1 and 2 present the first five modes of the normalized frequencies for and , respectively. The exact solutions  and the results obtained by the numerical approaches based on the CPT and the FSDT are also presented for comparison. The frequencies computed by the present method match well with the exact solutions and the solutions by the NURBS-based IGA based on the FSDT. The results from the present method are better than those obtained from the NURBS-based IGA based on the CPT.
Higher order basis functions can be easily obtained in the IGA, thus the formulation of CPT based on the IGA can be simplified. For nonhomogeneous thin plate such as FGM, the midplane displacements are neglected in the CPT, so the CPT is not suitable for analyzing obviously nonhomogeneous thin plate such as FGM thin plates with moderate gradient indices. To eliminate the effects of midplane displacements without introducing new unknown variables, the physical neutral surface is integrated into the CPT, that is, CPT_neu. The NURBS-based isogeometric analysis with the CPT_neu is developed to solve natural frequencies of FGM thin plates. Numerical simulations indicate that the accuracy of the CPT_neu is almost the same as that of the FSDT. The CPT_neu has one degree of freedom in one point, whereas the FSDT has five degrees of freedom in one point. Hence, the result of the CPT_neu is much more efficient than that of the FSDT.
This work was supported by the National Natural Science Foundation of China (Grant no. 51179063) and the Fundamental Research Funds for the Central Universities (Grant no. 2011B03014) and Jiangsu Province Graduate Students Research and Innovation Plan (Grant no. CXZZ13_023).
- S. Shojaee, E. Izadpanah, N. Valizadeh, and J. Kiendl, “Free vibration analysis of thin plates by using a NURBS-based isogeometric approach,” Finite Elements in Analysis and Design, vol. 61, pp. 23–34, 2012.
- C. H. Thai, H. Nguyen-Xuan, N. Nguyen-Thanh, T. H. Le, T. Nguyen-Thoi, and T. Rabczuk, “Static, free vibration, and buckling analysis of laminated composite Reissner-Mindlin plates using NURBS-based isogeometric approach,” International Journal For Numerical Methods in Engineering, vol. 91, no. 6, pp. 571–603, 2012.
- J. Kiendl, K.-U. Bletzinger, J. Linhard, and R. Wüchner, “Isogeometric shell analysis with Kirchhoff-Love elements,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 49–52, pp. 3902–3914, 2009.
- S. Abrate, “Functionally graded plates behave like homogeneous plates,” Composites Part B, vol. 39, no. 1, pp. 151–158, 2008.
- D.-G. Zhang and Y.-H. Zhou, “A theoretical analysis of FGM thin plates based on physical neutral surface,” Computational Materials Science, vol. 44, no. 2, pp. 716–720, 2008.
- T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, “Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement,” Computer Methods in Applied Mechanics and Engineering, vol. 194, no. 39–41, pp. 4135–4195, 2005.
- J. A. Cottrell, A. Reali, Y. Bazilevs, and T. J. R. Hughes, “Isogeometric analysis of structural vibrations,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 41–43, pp. 5257–5296, 2006.
- S. Shojaee, N. Valizadeh, E. Izadpanah, T. Bui, and T.-V. Vu, “Free vibration and buckling analysis of laminated composite plates using the NURBS-based isogeometric finite element method,” Composite Structures, vol. 94, no. 5, pp. 1677–1693, 2012.
- D. J. Benson, Y. Bazilevs, M. C. Hsu, and T. J. R. Hughes, “Isogeometric shell analysis: the Reissner-Mindlin shell,” Computer Methods in Applied Mechanics and Engineering, vol. 199, no. 5-8, pp. 276–289, 2010.
- J. Kiendl, K.-U. Bletzinger, J. Linhard, and R. Wüchner, “Isogeometric shell analysis with Kirchhoff-Love elements,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 49–52, pp. 3902–3914, 2009.
- Y. Bazilevs and T. J. R. Hughes, “NURBS-based isogeometric analysis for the computation of flows about rotating components,” Computational Mechanics, vol. 43, no. 1, pp. 143–150, 2008.
- Y. Bazilevs, V. M. Calo, Y. Zhang, and T. J. R. Hughes, “Isogeometric fluid-structure interaction analysis with applications to arterial blood flow,” Computational Mechanics, vol. 38, no. 4-5, pp. 310–322, 2006.
- Y. Bazilevs, V. M. Calo, T. J. R. Hughes, and Y. Zhang, “Isogeometric fluid-structure interaction: theory, algorithms, and computations,” Computational Mechanics, vol. 43, no. 1, pp. 3–37, 2008.
- C. V. Verhoosel, M. A. Scott, T. J. R. Hughes, and R. de Borst, “An isogeometric analysis approach to gradient damage models,” International Journal for Numerical Methods in Engineering, vol. 86, no. 1, pp. 115–134, 2011.
- W. A. Wall, M. A. Frenzel, and C. Cyron, “Isogeometric structural shape optimization,” Computer Methods in Applied Mechanics and Engineering, vol. 197, no. 33–40, pp. 2976–2988, 2008.
- S.-M. Hu, Y.-F. Li, T. Ju, and X. Zhu, “Modifying the shape of NURBS surfaces with geometric constraints,” CAD Computer Aided Design, vol. 33, no. 12, pp. 903–912, 2001.
- A. Hasani Baferani, A. R. Saidi, and E. Jomehzadeh, “An exact solution for free vibration of thin functionally graded rectangular plates,” Proceedings of the Institution of Mechanical Engineers, Part C, vol. 225, no. 3, pp. 526–536, 2011.
- J. N. Reddy, “Analysis of functionally graded plates,” International Journal for Numerical Methods in Engineering, vol. 47, no. 1–3, pp. 663–684, 2000.
- X. Zhao, Y. Y. Lee, and K. M. Liew, “Free vibration analysis of functionally graded plates using the element-free kp-Ritz method,” Journal of Sound and Vibration, vol. 319, no. 3–5, pp. 918–939, 2009.
- S. N. Atluri and T. Zhu, “A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics,” Computational Mechanics, vol. 22, no. 2, pp. 117–127, 1998.
- S. Hosseini-Hashemi, M. Fadaee, and S. R. Atashipour, “A new exact analytical approach for free vibration of ReissnerMindlin functionally graded rectangular plates,” International Journal of Mechanical Sciences, vol. 53, no. 1, pp. 11–22, 2011.