Composites Part B 163 (2019) 779–793

Contents lists available at ScienceDirect

Composites Part B
journal homepage: www.elsevier.com/locate/compositesb

Influence of material uncertainties on vibration and bending behaviour of T

skewed sandwich FGM plates
Sanjay Singh Tomara, Mohammad Talhab,∗
a
Aerospace Engineering Department, Indian Institute of Technology Kanpur, Uttar Pradesh, 208016, India
b
School of Engineering, Indian Institute of Technology Mandi, Himachal Pradesh, 175005, India

ARTICLE INFO ABSTRACT

Keywords: Present study aims to investigate the influence of material uncertainties on vibration and bending behaviour of
Functionally graded materials skewed sandwich FGM plates. Reddy's higher order shear deformation theory has been employed to model the
Material uncertainties displacement field. Variational approach has been used to derive the governing differential equations. Effect of
Perturbation technique material uncertainties in the formulation have been incorporated using first order perturbation technique
Monte Carlo simulation
(FOPT). An efficient stochastic finite element formulation (SFEM) have been used for the calculation of first and
Finite element method
second order statics of natural frequency and transverse deflection. Validation of the results have been per-
formed with the help of available literature and separately developed Monte Carlo formulation (MCS) algorithm.
A large number of examples have been solved to quantify the effect of uncertainties on the vibration and bending
characteristics of functionally graded skew sandwich plates.

