Document Type: Research Paper
Authors
^{1} Young Researchers and Elite Club,Khomeinishahr Branch, Islamic Azad University, Khomeinishahr, Iran
^{2} Department of Mechanical Engineering, University of Jiroft, Jiroft, Iran
^{3} Department of Water Engineering, Isfahan University of Technology, Isfahan, Iran
Abstract
Keywords

Mechanics of Advanced Composite Structures 4 (2017) 153168 

Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://MACS.journals.semnan.ac.ir 
Meshfree Dynamic Analyses of FGM Sandwich Plates Resting on A Pasternak Elastic Foundation
R. MoradiDastjerdi ^{a}^{*}, H. MomeniKhabisi ^{b}, R. Baghbani ^{c}
^{a }Young Researchers and Elite Club,Khomeinishahr Branch, Islamic Azad University, Khomeinishahr, Iran
^{b }Department of Mechanical Engineering, University of Jiroft, Jiroft, Iran
^{c }Department of Water Engineering, Isfahan University of Technology, Isfahan, Iran
Paper INFO 

ABSTRACT 
Paper history: Received 20170413 Revised 20170512 Accepted 20170603 
This study analyzes the free vibration, forced vibration, resonance, and stress wave propagation of orthotropic sandwich plates made of functionally graded materials (FGMs). Dynamic analyses are conducted using a meshfree method based on firstorder shear deformation theory and the shape functions constructed using moving least squares approximation. The sandwich plates are rested on a Pasternak elastic foundation and subjected to periodic or impact loading and essential boundary conditions, which are imposed through a transfer function method. The sandwich plates are assumed to be composed of a homogeneous orthotropic core and two orthotropic FGM face sheets made of two orthotropic materials. The volume fractions of the materials are varied smoothly along the thickness of the face sheets. The convergence and accuracy of the applied method are demonstrated, after which numerical analyses are conducted to investigate the effects of elastic foundation coefficients, material distributions, geometrical dimensions, timedependent loading, and boundary conditions on the vibrational and dynamic characteristics of the orthotropic FGM sandwich plates. 



Keywords: Wave propagation Sandwich plate Functional graded material Meshfree Firstorder shear deformation theory


