Isogeometric analysis for functionally graded microplates based on
modified couple stress theory
Hoang X. Nguyena , Tuan N. Nguyenb , M. Abdel-Wahabc , S.P.A. Bordasd , H.
Nguyen-Xuane , Thuc P. Voa,∗

arXiv:1604.00547v1 [] 2 Apr 2016

a

Department of Mechanical and Construction Engineering, Northumbria University, Newcastle upon Tyne
NE1 8ST, United Kingdom
b
Department of Computational Engineering, Vietnamese-German University, Binh Duong New City,
Vietnam
c
Soete Laboratory, Faculty of Engineering and Architecture, Ghent University, 9000 Ghent, Belgium
d
Faculte des Sciences, de la Technologie et de la Communication, University of Luxembourg, Luxembourg
e
Center for Interdisciplinary Research in Technology (CIRTech), Ho Chi Minh University of Technology
(HUTECH), Ho Chi Minh City 700000, Vietnam

Abstract
Analysis of static bending, free vibration and buckling behaviours of functionally graded
microplates is investigated in this study. The main idea is to use the isogeometric analysis
in associated with novel four-variable refined plate theory and quasi-3D theory. More importantly, the modified couple stress theory with only one material length scale parameter
is employed to effectively capture the size-dependent effects within the microplates. Meanwhile, the quasi-3D theory which is constructed from a novel seventh-order shear deformation
refined plate theory with four unknowns is able to consider both shear deformations and
thickness stretching effect without requiring shear correction factors. The NURBS-based
isogeometric analysis is integrated to exactly describe the geometry and approximately calculate the unknown fields with higher-order derivative and continuity requirements. The
convergence and verification show the validity and efficiency of this proposed computational
approach in comparison with those existing in the literature. It is further applied to study the
static bending, free vibration and buckling responses of rectangular and circular functionally
graded microplates with various types of boundary conditions. A number of investigations
are also conducted to illustrate the effects of the material length scale, material index, and
length-to-thickness ratios on the responses of the microplates.
Keywords: Functionally graded microplates, Modified couple stress theory, Refined plate
theory, Quasi-3D theory, Isogemetric analysis.

∗

Corresponding author
Email addresses: xuan.h.nguyen@northumbria.ac.uk (Hoang X. Nguyen),
ngx.hung@hutech.edu.vn (H. Nguyen-Xuan), thuc.vo@northumbria.ac.uk (Thuc P. Vo )

Preprint submitted to Elsevier

April 5, 2016

1. Introduction
Functionally graded materials (FGMs), which are initially proposed by group of Japanese
aerospace researchers, are composite materials formed of two or more constituent phases in
which material properties vary smoothly from one surface to the other. Consequently, FGMs
are able to avoid high interlaminar shear stresses, stress concentration and delamination
phenomenon which are often cited as shortcomings of laminated composite materials. A
FGM consisting of ceramic and metal possesses higher thermal resistance and better ductility
which are inherited from the ceramic and metal phases, respectively. Owning these striking
features, FGMs are applicable to various fields of engineering including aerospace, nuclear
power, chemistry and biology engineering. The FGMs are also widely studied within various
types of structure such as beams [1–4], plates [5–8], and shells [9–11].
Recent advances in technology lead to new industrial fields in which small-scale elements
are involved. Such elements have been applied in the micro- and nano-electro-mechanical
system [12, 13], actuators [14], space and bio-engineering [15]. These applications encourage
new research area that focuses on investigating and predicting the behaviours of such micro
structures. A number of attempts have been conducted to analyse characteristic of smallscale structures including experiment and computer simulation [16, 17]. Nevertheless, while
the former is mechanical and economical challenging due to the extremely small size of the
elements, the latter exponentially raises computational cost. Therefore, the mathematical
and mechanical models which account for the small-scale effects are needed. It is widely
known that the behaviour of small-scale structures is size-dependent. This is confirmed
following a number of theoretical and experimental works of Fleck et al. [18], Stolken and
Evans [19], and Lam et al. [20]. Having an experimental observation from bending test of
epoxy polymeric microbeams, Lam et al. [20] point out that the bending rigidity increases
2.4 times as a result of the reduction of beam thickness from 115 µm to 20 µm.
In order to take into account the size effects, based on the classical elasticity theories, a
number theories have been developed including nonlocal elasticity theory [21], strain gradient theory [22], and modified couple stress theory [23]. It is worth noting that the classical
elasticity is fundamentally founded by the introduction of the Hook’s Law in which the force
and the change in displacement are linearly related via the stiffness of the component where
the forces are applied. This principal of physics governs the linearly elastic behaviours of
materials. As physics and mechanics develop, stress and strain which are linearly connected
via the elastic modulus are cited as the primary measures of classical elasticity rather than
force and displacement counterparts. Aiming for better prediction of materials’ responses,
Mindlin and Tiersten [24] and Mindlin [25] have initially developed higher-order theories of
elasticity. Based on the employment of deformation metrics, those theories could be classified into two categories including strain gradient and couple stress theories. With regard
to the strain gradient theories, this concept was first developed by Fleck and Hutchinson
[22] in which the Mindlin’s theory [25] was extended. There are two components, which are
classified using the second-order deformation tensor that include stretch gradient tensor and
rotation gradient tensor [26]. Within the concept of couple stress theory, both strain and
curvature jointly govern the strength of the solid. In addition, while the antisymmetric part
2

of the second-order deformation gradients is served as rotation gradients, the symmetric
part is neglected. Based on the initial ideas of couple stress theory, a number of attempts
have been conducted to further develop the concepts that are applicable for size-dependent
problems. Yang et al. [23] proposed the equilibrium of moment of couples which was an
additional equilibrium relation that forces the couple stress tensor to be symmetric. Therefore, the deformation energy is only influenced by the symmetric part of rotation gradient
and the symmetric part of displacement gradient. In addition, instead of using two material
length scale parameters as need in the classical couple stress theory, this modified couple
stress theory (MCST) requires only one material length scale parameter to construct the
constitutive relation. Park and Gao [27] utilised the principle of minimum total potential energy to develop a variational formulation of the modified couple stress theory. This method
not only derives the equilibrium equations but also forms the determination of boundary
conditions which are not available in the Yang’s theory [23]. In recent years, the modified couple stress theory has been applied to study the various behaviours of small-scale
structures. For beam analysis, static bending, buckling, and vibration analysis have been
conducted using Euler-Bernoulli [28], Timoshenko [29–32], and higher-order beam theories
[33]. The modified couple stress theory was also applied in the small-scale plate analysis
in several ways, Tsiatas [34] initially employed MCST to investigate the static bending response of isotropic Kirchhoff microplates. Yin et al. [35] solely investigated the vibration
behaviour of Kirchhoff microplates using the standard separation of variables to derive the
closed-form solution for natural frequency. Bending and vibration behaviours of Mindlin
microplates were studied by Ma et al. [36] in which the thickness stretching effect was also
taken into account. With regard to the small-scale functionally graded (FG) structures, a
number of investigations have been conducted for micro FG beams and plates. Simsek et al.
[37] investigated the static bending of functionally graded microbeams based on the modified couple stress theory. Buckling behaviour of FG microbeams was studied by Nateghi
et al. [38]. Based on the modified couple stress theory, Reddy et al. [39–42] studied nonlinear behaviour of small-scale FG microplates with different geometry shapes. Thai and
his colleagues utilised the Navier’s approach to derive solutions for FG microplates in which
Kirchoff, Minlind and sinusoidal plate theories were used along with the MCST [43–45]. Ke
et al. [46, 47] employed p-version Ritz method and different quadrature method to solve
free vibration and bending, buckling problems of the rectangular and annular FG Mindlin
microplates based on MCST, respectively. Using modified couple stress theory, a refined
plate theory was utilised to predict the bending, buckling, and vibration behaviours of FG
microplates by He et al. [48] following the Navier’s approach. Most of these efforts followed
analytical approach which was able to solve some specific types of boundary condition.
A large proportion of the studies in small-scale FG structures employ the classical plate
theory (CPT) and the first-order shear deformation theory (FSDT). However, the classical plate theory (or the Kirchoff-Love theory), which neglects shear deformation, provides
acceptable solutions for thin plates, i.e. length-to-thickness ratios are quite large, only.
Meanwhile, the first-order shear deformation theory (or Reissner-Mindlin plate theory),
which accounts for transverse shear effects, is applicable for both thin and moderately thick
plates [49–51]. The shortcomings of the FSDT include inaccurate distribution of transverse
3

shear strain/stress and violation of traction free boundary conditions at top and bottom
surfaces. For this reason, shear correction factors are required to adjust the transverse shear
stress distribution. However, these factors vary from problems and one may find it difficult to choose an appropriate value. In order to avoid using shear correction factors, the
third-order shear deformation theory (TSDT) [52], higher-order shear deformation theory
(HSDT) [53], sinusoidal shear deformation theory (SSDT) [54], and refined plate theories
(RPT) [55] have been developed yielding more accurate and reliable results. The refined
plate theory is initially proposed by Senthilnathan et al. [55] by employing four unknowns
and is one variable less than that of TSDT. Shimpi et al. [56–58] then further developed
the RPT for isotropic and orthotropic plates by using only two variables. However, HSDT
and RPT require C1 -continuity of general displacements which cause significant challenge
to derive second derivative of deflection in the platform of finite element analysis (FEA)
where C0 elements are frequently used. In order to overcome the continuity issues, some C0
approximations [59] and Hermite interpolation functions with C1 elements [8] which involve
adding extra unknown variables of derivative of displacement are adopted.
Recently, a new numerical method so-called Isogeometric Analysis (IGA) which is able
to deal with C1 -continuity problem without using any additional variables or Hermite interpolation function has been introduced by Hughes and his co-workers [60]. This method
bridges the existing gap between the fields of finite element analysis and computer-aided design (CAD). The essential idea of the isogeometric analysis is that the basis functions which
are employed to exactly describe the geometry domain will also be used for approximations
of unknown fields. Since geometry is exactly modeled and the number of unknown terms is
not increased, it is expected that IGA would yield more accurate results with lower computational cost for HSDT and RPT problems in the comparison with regular FEA [61]. One
may find the guidance on computer implementation of IGA in the open literature [62–64].
IGA has been widely implemented in a number of linearly and non-linearly mechanical and
thermal problems such as static, free vibration, and buckling of laminated composite and
functionally graded plates with various plate theories including layerwise, FSDT, HSDT, and
RPT [53, 61, 65–73]. However, as far as author’s knowledge, there is no publication dealing
with analysis of small-scale plates based on modified couple stress theory and NURBS basis
functions.
In this study, in order to fill the aforementioned gap, the bending, free vibration and
buckling behaviours of FG microplates based on the modified couple stress theory and
four-variable refined plate and quasi-3D theories is investigated in conjunction with the
isogeometric approach. While the modified couple stress theory is employed to capture the
small-scale effects, the displacement fields of those microplates are expressed based on a novel
seventh-order RPT and quasi-3D theory. The mechanical behaviours of FG microplates are
then solved by the isogeometric analysis in which NURBS functions are simultaneously
used to exactly describe the geometry boundaries and construct the basis functions of the
approximations.
The outline of this study is expressed as follow. The next section presents theories which
are applicable for analysis of FG microplates including modified couple stress theory, RPT
and quasi-3D theory. In addition, a brief note on functionally graded material will also
4

be introduced in this section. Section 3 focuses on the isogeometric and NURBS-based
formulation of quasi-3D theory. The numerical examples which cover static bending, free
vibration and buckling analysis of rectangular and circular FG microplates with various
boundary conditions will be provided in Section 4. Finally, some conclusions will be drawn
in Section 5.
2. A novel theory for FG microplates
In this section, a brief review on the formulation of modified couple stress theory with
only one material length scale that accounts for size-dependent effects is firstly presented. It
is followed by the theoretical definition of functionally graded material of which microplates
are made. Displacement field of plates is then derived based on the four-variables refined
plate theory and quasi-3D plate theory where a novel seventh-order shear deformation theory
is proposed.
2.1. Modified couple stress theory
According to the modified couple stress theory which is proposed by Yang et al. [23], the
strain energy density w for linear isotropic material is a quadratic function of generalized
strains

1
(1)
w = λ (trεε)2 + µ ε : ε + l2 χ : χ ,
2
where λ and µ are Lame’s constants, µ is also known as shear modulus which is often denoted
as G, l represents material length scale parameter and the strain tensor ε and symmetric
curvature tensor χ are defined by
i
1h
T
ε=
∇u + (∇u) ,
(2a)
2
h
i
1
∇θ + (∇θ)T ,
(2b)
χ=
2
where u is displacement vector and the rotation vector θ is given by
1
θ = curl (u) .
(3)
2
The strain energy U restored in a deformed elastic body is then defined as
Z
Z
(4)
U = wdV = (σ : ε + m : χ) dV ,
V

V

where σ and m are the symmetric stress tensor and the deviatoric part of the couple stress
tensor, respectively. These components, σ and m, which are conjugated to the deformation
measures ε and χ, respectively, are given as
σ = λ (trεε) I + 2µεε,
m = 2µl2 χ,

(5a)
(5b)

where I is the identity matrix. Apparently, Eq. (5) shows there is only one material length
scale needed and the deviatoric couple stress tensor m is also symmetric, from which the
modified couple stress theory is formed.
5

2.2. Functionally graded material
The model configuration of the functional graded material which is made of metal and
ceramic is illustrated in Figure 1. There are several homogeneous models that are employed
to estimate the effective properties of the FGMs. According to the rule of mixtures, the
corresponding effective properties of the FGMs could be expressed as follow
Pe = Pm Vm + Pc Vc ,

(6)

where Pm and Pc are the material properties of metal and ceramic phases, respectively,
including Young’s modulus E, density ρ and Poisson’s ratio ν. Meanwhile, Vm and Vc
represent the volume fraction of metal and ceramic phases, respectively, which could be
defined as follow [8]
n

h
h
1 z
+
, Vm = 1 − Vc , − ≤ z ≤ ,
(7)
Vc (z) =
2 h
2
2
where n is the material index. This equation implies smooth variation in material properties
of FGM is guaranteed and governed by the material index n. As can be inferred from Eq.
(7), n = 0 leads to a FGM with fully homogeneous ceramic while the FGM becomes fully
metal as n goes toward ∞. Nevertheless, the rule of mixtures model fails to describe the
interaction between the material phases [74, 75]. Therefore, Mori-Tanaka scheme [76, 77]
is developed to fill this gap in estimation of FGM’s properties by integrating effective bulk
modulus Ke and shear modulus Ge which are given by
Vc
Ke − Km
=
,
m
Kc − Km
1 + Vm KKmc+−K
4
Gm

Vc
Ge − Gm
=
,
m
Gc − Gm
1 + Vm GGcm−G
+f1

3

where
f1 =

Gm (9Km + 8Gm )
.
6 (Km + 2Gm )

(8)

(9)

The effective Young’s modulus Ee and Poisson’s ratio νe are then defined as
Ee =

9Ke Ge
,
3Ke + Ge

νe =

3Ke − 2Ge
.
2 (3Ke + Ge )

(10)

