See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/358933027

FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations

for Membrane Effects

Article  in  Teknik Dergi · March 2022

DOI: 10.18400/tekderg.878982

CITATION READS

1 8

1 author:

Ülkü H. Çalık-Karaköse
Istanbul Technical University
7 PUBLICATIONS   22 CITATIONS   

SEE PROFILE

All content following this page was uploaded by Ülkü H. Çalık-Karaköse on 31 October 2022.

The user has requested enhancement of the downloaded file.

Teknik Dergi, 2022 11799-11822, Paper 659

FE Analysis of FGM Plates on Arbitrarily Orthotropic

Pasternak Foundations for Membrane Effects*

Ülkü Hülya ÇALIK-KARAKÖSE1

ABSTRACT
In this study, the finite element analysis of sigmoid functionally graded material (S-FGM)
plates resting on orthotropic Pasternak elastic foundations with different material angles is
presented. For modelling, SAP2000 software package is used in which the required
adjustments are made to obtain the expected behaviour of the plate-foundation system. The
plate is modelled using both solid elements and layered shell elements by defining a
number of solid elements and layers in the thickness direction having elastic properties
equivalent to the properties of the S-FGM plate. The interaction between the plate and the
foundation is provided by equalizing the vertical displacements of the plate and foundation
nodal points. The orthotropic Pasternak foundation is modelled using plane strain elements
with some adjustments to the elastic properties. The membrane effects of the simply
supported S-FGM plate on Pasternak foundation are considered by defining the edge
boundaries of the system as pinned supports. These effects are excluded by converting all
boundary nodes into roller supports except for one of the corner nodes of the plate and the
foundation due to the stability requirement. A number of verification examples are
performed to demonstrate the convenience and robustness of the proposed model. This
work can be easily extended to static and dynamic analyses of FGM plates with various
geometries resting on arbitrarily orthotropic Pasternak elastic foundations for further
studies.
Keywords: Functionally graded material plates, membrane effects, orthotropic Pasternak
foundation, material angle.

1. INTRODUCTION
Functionally graded materials (FGMs) are heterogeneous materials in which the material
properties are varied continuously in the thickness direction as presented by Suresh and
Mortenson [1], and are used to improve the strength and stiffness of structural elements in
many engineering applications. In laminated composites which are used to satisfy the
desired high-performance demands, stress singularities may occur at the interface of two

Note:
- This paper was received on February 12, 2021 and accepted for publication by the Editorial Board
on October 22, 2021.
- Discussions on this paper will be accepted by May 31, 2022.
 https://doi.org/10.18400/tekderg.878982

1 Department of Civil Engineering, Istanbul Technical University, Istanbul, Turkey

[email protected] - https://orcid.org/0000-0002-2944-7434
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

different materials which may cause de-bonding or cracking. This problem is eliminated by
the gradual change of the material properties of the functionally graded materials. Studies
on FGMs were first made for thermal barriers by Koizumi [2] and thermal mechanical
behaviour of FGMs was investigated in [3,4].
Reddy [5] presented theoretical formulation and finite element models of functionally
graded rectangular plates based on the third-order shear deformation theory. Cheng and
Batra [6] derived field equations for a functionally graded plate using both first-order and
third-order shear deformation theories.
Three different homogenization methods are used to describe the variation of material
properties; namely power law [7,8], exponential [9,10] and sigmoid functions. Since the
material properties vary continuously but rapidly and stress concentrations appear in one of
the interfaces in both the power-law and exponential functions, Chung and Chi [11] defined
a new volume fraction proposing an S-FGM, which was composed of two power-law
functions. It is demonstrated that the stress intensity factors of a cracked body are
significantly reduced by Chi and Chung [12] through the use of an S-FGM,. Chi and Chung
[13] derived theoretical formulations of a simply supported rectangular FGM plate of
medium thickness subjected to transverse loading and numerical solutions were obtained by
using the formulations derived by Chi and Chung [13], and also checked by the finite
element method in Chi and Chung [14]. The coupling effect of extension and bending in
functionally graded plates under transverse loading was investigated by Orakdöğen
et.al.[15] using the finite element method. Zenkour [16] introduced the bending analysis of
a simply supported functionally graded rectangular plate subjected to a transverse uniform
load. Elishakoff et.al.[17] presented the static response of a three-dimensional functionally
graded material rectangular plate with clamped boundary conditions and a power-law
distribution of the mechanical properties is adopted in the modelling. Swaminathan
et.al.[18] presented an extensive review of the various methods used in studying the static,
dynamic and stability behaviours of FGM plates where analytical and numerical methods
are both considered.
Plates resting on elastic foundations have been widely studied by many researchers to be
used in modelling various engineering applications. The simplest elastic foundation model
is the Winkler or one-parameter model by Winkler [19] in which the interaction between
the plate and the foundation is represented by a series of independent vertical linear elastic
springs. A more realistic representation of the elastic medium is made with two-parameter
or Pasternak model by Pasternak [20] which has a second dependent parameter that
introduces the shear interaction between springs avoiding the deflection discontinuity on
the interacted surface of the plate.
There exist many studies performing the static analysis of isotropic rectangular plates
resting on isotropic Pasternak foundations [21,22] and vibration analysis of this type of
plate-foundation system [23-26]. Shen [27] carried out the post-buckling analysis of a
simply supported, composite laminated rectangular plate resting on an isotropic two-
parameter elastic foundation where the plate is subjected to in-plane loading. Huang et.al.
[28] presented exact solutions for functionally graded thick plates based on the three-
dimensional theory of elasticity and examined interactions between the Winkler-Pasternak
elastic foundation and the plate demonstrating that elastic foundations have significant
effect on the mechanical behaviour of functionally graded thick plates. Lee et.al. [29]

11800
 Ülkü Hülya ÇALIK-KARAKÖSE