DOI: 10.22075/MACS.2017.11043.1107 
© 2017 Published by Semnan University Press. All rights reserved. 
A new class of materials known as functionally graded materials (FGMs) have attracted considerable attention as advanced structural materials for many structural members. FGMs are heterogeneous composite materials with gradient compositional variations in their constituents (e.g., metal and ceramic) from one surface of the materials to another. These variations result in continuously varying material properties. The materials are designed in such a way that they possess desirable properties for specific applications. The concept of FGMs, initially developed for superheatresistant materials for use in space planes or nuclear fusion reactors, is currently of interest to designers of functional materials for energy converters, dental and orthopedic implants, sensors, and thermogenerators [1]. Generally, numerous engineering problems can be modeled as thick plates on elastic foundations, such as the mat foundations of buildings, the pavement underneath roads, and the bases of heavy machines. The mechanical behavior of elastic foundations was widely discussed by Winkler [2] as a oneparameter model or linear model and by Pasternak [3] as a twoparameter model.
Twodimensional plate theories, including classical plate theory (CPT), firstorder shear deformation theory (FSDT), and higherorder shear deformation theories (HSDTs), are commonly used in plate analysis. CPT provides reasonable results for thin plates, but it neglects transverse shear deformation effects, underestimates deflections, and overestimates frequencies, as well as the buckling loads of moderately thick plates [4]. To overcome the limitations of CPT, researchers have developed many shear deformation plate theories that account for transverse shear deformation effects. The theories put forward by Reissner [5] and Mindlin [6] are known as FSDTs. FSDT provides a sufficiently accurate description of response for thin to moderately thick plates [7]. It is widely used in the finite element analyses of composite shells and plates because of its acceptable accuracy and computational complexity [9–11]. However, its performance strongly depends on shear correction factors, which are sensitive not only to material and geometric properties but also to loading and boundary conditions. To avoid the use of shear correction factors and include the actual crosssection warping of a plate into analyses, scholars extensively developed HSDTs, with consideration for the higherorder variations in inplane displacement throughout the thickness of plates [8].
Aiello and Ombres [12] presented an analytical approach to evaluating the buckling load of sandwich panels made of hybrid laminated faces and a transversely flexible core. The authors applied the a priori assumption of the displacement field throughout the thickness of the panels; that is, the authors superpositioned symmetric and antisymmetric components, aside from adopting a pure compressive mode. Ferreira et al. [13] presented the static deformations and free vibrations of shear flexible isotropic and laminated composite plates with an FSDT theory based on a highorder collocation method. They also analyzed isotropic and laminated plates by using Kansa’s nonsymmetric radial basis function collocation method, which is based on an HSDT [14]. Shariyat [15] investigated the effects of thermopiezoelasticity on dynamic buckling under suddenly applied thermal and mechanical loads for imperfect rectangular composite plates with surfacebonded or embedded piezoelectric sensors and actuators. The author developed a finite element formulation grounded in an HSDT and considered the temperature dependence of material properties. Other researchers used HSDT and the general von Karman equation as bases in analyzing the nonlinear vibration and dynamic response of an FGM plate with piezoelectric actuators [16] and surfacebonded piezoelectric fiberreinforced composite actuators [17]. Malekzadeh et al. [18] investigated the dynamic response of thick laminated annular sector plates with simply supported radial edges subjected to a radially distributed line load moving along the circumferential direction. The authors used a threedimensional (3D) hybrid method composed of a series solution, layerwise theory, and a differential quadrature method (DQM), in conjunction with the finite difference method. Zhang et al. [19] employed Reddy’s thirdorder shear deformation theory and the Galerkin procedure in determining the nonlinear dynamics and chaos that underlie a simply supported orthotropic FGM rectangular plate in a thermal environment. The researchers subjected the plate to parametric and external excitations. Dehghan and Baradaran [20] used a combination of finite element methods and DQMs to solve the eigenvalue (buckling and free vibration) equations of rectangular thick plates resting on Pasternak elastic foundations. Thai and Choi [21, 22] developed a refined plate theory for the free vibration and buckling analyses of FGM plates resting on an elastic foundation. Asemi and Shariyat [23] developed a highly accurate nonlinear 3D energybased finite element elasticity formulation for the buckling investigation of anisotropic FGM plates with arbitrary orthotropic directions. To achieve the most accurate results, the authors used a fully compatible Hermitian element with 168 degrees of freedom, which satisfies the continuity of strain and stress components at the mutual edges and nodes of the element a priori. Mansouri and Shariyat [24] conducted a thermomechanical buckling analysis of orthotropic auxetic plates (with negative Poisson ratios) resting on an elastic foundation in hygrothermal environments. Shariyat and Asemi [25] used a nonlinear finite element method and 3D elasticity theory to probe into the shear buckling of orthotropic heterogeneous FGM plates resting on a Winkler elastic foundation. Sofiyev et al. [26] presented analytical formulations and solutions for the stability analysis of heterogeneous orthotropic truncated conical shells subjected to external (lateral and hydrostatic) pressures with mixed boundary conditions. For this purpose, the researchers used Donnell shell theory. Furthermore, the vibrational behavior of single or multidirectional FGMs and functionally graded carbon nanotubereinforced composite (FGCNTRC) structures was investigated using 3D elasticity theory and generalized DQM [27–35].
Some forms of the meshfree method were used to analyze FGM structures. Qian et al. [36] inquired into the static, free, and forced vibrations of a thick rectangular FGM plate on the basis of higherorder shear and normal deformation theory and a meshless local Petrov–Galerkin (MLPG) method. Liew et al. [37] investigated the active control of laminated composite plates with piezoelectric sensor/actuator patches by employing an elementfree Galerkin (EFG) method and FSDT. They used a simple control algorithm to regulate the dynamic response of laminated plates with distributed sensor/actuator patches through a closed control loop. Lanhe et al. [38] examined the dynamic stability of thick FGM plates subjected to aerothermomechanical loads by using a moving least squares DQM. Rezaei Mojdehi et al. [39] carried out a 3D static and dynamic analysis of thick isotropic FGM plates on the basis of MLPG, which is used to construct 3D moving least squares shape functions. The static deformation, free vibration, and dynamic and stress wave propagation of FGM cylinders were analyzed using the same meshfree method adopted in the current research [40–42]. In these previous works, however, the structures examined were axisymmetric and isotropic cylinders. Lei et al. [43, 44] analyzed the buckling and free vibration of FGCNTRC plates by using the elementfree kpRitz method based on FSDT. MoradiDastjerdi et al. [45–47] carried out a dynamic analysis of functionally graded nanocomposite cylinders. The authors also analyzed the static, vibrational, and dynamic behaviors of functionally graded nanocomposite plates reinforced by wavy nanotubes through the same meshfree method adopted in the present study. In an absorbing work, Yaghoubshahi and Alinia [48] developed an EFG method based on HSDT to eliminate transverse shear locking in the analysis of laminated composite plates. The authors compared their results with those obtained using an EFG procedure based on FSDT. Finally, Zhang et al. [49, 50] proposed an elementfreebased improved moving least squaresRitz method and FSDT to study the buckling behavior of FGCNTRC plates resting on Winkler foundations. The authors also examined the nonlinear bending of the plates as they rested on a twoparameter elastic foundation.
As can be seen in the discussion above, no study has been devoted to the analysis of the free vibration, forced vibration, resonance, and stress wave propagation of orthotropic FGM sandwich plates subjected to periodic or impact loading. To address this deficiency, the present research develops a meshfree method on the basis of FSDT to investigate the dynamic behaviors of orthotropic FGM sandwich plates resting on a twoparameter Pasternak elastic foundation. In the meshfree method, shape functions constructed through a moving least squares method are used to approximate the displacement field in the weak form of a motion equation. A transfer function method is employed in the implementation of essential boundary conditions. The developed method does not increase calculations against EFG [41]. The orthotropic FGM sandwich plates are assumed to be composed of two orthotropic FGM face sheets and a homogeneous orthotropic core. The face sheets are assumed to be made of two orthotropic materials, whose volume fractions are varied smoothly along the thickness of the face sheets. The study also examines the effects of elastic foundation coefficients, material distributions, sandwich plate thickness, face sheet thickness, plate aspect ratio, and timedependent force and boundary conditions on the free vibration, forced vibration, resonance, and stress wave propagation of the plates.
Let us consider an orthotropic sandwich FGM plate of length a, width b, total thickness h, and face sheet thickness h_{f}. The material composition of the plate is shown in Fig. 1. The material properties of the face sheets are assumed to be graded along the thickness of the sheets. The profile of this variation exerts important effects on plate behavior. Several models have been proposed for variations in material properties. Among these, the volume fraction model is the most frequently used. In this model, the volume fraction and material properties of a sandwich plate are varied as follows:
(1) 