The variations in effective Young’s modulus of Al/Al2 O3 estimated by the rule of mixtures
and Mori-Tanaka scheme are depicted in Figure 2. As can be seen, corresponding to the
material index n, the effective property of FGM is varied continuously from the metal-rich
surface at the bottom to the ceramic-rich one at the top of the plate.
2.3. A novel seventh-order shear deformation plate theory
With regard to the plate theories, the third-order shear deformation model initially
proposed by Reddy [52] is widely considered as a reliable theory in which no shear correction
6

factor is required. In Reddy’s theory, the displacement field, for z ∈ [−h/2; h/2], is defined
as
u (x, y, z) = u0 (x, y) + zβx (x, y) + g (z) (βx (x, y) + w,x (x, y)) ,
v (x, y, z) = v0 (x, y) + zβy (x, y) + g (z) (βy (x, y) + w,y (x, y)) ,
w (x, y, z) = w0 (x, y) ,

(11a)
(11b)
(11c)

where comma notation represents derivative, g(z) = −4z 3 /(3h2 ) and the variables u0 =
[u0 v0 ]T , w0 , and β = [βx βy ]T are the membrane displacements, the transverse deflection
of the mid-plane surface, and the rotations, respectively. By making further assumptions,
w0 = wb + ws , βx = −wb,x , βy = −wb,y , to the Reddy’s theory which contains five unknowns,
Senthilnathan [55] initially proposes the four-variable refined plate theory which could be
expressed in the generalized form as follows
u (x, y, z) = u0 (x, y) − zwb,x (x, y) + g (z) ws,x (x, y) ,
v (x, y, z) = v0 (x, y) − zwb,y (x, y) + g (z) ws,y (x, y) ,
w (x, y, z) = wb (x, y) + ws (x, y) ,

(12a)
(12b)
(12c)

where wb and ws represent bending and shear components of transverse displacement,
respectively. The function g(z) = f (z) − z is employed to describe the distribution of
transverse shear strains and stresses through the plates’ thickness. It is necessary to have the
first derivative of f (z) satisfies the tangential zero value at z = ±h/2 such that the tractionfree condition at top and bottom surfaces is meet. Consequently, the shear correction factor
is no longer required for higher-order shear deformation theory and refined plate theory.
It should be noted that the higher-order shear deformation theory and refined plate
theory both fail to capture the thickness stretching effect of normal deformation (εz 6= 0)
due to the constant deflection through the plate thickness which can be inferred from Eq.
(12c). In order to bypass this shortcoming, a number of theories which consider the thickness
stretching effect have been developed [78–80]. Zenkour [81, 82] initially proposes the fourvariable quasi-3D plate theory accounting for both transverse shear and normal deformations
which could be alternatively expressed as follows
u (x, y, z) = u0 (x, y) − zwb,x (x, y) + f (z) ws,x (x, y) ,
v (x, y, z) = v0 (x, y) − zwb,y (x, y) + f (z) ws,y (x, y) ,
w (x, y, z) = wb (x, y) + φ (z) ws (x, y) .

(13a)
(13b)
(13c)

As can be observed, this quasi-3D model has a similar form to that of four-variable
refined plate theory shown in Eq. (12). Indeed, the displacement field based on the refined
plate theory could be readily obtained by degrading those of quasi-3D theory in which f (z)
and φ(z) are replaced by g(z) and 1, respectively.
A number of distribution functions, f (z) and φ(z), are available to analyse the FGM
plates based on higher-order shear deformation theory, refined plate theory, and quasi-3D
7

plate theory. One may find the general framework to construct such polynomial functions
in the recent work of Nguyen et al. [83]. In this study, a novel seventh-order function of
f (z) and corresponding function φ(z) are proposed for the four-variable refined plate theory
and quasi-3D theory. Those functions are presented in Table 1 along with other’s existing
in the literature.
According to the displacement field and the strain-displacement relation, which are presented in Eqs. (13) and (2a), respectively, the following strain expressions could be obtained
as
ε = ε 0 + zκb + f (z) κs ,
γ = [f 0 (z) + φ (z)] ε s ,

(14a)
(14b)

where







εx
u0,x
wb,xx
 , κb = −  wb,yy  ,
ε =  εy  , ε 0 =  v0,y
γxy
u0,y + v0,x
2wb,xy




ws,x
γxz
, εz = φ0 (z) ws .
, εs =
γ=
ws,y
γyz




ws,xx
κs =  ws,yy  ,
2ws,xy

(15)

Using Eqs. (13), (3), and (2b), rotation vector and curvature vector are respectively expressed by




θx
2wb,y − (f 0 − φ) ws,y
1
(16a)
θ =  θy  =  −2wb,x + (f 0 − φ) ws,x  ,
2
θz
v0,x − u0,y
 
  0
 


χb
χb0
f χb1
φχb3
χ =  χs  =  χs0  +  f 00 χs2  +  φ0 χs4  ,
(16b)
χzz
0
0
0
where









χxx
wb,xy
−2ws,xy
2ws,xy
1
 , χb1 = 1  2ws,xy
 , χb3 = 1  −2ws,xy
,
χb =  χyy  , χb0 =  −wb,xy
2
4
4
χxy
wb,yy − wb,xx
−ws,yy + ws,xx
ws,yy − ws,xx








1 −ws,y
1 ws,y
1 v0,xx − u0,xy
χxz
χs =
, χs0 =
, χs2 =
, χs4 =
.
χyz
4 v0,xy − u0,yy
4 ws,x
4 −ws,x
(17)
According to Eq. (5), the constitutive relations for classical and modified couple stress

8

theories could be presented in explicit form as follows

 
σ
Q11 Q12 Q13 0
0
0


x







σ
Q
Q
0
0
0


y 
22
23


 

σz
Q33 0
0
0
=

σ
Q
0
0

xy 
66








τ
Q55 0



 xz 

τyz
sym.
Q44


















εx
εy
εz
εxy
γxz
γyz










,

(18a)









mij = 2Ge l2 χij ,

(18b)

where, for the proposed quasi-3D theory (εz 6= 0), Qij are the three-dimensional elastic
constant which could be calculated as
(1 − νe ) Ee
,
(1 − 2νe ) (1 + νe )
νe E e
=
,
(1 − 2νe ) (1 + νe )
Ee
,
=
2 (1 + νe )

Q11 = Q22 = Q33 =
Q12 = Q13 = Q23
Q44 = Q55 = Q66

(19)

meanwhile, for the proposed refined plate theory (εz = 0), Qij are reduced plane-stress
elastic constant and are expressed as
Ee
,
1 − νe2
Ee νe
=
,
1 − νe2

Q11 = Q22 =
Q12 = Q21

Q44 = Q55 = Q66 =

(20)

Ee
,
2 (1 + νe )

Ee
.
2(1 + νe )
In this study, the weak form of the static bending, vibration, and buckling problems
are derived using the Hamilton’s principle and weak formulation. One could find those
straightforward procedures in the literature [61, 84, 85]. Firstly, the weak form of the static
bending of the couple-stress-based microplates subjected to transverse load loading q0 could
be expressed in the following compact form
Z
Z
Z
Z
b c
T
b
T
s
c T
δεεb D ε b dΩ + δεεs D ε s dΩ + δ (χb ) D c χb dΩ + δ (χcs )T Dsc χcs dΩ
and the shear modulus Ge =

Ω

Ω

Ω

Ω

Z 
=
Ω

9

 

,
h
δwb + φ
δws q0 dΩ
2

(21)

where the strain tensors and matrices in the first two terms in Eq. (21) related to classical
elastic theory are




ε0
A
B
E
X


 κb 
 B D F Yb 
ws,x
b




εb = 
,
ε
=
,
D
=
(22)
s
 E F H Ys  ,
κs 
ws,y
ws
X Yb Ys Z33
in which the material matrices are calculated by
Zh/2
(A, B, D, E, F, H) =




1, z, z 2 , f (z) , zf (z) , f 2 (z) Q̄dz,

(23a)

−h/2

X, Yb , Y


s

Zh/2
=

[φ0 (z) , zφ0 (z) , f (z) φ0 (z)] Q̃dz,

(23b)

−h/2

Zh/2
Z33 =

2

[φ0 (z)] Q33 dz,

(23c)

2

(23d)

−h/2

Zh/2

Ds =

[f 0 (z) + φ (z)] Q̂dz,

−h/2




Q11 Q12 0
Q̄ =  Q21 Q22 0  ,
0
0 Q66




Q13
Q̃ =  Q23  ,
0

and where the curvature tensors and matrices in the
representing the couple stress theory are




 c
χb0
χs0
A Bc
χcb =  χb1  , χcs =  χs2  , Dbc =  Bc Dc
χb3
χs4
Ec Fc


Q̂ =

Q44 0
0 Q55


,

(23e)

third and fourth terms in Eq. (21)

Ec
Fc  ,
Hc




X c Y c Tc
Dsc =  Yc Zc Vc  ,
Tc Vc Wc
(24)

in which the material matrices could be defined as
(Ac , Bc , Dc , Ec , Fc , Hc ) =

Zh/2 


2
1, f 0 (z) , [f 0 (z)] , φ (z) , f 0 (z) φ (z) , [φ (z)]2 Ḡdz, (25a)

−h/2

(Xc , Yc , Zc , Tc , Vc , Wc ) =

Zh/2 

2

2

1, f 00 (z) , [f 00 (z)] , φ0 (z) , f 00 (z) φ0 (z) , [φ0 (z)]



Ĝdz,

−h/2

(25b)
10

where




1 0 0
Ḡ = 2Ge l2  0 1 0  ,
0 0 1

Ĝ = 2Ge l

2



1 0
0 1


.

(26)

The weak form of the free vibration of the couple-stress-based microplates is briefly
expressed as
Z
Z
Z
Z
Z
c T
s c
c T
b c
T
s
T
b
¨
δε b D ε b dΩ + δεεs D ε s dΩ + δ (χb ) Dc χb dΩ + δ (χs ) Dc χs dΩ = δ ũT m̃ũdΩ,
where ũ = [u0 − wb,x ws,x v0 − wb,y ws,y wb ws
defined by




I0 0 0
I1 I2 I4
m̃ =  0 I0 0  in which I0 =  I2 I3 I5  ,
0 0 I1
I4 I5 I6
Zh/2
(I1 , I2 , I3 , I4 , I5 , I6 , I7 , I8 ) =

Ω

Ω

Ω

Ω

Ω

(27)
0] , and the mass matrix m̃ is
T




I1 I7 0
I1 =  I7 I8 0  ,
0 0 0

(28a)



ρ 1, z, z 2 , f (z) , zf (z) , f 2 (z) , φ (z) , φ2 (z) dz.

(28b)

−h/2

For buckling analysis, the weak form of the couple-stress-based microplates subjected to
in-plane loading is of the form
Z
Z
Z
Z
T
b
T
s
c T
b c
δεεb D εb dΩ + δεεs D εs dΩ + δ (χb ) Dc χb dΩ + δ (χcs )T Dsc χcs dΩ
Ω

Ω

Ω

Z
+

Ω

,

(29)

∇T δ [wb + φ (0) ws ] N0 ∇ [wb + φ (0) ws ] dΩ = 0

Ω
0
Nx0 Nxy
where ∇T = [∂/∂x ∂/∂y] and N0 =
0
Nxy Ny0
and matrix of pre-buckling loads, respectively.
T




are the transpose of gradient operator

3. FG microplate formulation based on NURBS basis functions
In this section, briefing review on NURBS which are served as the basis functions of
isogeometric analysis will be presented. It is followed by a novel NURBS-based formulation
for couple-stress microplate bending, free vibration and buckling that are relied on refined
plate theory and quasi-3D theory.
3.1. B-splines and NURBS basis functions
The starting point of the NURBS basis functions is non-decreasing knot vector Ξ =
{ξ1 , ξ2 , ..., ξn+p+1 } where the ith knot ξi ∈ R, n represents the number of basis functions, and
11

p denotes the polynomial order. The knot vectors could be either uniform if the knots are
equally spaced in the parameter space or open if its first and last knot’s values are repeated
p + 1 times. The knot spans which are bounded by knots defines element domain. The
B-spline basis functions that are constructed by the Cox-de Boor recursion formula, starting
with the zeroth order (p = 0) are defined by [60, 86]

1 if ξi ≤ ξ < ξi+1 ,
Ni,0 (ξ) =
(30a)
0 otherwise,
ξi+p+1 − ξ
ξ − ξi
Ni,p−1 (ξ) +
Ni+1,p−1 (ξ) , for p ≥ 1,
(30b)
Ni,p (ξ) =
ξi+p − ξi
ξi+p+1 − ξi+1
in which the fraction 0/0 is defined as zero. While the basis functions are smooth, e.g. C ∞
continuity, within this domain, they are C p−m continuity across the knots, where m is the
multiplicity of the knot. Therefore, for p ≥ 2, the basis functions would be of C1 continuity
at single knot and knot span as well. Two-dimensional B-splines are obtained by introducing
the second knot vector H = {η1 , η2 , ..., ηm+q+1 } and using tensor product of Ξ and H in the
parametric dimensions yielding
NA (ξ, η) = Ni,p (ξ) Mj,q (η) .

(31)

For the illustration purposes, Fig. 3 depicts the one- and two-dimensional B-spline basis
functions which are generated from the knot vector Ξ = {0, 0, 0, 0, 1/5, 2/5, 3/5, 3/5, 3/5, 4/5, 1, 1, 1, 1}
and its combination with the knot vector H = {0, 0, 0, 1/4, 1/2, 3/4, 1, 1, 1}, respectively.
The non-uniform rational B-splines (NURBS) basis functions are then further defined by
providing additional weight ζA to each control point given by [86]
RA (ξ, η) =

NA ζA
m×n
P

.

(32)

NÂ (ξ, η) ζÂ

Â

It is noted that the B-spline basis function is served as a special case of the NURBS function.
Indeed, if all the individual weights corresponding the control points are assigned an equal
constant, the NURBS function yields B-spline function.
3.2. A novel NURBS-based formulation of modified couple stress theory
By using the NURBS basis functions, the displacement variables u of a microplate could
be approximately calculated as follows
h

u (ξ, η) =

n×m
X

RA (ξ, η)qA ,

(33)

A=1

where n × m is the number of basis functions and qA = [u0A v0A wbA wsA ]T denotes the
vector of degrees of freedom associated with the control point A. By substituting the variable
12

approximations Eq. (33) into the strain-displacement relations Eq. (15), the in-plane and
shear strains could be obtained as follows
 T
ε0

κTb

κTs

εTs

T

=

n×m
Xh

T
(Bm
A)

Bb1
A

T

Bb2
A

T

(BsA )T

iT

qA ,

(34)

A=1

where




RA,x
0
0 0
0
m
b1



0
RA,y 0 0 , BA = − 0
BA =
RA,y RA,x 0 0
0



0 0 0 RA,xx
0 0
b2
s
BA =  0 0 0 RA,yy  , BA =
0 0
0 0 0 2RA,xy


0 RA,xx 0
0 RA,yy 0  ,
0 2RA,xy 0

0 RA,x
,
0 RA,y

(35)

and the curvatures could be obtained by substituting Eq. (33) into Eq. (17) :
 T
χb0 χTb1 χTb3
n×m
X   T
B̃b0
=
A