developed a refined higher order shear and normal deformation theory for the bending
analysis of power law, exponential and S-FGM plates resting on isotropic Pasternak
foundations. Vibration analysis of isotropic FGM plates on isotropic two-parameter
foundations is also performed by many researchers, [30-33]. Mansouri and Shariyat [34]
performed the analysis of orthotropic FGM plates on elastic foundations in their paper. In
[35,36], seismic behaviour of a 3D structural system resting on a two-parameter elastic
foundation was analysed using SAP2000. In the modelling of the foundation, the part of the
governing equation that belongs to Winkler behaviour is represented by elastic springs and
the part that belongs to elastic shear layer is taken into account by using four-noded
quadrilateral finite elements with four degrees of freedom.
There are also studies in the literature on the analysis of homogeneous plates and FGM
plates resting on orthotropic Pasternak foundations [37-42]. In this paper, the finite element
analysis of S-FGM plates resting on orthotropic Pasternak elastic foundations with different
material angles is presented and coupling effect of extension and bending is investigated.
SAP2000 software package, [43] is used in the implementation process and some
adjustments are made to the software for the expected behaviour of the plate-foundation
system. S-FGM plate is first modelled using layered shell elements by defining a number of
layers in the thickness direction having elastic properties equivalent to the properties of the
S-FGM plate. The plate is also modelled using solid elements. The membrane force-
bending moment coupling effect of the simply supported S-FGM plate is considered by
assigning the edge boundary nodes as pinned supports.
The interaction between the plate and the foundation is provided by equalizing the plate and
corresponding foundation nodal point vertical displacements. Orthotropic Pasternak
foundation is modelled using plane strain elements with some adjustments to the elastic
properties. After verifying the proposed model with a number of examples taken from the
literature, the influence of the material properties, membrane effects, number of layers used
for FGM plate and material angles of the orthotropic Pasternak foundation on deflections
and stresses are examined with numerical examples.

2. MATERIAL PROPERTIES OF S-FGM PLATES

In Figure 1, an S-FGM plate example is shown where z axis is in the thickness direction
and originated at the middle surface of the plate. The material properties vary continuously
and functionally based on the volume fraction of the sigmoid distribution.

x
S-FGM plate
E(z), (z)
h/2
y h/2
z

Figure 1 - An S-FGM plate example

11801
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

The volume fraction for S-FGM is defined using two power-law functions to ensure smooth
distribution of stresses among all interfaces of the plate, Chung and Chi [11]. The two
power-law functions are defined by

𝑉 (𝑧) = 1 − for 0≤𝑧≤ (1)

𝑉 (𝑧) = for − ≤𝑧≤0 (2)

where p is the material exponent. The elastic modulus of the S-FGM can be calculated
using

𝐸(𝑧) = 𝑉 (𝑧)(𝐸 − 𝐸 ) + 𝐸 for 0≤𝑧≤ (3)

𝐸(𝑧) = 𝑉 (𝑧)(𝐸 − 𝐸 ) + 𝐸 for − ≤𝑧≤0 (4)

Here, E1 and E2 are the elastic moduli of the bottom and top surfaces of the plate,
respectively. The Poisson’s ratio is assumed to be constant since it is stated by Delale and
Erdogan [9] that its effect on deformation is much less than that of the elastic modulus.

3. FINITE ELEMENT MODELLING

3.1. S-FGM Plate Model
The S-FGM plate is first modelled using shell finite elements for which the section type is
defined as “Shell-Layered/Nonlinear” in SAP2000, consisting of a number of layers in the
thickness direction. Each layer has a constant elastic modulus value which is pre-calculated
using the functions given in Eqs. [1]-[4] representing the continuous variation of the elastic
modulus of the plate starting from the bottom surface. The distance parameter “z” in the
functions indicates the ordinate of the middle section of each layer. In Figure 2, a segment
of an S-FGM plate and its equivalent representation with 8 layers is shown.
The shell element in the software is a four-noded finite element and its formulation
combines membrane and plate-bending behaviours. The element always activates all six
degrees of freedom which are 3 translations and 3 rotations at each of its connected joints.
When the element is used as a pure plate, it must be ensured that restraints or other supports
are provided to the degrees of freedom for in-plane translations and the rotation about the
normal.
The layered shell allows any number of layers to be defined in the thickness direction, each
with an independent location, thickness, behaviour, and material. For bending, a Mindlin
formulation is used which always includes transverse shear deformations. Membrane and
plate behaviours will be coupled if the layering is not symmetrical in the thickness direction
and the edge boundaries are assigned as pinned supports.

11802
 Ülkü Hülya ÇALIK-KARAKÖSE

The S-FGM plate is also modelled using solid elements where the elastic modulus of each
solid element in the thickness direction is equal to the elastic modulus of the corresponding
layer of the shell element.
The solid element used is an eight-noded finite element with three translational degrees of
freedom at each joint. Rotational degrees of freedom are not activated. This element
contributes stiffness to all of these translational degrees of freedom.
The membrane effects appear in the simply supported S-FGM plate when the edge
boundary nodes of the plate are assumed to be pinned supports (no horizontal
displacement) whereas these effects are excluded by converting all boundary nodes into
roller supports except at one of the corner nodes due to the stability requirement.
Layer
8

7
h/2

5
x,y
4

z 3
h/2

Figure 2 - A segment of an S-FGM plate and equivalent representation with 8 layers

3.2. Isotropic and Orthotropic Pasternak Foundation Models

Isotropic and orthotropic Pasternak foundations are also modelled using SAP2000. The
material symmetry type of the isotropic soil element is selected as “orthotropic” so as to
change the shear modulus values. All material properties are assigned to zero except the
shear moduli G13 and G23 that represent the second parameter of the Pasternak foundation.
The finite element to be used in modelling the Pasternak foundation is defined as “shell”
element, the section type is selected as “Plane-Strain” and the thickness is assigned a unit
value, [44,45].
The plane strain element is one type of area objects and is used to model plane-strain
behaviour in two-dimensional solids. It activates the three translational degrees of freedom
at each of its joints. Rotational degrees of freedom are not activated. The plane-strain
element models shear that is normal to the plane of the element, in addition to the in-plane
behaviour. Thus, stiffness is created for all three translational degrees of freedom.
In the case of plane strain, when the elastic moduli and Poisson’s ratios in three directions
are defined as zero and the in-plane shear moduli in the two directions (G13 and G23) are
assigned to non-zero values, only shear stresses in the thickness direction of the element
and end forces in the vertical direction appear. When the in-plane shear moduli are assigned
to the shear moduli of the two-parameter foundation, the resulting element stiffness matrix
is converted into the element stiffness matrix of the two-parameter orthotropic soil element.
If the in-plane shear moduli of the soil element are assigned to the same values, the
stiffness matrix of the isotropic soil element is obtained. The unknowns are only the
deflections of the nodal points.

11803
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

The first parameter of the Pasternak foundation which is represented by springs is taken
into account as “area springs” assigned to the surface of the foundation.
For the orthotropic foundation model with different material angles, shear moduli G13 and
G23 are assigned to different values and a desired material angle is selected.

3.3. Interaction between Plate and Pasternak Foundation