(2) 
where P is an indicator of the material properties of a plate; this indicator is used in lieu of modulus elasticity E, Poisson’s ratio , and density . Subscripts 1 and 2 represent the z=0 and z=h constituents, respectively; and n denotes the volume fraction exponent. An n equal to zero represents a sandwich plate made of two different materials in the face sheets and a core without any mixture of materials. By contrast, an infinite n indicates a homogeneous plate made of core constituents. Fig. 2 illustrates the variations in the material volume fraction of the FGM sandwich plate along the thickness of the plate at different volume fraction exponents.
On the basis of FSDT, displacement components can be defined as [4]
(3) 
where u, v, and w are displacements in the x, y, and z directions, respectively; u_{0}, v_{0}, and w_{0} denote the midplane displacements and the and rotations of the normal to the midplane directions about the yaxis and xaxis, respectively. Kinematic relations can be obtained as follows:
(4) 
where
(5) 
Figure 1 Schematic of orthotropic FGM sandwich plate resting on Pasternak elastic foundation
Figure 2 Variations in material volume fraction of orthotropic FGM sandwich plate along the plate thickness at different volume fraction exponents
The linear constitutive relations of a functionally graded plate can be written as
(6) 
in which , , , , and are the normal stress vector, shear stress vector, normal strain vector, shear strain vector, and engineering constants, respectively. denotes the transverse shear correction coefficient, which is set at for homogeneous materials and for FGMs (where for orthotropic FGM structures) [51].
In Eq. (6), as well,
(7) 
Considering the Pasternak foundation model, the total energy of the plate is expressed thus:
(8) 
where q(t) is the timedependent applied load, and k_{w} and k_{s} are the coefficients of Winkler and Pasternak (shear) foundations, respectively. If a foundation is modeled as a linear Winkler foundation, the coefficient k_{s} in Eq. (8) takes the value of zero.
4. MESHFREE NUMERICAL ANALYSES
In these analyses, the moving least squares shape functions introduced by Lancaster and Salkauskas [52] are used to approximate the displacement vector in the weak form of a motion equation. Displacement vector u can be approximated using the shape functions as follows:
(9) 
where N is the total number of nodes, denotes the virtual nodal values vector, and is the moving least squares shape function of the node located at X(x,y)=X_{i}. These variables are defined in the following manner:
(10) 
and
(11) 
In Eq. (11), is the cubic spline weight function, is the base vector, and is the moment matrix. These are defined thus:
(12) 

(13) 
Using the moving least squares shape function enables Eq. (4) to be written as
(14) 
in which
(15) 
For an elastic foundation, and can be defined as follows:
(16) 
Substituting Eqs. (6) and (14) into Eq. (8) yields

(17)
in which the components of extensional stiffness , bendingextensional coupling stiffness , bending stiffness , transverse shear stiffness , and and are introduced into a mass matrix. These variables are defined as
(18) 
and
(19) 
where , , and are the normal, coupled normalrotary, and rotary inertial coefficients, respectively. These coefficients are defined by
(20) 
The arrays of bendingextensional coupling stiffness matrix are zero for symmetric laminated composites.
Finally, through a derivative with respect to displacement vector , Eq. (17) can be expressed as
(21) 
in which , , and are the mass matrix, stiffness matrix, and force vector, respectively. These are defined as
(22) 

(23) 

(24) 
where , , and are the stiffness matrixes of extensional, bendingextensional, and bending modes, respectively. and are the stiffness matrixes that represent Winkler and Pasternak elastic foundations. They are defined in the following equations:
(25) 

(26) 
For numerical integration, the problem domain is discretized to a set of background cells with Gauss points inside each cell. Then, global stiffness matrix is obtained numerically by sweeping all the Gauss points.
The imposition of essential boundary conditions in the system of Eq. (21) is impossible because moving least squares shape functions do not satisfy the Kronecker delta property. As previously stated, this work uses a transfer function method in implementing essential boundary conditions. For this purpose, a transformation matrix is constructed by establishing a relationship between nodal displacement vector and virtual displacement vector .
(27) 
is the transformation matrix that is a 5N×5N matrix. For each node, this matrix is defined as

(28) 
where is an identity matrix of size 5. Using Eq. (27) rearranges the system of linear Eq. (21) to
(29) 
where
(30) 
Now, the essential B. Cs. can be easily enforced in the modified system of Eq. (29), as is possible with the finite element method. In this work, the Newmark (central difference) method is used for the solution of Eq. (29) in the time domain.
In this method, the displacement field is obtained in each time step as follows:
(31) 