T
χTs0 χTs2 χTs4
 T  T
B̃b1
B̃b3
A
A



B̃s0
A

T



B̃s2
A

T



B̃s4
A

T T

qA ,

(36)

A=1

where





0 0 0
−2RA,xy
1
,
0 0 0
2RA,xy
B̃b1
A =
4
0 0 0 (RA,xx − RA,yy )


1 −RA,xy RA,xx 0 0
(37)
s0
B̃A =
,
−R
R
0
0
4
A,yy
A,xy

0 0 0 RA,y
.
0 0 0 −RA,x


0 0
2RA,xy
0
1
0 0
−2RA,xy
0 ,
B̃b0
A =
2
0 0 (−RA,xx + RA,yy ) 0


0 0 0
2RA,xy
1
,
0 0 0
−2RA,xy
B̃b3
A =
4
0 0 0 (−RA,xx + RA,yy )



1
1 0 0 0 −RA,y
s2
s4
, B̃A =
B̃A =
4 0 0 0 RA,x
4

Substituting Eqs. (34) and (36) into Eqs. (21), (27), and (29), the matrix form of
the global equilibrium equations for static bending, free vibration, and buckling could be
respectively presented by
Kq = F,

K − ω M q = 0,
(K − λcr Kg ) q = 0,
2

(38a)
(38b)
(38c)

where the global stiffness matrix K is the summation of the stiffness matrices corresponding
to the classical theory Ks and the couple stress theory Kc , i.e. K = Ks +Kc . These matrices
13

are calculated as follow
 m T 
 B 

Z 
 Bb1  

Ks = 
 Bb2  

 z 

Ω
B

T 
Z  B̃b0 


Kc =  B̃b1
 b3 
B̃
Ω

 m 

A B E X
B 



 b1 

B D F Yb 
B
s T
s s

+
(B
)
D
B
(39a)
 dΩ,
E F H Ys  
Bb2 


 z 
X Yb Ys Z33
B

 s0 T  c

  s0 

b0
c
c
c
c
c
X Y T
A B E
 B̃ 
 B̃   B̃ 
 Yc Zc Vc  B̃s2 
Bc Dc Fc  B̃b1
+
B̃s2
 dΩ,
 s4 
 b3   s4 
c
c
c
c
c
c
T V W
E F H
B̃
B̃
B̃
(39b)


in which BzA = 0 0 0 RA . The load vector F is given by
Z
(40)
F = q0 RdΩ,
Ω

where R =



0 0 RA φ

h
2



RA

T

, the global mass matrix is computed by
Z
M = R̃T m̃R̃dΩ,

(41)

Ω

in which




RA 0
0
0
 R1 
R2
0 ,
, R1 =  0 0 −RA,x
R̃ =


R3
0 0
0
RA,x




0 RA
0
0
0 0 RA 0
0  , R3 =  0 0 0 R A  ,
R2 =  0 0 −RA,y
0 0
0
RA,y
0 0 0
0
the geometric stiffness matrix is given as
Z
Kg = (Bg )T N0 Bg dΩ,

(42)

(43)

Ω

where
g

B =



0 0 RA,x φ (0) RA,x
0 0 RA,y φ (0) RA,y


,

(44)

and ω and λcr represent the natural frequency and the critical buckling value, respectively.
As can be observed from Eq. (39), by introducing the distribution function f (z), the fourvariable refined plate theory and quasi-3D theory do not require any shear correction factor,
which is usually needed if the first-order shear deformation theory is applied, to describe
the transverse shear stresses satisfying traction-free conditions. In addition, the expressions
14

of B and B̃ matrices in Eqs. (35) and (37), respectively, show the employment of the
second-order derivatives of the approximation functions RA . Consequently, C1 -continuity
approximations are required. This requirement may cause severe issues in finite element
analysis which could only be solved by increasing the degrees of freedom resulting in higher
number of variables and computational cost. However, within the platform of isogeometric
analysis in which NURBS basis functions are employed, the C1 -continuity requirement could
be naturally satisfied for p ≥ 2 since the basis functions are Cp−1 continuous across the knot
spans, i.e. elements.
4. Numerical examples and discussion
In this section, the convergence and verification studies are firstly conducted to demonstrate the validity and accuracy of the novel approaches presented in Section 2 and 3. In
order to illustrate the efficiency of the IGA approach in dealing with modified couple stress
theory, this section is then continued by the computational analysis of FG rectangular and
circular microplates with various types of boundary condition for static bending, free vibration and buckling problems. In these investigations, the FG microplates made of mixtures
of metal and ceramic whose material properties presented in Table 2 are used. Throughout
the numerical examples, unless otherwise specified, the material length scale l is chosen as
17.6 × 10−6 m which was suggested by Lam et al. [20]. There are two types of boundary
condition considered as follows
Simply supported (S)
v0 = wb = ws = 0 at x = 0, a
u0 = wb = ws = 0 at y = 0, b
Clamped (C)
u0 = v0 = wb = ws = 0 and wb,x = wb,y = wb,z = ws,x = ws,y = ws,z = 0
It should be noted that, within the IGA platform, while the homogeneous boundary
conditions corresponded to displacement itself, e.g. u0 , v0 , wb , ws , are easily treated in a
similar way to the traditional finite element method, those require first derivative of displacement, e.g. wb,x , wb,y , wb,z , ws,x , ws,y , ws,z , could be done by assigning zero values for
all displacements of control points which are directly related to clamped edges and their
adjacent points [85, 87].
4.1. Convergence and verification studies
In order to evaluate the convergence and reliability of the novel shape functions and
IGA approach proposed in Section 2 and 3, the MCST-based size-dependent analysis of
homogeneous fully simply-supported (SSSS) square microplate which is shown in Fig. 5 is
conducted using RPT model. Moderately thick plates (a/h = 20) with four different values
of material length scale ratio (l/h = 0, 0.2, 0.6, 1) are investigated. For each case, eight
different finite element meshes are analysed to study the convergence rate of the proposed
IGA approach. As can be seen from Table 3, while fast convergence of the analysis for
polynomial order p = 3 and p = 4 is archived, solutions using p = 2 experience relatively
slower convergence rate toward analytical solutions reported by Thai and Kim [43]. Fig. 4
15

presents the convergence study of non-dimensional central deflection of homogeneous square
microplates with respect to the analytical solutions [43] for the case of l/h = 0.2. Based
on the convergence study, the cubic (p = 3) NURBS element mesh of 11 × 11 is relatively
sufficient for all analysis cases. Therefore, this mesh whose geometry shown in Fig. 5 will
be used throughout the next examples unless otherwise specified.
Further investigation on the accuracy of the proposed method is conducted using FG
plates made of alumina and aluminum (Al/Al2 O3 ). In this case, without considering couple
stress effects, proposed RPT (εz = 0) and quasi-3D (εz 6= 0) theories are applied in analysis
of SSSS square plates using the rule of mixtures scheme. The plates are subjected
to
uni πy

 πx 
sin
,
formly and sinusoidally distributed loads which are defined as q0 and q0 sin
a
a
respectively. As can be observed in Table 4, the present results are in good agreement with
those available in published works using various 2-D and quasi-3D theories. Apparently,
these investigations have confirmed the validity and reliability of the present study.
4.2. Static bending analysis
In this section, the static bending analysis of FG microplates based on modified couple
stress theory will be investigated. The SSSS square microplates whose materials follow
the rule of mixtures. The aspect ratio a/h, material length scale ratio l/h, and material
index n are taken into account. Table 5 presents the comparison of non-dimensional central
deflection of SSSS square plate with those of Thai and Kim [43]. While the results generated
from the proposed RPT theory are in very good agreement, quasi-3D theory yields slightly
different responses in term of displacement. This is due to the consideration of the thickness
stretching effect in the quasi-3D theory.
The bending responses of CCCC square Al/Al2 O3 microplates under sinusoidally and
uniformly distributed loads are further studied and presented in Table 6. It is noted that
since there is no study on the static behaviours of CCCC microplates using MCST are
reported in the literature, the results are compared with those generated from Reddy’s
HSDT model [52] using the proposed IGA approach. As can be seen, the results based on
Reddy’s model yields excellent agreement with proposed RPT-based solutions. The effects
of material index n and material length scale l on the central displacement of CCCC square
Al/Al2 O3 plate are depicted in Fig. 6 in which data are generated by the proposed RPT
and quasi-3D theories for three different ratios l/h of 0, 0.4 and 1.0. As can be observed, the
increase of material index n leads to a rise of plates’ central deflection due to the decrease in
plates’ stiffness. On the contrary, the growth of material length scale ratio l/h is followed by
a decline of displacement due to the increase in magnitude of stiffness matrix. In other word,
for specific material length scale l, the thinner the microplates are, the higher plates’ stiffness
would be. It can be further observed that the discrepancy in term of central deflection by
the proposed RPT and quasi-3D is significantly decreased as l/h increase and vanishes when
l/h = 1.0. Fig. 7 depicts the deformed configurations of the Al/Al2 O3 square microplates
with various boundary conditions subjected to sinusoidally distributed load in which a/h, l/h
and n are equal to 5, 0.4 and 10, respectively. It should be noted that the deformed shape
of the microplates are scaled up for illustration purpose.
16

4.3. Free vibration analysis
In this part, the free vibration analysis of FG microplates based on modified couple
stress theory will be discussed. The proposed Quasi-3D model is initially tested for normalscale SSSS Al/ZrO2 -1 plates with various theories taking into account the normal shear
deformation. It can be seen that the comparison results which are given in Table 7 agree
well with those of other published works. The proposed RPT and quasi-3D models using
IGA are further tested for homogeneous square microplates. The results which are shown in
Table 8 are compared with analytical solutions generated from CPT by Yin et al. [35] and
TSDT by Thai and Kim [43] . As can be observed, while the proposed RPT model yields
slight discrepancy with respect to Yin et al.’s [35] due to their CPT assumption neglecting
shear deformations, it shows excellent agreement with Thai and Kim’s [43], especially as
plates become thinner, i.e. a/h ratio is relatively large. Due to the consideration of thickness
stretching effects, the proposed quasi-3D model gives slightly different results in comparison
with other results.
Table 9 presents non-dimensional natural frequency of SSSS Al/Al2 O3 square microplates.
The results are compared with those of Thai and Kim [43] in which analytical approach
based on TSDT model is employed. Thick (a/h = 5), moderately thick (a/h = 20) and thin
(a/h = 100) microplates are considered. The results reveal good agreement between RPT
model and TSDT [43], especially when the material length scale ratio l/h is small, e.g. 0
or 0.2. On the contrary, the discrepancy becomes larger as l/h gets closer to 1. However,
this phenomenon just happens for thick plates and tends to be vague as the plates become
thinner. Meanwhile, the quasi-3D gives slightly different results compared to that of RPT
model. General observation from Table 9 reveals that the higher material length scale ratio
is chosen, the bigger natural frequencies of the plates could be obtained as a result of the
increase in plate’s stiffness. Fig. 8 provides a closer look at the effects of material index
n and material length scale l on the plates’ natural frequencies which are computed using proposed RPT and quasi-3D models. Similar to the previous case of bending analysis,
the plate’s stiffness decreases as a result of rising in material index n and falling in material length scale ratio l/h which leads to the decrease of natural frequency of the plates.
In addition, the discrepancy in term of frequency results predicted by the proposed RPT
and quasi-3D becomes less significant as l/h gets bigger and is almost zero difference when
l/h = 1.0. The first six natural frequencies of Al/Al2 O3 square microplates with different
types of boundary conditions are given in Table 10 in which the results are generated for
n = 1 and l/h = 0.2. The present quasi-3D results show good agreement with those of
Zenkour’s quasi-3D model [81] using the proposed IGA approach. The first six mode shapes
corresponding to the quasi-3D vibration analysis of CCCC microplates with a/h = 10 are
presented in Fig. 9.
Next, free vibration of circular plates whose geometry configuration and mesh are shown
in Fig. 10 will be investigated. Since there is no publication on the vibration behaviours of
FG circular microplates based on the modified couple stress theory is found, the investigation
of circulars could be served as benchmark examples. Table 11 presents the fundamental
natural frequencies which are firstly tested for homogeneous plates. Natural frequencies of
plates without considering size-dependent effects (l = 0) are compared with results reported
17

by Mohammandi et al. [88] and Nguyen et al. [85]. The results show very good agreement
between the theories, especially proposed RPT and Nguyen et al.’s [85] which also uses
another polynomial-based RPT model. Table 12 presents the first six natural frequencies
of Al/Al2 O3 circular plates with simple and clamped supports. The plates’ thickness h and
material index n are set as 0.2R and 1, respectively. As previous case of square microplates,
Zenkour’s quasi-3D model [81] using proposed IGA approach is added for reference purpose
since there is no study can be found for this analysis problem. The six fundamental mode
shapes of clamped plates with l/h = 0.2 using quasi-3D model are printed in Fig. 11.
4.4. Buckling analysis
In this stage, buckling behaviours of square and circular FG microplates will be discussed.
In order to verify the proposed method and models in dealing with buckling analysis, first
attempt is to calculate the critical buckling load of SSSS FG microplates bearing biaxial
load. The results are compared with analytical solutions based on CPT and FSDT by Thai
and Choi [44] and RPT by He et al. [48] for which material properties are E1 = 14.4GPa,
ρ1 = 12.2 × 103 kg/m3 , E2 = 1.44GPa, ρ2 = 1.22 × 103 kg/m3 , ν1 = ν2 = 0.38. As can be seen
in Table 13, although there is discrepancy for the case of thick microplates (a/h = 5), the
results predicted by proposed RPT theory are in good agreement with those calculated by
FSDT [44] and RPT [48] as plates become thinner. On the other hand, while CPT-based
solutions are of remarkable difference to that of FDST and RPTs due to the ignorance of
shear deformations, especially for thick plates, the proposed quasi-3D theory which takes into
account normal deformation yields similar results with respect to those of shear deformable
theories of FSDT and RPTs.
Table 14 presents the biaxial buckling analysis results of Al/Al2 O3 square microplates.
The results are calculated based on the proposed RPT and quasi-3D theories which could
be served as benchmark examples for future references since no report in the similar case
is found. The first six non-dimensional biaxial buckling loads of Al/Al2 O3 square plates
are reported in Table 15 for n = 10 and l/h = 0.2. While the results generated from the
proposed quasi-3D and IGA-based Zenkour’s quasi-3D theories [81] are relatively close to
each other for both type of boundary conditions, the RPT’s show remarkable discrepancy
to the other theories, especially for the CCCC plates. This is due to the consideration of
normal deformation of the quasi-3D theories. For a/h = 5, six buckling mode shapes of
CCCC plates based on the proposed quasi-3D theory which are scaled up for illustration
purposes are presented in Fig. 12.
Finally, this section ends up with a number of investigations on buckling behaviours of
circular FG microplates. Table 16 presents the results of critical buckling load of CCCC
Al/ZrO2 -2 circular plates without considering couple stress effects. It should be noted that,
for this particular attempt of comparison purpose, the material volumes are defined as
Vm = (0.5 − z/h)n and Vc = 1 − Vm [61, 68, 89, 90]. The comparison reveals that the results
generated from proposed RPT are in good agreement with those of other shear deformation
theories even with relatively thick and thick plates. The proposed quasi-3D yields slightly
different results in all cases.
18