The plate-foundation interaction is provided by matching the vertical displacements of the
plate and the corresponding foundation nodal points. “Equal constraint” property is used for
this purpose where the number of constraints is equal to the number of nodes of the plate.
Each constraint is then assigned to the respective plate and foundation nodes.
The S-FGM plate model of the system is created using both layered shell and solid
elements as shown in Figure 3.
The membrane effects are excluded by releasing the in-plane translational degrees of
freedom of both the plate and foundation edge nodes. Note that the in-plane translations of
one of the corner nodes of the plate and the foundation are restrained for the stability
requirement.

a) b)
Figure 3 - Plate-foundation finite element model with membrane effects
a) using layered shell finite elements for the plate b) using solid finite elements for the plate

4. VERIFICATION EXAMPLES
4.1. Simply supported S-FGM Plate
A simply supported square S-FGM plate which has been solved both analytically by Chi
and Chung [14] and numerically by [14,15] is analysed and the results are compared only
with the results of Orakdöğen et.al.[15] because membrane effect analysis was not carried
out by Chi and Chung [14]. The thickness of the plate is h =0.02 m, the length is a=1 m
and the Poisson’s ratio is υ= 0.3. The elastic modulus at the bottom surface is E1=25x106
kN/m2 and the material exponent is p=2. The plate is subjected to a uniformly distributed

11804
 Ülkü Hülya ÇALIK-KARAKÖSE

load of q=100 kN/m2. The problem is solved for two different elastic modulus ratios,
E1/E2= 30 and E1/E2=2 where E2 indicates the elastic modulus of the top surface.
First, the plate is modelled using 8 layered shell finite elements and (4x4), (8x8) and
(20x20) meshes for E1/E2= 30, including membrane effects. This mesh refinement shows
that convergence to the deflections along the x direction (at y=0.5 m) obtained by
Orakdöğen et.al.[15] is achieved using an (8x8) mesh, Figure 4. However, due to one-to-
one comparison of nodal values, the same mesh discretization with Orakdöğen et.al.[15]
which is a (20x20) mesh has been used in the analyses.
The same example is also solved using 2 and 4 layered shell elements to see if the results
change with the number of layers used in the modelling. It is observed that for the case
E1/E2=2, the dimensionless deflections are almost the same for all number of layers both for
the plate with and without membrane effects, Figure 5.
x/a
0,50 0,60 0,70 0,80 0,90 1,00
0,00
0,05
0,10
0,15
w/h

0,20
(4x4 mesh)
0,25 (8x8 mesh)
(20x20 mesh)
0,30
[15]
0,35

Figure 4 - Variation of dimensionless deflections with mesh refinement (membrane effects

are included)
x/a x/a
0,50 0,60 0,70 0,80 0,90 1,00 0,5 0,6 0,7 0,8 0,9 1
0,00 0,00
0,02 0,02
0,04 0,04
0,06 0,06
0,08 0,08
w/h

w/h

0,10 0,10
0,12 0,12
2 layers
0,14 2 layers 0,14
0,16 0,16 4 layers
4 layers
0,18 0,18 8 layers
8 layers 0,20
0,20

a) b)
Figure 5 - Comparison of dimensionless deflections for different number of layers
(E1/E2=2)
a) membrane effects are included b) membrane effects are excluded

11805
 FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

On the other hand, the dimensionless deflections of the plate with elastic modulus ratio
E1/E2=30 increase as the number of layers used increases and this difference becomes much
larger when the membrane effects are excluded as shown in Figure 6.
The dimensionless deflections of the S-FGM plate along the x direction (at y=0.5 m) for the
two elastic modulus ratios and for the cases with and without membrane effects (m.e.) are
obtained and compared with the solution given by Orakdöğen et.al.[15], Figure 7. Note that
8 layers are used in the modelling while the number of layers used is 10 in Orakdöğen
et.al.[15]. It is verified that the results obtained in this study have very good agreement with
the reference solution.

x/a x/a
0,50 0,60 0,70 0,80 0,90 1,00 0,50 0,60 0,70 0,80 0,90 1,00
0,00 0,00
0,05 0,05
0,10 0,10
0,15 0,15
0,20 0,20
w/h

0,25
w/h

0,25
0,30 2 layers
0,30 2 layers
0,35 0,35
4 layers 4 layers
0,40 0,40
8 layers 0,45 8 layers
0,45
0,50 0,50

a) b)
Figure 6- Comparison of dimensionless deflections for different number of layers
(E1/E2=30)
a) membrane effects are included b) membrane effects are excluded

x/a x/a
0,50 0,60 0,70 0,80 0,90 1,00 0,50 0,60 0,70 0,80 0,90 1,00
0,00 0,00
0,05 0,02
0,10 0,04
0,15 0,06
0,20 0,08
w/h
w/h

0,25 0,10
0,30 0,12
m.e. included 0,14 m.e. included
0,35
0,40 [15] m.e. included 0,16 [15] m.e. included

0,45 m.e. excluded 0,18 m.e. excluded

0,50 [15] m.e. excluded 0,20 [15] m.e. excluded

a) b)
Figure 7 - Comparison of dimensionless deflections a) E1/E2=30 b) E1/E2=2

11806
 Ülkü Hülya ÇALIK-KARAKÖSE

The example is modelled and solved using solid finite elements as well and the
dimensionless deflections are compared with the values obtained using the model with
layered shell elements. It is observed that solid elements can also be used for the modelling
since all results match, Figure 8. It is also obtained that the mid-point deflections increase
1.8 times and 2.5 times for the cases with and without membrane effects, respectively as the
elastic modulus decreases with a larger ratio from the bottom to the top surface of the plate.
x/a x/a
0,50 0,60 0,70 0,80 0,90 1,00 0,50 0,60 0,70 0,80 0,90 1,00
0,00 0,00
0,05 0,05
0,10 0,10
0,15 0,15
0,20 0,20
w/h

0,25 w/h 0,25

0,30 8 layered shell , E1/E2=2 0,30
0,35 8 layered solid , E1/E2=2 0,35 8 layered shell , E1/E2=2
0,40 8 layered shell , E1/E2=30 0,40 8 layered solid , E1/E2=2
0,45 8 layered solid , E1/E2=30 0,45 8 layered shell , E1/E2=30
0,50 0,50 8 layered solid , E1/E2=30

a) b)
Figure 8- Comparison of dimensionless deflections using shell and 3D solid elements
a) membrane effects are included b) membrane effects are excluded

-0,5
(4x4) mesh , m.e. included
(8x8) mesh , m.e. included
-0,4
(20x20) mesh , m.e. included
(4x4) mesh , m.e. excluded
-0,3
(8x8) mesh , m.e. excluded
(20x20) mesh , m.e. excluded
-0,2

-0,1
-800 200 1200
0,0
z/h

0,1

0,2

0,3

0,4

0,5 x/q

Figure 9 - Variation of dimensionless stresses with mesh refinement

