Next Article in Journal
Application of the ADMM Algorithm for a High-Dimensional Partially Linear Model
Next Article in Special Issue
A Closed-Form Parametrization and an Alternative Computational Algorithm for Approximating Slices of Minkowski Sums of Ellipsoids in R3
Previous Article in Journal
Collaborative Screening of COVID-19-like Disease from Multi-Institutional Radiographs: A Federated Learning Approach
Previous Article in Special Issue
Numerical Investigation of Heat Transfer Enhancement in a Microchannel with Conical-Shaped Reentrant Cavity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Approach for Mechanical Behaviour Analysis of FGM Plates on Elastic Foundation

by
Fatima Zohra Zaoui
1,
Djamel Ouinas
1,
Belkacem Achour
2,*,
Mabrouk Touahmia
2,
Mustapha Boukendakdji
2,
Enamur R. Latifee
2,
Ahmed A. Alawi Al-Naghi
2 and
Jaime Aurelio Viña Olay
3
1
Laboratory of Science and Technology Environment and Valorization, Faculty of Sciences and Technology, Ibn Badis University, Mostaganem 27000, Algeria
2
Civil Engineering Department, University of Ha’il, Ha’il 55476, Saudi Arabia
3
Department of Materials Science and Metallurgical Engineering, University of Oviedo, Viesques Campus, 33203 Gijón, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(24), 4764; https://doi.org/10.3390/math10244764
Submission received: 17 October 2022 / Revised: 28 November 2022 / Accepted: 6 December 2022 / Published: 15 December 2022
(This article belongs to the Special Issue Mathematical Problems in Mechanical Engineering)

Abstract

:
This paper presents the flexural analysis of functionally graded plates resting on elastic foundations using new two-dimensional (2D) and quasi-three-dimensional (quasi-3D) higher order shear deformation theories. The main interesting feature of this theory is that it proposes a new displacement field with undetermined integral variables which involves only five unknown functions, unlike other shear and normal deformation theories, hence making it easier to use. A parabolic transverse shear deformation shape function satisfying the zero shear stress conditions on the plate outer surfaces is considered. The elastic foundation follows the Pasternak mathematical model. The material properties change continuously across the thickness of the FG plate using different distributions: power law, exponential, and Mori–Tanaka models. The governing equations of FG plates subjected to sinusoidal and uniformly distributed loads are established through the principle of virtual works and then solved via Navier’s procedure. In this work, a detailed discussion on the influence of material composition, geometric parameters, stretching effect, and foundation parameters on the deflection, axial displacements, and stresses is given, and the obtained results are compared with those published in previous works to demonstrate the accuracy and the simplicity of the present formulations. The different obtained results were found to be in good agreement with the available solutions of other higher-order theories. The proposed model is able to represent the cross section warping in the deformed shape and to demonstrate the validity and efficiency of the approach, the findings reported herein prove that this theory is capable of predicting displacements and stresses more accurately than other theories, as its results are closer when compared to numerical methods reported in other literatures.

1. Introduction

Functionally graded materials (FGMs) are a new type of advanced composite materials made up of a combination of ceramics and metals that is designed to endure high thermal environment and to provide stiffness to the host structure [1]. Their mechanical properties are assumed to change following a smooth and continuous distribution from the lower to the upper surface of the structure and thus eliminate stress concentration found in laminate composites and maintain structural integrity at a desirable level [2]. Due to this feature, FGMs have been explored in many engineering applications, such as aerospace, aircraft, automobile, and defence industries and most recently, electronic and biomedical applications [3,4].
In recent years, due to widespread applications of FGMs in various engineering fields, many theories have been developed in order to analyse bending, buckling and dynamic behaviours of different FG structures [5,6,7,8,9]. The use of plates resting on elastic foundations is very popular in structural engineering, such as reinforced concrete pavements, airport runways, storage tanks foundations, swimming pools, retaining walls, and concrete footings [10,11,12,13]. The most famous model of elastic foundations is known as the Winkler–Pasternak model. This contains, of course, two parameters—Winkler’s transverse stiffness coefficient and Pasternak’s shear stiffness coefficient. The first parameter of the Pasternak model represents the modelling of the elastic medium as a series of closely spaced, mutually independent vertical linear elastic springs, while the second parameter represents the transverse shear stress due to the interaction of the shear deformation of the surrounding elastic medium [14,15,16,17,18].
The bending problem of transverse load acting on an isotropic inhomogeneous rectangular plate using both two-dimensional (2D) trigonometric and three-dimensional (3D) elasticity solutions was analysed by Zenkour [19]. Additionally, the stress and displacement response of the FG plates resting on Pasternak elastic foundations have been analysed under uniform thermal and mechanical loading by Zenkour [20]. The effects of thickness stretching in functionally graded plates and shells were studied in a paper by Carrera et al. [21], in which variable plate/shell models are implemented according to Carrera’s unified formulation. Quasi-3D sinusoidal and hyperbolic shear deformation theories to study the static and free vibration responses of FG plates were proposed by Neves et al. [22,23]. Based on a simple refined higher-order shear deformation theory with two variables, Thai et al. [24] examined the bending, buckling and free vibration of thick FG plates resting on elastic foundation. The buckling response of a composite laminated material square plate having an elliptical notch was studied by Ouinas and Achour [25] using the finite element method. Mantari and Soares [26] developed an analytical solution to the static analysis of functionally graded plates by using a new trigonometric higher-order theory in which the stretching effect is included. A 3D elasticity solution for bending thick FG plates was proposed by Zhang et al. [27] using a hybrid semi-analytical approach—the state-space-based differential quadrature method. The influence of temperature and moisture on the bending behaviour of FG plates resting on elastic foundations using refined shear and normal deformations plate theory were investigated by Al Khateeb and Zenkour [28]. Thai and Choi [29] used a zeroth-order shear deformation theory for bending and vibration analyses of functionally graded plates resting on elastic foundation. A refined higher-order shear and normal deformation theory for E-, P-, and S-FGM plates on Pasternak elastic foundation was presented by Lee et al. [30]. Akavci and Tanrikulu [31] presented 2D and quasi-3D shear deformation theories for bending and free vibration analysis of single-layer FG plates using a new hyperbolic shape function. Mantari and Granados [32] developed an original first shear deformation theory for bending and free vibration analysis of functionally graded plates resting on elastic foundation. Houari et al. [33] developed a new simple and accurate three-unknowns sinusoidal shear deformation theory for the bending and vibration analysis of FG plates. Recently, a bending analysis of different material distributions of functionally graded beam was established by Aldousari [34], in which two symmetric and anti-symmetric functions were developed and their effects on the static deflection and bending stresses were compared with classical power law distribution. A non-polynomial four-refined-shear deformation theory for free vibration analysis of FG plates on elastic foundation was performed by Meftah et al. [35]. Zaoui et al. [36] and Guerroudj et al. [37] presented hybrid quasi-3D shear deformation theories to study the free vibration of functionally graded plates. A new simple shear deformation theory with three unknown was developed by Amar et al. [38] for static analysis of FG plates on elastic foundation. Younsi et al. [39] developed a non-polynomial 2D and quasi-3D theory to investigate static and dynamic responses of FG plates. Belkhodja et al. [40] investigated the flexion, free vibrations and buckling of FGMs plate with simply supported edges using a new exponential-trigonometric shear function. Thermal vibration of functionally graded graphene platelets reinforced composite (FG-GPLRC) annular plate resting on an elastic foundation under the mechanical load framework of higher order shear deformation (HSDT) was analysed by Wang et al. [41]. Many researchers [42,43,44,45] have studied the flows of MHD hybrid nanofluid, while Qureshi et al. [46] have specially investigated the dispersion of polymer matrix/CNT nanocomposite material through permeable surfaces subjected to the magnetized hybrid nanofluids. The deduced mathematical model for nanocomposites based on PDEs is beneficial for nanotechnology. The analytical solutions for the temperature and velocity fields are obtained using the Laplace transform. The discrete fractional proportional operator, the Caputo and Caputo–Fabrizio fractional operators have been studied in references [47,48]. Vieru et al. [49] used the fractional time derivative of Caputo–Fabrizio to develop a new mathematical model. The latter is based on a fractional constitutive equation of the clean heat flux to describe the memory effects. Saffari et al. [50] used a size-dependent non-local strain gradient theory to analyse free and forced vibration responses of fluid-conveying viscoelastic nanotubes sitting on nonlinear Visco–Pasternak foundations under magnetic fields. Alazwari et al. [51] used the Visco–Winkler–Pasternak elastic foundation model in conjunction with a quasi-3D refined theory for the first time to study the vibration response of FG plates.
The objective of this work is to study the static behaviour of FG plates simply supported and resting on elastic foundations of the Winkler–Pasternak type by applying modified 2D and quasi-3D shear strain theories developed by Zaoui et al. [52], for which only the problem of free vibration of FGM plates resting on elastic supports was assessed. The most interesting feature of these theories is that they have a new displacement field that contains fewer unknowns compared to other theories such as Neves et al. [22,23]. Indeterminate integrals have been used in the proposed field of displacements in which the transverse shear and thickness stretching effects are considered in quasi-3D theory while ignored in the 2D model. Moreover, these theories account for a parabolic variation of transverse shear stresses across the thickness to satisfy shear stress free boundary conditions without including shear correction factor. The principle of virtual work is applied to establish the fundamental governing equations of FG plates subjected to sinusoidal and uniformly distributed loads. Closed-form solutions of deflections for simply supported plates are derived via Navier’s method, and the found results are compared to those published in the literature to demonstrate accuracy and simplicity of the proposed theories. The effects of power law index, slenderness ratio, side-to-thickness ratio, normal strain, Winkler–Pasternak parameters, and different rules of mixture on the bending response of FG plates have been carried out. The findings presented here demonstrate that the suggested model can accurately describe the cross section warping in the deformed shape and that it can estimate displacements and stresses of FG plates more precisely than other shear deformation theories.

2. Mathematical Formulation

2.1. Material Properties

In this work, a functionally graded rectangular plate, having uniform thickness h , length a , width b , is considered as seen in Figure 1. The material properties of FGM are assumed to vary gradually through the thickness according to the following rules of mixture.

2.1.1. The Power Law (P-FGM) Variation

The volume fraction of the P-FGM plate is supposed to change continuously across the thickness of the plate according to the power law variation [53], as given in Equation (1):
P ( z ) = P m + ( P c P m )   ( 1 2 + z h ) k

2.1.2. The Exponential (E-FGM) Variation

The volume fraction of the E-FGM plate is considered to vary continuously across the thickness direction of the plate according to the exponential distribution [30,54] as follows:
P ( z ) = P m e p ( z + h / 2 ) , p = 1 h ln ( P c P m )

2.1.3. The Mori–Tanaka Homogenization Model

For the Mori–Tanaka scheme [29,39], the volume fraction of the FGM plate is defined in the equations below:
P ( z ) = P m + ( P c P m ) V c 1 + V m ( P c P m 1 ) 1 + ν 3 3 ν ,  
V c = ( 1 2 + z h ) k , V m + V c = 1
where P represents the effective material property, such as Young’s modulus E . P m   P m , P c , V m and V c denote material properties and the volume fraction of the upper (metal) and the lower (ceramic) faces of the plate, respectively. k is the power law index and P is the volume fraction exponent. Since the effects of Poisson’s ratio ( ν ) on the response of FG plates are very small, they are supposed to be constant for all FGM graded layers.