Table 17 presents the non-dimensional critical buckling load of Al/Al2 O3 circular microplates with various boundary conditions based on the proposed RPT and quasi-3D theories. The difference between the two theories is relatively small for both case of boundary
conditions. The effects of material index n and material length scale ratio l/h on the critical buckling load of simply and clamped support of Al/Al2 O3 circular microplates with
h/R = 0.2 are illustrated in Fig. 13. The first six buckling loads of Al/Al2 O3 circular microplates with various aspects ratios for n = 1 and l/h = 0.6 along with the results generated
from Zenkour’s quasi-3D theory [81] using proposed IGA approach are reported in Table 18.
Six scaled buckling mode shapes corresponding to simply-supported circular microplates for
h/R = 0.2 based on proposed quasi-3D theory are presented in Fig. 14.
5. Conclusions
In this study, a novel computational approach, which employs the modified couple stress
theory, the four-variable refined plate theory and quasi-3D, and the NURBS-based isogemetric analysis, has been presented to investigate the static bending, free vibration, and
buckling behaviours of functionally graded microplates with various geometric domains and
boundary conditions. Within this proposed approach, the mathematical model governing
the behaviours of the plates is constructed based on the modified couple stress theory with
only one material length scale which efficiently taken into account the size-dependent effects
of small-scale structures and the novel seventh-order refined plate theory and quasi-3D theory with only four unknowns. This proposed quasi-3D not only considers shear deformations
but also is able to well capture the thickness stretching effect which is neglected by other
classical, higher-order shear deformation, and refined plate theories. The NURBS-based isogeometric analysis is incorporated within the computational process to exactly describe the
geometry an approximately construct the unknown field in which higher-order derivatives
and continuity elements are readily satisfied. A large number of investigations including some
benchmark examples have confirmed the validity and efficiency of the proposed approach.
The obtained results also reveal that the increase of the material length scale parameter
rises the microplates’ stiffness which results in the decline in the central displacement and
increase in the natural frequency and buckling load. On the contrary, the growth of the
material index leads to completely reverse trends in which microplates’ stiffness are fallen.
Meanwhile, the numerical results that are generated by the proposed refined plate theory
and quasi-3D theory are slightly different due to the consideration of the thickness stretching
effect. However, these discrepancies become gradually vague as the microplates get thinner.
References
[1] L. C. Trinh, T. P. Vo, A. I. Osofero, J. Lee, Fundamental frequency analysis of functionally graded
sandwich beams based on the state space approach, Composite Structuresdoi:10.1016/j.compstruct.
2015.11.010.
[2] T. P. Vo, H.-T. Thai, T.-K. Nguyen, A. Maheri, J. Lee, Finite element model for vibration and buckling of functionally graded sandwich beams based on a refined shear deformation theory, Engineering
Structures 64 (2014) 12–22. doi:10.1016/j.engstruct.2014.01.029.

19

[3] T. P. Vo, H.-T. Thai, T.-K. Nguyen, F. Inam, Static and vibration analysis of functionally graded
beams using refined shear deformation theory, Meccanica 49 (1) (2013) 155–168. doi:10.1007/
s11012-013-9780-1.
[4] M. Asghari, M. T. Ahmadian, M. H. Kahrobaiyan, M. Rahaeifard, On the size-dependent behavior
of functionally graded micro-beams, Materials & Design 31 (5) (2010) 2324–2329. doi:10.1016/j.
matdes.2009.12.006.
[5] V.-H. Nguyen, T.-K. Nguyen, H.-T. Thai, T. P. Vo, A new inverse trigonometric shear deformation
theory for isotropic and functionally graded sandwich plates, Composites Part B: Engineering 66 (2014)
233–246. doi:10.1016/j.compositesb.2014.05.012.
[6] H.-T. Thai, T. P. Vo, A new sinusoidal shear deformation theory for bending, buckling, and vibration
of functionally graded plates, Applied Mathematical Modelling 37 (5) (2013) 3269–3281. doi:10.1016/
j.apm.2012.08.008.
[7] A. M. Zenkour, Generalized shear deformation theory for bending analysis of functionally graded plates,
Applied Mathematical Modelling 30 (1) (2006) 67 – 84. doi:http://dx.doi.org/10.1016/j.apm.
2005.03.009.
[8] J. N. Reddy, Analysis of functionally graded plates, International Journal for Numerical Methods in
Engineering 47 (1-3) (2000) 663–684. doi:10.1002/(SICI)1097-0207(20000110/30)47:1/3<663::
AID-NME787>3.0.CO;2-8.
[9] J. L. Mantari, Refined and generalized hybrid type quasi-3d shear deformation theory for the bending
analysis of functionally graded shells, Composites Part B: Engineering 83 (2015) 142–152. doi:10.
1016/j.compositesb.2015.08.048.
[10] J. Torabi, Y. Kiani, M. R. Eslami, Linear thermal buckling analysis of truncated hybrid FGM conical
shells, Composites Part B: Engineering 50 (2013) 265–272. doi:10.1016/j.compositesb.2013.02.
025.
[11] F. Tornabene, Free vibration analysis of functionally graded conical, cylindrical shell and annular plate
structures with a four-parameter power-law distribution, Computer Methods in Applied Mechanics and
Engineering 198 (3740) (2009) 2911–2935. doi:10.1016/j.cma.2009.04.011.
[12] Y. Fu, H. Du, W. Huang, S. Zhang, M. Hu, TiNi-based thin films in MEMS applications: a review,
Sensors and Actuators A: Physical 112 (23) (2004) 395–408. doi:10.1016/j.sna.2004.02.019.
[13] Z. Lee, C. Ophus, L. M. Fischer, N. Nelson-Fitzpatrick, K. L. Westra, S. Evoy, V. Radmilovic, U. Dahmen, D. Mitlin, Metallic NEMS components fabricated from nanocomposite AlMo films, Nanotechnology 17 (12) (2006) 3063. doi:10.1088/0957-4484/17/12/042.
[14] R. H. Baughman, C. Cui, A. A. Zakhidov, Z. Iqbal, J. N. Barisci, G. M. Spinks, G. G. Wallace, A. Mazzoldi, D. D. Rossi, A. G. Rinzler, O. Jaschinski, S. Roth, M. Kertesz, Carbon Nanotube Actuators,
Science 284 (5418) (1999) 1340–1344. doi:10.1126/science.284.5418.1340.
[15] K.-T. Lau, H.-Y. Cheung, J. Lu, Y.-S. Yin, D. Hui, H.-L. Li, Carbon Nanotubes for Space and BioEngineering Applications, Journal of Computational and Theoretical Nanoscience 5 (1) (2008) 23–35.
doi:10.1166/jctn.2008.003.
[16] S. J. V. Frankland, A. Caglar, D. W. Brenner, M. Griebel, Molecular Simulation of the Influence of
Chemical Cross-Links on the Shear Strength of Carbon NanotubePolymer Interfaces, The Journal of
Physical Chemistry B 106 (12) (2002) 3046–3048. doi:10.1021/jp015591+.
[17] K. M. Liew, X. Q. He, C. H. Wong, On the study of elastic and plastic properties of multi-walled
carbon nanotubes under axial tension using molecular dynamics simulation, Acta Materialia 52 (9)
(2004) 2521–2527. doi:10.1016/j.actamat.2004.01.043.
[18] N. A. Fleck, G. M. Muller, M. F. Ashby, J. W. Hutchinson, Strain gradient plasticity: Theory and
experiment, Acta Metallurgica et Materialia 42 (2) (1994) 475–487. doi:10.1016/0956-7151(94)
90502-9.
[19] J. S. Stlken, A. G. Evans, A microbend test method for measuring the plasticity length scale, Acta
Materialia 46 (14) (1998) 5109–5115. doi:10.1016/S1359-6454(98)00153-0.
[20] D. Lam, F. Yang, A. Chong, J. Wang, P. Tong, Experiments and theory in strain gradient elasticity,
Journal of the Mechanics and Physics of Solids 51 (8) (2003) 1477 – 1508. doi:http://dx.doi.org/

20

10.1016/S0022-5096(03)00053-X.
[21] A. C. Eringen, Nonlocal polar elastic continua, International Journal of Engineering Science 10 (1)
(1972) 1–16. doi:10.1016/0020-7225(72)90070-5.
[22] N. A. Fleck, J. W. Hutchinson, A phenomenological theory for strain gradient effects in plasticity,
Journal of the Mechanics and Physics of Solids 41 (12) (1993) 1825–1857. doi:10.1016/0022-5096(93)
90072-N.
[23] F. Yang, A. C. M. Chong, D. C. C. Lam, P. Tong, Couple stress based strain gradient theory for
elasticity, International Journal of Solids and Structures 39 (10) (2002) 2731–2743. doi:10.1016/
S0020-7683(02)00152-X.
[24] R. D. Mindlin, H. F. Tiersten, Effects of couple-stresses in linear elasticity, Archive for Rational Mechanics and Analysis 11 (1) (1962) 415–448. doi:10.1007/BF00253946.
[25] R. D. Mindlin, Second gradient of strain and surface-tension in linear elasticity, International Journal
of Solids and Structures 1 (4) (1965) 417–438. doi:10.1016/0020-7683(65)90006-5.
[26] R. Ansari, R. Gholami, M. Faghih Shojaei, V. Mohammadi, S. Sahmani, Size-dependent bending,
buckling and free vibration of functionally graded timoshenko microbeams based on the most general
strain gradient theory, Composite Structures 100 (2013) 385–397. doi:10.1016/j.compstruct.2012.
12.048.
[27] S. K. Park, X.-L. Gao, Variational formulation of a modified couple stress theory and its application
to a simple shear problem, Zeitschrift fr angewandte Mathematik und Physik 59 (5) (2007) 904–917.
doi:10.1007/s00033-006-6073-8.
[28] S. K. Park, X.-L. Gao, BernoulliEuler beam model based on a modified couple stress theory, Journal
of Micromechanics and Microengineering 16 (11) (2006) 2355. doi:10.1088/0960-1317/16/11/015.
[29] W. Xia, L. Wang, L. Yin, Nonlinear non-classical microscale beams: Static bending, postbuckling and
free vibration, International Journal of Engineering Science 48 (12) (2010) 2044–2053. doi:10.1016/
j.ijengsci.2010.04.010.
[30] L.-L. Ke, Y.-S. Wang, Size effect on dynamic stability of functionally graded microbeams based
on a modified couple stress theory, Composite Structures 93 (2) (2011) 342–350. doi:10.1016/j.
compstruct.2010.09.008.
[31] C. M. C. Roque, D. S. Fidalgo, A. J. M. Ferreira, J. N. Reddy, A study of a microstructure-dependent
composite laminated Timoshenko beam using a modified couple stress theory and a meshless method,
Composite Structures 96 (2013) 532–537. doi:10.1016/j.compstruct.2012.09.011.
[32] H.-T. Thai, T. P. Vo, T.-K. Nguyen, J. Lee, Size-dependent behavior of functionally graded sandwich
microbeams based on the modified couple stress theory, Composite Structures 123 (2015) 337–349.
doi:10.1016/j.compstruct.2014.11.065.
[33] M. Salamat-talab, A. Nateghi, J. Torabi, Static and dynamic analysis of third-order shear deformation
FG micro beam based on modified couple stress theory, International Journal of Mechanical Sciences
57 (1) (2012) 63–73. doi:10.1016/j.ijmecsci.2012.02.004.
[34] G. C. Tsiatas, A new Kirchhoff plate model based on a modified couple stress theory, International
Journal of Solids and Structures 46 (13) (2009) 2757–2764. doi:10.1016/j.ijsolstr.2009.03.004.
[35] L. Yin, Q. Qian, L. Wang, W. Xia, Vibration analysis of microscale plates based on modified couple
stress theory, Acta Mechanica Solida Sinica 23 (5) (2010) 386 – 393. doi:http://dx.doi.org/10.
1016/S0894-9166(10)60040-7.
[36] H. M. Ma, X.-L. Gao, J. N. Reddy, A non-classical Mindlin plate model based on a modified couple
stress theory, Acta Mechanica 220 (1-4) (2011) 217–235. doi:10.1007/s00707-011-0480-4.
[37] M. Simsek, T. Kocaturk, S. D. Akbas, Static bending of a functionally graded microscale Timoshenko
beam based on the modified couple stress theory, Composite Structures 95 (2013) 740–747. doi:
10.1016/j.compstruct.2012.08.036.
[38] A. Nateghi, M. Salamat-talab, J. Rezapour, B. Daneshian, Size dependent buckling analysis of functionally graded micro beams based on modified couple stress theory, Applied Mathematical Modelling
36 (10) (2012) 4971–4987. doi:10.1016/j.apm.2011.12.035.
[39] J. N. Reddy, J. Kim, A nonlinear modified couple stress-based third-order theory of functionally graded

21