(32) 
where t shows the times, and dt is the size of a time step.
By solving Eq. (29), the time history of the displacement field and then stress wave propagation can be derived with the consideration of type of periodic or impact loading. In the absence of external forces, Eq. (29) is simplified in this manner:
(33) 
Thus, the natural frequencies and mode shapes of the plates are determined by solving this eigenvalue problem.
5. RESULTS AND DISCUSSION
This section discusses numerical examples of the free vibration, forced vibration, resonance, and stress wave propagation behaviors of the examined orthotropic FGM sandwich plates. As mentioned earlier, the plates are subjected to periodic or impact loading, and the developed meshfree approach and the Newmark method are used. The frequencies of the sandwich plates are derived and investigated through the consideration of their resonance behaviors.
First, the convergence and accuracy of the meshfree method in determining the vibrational behaviors of the plates are examined by a comparison between the results and those reported in the literature. Second, the meshfree results on the vibrational and dynamic characteristics of the plates are presented.
In simulations, the orthotropic FGM sandwich plates are assumed to be made of a homogeneous glassepoxy core and two FGM face sheets. In the face sheets, the material volume fractions are varied from those of the glassepoxy at the interfaces of the core and from those of the graphiteepoxy at the faces of the sandwich plates. The modifications are conducted using Eq. (1) (Fig. 2). The material properties of the glassepoxy and graphiteepoxy are listed in Table 1. Note that transverse isotropic materials are a special class of orthotropic materials that have the same properties in one plane (e.g., the xy plane) and different properties in the direction normal to the aforementioned plane (e.g., the zaxis). Thus, glassepoxy and graphiteepoxy are transverse isotropic materials. In all the examples of the orthotropic sandwich plates, the foundational parameters are presented in the nondimensional forms K_{w}=k_{w}a^{4}/D and K_{s}=k_{s}a^{2}/D, in which D=E_{1}h^{3}/12(1υ_{12}^{2}) is a reference bending rigidity of the plate and is based on the mechanical properties of graphiteepoxy. The nondimensional deflections (for forced and dynamic analyses) and natural frequencies of the orthotropic sandwich plates are based on the mechanical properties of graphiteepoxy and glassepoxy, respectively. The deflection and frequency are defined as [13]
(34) 

(35) 
where q_{0} is the value of the amplitude of timedependent applied load, and w denotes the central deflection of the plates. The following sets of support conditions are employed to compute the desired responses:
(a) Simply supported conditions
at x=0 and a: v=w=θ_{y}=0 and at y=0 and b: u=w=θ_{x}=0
(b) Clamped conditions:
at x=0 and a: and at y=0 and b: u=v=w=θ_{x}=θ_{y}=0
Table 1 Material properties of the examined orthotropic (transverse isotropic) materials [53]
Materials 
E_{1} (GPa) 
E_{2} (GPa) 
E_{3} (GPa) 
υ_{23} 
υ_{31} 
υ_{12} 
G_{23} (GPa) 
G_{31} (GPa) 
G_{12} (GPa) 
ρ (kg/m^{3}) 
Graphiteepoxy 
155 
12.1 
12.1 
0.458 
0.248 
0.248 
3.2 
4.4 
4.4 
1500 
Glassepoxy 
50 
15.2 
15.2 
0.428 
0.254 
0.254 
3.28 
4.7 
4.7 
1800 
5.1 Validation of models
To investigate the convergence and accuracy of the developed method, let us consider a simply supported FGM square plate, as was done by Thai and Choi [21]. Fig. 3 depicts the convergence of the proposed meshfree method in the nondimensional fundamental frequencies of the plates resting on WinklerPasternak elastic foundations at h/a=0.2, K_{w}=100 and K_{s}=100, and volume fraction exponents of n=0 and n=1. These values reveal that using only a 5×5 node arrangement enables the proposed method to achieve very good accuracy and agreement with the results of Thai and Choi [21] for the homogeneous (n=0) and FGM (n=1) plates. The nondimensional fundamental frequencies of these plates are presented in Table 2 for various values of h/a (0.05, 0.1, and 0.2) and elastic foundation coefficients. The table shows that the proposed method exhibits very good accuracy and agreement with previously reported results, especially with respect to thin plates.
5.2 Free vibration of orthotropic FGM sandwich plates
Clamped orthotropic FGM sandwich plates are examined to investigate the effects of material distributions, elastic foundation coefficients, and geometric dimensions on the natural frequencies of the plates. Table 3 lists the nondimensional fundamental frequencies of the plates at different elastic foundation coefficients (K_{w} and K_{s}), aspect ratios (b/a=1 and 3), plate thickness ratios (h/a=0.1 and 0.2), face sheet thickness ratios (h_{f}/h=0.1 and 0.2), and volume fraction exponents (n=0, 0.1, 1, 10, and 100). The frequency parameter decreases with an increase in b/a from 1 to 3 because the plate manners are near the beam manners. An increase in h/a and h_{f}/h and a decrease in volume fraction exponent increase the frequency parameter. The values of h/a significantly affect the frequency parameter of orthotropic FGM sandwich plates. Finally, the elastic foundation also increases the frequency parameter.
(a)