1. Introduction manufacturing with full accuracy. Hence the inclusion of the deviation
of material and geometric properties in the structures are inherent. This
Laminated composites have found their wide application in various stochastic nature of material properties affects structural response of
engineering field's. Basically laminated composite structures constitute the system. Hence it becomes mandatory to involve such behaviour in
of various isotropic layers stacked together by means of some bonding the analysis. In recent years, this field have attracted the various re-
materials. The major drawback in this construction is the presence of searchers across the globe and as a result of which a large number of
stress concentrations at the interface of the various layers. Due to which literature have been available in this field.
laminated composites go through the delamination, cracking and other Nakagiri et al. [4] studied vibration eigenvalue problem for FRP
damages. This limitation of the laminated composites can be overcome laminated plates with the help of second order perturbation method
with the help of functionally graded material (FGM's) [1]. FGM is an while considering various stacking sequence and layer thickness of the
advance class of materials which have excellent thermo-mechanical laminate. Oh and Librescu [5] studied the effect of uncertainties in the
properties which enables them to maintain their structural integrity in properties on the vibration and reliability of cantilever beam. In order
higher thermal loading conditions. Need of such material arises since to account the stochastic behaviour of uncertainties first and second
the second world war when scientist started searching a material which order methods had been employed on the eigenvalues. Vanini and
can withstand higher temperature conditions [2]. They are usually Mariani [6] investigated the effect of uncertainties on the vibration
composed of metal and ceramics. These constituent materials are behaviour of composite plates resting on elastic foundation by solving
graded in a fashion following various laws such as power law, sigmoid the vibration problem in both deterministic and stochastic cases. Raj
law, exponential law etc. The graded variation of FGM structures pre- et al. [7] calculated the static response of composite plates with the
vents them from delamination unlike the laminated composite mate- help of MCS and finite element methodology. With the help of this
rials [3]. analytical investigations authors predicted the behaviour of normally
The process of fabrication of FGM is a complex and challenging as a distributed parameters. Manohar and Ibrahim [8] reviewed the various
large number of design variables are involved which enables the work published in the area of vibration analysis of strictures having
smooth variation of material properties in the structures. In practical material uncertainties. Author covered various algorithm and methods
cases it is very difficult to have complete control on each step of such as finite element, boundary element method along with various

∗
Corresponding author.
E-mail address: [email protected] (M. Talha).

https://doi.org/10.1016/j.compositesb.2019.01.035
Received 23 April 2018; Received in revised form 2 January 2019; Accepted 3 January 2019
Available online 10 January 2019
1359-8368/ © 2019 Elsevier Ltd. All rights reserved.
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

statistical techniques such as MCS and statistical energy analysis. employed for the validation purpose. The effect of material un-
Singh et al. [9] investigated the vibration response of laminated certainties on the deflection statics of sandwich composite plates have
composite plates having material uncertainties. They used HSDT to been investigated by Pandit et al. [26]. Authors developed an improved
account the effect of shear deformation and rotary inertia of the plates. theory to calculate the structural response of sandwich plates. Theory
FOPT had been used to model the stochastic nature of the system. The was assumed to follow the stress continuity at the layer interfaces.
results obtained have been validated with MCS and compared with Shahabian and Hosseini [27] studied the dynamic response of FG
various plate models. Same group of authors performed another study hollow cylinder having the material uncertainties. Authors employed
[10] on the stochastic vibration of laminated composite cylindrical Navier solution to calculate the displacement at each point. Whereas
panel's. Considering uncertainties in the material properties. They Newmark method with MCS had been employed to calculate mean and
employed the effect of transverse and rotary inertia in the calculation variance of the deflection statics. Same group of authors [28] per-
whereas FOPT had been employed to account the effect of uncertainties formed one more study on the stress calculation on the FG hollow cy-
in the system. linder under the dynamic shock loading. The material properties of
Onkar and Yadav [11] studied deflection response of laminated Al/ Al2 O3 were considered as basic random variables following Gaussian
composite plates under random transverse loading conditions. They distribution. Chandrashekhar and Ganguli [29] used MCS with Latin
employed laminated plate theory along with von-karman assumptions hypercube sampling algorithm to perform the non-linear vibration
for the formulation. FOPT had been used to calculate the transverse analysis of laminated composite and sandwich composite plates. Pandit
deflection of the problem whereas MCS had been used for the validation et al. [30] proposed an improved higher order zigzag theory and used it
purpose. Same group of authors performed two more studies on lami- to calculate the vibration behaviour of soft core sandwich composite
nated composites. Firstly they [12] calculated the non-linear vibration plates. Later same group of author's [31] performed another study using
response of laminated composite plates. Later they [13] calculated the improved zigzag theory calculating the vibration response of sandwich
forced vibration response of laminated composite plates. They used composite plates. Theory relaxes the assumption of layerwise theory
Hamilton's principle to calculate the governing equation's for composite while considering the advantages of single later theories. They used
plates. Closed form solution along with FOPT had been used to calcu- FOPT to account the effect of material uncertainties.
late the variance of various responses. Lal and Singh [32] analysed the bending response of laminated
Yang et al. [14] performed stochastic bending response of FGM composite plates under thermo-mechanical type of loading conditions.
plates having uncertain material properties. They used HSDT in con- They employed a C 0 continuous finite element model along with FOPT
junction with FOPT to calculate the first and second order moment of to calculate deflection statics of the plates loaded under transverse
deflection statics. Kitipornchai et al. [15] studied the vibration response loading conditions. Same group of authors [33] performed another
of FG laminates having uncertain material properties in the thermal study on the spherical shell panel considering material uncertainties.
environment. They used HSDT for the formulation whereas, semi-ana- They used HSDT in Von-Karman sense for formulation. Shell was as-
lytical approach to solve the governing differential equations. FOPT had sumed to be under hygro-thermo-mechanical loading. Giunta et al. [34]
been used to calculate the variance of the eigen solution of the problem. performed free vibration analysis on the composite plates having ma-
Giunta and Carrera [16] proposed a Newton series approximation in the terial and geometric uncertainties. Author's used classical and two di-
finite element environment analysis. In order to prove the improved mensional model via Carrera's unified formulation to formulate the
computational efficiency with the presented formulation author solved eigenvalue problem. Monte Carlo method was used to incorporate the
the truss problem under static loading conditions. Tripathi et al. [17] uncertainties. They prove that the frequency is not sensitive for the
analysed the stochastic vibration response of laminated composite materialistic uncertainty as compared to the geometric uncertainty.
conical shell structures considering the material properties as basic Jagtap et al. [35] analysed the non-linear stochastic bending response
basic random variables. Effect of random material properties and of FGM plates in thermal environment. Authors used a C 0 continuous
random excitation on the reliability and vibration response of rotating HSDT in conjunction with Von-Karman nonlinear assumptions.
beam had been studied by Hosseini and Khadem [18]. Thermal analysis had been done while considering temperature de-
Gao [19] performed vibration analysis on the truss structure using pendent and independent material properties. The noninear vibration
random factor and interval factor methods. Author considered the fre- analysis of composite plates laminated with piezoelectric layer having
quency response and mode shapes to be random variable in analysis. material uncertainty had been dome by Lal et al. [36].
Algebra synthesis method for random variables had been employed to Sepahvand et al. [37] used polynomial chaos free expansion to
derive the first and second order statics of the system as Rayleigh calculate the effect of material uncertainty in the elastic modulus on the
quotient. Lal et al. [20] investigated the stochastic bending behaviour free vibration response of orthotropic plates. Shegokar and Lal [38]
of laminated composite plates resting on two parameter Pasternak performed a study on the deflection statics of the FGM beam beam
foundation. They used a C 0 continuous finite element model in con- bonded with piezoelectric layers having material uncertainties. Plate
junction with FOPT to calculate the mean and variance of deflection was assumed to be under thermo-electro-mechanical loading. They
statics. Same group of authors performed few more studies. Firstly they employed a direct iteration based FOPT to calculate the second order
[21] calculated the vibraiton response of laminated composite plates statics of the problem. Chang [39] investigated the non-linear dynamic
resting on elatic foundation. Later they performed two nonlinear stu- behaviour of single walled carbon nano-tube using SFEM technique
dies for calculation the nonlinear vibration [22] and nonlinear bending considering the effect of uncertainty in elastic moduli and mass density.
[23]of laminated composite plates resting on elastic. In both studies Kumar et al. [40] studied the stochastic flexural response of lami-
author used Von-Karman non-linear strain assumptions to account the nated composite plates resting on elastic foundation having cubic
non-linearity in the system. Plate was assumed to be under uniform nonlinearity. They modeled various material properties, lateral load,
pressure and thermal loading. hygroscopic expansion coefficient as random variable. Dash and Singh
Shaker et al. [24] used stochastic FEM methodology to calculate the [41] calculated the stochastic banding response of laminated composite
vibration response of FGM plates. They used first and second order plates. They used used a shear and normal deformation theory to ac-
reliability approaches to calculate the variance of the eigenvalue solu- count the effect of large deformation's. Green-Lagrange nonlinear as-
tions of the vibration problem. Singh et al. [25] presented a study on sumptions have been considered to account the non-linearity in system.
the non-linear bending response of laminated composite plates resting A C 0 continuous finite element formulation along with perturbation
on two parameter elastic foundation. Authors employed HSDT in von- technique had been employed to calculate deflection statics of the
karman sense for the formulation of the problem. Plate was assumed to composite plates. Talha and Singh [42] used an improved kinematics
be under the uniformly distributed loading conditions, MCS have been for the stochastic vibration analysis of FGM plates. Theory considers the

780
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 1. Geometry of plate in Cartesian coordinate system [59].

effect transverse and normal shear deformation's. In order to consider has been used to quantify the uncertainties. They proved that the
the effect of stochastic nature in the formulation material properties sources of uncertainties effects the low velocity impact of the FGM
such as elastic modulus, volume fraction index of FGM plates had been plates.
considered as basic random variables. The present analysis aims to investigate stochastic bending and free
Xu et al. [43] employed random factor method to calculate sto- vibration response of functionally graded skew sandwich plates. The
chastic dynamic characteristics of functionally graded beam and formulations have been performed with Reddys higher order shear
showed the characteristic effect of constituent volume fraction index on deformation theory (HSDT) by employing isoparametric C 0 continuous
the dynamic characteristics distribution. Asnafi and Abedi [44] in- finite element. The governing equation's have been derived using var-
vestigated the dynamic stability and bifurcation of FGM loads under iational approach. FOPT has been used to account the effect of material
stochastic random loading conditions. The properties of the plates were uncertainties in the functionally graded sandwich plate. Comparison
assumed to be varying in accordance with power law, exponential law and convergence studies have been performed to show the effectiveness
and sigmoid law. They employed probability density function to eval- and reliability of the present formulation. Separately developed Monte
uate the instability region and the response of the FGM plates had been Carlo simulation methods for validation purposes. Various results have
solved using exact solution of Fokker Planck Kolmogorov equation. been shown to demonstrate the effect of skew angles, aspect ratio,
Chakraborty et al. [45] used PCFE (polynomial correlated function thickness ratios, volume fraction index, and boundary conditions on the
expansion) approach to perform the free vibration stochastic analysis second order statics of defection and fundamental natural frequencies.
on the composite plates. Garca-Macas et al. [46] studied the free vi-
bration response of CNT reinforced FGM plates. The analysis had been 2. Material properties and plate geometry definition
performed with finite element methodology and probability theory was
used to incorporate the uncertainties in the material properties. Consider a rectangular sandwich FGM plate shown in Fig. 1 with
Mukhopadhyay et al. [47] calculated the effect of noise on sto- sides a, b and height h respectively. Assuming that the mid-plane of
chastic free vibration behaviour composite laminated shallow shell plate coincides with the x-y plane.
structures. Authors used surrogate based methodology to quantify The gradation of FGM layer follows the power law distribution. The
material uncertainties. Dey et al. [48] presented a Kriging model to effective material properties are calculated using rule of mixtures
analyse the free vibration response of composite shallow doubly curved (Vc + Vm = 1),
shells. They used Mindlin's first order plate theory for formulation.
Peff = Pc Vc + Pm Vm (1)
Later same group of authors performed many studies which can be seen
in Refs. [49–53]. Corradi et al. [54] performed a large number of ex- Where, Peff represents the effective material properties such as young's
periments to calculate capacity and strength uncertainty of various modulus, Poisson's ratio, mass density and coefficient of elastic ex-
types of wood beams. A special attention had been given to timber pansion, whereas Vc , Vm represents the volume fraction of the ceramic
beams as they are commonly used for construction purpose. Tomar and metal, respectively.
et al. [55] presented a critical review on the non-deterministic meth- Plate having FGM core with homogeneous face sheet shown in Fig. 2
odologies for the analysis of the composite structures. They performed a considered to follow given volume fractions [61],
thorough survey on the various on the research work which had been
Vc = 1, z [h1, h/2]
performed in the last couple of centuries in the stochastic domain.
( ) ,z
n
Wu et al. [56] performed the static analysis on the Euler Bernoulli Vc =
z h2
[h2, h1], hf = h1 + h2
hf
FGM beam while considering the uncertainties in various governing
parameters such as geometry, material properties and loading's. Per-
VC = 0, z [ h /2, h2] (2)
turbation theory along with finite element methodology had been used In order to consider the effect of thermal environment material
to model the uncertainties and calculating second order deflection properties are assumed to be the function of temperature.
statics of the beam. Karsh et al. [57] studied the stochastic dynamic
response of the twisted FGM plates using neural network based ap- P (t ) = P0 (P 1 T 1 + 1 + P1 T + P2 T 2 + P3 T 3), (3)
proach integrated with finite element formulation. They considered the Where, P0, P 1, P1, P2, P3 represents coefficient of temperature and are
uncertainties in the material and geometric properties of the plates and unique to the type of materials. A linear temperature distribution has
showed that the uncertainties significantly affect the dynamic beha- been employed across the thickness of the plate which is written as,
viour of the functionally graded twisted plates. Same group of authors
presented another study [58] on the FGM plates considering the low z 1
T (z ) = Tb + T +
velocity response. They employed power law to calculate the material h 2 (4)
properties whereas the modified Hertzian contact law had been used to Where, T = Tt Tb , Tt and Tb represents the temperature at top and
calculate the contact forces. Fuzzy logic based non-probabilistic method bottom face of the functionally graded sandwich plate.

781
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 2. Cross section of Sandwich FGM plate with FGM core and homogeneous face sheets [60].

W0 W0
x = 0, y =0
x y (7)

Final displacement vector at a node ‘i’ can be written as,

{ i} = {U0, V0, W0, T

x, y, x, y} (8)

3.2. Strain displacement relationship

Fig. 3. Nide noded isoparametric element in local coordinate and global co-
ordinate system. The strain displacement relationship for sandwich FGM plate can be
represented as,
9 9 9
{ }= Ni { }i; x= Ni x i ; y= Ni yi ={} { th} (9)
i=1 i=1 i=1 (14)
Where, ε and t represents the linear strain vector and thermal strain
3. Mathematical formulation vector.
0
3.1. Displacement field x
1 1
1 0 3
1
0
y 2 1 0 3
2 2 2
The origin of Cartesian coordinate system is assumed to be lying on = yz = 3
0
+z 0 + z2 1 + z3 0 ;
the edge of the sandwich plate shown in Fig. 1. Reddy's higher order zx 0
4
0 2
2 0
shear deformation theory has been adopted for the mathematical for- xy
0
1
0 3
5
3 3
(10)
mulation. As this theory includes the effect of transverse shear de-
formation and rotary inertia. According to Reddy's theory displacement Where
field at a point can be written in terms of third order polynomial,
0 U0 0 V0 W0 W0 V0 U0
1 = , 2 = , 30 = y + , 40 = x + , 50 = + ,
U U0 x x x
x y y x x y

V = V0 + z y + z2 y + z3 y
1
1 = x
, 1
2 =
y
, 1
3 = x
+
y
, 2
1 = 3C2 ( + y ),
2
2
x y y x y
W W0 0 0 0 (5)
= 3C2 ( x + x ),
Where, U,V and W represents inplane and transverse displacement at a
point, U0, V0, W0 are the displacements corresponding to mid-plane. y
3
1 = C2 ( x
x
+ x
x
), 3
1 = C2 ( y
y
+
y
y ), 3
1

and x represents the rotation of the normal to the mid-plane about y

and x axis, x , y , x , y are the higher order terms of Taylor's series = C2 ( y
x
+ x
y
+ x
y
+ x
y )
expansion. Value of the higher order terms are calculated with the help
of zero transverse stress conditions at top and bottom surface of the Thermal strain vector can be represented as,
plate. After applying these boundary condition following C1 continuous
1
displacement field is obtained, 2
W0 th = 0 T
U U0 x x 0
V = V0 + f1 (z ) y + f2 (z ) W0 0 (11)
y
W W0 0 0 (6) Where, 1 and 2 represents the coefficient of thermal expansion in in-
plane directions respectively.
Where, f1 (z ) = C1 z C2 z 3,
f2 (z ) = C3 C1 = 1, C2 = C3 = z 3, 4/3h2
The complexities associated with the analysis of C1 continuous finite
element is well known. To overcome these difficulty a C 0 continuitous 3.3. Constitutive relationship
displacement field have been prepared with panalty approch, due to
this some artificial constants are imposed to the displacement field The constitutive relation depicts the relationship among stress and
[62], strains of the sandwich FGM plate and is written as,