plates, Composite Structures 94 (3) (2012) 1128–1143. doi:10.1016/j.compstruct.2011.10.006.
[40] J. N. Reddy, J. Berry, Nonlinear theories of axisymmetric bending of functionally graded circular
plates with modified couple stress, Composite Structures 94 (12) (2012) 3664–3668. doi:10.1016/j.
compstruct.2012.04.019.
[41] J. Kim, J. N. Reddy, A general third-order theory of functionally graded plates with modified couple
stress effect and the von karman nonlinearity: theory and finite element analysis, Acta Mechanica
226 (9) (2015) 2973–2998. doi:10.1007/s00707-015-1370-y.
[42] J. N. Reddy, J. Romanoff, J. A. Loya, Nonlinear Finite Element Analysis of Functionally Graded
Circular Plates with Modified Couple Stress Theory, European Journal of Mechanics - A/Solidsdoi:
10.1016/j.euromechsol.2015.11.001.
[43] H.-T. Thai, S.-E. Kim, A size-dependent functionally graded reddy plate model based on a modified
couple stress theory, Composites Part B: Engineering 45 (1) (2013) 1636 – 1645. doi:http://dx.doi.
org/10.1016/j.compositesb.2012.09.065.
[44] H.-T. Thai, D.-H. Choi, Size-dependent functionally graded kirchhoff and mindlin plate models based
on a modified couple stress theory, Composite Structures 95 (2013) 142 – 153. doi:http://dx.doi.
org/10.1016/j.compstruct.2012.08.023.
[45] H.-T. Thai, T. P. Vo, A size-dependent functionally graded sinusoidal plate model based on a modified
couple stress theory, Composite Structures 96 (2013) 376–383. doi:10.1016/j.compstruct.2012.09.
025.
[46] L.-L. Ke, Y.-S. Wang, J. Yang, S. Kitipornchai, Free vibration of size-dependent Mindlin microplates
based on the modified couple stress theory, Journal of Sound and Vibration 331 (1) (2012) 94–106.
doi:10.1016/j.jsv.2011.08.020.
[47] L.-L. Ke, J. Yang, S. Kitipornchai, M. A. Bradford, Bending, buckling and vibration of size-dependent
functionally graded annular microplates, Composite Structures 94 (11) (2012) 3250–3257. doi:10.
1016/j.compstruct.2012.04.037.
[48] L. He, J. Lou, E. Zhang, Y. Wang, Y. Bai, A size-dependent four variable refined plate model for
functionally graded microplates based on modified couple stress theory, Composite Structures 130
(2015) 107 – 115. doi:http://dx.doi.org/10.1016/j.compstruct.2015.04.033.
[49] L. Della Croce, P. Venini, Finite elements for functionally graded ReissnerMindlin plates, Computer
Methods in Applied Mechanics and Engineering 193 (911) (2004) 705–725. doi:10.1016/j.cma.2003.
09.014.
[50] H. Nguyen-Xuan, G. R. Liu, C. Thai-Hoang, T. Nguyen-Thoi, An edge-based smoothed finite element
method (ES-FEM) with stabilized discrete shear gap technique for analysis of ReissnerMindlin plates,
Computer Methods in Applied Mechanics and Engineering 199 (912) (2010) 471–489. doi:10.1016/
j.cma.2009.09.001.
[51] B. Zhang, Y. He, D. Liu, Z. Gan, L. Shen, A non-classical Mindlin plate finite element based on
a modified couple stress theory, European Journal of Mechanics - A/Solids 42 (2013) 63–80. doi:
10.1016/j.euromechsol.2013.04.005.
[52] J. N. Reddy, A Simple Higher-Order Theory for Laminated Composite Plates, Journal of Applied
Mechanics 51 (4) (1984) 745–752. doi:10.1115/1.3167719.
[53] H. Nguyen-Xuan, C. H. Thai, T. Nguyen-Thoi, Isogeometric finite element analysis of composite sandwich plates using a higher order shear deformation theory, Composites Part B: Engineering 55 (2013)
558–574. doi:10.1016/j.compositesb.2013.06.044.
[54] H. Arya, R. P. Shimpi, N. K. Naik, A zigzag model for laminated composite beams, Composite Structures 56 (1) (2002) 21–24. doi:10.1016/S0263-8223(01)00178-7.
[55] N. R. Senthilnathan, S. P. Lim, K. H. Lee, S. T. Chow, Buckling of Shear-Deformable Plates, AIAA
Journal 25 (9) (1987) 1268–1271. doi:10.2514/3.48742.
[56] R. P. Shimpi, Refined plate theory and its variants, AIAA journal 40 (1) (2002) 137–146.
[57] R. P. Shimpi, H. G. Patel, A two variable refined plate theory for orthotropic plate analysis, International Journal of Solids and Structures 43 (2223) (2006) 6783–6799. doi:10.1016/j.ijsolstr.2006.
02.007.

22

[58] R. P. Shimpi, H. G. Patel, Free vibrations of plate using two variable refined plate theory, Journal of
Sound and Vibration 296 (45) (2006) 979–999. doi:10.1016/j.jsv.2006.03.030.
[59] C. A. Shankara, N. G. R. Iyengar, A C0 element for the free vibration analysis of laminated composite
plates, Journal of Sound and Vibration 191 (5) (1996) 721–738. doi:10.1006/jsvi.1996.0152.
[60] T. Hughes, J. Cottrell, Y. Bazilevs, Isogeometric analysis: Cad, finite elements, nurbs, exact geometry
and mesh refinement, Computer Methods in Applied Mechanics and Engineering 194 (3941) (2005)
4135 – 4195. doi:http://dx.doi.org/10.1016/j.cma.2004.10.008.
[61] H. Nguyen-Xuan, L. V. Tran, C. H. Thai, S. Kulasegaram, S. Bordas, Isogeometric analysis of functionally graded plates using a refined plate theory, Composites Part B: Engineering 64 (2014) 222 –
234. doi:http://dx.doi.org/10.1016/j.compositesb.2014.04.001.
[62] A. V. Vuong, C. Heinrich, B. Simeon, ISOGAT: A 2d tutorial MATLAB code for Isogeometric Analysis,
Computer Aided Geometric Design 27 (8) (2010) 644–655. doi:10.1016/j.cagd.2010.06.006.
[63] C. de Falco, A. Reali, R. Vzquez, GeoPDEs: A research tool for Isogeometric Analysis of PDEs,
Advances in Engineering Software 42 (12) (2011) 1020–1034. doi:10.1016/j.advengsoft.2011.06.
010.
[64] V. P. Nguyen, C. Anitescu, S. P. A. Bordas, T. Rabczuk, Isogeometric analysis: An overview and
computer implementation aspects, Mathematics and Computers in Simulation 117 (2015) 89–116. doi:
10.1016/j.matcom.2015.05.008.
[65] D. Wang, W. Liu, H. Zhang, Superconvergent isogeometric free vibration analysis of EulerBernoulli
beams and Kirchhoff plates with new higher order mass matrices, Computer Methods in Applied
Mechanics and Engineering 286 (2015) 230–267. doi:10.1016/j.cma.2014.12.026.
[66] L. Beiro da Veiga, A. Buffa, C. Lovadina, M. Martinelli, G. Sangalli, An isogeometric method for the
ReissnerMindlin plate bending problem, Computer Methods in Applied Mechanics and Engineering
209212 (2012) 45–53. doi:10.1016/j.cma.2011.10.009.
[67] C. H. Thai, H. Nguyen-Xuan, N. Nguyen-Thanh, T.-H. Le, T. Nguyen-Thoi, T. Rabczuk, Static, free
vibration, and buckling analysis of laminated composite ReissnerMindlin plates using NURBS-based
isogeometric approach, International Journal for Numerical Methods in Engineering 91 (6) (2012) 571–
603. doi:10.1002/nme.4282.
[68] L. V. Tran, A. Ferreira, H. Nguyen-Xuan, Isogeometric analysis of functionally graded plates using
higher-order shear deformation theory, Composites Part B: Engineering 51 (2013) 368 – 383. doi:
http://dx.doi.org/10.1016/j.compositesb.2013.02.045.
[69] C. H. Thai, A. J. M. Ferreira, E. Carrera, H. Nguyen-Xuan, Isogeometric analysis of laminated composite and sandwich plates using a layerwise deformation theory, Composite Structures 104 (2013)
196–214. doi:10.1016/j.compstruct.2013.04.002.
[70] L. V. Tran, C. H. Thai, H. Nguyen-Xuan, An isogeometric finite element formulation for thermal
buckling analysis of functionally graded plates, Finite Elements in Analysis and Design 73 (2013) 65–
76. doi:10.1016/j.finel.2013.05.003.
[71] P. Phung-Van, M. Abdel-Wahab, K. M. Liew, S. P. A. Bordas, H. Nguyen-Xuan, Isogeometric analysis of
functionally graded carbon nanotube-reinforced composite plates using higher-order shear deformation
theory, Composite Structures 123 (2015) 137–149. doi:10.1016/j.compstruct.2014.12.021.
[72] L. V. Tran, J. Lee, H. Nguyen-Van, H. Nguyen-Xuan, M. A. Wahab, Geometrically nonlinear isogeometric analysis of laminated composite plates based on higher-order shear deformation theory, International
Journal of Non-Linear Mechanics 72 (2015) 42–52. doi:10.1016/j.ijnonlinmec.2015.02.007.
[73] C. H. Thai, H. Nguyen-Xuan, S. P. A. Bordas, N. Nguyen-Thanh, T. Rabczuk, Isogeometric Analysis
of Laminated Composite Plates Using the Higher-Order Shear Deformation Theory, Mechanics of
Advanced Materials and Structures 22 (6) (2015) 451–469. doi:10.1080/15376494.2013.779050.
[74] S. S. Vel, R. C. Batra, Exact Solution for Thermoelastic Deformations of Functionally Graded Thick
Rectangular Plates, AIAA Journal 40 (7) (2002) 1421–1433. doi:10.2514/2.1805.
[75] L. F. Qian, R. C. Batra, L. M. Chen, Static and dynamic deformations of thick functionally
graded elastic plates by using higher-order shear and normal deformable plate theory and meshless local PetrovGalerkin method, Composites Part B: Engineering 35 (68) (2004) 685–697. doi:

23

10.1016/j.compositesb.2004.02.004.
[76] T. Mori, K. Tanaka, Average stress in matrix and average elastic energy of materials with misfitting
inclusions, Acta Metallurgica 21 (5) (1973) 571–574. doi:10.1016/0001-6160(73)90064-3.
[77] Y. Kiani, M. R. Eslami, Thermal Postbuckling of Imperfect Circular Functionally Graded Material
Plates: Examination of Voigt, MoriTanaka, and Self-Consistent Schemes, Journal of Pressure Vessel
Technology 137 (2) (2014) 021201–021201. doi:10.1115/1.4026993.
[78] E. Carrera, S. Brischetto, M. Cinefra, M. Soave, Effects of thickness stretching in functionally graded
plates and shells, Composites Part B: Engineering 42 (2) (2011) 123 – 133. doi:http://dx.doi.org/
10.1016/j.compositesb.2010.10.005.
[79] J. L. Mantari, C. Guedes Soares, A novel higher-order shear deformation theory with stretching effect
for functionally graded plates, Composites Part B: Engineering 45 (1) (2013) 268–281. doi:10.1016/
j.compositesb.2012.05.036.
[80] C. H. Thai, A. M. Zenkour, M. Abdel Wahab, H. Nguyen-Xuan, A simple four-unknown shear and
normal deformations theory for functionally graded isotropic and sandwich plates based on isogeometric
analysis, Composite Structures 139 (2016) 77–95. doi:10.1016/j.compstruct.2015.11.066.
[81] A. M. Zenkour, A simple four-unknown refined theory for bending analysis of functionally graded plates,
Applied Mathematical Modelling 37 (2013) 9041 – 9051. doi:http://dx.doi.org/10.1016/j.apm.
2013.04.022.
[82] A. M. Zenkour, Bending analysis of functionally graded sandwich plates using a simple four-unknown
shear and normal deformations theory, Journal of Sandwich Structures and Materials 15 (6) (2013)
629–656. doi:10.1177/1099636213498886.
[83] T. N. Nguyen, C. H. Thai, H. Nguyen-Xuan, On the general framework of high order shear deformation
theories for laminated composite plate structures: A novel unified approach, International Journal of
Mechanical Sciencesdoi:10.1016/j.ijmecsci.2016.01.012.
[84] J. N. Reddy, Theory and Analysis of Elastic Plates and Shells, 2nd Edition, CRC Press, 2006.
[85] N.-T. Nguyen, D. Hui, J. Lee, H. Nguyen-Xuan, An efficient computational approach for size-dependent
analysis of functionally graded nanoplates, Computer Methods in Applied Mechanics and Engineering
297 (2015) 191 – 218. doi:http://dx.doi.org/10.1016/j.cma.2015.07.021.
[86] J. A. Cottrell, T. J. R. Hughes, Y. Bazilevs, Isogeometric Analysis: Toward Integration of CAD and
FEA, 1st Edition, Wiley Publishing, 2009.
[87] F. Auricchio, L. B. da Veiga, A. Buffa, C. Lovadina, A. Reali, G. Sangalli, A fully locking-free isogeometric approach for plane linear elasticity problems: A stream function formulation, Computer
Methods in Applied Mechanics and Engineering 197 (14) (2007) 160 – 172. doi:http://dx.doi.org/
10.1016/j.cma.2007.07.005.
[88] M. Mohammadi, M. Ghayour, A. Farajpour, Free transverse vibration analysis of circular and annular
graphene sheets with various boundary conditions using the nonlocal continuum plate model, Composites Part B: Engineering 45 (1) (2013) 32 – 42. doi:http://dx.doi.org/10.1016/j.compositesb.
2012.09.011.
[89] L. Ma, T. Wang, Relationships between axisymmetric bending and buckling solutions of {FGM} circular
plates based on third-order plate theory and classical plate theory, International Journal of Solids and
Structures 41 (1) (2004) 85 – 101. doi:http://dx.doi.org/10.1016/j.ijsolstr.2003.09.008.
[90] A. Saidi, A. Rasouli, S. Sahraee, Axisymmetric bending and buckling analysis of thick functionally
graded circular plates using unconstrained third-order shear deformation plate theory, Composite Structures 89 (1) (2009) 110 – 119. doi:http://dx.doi.org/10.1016/j.compstruct.2008.07.003.
[91] H.-T. Thai, S.-E. Kim, A simple quasi-3d sinusoidal shear deformation theory for functionally graded
plates, Composite Structures 99 (2013) 172–180. doi:10.1016/j.compstruct.2012.11.030.
[92] A. Neves, A. Ferreira, E. Carrera, C. Roque, M. Cinefra, R. Jorge, C. Soares, A quasi-3d sinusoidal shear
deformation theory for the static and free vibration analysis of functionally graded plates, Composites
Part B: Engineering 43 (2) (2012) 711 – 725. doi:http://dx.doi.org/10.1016/j.compositesb.
2011.08.009.
[93] J. L. Mantari, C. G. Soares, A quasi-3d tangential shear deformation theory with four un-

24

[94]

[95]

[96]

[97]

[98]

knowns for functionally graded plates, Acta Mechanica 226 (3) (2014) 625–642. doi:10.1007/
s00707-014-1192-3.
S. Akavci, A. Tanrikulu, Static and free vibration analysis of functionally graded plates based on a
new quasi-3d and 2d shear deformation theories, Composites Part B: Engineering 83 (2015) 203 – 215.
doi:http://dx.doi.org/10.1016/j.compositesb.2015.08.043.
S. S. Vel, R. Batra, Three-dimensional exact solution for the vibration of functionally graded rectangular
plates, Journal of Sound and Vibration 272 (35) (2004) 703 – 730. doi:http://dx.doi.org/10.1016/
S0022-460X(03)00412-7.
H. Matsunaga, Free vibration and stability of functionally graded plates according to a 2-d higher-order
deformation theory, Composite Structures 82 (4) (2008) 499 – 512. doi:http://dx.doi.org/10.1016/
j.compstruct.2007.01.030.
Z. Belabed, M. S. A. Houari, A. Tounsi, S. Mahmoud, O. A. Bg, An efficient and simple higher order
shear and normal deformation theory for functionally graded material (fgm) plates, Composites Part B:
Engineering 60 (2014) 274 – 283. doi:http://dx.doi.org/10.1016/j.compositesb.2013.12.057.
F. Alijani, M. Amabili, Effect of thickness deformation on large-amplitude vibrations of functionally
graded rectangular plates, Composite Structures 113 (2014) 89 – 107. doi:http://dx.doi.org/10.
1016/j.compstruct.2014.03.006.

25

Table 1: Various forms of distribution function used for HSDT, RPT, and quasi-3D theories
Theory
εz
HSDT [52]
=0
RPT [61]
=0
Quasi-3D [91]
6= 0
Quasi-3D [81]
6= 0
Quasi-3D [85]
6= 0
Present RPT
=0
Present quasi-3D 6= 0

f (z)
3
z − 43 hz 2

π
arctan sin
z
h

h
−z
sin πz
π
h 

z
4z 3
1
h sinh h − 3h
2 cosh 2
9π 3
28π 5
π
z − 5h
3 z + 25h5 z
h
3
6z 5
8z 7
−8z + 10z
+ 5h
4 + 7h6
h2
3
6z 5
8z 7
−8z + 10z
+ 5h
4 + 7h6
h2

26

φz
f 0 (z) + 1
1 0
f (z)
12
1 0
f
(z)
8
3 0
f (z)
20