(b)

Figure 3 Convergence of nondimensional fundamental frequency ( ) of the FGM plate at (a) n=0; (b) n=1 and h/a=0.2, K_{w}=100, K_{s}=100 for different numbers of nodes in each direction
Table 2 Comparison of normalized fundamental frequencies, , in simply supported square FGM plates

Table 3 Nondimensional fundamental frequencies, , in clamped orthotropic FGM sandwich plates
* 2^{nd} frequency: 0.1090, 3^{rd} frequency: 0.1430 ** 2^{nd} frequency: 0.1693, 3^{rd} frequency: 0.1927

5.3 Forced vibration of orthotropic FGM sandwich plates
Square orthotropic FGM sandwich plates are subjected to periodic uniform pressure loading at the top face of the sandwich plates as follows:
(36) 
where is the frequency of timedependent applied load and is equal to in the succeeding simulation.
The effects of loading frequency on the time history of the central deflection of the orthotropic FGM sandwich plates are investigated. Let us consider clamped square plates with h/a=0.1, h_{f}/h=0.1, and n=1 and subjected to timedependent pressure at loading frequencies of , , and . Figs. 4a and 4b show the time histories of the normalized central deflections of the plates without an elastic foundation and with resting on the Pasternak foundation (K_{w}=100 and K_{s}=10), respectively. The figures indicate that the severity of the elastic foundation decreases the central deflections of the sandwich plates and that the amplitudes of the deflections are reduced by a decrease in loading frequency.
(a)

(b)

Figure 4 Time history of normalized central deflections ( ) of clamped orthotropic FGM sandwich plates at (a) K_{w}=0, K_{s}=0; (b) K_{w}=100, K_{s}=10, and b/a=1, h/a=0.1, h_{f}/h=0.1, and n=1
Now, let us consider square clamped sandwich plates subjected to periodic pressure ( ) and resting on the elastic foundation (K_{w}=100 and K_{s}=10) to delve into the effects of the plates’ geometric dimensions on the time history of deflection. Fig. 5 illustrates the time history of the plates’ normalized central deflections at h/a=0.1 and 0.2, h_{f}/h=0.1 and 0.2, and n=0, 0.1, 1, and 10. The amplitudes of the deflections decline with decreasing volume fraction exponent. Such amplitudes also decrease with increasing h/a and h_{f}/h. The thickness of the plates (h/a) is higher than that of the face sheets (h_{f}/h).
Finally, the effects of essential boundary conditions on the forced vibration behaviors of orthotropic FGM sandwich plates subjected to a periodic load of are investigated. Fig. 6 show the time history of the normalized central deflections of the plates at h/a=0.1, h_{f}/h=0.1, K_{w}=100 and K_{s}=10, and n=0, 0.1, 1, and 10 for the boundary conditions of CSCS and CFCF, where C, S, and F denote clamped, simply supported, and free edges, respectively. Comparing Figs. 6 and 5a indicates that the clamped and CFCF sandwich plates have the smallest and largest deflection amplitudes, respectively.
5.4 Resonance behavior of orthotropic FGM sandwich plates
The resonance of orthotropic FGM sandwich plates are also explored. Let us consider sandwich plates subjected to periodic uniform loading at the top face of the sandwich plates, as with Eq. (36), in which loading frequency is equal to the n^{th} natural frequency of the sandwich plates reported in Table 3.
Figs. 7a and 7b illustrate the time histories of the normalized central deflections of these sandwich plates without an elastic foundation and with resting on the Pasternak foundation (K_{w}=100 and K_{s}=10), respectively. These sandwich plates are clamped at h/a=0.1, h_{f}/h=0.1, and n=1 and subjected to periodic loading with loading frequencies equal to the first, second, and third frequencies of the sandwich plates (i.e., ). This loading leads to a divergence in the amplitudes of deflection in the first mode despite the elastic foundation’s reduction of deflection amplitudes.
Fig. 8 presents the first mode of resonance in the clamped orthotropic sandwich plates at b/a=1, h/a=0.1 and 0.2, h_{f}/h=0.1 and 0.2, and n=0, 0.1, 1, and 10. An increase in h/a elevates the amplitudes of deflection and decreases the periods of their vibration because of the increase in frequency loading. The volume fraction exponent exerts a stronger effect on the periods of vibrations through an increase in h_{f}/h.
5.5 Dynamic behavior of orthotropic FGM sandwich plates
Now, the dynamic behaviors of orthotropic FGM sandwich plates resting on an elastic foundation are determined by considering sandwich plates subjected to an impact uniform pressure load at the top face of the plates as follows:
(35) 
Let us consider clamped orthotropic FGM sandwich plates subjected to an impact load, as with Eq. (37), and resting on the Pasternak elastic foundation at b/a=1, h/a=0.1, h_{f}/h=0.1, K_{w}=100, and K_{s}=10 (a=1 m, q_{0}=0.1 MPa). Fig. 9 shows the time history of stress wave propagation ( , , , , and ) at the top face and central deflection of the sandwich plates for n=0.1, 1, and 10. The amplitudes of normal stresses are higher than those of shear stresses. An increasing volume fraction exponent leads to a rise in the amplitudes of stresses and a decline in wave speed. After load is eliminated, the sandwich plates show harmonic vibration with amplitudes less than those of the same forced vibration (in Fig. 5a).
(a)