782
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

ne ne
xx Q11 Q12 0 0 0 1
U= U (e ) = { } e (T ) [K (e) ]{ } e
yy Q12 Q22 0 0 0 2 (19)
e=1 e=1
yz = 0 0 Q44 0 0 { l t }or{ i} = [Q¯ ij]{ i}
zx 0 0 0 Q55 0 Where, ‘ne’ represent the number of elements present in the domain of
the plate and K (e) be the element stiffness matrix for the plate can be
xy 0 0 0 0 Q66 (12)
written as,
Where , [Q¯ ] and [ ] represents the stress vector, reduced elastic coefficient
matrix and strain vector respectively. The terms in the elastic coefficient matrix K (e) = [B]Ti [D][B ]i (20)
can be written as.
Where,
E (z , T ) (z , T ) E (z , T )
Q11 = Q22 = , Q12 = , Q44 = Q55 = Q66 D = [T ]Ti [Qij][T ]i
1 (z , T )2 1 (z , T ) 2
E (z , T )
=
2(1 + (z , T )) (13) 4.1.2. Kinetic energy
Material properties in matrix [Qij ] are graded according to equations The kinetic energy of the plate is given by,
(1)–(3) respectively. NL
1 1 hk + 1
T= {X}T {X} dv = {X}T {X} dz dA
2 2 A
k=1
hk
(21)
4. Finite element implementation v

Where, ρ and X represents the mass density and the velocity vector for
In the present finite element formulation C 0 continuous 9 noded the plate respectively. The velocity vector can be written as,
isoparametric Lagrangian element has been employed as shown in
{X} = [N¯ ]{ } (22)
Fig. 3. Each node of the element contains 7 degree of freedom (DOF)
per node and hence 63 DOF per element. The generalized displacement Where, [N̄] is the function of the thickness coordinate can be written as,
vector and element geometry at a point can be written as [63],Where, i
represents the displacement vector at ith node and x i , yi are the Cartesian 1 0 0 f1 (z ) 0 f2 (z ) 0
coordinates at ith node. Ni represents the interpolation functions which [N̄] = 0 1 0 0 f1 (z ) 0 f1 (z )
can be written as [64], 0 0 1 0 0 0 0 (23)

N1 = 4 (
1 2 )( 2 ), N2 = 4 (
1 2 + )( 2 ) The kinetic energy of FGM laminate will be summation of each
1 2 1 2
element. By applying the finite elements approximation the global ki-
N3 = ( + )( 2 + ), N4 = ( )( 2 + )
4 4 netic energy expression can be written as,
1 2 )( 2 1 2 2)
N5 = (1 ), N6 = ( + )(1 ne
2 2 1 T
1 2 )( 2 1 2 2)
T= { } (e) [M ](e) { } (e)
N7 = 2
(1 + ), N8 = 2
( )(1 2 e=1 (24)
N9 = (1 2 )(1 2)
(15) Where, m represents the elemental mass matrix for element ‘i’.

4.1. Energy calculation 4.2. Work done due to transverse load

4.1.1. Strain energy The work potential generated due to distributed force acting on the
The strain energy for FGM plate element ‘i’ is given by Ref. [65], top surface of the plate,

1 1 Wp = { }T {q} dxdy
Ue = { }Ti { }i dv = { }Ti [Qij]{ }i dv
2 v 2 v (16) A

Sandwich FGM plate consist of more than one number of layers, due Where, ‘q’ represents distributed load vector acting on top surface of the
to which for the calculation of strain energy material property matrix plate.
[Q] has to be calculated for each layer of sandwich plate and summed
for all ‘NL’ nember of layers. New strain energy expression can be 4.3. Skew boundary transformation
written as,
Considering a skew plate shown in Fig. 4 whose boundaries are
1
NL
hk + 1 supported on the two adjacent edges, it is not necessary condition that
Ue = { l}Ti [Qij]k { l}i dz dA the edges of the boundaries are parallel to the axis of the global co-
2 A hk
k=1 (17) ordinates. This fact leads us to define of the nodal displacements in
terms of local coordinates and then transform it in terms of the global
[T ] is the matrix of z extracted from strain terms,
coordinates in order to define the boundary conditions. The transfor-
{ }i = [T ]{¯}i = [T ][B]i { }i (18) mation of the coordinates followed by the simple rules expressed as
[66],
Where, [T] is the matrix made by extracting z terms from strain matrix
and is given by, i = Tg l
i (25)

1 0 0 0 0 z 0 0 0 0 z30 0 Where, i , are the generalized and global displacement vectors at the
l
i
0 1 0 0 0 0 z 0 0 0 0 z3 0 respective i node. These can be defined as,
T= 0 0 1 0 0 0 0 0 z2 0 0 0 0 T
{ i} = {U0, V0, W0, x, y, x, y}
0 0 0 1 0 0 0 0 0 z2 0 0 0 T
0 0 0 0 1 0 0 z 0 0 0 0 z3 { li} = U0l , V0l , W0l , l l l l
x, y, x, y
(26)
The overall strain energy of the FGM sandwich plate will be sum-
mations of the stain energy of each individual element. The transformation matrix at node i can be written as,

783
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 4. Orientation of skew edges of plate. Fig. 6. Comparison of COV ( 2 ) of Al FGM Al2 O3 sandwich plate with MCS
considering Ec as random variable.

Table 1
Material properties of various constituent materials.
Material E (N /m2) Properties ν ρ(kg/m3) Table 3
Variation of non-dimensional central deflection of Al FGM Al2 O3 sandwich
Al 70 × 109 0.30 2707 plate with n = 4.
Al2 O3 380 × 109 0.30 3800
Mesh size Thickness ratio (a/h)
ZrO2 151 × 109 0.30 3000
Ti 6Al 4V 105.7 × 109 0.298 4429 4 10

4 1.538 1.1666
9 1.1968 0.9146
Table 2
16 1.0651 0.8239
Variation of non-dimensional frequency parameter of Al FGM Al2 O3
25 1.0656 0.8225
sandwich plate with n = 1, a/h = 10. 36 1.0596 0.8205
Mesh size Frequency parameter Neves et al. [69] 1.1108 0.8700

1:1:2 1:2:1 (%) differ. 4.6095 5.6857

4 1.4024 1.4966
9 1.3596 1.4538
16 1.3485 1.4429
25 1.3446 1.4391
36 1.3429 1.4375
El Meiche et al. [68] 1.3332 1.4394

(%) Difference 0.7280 0.1304

Fig. 7. Comparison of COV of transverse deflection of Ni / Al2 O3 FGM plate.

Fig. 5. Comparison of COV of square of natural frequency Al/ ZrO2 FGM plate.

784
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

elemental stiffness and mass matrices as,

KT = TTf K (e) Tf (29)

MT = TTf M (e) Tf (30)

Where, MT , KT represents the transformed global mass and stiffness

matrices.

4.4. Governing equations

In order derive the governing equations of plate variational prin-

ciple have been adopted. The generalized governing differential equa-
tion for combined bending and vibration case can be written as,

[MT ]{q¨} + [KT ]{q} = {F } (31)

Assumption has been made that the system is vibrating in the

principle mode, the eigenvalue problem for bending and vibration re-
sponse of functionally graded skew sandwich plate can be written as,
Fig. 8. Comparison of COV of transverse deflection of Al FGM Al2 O3
[KT ]{q} = [MT ]{q} + {F } (32)
sandwich plate with MCS and FOPT.
Where, = 2 , ω represents the natural frequency of the functionally
c s 0 0 0 0 0 graded skew sandwich plate, q represents the global displacement
s c 0 0 0 0 0 vector, whereas F represents the global transverse load vector.
0 0 1 0 0 0 0
Tg = 0 0 0 c s 0 0
0 0 0 s c 0 0 4.5. Perturbation approach
0 0 0 0 0 c s
0 0 0 0 0 s c (27) In the present study in order to account the effect of material un-
certainties in the vibration and bending response of functionally graded
Where, c = cos( ) and s = sin( ), represent the skew angle for the skew sandwich plates first order perturbation approach have been
plate. From the above transformation matrix it may be noted that for adopted. The basic material properties of skew sandwich plates are
the various nodes which are not present on the skew edges transfor- considered as basic input random variables. In the framework of per-
mation matrix have nonzero value only for the principle diagonal ele- turbation technique, let following random field variable's have been
ments, which are equal to unity. Thus the transformation matrix for a assumed,
complete element the element transformation matric can be written as,
(x , y ) = { 1 (x , y ), 2 (x , y ), 3 (x , y ), ………… m (x , y )} (33)
Tf = Tg [I9] (28)
Here elastic constants such as youngs modulus, poisons ratio and
Where, I9 represents the identity matrix of order 9 × 9, the order of the volume fraction index are assumed as independent random vari-
identity matrix will be same as that of number of nodes present in the able. Using finite element considerations, variable's in the random field
element. Elemental transformation matrix will be used to transform the can be written as,

Fig. 9. Variation of COV of square of natural frequency with thickness ratio(a/ h ) considering (a) Ec (b) Em as random variable.

785
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 10. Variation of COV of square of natural frequency with volume fraction index(n) considering (a) Ec and Em (b) c and m as random variable.

nn m m
0 0 r r
= Ni i Var ( i ) = i, j i, k Cov ( j, k)
i=1 (34) j =1 k=1 (40)
As, material properties have been considered as random variables. where, Cov ( represents the coefficient of variance between
r
j,
r
k)
r
j
Hence vectors and matrices in eq. (3.38) will be function of and will and rk respectively.
be random in nature.
m m
[KT ] = [KT ( )] Var (q) = diag (qi, j )(qi, j )T Exp r
j,
r
j
{qi } = {qi ( )} j =1 j=1 (41)
{ } = { ( )} (35)
Where, Exp ( r
j,
r
j) represents the expectation matrix which can be
In FOPT random variable are expanded using the Taylor series ex- calculated as,
pansion and upto first order terms of expansions are retained. Here, r r
Exp ( j, k) = [ ][ ][ ]
variables in the random field are expanded about mean values 0 of
random variables, and zero mean random part r [67]. and, [ ] and [ ] are standard deviation and correlation coefficient
matrices respectively.
[KT ] = [KT0] + [KTr ]
{qi} = {qi0} + {qir } 0 ... 0 1 12 ... 1k

{ i} = { i0} + { ir } (36) 0 ... 0 21 1 ... 2k

[ ]= . . ... . , [ ] = . . ... .
. . ... . . . ... .
Considering the variation in random part be very small as compared 0 0 ... k1 k2 ... 1
to mean value. Using the above equation and putting in equation (32)
while collecting terms with same order and neglecting terms above first
order following two equations can be derived, 5. Results and discussions

• Zeroth order equation In this section, stochastic vibration and bending response of func-
[KT0]{qi0} 0 0 tionally graded skew sandwich plates various numerical results have
= {F } + i [MT ]{qi } (37)
been solved. This section have been divided into two major parts
namely convergence and comparison study and the parametric studies.
• First order equation Stochastic results have been demonstrated in terms of coefficient of
[KT0]{qir } + [KTr ]{qi0} = {F } + 0 r
i [MT ]{qi } + r 0
i [MT ]{qi } (38) variance (COV=SD/mean) of square of natural frequencies and de-
flection statics of skew sandwich plates against COV of random vari-
Using first order Taylor series expansion, for small values of random ables (R.V.).
part, terms [KTr ], ir , {qir } can be written as,

i
r
=
m 0 r 5.1. Convergence and comparison studies
j = 1 i, d j
m
qir = q0 r
j = 1 i, j j In this section accuracy and applicability of the present formulation
m
[KTr ]i = [KT0] rj has been demonstrated through various examples. The current ap-
j=1 i, j (39)
proach has been validated against the results available in the published
In the current case, results for vibration and bending case are cal- literature and using Monte Carlo simulation (MCS) technique. This
culated separately. The variance for eigenvalues of frequency and de- section aims to demostrate appropriate mesh size which can be used to
flection can be written as, obtain the results.

786
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 11. Variation of COV of square of natural frequency with skew angle (ψ) considering (a) Ec (b) Em (c) c (d) m (e) Ec , Em (f) Ec ,Em , c , m ,n as random variable.

5.1.1. Vibration analysis fractions. Author's used a new hyperbolic shear deformation theory
Example 1: This example presents the validation study for fre- along with the Navier's closed form solution to calculate vibration re-
quency parameter ( ¯ ) of sandwich FGM plates simply supported (SSSS) sponse of the plate. The elastic modulus of metal ( Al ) and ceramic
at all edges have been validated against the results of El Meiche et al. ( Al2 O3 ) have been considered to be 70GPa and 380Gpa , respectively.
[68]. Table 2 shows that the results obtained by the current formulation The material properties are also given in Table 1. The thickness ratio
agrees well with the literature El Meiche et al. [68] at various volume (a/ h ) has been taken as 10, frequency parameter ( ¯ ) used in study can

787
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 12. Variation of COV of square of natural frequency considering (a) all properties as R.V. (b) Various boundary conditions.

be written as: Example:2 This example presents the validation of deflection statics
of Ni /Al2 O3 FGM plate having volume fraction index(n) and thickness
¯ = (b2 / h) 0 / E0 (42) ratio(a/ h ) to be 2 and 10 respectively. Elastic modulus of metal (Em ) is
where, ω, b and h are natural frequency, width and thickness of plate. 0 considered to be independent random variable. Fig. 7 shows that results
and E0 are constant terms in frequency parameters having values obtained with current formulation agrees fairly well with the literature
1Kg / m3 and 1GPa , respectively. [14]. Yang et al. [14] used semi-analytical method based HSDT to
Example:2 In this example validation of coefficient of variation calculate transverse deflection statics of FGM plate.
(COV) of square of natural frequency ( 2 ) of square Al/ ZrO2 FGM plate Example:3 In this example comparative study have been performed
with Shaker et al. [24] has been performed. The volume fraction index between transverse deflection statics of Al FGM Al2 O3 sandwich
(n) and thickness ratio (a/ h ) in the study had been considered as 2 and FGM plate obtained with FOPT and MCS. Plate is assumed to be under
10, respectively. Plate is assumed to be simply supported at all edges. uniformly distributed loading and simply supported boundary condi-
The elastic modulus of metal and ceramic have been considered to be tions. Elastic modulus of ceramic(Ec ) has been considered to be random
70 GPa and 151 Gpa, respectively. Shaker et al. [24] employed first and variable. The volume fraction index(n) and thickness ratio(a/ h ) to be 2
second order reliability method whereas in present study FOPT has and 10 respectively. Fig. 8 shows that results obtained with FOPT
been used to obtain the second order statics of vibration response of the agrees with MCS and thus current formulation can be used to produce
FGM plates. Fig. 5 shows that results obtained from the present for- results for sandwich FGM plate.
mulation agrees well while considering elastic modulus (Ec and Em ) as
random variables (R. V .). 5.2. Parametric studies
Example:3 In this example coefficient of variation (COV) of square
of natural frequency ( 2 ) of simply supported square Objective of the parametric studies is to represent the effect of
Al FGM Al2 O3 sandwich FGM plate is validated with in- dispersion of various material properties and random variables on the
dependently developed Monte Carlo simulation (MCS) and FOPT ap- COV of square of natural frequency( 2 ) and transverse deflection(w) on
proach. Elastic modulus of ceramic (Ec ) has been considered as random Al FGM Al2 O3 skew sandwich plate. The material properties in the
variable (R.V.). Thickness ratio (a/ h ) and volume fraction index (n)of study have been assumed to be varying from 0% to 15%. The thickness
simply supported plates has been considered as 10 and 2, respectively. ratio of each layer is consider to be in 1:1:2 ratio moving from top to
The elastic modulus of metal and ceramic have been taken to be 70 GPa bottom layer. The material properties used in the analysis are as fol-
and 380 Gpa, respectively. The results obtained MCS and FOPT matches lows:
well as shown in Fig. 6. Since, the MCS is computationally expensive,
therefore FOPT approach has been used to obtain the results for sto- • Aluminium(Al): Elastic modulus (E ) = 70 GPa, Mass density
m
chastic vibration response of sandwich FGM plates. ( ) = 2707 Kg / m , Poisson's ratio ( ) = 0.3
3

• Aluminium Oxide ( Al O ): Elastic modulus (E ) = 380 GPa, Mass

m m
2 3 c

5.1.2. Bending analysis density ( c ) = 3800 Kg / m3 , Poisson's ratio ( c ) = 0.3

Example:1 In this example the non-dimensional central deflection
of Al FGM Al2 O3 sandwich FGM plate has been validated with 5.2.1. Vibration analysis
Neves et al. [69]. Validation has been performed while considering Example:1 This example presents the variation in COV of square of
different values of thickness ratios(a/ b ) whereas, volume fraction index natural frequency ( 2 ) of Al FGM Al2 O3 skew sandwich plate with
(n) has been considered to be 4. Neves et al. [69] used quasi 3D HSDT thickness ratio (a/ h ) by considering elastic modulus (Ec and Em ) as
while incorporating collocation method having radial basis function. random variables. Plate has been considered to be simply supported at
Table 3 shows comparison of results from both the methodologies. It all edges. The volume fraction index (n) and skew angle (ψ) of square
can be seen that at 6 × 6 mesh the results obtained with current for- sandwich plate is taken to be 2 and 100 , respectively. It is observed from
mulation agrees well with the published literature [69]. Fig. 9 that the values of COV of square of natural frequency ( 2 )

788
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 13. Variation of COV of transverse displacement with thickness ratio (a/ h ) considering (a) Ec (b) Em (c) c (d) m (e) n (f) Ec ,Em as random variable.

increases with increase in thickness ratio. This is more predominant in Em ) and poisson's ratio's ( c and m ) have been demonstrated. The
case of Ec as compared to Em . This is due to the fact that the value of thickness ratio (a/ h ) and skew angle (ψ) of skew sandwich plat are
elastic modulus of ceramic is higher than that of metal. considered to be 10 and 100 , respectively. Fig. 10 shows that COV of 2
Example:2 In this example the variation of COV of square of natural varies in cases where Em and m are varying. This is mainly due to the
frequency ( 2 ) of simply supported Al FGM Al2 O3 skew sandwich fact that on increasing the volume fraction index (n) of the plate, metal
plate with volume fraction index (n) considering elastic moduli (Ec and fraction in FGM layer increases and hence the variation in the COV can

789
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 14. Variation of COV of transverse displacement with thickness ratio (a/ h ) considering (a) Ec (b) Em (c) c (d) m as random variable.

be observed. Fig. 12a that with increase in number of independent random variables
Example:3 In this example various results for COV of square of (R.V.) COV of natural frequency ( 2 ) increases. This shows that the
natural frequency ( 2 ) of Al FGM Al2 O3 skew sandwich FGM effect of randomness in material properties have a significant effect on
plates with skew angle (ψ) has been demonstrated while varying all the sensitivity of natural frequency of plate. It is also observed from
material properties (Ec , Em , c , m , n) one by one. The thickness ratio Fig. 12b that the value of COV of natural frequency ( 2 ) is comes out to
(a/ h ) and volume fraction index (n) of the simply supported plates have be maximum in case of SFSF boundary conditions. This can be under-
been considered to be 2 and 10, respectively. It can be observed from stood as the this condition offers least constrained condition among
Fig. 11 that in all the cases the value of COV of 2 increases with in- various boundary conditions used, where as the value is minimum in
crease of skew angle (ψ). This is due to the increase in mean natural case of CCCC boundary conditions as it offers highly constrained con-
frequency of the system with skew angle. The numeric values of COV, dition for the plate.
2 is more in case of elastic moduli (E and E ) is more than poissons
c m
ratios ( c and m ) as the numeric mean values of elastic moduli are more 5.2.2. Bending analysis
than poissons ratios. Due to which small deviation in elastic modulus In this section various results have been shown while considering
largely influences the dispersion behaviour of natural frequency of square Al FGM Al2 O3 skew sandwich plates under simply sup-
system. ported boundary condition and with uniform loading conditions of
Example:4 In this examples, two special cases have been demon- 100 MPa.
strated. Firstly by considering the effect of increase of random material Example:1 This examples presents the dispersion of transverse
properties one by one on COV of square of natural frequency ( 2 ) later displacement (w) of square Al FGM Al2 O3 skew sandwich plate
by considering various boundary conditions. It can be observed from with thickness ratio (a/ h ) while considering all properties (Ec , Em , c ,

790
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

Fig. 15. Variation of COV of transverse displacement with volume fraction index (n) considering (a) Ec (b) Em (c) c (d) m as random variable.

Fig. 16. Variation of COV of square of natural frequency considering all properties as R.V.

m ,n) as random variables shown in Fig. 13. Plates are assumed to be in case of Em and m on the contrary COV decreases with increase in
under simply supported boundary condition. The volume fraction index volume fraction index (n) while considering Ec and c as random vari-
(n) and skew angle (ψ) are considered as 2 and 100 , respectively. It has ables. This is mainly due to the fact that with the increase in volume
been observed that COV (S. D. / Mean ) of transverse displacement (w) fraction index (n) metal volume fraction in FGM layer of sandwich plate
decreases with the increase in thickness ratio (a/ h ) while considering Ec increases. Due to which the sensitivity of transverse deflection (w) can
and c as random variable. While contrary COV of transverse dis- be easily observed.
placement (w) decreases with increase in Em , n and m . In general it can Example:3 This example presents the effect of dispersion of volume
be observed that COV of transverse displacement (w) is most sensitive fraction index(n) on COV of transverse displacement(w) of
to ceramic elastic modulus (Ec ) and metal poisson's ratio ( m ). Al FGM Al2 O3 skew sandwich plate considering more than one
Example:2 In this example the effect of variation of volume fraction material properties as random variable. The Thickness ratio and skew
index (n) on COV (S. D. / Mean ) of transverse displacement (w) have angle of the plate have been considering to be 10 and 100 , respectively.
been observed considering Ec , Em , c , m as independent random vari- Plate is assumed to be under simply supported boundary conditions.
ables. The thickness ratio (a/ h ) and skew angle (ψ) of the plate have From Fig. 15 it can be observed that with increase in volume fraction
been considered to be 10 and 100 , respectively. The variation has been index(n) COV of transverse displacement(w) increases with an excep-
observed at n = 1, 2, 5 and 10, respectively. It has been observed from tion at n = 1. It can be understood as elastic modulus of ceramic(Ec )
Fig. 14 that dispersion of COV of transverse displacement (w) increases and metal(Em ) being considered to be random. With results in increase

791
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

in COV of transverse displacement(w) due to (Em ) and decrease sue to References

(Ec ). Hence due to combined effect of both random variable such in-
teresting result has been obtained. [1] Zidi M, Tounsi A, Houari MSA, Adda Bedia EA, Anwar Bég O. Bending analysis of
Example:4 This example shows the effect of increase of random FGM plates under hygro-thermo-mechanical loading using a four variable refined
plate theory. Aero Sci Technol 2014;34(1):24–34. https://doi.org/10.1016/j.ast.
variable on the dispersion of transverse displacement (w) of 2014.02.001.
Al FGM Al2 O3 skew sandwich plate. Volume fraction index (n), [2] Koizumi M. FGM activities in Japan. Compos B Eng 1997;28(1–2):1–4. https://doi.
thickness ratio (a/ h ) and skew angle (ψ) of plate have been considered org/10.1016/S1359-8368(96)00016-9.
[3] Gupta A, Talha M. Recent development in modeling and analysis of functionally
to be 2, 10 and 100 , respectively. It has been observed from Fig. 16 that graded materials and structures. Prog Aero Sci 2015;79:1–14. https://doi.org/10.
COV of transverse displacement (w) increases with increase in number 1016/j.paerosci.2015.07.001.
of independent random variable. This shows that all the material [4] Nakagiri S, Takabatake H, Tani S. Uncertain eigenvalue analysis of composite la-
minated plates by the stochastic finite element method. J. Eng. Ind. 1987;109(1):9.
properties leads to the increase in the sensitivity of the transverse dis- https://doi.org/10.1115/1.3187096.
placement (w). [5] Oh D, Librescu L. Free vibration and reliability of composite cantilevers featuring
This example shows the effect of increase of random variable on the uncertain properties. Reliab Eng Syst Saf 1997;56(3):265–72. https://doi.org/10.
1016/S0951-8320(96)00038-5.
dispersion of transverse displacement(w) of Al FGM Al2 O3 skew
[6] Venini P, Mariani C. Computational structures technology free vibrations of un-
sandwich plate. Volume fraction index(n), thickness ratio(a/ h ) and certain composite plates via stochastic Rayleigh-Ritz approach. Comput Struct
skew angle(ψ) of plate have been considered to be 2,10 and 100 , re- 1997;64(1):407–23https://doi.org/10.1016/S0045-7949(96)00161-7.
spectively. It can be observed from Fig. 16 that COV of transverse [7] Navaneetha Raj B, Iyengar N, Yadav D. Response of composite plates with random
material properties using FEM and Monte Carlo simulation. Adv Compos Mater
displacement(w) increases with increase in number of independent 1998;7(3):219–37. https://doi.org/10.1163/156855198X00165.
random variable. This shows that all the material properties leads to the [8] Manohar CS, Ibrahim RA. Progress in structural dynamics with stochastic parameter
increase in the sensitivity of the transverse displacement(w). variations: 1987-1998. Appl Mech Rev 1999;52(5):177. https://doi.org/10.1115/1.
3098933.
[9] Singh BN, Yadav D, Iyengar NGR. Natural frequencies of composite plates with
6. Conclusion random material properties using higher-order shear deformation theory. Int J
Mech Sci 2001;43(10):2193–214. https://doi.org/10.1016/S0020-7403(01)
00046-7.
Stochastic vibration and bending response of functionally graded [10] Singh BN, Yadav D, Iyengar NGR. Free vibration of composite cylindrical panels
skew sandwich plates have been studied while considering the un- with random material properties. Compos Struct 2002;58(4):435–42. https://doi.
certainty in material properties. Reddy's HSDT has been used to model org/10.1016/S0263-8223(02)00133-2.
[11] Onkar AK, Yadav D. Non-linear response statistics of composite laminates with
the displacement field. Variational principle has been used to derive random material properties under random loading. Compos Struct
governing differential equations. Finite element methodology in con- 2003;60(4):375–83. https://doi.org/10.1016/S0263-8223(03)00049-7.
junction with FOPT have been used to calculate the first and second [12] Onkar AK, Yadav D. Non-linear free vibration of laminated composite plate with
random material properties. J Sound Vib 2004;272(3–5):627–41. https://doi.org/
order statics of natural frequency and transverse displacement of skew
10.1016/S0022-460X(03)00387-0.
sandwich plate. Convergence and validation studies have been per- [13] Onkar AK, Yadav D. Forced nonlinear vibration of laminated composite plates with
formed to show the reliability of present formulation. Various numer- random material properties. Compos Struct 2005;70(3):334–42. https://doi.org/10.
ical results have been shown to observe the effect of various random 1016/j.compstruct.2004.08.037.
[14] Yang J, Liew KM, Kitipornchai S. Stochastic analysis of compositionally graded
variables and system parameters on the vibration and deflection statics plates with system randomness under static loading. Int J Mech Sci
of functionally graded skew sandwich plates. The list of conclusions 2005;47(10):1519–41. https://doi.org/10.1016/j.ijmecsci.2005.06.006.
that can be drawn from current study are as follows: [15] Kitipornchai S, Yang J, Liew KM. Random vibration of the functionally graded la-
minates in thermal environments. Comput Methods Appl Mech Eng

• COV of natural frequency( ) and transverse displacement(w) varies

2006;195(9–12):1075–95. https://doi.org/10.1016/j.cma.2005.01.016.
2
[16] Giunta G, Carrera E. Stochastic static analyses of fe models by means of newton's
linearly with COV of random material properties. series expansions. III European Conference on Computational Mechanics.

• Effect of increase of thickness ratio(a/h) on COV of natural fre- Dordrecht: Springer Netherlands; 2006. p. 293. ISBN 978-1-4020-5370-2;.
[17] Tripathi V, Singh B, Shukla K. Free vibration of laminated composite conical shells
quency( 2 ) has been comes out to be more while considering Ec as with random material properties. Compos Struct 2007;81(1):96–104. https://doi.
random variable compared to Em . org/10.1016/j.compstruct.2006.08.002.

• Effect of increase of volume fraction index(n) on COV of natural [18] Hosseini S, Khadem S. Vibration and reliability of a rotating beam with random
properties under random excitation. Int J Mech Sci 2007;49(12):1377–88. https://
frequency( 2 ) have been comes out to be maximum while con- doi.org/10.1016/j.ijmecsci.2007.04.008.
sidering Em and m as random variable whereas effect is minimum [19] Gao W. Natural frequency and mode shape analysis of structures with uncertainty.
Mech Syst Signal Process 2007;21(1):24–39. https://doi.org/10.1016/j.ymssp.
for Ec and c , respectively.
•
2006.05.007.
The dispersion of elastic modulus of ceramic(Ec ) have maximum [20] Lal A, Singh BN, Kumar R. Static response of laminated composite plates resting on
effect on the sensitivity of natural frequency of skew sandwich FGM elastic foundation with uncertain system properties. J Reinforc Plast Compos
plate on contrary it found to be minimum in case of m . 2007;26(8):807–29. https://doi.org/10.1177/0731684407076730.

•
[21] Lal A, Singh BN, Kumar R. Natural frequency of laminated composite plate resting
COV of natural frequency( 2 ) increases with increase in number of on an elastic foundation with uncertain system properties. Struct Eng Mech
independent random variables considering other parameter to be 2007;27(2):199–222. https://doi.org/10.12989/sem.2007.27.2.199.
constant. [22] Lal A, Singh BN, Kumar R. Nonlinear free vibration of laminated composite plates

•
on elastic foundation with random system properties. Int J Mech Sci
As plate thickness ratio(a/h) the sensitivity of transverse deflection 2008;50(7):1203–12. https://doi.org/10.1016/j.ijmecsci.2008.04.002.
(w) increases while considering Em and m as R.V. whereas, on [23] Lal A, Singh BN, Kumar R. Stochastic nonlinear bending response of laminated
contrary sensitivity decreases while considering Em , m and n as composite plates with system randomness under lateral pressure and thermal
loading. Arch Appl Mech 2011;81(6):727–43. https://doi.org/10.1007/s00419-
random variable.
•
010-0442-7.
Sensitivity of transverse deflection(w) increases while considering [24] Shaker A, Abdelrahman W, Tawfik M, Sadek E. Stochastic Finite element analysis of
Em and m as random variable whereas, on contrary it decreases the free vibration of functionally graded material plates. Comput Mech
2008;41(5):707–14. https://doi.org/10.1007/s00466-007-0226-2.
while considering Ec as random variable.
•
[25] Singh BN, Lal A, Kumar R. Nonlinear bending response of laminated composite
COV of transverse deflection(w) increases with the increase in plates on nonlinear elastic foundation with uncertain system properties. Eng Struct
number of independent random variables. 2008;30(4):1101–12. https://doi.org/10.1016/j.engstruct.2007.07.007.
[26] Pandit M, Singh BN, Sheikh A. Stochastic perturbation-based finite element for
deflection statistics of soft core sandwich plate with random material properties. Int
Appendix A. Supplementary data J Mech Sci 2009;51(5):363–71. https://doi.org/10.1016/j.ijmecsci.2009.03.003.
[27] Shahabian F, Hosseini SM. Stochastic dynamic analysis of a functionally graded
Supplementary data to this article can be found online at https:// thick hollow cylinder with uncertain material properties subjected to shock loading.
Mater Des 2010;31(2):894–901. https://doi.org/10.1016/j.matdes.2009.07.040.
doi.org/10.1016/j.compositesb.2019.01.035.

792
S.S. Tomar, M. Talha Composites Part B 163 (2019) 779–793

[28] Hosseini SM, Shahabian F. Reliability of stress field in AlAl2O3 functionally graded [49] Dey S, Mukhopadhyay T, Adhikari S. Stochastic free vibration analysis of angle-ply
thick hollow cylinder subjected to sudden unloading, considering uncertain me- composite plates - A RS-HDMR approach. Compos Struct 2015;122:526–36. https://
chanical properties. Mater Des 2010;31(8):3748–60. https://doi.org/10.1016/j. doi.org/10.1016/j.compstruct.2014.09.057.
matdes.2010.03.016. [50] Dey S, Mukhopadhyay T, Sahu S, Li G, Rabitz H, Adhikari S. Thermal uncertainty
[29] Chandrashekhar M, Ganguli R. Nonlinear vibration analysis of composite laminated quantification in frequency responses of laminated composite plates. Compos B Eng
and sandwich plates with random material properties. Int J Mech Sci 2015;80:186–97. https://doi.org/10.1016/j.compositesb.2015.06.006.
2010;52(7):874–91. https://doi.org/10.1016/j.ijmecsci.2010.03.002. [51] Dey S, Mukhopadhyay T, Khodaparast HH, Adhikari S. Fuzzy uncertainty propa-
[30] Pandit MK, Singh BN, Sheikh AH. Stochastic free vibration response of soft core gation in composites using Gram-Schmidt polynomial chaos expansion. Appl Math
sandwich plates using an improved higher-order zigzag theory. J Aero Eng Model 2016;40(7–8):4412–28. https://doi.org/10.1016/j.apm.2015.11.038.
2010;23(1):14–23. https://doi.org/10.1061/(ASCE)0893-1321(2010)23:1(14). [52] Dey S, Naskar S, Mukhopadhyay T, Gohs U, Spickenheuer A, Bittrich L, et al.
[31] Pandit MK, Singh BN, Sheikh AH. Vibration of sandwich plates with random ma- Uncertain natural frequency analysis of composite plates including effect of noise A
terial properties using improved higher-order zig-zag theory. Mech Adv Mater polynomial neural network approach. Compos Struct 2016;143:130–42. https://
Struct 2010;17(7):561–72. https://doi.org/10.1080/15376490903398698. doi.org/10.1016/j.compstruct.2016.02.007.
[32] Lal A, Singh BN. Effect of random system properties on bending response of thermo- [53] Dey S, Mukhopadhyay T, Adhikari S. Metamodel based high-fidelity stochastic
mechanically loaded laminated composite plates. Appl Math Model analysis of composite laminates: A concise review with critical comparative as-
2011;35(12):5618–35. https://doi.org/10.1016/j.apm.2011.05.014. sessment. 2017. https://doi.org/10.1016/j.compstruct.2017.01.061.
[33] Lal A, Singh BN, Anand S. Nonlinear bending response of laminated composite [54] Corradi M, Borri A, Righetti L, Speranzini E. Uncertainty analysis of FRP reinforced
spherical shell panel with system randomness subjected to hygro-thermo-mechan- timber beams. Compos B Eng 2017;113:174–84. https://doi.org/10.1016/j.
ical loading. Int J Mech Sci 2011;53(10):855–66. https://doi.org/10.1016/j. compositesb.2017.01.030.
ijmecsci.2011.07.008. [55] Tomar SS, Zafar S, Talha M, Gao W, Hui D. State of the art of composite structures
[34] Giunta G, Carrera E, Belouettar S. Free vibration analysis of composite plates via in non-deterministic framework: A review. Thin-Walled Struct 2018;132:700–16.
refined theories accounting for uncertainties. Shock Vib 2011;18(4):537–54. [56] Wu D, Gao W, Hui D, Gao K, Li K. Stochastic static analysis of Euler-Bernoulli type
[35] Jagtap KR, Lal A, Singh BN. Stochastic nonlinear bending response of functionally functionally graded structures. Compos B Eng 2018;134:69–80. https://doi.org/10.
graded material plate with random system properties in thermal environment. Int J 1016/j.compositesb.2017.09.050.
Mech Mater Des 2012;8(2):149–67. https://doi.org/10.1007/s10999-012-9183-9. [57] Karsh P, Mukhopadhyay T, Dey S. Stochastic dynamic analysis of twisted func-
[36] Lal A, Choksi P, Singh BN. Stochastic nonlinear free vibration analysis of piezo- tionally graded plates. Compos B Eng 2018;147:259–78. https://doi.org/10.1016/
laminated composite conical shell panel subjected to thermoelectromechanical J.COMPOSITESB.2018.03.043.
loading with random material properties. J Appl Mech 2012;79(6):061008https:// [58] Karsh P, Mukhopadhyay T, Dey S. Stochastic low-velocity impact on functionally
doi.org/10.1115/1.4006765. graded plates: Probabilistic and non-probabilistic uncertainty quantification.
[37] Sepahvand K, Marburg S, Hardtke HJ. Stochastic free vibration of orthotropic plates Compos B Eng 2019;159:461–80. https://doi.org/10.1016/J.COMPOSITESB.2018.
using generalized polynomial chaos expansion. J Sound Vib 2012;331(1):167–79. 09.066.
https://doi.org/10.1016/j.jsv.2011.08.012. [59] Tomar SS, Talha M. Thermo-mechanical buckling analysis of functionally graded
[38] Shegokar NL, Lal A. Stochastic nonlinear bending response of piezoelectric func- skew laminated plates with initial geometric imperfections. Int. J. Appl. Mech.
tionally graded beam subjected to thermoelectromechanical loadings with random 2018;0(0):1850014. https://doi.org/10.1142/S175882511850014X.
material properties. Compos Struct 2013;100:17–33. https://doi.org/10.1016/j. [60] Tomar SS, Talha M. Large amplitude vibration analysis of functionally graded la-
compstruct.2012.12.032. minated skew plates in thermal environment. Mech Adv Mater Struct
[39] Chang TP. Stochastic FEM on nonlinear vibration of fluid-loaded double-walled 2017;0(0):1–14. https://doi.org/10.1080/15376494.2017.1400619.
carbon nanotubes subjected to a moving load based on nonlocal elasticity theory. [61] Tomar SS, Talha M. On the flexural and vibration response of imperfection sensitive
Compos B Eng 2013;54:391–9. https://doi.org/10.1016/j.compositesb.2013.06. higher order skewed fgm sandwich plates in thermal environment. Proc IME C J
012. Mech Eng Sci 2018;0(0). https://doi.org/10.1177/0954406218766959.
[40] Kumar R, Patil HS, Lal A. Nonlinear flexural response of laminated composite plates [62] Shankara C, Iyengar N. A c0 element for the free vibration analysis of laminated
on a nonlinear elastic foundation with uncertain system properties under lateral composite plates. J Sound Vib 1996;191(5):721–38. https://doi.org/10.1006/jsvi.
pressure and hygrothermal loading: micromechanical model. J Aero Eng 1996.0152.
2014;27(3):529–47. https://doi.org/10.1061/(ASCE)AS.1943-5525.0000219. [63] Gupta A, Talha M. An assessment of a non-polynomial based higher order shear and
[41] Dash P, Singh BN. Static response of geometrically nonlinear laminated composite normal deformation theory for vibration response of gradient plates with initial
plates having uncertain material properties. Mech Adv Mater Struct geometric imperfections. Compos B Eng 2016;107:141–61. https://doi.org/10.
2015;22(4):269–80. https://doi.org/10.1080/15376494.2012.736056. 1016/j.compositesb.2016.09.071.
[42] Talha M, Singh BN. Stochastic vibration characteristics of finite element modelled [64] Gupta A, Talha M. Nonlinear flexural and vibration response of geometrically im-
functionally gradient plates. Compos Struct 2015;130:95–106. https://doi.org/10. perfect gradient plates using hyperbolic higher-order shear and normal deformation
1016/j.compstruct.2015.04.030. theory. Compos B Eng 2017;123:241–61. https://doi.org/10.1016/j.compositesb.
[43] Xu Y, Qian Y, Chen J, Song G. Stochastic dynamic characteristics of FGM beams 2017.05.010.
with random material properties. Compos Struct 2015;133:585–94. https://doi. [65] Gupta A, Talha M, Singh BN. Vibration characteristics of functionally graded ma-
org/10.1016/j.compstruct.2015.07.057. terial plate with various boundary constraints using higher order shear deformation
[44] Asnafi A, Abedi M. A complete analogical study on the dynamic stability analysis of theory. Compos B Eng 2016;94:64–74. https://doi.org/10.1016/j.compositesb.
isotropic functionally graded plates subjected to lateral stochastic loads. Acta Mech 2016.03.006.
2015;226(7):2347–63. https://doi.org/10.1007/s00707-015-1321-7. [66] Garg AK, Khare RK, Kant T. Free vibration of skew fiber-reinforced composite and
[45] Chakraborty S, Mandal B, Chowdhury R, Chakrabarti A. Stochastic free vibration sandwich laminates using a shear deformable finite element model. J Sandw Struct
analysis of laminated composite plates using polynomial correlated function ex- Mater 2006;8(1):33–53. https://doi.org/10.1177/1099636206056457.
pansion. Compos Struct 2016;135:236–49. https://doi.org/10.1016/j.compstruct. [67] Kleiber M, Hien T. The stochastic finite element method (basic perturbation tech-
2015.09.044. nique and computer implementation). Appl Stoch Model Data Anal 1994;10(4).
[46] García-Macías E, Castro-Triguero R, Friswell MI, Adhikari S, Sáez A. Metamodel- https://doi.org/10.1002/asm.3150100412. 297–297.
based approach for stochastic free vibration analysis of functionally graded carbon [68] El Meiche N, Tounsi A, Ziane N, Mechab I, AddaBedia EA. A new hyperbolic shear
nanotube reinforced plates. Compos Struct 2016;152:183–98. https://doi.org/10. deformation theory for buckling and vibration of functionally graded sandwich
1016/j.compstruct.2016.05.019. plate. Int J Mech Sci 2011;53(4):237–47. https://doi.org/10.1016/j.ijmecsci.2011.
[47] Mukhopadhyay T, Naskar S, Dey S, Adhikari S. On quantifying the effect of noise in 01.004.
surrogate based stochastic free vibration analysis of laminated composite shallow [69] Neves A, Ferreira A, Carrera E, Cinefra M, Roque C, Jorge R, et al. Static, free
shells. Compos Struct 2016;140:798–805. https://doi.org/10.1016/j.compstruct. vibration and buckling analysis of isotropic and sandwich functionally graded
2015.12.037. plates using a quasi-3D higher-order shear deformation theory and a meshless
[48] Dey S, Mukhopadhyay T, Adhikari S. Stochastic free vibration analyses of composite technique. Compos B Eng 2013;44(1):657–74. https://doi.org/10.1016/j.
shallow doubly curved shells A Kriging model approach. Compos B Eng compositesb.2012.01.089.
2015;70:99–112. https://doi.org/10.1016/j.compositesb.2014.10.043.

793