11807
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

-0,5

-0,4
E1/E2=2
[15] E1/E2=2
-0,3
E1/E2=30
[15] E1/E2=30
-0,2

-0,1
-800 200 1200
0,0
z/h

0,1

0,2

0,3

0,4

0,5 x/q

Figure 10- Comparison of dimensionless stresses (membrane effects are excluded)

-0,5 -0,5

-0,4 8 layered shell , m.e. included -0,4 8 layered shell , m.e. included
8 layered solid , m.e. included
8 layered solid , m.e. included
-0,3 8 layered shell , m.e. excluded -0,3
8 layered shell , m.e. excluded
8 layered solid , m.e. excluded
-0,2 -0,2 8 layered solid , m.e. excluded

-0,1 -0,1
-800 200 1200 -800 200 1200
0,0 0,0
z/h

z/h

0,1 0,1

0,2 0,2

0,3 0,3

0,4 0,4

0,5 x/q 0,5 x/q

a) b)
Figure 11- Comparison of dimensionless stresses using layered shell and 3D solid finite
elements a) E1/E2=2 b) E1/E2=30

11808
 Ülkü Hülya ÇALIK-KARAKÖSE

A mesh refinement study using (4x4), (8x8) and (20x20) meshes is also carried out for the
dimensionlessx stress distributions through the thickness of the plate centre and results are
given comparatively in Figure 9. Note that 8 layered shell finite elements are used in the
models and the elastic modulus ratio is E1/E2=30. It is observed that all values match except
for (4x4) mesh where a negligible difference occurs when membrane effects are included.
The dimensionlessx stress distributions through the thickness of the plate centre are
obtained using a (20x20) mesh and 8 layered shell finite elements for different elastic
modulus ratios and compared with the values given in [15]. It is observed that the results
have very good agreement with the reference solution, as shown in Figure 10.
The dimensionless stress distributions are also verified with the results obtained using the
3D plate model for the cases with and without membrane effects. It is observed that the
dimensionless stress values decrease with the inclusion of membrane effects and this
decrement is more apparent when E1/E2 ratio increases. It is also seen that the neutral axis
moves towards the mid-section of the plate when membrane effects are included, as in
Figure 11.

4.2. Homogeneous Plate on Isotropic Pasternak Foundation

Static analysis of a rectangular plate resting on an isotropic Pasternak foundation subjected
to a uniform loading and a central point load is performed and compared with the reference
solutions given in [21,22] to verify the proposed plate-foundation model. The geometrical
and loading properties of the plate-foundation system are shown in Figure 12.
9.144

P=133.34 kN
12.192

plate q=23.94 kN/m2

9.144 m

foundation

9.144 m 9.144 9.144

Figure 12 - Plan and vertical section of plate-foundation system

The plate thickness is h=0.1524 m and the material properties of the plate are;
E=2.0685x106 kN/m2 and υ=0.2. The foundation is extended in the in-plane directions and
the plate rests freely on the foundation without assigning any boundary conditions to the
edges. Here, the plate is modelled using thin plate finite elements.
The four-noded plate finite element used for modelling homogeneous plates has 3 degrees
of freedom at each node, one translational degree of freedom in the direction normal to the
plane of the element and 2 out-of-plane rotational degrees of freedom. In the thick plate

11809
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

finite element formulation, the effects of transverse shear deformation are included whereas
in the thin plate finite element formulation, transverse shearing deformation is neglected.
The midpoint deflections of the plate under uniformly distributed loading and under central
point load for different soil depths are given in Table 1 and Table 2, respectively. Here, kw
is the first parameter (the coefficient of subgrade reaction), kp is the second parameter (the
shear modulus) of the Pasternak foundation.
The results are in good agreement with the reference solutions indicating that the proposed
model fully represents the expected behaviour of the isotropic plate-foundation system.

Table 1 - Midpoint deflections of the plate under uniformly distributed loading for different
soil depths
H(m) 3.048 6.096 9.144 15.24
3
kw (kN/m ) 27206.59 13757.63 9377.96 5964.21
2
kp (kN/m ) 26905.91 50410.87 70586.50 104664.46
Ref.[22] 0.08720 0.15260 0.18930 0.22120
Ref.[21] 0.08530 0.15240 0.18900 0.20700
w x10-2(m)
Present
0.08713 0.15300 0.18900 0.22100
study

Table 2 - Midpoint deflections of the plate under central point load for different soil depths
H(m) 3.048 6.096 9.144 15.24
3
kw (kN/m ) 31898.08 24256.05 23737.98 23710.59
2
kp (kN/m ) 18913.08 23597.82 24035.56 24060.66
Ref.[22] 0.08180 0.08450 0.08460 0.08460
Ref.[21] 0.04800 0.09750 0.09750 0.09750
w x10-2 (m)
Present
0.07617 0.08052 0.08076 0.08077
study

4.3. Simply Supported Homogeneous Plate on Isotropic Pasternak Foundation

Example 1:
A simply supported homogeneous square plate resting on an isotropic Pasternak foundation
given in [28,46] is modelled using thin plate finite elements and solved for different soil
parameters. The plate is subjected to a uniformly distributed load of q=E/105 kN/m2. Length
to height ratio of the plate is a/h =100 and the Poisson’s ratio is υ=0.3.
Dimensionless central deflections are obtained and compared with the reference solutions
as given in Table 3. Here, the dimensionless central deflection is defined as w’ = 103 D w /
qa4 where 𝐷 = , Kw and Kp are the soil parameters in non-dimensional form
( )

11810
 Ülkü Hülya ÇALIK-KARAKÖSE

expressed as Kw =kw a4 / D and Kp =kp a2/ D. Results obtained in this study are very close to
the reference values, and it is observed that the deflections decrease with increasing shear
modulus as expected.

Table 3 - Dimensionless central deflections of a uniformly loaded square plate on

Pasternak foundation
Dimensionless central deflection
Kw Kp (w')
present Ref. [28] Ref. [46]
1 1 3.8536 3.8546 3.8530
1 81 0.7690 0.7630 0.7630
1 625 0.1162 0.1153 0.1150

Example 2:
A simply supported homogeneous square plate resting on two different foundation types,
Winkler and isotropic Pasternak foundations, is modelled using thick plate finite elements
and the deflection and moment values obtained at the plate centre are compared with those
given in [24]. The plate is subjected to a uniformly distributed load of q=30 kN/m2. Plate
length is a=10 m and the thickness is h=0.5m. elastic modulus and Poisson’s ratio are
E=200x106 kN/m2 and =0.167, respectively. The coefficient of subgrade reaction is
kw=5620 kN/m3 and the shear modulus of the Pasternak foundation is kp=120000 kN/m2.
Note that the shear moduli of the soil element are assigned to zero values to represent the
Winkler foundation.
The results given in Table 4 show that the central deflection and moment values obtained in
this study are very close to the reference values where the relative errors are between 2.5 - 3
%.

Table 4 - Central deflections and moments of a uniformly loaded square plate on Winkler
and Pasternak foundations
Winkler Pasternak
Present Present
Ref.[24] study Ref.[24] study
w x 104 (m) 5.329 5.498 4.187 4.298
M (kNm/m) 120.250 124.125 92.730 95.281

4.4. Simply Supported Homogeneous Plate on Orthotropic Pasternak Foundation

The dimensionless angular frequencies of the first 5 modes of a square homogeneous plate
resting on an orthotropic foundation with different material angles are obtained and
compared with the solution given by Kutlu [37] to verify the orthotropic foundation model.

11811
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

The plate is modelled using thick plate finite elements. Material and geometrical properties
of the homogeneous plate are E=25x106 kN/m2, υ=0.3, a=10 m and h=1 m. The
dimensionless coefficient of subgrade reaction and shear moduli are Kw=100, Kpx=10 and
Kpy =70, respectively.
The coefficient of subgrade reaction of the soil is kw =22893.773 kN/m3 and shear moduli
of the soil are kpx = 228937.73 kN/m2 and kpy =1602564.1 kN/m2. Dimensionless angular
frequencies are calculated using 𝜔 = where = 2550 kg/m3.

The dimensionless angular frequencies of the system obtained for 7 different material
angles are given comparatively in Table 5. It is seen that the values are very close to each
other which verifies the orthotropic Pasternak foundation model. Note that the
dimensionless angular frequency values are the same for complementary angles due to the
symmetry of the system in both geometry and loading.

Table 5 - Dimensionless angular frequencies of the first 5 modes for different material
angles
θ θ θ θ θ θ θ
Ref.[37] 3.5716 3.5616 3.5413 3.5310 3.5413 3.5616 3.5716
mode 1 Present
study 3.5371 3.5276 3.5082 3.4984 3.5082 3.5276 3.5371
Ref.[37] 5.7530 5.8207 5.9895 6.1075 5.9895 5.8207 5.7530
mode 2 Present
study 5.6807 5.7510 5.9266 6.0514 5.9266 5.7510 5.6807
Ref.[37] 7.1451 7.0716 6.8915 6.7678 6.8915 7.0716 7.1451
mode 3 Present
study 7.0546 6.9811 6.7995 6.6716 6.7995 6.9811 7.0546
[Ref.37] 9.0773 8.9905 8.8897 8.8529 8.8897 8.9905 9.0773
mode 4 Present
study 8.7947 8.7332 8.6499 8.6178 8.6499 8.7332 8.7947
Ref.[37] 9.5259 9.7388 10.1921 10.6600 10.1921 9.7388 9.5259
mode 5 Present
study 9.4438 9.6309 10.0631 10.5428 10.0631 9.6309 9.4438

5. SIMPLY SUPPORTED S-FGM PLATE ON ORTHOTROPIC PASTERNAK

FOUNDATION
The finite element analysis of a simply supported S-FGM square plate resting on an
orthotropic Pasternak foundation is performed next, using layered shell finite elements. The
elastic modulus at the bottom surface of the plate is E1= 25 x106 kN/m2, the Poisson’s ratio
is υ=0.3, the thickness is h=1 m and the length is a=10 m. The plate is subjected to a
uniformly distributed load of q=100 kN/m2. The dimensionless coefficient of subgrade

11812
 Ülkü Hülya ÇALIK-KARAKÖSE

reaction and shear moduli are Kw=100, Kpx=10 and Kpy=70, respectively where the material
angle is 0o.
The deflections of the S-FGM plate along the x direction (at y=5 m) for the cases with and
without membrane effects are obtained for E1/E2=2 and E1/E2=30. The plate is modelled
with 2 and 8 layers in order to observe the effect of the number of layers used defining the
variation of the elastic modulus of the S-FGM plate on deflections and it is seen that for all
cases, the deflection values do not change with the number of layers used, Figures 13-14.
Comparison of the deflections according to the inclusion and exclusion of membrane
effects is given in Figure 15 and it is observed that the deflection difference becomes more
apparent as the elastic modulus ratio increases or in other words, as the elastic modulus
decreases with a larger ratio from the bottom to the top surface of the plate.

x/a x/a
0,50 0,60 0,70 0,80 0,90 1,00 0,50 0,60 0,70 0,80 0,90 1,00
0,00 0,00

1,00 1,00
w x 10-4 (m)

2,00 2,00
w x 10-4 (m)

3,00 3,00

4,00 4,00
(E1/E2=2) 2 layers
5,00 5,00 (E1/E2=2) 2 layers
(E1/E2=2) 8 layers
(E1/E2=2) 8 layers
6,00 6,00

a) b)
Figure 13 - Comparison of deflections for different number of layers (E1/E2=2)
a) membrane effects are included b) membrane effects are excluded

x/a x/a
0,50 0,60 0,70 0,80 0,90 1,00 0,50 0,60 0,70 0,80 0,90 1,00
0,00 0,00
1,00 1,00
2,00
w x 10-4 (m)

2,00
w x 10-4 (m)

3,00 3,00
4,00 4,00
5,00 (E1/E2=30) 2 layers
5,00 (E1/E2=30) 2 layers
6,00 (E1/E2=30) 8 layers 6,00 (E1/E2=30) 8 layers
7,00 7,00

a) b)
Figure 14 - Comparison of deflections for different number of layers (E1/E2=30)
a) membrane effects are included b) membrane effects are excluded

11813
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

x/a
0,50 0,60 0,70 0,80 0,90 1,00
0,00

1,00

2,00
w x 10-4 (m)

3,00

4,00

5,00 E1/E2=30 (m.e. included)

E1/E2=30 (m.e. excluded)
6,00
E1/E2=2 (m.e. included)
7,00 E1/E2=2 (m.e. excluded)

Figure 15 - Comparison of deflections using 8 layers

The stress distributions through the thickness of the plate centre are investigated according
to the elastic modulus ratio and membrane effects using two layers, Figure 16. x and y
stresses are analogous but their magnitudes differ from each other due to the orthotropic
property of the soil material. The stress distributions for E1/E2=1 corresponding to a
homogeneous plate material are linear and the graphs for x and y gradually turn into
curved shapes towards the top surface, as E1/E2 ratio increases. Besides, stress distributions
diverge from each other for the cases with and without membrane effects as E1/E2 ratio
increases. Note that the stress values are obtained at the top and bottom levels of each layer
to be able to capture 3 stress values for the 2 layered solution.
E1/E2=1 (m.e. included)
E1/E2=1 (m.e. included)
-0,5 E1/E2=1 (m.e. excluded) -0,5
E1/E2=1 (m.e. excluded)
E1/E2=2 (m.e. included)
-0,4 E1/E2=2 (m.e. included)
E1/E2=2 (m.e. excluded) -0,4
E1/E2=2 (m.e. excluded)
E1/E2=10 (m.e. included) E1/E2=10 (m.e. included)
-0,3 -0,3
E1/E2=10 (m.e. excluded) E1/E2=10 (m.e. excluded)
E1/E2=30 (m.e. included) E1/E2=30 (m.e. included)
-0,2 -0,2
E1/E2=30 (m.e. excluded) E1/E2=30 (m.e. excluded)
-0,1 -0,1
-800 -600 -400 -200 0 200 400 600 800 -800 -600 -400 -200 0 200 400 600 800
0,0 0,0
z/h

z/h

0,1 0,1

0,2 0,2

0,3 0,3

0,4 0,4

0,5 0,5
x y

Figure 16 - x and y stress distributions for different E1/E2 ratios

11814
 Ülkü Hülya ÇALIK-KARAKÖSE

Then, the stress distributions are recalculated for different number of layers used in
modelling the S-FGM plate and shown comparatively in Figure 17. x and y stresses
obtained for 4 and 8 layers are almost the same while they differ slightly near the top and
bottom surfaces when 2 layers are used.
-0,5 -0,5

-0,4 E1/E2=30 (2 layers) -0,4 E1/E2=30 (2 layers)

E1/E2=30 (4 layers) E1/E2=30 (4 layers)
-0,3 -0,3
E1/E2=30 (8 layers) E1/E2=30 (8 layers)

-0,2 -0,2

-0,1 -0,1
-400 -200 0 200 400 600 800
-400 -200 0 200 400 600 800
0,0
0,0
z/h

0,1 z/h 0,1

0,2
0,2
0,3
0,3
0,4
0,4
0,5
0,5
x y

Figure 17 - x and y stress distributions for different number of layers (membrane effects
are included)

In addition, two alternate arrangements of the S-FGM plate on the orthotropic Pasternak
foundation are considered and solved for the stresses through the thickness of the plate for
comparison. The x and y stress components in Figures 18a-19a belong to the arrangement
where the surface with larger elastic modulus is attached to the elastic foundation whereas
Figures 18b-19b indicate the stress distributions of the plate in the reverse arrangement. It
is observed that x and y stresses start diverging towards the top surface as the elastic
modulus decreases rapidly and this situation is the reverse for the other case, Figures 18-19.
It is also obtained that zero x and y stresses which are at the mid-section of the
homogeneous plate, are below and above the mid-section for the two arrangements,
respectively due to the heterogeneity of the material.
x and y stress components are obtained for different material angles ranging from 0o to
90o at intervals of 15o and given comparatively in Figure 20. The elastic modulus ratio of
the S-FGM plate is taken as E1/E2=30 and the system is analysed for the case where
membrane effects are included. x and y stresses have the maximum divergence at the
bottom surface and after converging, they reach the top surface almost with a constant
difference.
Finally, the midpoint deflections are investigated for three different cases in which the non-
dimensional soil parameters are; Kw=100 Kpx=10 Kpy=70 (case 1) ; Kw=100 Kpx=10
Kpy=140 (case 2) ; Kw=100 Kpx=10 Kpy=280 (case 3) . Coefficient of subgrade reaction
and one of the shear moduli of the foundation are kept constant while the other shear
modulus is doubled in each case.

11815
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

-0,5 -0,5
E1/E2=1 -0,4 E2/E1=1
-0,4
E1/E2=2 E2/E1=2
-0,3 E1/E2=4 -0,3 E2/E1=4
E1/E2=10 E2/E1=10
-0,2 E1/E2=30 -0,2 E2/E1=30

-0,1 -0,1
-800 -600 -400 -200 0 200 400 600 800 -800 -600 -400 -200 0 200 400 600 800
0,0 0,0

z/h
z/h

0,1 0,1

0,2 0,2

0,3 0,3

0,4 0,4

0,5 0,5
x x

Figure 18 - x stress distributions for different E1/E2 and E2/E1 ratios (membrane effects
are included)

-0,5 -0,5
E1/E2=1 E2/E1=1
-0,4 -0,4
E1/E2=2 E2/E1=2
-0,3 E1/E2=4 -0,3 E2/E1=4
E1/E2=10 E2/E1=10
-0,2 E1/E2=30 -0,2 E2/E1=30

-0,1 -0,1
-800 -600 -400 -200 0 200 400 600 800 -800 -600 -400 -200 0 200 400 600 800
0,0 0,0
z/h

z/h

0,1 0,1

0,2 0,2

0,3 0,3

0,4 0,4

0,5 0,5
y y

Figure 19 - y stress distributions for different E1/E2 and E2/E1 ratios (membrane effects
are included)

Two different elastic modulus ratios for the plate and seven material angles for the
foundation are considered for all cases. It is observed in Figure 21 that the midpoint

11816
 Ülkü Hülya ÇALIK-KARAKÖSE

deflections decrease with increasing shear modulus. The maximum midpoint deflection
value is achieved when the material angle is 45o and the remaining values are in descending
order and are the same for complementary angles. Besides, the midpoint deflections for
E1/E2=30 are larger than the values obtained for E1/E2=2. In other words, midpoint
deflections increase for “softer” top surfaces. This increment becomes more apparent when
the membrane effects are excluded and with decreasing shear modulus, Figures 21-22.

0 degree 0 degree
-0,5 -0,5
15 degrees 15 degrees
-0,4 30 degrees -0,4 30 degrees
45 degrees 45 degrees
-0,3 -0,3 60 degrees
60 degrees
75 degrees 75 degrees
-0,2 -0,2
90 degrees 90 degrees
-0,1 -0,1
-400 -200 0 200 400 600 800 -400 -200 0 200 400 600 800
0,0 0,0
z/h

z/h

0,1 0,1

0,2 0,2

0,3 0,3

0,4 0,4

0,5 0,5
x y

Figure 20 - x and y stress distributions for different material angles (E1/E2=30,

membrane effects are included)

material angle material angle

E1/E2=2 E1/E2=2
0 15 30 45 60 75 90 (case 1) 0 15 30 45 60 75 90 (case 1)
0,00 E1/E2=30 0,00 E1/E2=30
(case 1) (case 1)
1,00 1,00
E1/E2=2 E1/E2=2
w x 10-4 (m)

(case 2) 2,00
w x 10-4 (m)

2,00 (case 2)
E1/E2=30 E1/E2=30
3,00 3,00
(case 2) (case 2)
4,00 E1/E2=2 4,00
E1/E2=2
(case 3) (case 3)
5,00 5,00
E1/E2=30
E1/E2=30
6,00 (case 3) 6,00 (case 3)

7,00 7,00

a) b)
Figure 21 - Mid-point deflections for different material angles
a) membrane effects are included b) membrane effects are excluded

11817
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

material angle
0 15 30 45 60 75 90
0,00

1,00

2,00
w x 10-4 (m)

m.e. included (case 1)

3,00
m.e. excluded (case 1)
4,00 m.e. included (case 2)

5,00 m.e. excluded (case 2)

m.e. included (case 3)
6,00
m.e. excluded (case 3)
7,00

Figure 22 - Mid-point deflections for different material angles (E1/E2=30)

6. CONCLUSIONS
In this study, the finite element analysis of S-FGM plates with and without membrane
effects, resting on orthotropic Pasternak elastic foundations with different material angles is
carried out using SAP2000 software package. It is demonstrated with the verification
examples that the proposed model is convenient and robust. This work can be easily
extended to static and dynamic analyses of FGM plates with various geometries resting on
arbitrarily orthotropic Pasternak foundations for further studies.
The results obtained in the present study can be summarised as follows:
1) For small elastic modulus ratios (e.g. E1/E2=2), the deflections of the S-FGM plate both
with and without Pasternak foundation do not change with the number of layers used for
the plate.
2) For large elastic modulus ratios (e.g. E1/E2=30), the deflections of the plate on Pasternak
foundation do not change with the number of layers used whereas the deflections of the
plate without Pasternak foundation increase at a lower rate as the number of layers
increases.
3) The deflections increase when the membrane effects are excluded and the deflection
difference between the cases with and without membrane effects becomes more apparent as
the elastic modulus decreases with a larger ratio from the bottom to the top surface of the
plate.
 Stresses x and y through the thickness of the plate centre are analogous but their
magnitudes differ from each other due to the orthotropic property of the soil material. All
stress distributions for a homogeneous plate material are linear while the distributions for
x and y gradually turn into curved shapes towards the top surface as E1/E2 ratio increases.

11818
 Ülkü Hülya ÇALIK-KARAKÖSE

 The stress distributions diverge for the cases with and without membrane effects as
E1/E2 ratio increases.
 Stresses x and y obtained for 4 and 8 layers are almost the same while they differ
slightly near the top and bottom surfaces when 2 layers are used.
 Stresses x and y start diverging towards the top surface as the elastic modulus
decreases with a larger rate for the arrangement of the S-FGM plate where the surface with
larger elastic modulus is attached to the orthotropic Pasternak foundation. This divergence
is towards the bottom surface for the opposite positioning of the plate. Zero x and y
stresses which are at the mid-section of the homogeneous plate, are below and above the
mid-section for the two arrangements, respectively due to the heterogeneity of the material.
8) For different material angles, x and y stresses have the maximum divergence at the
bottom surface and after converging, they reach the top surface almost with a constant
difference.
9) The midpoint deflections decrease with increasing shear modulus of the foundation. The
maximum midpoint deflection value is achieved when the material angle is 45o and the
remaining values are in descending order and are the same for complementary angles.
10) The midpoint deflections increase for softer top surfaces and this increment becomes
more apparent with decreasing shear modulus and when the membrane effects are
excluded.

Acknowledgements
Engin Orakdöğen and Kutlu Darılmaz are gratefully acknowledged for their fruitful
suggestions on this study.

References
[1] Suresh S., Mortenson A., Fundamentals of Functionally Graded Materials: Processing
and Thermomechanical Behavior of Graded Metals and Metal-Ceramic Composites.
IOM Communications Ltd, London, UK, 1998.
[2] Koizumi M., FGM activities in Japan. Composites Part B: Eng 28(1-2), 1-4, 1997.
[3] Tanigawa Y., Some basic thermoelastic problems for nonhomogeneous structural
materials. Appl. Mech. Rev. 48(6), 287-300, 1995.
[4] Suresh S., Mortensen A., Functionally graded metals and metalceramic composites 2:
thermomechanical behaviour. Int. Mater. Rev. 42(3), 85-116, 1997.
[5] Reddy J.N., Analysis of functionally graded plates. Int. J. Numer. Methods Eng., 47,
663-684, 2000.
[6] Cheng Z.Q., Batra R.C. Deflection relationships between the homogenous Kirchhoff
plate theory and different functionally graded plate theories. Archives of Mechanics
52, 143-158, 2000.

11819
FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

[7] Bao G., Wang L., Multiple cracking in functionally graded ceramic/metal coatings.
Int. J. Solids Struct., 32, 2853-2871, 1995.
[8] Jin Z.H., Paulino G.H., Transient thermal stress analysis of an edge crack in a
functionally graded material. Int. J. Fracture, 107, 73-98, 2001.
[9] Delale F., Erdogan F., The crack problem for a nonhomogeneous plane. ASME J.
Appl. Mech., 50, 609-614, 1983.
[10] Erdogan F., Wu B.H., Crack problems in FGM layers under thermal stresses. J.
Therm. Stress., 19, 237-265, 1996.
[11] Chung Y.L., Chi S.H., The residual stress of functionally graded materials. J. Chin.
Inst. Civ. Hydraul. Eng., 13, 1-9, 2001.
[12] Chi S.H., Chung Y.L., Cracking in sigmoid functionally graded coating. J. Mech., 18,
41-53, 2002.
[13] Chi S.H., Chung Y.L., Mechanical behavior of functionally graded material plates
under transverse load-Part I: Analysis. Int. J. of Solids and Str., 43, 3657-3674, 2006a.
[14] Chi S.H., Chung Y.L., Mechanical behavior of functionally graded material plates
under transverse load-Part II: Numerical results. Int. J. of Solids and Str., 43, 3675-
3691, 2006b.
[15] Orakdöğen E., Küçükarslan S., Sofiyev A., Omurtag M.H., Finite element analysis of
functionally graded plates for coupling effect of extension and bending. Meccanica,
45, 63-72, 2010.
[16] Zenkour A.M., Generalized shear deformation theory for bending analysis of
functionally graded plates. Appl. Math. Model., 30, 67-84, 2006.
[17] Elishakoff I., Gentilini C., Viola E., Three-dimensional analysis of an all-round
clamped plate made of functionally graded materials. Acta Mech., 180(1-4), 21-36,
2005.
[18] Swaminathan K., Naveenkumar D.T., Zenkour A.M., Carrera E., Stress, vibration and
buckling analyses of FGM plates-A state-of-the-art Review. Composite Structures,
120, 10-31, 2015.
[19] Winkler E. Theory of Elasticity and Strength of Materials. Dominicus, Prague, 1867.
[20] Pasternak P.L., On a new method of analysis of an elastic foundation by means of two
foundation constants. Cosudarstrennoe Izdatelstvo Literaturi po Stroitelstvu i
Arkhitekture, Moscow, USSR, 1-56, 1954.
[21] Vallabhan C.V.G., Straughan W.T., Das Y.C., Refined Model of Analysis of Plates of
Elastic Foundations. Journal of Engineering Mechanics, 117(12), 2830-2844, 1994.
[22] Celik M., Saygun A., A method for the analysis of plates on a two-parameter
foundation. Int. J. of Solids and Str., 36, 2891-2915, 1999.
[23] Xiang Y., Wang C.M., Kitipornchai S., Exact vibration solution for initially stressed
Mindlin plates on Pasternak foundation. Int. J. Mech. Sci., 36(4), 311-316, 1994.

11820
 Ülkü Hülya ÇALIK-KARAKÖSE

[24] Omurtag M.H., Ozutok A., Akoz A.Y., Free vibration analysis of Kirchhoff plates
resting on elastic foundation by mixed finite element formulation based on Gateaux
differential. Int. J. Numer. Methods Eng., 40(2), 295-317, 1997.
[25] Matsunaga H. Vibration and stability of thick plates on elastic foundations. ASCE J.
Eng. Mech., 126(1), 27-34, 2000.
[26] Zhou D., Cheung Y.K., Lo S.H., Au F.T.K., Three-dimensional vibration analysis of
rectangular thick plates on Pasternak foundations. Int. J. Numer. Methods. Eng.,
59(10), 1313-1334, 2004.
[27] Shen H.S., Postbuckling analysis of composite laminated plates on two-parameter
elastic foundations. Int. J. Mech. Sci., 37(12), 1307-131, 1995.
[28] Huang Z.Y., Lü C.F., Chen W.Q., Benchmark solutions for functionally graded thick
plates resting on Winkler–Pasternak elastic foundations. Composite Structures, 85,
95-104, 2008.
[29] Lee W.H., Han S.C., Park W.T., A refined higher order shear and normal deformation
theory for E-,P-, and S-FGM plates on Pasternak elastic foundation. Composite
Structures, 122, 330-342, 2015.
[30] Tajeddini V., Ohadi A., Sadighi M., Three-dimensional free vibration of variable
thickness thick circular and annular isotropic and functionally graded plates on
Pasternak foundation. Int. J. Mech. Sci., 53(4), 300-308, 2011.
[31] Baferani A.H., Saidi A.R., Ehteshami H., Accurate solution for free vibration analysis
of functionally graded thick rectangular plates resting on elastic foundation.
Composite Structures, 93(7), 1842-1853, 2011.
[32] Malekzadeh P., Golbahar Haghighi M.R., Atashi M.M., Free vibration analysis of
elastically supported functionally graded annular plates subjected to thermal
environment. Meccanica, 46, 893-913, 2011.
[33] Thai H.T., Choi D.H., A refined shear deformation theory for free vibration of
functionally graded plates on elastic foundation. Composites Part B, 43, 2335-2347,
2012.
[34] Mansouri M.H., Shariyat M., Differential quadrature thermal buckling analysis of
general quadrilateral orthotropic auxetic FGM plates on elastic foundations. Thin-
Walled Structures, 112, 194-207, 2017.
[35] Hamarat M.A., Dynamic analysis of structures resting on two parameter elastic
foundation. MSc Dissertation, Istanbul Technical University, Istanbul, 2012.
[36] Hamarat M.A., Çalık-Karaköse Ü.H., Orakdöğen E., Seismic Analysis of Structures
Resting on Two Parameter Elastic Foundation. 15th World Conference on Earthquake
Engineering. Lisbon, Portugal, September, 2012.

11821
 FE Analysis of FGM Plates on Arbitrarily Orthotropic Pasternak Foundations for …

[37] Kutlu A., Analysis of free vibrations of Mindlin plates resting on arbitrarily
orthotropic Pasternak foundation with mixed finite elements. MSc Dissertation,
Istanbul Technical University, Istanbul, 2007.
[38] Kutlu A., Omurtag M.H., Large deflection bending analysis of elliptic plates on
orthotropic elastic foundation with mixed finite element method, International Journal
of Mechanical Sciences, 65, 64–74, 2012.
[39] Kutlu A., Ugurlu B., Omurtag M.H., Ergin A., Dynamic response of Mindlin plates
resting on arbitrarily orthotropic Pasternak foundation and partially in contact with
fluid. Ocean Engineering, 42,112-125, 2012.
[40] Kolahchi R., Safari M., Esmailpour M., Dynamic stability analysis of temperature-
dependent functionally graded CNT-reinforced visco-plates resting on orthotropic
elastomeric medium. Composite Structures, 150, 255-265, 2016.
[41] Heydari M. M., Bidgoli A. H., Golshani H. R., Beygipoor G., Davoodi, A., Nonlinear
bending analysis of functionally graded CNT-reinforced composite Mindlin
polymeric temperature-dependent plate resting on orthotropic elastomeric medium
using GDQM. Nonlinear Dynamics, 79(2), 1425- 1441, 2015.
[42] Arani A. G., Cheraghbak A., Kolahchi R., Dynamic buckling of FGM viscoelastic
nanoplates resting on orthotropic elastic medium based on sinusoidal shear
deformation theory. Structural engineering and mechanics: An international journal,
60(3), 489-505, 2016.
[43] SAP2000, v18., Integrated Finite Elements Analysis and Design of Structures,
Computers and Structures. p. Inc, Berkeley, CA, 2016.
[44] Elmacı Ş., Spectral analysis and determination of free vibration characteristics of
quadratic plates resting on arbitrarily orthotropic two parameter elastic foundation.
MSc Dissertation, Istanbul Technical University, Istanbul, 2019.
[45] Aykılıç B., Spectral analysis of circular and elliptic plates resting on arbitrarily
orthotropic Pasternak type foundation. MSc Dissertation, Istanbul Technical
University, Istanbul, 2019.
[46] Lam K.Y., Wang C.M., He X.Q., Canonical exact solution for Levy-plates on two
parameter foundation using Green’s functions. Eng. Struct., 22(4),364-378, 2000.

11822

View publication stats