(b)

(c)

(d)

Figure 5 Time history of normalized central deflections ( ) of clamped orthotropic FGM sandwich plates at (a) h/a=0.1, h_{f}/h=0.1; (b) h/a=0.1, h_{f}/h=0.2; (c) h/a=0.2, h_{f}/h=0.1; (d) h/a=0.1, h_{f}/h=0.1, and b/a=1, K_{w}=100, K_{s}=10, and
(a)

(b)

Figure 6 Time history of normalized central deflections ( ) of (a) CSCS, (b) CFCF orthotropic FGM sandwich plates at b/a=1, h/a=0.1, h_{f}/h=0.1, K_{w}=100, K_{s}=10, and
(a)

(b)

Figure 7 Time history of normalized central deflections ( ) of orthotropic FGM sandwich plates at (a) K_{w}=0, K_{s}=0; (b) K_{w}=100, K_{s}=10, and b/a=1, h/a=0.1, h_{f}/h=0.1, and n=1
(a)

(b)

(c)

(d)

Figure 8 Time history of normalized central deflections ( ) of clamped orthotropic FGM sandwich plates at (a) h/a=0.1, h_{f}/h=0.1; (b) h/a=0.1, h_{f}/h=0.2; (c) h/a=0.2, h_{f}/h=0.1; (d) h/a=0.1, h_{f}/h=0.1, and b/a=1, K_{w}=100, K_{s}=10, and
(a)

(b)

(c)

(d)

(e)

(f)

Figure 9 Time history of (a) , (b) , (c) , (d) , (e) , and (f) at the top face of clamped orthotropic FGM sandwich plates at b/a=1, h/a=0.1, h_{f}/h=0.1, K_{w}=100, and K_{s}=10
(a)

(b)

(c)

(d)

(e)

(f)

Figure 10 Time history of (a) , (b) , (c) , (d) , (e) , and (f) at the top face of clamped orthotropic FGM sandwich plates at b/a=1, h/a=0.1, h_{f}/h=0.2, K_{w}=100, and K_{s}=10
Let us consider the same sandwich plates but with h_{f}/h=0.2 as the variable. Fig. 10 displays the stress wave propagation and central vibration of the plates. The comparison of Figs. 10 and 9 indicates that the plates exhibit almost the same dynamic behaviors but that the volume fraction exponent exerts a stronger effect on the plates with thicker face sheets.
Finally, let us reexamine the first model of dynamic analysis but with h/a=0.2 as the variable. The time history of stress wave propagation and the central vibration of the orthotropic FGM sandwich plates are illustrated in Fig. 11. The comparison of Figs. 9 and 11 demonstrates that an increase in the thickness of the sandwich plates leads to a decrease in the amplitudes of stresses and vibrations and the speeds of wave propagation. Fig. 12 shows the time history of the normal stresses of the sandwich plate with n=10 imposed at the top, middle, and bottom planes (z=−h, 0, h). The sandwich plate at z=0 senses almost no normal stresses, whereas the top and bottom planes sense the same inplane stresses but with a different sign.
(a)

(b)

(c)

(d)

Figure 11 Time history of (a) , (b) , (c) , and (d) at the top face of clamped orthotropic FGM sandwich plates at b/a=1, h/a=0.2, h_{f}/h=0.1, K_{w}=100, and K_{s}=10
(a)

(b)