2.2. Kinematics

On the basis of higher-order shear deformation theories and considering the stretching effect [52,55], the supposed displacement field of the plate can be expressed as
u ( x , y , z , t ) = u 0 ( x , y , t ) z w 0 x + k 1   f ( z ) θ   ( x , y , t )   d x v ( x , y , z , t ) = v 0 ( x , y , t ) z w 0 y + k 2   f ( z ) θ   ( x , y , t )   d y w ( x , y , z , t ) = w 0 ( x , y , t ) + g ( z ) ϕ z ( x , y , t )
where u 0 , v 0 , w 0 , θ and ϕ z are the five unknown displacement functions of the middle surface of the plate. Note that the integrals do not have limits. The constants k 1 and k 2 depend on the geometry. f ( z ) is the shape function which represents the variation of the transverse shear strains within the thickness. Note that g ( z ) = 0 for 2D analysis.
The shape function is chosen according to Zaoui et al. [9,52] as
f ( z ) = π h π 4 + h 4 e ( h z / π ) ( π 2 sin   ( π z h ) + h 2 cos   ( π z h ) ) π h 3 π 4 + h 4
and
g ( z ) = d f d z
The general strain–displacement relations can be defined from Equation (4) by the application of the linear small-strain elasticity theory as follows
{ ε x ε y γ x y } = { ε x 0 ε y 0 γ x y 0 } + z { k x b k y b k x y b } + f ( z ) { k x s k y s k x y s }
{ γ y z γ x z } = g ( z ) { γ y z 0 γ x z 0 } , ε z = g ( z )   ε z 0
where
{ ε x 0 ε y 0 γ x y 0 } = { u 0 x v 0 x u 0 y + v 0 x } ,   { k x b k y b k x y b } = { 2 w 0 x 2 2 w 0 y 2 2 2 w 0 x y }
{ k x s k y s k x y s } = { k 1 θ k 2 θ k 1 y θ   d x + k 2 x θ   d y } ,
{ γ y z 0 γ x z 0 } = { k 2 θ   d y + ϕ z y k 1 θ   d x + ϕ z x } ,
ε z 0 = ϕ z   and   g ( z ) = d g ( z ) d z
It can be seen from Equation (6) that the transverse shear strains ( γ x z , γ y z ) are equal to zero at the top ( z = h / 2 ) and bottom ( z = h / 2 ) surfaces of the plate. A shear correction coefficient is, hence, not required.
The integrals employed in the above expressions are solved by a Navier’s type solution and can be described as
y θ   d x = A 2 θ x y ,   x θ   d y = B 2 θ x y
θ   d x = A θ x ,   θ   d y = B θ y
where the coefficients A and B —which are defined according to the type of solutions employed, k 1 and k 2 —are expressed as follows:
A = 1 α 2 ,   B = 1 β 2 ,   k 1 = α 2 ,   k 2 = β 2
Note that α and β are defined in Equation (24).
The stress–strain relationships for FG plate according to three-dimensional (3D) elasticity can be written as
{ σ x σ y σ z τ y z τ x z τ x y } = [ Q 11 Q 12 Q 13 0 0 0 Q 12 Q 22   Q 23 0 0 0 Q 13 Q 23 Q 33 0 0 0 0 0 0 Q 44 0 0 0 0 0 0 Q 55 0 0 0 0 0 0 Q 66 ] { ε x ε y ε z γ y z γ x z γ x y }
The Q i j expressions in terms of engineering constants depend on the normal strain ε z .
  • If the εz ≠ 0 then Qij are
    Q 11 = Q 22 = Q 33 = ( 1 ν ) E ( z ) ( 1 2 ν ) ( 1 + ν )
    Q 12 = Q 13 = Q 23 = ν   E ( z ) ( 1 2 ν ) ( 1 + ν )
    Q 44 = Q 55 = Q 66 = E ( z ) 2 ( 1 + ν )
  • If the εz = 0, then Qij are
    Q 11 = Q 22 = E ( z ) ( 1 ν 2 )
    Q 12 = ν   E ( z ) ( 1 ν 2 )
    Q 44 = Q 55 = Q 66 = E ( z ) 2 ( 1 + ν )

2.3. Equilibrium Equations and Stress Components

The equilibrium equations of functionally graded plate subjected to mechanical loadings can be derived by applying the principle of virtual works [56], which can be expressed in analytical form as
0 = δ U + δ V p + δ V e
where δ U is the virtual strain energy, δ V p is the virtual potential energy of applied distributed transverse loads, and δ V e is the virtual strain energy of the elastic foundation.
The virtual strain energy of the plate is defined by
δ   U = V [ σ x   δ   ε x + σ y   δ   ε y + σ z   δ   ε z + τ x y   δ   γ x y + τ y z   δ   γ y z + τ x z   δ   γ x z ]   d V     = A [ N x   δ   ε x 0 + N y   δ   ε y 0 + N z   δ   ε z 0 + N x y   δ   γ x y 0 + M x b   δ   k x b + M y b   δ   k y b     + M x y b   δ   k x y b + M x s   δ   k x s + M y s   δ   k y s + M x y s   δ   k x y s + Q y z s   δ   γ y z 0 + Q x z s   δ   γ x z 0 ]   d A
where A is the top surface and the stress resultants N, M and Q are given by
( N i , M i b , M i s ) = h / 2 h / 2 ( 1 , z , f ) σ i d z ,   ( i = x , y , x y )
N z = h / 2 h / 2 g ( z )    σ z   d z  
( Q x z s , Q y z s ) = h / 2 h / 2 g ( τ x z , τ y z ) d z
The virtual potential energy of the applied external loads to the plate can be expressed as
δ V p = A q δ   w 0   d   A
where q is the external load applied to the plate.
The virtual strain energy of the foundation can be written as
δ V e = A [ K w w 0   δ   w 0 + K s (   w 0 x δ   w 0 x +   w 0 y δ   w 0 y ) ]   d A
where K w and K s are the transverse and shear stiffness coefficients of the foundation, respectively.
The governing equations of equilibrium can be obtained by substituting Equations (14), (16) and (17) into Equation (13), integrating the displacement gradients by parts, and setting the coefficients δ u 0 , δ v 0 , δ w 0 , δ θ and δ ϕ z to zero separately, as follows:
δ u 0 : N x x + N x y y = 0
δ v 0 : N y y + N x y x =   0
δ w 0 : 2 M x b x 2 + 2 M y b y 2 + 2 2 M x y b x y K w w 0 + K s 2 w 0 =   q
δ θ : k 1 M x s k 2 M y s ( k 1 A + k 2 B ) 2 M x y s x y + k 1 A S x z s x + k 2 B S y z s y = 0
δ ϕ z : S x z s x + S y z s y N z = 0
where 2 = ( 2 / x 2 ) + ( 2 / y 2 ) is the Laplacian operator in a 2D Cartesian coordinate system.
The stresses and moment resultants which appeared in Equation (15) can be written in terms of generalized displacements ( u 0 , v 0 , w 0 , θ , ϕ z ) as:
{ N x N y N x y M x b M y b M x y b M x s M y s M x y s N z } = [ A 11 A 12 0 B 11 B 12 0 B 11 s B 12 s 0 X 13 A 12 A 22 0 B 12 B 22 0 B 12 s B 22 s 0 X 23 0 0 A 66 0 0 B 66 0 0 B 66 s 0 B 11 B 12 0 D 11 D 12 0 D 11 s D 12 s 0 Y 13 B 12 B 22 0 D 12 D 22 0 D 12 s D 22 s 0 Y 23 0 0 B 66 0 0 D 11 0 0 D 66 s 0 B 11 s B 12 s 0 D 11 s D 12 s 0 H 11 s H 12 s 0 Y 13 s B 12 s B 22 s 0 D 12 s D 22 s 0 H 12 s H 22 s 0 Y 23 s 0 0 B 66 s 0 0 D 66 s 0 0 H 66 s 0 X 13 X 23 0 Y 13 Y 23 0 Y 13 s Y 23 s 0 Z 33 ] { ε x 0 ε y 0 γ x y 0 k x b k y b k x y b k x s k y s k x y s ε z 0 }
{ S y z s S x z s } = [ A 44 s 0 0 A 55 s ]   { γ y z 0 γ x z 0 }
where the stiffness components and inertias are given as:
( A i j , A i j s , B i j , D i j , B i j s , D i j s , H i j s ) = h / 2 h / 2 Q i j ( 1 , g 2 ( z ) , z , z 2 , f ( z ) , z   f ( z ) , f 2 ( z ) ) d z
( X i j ,   Y i j , Y i j s , Z i j ) = h / 2 h / 2 ( 1 , z , f ( z ) , g ( z ) ) g ( z )   Q i j d z
( I 0 ,    I 1 ,   I 2 ,   J 1 ,    J 2 ,   J 0 ,   K 0 ,   K 2 ) = h / 2 h / 2 ( 1 ,   z ,   z 2 ,   f ,   z f ,   g ,   g 2 , f 2 ) ρ ( z )     d z  
and ρ ( z ) is the mass density.
The in-plane normal and shear stresses ( σ x , σ y and τ x y ) can be determined accurately by the constitutive Equation (10) for FG plates. However, if the transverse normal and shear stresses ( σ z , τ y z and τ x z ) are computed from the constitutive Equation (10), they may not respect the boundary conditions at the upper and lower surfaces of the plate. Therefore, these stresses are determined by integrating the equilibrium equations of 3D elasticity with respect to thickness coordinates as:
τ x z = h / 2 z ( σ x x + τ x y y )   d z + C 1 ( x , y )
τ y z = h / 2 z ( τ x y x + σ y y )   d z + C 2 ( x , y )
σ z = h / 2 z ( h / 2 z [ 2 σ x x 2 + 2 2 τ x y x y + 2 σ y y 2 ] d z )   d z + C 3 ( x , y ) z + C 4 ( x , y )
where C i ’s (i = 1, 4) are constants and are determined by the following boundary conditions at the upper and lower surfaces of the plate.
τ x z | z = ± h / 2 = 0 ,   τ y z | z = ± h / 2 = 0 , σ z | z = h / 2 = q ( x , y ) , σ z | z = h / 2 = 0

3. Closed-Form Solutions for Simply Supported FG Plates