Table 2: Material properties of FG plates
Property

Material
Al
Al2 O3
E (GPa)
70
380
ν
0.3
0.3
3
ρ (kg/m ) 2707 3800

27

ZrO2 -1 ZrO2 -2
200
151
0.3
0.3
3000
3000

10Eh3
w (a/2, a/2, 0) of
q0 L 4
SSSS homogeneous square plate subjected to sinusoidally distributed load, a/h = 20
Table 3: Convergence of non-dimensional central deflection w̄ =

l/h

p Element Mesh
3×3
5×5

Analytical [43]
7×7

9×9

11×11

13×13

15×15

17×17

0

0.2842
2
3
4

0.2734
0.2823
0.2843

0.2799
0.2841
0.2842

0.2819
0.2842
0.2842

0.2828
0.2842
0.2842

0.2833
0.2842
0.2842

0.2835
0.2842
0.2842

0.2837
0.2842
0.2842

0.2838
0.2842
0.2842

2
3
4

0.2346
0.2415
0.2432

0.2397
0.2430
0.2431

0.2413
0.2431
0.2431

0.2420
0.2431
0.2431

0.2424
0.2431
0.2431

0.2426
0.2431
0.2431

0.2427
0.2431
0.2431

0.2428
0.2431
0.2431

0.2

0.2430

0.6

0.1124
2
3
4

0.1098
0.1120
0.1127

0.1115
0.1127
0.1127

0.1121
0.1127
0.1127

0.1123
0.1127
0.1127

0.1125
0.1127
0.1127

0.1125
0.1127
0.1127

0.1126
0.1127
0.1127

0.1126
0.1127
0.1127

2
3
4

0.0532
0.0541
0.0544

0.0539
0.0544
0.0544

0.0541
0.0544
0.0544

0.0542
0.0544
0.0544

0.0543
0.0544
0.0544

0.0543
0.0544
0.1127

0.0543
0.0544
0.1127

0.0543
0.0544
0.1127

1

0.0542

28

Table 4: Comparison of non-dimensional deflection w̄ =
Al/Al2 O3 square plates (rule of mixtures scheme)
n

Theory

10Ec h3
w (a/2, b/2, 0) of SSSS
q0 L4

εz = 0
a/h = 4 a/h = 10 a/h = 100

Sinusoidally distributed load
1 Zenkour [81]
Neves et al. [92]
Mantari and Soares [93]
Carrera et al. [78]
0.7289
Akavci and Tanrikulu [94] 0.7282
Present†
0.7284
4 Zenkour [81]
Neves et al. [92]
Mantari and Soares [93]
Carrera et al. [78]
1.1673
Akavci and Tanrikulu [94] 1.1613
Present†
1.1590
10 Zenkour [81]
Neves et al. [92]
Mantari and Soares [93]
Carrera et al. [78]
1.3925
Akavci and Tanrikulu [94] 1.3917
Present†
1.3902
Uniformly distributed load
1 Zenkour [7]
Akavci and Tanrikulu [94] Present†
1.1319
4 Zenkour [7]
Akavci and Tanrikulu [94] Present†
1.7928
10 Zenkour [7]
Akavci and Tanrikulu [94] Present†
2.1432
†
Proposed RPT and quasi-3D models are

0.5890
0.5889
0.5889
0.8828
0.8818
0.8813
1.0090
1.0089
1.0086

0.5625
0.5625
0.5625
0.8286
0.8287
0.8287
0.9361
0.9362
0.9362

0.9287
0.9288
0.9288
0.8904
1.3890
1.3888
1.3882
1.3116
1.5876
1.5875
1.5870
1.4818
used for the case of εz = 0

29

εz 6= 0
a/h = 4 a/h = 10 a/h = 100
0.6828
0.6997
0.693
0.7171
0.6908
0.6935
1.1001
1.1178
1.085
1.1585
1.0983
1.0868
1.3391
1.3490
1.308
1.3745
1.3352
1.3116

0.5592
0.5845
0.569
0.5875
0.5691
0.5691
0.8404
0.8750
0.838
0.8821
0.8417
0.8392
0.9806
0.8750
0.972
1.0072
0.9818
0.9748

0.5624
0.5624
0.545
0.5625
0.5457
0.5460
0.7933
0.8286
0.793
0.8286
0.7925
0.7933
0.9140
0.8286
0.911
0.9361
0.9141
0.9132

0.8977
1.0788
0.8978
0.8642
1.3259
1.6827
1.3223
1.2556
1.5453
2.0244
1.5347
1.4454
and εz 6= 0, respectively.

10Ec h3
w (a/2, b/2, 0) of SSSS Al/Al2 O3 square microplates
q0 L4
subjected to sinusoidally distributed load (rule of mixtures scheme)

Table 5: Non-dimensional deflection w̄ =

a/h
5

20

100

l/h
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1

n=0
RPT
0.3433
0.2898
0.1975
0.1292
0.0871
0.0614
0.2842
0.2431
0.1695
0.1127
0.0767
0.0544
0.2804
0.2401
0.1677
0.1116
0.0760
0.0539

Quasi-3D
0.3360
0.2853
0.1965
0.1296
0.0879
0.0623
0.2836
0.2427
0.1694
0.1127
0.0768
0.0544
0.2803
0.2400
0.1677
0.1116
0.0760
0.0539

TSDT[43]
0.3433
0.2875
0.1934
0.1251
0.0838
0.0588
0.2842
0.2430
0.1693
0.1124
0.0765
0.0542
0.2804
0.2401
0.1677
0.1116
0.0760
0.0539

n=1
RPT
0.6688
0.5505
0.3601
0.2288
0.1517
0.1060
0.5689
0.4739
0.3157
0.2029
0.1352
0.0947
0.5625
0.4689
0.3128
0.2012
0.1341
0.0939

30

Quasi-3D
0.6401
0.5321
0.3537
0.2274
0.1520
0.1069
0.5516
0.4619
0.3105
0.2008
0.1343
0.0943
0.5460
0.4574
0.3076
0.1990
0.1332
0.0934

TSDT[43]
0.6688
0.5468
0.3535
0.2224
0.1464
0.1017
0.5689
0.4737
0.3153
0.2025
0.1349
0.0944
0.5625
0.4689
0.3128
0.2011
0.1341
0.0939

n = 10
RPT
1.2271
1.0400
0.7140
0.4694
0.3174
0.2242
0.9537
0.8313
0.6001
0.4102
0.2842
0.2038
0.9362
0.8176
0.5925
0.4062
0.2820
0.2024

Quasi-3D
1.1663
1.0019
0.7043
0.4711
0.3220
0.2289
0.9280
0.8120
0.5906
0.4061
0.2825
0.2031
0.9132
0.8001
0.5833
0.4018
0.2799
0.2014

TSDT[43]
1.2276
1.0247
0.6908
0.4514
0.3052
0.2158
0.9538
0.8303
0.5986
0.4090
0.2834
0.2033
0.9362
0.8176
0.5925
0.4061
0.2820
0.2024

Table 6: Non-dimensional deflection w̄ =
(rule of mixtures scheme)
a/h

n=0
RPT
Quasi-3D
Sinusoidally distributed load
5
0
0.1601 0.1359
0.2 0.1378 0.1197
0.4 0.0974 0.0883
0.6 0.0655 0.0616
0.8 0.0450 0.0435
1
0.0321 0.0316
20
0
0.1035 0.0950
0.2 0.0919 0.0849
0.4 0.0688 0.0645
0.6 0.0485 0.0462
0.8 0.0343 0.0331
1
0.0250 0.0243
100 0
0.0999 0.0955
0.2 0.0889 0.0853
0.4 0.0668 0.0646
0.6 0.0473 0.0461
0.8 0.0336 0.0329
1
0.0244 0.0241
Uniformly
5
0
0.2
0.4
0.6
0.8
1
20
0
0.2
0.4
0.6
0.8
1
100 0
0.2
0.4
0.6
0.8
1

10Ec h3
w (a/2, b/2, 0) of CCCC Al/Al2 O3 square microplates
q0 L4

l/h

distributed load
0.2239 0.1860
0.1924 0.1641
0.1358 0.1215
0.0914 0.0851
0.0627 0.0602
0.0447 0.0439
0.1436 0.1300
0.1275 0.1164
0.0956 0.0887
0.0674 0.0637
0.0478 0.0457
0.0348 0.0336
0.1384 0.1314
0.1232 0.1174
0.0927 0.0891
0.0657 0.0637
0.0467 0.0456
0.0340 0.0334

IGA-Reddy

n=1
RPT

IGA-Reddy

n = 10
RPT
Quasi-3D

Quasi-3D

IGA-Red

0.1601
0.1377
0.0973
0.0655
0.0449
0.0320
0.1035
0.0919
0.0688
0.0485
0.0343
0.0250
0.0999
0.0889
0.0668
0.0473
0.0336
0.0244

0.3021
0.2555
0.1751
0.1151
0.0779
0.0551
0.2065
0.1797
0.1294
0.0882
0.0611
0.0438
0.2003
0.1746
0.1262
0.0863
0.0599
0.0430

0.2554
0.2214
0.1586
0.1081
0.0752
0.0542
0.1863
0.1638
0.1204
0.0837
0.0587
0.0425
0.1872
0.1643
0.1204
0.0835
0.0584
0.0422

0.3020
0.2554
0.1751
0.1151
0.0779
0.0551
0.2065
0.1797
0.1294
0.0882
0.0611
0.0438
0.2003
0.1746
0.1262
0.0863
0.0599
0.0430

0.6111
0.5178
0.3568
0.2358
0.1602
0.1136
0.3505
0.3150
0.2419
0.1746
0.1258
0.0926
0.3336
0.3013
0.2336
0.1701
0.1232
0.0910

0.4740
0.4199
0.3132
0.2204
0.1559
0.1134
0.3150
0.2857
0.2237
0.1647
0.1204
0.0896
0.3131
0.2842
0.2228
0.1640
0.1199
0.0891

0.6113
0.5183
0.3575
0.2364
0.1606
0.1139
0.3505
0.3151
0.2420
0.1747
0.1258
0.0926
0.3336
0.3013
0.2336
0.1701
0.1233
0.0910

0.2238
0.1924
0.1358
0.0913
0.0627
0.0447
0.1436
0.1275
0.0956
0.0674
0.0478
0.0348
0.1384
0.1232
0.0927
0.0657
0.0467
0.0340

0.4220
0.3566
0.2443
0.1606
0.1087
0.0769
0.2864
0.2493
0.1797
0.1227
0.0850
0.0610
0.2775
0.2421
0.1751
0.1199
0.0833
0.059831

0.3500
0.3040
0.2186
0.1495
0.1041
0.0752
0.2553
0.2248
0.1657
0.1156
0.0813
0.0589
0.2577
0.2265
0.1664
0.1155
0.0810
0.0586

0.4219
0.3565
0.2442
0.1605
0.1087
0.0769
0.2864
0.2493
0.1797
0.1227
0.0850
0.0610
0.2775
0.2421
0.1751
0.1199
0.0833
0.0598

0.8557
0.7233
0.4976
0.3288
0.2234
0.1584
0.4863
0.4372
0.3360
0.2428
0.1750
0.1289
0.4622
0.4176
0.3242
0.2363
0.1713
0.1266

0.6408
0.5701
0.4286
0.3034
0.2154
0.1570
0.4307
0.3912
0.3072
0.2269
0.1664
0.1240
0.4310
0.3915
0.3075
0.2268
0.1661
0.1236

0.8559
0.7239
0.4985
0.3295
0.2239
0.1588
0.4863
0.4373
0.3361
0.2429
0.1751
0.1289
0.4622
0.4176
0.3242
0.2363
0.1713
0.1266]

a2
Table 7: Comparison of non-dimensional natural frequencies ω̄ = ω
h
Al/ZrO2 -1 plates (Mori-Tanaka scheme)
Theory
Vel and Batra [95]
Matsunaga [96]
Neves et al. [92]
Belabed et al. [97]
Alijani and Amabili [98]
Akavci and Tanrikulu [94]
Present

n=1
a/h = 5
5.4806
5.7123
5.4825
5.4800
5.4796
5.4829
5.5172

a/h = 10
5.9609
6.1932
5.9600
5.9700
5.9578
5.9676
6.0023

32

a/h = 20
6.1076
6.3390
6.1200
6.1200
6.1040
6.1160
6.1505

a/h = 5
n=2
5.4923
5.6599
5.4950
5.5025
5.4919
5.5064
5.5324

r

ρm
of SSSS
Em

n=3
5.5285
5.6757
5.5300
5.5350
5.5279
5.5388
5.5642

n=5
5.5632
5.7020
5.5625
5.5625
5.5633
5.5644
5.5886

a2
Table 8: Comparison of non-dimensional natural frequencies ω̄ = ω
h
geneous microplates

ρ
of SSSS homoE

a/h Theory
5

20

100

l/h
0
CPT [35]
5.9734
TSDT [43]
5.2813
RPT (Present)
5.2813
Quasi-3D (Present) 5.3090
CPT [35]
5.9734
TSDT [43]
5.9199
RPT (Present)
5.9199
Quasi-3D (Present) 5.9235
CPT [35]
5.9734
TSDT [43]
5.9712
RPT (Present)
5.9712
Quasi-3D (Present) 5.9723

r

0.2
6.4556
5.7699
5.7496
5.7622
6.4556
6.4027
6.4009
6.4030
6.4556
6.4535
6.4534
6.4544

33

0.4
7.7239
7.0330
6.9667
6.9438
7.7239
7.6708
7.6646
7.6633
7.7239
7.7217
7.7215
7.7222

0.6
9.4673
8.7389
8.6191
8.5509
9.4673
9.4116
9.4005
9.3952
9.4673
9.4651
9.4646
9.4650

0.8
11.4713
10.6766
9.8943
9.8943
11.4713
11.4108
11.3945
11.3854
11.4713
11.4689
11.4682
11.4683

1
13.6213
12.7408
9.9791
9.9791
13.6213
13.5545
13.5330
13.5202
13.6213
13.6186
13.6178
13.6177

a2
Table 9: Non-dimensional natural frequency ω̄ = ω
h
mixtures scheme)
a/h
5

20

100

l/h
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1

n=0
RPT
5.2813
5.7496
6.9667
8.6191
9.8943
9.9791
5.9199
6.4009
7.6646
9.4005
11.3945
13.5330
5.9712
6.4534
7.7215
9.4646
11.4682
13.6178

Quasi-3D TSDT[43]
5.3090
5.2813
5.7622
5.7699
6.9438
7.0330
8.5509
8.7389
9.8943
10.6766
9.9791
12.7408
5.9235
5.9199
6.4030
6.4027
7.6633
7.6708
9.3952
9.4116
11.3854
11.4108
13.5202
13.5545
5.9723
5.9712
6.4544
6.4535
7.7222
7.7217
9.4650
9.4651
11.4683
11.4689
13.6177
13.6186

r

ρc
of SSSS Al/Al2 O3 square plates (rule of
Ec

n=1
RPT
Quasi-3D
4.0781 4.1521
4.4959 4.5542
5.5620 5.5865
6.9822 6.9681
8.2313 8.2313
8.3019 8.3019
4.5228 4.5919
4.9556 5.0179
6.0714 6.1203
7.5739 7.6107
9.2768 9.3042
11.0882 11.1082
4.5579 4.6263
4.9922 5.0546
6.1124 6.1635
7.6220 7.6630
9.3339 9.3673
11.1554 11.1832

34

TSDT[43]
4.0781
4.5094
5.6071
7.0662
8.7058
10.4397
4.5228
4.9568
6.0756
7.5817
9.2887
11.1042
4.5579
4.9922
6.1126
7.6224
9.3344
11.1560

n = 10
RPT
3.2519
3.5312
4.2584
5.2471
5.8571
5.9073
3.7623
4.0299
4.7428
5.7369
6.8914
8.1384
3.8058
4.0724
4.7837
5.7778
6.9341
8.1842

Quasi-3D
3.3126
3.5740
4.2627
5.2115
5.8571
5.9073
3.8129
4.0761
4.7794
5.7640
6.9106
8.1510
3.8533
4.1168
4.8215
5.8090
6.9600
8.2060

TSDT[4
3.2514
3.5548
4.3200
5.3335
6.4759
7.6895
3.7622
4.0323
4.7488
5.7453
6.9013
8.1494
3.8058
4.0725
4.7840
5.7782
6.9345
8.1846

a2
Table 10: The first six non-dimensional natural frequencies ω̄ = ω
h
square plates (Mori-Tanaka scheme)
BC
SSSS

a/h Theory
5

10

100

CCCC

5

10

100

IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)

Mode
1
7.9366
7.8883
7.7844
8.5940
8.5607
8.4401
8.8399
8.8390
8.7127
12.7213
13.1029
12.1531
15.2379
15.4413
14.3638
16.0795
16.0540
15.5212

35

2
13.8049
13.8049
13.8049
20.3230
20.2031
19.9188
21.7851
21.7802
21.4598
22.6661
22.9300
22.0683
28.9665
29.3267
27.8276
32.6473
32.6082
31.5708

3
13.8049
13.8049
13.8049
20.3230
20.2031
19.9188
21.7851
21.7802
21.4598
22.6661
22.9300
22.0683
28.9665
29.3267
27.8276
32.6473
32.6082
31.5708

r

ρm
of Al/Al2 O3
Em

4
17.3259
17.2045
16.9943
27.5893
27.5893
27.5893
35.3333
35.3221
34.8171
27.9021
27.8791
26.3182
41.2221
41.5955
39.9305
48.5912
48.5477
47.0903

5
17.3259
17.2045
16.9943
27.5893
27.5893
27.5893
43.0764
43.0573
42.4074
27.9021
27.8791
26.3182
47.8649
48.2657
46.0435
58.0858
58.0353
56.2566

6
19.5422
19.5422
19.5422
31.7475
31.5547
31.1386
43.0764
43.0573
42.4074
31.1450
31.3005
30.5398
47.9221
48.4504
46.4772
58.3951
58.3402
56.5207

r
Table 11: The first six non-dimensional natural frequencies ω̄ = ωR

2

circular microplates
l/h Theory
Simple support
0
Mohammadi et al. [88]
Nguyen et al. [85]
RPT (Present)
Quasi-3D (Present)
0.2 RPT (Present)
Quasi-3D (Present)
0.4 RPT (Present)
Quasi-3D (Present)
0.6 RPT (Present)
Quasi-3D (Present)
0.8 RPT (Present)
Quasi-3D (Present)
1
RPT (Present)
Quasi-3D (Present)
Clamped support
0
Mohammadi et al. [88]
Nguyen et al. [85]
RPT (Present)
Quasi-3D (Present)
0.2 RPT (Present)
Quasi-3D (Present)
0.4 RPT (Present)
Quasi-3D (Present)
0.6 RPT (Present)
Quasi-3D (Present)
0.8 RPT (Present)
Quasi-3D (Present)
1
RPT (Present)
Quasi-3D (Present)
Eh3
†
D=
12(1 − ν 2 )

ρh
of homogeneous
D†

Mode
1

2

3

4

5

6

4.9345
4.9304
4.9304
4.9385
4.9925
5.0024
5.1213
5.1365
5.2422
5.2649
5.3324
5.3642
5.3954
5.4379

13.8981
13.8587
13.8591
13.8701
14.5095
14.5206
16.2743
16.2857
18.8078
18.8192
21.8342
21.8450
25.1769
25.1869

25.6132
25.4798
25.4799
25.4983
26.5426
26.5613
29.4369
29.4529
33.6370
33.6497
38.6993
38.7080
44.3292
44.3335

29.7198
29.5390
29.5456
29.5691
31.1786
31.1981
35.5996
35.6092
41.9096
41.9065
49.3881
49.3713
57.5842
57.5538

39.9571
39.6331
39.6518
39.6881
41.8855
41.9148
47.7406
47.7542
55.8684
55.8623
65.3972
65.3694
75.8223
75.7721

48.4788
48.0046
48.0402
48.0906
50.8427
50.8808
58.5547
58.5621
69.6631
69.6331
82.8177
82.7497
97.1772
97.0720

10.2158
10.1839
10.1842
10.4466
10.8087
11.0612
12.4963
12.7255
14.8897
15.0927
17.7051
17.8850
20.7715
20.9330

21.2604
21.1433
21.1459
21.6458
22.4449
22.9236
25.9527
26.3811
30.9242
31.2968
36.7699
37.0938
43.1360
43.4206

34.8772
34.5892
34.5885
35.2774
36.3961
37.1800
41.2953
41.9956
48.3696
48.9732
56.8123
57.3265
66.1009
66.5388

39.7706
39.3624
39.3832
40.2833
41.8103
42.6664
48.3406
49.0934
57.5674
58.2073
68.4011
68.9439
80.1939
80.6591

51.0295
50.4385
50.4865
51.5045
53.5802
54.5550
61.7419
62.6149
73.1139
73.8637
86.4224
87.0506
100.9243
101.4425

60.8290
59.9580
60.0416
61.3186
63.7743
64.9738
73.9729
74.9922
88.5334
89.3634
105.6689
106.3411
124.3039
124.8498

36

Table 12: The first six non-dimensional natural frequencies ω̄ = ωR

2

circular microplates (Mori-Tanaka scheme)
l/h Theory
Simple support
0
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.2 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.4 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.6 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.8 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
1
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
Clamped support
0
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.2 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.4 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.6 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
0.8 IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
1
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
Ec h3
†
Dc =
12(1 − νc2 )

r

ρc h
Dc†

of Al/Al2 O3

Mode
1

2

3

4

5

6

3.4629
3.4132
3.3572
3.5224
3.4675
3.4118
3.6392
3.5698
3.5114
3.7577
3.6687
3.6018
3.8632
3.7546
3.6721
3.9548
3.8303
3.9557

8.9254
8.8258
8.6722
9.3961
9.2928
9.1595
10.6411
10.5288
10.4415
12.2959
12.1873
12.1378
13.2220
13.2030
13.1245
13.4288
13.4243
13.5531

8.9254
8.8258
8.6722
9.3961
9.2928
9.1595
10.6411
10.5288
10.4415
12.2959
12.1873
12.1378
13.2220
13.2030
13.1245
13.4288
13.4243
13.5531

12.9890
12.9916
12.9121
13.1894
13.1928
13.1037
13.3115
13.3143
13.2326
13.5601
13.5445
13.4872
14.8766
14.7552
14.7850
17.1268
16.9751
17.4735

12.9890
12.9916
12.9121
13.1894
13.1928
13.1037
13.3115
13.3143
13.2326
13.5601
13.5445
13.4872
14.8766
14.7552
14.7850
17.1268
16.9751
17.4735

15.4780
15.3331
15.0490
16.2612
16.1022
15.8706
18.2644
18.1203
18.0361
18.4956
18.4939
18.4475
18.6662
18.6655
18.6192
18.8623
18.8618
18.8633

6.7718
6.8745
6.3384
7.2163
7.3195
6.8208
8.4091
8.5121
8.0909
10.0860
10.1880
9.8463
12.0447
12.1441
11.8741
14.1674
14.2615
14.0568

12.9968
13.1770
12.4133
13.9289
14.0950
13.4172
16.4035
16.5390
16.0391
19.8463
19.9431
19.6321
23.8379
23.8850
23.7016
26.0476
25.9953
37
24.4150

12.9968
13.1770
12.4133
13.9289
14.0950
13.4172
16.4035
16.5390
16.0391
19.8463
19.9431
19.6321
23.8379
23.8850
23.7016
26.0476
25.9953
24.4150

19.8267
20.0202
19.0879
21.1653
21.3095
20.5401
24.6620
24.7273
24.2526
25.9424
25.8948
24.3655
25.9920
25.9507
24.4627
28.1435
28.1492
28.2217

19.8608
20.0388
19.1363
21.5889
21.7503
20.9799
25.9011
25.8504
24.3345
25.9424
25.8948
24.3655
25.9920
25.9507
24.4627
28.1435
28.1492
28.2217

21.9673
22.2840
20.9877
23.7182
23.9885
22.8902
25.9011
25.8504
24.3345
28.6525
28.6525
28.6525
28.9380
28.9380
28.9380
29.3010
29.3010
29.3010

Table 13: Comparison of non-dimensional critical buckling load P̄cr =
(rule of mixtures scheme)
l/h Theory
0

CPT [44]
FSDT [44]
RPT [48]
RPT†
Quasi-3D†
0.2 CPT [44]
FSDT [44]
RPT [48]
RPT†
Quasi-3D†
0.4 CPT [44]
FSDT [44]
RPT [48]
RPT†
Quasi-3D†
0.6 CPT [44]
FSDT [44]
RPT [48]
RPT†
Quasi-3D†
0.8 CPT [44]
FSDT [44]
RPT [48]
RPT†
Quasi-3D†
1
CPT [44]
FSDT [44]
RPT [48]
RPT†
Quasi-3D†
†
Present

a/h = 5
n=0
19.2255
15.3228
15.3322
15.3321
15.3629
22.0863
17.6150
18.0422
17.8878
17.7286
30.6685
24.2899
26.1539
25.5457
24.8060
44.9723
34.7856
39.6393
38.2867
36.5415
64.9976
48.2915
58.4862
56.0961
52.8623
90.7444
63.8913
82.6938
78.9675
73.6925

n=1
8.2145
6.8576
6.8611
6.8610
7.3905
9.7879
8.1715
8.3399
8.2820
8.7153
14.5082
11.9922
12.7754
12.5322
12.6741
22.3753
17.9838
20.1658
19.5858
19.2256
33.3892
25.6654
30.5105
29.4240
28.3151
47.5499
34.4981
43.8094
42.0388
39.8872

n = 10
3.8359
2.9979
2.7672
2.7702
3.0118
4.3560
3.4076
3.3619
3.2917
3.4728
5.9164
4.6013
5.0407
4.8371
4.8557
8.5171
6.4804
7.7001
7.3772
7.1597
12.1581
8.9020
11.3322
10.9005
10.3843
16.8393
11.7042
15.9522
15.4071
14.5287

a/h = 10
n=0
19.2255
18.0746
18.0754
18.0756
18.1561
22.0863
20.7607
20.9025
20.8497
20.8583
30.6685
28.7478
29.3735
29.1700
28.9624
44.9723
41.8271
43.4732
43.0329
42.4620
64.9976
59.6657
63.1958
62.4358
61.3467
90.7444
81.8269
88.5416
87.3775
85.6043

38

n=1
8.2145
7.8273
7.8276
7.8277
8.5396
9.7879
9.3241
9.3767
9.3581
10.0344
14.5082
13.7742
14.0232
13.9459
14.5168
22.3753
21.0597
21.7657
21.5846
21.9814
33.3892
30.9928
32.6036
32.2693
32.4199
47.5499
43.3274
46.5372
45.9981
45.8223

Pcr a2
of square FG microplates
E2 h3

n = 10
3.8359
3.5853
3.4969
3.4982
3.8921
4.3560
4.0710
4.0513
4.0246
4.3958
5.9164
5.5151
5.6631
5.5925
5.9066
8.5171
7.8802
8.2906
8.1871
8.4246
12.1581
11.1065
11.9349
11.8036
11.9498
16.8393
15.1152
16.6033
16.4431
16.4819

a/h = 20
n=0
19.2255
18.9243
18.9243
18.9244
18.9675
22.0863
21.7387
21.7771
21.7628
21.7852
30.6685
30.1625
30.3324
30.2773
30.2381
44.9723
44.1369
44.5855
44.4673
44.3258
64.9976
63.5656
64.5348
64.3321
64.0474
90.1804
89.8715
89.4018

n=1
8.2145
8.1142
8.1142
8.1143
8.8639
9.7879
9.6675
9.6815
9.6766
10.4160
14.5082
14.3167
14.3832
14.3626
15.0722
22.3753
22.0292
22.2188
22.1708
22.8320
33.3892
32.7517
33.1882
33.0999
33.6948
47.2914
47.1494
47.6596

n = 10
3.8359
3.7700
3.7450
3.7454
4.1850
4.3560
4.2809
4.2752
4.2677
4.7009
5.9164
5.8102
5.8505
5.8312
6.2486
8.5171
8.3472
8.4589
8.4312
8.8281
12.1581
11.8745
12.1011
12.0666
12.4394
16.7793
16.7376
17.0825

Table 14: Non-dimensional critical buckling load P̄cr =
(Mori-Tanaka scheme)
BC

a/h

l/h

SSSS

5

0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1

n=0
RPT
87.4747
103.6359
152.0215
232.4351
344.7280
488.8378
105.6668
123.5339
177.1295
266.4425
391.4648
552.1931
107.0958
125.0926
179.0825
269.0653
395.0406
557.0082

0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1

20

100

CCCC

5

20

100

†

Dm =

Em h3
2)
12(1 − νm

Pcr a2
†
Dm

Quasi-3D
86.5475
101.8797
147.6312
223.1633
327.6683
460.3907
105.6221
123.4102
176.7715
265.6972
390.1733
550.1813
107.1271
125.1206
179.1009
269.0682
395.0223
556.9633

n=1
RPT
35.0795
42.3816
64.2668
100.6957
151.6424
217.0988
42.0033
49.9504
73.7907
113.5226
169.1451
240.6580
42.5423
50.5370
74.5212
114.4947
170.4576
242.4098

178.2578 188.3478
206.9297 215.7331
292.6287 295.5433
435.0030 424.8837
633.8611 601.6084
889.1171 823.8983
273.9507 288.6976
308.7803 323.9291
413.0833 429.0109
586.6304 603.0356
829.3607 845.5938
1141.2885 1156.5929
283.7646 292.2008
319.1326 327.9091
425.0550 434.6081
601.2918 611.6685
847.7694 858.8291
1164.4998 1176.0787

72.2150
85.4028
124.8630
190.5274
282.3839
400.4310
109.0237
124.5481
171.0330
248.3764
356.5558
495.5793
112.7294
128.4414
175.4917
253.7742
363.2623
503.9643
39

of Al/Al2 O3 square microplates

Quasi-3D
35.6610
42.6185
63.3902
97.7189
145.2904
205.8260
43.4065
51.3272
75.0881
114.6856
170.1141
241.3659
44.0165
52.0102
75.9913
115.9596
171.9154
243.8584

n = 10
RPT
21.7236
25.4909
36.7873
55.6034
81.9336
115.7764
27.5182
31.5521
43.6520
63.8152
92.0401
128.3267
27.9978
32.0486
44.2009
64.4546
92.8097
129.2661

Quasi-3D
21.7189
25.2307
35.7417
53.1891
77.4961
108.5941
27.8578
31.8705
43.9085
63.9707
92.0558
128.1618
28.3566
32.4065
44.5563
64.8061
93.1557
129.6052

78.0421
90.6731
127.4320
187.0650
268.7548
371.8457
117.3906
133.0578
179.8054
257.2793
365.3387
503.9582
119.3117
135.1344
182.4322
260.9696
370.6627
511.5167

42.1220
49.1101
69.9984
104.7011
153.1709
215.3906
70.8926
78.8262
102.5575
141.9833
197.0586
267.7793
74.1644
82.1289
105.9838
145.6726
201.1735
272.4865

45.4860
52.1032
71.1001
101.8084
144.0664
197.7752
75.2177
83.2085
107.0378
146.4815
201.4295
271.8584
77.0743
85.1055
109.1182
148.9836
204.6383
276.0730

Table 15: The first six non-dimensional buckling loads P̄ =
croplates, n = 10, l/h = 0.2 (Mori-Tanaka scheme)
BC
SSSS

a/h Theory
5

10

20

100

CCCC

5

10

20

100

IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)

Mode
1
25.4064
25.2307
25.4909
30.3554
30.2890
30.1012
31.8815
31.8705
31.5521
32.3844
32.4065
32.0486
50.6557
52.1032
49.1101
74.8665
75.4661
70.1001
83.8952
83.2085
78.8262
85.2475
85.1055
82.1289

40

2
47.0158
46.2828
47.9306
67.5835
67.2254
67.3570
76.0319
75.8988
75.2628
79.0722
79.1198
78.2296
66.0350
66.3156
66.7907
114.2196
115.7598
110.6293
140.5084
140.4610
133.9531
149.0723
148.8861
143.8128

P a2
of Al/Al2 O3 square miDm

3
47.0158
46.2828
47.9306
67.5835
67.2254
67.3570
76.0319
75.8988
75.2628
79.0722
79.1198
78.2296
66.0350
66.3156
66.7907
114.2196
115.7598
110.6293
140.5084
140.4610
133.9531
149.0723
148.8861
143.8128

4
62.7947
61.4689
64.7895
101.6321
100.9300
101.9693
121.4313
121.1662
120.4131
129.2976
129.3662
127.9519
78.0114
77.2179
80.3088
146.4170
147.8696
144.4919
191.9601
192.6237
185.0964
210.8794
210.7133
203.9109

5
66.4706
64.8256
68.7456
115.5866
114.4706
116.1512
143.5095
143.0021
142.2468
155.1773
155.2457
153.4678
80.2597
78.9774
82.9464
160.1781
161.3725
158.8413
219.4306
220.3776
211.8925
246.4373
246.2408
237.9639

6
66.4706
64.8256
68.7456
115.5866
114.4706
116.1512
143.5095
143.0021
142.2468
155.1773
155.2457
153.4678
84.0819
84.1043
85.5400
172.6816
174.2442
167.4729
239.0973
239.1104
227.4324
267.8293
267.4391
257.7571

Pcr R2
Table 16: Comparison of non-dimensional critical buckling load P¯cr =
of CCCC
Dm
Al/ZrO2 -2 circular plates (rule of mixtures scheme)
n
0

Theory

TSDT [89]
UTSDT [90]
TSDT [68]
RPT [61]
RPT (Present)
Quasi-3D (Present)
0.5 TSDT [89]
UTSDT [90]
TSDT [68]
RPT [61]
RPT (Present)
Quasi-3D (Present)
2
TSDT [89]
UTSDT [90]
TSDT [68]
RPT [61]
RPT (Present)
Quasi-3D (Present)
5
TSDT [89]
UTSDT [90]
TSDT [68]
RPT [61]
RPT (Present)
Quasi-3D (Present)
10 TSDT [89]
UTSDT [90]
TSDT [68]
RPT [61]
RPT (Present)
Quasi-3D (Present)

h/R
0.1
14.089
14.089
14.1089
14.2023
14.0932
14.8264
19.411
19.413
19.4391
19.5663
19.4169
20.5166
23.074
23.075
23.1062
23.2592
23.0809
24.4332
25.439
25.442
25.4743
25.6418
25.4469
26.8812
27.133
27.131
27.1684
27.3429
27.1395
28.6197

41

0.2
12.574
12.575
12.5914
12.7281
12.5776
13.4557
17.311
17.310
17.3327
17.5180
17.3133
18.6074
20.803
20.805
20.8319
21.0569
20.8088
22.3510
22.971
22.969
22.9992
23.2426
22.9742
24.6195
24.423
24.422
24.4542
24.6994
24.4287
26.1483

0.25
11.638
11.639
11.6540
11.8143
11.6409
12.4564
16.013
16.012
16.0334
16.2506
16.0153
17.2206
19.377
19.378
19.4033
19.6687
19.3812
20.8035
21.414
21.412
21.4407
21.7268
21.4168
22.9303
22.725
22.725
22.7536
23.0389
22.7297
24.3140

0.3
10.670
10.670
10.6842
10.8666
10.6719
11.3775
14.672
14.672
14.6910
14.9381
14.6740
15.7222
17.882
17.881
17.9060
18.2099
17.8848
19.1161
19.780
19.778
19.8043
20.1313
19.7813
21.0878
20.948
20.949
20.9750
21.2986
20.9524
22.3196

Pcr R2
Table 17: Non-dimensional critical buckling load P¯cr =
of Al/Al2 O3 circular miDm
croplates (Mori-Tanaka scheme)
h/R

l/h

n=0
RPT
Simple support
0.1
0
22.5182
0.2 23.0489
0.4 24.1022
0.6 25.0682
0.8 25.7985
1
26.3291
0.2
0
21.7456
0.2 22.2719
0.4 23.2739
0.6 24.2066
0.8 24.9581
1
25.5480
0.3
0
20.5707
0.2 21.0993
0.4 22.0854
0.6 23.0439
0.8 23.8746
1
24.5707
Clamped support
0.1
0
76.5059
0.2 86.3696
0.4 115.9595
0.6 165.2737
0.8 234.3115
1
323.0727
0.2
0
68.2782
0.2 77.6521
0.4 105.7687
0.6 152.6199
0.8 218.1997
1
302.5060
0.3
0
57.9331
0.2 66.6302
0.4 92.7032
0.6 136.1160
0.8 196.8399
1
274.8619

Quasi-3D

n=1
RPT

Quasi-3D

n = 10
RPT

Quasi-3D

22.6953
23.2627
24.4172
25.5337
26.4629
27.2420
22.0263
22.5928
23.7059
24.8190
25.8386
26.7897
20.8606
21.4128
22.4669
23.5641
24.6485
25.7228

9.5368
9.8179
10.3474
10.8125
11.1543
11.3981
9.2059
9.4873
9.9893
10.4423
10.8003
11.0777
8.7033
8.9860
9.4810
9.9511
10.3517
10.6823

9.7960
10.0793
10.6292
11.1391
11.5512
11.8897
9.4510
9.7327
10.2520
10.7545
11.2074
11.6254
8.9036
9.1777
9.6686
10.1693
10.6610
11.1443

5.9574
6.0876
6.3543
6.6130
6.8182
6.9718
5.6942
5.8318
6.1019
6.3676
6.5891
6.7647
5.3041
5.4526
5.7409
6.0352
6.2927
6.5055

6.0309
6.1617
6.4360
6.7163
6.9619
7.1752
5.7510
5.8833
6.1509
6.4377
6.7170
6.9866
5.3287
5.4621
5.7276
6.0281
6.3431
6.6623

80.4859
90.4095
120.0466
169.2564
237.9906
326.2293
73.0452
82.3859
110.0650
155.7501
219.2330
300.3159
61.7636
70.2136
95.0207
135.6444
191.6493
262.6208

30.4539
34.8524
48.0474
70.0379
100.8235
140.4042
27.3245
31.5460
44.2098
65.3154
94.8629
132.8523
23.3438
27.3083
39.2004
59.0182
86.7603
122.4264
42

32.7398
37.1578
50.3612
72.3004
102.9597
142.3316
29.8617
34.0702
46.5462
67.1544
95.8167
132.4584
25.4648
29.3268
40.6571
59.2280
84.8742
117.4429

19.7743
22.0235
28.7616
39.9728
55.6486
75.7879
17.1442
19.3412
25.9122
36.8218
52.0495
71.5905
14.0464
16.1359
22.3875
32.7704
47.2643
65.8628

20.9413
23.1967
29.9236
41.0821
56.6598
76.6551
18.6240
20.7979
27.1925
37.7110
52.3297
71.0373
15.3047
17.3213
23.1508
32.6634
45.8312
62.6312

Table 18: The first six non-dimensional buckling loads P̄ =
n = 1, l/h = 0.6 (Mori-Tanaka scheme)
BC

h/R

Theory

Simple support

0.1

IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)
IGA-Zenkour
Quasi-3D (Present)
RPT (Present)

0.2

0.3

Clamped support 0.1

0.2

0.3

Mode
1
11.4689
11.1391
10.8125
11.3113
10.7545
10.4423
10.8014
10.1693
9.9511
72.6509
72.3004
70.0379
66.8949
67.1544
65.3154
58.9576
59.2280
59.0182

43

P R2
of Al/Al2 O3 circular microplates,
Dm

2
54.8031
54.0884
53.6215
50.7663
49.5358
50.2378
45.2015
43.6917
45.6739
126.0657
125.7477
123.1185
108.9270
108.3685
109.4813
90.1444
88.4285
93.6789

3
54.8031
54.0884
53.6215
50.7663
49.5358
50.2378
45.2015
43.6917
45.6739
126.0657
125.7477
123.1185
108.9270
108.3685
109.4813
90.1444
88.4285
93.6789

4
95.3143
93.7842
93.3944
83.1176
80.6606
83.3889
69.1649
66.3080
71.5215
172.3464
171.5372
169.2944
140.4481
138.0488
143.4150
110.5155
106.6289
117.0824

5
123.2488
122.0854
122.3489
107.1987
104.8026
109.1766
89.4125
86.1911
93.8041
205.4419
205.0589
203.3425
166.8569
164.4947
171.7706
131.5270
127.0395
140.3857

6
129.6396
128.1298
128.3281
110.9343
108.1436
112.9355
90.9795
87.4685
95.5819
222.4846
221.0537
218.6139
175.2599
172.3137
179.5962
134.2574
129.6436
143.1361

z
Ceramic

h/2

h/2

b

y

x
Metal
a

Figure 1: The functionally graded microplate model.

44

Figure 2: The effective modulus of Al/Al2 O3 plates according to the rule of mixtures scheme
(in solid lines) and Mori-Tanaka scheme (in dash lines).

45

1

0
0 ,0 ,0 ,0

1/5

3/5,3/5,3/5

2/5

4/5

1 ,1 ,1 ,1

1 2 3 3 3 4
(a) Cubic basis functions corresponding to Ξ = {0, 0, 0, 0, , , , , , , 1, 1, 1, 1}.
5 5 5 5 5 5

(b) Bivariate B-spline basis functions.

Figure 3: One- and two-dimensional B-spline basis functions.

46

0 .2 4 6
0 .2 4 4
0 .2 4 2
A n
p =
p =
p =

w

0 .2 4 0

a ly tic a l
2
3
4

0 .2 3 8
0 .2 3 6
0 .2 3 4
0 .2 3 2
0 .2 3 0
3

5

7

9

1 1

1 3

1 5

1 7

N u m b e r o f e le m e n ts p e r e d g e

Figure 4: Convergence study of non-dimensional central deflection of homogeneous square
microplates.

47

h
q0
a

z
y
x
a
(a) Geometric configuration.

(b) Control point net and 11×11 cubic elements.

Figure 5: Geometry and element mesh of a square microplate.

48

1 .0

0 .8

l/h = 0
l/h = 0.4
l/h = 1.0

Quasi-3D
RPT

w

0 .6

0 .4

0 .2

0 .0
0

2

4

6

8

1 0

n
Figure 6: Effects of material index n and material length scale ratio l/h on the central
deflection of CCCC Al/Al2 O3 square microplates subjected to uniformly distributed load,
a/h = 5 (rule of mixtures scheme).

49

(a) CCCC.

(b) SSSS.

(c) SCSC.

(d) SFSF.

Figure 7: Deformed configuration of Al/Al2 O3 square microplates.

50

2 0
Quasi-3D
RPT

l/h = 0
l/h = 0.4
l/h = 1.0

ω

1 5

1 0

5

0
0

2

4

6

8

1 0

n

Figure 8: Effects of material index n and material length scale ratio l/h on the natural
frequency of CCCC Al/Al2 O3 square microplates, a/h = 5 (rule of mixtures scheme).

51

(a) ω̄1 = 15.4413.

(b) ω̄2 = 29.3267.

(c) ω̄3 = 29.3267.

(d) ω̄4 = 41.5955.

(e) ω̄5 = 48.2657.

(f) ω̄6 = 48.4504.

Figure 9: The first six free vibration mode shapes of Al/Al2 O3 CCCC square microplates.
52

z
y
x
R

h

(a) Geometric configuration.

(b) Control point net and 11×11 cubic elements.

Figure 10: Geometry and element mesh of a circular microplate.

53

(a) ω̄1 = 7.3195.

(b) ω̄2 = 14.0950.

(c) ω̄3 = 14.0950.

(d) ω̄4 = 21.3095.

(e) ω̄5 = 21.7503.

(f) ω̄6 = 23.9885.

Figure 11: The first six free vibration mode shapes of clamped Al/Al2 O3 circular microplates.
54

(a) P̄1 = 52.1032.

(b) P̄2 = 66.3156.

(c) P̄3 = 66.3156.

(d) P̄4 = 77.2179.

(e) P̄5 = 78.9774.

(f) P̄6 = 84.1043.

Figure 12: The first six buckling mode shapes of Al/Al2 O3 CCCC square microplates.
55

3 0
Quasi-3D
RPT

l/h = 0
l/h = 0.4
l/h = 1.0

2 4

P

c r

1 8

1 2

6

0
0

2

4

6

8

1 0

n

(a) Simply support.

3 5 0
Quasi-3D
RPT

l/h = 0
l/h = 0.4
l/h = 1.0

3 0 0
2 5 0

P

c r

2 0 0
1 5 0
1 0 0
5 0
0
0

2

4

6

8

1 0

n

(b) Clamped support.

Figure 13: Effects of material index n and material length scale ratio l/h on the critical
56
buckling load of Al/Al2 O3 circular microplates, h/R = 0.2 (Mori-Tanaka scheme).

(a) P̄1 = 10.7545.

(b) P̄2 = 49.5358.

(c) P̄3 = 49.5358.

(d) P̄4 = 80.6606.

(e) P̄5 = 104.8026.

(f) P̄6 = 108.1436.

Figure 14: The first six buckling mode shapes of simply-supported Al/Al2 O3 circular microplates.
57