Figure 12 Time history of (a) and (b) at the top, middle, and bottom planes of clamped orthotropic FGM sandwich plates at n=10, b/a=1, h/a=0.1, h_{f}/h=0.2, K_{w}=100, and K_{s}=10
6. CONCLUSION
In this study, the vibration, resonance, and dynamic behaviors of orthotropic FGM sandwich plates resting on a Pasternak elastic foundation are analyzed using a developed meshfree method based on FSDT and moving least squares shape functions. Essential boundary conditions are imposed via a transfer function method. The sandwich plates are assumed to be composed of a homogeneous orthotropic core and two orthotropic FGM face sheets made of two orthotropic materials. Numerical examples are provided to illuminate the vibrational and dynamic characteristics of the plates at different elastic foundation coefficients, material distributions, geometrical dimensions, timedepended loading levels, and boundary conditions. The primary findings are summarized as follows:
References
[1] Koizumi M. The concept of FGM. Ceram Trans Func Grad Mater 1993; 34: 3–10.
[2] Winkler E. 1867, Die Lehre von der Elasticitaet und Festigkeit, Dominicus, Prag.
[3] Pasternak PL. 1954, On a new method of analysis of an elastic foundation by means of two foundation constants (in Russian), Gosudarstrennoe Izdatelstvo Literaturi po Stroitelstvu i Arkhitekture, Moscow, USSR.
[4] Reddy JN. 2004, Mechanics of Laminated Composite Plates and Shells: Theory and Analysis, CRC.
[5] Reissner E. The effect of transverse shear deformation on the bending of elastic plates, J Applied Mech, 1945; 12: 69–72.
[6] Mindlin RD. Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates, J Applied Mech, 1951; 18: 31–38.
[7] Thai HT, Choi DH. An efficient and simple refined theory for buckling analysis of functionally graded plates. Applied Math Modelling, 2012; 36: 1008–1022.
[8] Reddy JN. A simple higher order theory for laminated composite plates. J Applied Mech, 1984; 51: 745752.
[9] Kundalwal, SI, Kumar, RS, Ray, MC, Smart damping of laminated fuzzy fiber reinforced composite shells using 1–3 piezoelectric composites, Smart Mater Struct, 2013: 22, 105001.
[10] Kundalwal, SI, Meguid, SA, Effect of carbon nanotube waviness on active damping of laminated hybrid composite shells, Acta Mech, 2015; 226, 2035–2052.
[11] Kundalwal, SI, Ray MC, Smart damping of fuzzy fiber reinforced composite plates using 13 piezoelectric composites, J Vib Control, 2016; 22: 15261546.
[12] Aiello MA, Ombres L. Buckling Load Design of Sandwich Panels Made with Hybrid Laminated Faces and Transversely Flexible Core. J Sandwich Struct Mater, 2007; 9: 467485.
[13] Ferreira AJM, Castro LMS, Bertoluzza S. A high order collocation method for the static and vibration analysis of composite plates using a firstorder theory. Compos Struct, 2009; 89: 424–432.
[14] Ferreira AJM, Roque CMC, Martins PALS. Analysis of composite plates using higherorder shear deformation theory and a finite point formulation based on the multiquadric radial basis function method. Composites Part B 2003; 34: 627–36.
[15] Shariyat M. Dynamic buckling of imperfect laminated plates with piezoelectric sensors and actuators subjected to thermoelectromechanical loadings, considering the temperaturedependency of the material properties. Compos Struct, 2009; 88: 228–239.
[16] Huang XL, Shen HS. Vibration and dynamic response of functionally graded plates with piezoelectric actuators in thermal environments. J Sound Vibr, 2006; 289: 25–53.
[17] Xia XK, Shen HS. Nonlinear vibration and dynamic response of FGM plates with piezoelectric fiber reinforced composite actuators, Compos Struct, 2009; 90: 254–262.
[18] Malekzadeh P, Golbahar Haghighi MR, Gholami M, Dynamic response of thick laminated annular sector plates subjected to moving load, Compos Struct, 2010; 92: 155–163.
[19] Zhang W, Yang J, Hao Y. Chaotic vibrations of an orthotropic FGM rectangular plate based on thirdorder shear deformation theory, Nonlinear Dyn, 2010; 59: 619–660
[20] Dehghan M, BaradaranGH. Buckling and free vibration analysis of thick rectangular plates resting on elastic foundation using mixed finite element and differential quadrature method, Applied Math Comput, 2011; 218: 2772–2784.
[21] Thai HT, Choi DH. A refined plate theory for functionally graded plates resting on elastic foundation, Compos Sci Technol, 2011; 71: 1850–1858.
[22] Thai HT, Choi DH. An efficient and simple refined theory for buckling analysis of functionally graded plates, Applied Math Modelling, 2012; 36: 1008–1022.
[23] Asemi K, Shariyat M. Highly accurate nonlinear threedimensional finite element elasticity approach for biaxial buckling of rectangular anisotropic FGM plates with general orthotropy directions, Compos Struct, 2013; 106: 235249.
[24] Mansouri MH, Shariyat M. Biaxial thermomechanical buckling of orthotropic auxetic FGM plates with temperature and moisture dependent material properties on elastic foundations, Composites Part B 2015; 83: 88104.
[25] Shariyat M, Asemi K. Threedimensional nonlinear elasticitybased 3D cubic Bspline finite element shear buckling analysis of rectangular orthotropic FGM plates surrounded by elastic foundations, Composites Part B 2014; 56: 934947.
[26] Sofiyev AH, Huseynov SE, Ozyigit P, Isayev FG. The effect of mixed boundary conditions on the stability behavior of heterogeneous orthotropic truncated conical shells, Meccanica 2015; 50: 21532166.
[27] Tahouneh, V, EskandariJam, J, A Semianalytical Solution for 3D Dynamic Analysis of Thick Continuously Graded Carbon Nanotubereinforced Annular Plates Resting on a Twoparameter Elastic Foundation, Mech Adv compos struct, 2014; 1, 113130.
[28] Tahouneh, V, Naei, MH, SemiAnalytical Solution for Free Vibration Analysis of Thick Laminated Curved Panels with PowerLaw Distribution FG Layers and Finite Length Via TwoDimensional GDQ Method, J solid mech, 2016; 8, 334347.
[29] Tahouneh, V, Naei, MH, The effect of multidirectional nanocomposite materials on the vibrational response of thick shell panels with finite length and rested on twoparameter elastic foundations, Int J Adv Struct Eng, 2016; 8, 1128.
[30] Tahouneh, V, Using an equivalent continuum model for 3D dynamic analysis of nanocomposite plates, Steel Compos Struct, 2016; 20, 623649.
[31] Tahouneh, V, Naei, MH, Free vibration and vibrational displacements analysis of thick elastically supported laminated curved panels with powerlaw distribution functionally graded layers and finite length via 2D GDQ method, JSand Struct Mater, 2016; 18, 263293.
[32] Tahouneh, V, Naei, MH, Using EshelbyMoriTanaka scheme for 3D free vibration analysis of sandwich curved panels with functionally graded nanocomposite face sheets and finite length, Polym Compos, 2016; DOI: 10.1002/pc. 23929.
[33] Tahouneh, V, Yas, MH, Tourang, H, Kabirian, M, Semianalytical solution for threedimensional vibration of thick continuous grading fiber reinforced (CGFR) annular plates on Pasternak elastic foundations with arbitrary boundary conditions on their circular edges, Meccanica, 2013; 48, 13131336.
[34] Tahouneh, V, Yas, MH, Influence of equivalent continuum model based on the EshelbyMoriTanaka scheme on the vibrational response of elastically supported thick continuously graded carbon nanotubereinforced annular plates, Polym Compos, 2014; 35, 16441661.
[35] Tahouneh, V, Yas, MH, Semianalytical solution for threedimensional vibration analysis of thick multidirectional functionally graded annular sector plates under various boundary conditions, J Eng Mech, 2013, 140, 3146.
[36] Qian LF, Batra RC, Chen LM. Static and dynamic deformations of thick functionally graded elastic plates by using higherorder shear and normal deformable plate theory and meshless local Petrov–Galerkin method, Compos Part B 2004; 35: 685–697.
[37] Liew KM, He XQ, Tan MJ, Lim HK. Dynamic analysis of laminated composite plates with piezoelectric sensor/actuator patches using the FSDT meshfree method, Inter J Mech Sci, 2004; 46: 411–431.
[38] Lanhe W, Hongjun W, Daobin W, Dynamic stability analysis of FGM plates by the moving least squares differential quadrature method, Compos Struct, 2007; 77: 383–394.
[39] Rezaei Mojdehi A, Darvizeh A, Basti A, Rajabi H. Three dimensional static and dynamic analysis of thick functionally graded plates by the meshless local Petrov–Galerkin (MLPG) method, Eng Anal Bound Elements, 2011; 35: 1168–1180
[40] Foroutan M, MoradiDastjerdi R, SotoodehBahreini R. Static analysis of FGM cylinders by a meshfree method, Steel Compos Struct, 2012; 12: 111.
[41] Mollarazi HR, Foroutan M, MoradiDastjerdi R. Analysis of free vibration of functionally graded material (FGM) cylinders by a meshless method, J Compos Mater 2012; 46: 507–15.
[42] Foroutan M, MoradiDastjerdi R. Dynamic analysis of functionally graded material cylinders under an impact load by a meshfree method, Acta Mech, 2011; 219: 28190.
[43] Lei ZX, Liew KM, Yu JL. Buckling analysis of functionally graded carbon nanotubereinforced composite plates using the elementfree kpRitz method, Compos Struct, 2013; 98:160–168.
[44] Lei ZX, Liew KM, Yu JL. Free vibration analysis of functionally graded carbon nanotubereinforced composite plates using the elementfree kpRitz method in thermal environment, Compos Struct, 2013; 106: 128–138.
[45] MoradiDastjerdi R, MomeniKhabisi H, Dynamic analysis of functionally graded nanocomposite plates reinforced by wavy carbon nanotube. Steel Compos Struct, 2016; 22: 277299.
[46] MoradiDastjerdi R, Pourasghar A. Dynamic analysis of functionally graded nanocomposite cylinders reinforced by wavy carbon nanotube under an impact load. J Vib Control 2016; 22: 10621075.
[47] MoradiDastjerdi R, Payganeh G, Rajabizadeh Mirakabad S, Jafari Mofrad Taheri M, Static and Free Vibration Analyses of Functionally Graded Nanocomposite Plates Reinforced by Wavy Carbon Nanotubes Resting on a Pasternak Elastic Foundation. Mech Adv Compos Struct, 2016; 3: 123135.
[48] Yaghoubshahi M, Alinia MM. Developing an element free method for higher order shear deformation analysis of plates, ThinWalled Struct, 2015; 94: 225–233.
[49] Zhang LW, Lei ZX, Liew KM, An elementfree IMLSRitz framework for buckling analysis of FG–CNT reinforced composite thick plates resting on Winkler foundations, Eng Anal Bound Elements, 2015; 58: 7–17.
[50] Zhang LW, Song ZG, Liew KM, Nonlinear bending analysis of FGCNT reinforced composite thick plates resting on Pasternak foundations using the elementfree IMLSRitz method, Compos Struct, 2015; 128: 165–175.
[51] Efraim E, Eisenberger M, Exact vibration analysis of variable thickness thick annular isotropic and FGM plates. J Sound Vib, 2007; 299: 720–38.
[52] Lancaster P, Salkauskas K, Surface Generated by Moving Least Squares Methods, Math Comput,1981; 37: 14158.
[53] Hyer MW. 1998, Mechanics of composite materials. McGrawHill.
[54] Baferani AH, Saidi AR, Ehteshami H, Accurate solution for free vibration analysis of functionally graded thick rectangular plates resting on elastic foundation. Compos Struct, 2011; 93: 1842–53.