The Navier’s procedure [18,57], based on double Fourier series, is employed herein to define the closed-form solution of the partial differential equations (Equation (18)) for which the displacement variables satisfying the boundary conditions can be given as
{ u 0 v 0 w 0 θ ϕ z } = m = 1 n = 1 { U m n cos ( α   x ) sin ( β   y ) V m n sin ( α   x ) cos ( β   y ) W m n sin ( α   x ) sin ( β   y ) X m n sin ( α   x ) sin ( β   y ) ϕ m n sin ( α   x ) sin ( β   y ) }
where ( U m n , V m n , W m n , X m n , ϕ m n ) are unknown functions to be determined. α and β are expressed as:
α = m π / a   and   β = n π / b
The transverse distributed load q(x, y) is also expanded in a double Fourier series as
q ( x , y ) = m = 1 n = 1 q m n sin ( m π a x ) sin ( n π b y )
The coefficients qmn are given below for some general loadings:
  • For uniformly distributed load
    q m n = { 16 q 0 m n π 2 , m , n = 1 , 3 , 5 , . . 0 m , n = 2 , 4 , 6 , . .
  • For sinusoidal distributed load
    q m n = q 0
    in which q0 is the intensity of the load.
Substituting Equations (19) and (23) into Equation (18), closed-form solutions of static problem of FG plate can be obtained as follows:
[ S 11 S 12 S 13 S 14 S 15 S 12 S 22 S 23 S 24 S 25 S 13 S 23 S 33 S 34 S 35 S 14 S 24 S 34 S 44 S 45 S 15 S 25 S 35 S 45 S 55 ] { U m n V m n W m n X m n ϕ m n } = { 0 0 q m n 0 0 }
in which Sij are the stiffness matrix’s coefficients and are defined as follows:
S 11 = α 2 B 11 + β 2 A 66 S 12 = α β ( A 12 + A 66 ) S 13 = α 3 B 11 α β 2 ( B 12 + 2 B 66 ) S 14 = α ( k 1 B 11 s + k 2 B 12 s ) + α β 2 B 66 s ( k 1 A + k 2 B ) S 15 = α   X 13 S 22 = α 2 A 66 + β 2 A 22 S 23 = β 3 B 22 α 2 β ( B 12 + 2 B 66 ) S 24 = β ( k 1 B 12 s + k 2 B 22 s ) + α 2 β ( k 1 A + k 2 B ) B 66 s S 25 = β   X 23 S 33 = α 4 D 11 + β 4 D 22 + 2 α 2 β 2 ( D 12 + 2 D 66 ) + K w + K s ( α 2 + β 2 ) S 34 = α 2 k 1 D 11 s + ( k 2 α 2 + k 1 β 2 ) D 12 s + β 2 k 2 D 22 s 2 α 2 β 2 ( k 1 A + k 2 B ) D 66 s S 35 = α 2   Y 13 + β 2   Y 23 S 44 = k 1 2 H 11 s + k 2 2 H 22 s + 2 k 1 k 2 H 12 s + α 2 β 2 ( k 1 A + k 2 B ) 2 H 66 s + α 2 ( k 1 A ) 2 A 55 s + β 2 ( k 2 B ) 2 A 44 s S 45 = k 1 Y 13 s + k 2 Y 23 s + α 2   k 1 A A 55 s + β 2   k 2 B A 44 s S 55 = α 2   A 55 s + β 2   A 44 s + Z 33

4. Results and Discussion

4.1. Bending Analysis of Simply Supported FG Plates

In this section, different numerical examples are presented and compared to the results of various 2D, 3D, and quasi-3D shear deformation theories to check the accuracy of the present theory in investigating the bending behaviour of simply supported FG plates. The mechanical characteristics of the metal and ceramics used in this study are given in Table 1. In the calculations, both homogeneous isotropic plates and FGPs are studied. Additionally, an analysis of the parameters was performed. For convenience, the following non-dimensional displacements and stresses are used in presenting the numerical results in graphical and tabular forms:
w ¯ = 10 h 3 E c a 4 q 0 w ( a 2 , b 2 , z ) , σ ¯ x = h a q 0 σ x ( a 2 , b 2 , z ) , σ ¯ y = h a q 0 σ y ( a 2 , b 2 , z ) , τ ¯ x y = h a q 0 τ x y ( 0 , 0 , z ) ,   τ ¯ x z = h a q 0 τ x z ( 0 , b 2 , z ) ,   τ ¯ y z = h a q 0 τ y z ( a 2 , 0 , z ) u ˜ = E c h q 0 u   ( 0 , b 2 , z ) , w ˜ = E c h q 0 w ( a 2 , b 2 , z ) , σ ˜ x = σ x ( a 2 , b 2 , z ) q 0 , σ ˜ z = σ z ( a 2 , b 2 , z ) q 0 , τ ˜ x y = τ x y ( 0 , 0 , z ) q 0 ,   τ ˜ x z = τ x z ( 0 , b 2 , z ) q 0 u ^ = G 1 h q 0 u   ( 0 , b 2 , z ) , w ^ = G 1 h q 0 w ( a 2 , b 2 , z ) ,

4.1.1. Functionally Graded Plates (P-FGM)

In this section, non-dimensional displacements and stresses of an A l / A l 2 O 3 FG square plate under uniformly and sinusoidal distributed loads for various values of the power law index are calculated and presented in Table 2 and Table 3.
In Table 2, the computed results of non-dimensional deflection and stress components of moderately thick square FG plate under uniform load are compared to those given by quasi-3D and 2D shear deformation theories by Akavci and Tanrikulu [31] and Younsi et al. [39]. From this table, it can be observed that the results of the proposed 2D theory agree well with those of Akavci and Tanrikulu [31] and Younsi et al. [39] in all cases. Additionally, by comparing the obtained values of the present quasi-3D theory with those given by the other quasi-3D theories, an excellent correlation can be noted.
Table 3 presents the non-dimensional normal stress ( σ ¯ x ) and transverse displacement ( w ¯ ) of square thin and thick plates subjected to sinusoidal load for three different power law indices (k). The determined results are compared with those generated by Carerra et al. [21] based on fourth-order variations of both in-plane and transverse displacements across the thickness, quasi-3D sinusoidal and hyperbolic shear deformation theories of Neves et al. [22,23], quasi-3D shear deformation theories of Hebali et al. [2], and 2D and quasi-3D hyperbolic shear deformation theories of Akavci and Tanrikulu [31] and Younsi et al. [39]. A good agreement between the results is found for both thin and thick plates. Additionally, the present model provides a good prediction of both displacement and stress even in thick FG plates, where the stretching effects are more pronounced.
Moreover, it should be noticed that the number of unknown variables used in the present formulation is just four in the 2D theory and five in the quasi-3D model, while nine unknown functions were used in hyperbolic shear deformation theories of Neves et al. [22,23]. It can be concluded that the present theory is not only accurate but also comparatively simple in predicting the bending response of simply supported FGM plates.
The stresses and displacements variations through the thickness of Al/Al2O3 FG square plate subjected to sinusoidal load are plotted in Figure 2. The results are shown as compared to the quasi-3D shear deformation theory of Akavci and Tanrikulu [31] for different values of the power law index (k). From this figure, excellent agreement between the obtained results and those computed by Akavci and Tanrikulu can be observed [31]. It is important to note that the in-plane stresses ( σ ¯ x ) and ( τ ¯ x y ) exhibit linear variation through the thickness of homogeneous plate, while variation is parabolic for FG plates. The figure also shows that the deflection ( w ¯ ) and in-plane stresses— ( σ ¯ x ) and ( σ ¯ z ) , respectively—increase and that the shear stresses ( τ ¯ x y ) and ( τ ¯ x z ) decrease with the increasing value of material index (k).
In the middle of the FG plate, the axial stresses are negative in the lower half of the metal for all values of the power index k and then become positive in the upper part starting from the geometric ratio z/h > 0.17. Beyond this value, the stresses tend to increase with the increase in the power index, and the most important stress is found for k = 10. An inverse behaviour is observed for the tangential stresses in the plane; they are positive and are obtained in the lower face of the plate, whereas the minimal tangential stresses are obtained on the upper face of the plate FG. The importance of these stresses depends on the quality of the air of FGM materials.
The maximum transverse shear stresses are located in the median plane for homogeneous plates, and they tend to shift slightly towards the upper surface with respect to the median plane with the importance of k. The maximum transverse shear stresses will possess an asymmetric characteristic through the thickness of the FG plate. In this case and for an index k = 1 the maximum axial and shear stresses in the plan are found respectively 12 and 6 times more important than the maximum transverse shear stress. Regarding the deflection ( w ¯ ) , it is found to be almost stable through the thickness of the plate and increases significantly proportionally with the index k.
Figure 3 presents a 3D interaction diagram of the power-law index (k), side-to-thickness ratio (a/h), and centre deflection ( w ¯ ) using the proposed 2D and quasi-3D theories. It can be seen from this figure that the centre deflection increases with the increase of the power law index and decreases with the increasing of the thickness ratio. It is due to the fact that a higher value of (k) corresponds to lower value of volume fraction of the ceramic phase and thus leads to the decrease of the value of the elasticity modulus, which makes the plate softer. The figure shows also that the centre deflections computed from the present 2D theory, which neglects the thickness stretching effect, are higher than those calculated from quasi-3D theory.

4.1.2. Exponentially Graded Plates (E−FGM)

Computations in this section are carried out for a simply supported E-FGM plate. The material properties of the E-FGM plate are defined by the exponential function given in Equation (2). The non-dimensional stress and displacements of the E-FGM plate are calculated and compared with the results of different HSDTs for different loadings.
The centre deflections ( w ¯ ) in-plane and transverse shear stresses ( σ ¯ x , τ ¯ x y ) of Al/Al2O3 plates under sinusoidal loads are calculated for different values of aspect ratio (b/a), thickness ratio (a/h), and exponent values (k) in Table 4, Table 5, Table 6 and Table 7. The central deflections of the very thick E-FGM plates are analysed in Table 4. The obtained predictions are compared with the quasi-3D sinusoidal and exact 3D elasticity theories of Zenkour [19], 2D and quasi-3D trigonometric models of Mantari and Soares [26], and the quasi-3D and 2D shear deformation theories by Akavci and Tanrikulu [31], Younsi et al. [39]. Since the proposed and other quasi-3D models include the thickness-stretching influence, they lead to results close to each other, whereas 2D HSDTs overestimate the deflections due to omitting the thickness-stretching influence. In Table 5 and Table 6, the computed non-dimensional stresses are provided as compared with those given by Younsi et al. [39], Akavci and Tanrikulu [31] based on quasi-3D and 2D hyperbolic theories and Mantari and Soares [26] using 2D and quasi-3D trigonometric theories. It can be seen from the table that an excellent agreement is achieved between the results of present theory and those of other theories. Table 4, Table 5, Table 6 and Table 7 also demonstrate that deflection ( w ¯ ) and transverse shear stress ( τ ¯ x z ) decrease and that in-plane stress ( σ ¯ x ) increases with an increase in the exponent (k).
In Figure 4, the variations of non-dimensional displacements and stresses according to the thickness of an E-FGM plate subjected to sinusoidal loading for different Ec/Em ratios are displayed using the present theory, including the thickness-stretching effect and compared to the quasi-3D theory of Younsi et al. [39]. A very good accuracy between the solutions is observed. It can be seen also from these results that the non-dimensional displacements increase with the increasing of Ec/Em. In addition, it can be deduced that Ec/Em ratios affect considerably the non-dimensional stresses.
It should be noted that the normal stress ( σ ¯ x ) increases in importance in the two extreme surfaces, whether in the ceramic or in the metal depending on the couple Ec/Em When Ec/Em = 0.1, the compressive stress in the metal is on the order of 4.44 times greater than that of ceramic traction. An opposite behaviour occurs in the case that Ec/Em = 10, in which case the compressive stress in the metal is 4.44 times less than that of the traction in the upper layer of the ceramic. When Ec/Em = 1, the normal stress is asymmetrical on both sides of the middle layer of the plate. The maximum shear stresses ( τ ˜ x y ) are in the external layers and they are half in comparison with the normal stresses ( σ ¯ x ) It is also noted that the maximum shearing stresses through the thickness ( τ ˜ x z ) are four times less than the shear stresses ( τ ˜ x y ) .

4.1.3. Comparative Study

In this section, parametric studies have been presented to evaluate the effect of the power law index (k) and side-to-thickness ratio (a/h) on the bending of functionally graded plates using three rules of mixture (P-FGM, E-FGM, Mori–Tanaka models). Figure 5 illustrates the variation of non-dimensional deflection ( w ¯ ) with respect to exponent index (k) and slenderness ratio (a/h) for P-FGM, E-FGM, Mori–Tanaka homogenization scheme. It can be noted from this figure, that when the exponent index increases, deflections of plates using E-FGM decrease whereas the deflections increase when P-FGM and Mori–Tanaka models are used. Such behaviour is due to the fact that the increase in the power law index will increase the stiffness of the EGM plate and thus lead to a reduction of the transverse displacement.
Figure 6 is devoted to present the influence of side-to-thickness ratio on the non-dimensional deflection of a simply supported plate with different models. As can be observed from this figure, an increase in side-to-thickness ratio causes a reduction in the magnitude of deflection. This means that the effect of shear deformation is significant when plates are thick (a/h ≤ 5) and negligible for thin plates. Furthermore, it is apparent that plates made with exponential function have the lowest deflection than the other functions.
An additional comparative study for evaluating the dimensionless centre deflections of the Mori–Tanaka homogenization scheme and P-FGM and E-FGM plates subjected to sinusoidal load is also carried out. In Figure 7, it is shown that the deflection at the centre of plates with Mori–Tanaka homogenization model is larger than those of P-FGM and E-FGM plates. The E-FGM plate which has the smallest deflection is stiffer than the other FGM plates.

4.2. Analysis of FG Plates on Elastic Foundation

This section aims to demonstrate the accuracy of the proposed new models in predicting the static response of FG plates resting on elastic foundation. The following relations of non-dimensional displacements, stresses and foundation parameters in the bending problem are used:
u = 100 D c q 0 a 4 u    ( 0 , b 2 , h 2 ) , v = 100 D c q 0 a 4 v    ( b 2 ,   0 ,   h 2 ) , w = 100 D c q 0 a 4 w   ( a 2 , b 2 ) , σ x = h 2 q 0 a 2 σ x ( a 2 , b 2 , h 2 ) σ y = h 2 q 0 a 2 σ y ( a 2 , b 2 , h 2 ) , σ x y = h 2 q 0 a 2 σ x y ( 0 ,   0 ,   h 2 ) K 0 = K w a 4 E 0 h 3 , J 0 = K s a 2 ν   E 0 h 3 , E 0 = 1 G P a K ¯ w = K w a 4 D c , K ¯ s = K s a 2 D c , D c = E c h 3 12 ( 1 υ 2 )
In the first example, isotropic square plates on elastic foundation are analysed to verify the accuracy of the proposed theories. The dimensionless deflection of a homogeneous square plate subjected to uniform loads is presented in Table 8 for two values of the side-to-thickness ratio (a/h) and different values of the foundation parameters ( K ¯ w , K ¯ s ) . The obtained results are compared to those given by Thai and Choi [29] using the differential quadrature method and those of Al Khateeb and Zenkour [28] based on refined shear deformation plate theory. This table proves that the computed results are in excellent agreement with those reported by the other theories of Thai and Choi [29] and Al Khateeb and Zenkour [28] for all values of side-to-thickness ratio and foundation parameters.
The next example is implemented for Al/Al2O3 moderately thick rectangular plates on elastic foundation. Table 9 presents a comparison of non-dimensional displacements and stresses of FG plate for various exponent values (k) and foundation parameters (K0, J0). derived from the present theory, the sinusoidal shear deformation theory of Zenkour [20], and the 2D zeroth-order shear deformation theory of Thai and Choi [29]. It is clear that the present results are found to be closer with other theories for the different parameters. It should be noted that the present theory is not only efficient but more accurate in predicting the bending behaviour of FG plates resting on elastic foundation. The table exhibits that axial displacements, deflection, and stresses decrease with the increase in foundation parameters. In addition, it is apparent that the computed values from the present quasi-3D theory which takes into account the thickness stretching effect are smaller than those calculated from 2D theory.
To check the effect of Winkler’s and Pasternak’s foundation parameters ( K ¯ w , K ¯ s ) on the bending behaviour of FG plates subjected to sinusoidal loads based on the proposed 2D and quasi-3D theories, Figure 8 depicts the variation in non-dimensional deflection w ¯ of Al/Al2O3 FG square plates versus Winkler ( K ¯ w ) and Pasternak ( K ¯ s ) parameters. As can be seen from this figure, the dimensionless deflection diminishes when the foundation parameters increase. Compared to the Winkler parameter ( K ¯ w ) , the Pasternak foundation parameter ( K ¯ s ) has a more significant effect on decreasing the dimensionless deflection. Additionally, the deflection computed by present quasi-3D theory is lower than that obtained from 2D theory, which means the thickness-stretching effect is more significant when the thickness of the plate is higher, but it is negligible with reducing the thickness of the plate. This indicates that the stretching effect is pronounced and must be taken into account in the modelling of thick plate.

5. Conclusions

The aim of this work is to develop a mathematical approach in order to study the mechanical behaviour of functionally graded plates on elastic foundation by applying the 2D and quasi-3D shear strain theories proposed by Zaoui et al. (2019) regarding free vibration analysis. The theory, in which only five unknowns are used, takes into account a parabolic distribution of the transverse shear stress using a new shear shape function which satisfies the boundary conditions with no tension or compression on the bottom and top surfaces of the plate. Moreover, the shear correction factor has not been included. The principle of virtual work was applied to determine the equilibrium equations. The analytical solutions of simply supported plates are obtained using the double Fourier series and are compared to previously published solutions in order to demonstrate the accuracy of the proposed theory. The effects of a large number of parameters, such as the transverse normal strain, the form ratio of the plate, the length/thickness ratio, the volume fraction distributions, the various mixing laws, and the boundary conditions on the deflection, axial displacements and the various stresses are examined in detail. From the results analysed above, the following main points can be drawn:
  • Through all the comparative investigations, it can be seen that the proposed theory proves good agreement with the results of other 2D and quasi-3D HSDTs.
  • The proposed quasi-3D HSDT contains five unknowns but provides results comparable to those generated by other published quasi-3D theories with greater numbers of unknowns (e.g., quasi-3D theories of Neves et al. [22,23] with nine unknowns and Akavci and Tanrikulu [31] with six unknowns).
  • The power law index has the effect of increasing deflections of plate. The influence of the value of material index more than 5 is negligible.
  • The small difference between the present 2D and quasi-3D shear deformation results is due to neglecting the thickness stretching effect.
  • Results demonstrate that the plate becomes stiffer when the effects of normal deformations are considered, consequently leading to a decrease in deflection and an increase in stresses.
  • The thickness stretching effect is more pronounced for thick plates and it needs to be taken into account in the modelling.
  • The decreasing effect of foundation parameters on the deflection leads to an increase in the FG plate stiffness.
Consequently, the present method yields more accurate results for plate bending problems, as they are closer when compared to those reported in other literatures using numerical methods. Finally, it can be concluded that the formulation developed herein can serve as a reference for future research and could be extended to the analysis of classical and thermal buckling of FG plates by introducing various combined boundary conditions. Additionally, the present method can be combined with the nonlocal continuum mechanics to model the nano and micro material mechanical behaviours.

Author Contributions

Methodology and supervision, D.O.; project administration, writing—review and editing, B.A.; software and writing—original draft preparation, F.Z.Z.; visualization, M.B.; validation, J.A.V.O. and A.A.A.A.-N.; investigation, E.R.L. and M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Deanship of Scientific Research at the University of Hail, Saudi Arabia, under the project number RG-20 098.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their deepest gratitude to the Deanship of Scientific Research and to the College of Engineering at the University of Hail, for providing necessary support for conducting this research, through the project number RG-20 098.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

ParametersDescriptionUnitsParametersDescriptionUnits
aLengthmpPower law index/
bWidthmσijStress tensorN m−2
hThicknessmεijStrain tensor
EYoung’s modulusGPaQijEngineering constantsGPa
vPoisson ratio/u, v, wDisplacements in x, y, z axis.m
ρMass DensityKg/m3qMechanical loadN
PMaterial property/f(z)Shape function/
VVolume fraction/m, nmodes/

References

  1. Jha, D.K.; Kant, T.; Singh, R.K. Higher Order Shear and Normal Deformation Theory for Natural Frequency of Functionally Graded Rectangular Plates. Nucl. Eng. Des. 2012, 250, 8–13. [Google Scholar] [CrossRef]
  2. Hebali, H.; Tounsi, A.; Houari, M.S.A.; Bessaim, A.; Bedia, E.A.A. New Quasi-3D Hyperbolic Shear Deformation Theory for the Static and Free Vibration Analysis of Functionally Graded Plates. J. Eng. Mech. 2014, 140, 374–383. [Google Scholar] [CrossRef]
  3. Nguyen, T.-K. A Higher-Order Hyperbolic Shear Deformation Plate Model for Analysis of Functionally Graded Materials. Int. J. Mech. Mater. Des. 2015, 11, 203–219. [Google Scholar] [CrossRef]
  4. Abdulrazzaq, M.A.; Fenjan, R.M.; Ahmed, R.A.; Faleh, N.M. Thermal Buckling of Nonlocal Clamped Exponentially Graded Plate According to a Secant Function Based Refined Theory. Steel Compos. Struct. 2020, 35, 147–157. [Google Scholar]
  5. Alijani, F.; Amabili, M. Effect of Thickness Deformation on Large-Amplitude Vibrations of Functionally Graded Rectangular Plates. Compos. Struct. 2014, 113, 89–107. [Google Scholar] [CrossRef]
  6. Arefi, M.; Allam, M.N.M. Nonlinear Responses of an Arbitrary FGP Circular Plate Resting on the Winkler-Pasternak Foundation. Smart Struct. Syst. 2015, 16, 81–100. [Google Scholar] [CrossRef]
  7. Amir, S.; Arshid, E.; Rasti-Alhosseini, S.M.A.; Loghman, A. Quasi-3D Tangential Shear Deformation Theory for Size-Dependent Free Vibration Analysis of Three-Layered FG Porous Micro Rectangular Plate Integrated by Nano-Composite Faces in Hygrothermal Environment. J. Therm. Stresses 2020, 43, 133–156. [Google Scholar] [CrossRef]
  8. Azadi, M. Free and Forced Vibration Analysis of FG Beam Considering Temperature Dependency of Material Properties. J. Mech. Sci. Technol. 2011, 25, 69–80. [Google Scholar] [CrossRef]
  9. Zaoui, F.Z.; Ouinas, D.; Tounsi, A.; Viña Olay, J.A.; Achour, B.; Touahmia, M. Fundamental Frequency Analysis of Functionally Graded Plates with Temperature-Dependent Properties Based on Improved Exponential-Trigonometric Two-Dimensional Higher Shear Deformation Theory. Arch. Appl. Mech. 2021, 91, 859–881. [Google Scholar] [CrossRef]
  10. Hosseini-Hashemi, S.; Rokni Damavandi Taher, H.; Akhavan, H.; Omidi, M. Free Vibration of Functionally Graded Rectangular Plates Using First-Order Shear Deformation Plate Theory. Appl. Math. Model. 2010, 34, 1276–1291. [Google Scholar] [CrossRef]
  11. Golmakani, M.E.; Alamatian, J. Large Deflection Analysis of Shear Deformable Radially Functionally Graded Sector Plates on Two-Parameter Elastic Foundations. Eur. J. Mech.-A/Solids 2013, 42, 251–265. [Google Scholar] [CrossRef]
  12. Duc, N.D.; Quan, T.Q. Nonlinear Dynamic Analysis of Imperfect Functionally Graded Material Double Curved Thin Shallow Shells with Temperature-Dependent Properties on Elastic Foundation. J. Vib. Control 2015, 21, 1340–1362. [Google Scholar] [CrossRef]
  13. Khalili, A.; Vosoughi, A.R. An Approach for the Pasternak Elastic Foundation Parameters Estimation of Beams Using Simulated Frequencies. Inverse Probl. Sci. Eng. 2018, 26, 1079–1093. [Google Scholar] [CrossRef]
  14. Meksi, A.; Benyoucef, S.; Houari, M.S.A.; Tounsi, A. A Simple Shear Deformation Theory Based on Neutral Surface Position for Functionally Graded Plates Resting on Pasternak Elastic Foundations. Struct. Eng. Mech. Int. J. 2015, 53, 1215–1240. [Google Scholar] [CrossRef]
  15. Shahbaztabar, A.; Ranji, A.R. Effects of In-Plane Loads on Free Vibration of Symmetrically Cross-Ply Laminated Plates Resting on Pasternak Foundation and Coupled with Fluid. Ocean Eng. 2016, 115, 196–209. [Google Scholar] [CrossRef]
  16. Radaković, A.; Čukanović, D.; Bogdanović, G.; Blagojević, M.; Stojanović, B.; Dragović, D.; Manić, N. Thermal Buckling and Free Vibration Analysis of Functionally Graded Plate Resting on an Elastic Foundation According to High Order Shear Deformation Theory Based on New Shape Function. Appl. Sci. 2020, 10, 4190. [Google Scholar] [CrossRef]
  17. Zenkour, A.M.; El-Shahrany, H.D. Hygrothermal Forced Vibration of a Viscoelastic Laminated Plate with Magnetostrictive Actuators Resting on Viscoelastic Foundations. Int. J. Mech. Mater. Des. 2021, 17, 301–320. [Google Scholar] [CrossRef]
  18. Zaoui, F.Z.; Ouinas, D.; Achour, B.; Tounsi, A.; Latifee, E.R.; Al-Naghi, A.A.A. A Hyperbolic Shear Deformation Theory for Natural Frequencies Study of Functionally Graded Plates on Elastic Supports. J. Compos. Sci. 2022, 6, 285. [Google Scholar] [CrossRef]
  19. Zenkour, A.M. Benchmark Trigonometric and 3-D Elasticity Solutions for an Exponentially Graded Thick Rectangular Plate. Arch. Appl. Mech. 2007, 77, 197–214. [Google Scholar] [CrossRef]
  20. Zenkour, A.M. The Refined Sinusoidal Theory for FGM Plates on Elastic Foundations. Int. J. Mech. Sci. 2009, 51, 869–880. [Google Scholar] [CrossRef]
  21. Carrera, E.; Brischetto, S.; Cinefra, M.; Soave, M. Effects of Thickness Stretching in Functionally Graded Plates and Shells. Compos. Part B Eng. 2011, 42, 123–133. [Google Scholar] [CrossRef] [Green Version]
  22. Neves, A.M.A.; Ferreira, A.J.M.; Carrera, E.; Roque, C.M.C.; Cinefra, M.; Jorge, R.M.N.; Soares, C.M.M. A Quasi-3D Sinusoidal Shear Deformation Theory for the Static and Free Vibration Analysis of Functionally Graded Plates. Compos. Part B Eng. 2012, 43, 711–725. [Google Scholar] [CrossRef]
  23. Neves, A.M.A.; Ferreira, A.J.M.; Carrera, E.; Cinefra, M.; Roque, C.M.C.; Jorge, R.M.N.; Soares, C.M.M. A Quasi-3D Hyperbolic Shear Deformation Theory for the Static and Free Vibration Analysis of Functionally Graded Plates. Compos. Struct. 2012, 94, 1814–1825. [Google Scholar] [CrossRef]
  24. Thai, H.-T.; Park, M.; Choi, D.-H. A Simple Refined Theory for Bending, Buckling, and Vibration of Thick Plates Resting on Elastic Foundation. Int. J. Mech. Sci. 2013, 73, 40–52. [Google Scholar] [CrossRef]
  25. Ouinas, D.; Achour, B. Buckling Analysis of Laminated Composite Plates [(θ/−θ)] Containing an Elliptical Notch. Compos. Part B Eng. 2013, 55, 575–579. [Google Scholar] [CrossRef]
  26. Mantari, J.L.; Guedes Soares, C. A Novel Higher-Order Shear Deformation Theory with Stretching Effect for Functionally Graded Plates. Compos. Part B Eng. 2013, 45, 268–281. [Google Scholar] [CrossRef]
  27. Zhang, H.; Jiang, J.; Zhang, Z. Three-Dimensional Elasticity Solutions for Bending of Generally Supported Thick Functionally Graded Plates. Appl. Math. Mech.-Engl. Ed. 2014, 35, 1467–1478. [Google Scholar] [CrossRef]
  28. Al Khateeb, S.A.; Zenkour, A.M. A Refined Four-Unknown Plate Theory for Advanced Plates Resting on Elastic Foundations in Hygrothermal Environment. Compos. Struct. 2014, 111, 240–248. [Google Scholar] [CrossRef]
  29. Thai, H.-T.; Choi, D.-H. Zeroth-Order Shear Deformation Theory for Functionally Graded Plates Resting on Elastic Foundation. Int. J. Mech. Sci. 2014, 78, 35–43. [Google Scholar] [CrossRef]
  30. 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. Compos. Struct. 2015, 122, 330–342. [Google Scholar] [CrossRef]
  31. Akavci, S.S.; Tanrikulu, A.H. Static and Free Vibration Analysis of Functionally Graded Plates Based on a New Quasi-3D and 2D Shear Deformation Theories. Compos. Part B Eng. 2015, 83, 203–215. [Google Scholar] [CrossRef]
  32. Mantari, J.L.; Granados, E.V. An Original FSDT to Study Advanced Composites on Elastic Foundation. Thin-Walled Struct. 2016, 107, 80–89. [Google Scholar] [CrossRef]
  33. Houari, M.S.A.; Tounsi, A.; Bessaim, A.; Mahmoud, S.R. A New Simple Three-Unknown Sinusoidal Shear Deformation Theory for Functionally Graded Plates. Steel Compos. Struct. 2016, 22, 257–276. [Google Scholar] [CrossRef]
  34. Aldousari, S.M. Bending Analysis of Different Material Distributions of Functionally Graded Beam. Appl. Phys. A 2017, 123, 296. [Google Scholar] [CrossRef]
  35. Meftah, A.; Bakora, A.; Zaoui, F.Z.; Tounsi, A.; Bedia, E.A.A. A Non-Polynomial Four Variable Refined Plate Theory for Free Vibration of Functionally Graded Thick Rectangular Plates on Elastic Foundation. Steel Compos. Struct. 2017, 23, 317–330. [Google Scholar] [CrossRef]
  36. Zaoui, F.Z.; Tounsi, A.; Ouinas, D. Free Vibration of Functionally Graded Plates Resting on Elastic Foundations Based on Quasi-3D Hybrid-Type Higher Order Shear Deformation Theory. Smart Struct. Syst. 2017, 20, 509–524. [Google Scholar]
  37. Guerroudj, H.Z.; Yeghnem, R.; Kaci, A.; Zaoui, F.Z.; Benyoucef, S.; Tounsi, A. Eigenfrequencies of Advanced Composite Plates Using an Efficient Hybrid Quasi-3D Shear Deformation Theory. Smart Struct. Syst. 2018, 22, 121–132. [Google Scholar]
  38. Amar, L.H.H.; Kaci, A.; Yeghnem, R.; Tounsi, A. A New Four-Unknown Refined Theory Based on Modified Couple Stress Theory for Size-Dependent Bending and Vibration Analysis of Functionally Graded Micro-Plate. Steel Compos. Struct. 2018, 26, 89–102. [Google Scholar] [CrossRef]
  39. Younsi, A.; Tounsi, A.; Zaoui, F.Z.; Bousahla, A.A.; Mahmoud, S.R. Novel Quasi-3D and 2D Shear Deformation Theories for Bending and Free Vibration Analysis of FGM Plates. Geomech. Eng. 2018, 14, 519–532. [Google Scholar]
  40. Belkhodja, Y.; Ouinas, D.; Zaoui, F.Z.; Fekirini, H. An Exponential-Trigonometric Higher Order Shear Deformation Theory (HSDT) for Bending, Free Vibration, and Buckling Analysis of Functionally Graded Materials (FGMs) Plates. Adv. Compos. Lett. 2020, 29, 096369351987573. [Google Scholar] [CrossRef]
  41. Wang, Y.; Zeng, R.; Safarpour, M. Vibration Analysis of FG-GPLRC Annular Plate in a Thermal Environment. Mech. Based Des. Struct. Mach. 2022, 50, 352–370. [Google Scholar] [CrossRef]
  42. Shah, N.A.; Wakif, A.; El-Zahar, E.R.; Ahmad, S.; Yook, S.-J. Numerical Simulation of a Thermally Enhanced EMHD Flow of a Heterogeneous Micropolar Mixture Comprising (60%)-Ethylene Glycol (EG), (40%)-Water (W), and Copper Oxide Nanomaterials (CuO). Case Stud. Therm. Eng. 2022, 35, 102046. [Google Scholar] [CrossRef]
  43. Sajjan, K.; Shah, N.A.; Ahammad, N.A.; Raju, C.S.K.; Kumar, M.D.; Weera, W.; Sajjan, K.; Shah, N.A.; Ahammad, N.A.; Raju, C.S.K.; et al. Nonlinear Boussinesq and Rosseland Approximations on 3D Flow in an Interruption of Ternary Nanoparticles with Various Shapes of Densities and Conductivity Properties. AIMS Math. 2022, 7, 18416–18449. [Google Scholar] [CrossRef]
  44. Raza, Q.; Qureshi, M.Z.A.; Khan, B.A.; Kadhim Hussein, A.; Ali, B.; Shah, N.A.; Chung, J.D. Insight into Dynamic of Mono and Hybrid Nanofluids Subject to Binary Chemical Reaction, Activation Energy, and Magnetic Field through the Porous Surfaces. Mathematics 2022, 10, 3013. [Google Scholar] [CrossRef]
  45. Sabu, A.S.; Wakif, A.; Areekara, S.; Mathew, A.; Shah, N.A. Significance of Nanoparticles’ Shape and Thermo-Hydrodynamic Slip Constraints on MHD Alumina-Water Nanoliquid Flows over a Rotating Heated Disk: The Passive Control Approach. Int. Commun. Heat Mass Transf. 2021, 129, 105711. [Google Scholar] [CrossRef]
  46. Qureshi, M.Z.A.; Faisal, M.; Raza, Q.; Ali, B.; Botmart, T.; Shah, N.A.; Qureshi, M.Z.A.; Faisal, M.; Raza, Q.; Ali, B.; et al. Morphological Nanolayer Impact on Hybrid Nanofluids Flow Due to Dispersion of Polymer/CNT Matrix Nanocomposite Material. AIMS Math. 2023, 8, 633–656. [Google Scholar] [CrossRef]
  47. Mahsud, Y.; Shah, N.A.; Vieru, D. Influence of Time-Fractional Derivatives on the Boundary Layer Flow of Maxwell Fluids. Chin. J. Phys. 2017, 55, 1340–1351. [Google Scholar] [CrossRef]
  48. Elnaqeeb, T.; Shah, N.A.; Mirza, I.A. Natural Convection Flows of Carbon Nanotubes Nanofluids with Prabhakar-like Thermal Transport. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  49. Vieru, D.; Fetecau, C.; Shah, N.A.; Yook, S.-J. Unsteady Natural Convection Flow Due to Fractional Thermal Transport and Symmetric Heat Source/Sink. Alex. Eng. J. 2022; in press. [Google Scholar] [CrossRef]
  50. Saffari, P.R.; Fakhraie, M.; Roudbari, M.A. Nonlinear Vibration of Fluid Conveying Cantilever Nanotube Resting on Visco-Pasternak Foundation Using Non-Local Strain Gradient Theory. Micro Nano Lett. 2020, 15, 181–186. [Google Scholar] [CrossRef]
  51. Alazwari, M.A.; Zenkour, A.M. A Quasi-3D Refined Theory for the Vibration of Functionally Graded Plates Resting on Visco-Winkler-Pasternak Foundations. Mathematics 2022, 10, 716. [Google Scholar] [CrossRef]
  52. Zaoui, F.Z.; Ouinas, D.; Tounsi, A. New 2D and Quasi-3D Shear Deformation Theories for Free Vibration of Functionally Graded Plates on Elastic Foundations. Compos. Part B Eng. 2019, 159, 231–247. [Google Scholar] [CrossRef]
  53. Zaoui, F.Z.; Hanifi, H.A.L.; Younsi, A.; Meradjah, M.; Tounsi, A.; Ouinas, D. Free Vibration Analysis of Functionally Graded Beams Using a Higher-Order Shear Deformation Theory. Math. Model. Eng. Probl. 2017, 4, 7–12. [Google Scholar] [CrossRef]
  54. Meradjah, M.; Bouakkaz, K.; Zaoui, F.Z.; Tounsi, A. A Refined Quasi-3D Hybrid-Type Higher Order Shear Deformation Theory for Bending and Free Vibration Analysis of Advanced Composites Beams. Wind. Struct. 2018, 27, 269–282. [Google Scholar]
  55. Zaoui, F.Z.; Tounsi, A.; Ouinas, D.; Viña Olay, J.A. A Refined HSDT for Bending and Dynamic Analysis of FGM Plates. Struct. Eng. Mech. Int’l J. 2020, 74, 105–119. [Google Scholar]
  56. Zenkour, A.M.; Hafed, Z.S.; Radwan, A.F. Bending Analysis of Functionally Graded Nanoscale Plates by Using Nonlocal Mixed Variational Formula. Mathematics 2020, 8, 1162. [Google Scholar] [CrossRef]
  57. Rachid, A.; Ouinas, D.; Lousdad, A.; Zaoui, F.Z.; Achour, B.; Gasmi, H.; Butt, T.A.; Tounsi, A. Mechanical Behavior and Free Vibration Analysis of FG Doubly Curved Shells on Elastic Foundation via a New Modified Displacements Field Model of 2D and Quasi-3D HSDTs. Thin-Walled Struct. 2022, 172, 108783. [Google Scholar] [CrossRef]
  58. Vaghefi, R.; Baradaran, G.H.; Koohkan, H. Three-Dimensional Static Analysis of Thick Functionally Graded Plates by Using Meshless Local Petrov–Galerkin (MLPG) Method. Eng. Anal. Bound. Elem. 2010, 34, 564–573. [Google Scholar] [CrossRef]
Figure 1. Geometry and coordinates of the considered FGM plate on elastic foundation.
Figure 1. Geometry and coordinates of the considered FGM plate on elastic foundation.
Mathematics 10 04764 g001
Figure 2. Distributions of the non-dimensional displacements ( u ¯ , w ¯ ) and stresses ( σ ¯ x , σ ¯ z , τ ¯ x y , τ ¯ x z ) through the thickness of square P−GM plate (a/h = 10) [31].
Figure 2. Distributions of the non-dimensional displacements ( u ¯ , w ¯ ) and stresses ( σ ¯ x , σ ¯ z , τ ¯ x y , τ ¯ x z ) through the thickness of square P−GM plate (a/h = 10) [31].
Mathematics 10 04764 g002
Figure 3. Effect of the power law index k and side-to-thickness ratio a/h on the non-dimensional deflection w ¯ of simply supported Al/Al2O3 P−FGM plates.
Figure 3. Effect of the power law index k and side-to-thickness ratio a/h on the non-dimensional deflection w ¯ of simply supported Al/Al2O3 P−FGM plates.
Mathematics 10 04764 g003
Figure 4. Distributions of the non-dimensional displacements ( u ˜ , w ˜ ) and stresses ( σ ˜ x , σ ˜ z , τ ˜ x y , τ ˜ x z ) through the thickness of square E-FGM plate (a/h = 10) [39].
Figure 4. Distributions of the non-dimensional displacements ( u ˜ , w ˜ ) and stresses ( σ ˜ x , σ ˜ z , τ ˜ x y , τ ˜ x z ) through the thickness of square E-FGM plate (a/h = 10) [39].
Mathematics 10 04764 g004
Figure 5. Non-dimensional deflection variation w ¯ of different models of FG plates for various values of power-law index (k) and side-to-thickness ratio (a/h).
Figure 5. Non-dimensional deflection variation w ¯ of different models of FG plates for various values of power-law index (k) and side-to-thickness ratio (a/h).
Mathematics 10 04764 g005
Figure 6. Effect of the side-to-thickness ratio (a/h) on the non-dimensional deflection of different FG plates.
Figure 6. Effect of the side-to-thickness ratio (a/h) on the non-dimensional deflection of different FG plates.
Mathematics 10 04764 g006
Figure 7. Comparison of the deflection of the Mori–Tanaka model and P-FGM and E-FGM plates (a/h = 4, k = 2).
Figure 7. Comparison of the deflection of the Mori–Tanaka model and P-FGM and E-FGM plates (a/h = 4, k = 2).
Mathematics 10 04764 g007
Figure 8. Effect of foundation parameters ( K ¯ w , K ¯ s ) on the non-dimensional deflection w ¯ of FG square plates under sinusoidal load (k =1, a/h = 10).
Figure 8. Effect of foundation parameters ( K ¯ w , K ¯ s ) on the non-dimensional deflection w ¯ of FG square plates under sinusoidal load (k =1, a/h = 10).
Mathematics 10 04764 g008
Table 1. Material properties used in the FG plates [52,55].
Table 1. Material properties used in the FG plates [52,55].
MaterialProperties
E (GPa)Vρ (kg/m3)
Aluminum (Al)700.32702
Alumina (Al2O3)3800.33800
Zirconia (ZrO2)2000.35700
Table 2. The non-dimensional displacement and stress components of an Al/Al2O3 FG square plate subjected to uniformly distributed load (a/h = 10).
Table 2. The non-dimensional displacement and stress components of an Al/Al2O3 FG square plate subjected to uniformly distributed load (a/h = 10).
kTheory ε z w ¯ ( 0 ) σ ¯ x ( h 2 ) σ ¯ y ( h 3 ) τ ¯ x z ( 0 ) τ ¯ y z ( h 6 ) τ ¯ x y ( h 3 )
0Akavci and Tanrikulu [31]= 00.46652.89091.91030.49880.43631.2857
Akavci and Tanrikulu [31]≠ 00.46352.99811.89250.47820.43151.2578
Younsi et al. [39]= 00.46652.89131.91020.50430.43671.2855
Younsi et al. [39]≠ 00.46372.99191.89320.50420.43171.2585
Present study= 00.46652.89121.91020.50430.43691.2856
Present study≠ 00.46253.07291.87560.47610.43071.2548
1Akavci and Tanrikulu [31]= 00.92884.47072.16930.49880.53641.1141
Akavci and Tanrikulu [31]≠ 00.89774.61102.08220.47820.51191.0211
Younsi et al. [39]= 00.92874.47132.16920.50430.53701.1141
Younsi et al. [39]≠ 00.89804.60052.08320.47910.51211.0225
Present study= 00.92874.47132.16920.50420.53721.1141
Present study≠ 00.89614.73792.05780.47610.51141.0206
2Akavci and Tanrikulu [31]= 01.19405.22482.03420.45810.56430.9909
Akavci and Tanrikulu [31]≠ 01.13765.38251.92570.45240.50810.8921
Younsi et al. [39]= 01.19405.22562.03400.46370.56570.9908
Younsi et al. [39]≠ 01.13805.37261.92810.45320.50820.8926
Present study= 01.19405.22552.03400.46360.56580.9908
Present study≠ 01.13525.52321.89720.45050.50740.8902
4Akavci and Tanrikulu [31]= 01.38885.88551.72050.40900.52531.0305
Akavci and Tanrikulu [31]≠ 01.32596.03821.60620.43580.48040.9274
Younsi et al. [39]= 01.38905.88661.72020.41510.52781.0303
Younsi et al. [39]≠ 01.32626.03011.61010.43650.48060.9279
Present study= 01.38895.88651.72020.41490.52791.0303
Present study≠ 01.32376.19201.57440.43410.47970.9256
10Akavci and Tanrikulu [31]= 01.58757.36171.28280.44360.41591.0705
Akavci and Tanrikulu [31]≠ 01.54537.51231.20160.43320.45610.9860
Younsi et al. [39]= 01.58757.36281.28250.44950.41741.0703
Younsi et al. [39]≠ 01.54547.50641.20590.43390.45620.9862
Present study= 01.58757.36281.28250.44950.41761.0703
Present study≠ 01.54367.69141.17240.43140.45540.9852
Table 3. Non-dimensional displacement and stress of an Al/Al2O3 FG square plate subjected to sinusoidal load.
Table 3. Non-dimensional displacement and stress of an Al/Al2O3 FG square plate subjected to sinusoidal load.
kTheory ε z σ ¯ x ( h / 3 ) w ¯ ( 0 )
a/h = 4a/h = 10a/h = 100a/h = 4a/h = 10a/h = 100
1Carrera et al. [21]≠ 00.62211.506414.96900.71710.58750.5625
Neves et al. [22]≠ 00.59251.494514.96900.69970.58450.5624
Neves et al. [23]≠ 00.59101.491714.94400.70200.58680.5648
Hebali et al. [2]≠ 00.59521.495414.96300.69100.56860.5452
Akavci and Tanrikulu [31]= 00.58061.489514.96700.72820.58890.5625
Akavci and Tanrikulu [31]≠ 00.57541.432214.30600.69080.56910.5457
Younsi et al. [39]= 00.58081.489614.96750.72830.58890.5625
Younsi et al. [39]≠ 00.57581.433014.31350.69100.56920.5459
Present study= 00.58031.489414.96750.72800.58890.5625
Present study≠ 00.57051.415714.13300.68960.56800.5447
4Carrera et al. [21]≠ 00.48771.197111.92301.15850.88210.8286
Neves et al. [22]≠ 00.44041.178311.93201.11780.87500.8286
Neves et al. [23]≠ 00.43401.159311.73801.10950.86980.8241
Hebali et al. [2]≠ 00.45071.177911.87101.09640.84130.7926
Akavci and Tanrikulu [31]= 00.44311.178711.92001.16130.88180.8287
Akavci and Tanrikulu [31]≠ 00.42471.101711.08801.09830.84170.7925
Younsi et al. [39]= 00.44371.178911.92091.16090.88170.8287
Younsi et al. [39]≠ 00.42601.104511.11521.09820.84190.7928
Present study= 00.44241.178311.92081.16180.88180.8287
Present study≠ 00.41811.080210.86331.09700.84030.7910
10Carrera et al. [21]≠ 00.39650.89658.90771.37451.00720.9361
Neves et al. [22]≠ 00.32271.178311.93201.34900.87500.8286
Neves et al. [23]≠ 00.31080.84678.60131.33270.98860.9228
Hebali et al. [2]≠ 00.33250.88898.99771.33330.97910.9114
Akavci and Tanrikulu [31]= 00.32420.87788.90591.39171.00890.9362
Akavci and Tanrikulu [31]≠ 00.30950.82298.31851.33520.98180.9141
Younsi et al. [39]= 00.32480.87808.90591.39151.00880.9362
Younsi et al. [39]≠ 00.31090.82598.34731.33530.98190.9141
Present study= 00.32350.87758.90591.39171.00890.9362
Present study≠ 00.30330.80318.11181.33330.98070.9130
Table 4. Non-dimensional deflection w ¯ ( 0 ) = 10 h 3 E 0 a 4 q 0 w ( a 2 , b 2 , 0 ) of E-FGM plates subjected to sinusoidal distributed load (a/h = 2).
Table 4. Non-dimensional deflection w ¯ ( 0 ) = 10 h 3 E 0 a 4 q 0 w ( a 2 , b 2 , 0 ) of E-FGM plates subjected to sinusoidal distributed load (a/h = 2).
b/aTheoryεzk
0.10.30.50.711.5
1Zenkour [19]≠ 00.57690.52470.47660.43240.37260.2890
Zenkour [19]= 00.57300.51800.46780.42210.36110.2771
Mantari and Soares [26] ≠ 00.57780.52240.47170.42560.36480.2793
Mantari and Soares [26]= 00.63620.57510.51940.46870.40170.3079
Akavci and Tanrikulu [31]= 00.63510.57410.51850.46790.40040.3075
Akavci and Tanrikulu [31]≠ 00.57500.51980.46940.42360.36240.2780
Younsi et al. [39]= 00.63550.57450.51890.46830.40070.3077
Younsi et al. [39]≠ 00.57580.52050.47010.42420.36290.2784
Present study= 00.63430.57340.51790.46740.40000.3072
Present study≠ 00.57310.51810.46790.42220.36120.2771
2Zenkour [19]≠ 01.19441.08590.98640.89520.77260.6017
Zenkour [19]= 01.18791.07390.97000.87540.74930.5757
Mantari and Soares [26]≠ 01.19401.07940.97500.87990.75370.5786
Mantari and Soares [26]= 01.27761.15531.04410.94300.80920.6237
Akavci and Tanrikulu [31]= 01.27631.15411.04310.94220.80790.6234
Akavci and Tanrikulu [31]≠ 01.19381.07650.97230.87750.75110.5771
Younsi et al. [39]= 01.27681.15461.04350.94260.80820.6236
Younsi et al. [39]≠ 01.19171.07740.97310.87820.75170.5775
Present study= 01.27531.15321.04230.94150.80740.6231
Present study≠ 01.18801.07400.97010.87550.74940.5758
3Zenkour [19]≠ 01.44291.31161.91121.08110.93330.7275
Zenkour [19]= 01.43541.29771.17221.05790.90560.6961
Mantari and Soares [26]≠ 01.44211.30371.17761.06270.91040.6992
Mantari and Soares [26]= 01.53401.38731.25401.13290.97250.7506
Akavci and Tanrikulu [31]= 01.53271.38611.25301.13200.97120.7503
Akavci and Tanrikulu [31]≠ 01.43861.30051.17481.06020.90760.6976
Younsi et al. [39]= 01.53321.38661.25341.13240.97150.7504
Younsi et al. [39]≠ 01.43961.30151.17561.06100.90820.6981
Present study= 01.53161.38521.25211.13130.97060.7499
Present study≠ 01.43541.29771.17221.05790.90570.6961
Table 5. Non-dimensional stress σ ¯ x ( h 2 ) = h 2 a 2 q 0 σ x ( a 2 , b 2 , h 2 ) of E-FGM plates subjected to sinusoidal distributed load (a/h = 10).
Table 5. Non-dimensional stress σ ¯ x ( h 2 ) = h 2 a 2 q 0 σ x ( a 2 , b 2 , h 2 ) of E-FGM plates subjected to sinusoidal distributed load (a/h = 10).
b/aTheoryεzk
0.10.30.50.711.522.53
1Mantari and Soares [26]≠ 00.21960.23450.25030.26710.29440.34600.40650.47750.5603
Mantari and Soares [26]= 00.20620.22040.23550.25150.27740.32640.38350.45020.5278
Akavci and Tanrikulu [31]= 00.20630.22050.23560.25160.27760.32660.38380.45040.5281
Akavci and Tanrikulu [31]≠ 00.21420.22850.24380.26010.28660.33700.39640.46640.5485
Younsi et al. [39]= 00.20630.22050.23550.25160.27750.32650.38370.45040.5279
Younsi et al. [39]≠ 00.21370.22800.24330.25950.28600.33630.39570.46570.5478
Present study= 00.20630.22050.23560.25170.27760.32660.38380.45050.5282
Present study≠ 00.21950.23440.25020.26700.29430.43590.40640.47730.5602
2Mantari and Soares [26]≠ 00.45520.48670.52000.55540.61260.72010.84490.98981.1580
Mantari and Soares [26]= 00.43500.46490.49660.53030.58500.68810.80850.94901.1125
Akavci and Tanrikulu [31]= 00.43510.46500.49680.53050.58520.68840.80880.94931.1129
Akavci and Tanrikulu [31]≠ 00.44660.47730.50980.54430.60020.70580.82890.97251.1397
Younsi et al. [39]= 00.43510.46500.49670.53050.58510.68830.80870.94921.1128
Younsi et al. [39]≠ 00.44590.47650.50900.54350.59930.70480.82780.82781.1388
Present study= 00.43510.46500.49680.53060.58520.68840.80890.94941.1131
Present study≠ 00.45510.48650.51990.55530.61240.71990.84470.98971.1579
3Mantari and Soares [26]≠ 00.55140.58960.63020.67330.74270.87301.02401.19901.4017
Mantari and Soares [26]= 00.52880.56510.60370.64470.71120.83650.98281.15361.3523
Akavci and Tanrikulu [31]= 00.52900.56530.60390.64490.71140.83680.98321.15401.3528
Akavci and Tanrikulu [31]≠ 00.54180.57910.61870.66080.72890.85701.00611.17971.3813
Younsi et al. [39]= 00.52890.56520.60380.64490.71130.83670.98311.15381.3527
Younsi et al. [39]≠ 00.54100.57830.61790.65990.72790.85591.00501.17861.3803
Present study= 00.52900.56530.60390.64500.71140.83680.98331.15411.3529
Present study≠ 00.55120.58950.63000.67310.74250.87281.02381.19881.4016
Table 6. Non-dimensional stress τ ¯ x z ( 0 ) = h a   q 0 τ x z ( 0 , b 2 , 0 ) of E-FGM plates subjected to sinusoidal distributed load (a/h = 10).
Table 6. Non-dimensional stress τ ¯ x z ( 0 ) = h a   q 0 τ x z ( 0 , b 2 , 0 ) of E-FGM plates subjected to sinusoidal distributed load (a/h = 10).
b/aTheoryεzk
0.10.30.50.711.522.53
1Mantari and Soares [26]≠ 00.24540.24500.24420.24300.24050.23440.22630.21620.2045
Mantari and Soares [26]= 00.23800.23760.23680.23560.23300.22680.21850.20940.1985
Akavci and Tanrikulu [31]= 00.24340.24300.24220.24100.23850.23240.22420.21400.2023
Akavci and Tanrikulu [31]≠ 00.23670.23640.23590.23530.23380.23000.22490.21820.2102
Younsi et al. [39]= 00.24160.24120.24040.23920.23660.23050.22220.21210.2003
Younsi et al. [39]≠ 00.23710.23690.23640.23570.23420.23040.22520.21860.2105
Present study= 00.24610.24570.24490.24370.24120.23510.22690.21680.2051
Present study≠ 00.23570.23540.23500.23430.23280.22910.22400.21740.2094
2Mantari and Soares [26]≠ 00.39270.39210.39080.38890.38490.37520.36210.34600.3273
Mantari and Soares [26]= 00.38100.38030.37900.37700.37300.36300.34970.33440.3165
Akavci and Tanrikulu [31]= 00.38960.38890.38770.38570.38170.37190.35880.34250.3237
Akavci and Tanrikulu [31]≠ 00.37900.37870.37790.37680.37440.36840.36020.34960.3368
Younsi et al. [39]= 00.38670.38600.38470.38280.37870.36890.35570.33940.3206
Younsi et al. [39]≠ 00.37970.37930.37860.37740.37500.36910.36080.35010.3373
Present study= 00.39390.39330.39200.39010.38600.37630.36320.34700.3282
Present study≠ 00.37740.37700.37630.37520.37280.36690.35870.34820.3355
3Mantari and Soares [26]≠ 00.44180.44110.43960.43750.43300.42210.40740.38930.3683
Mantari and Soares [26]= 00.42860.42790.42640.42420.41960.40840.39340.37610.3558
Akavci and Tanrikulu [31]= 00.43830.43760.43610.43400.42940.41850.40360.38540.3642
Akavci and Tanrikulu [31]≠ 00.42650.42610.42520.42390.42120.41460.40530.39340.3789
Younsi et al. [39]= 00.43500.43430.43280.43070.42610.41510.40020.38190.3607
Younsi et al. [39]≠ 00.42730.42680.42600.42470.42200.41530.40590.39400.3795
Present study= 00.44320.44250.44100.43890.43430.42340.40860.39040.3693
Present study≠ 00.42460.42420.42340.42210.41940.41280.40360.39180.3775
Table 7. Non-dimensional central displacement w ¯ ( 0 ) = G ( h ) w / h   q 0 and in-plane normal stress σ ¯ x ( 0 ) = σ x ( 0 ) / q 0 of E-FGM plates subjected to uniformly distributed load.
Table 7. Non-dimensional central displacement w ¯ ( 0 ) = G ( h ) w / h   q 0 and in-plane normal stress σ ¯ x ( 0 ) = σ x ( 0 ) / q 0 of E-FGM plates subjected to uniformly distributed load.
h/aQuantityTheoryE0/E1
0.10.51210
0.2 w ¯ Vaghefi et al. [58] (BEM)4.09168.975112.599017.664039.0600
Vaghefi et al. [58] (FEM)4.12159.004712.613017.711039.1550
Akavci and Tanrikulu [31] ( ε z 0 )3.83338.872412.597017.744038.3330
Younsi et al. [39] ( ε z 0 )3.83458.875612.602517.751138.3451
Present study ( ε z = 0)4.10119.108712.865318.217141.0098
Present study ( ε z 0 )3.82658.856012.574017.712338.2668
σ ¯ x
(−h/2)
Vaghefi et al. [58] (BEM)−15.356−9.2902−7.4462−5.9410−3.4665
Vaghefi et al. [58] (FEM)−15.403−9.2995−7.4588−5.9591−3.4805
Akavci and Tanrikulu [31] ( ε z 0 )−16.3220−9.6545−7.6944−6.1109−3.4530
Younsi et al. [39] ( ε z 0 )−16.2898−9.6313−7.6770−6.0994−3.4504
Present study ( ε z = 0)−15.6820−9.2913−7.3718−5.8141−3.2271
Present study ( ε z 0 )−16.6927−9.8955−7.8723−6.2318−3.4891
0.3 w ¯ Vaghefi et al. [58] (BEM)0.97072.13782.98534.12088.7134
Vaghefi et al. [58] (FEM)0.97322.14072.97924.13338.7293
Akavci and Tanrikulu [31] ( ε z 0 )0.89232.08342.96024.16698.9229
Younsi et al. [39] ( ε z 0 )0.89252.08432.96154.16858.9253
Present study ( ε z = 0)0.96022.17723.08224.35439.6015
Present study ( ε z 0 )0.89082.07982.95494.15958.9080
σ ¯ x Vaghefi et al. [58] (BEM)−7.223−4.3084−3.4496−2.7499−1.6449
Vaghefi et al. [58] (FEM)−7.2639−4.3378−3.4681−2.7673−1.6499
Akavci and Tanrikulu [31] ( ε z 0 )−7.6576−4.5062−3.5748−2.8235−1.5731
Younsi et al. [39] ( ε z 0 )−7.6386−4.4941−3.5659−2.8175−1.5715
Present study ( ε z = 0)−7.2499−4.2796−3.3846−2.6589−1.4605
Present study ( ε z 0 )−7.7999−4.5974−3.6421−2.8693−1.5869
Table 8. Comparison of the dimensionless deflection w of isotropic square plate subjected to uniformly distributed load.
Table 8. Comparison of the dimensionless deflection w of isotropic square plate subjected to uniformly distributed load.
K ¯ w K ¯ s a/h = 10a/h = 200
Thai et al. [24]Al Khateeb and Zenkour [28]Present 2DPresent Quasi-3DThai et al. [24]Al Khateeb and Zenkour [28]Present 2DPresent Quasi-3D
153.34553.180683.34523.33023.22003.219593.22003.2117
102.75042.619772.75032.74522.66842.668092.66842.6628
152.33312.22532.33302.33292.27632.276022.27632.2722
202.02441.933042.02432.02701.98341.983171.98341.9803
3452.84212.706992.84202.83582.75522.754852.75522.7491
102.39832.287652.39822.39772.33902.338662.33892.3346
152.07301.979632.07292.07542.03062.030372.03062.0274
201.82441.743941.82441.82861.79321.792981.79321.7907
5451.37851.323441.37841.38541.36881.368641.36881.3674
101.26151.211691.26141.26841.25431.254121.25421.2531
151.16271.117251.16271.16941.15721.157111.15721.1562
201.07821.036381.07821.08471.07401.073891.07401.0732
Table 9. Comparisons of non-dimensional displacements and stresses of simply supported FG rectangular plate (Al/Al2O3) resting on elastic foundation under uniform loads (b = 3a, a/h = 10).
Table 9. Comparisons of non-dimensional displacements and stresses of simply supported FG rectangular plate (Al/Al2O3) resting on elastic foundation under uniform loads (b = 3a, a/h = 10).
kK0J0Method u v w σ x σ y σxy
000Zenkour [20]0.19720.10221.25830.71620.24480.2893
Thai and Choi [29]0.19710.10221.25830.71600.24470.2890
Present 2D0.19720.10211.25820.71590.24420.2869
Present Quasi-3D0.19530.10091.25030.87550.40740.2823
1000Zenkour [20]0.19220.10031.22590.69700.23760.2843
Thai and Choi [29]0.19220.10031.22600.69690.23750.2840
Present 2D0.19220.10021.22590.69670.23700.2819
Present Quasi-3D0.19050.09911.21860.85230.39590.2774
0100Zenkour [20]0.18300.09671.16620.66190.22450.2746
Thai and Choi [29]0.18300.09671.16620.66180.22450.2744
Present 2D0.18300.09661.16610.66160.22390.2723
Present Quasi-3D0.18150.09551.15990.80970.37490.2681
100100Zenkour [20]0.17870.09511.13820.64530.21840.2702
Thai and Choi [29]0.17870.09511.13820.64520.21830.2700
Present 2D0.17870.09501.13810.64510.21780.2679
Present Quasi-3D0.17720.09391.13230.78950.36500.2638
0.500Zenkour [20]0.34920.18101.93440.23370.07990.0941
Thai and Choi [29]0.34910.18091.93450.23370.07990.0941
Present 2D0.34920.18071.93430.23360.07970.0934
Present Quasi-3D0.33460.17291.89950.27630.12860.0891
1000Zenkour [20]0.33580.17591.85900.22420.07630.0916
Thai and Choi [29]0.33580.17581.85900.22420.07630.0916
Present 2D0.33580.17561.85890.22420.07610.0910
Present Quasi-3D0.32210.16811.82710.26530.12310.0867
0100Zenkour [20]0.31200.16651.72480.20750.07010.0871
Thai and Choi [29]0.31190.16651.72480.20750.07010.0870
Present 2D0.31200.16631.72470.20740.06990.0864
Present Quasi-3D0.29970.15931.69800.24580.11360.0824
100100Zenkour [20]0.30130.16231.66400.19990.06730.0850
Thai and Choi [29]0.30120.16231.66400.19990.06730.0850
Present 2D0.30130.16211.66390.19980.06710.0843
Present Quasi-3D0.28960.15541.63940.23690.10920.0805
100Zenkour [20]0.48550.25152.51330.32500.11110.1307
Thai and Choi [29]0.48540.25152.51340.32500.11110.1306
Present 2D0.48540.25122.51320.32490.11080.1298
Present Quasi-3D0.45440.23472.42870.37520.17460.1209
1000Zenkour [20]0.46170.24242.38740.30810.10470.1263
Thai and Choi [29]0.46160.24242.38750.30800.10470.1262
Present 2D0.46160.24212.38730.30800.10450.1254
Present Quasi-3D0.43290.22652.31150.35630.16520.1170
0100Zenkour [20]0.42040.22622.17020.27910.09400.1183
Thai and Choi [29]0.42030.22612.17030.27910.09400.1182
Present 2D0.42040.22582.17010.27900.09380.1174
Present Quasi-3D0.39550.21182.10830.32380.14930.1097
100100Zenkour [20]0.40230.21912.07460.26630.08930.1148
Thai and Choi [29]0.40220.21902.07460.26630.08930.1148
Present 2D0.40230.21882.07450.26620.08900.1139
Present Quasi-3D0.37900.20542.01840.30930.14220.1066
200Zenkour [20]0.65650.34013.22670.43960.15020.1766
Thai and Choi [29]0.65640.34003.22660.43950.15020.1766
Present 2D0.65650.33973.22660.43940.14990.1755
Present Quasi-3D0.60250.31133.07060.49750.23150.1604
1000Zenkour [20]0.61570.32453.02190.41060.13940.1690
Thai and Choi [29]0.61560.32443.02180.41050.13940.1690
Present 2D0.61570.32413.02180.41040.13900.1679
Present Quasi-3D0.56690.29772.88540.46610.21600.1538
0100Zenkour [20]0.54760.29752.68140.36280.12170.1557
Thai and Choi [29]0.54750.29742.68140.36280.12170.1557
Present 2D0.54750.29712.68130.36270.12140.1546
Present Quasi-3D0.50690.27402.57480.41410.19040.1421
100100Zenkour [20]0.51870.28612.53640.34230.11420.1502
Thai and Choi [29]0.51860.28602.53640.34230.11420.1501
Present 2D0.51870.28562.53630.34220.11380.1490
Present Quasi-3D0.48130.26392.44150.39160.17940.1372
500Zenkour [20]0.78050.40453.85170.52240.17850.2104
Thai and Choi [29]0.78020.40433.85060.52230.17850.2103
Present 2D0.78040.40403.85160.52220.17810.2089
Present Quasi-3D0.71980.37213.68930.59410.27650.1919
1000Zenkour [20]0.72320.38253.56290.48160.16330.1997
Thai and Choi [29]0.72300.38243.56200.48160.16330.1996
Present 2D0.72310.38203.56280.48160.16290.1982
Present Quasi-3D0.66930.35283.42470.54960.25440.1825
0100Zenkour [20]0.63050.34563.09790.41680.13940.1815
Thai and Choi [29]0.63040.34553.09720.41680.13940.1814
Present 2D0.63040.34513.09780.41660.13900.1800
Present Quasi-3D0.58680.31992.99480.47810.21940.1663
100100Zenkour [20]0.59230.33042.90520.38970.12940.1741
Thai and Choi [29]0.59220.33032.90460.38970.12940.1740
Present 2D0.59220.32982.90500.38950.12900.1726
Present Quasi-3D0.55240.30632.81530.44810.20470.1597
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zaoui, F.Z.; Ouinas, D.; Achour, B.; Touahmia, M.; Boukendakdji, M.; Latifee, E.R.; Al-Naghi, A.A.A.; Viña Olay, J.A. Mathematical Approach for Mechanical Behaviour Analysis of FGM Plates on Elastic Foundation. Mathematics 2022, 10, 4764. https://doi.org/10.3390/math10244764

AMA Style

Zaoui FZ, Ouinas D, Achour B, Touahmia M, Boukendakdji M, Latifee ER, Al-Naghi AAA, Viña Olay JA. Mathematical Approach for Mechanical Behaviour Analysis of FGM Plates on Elastic Foundation. Mathematics. 2022; 10(24):4764. https://doi.org/10.3390/math10244764

Chicago/Turabian Style

Zaoui, Fatima Zohra, Djamel Ouinas, Belkacem Achour, Mabrouk Touahmia, Mustapha Boukendakdji, Enamur R. Latifee, Ahmed A. Alawi Al-Naghi, and Jaime Aurelio Viña Olay. 2022. "Mathematical Approach for Mechanical Behaviour Analysis of FGM Plates on Elastic Foundation" Mathematics 10, no. 24: 4764. https://doi.org/10.3390/math10244764

APA Style

Zaoui, F. Z., Ouinas, D., Achour, B., Touahmia, M., Boukendakdji, M., Latifee, E. R., Al-Naghi, A. A. A., & Viña Olay, J. A. (2022). Mathematical Approach for Mechanical Behaviour Analysis of FGM Plates on Elastic Foundation. Mathematics, 10(24), 4764. https://doi.org/10.3390/math10244764

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop