Next Article in Journal
Monitoring the Self-healing of Concrete from the Ultrasonic Pulse Velocity
Next Article in Special Issue
Optical and Mechanical Properties of Highly Transparent Glass-Flake Composites
Previous Article in Journal
Carbon Fiber Reinforced Thermoplastics Molding by Using Direct Resistance Heating to Carbon Nanofilaments Grafted Carbon Fiber
Previous Article in Special Issue
Thermal Buckling Behaviour of Thin and Thick Variable-Stiffness Panels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bending, Free Vibration, and Buckling Analysis of Functionally Graded Porous Micro-Plates Using a General Third-Order Plate Theory

1
Department of Civil and Construction Engineering, Western Michigan University, Kalamazoo, MI 49008, USA
2
Department of Mechanical and Aerospace Engineering, Western Michigan University, Kalamazoo, MI 49008, USA
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2019, 3(1), 15; https://doi.org/10.3390/jcs3010015
Submission received: 15 January 2019 / Revised: 27 January 2019 / Accepted: 28 January 2019 / Published: 1 February 2019
(This article belongs to the Special Issue Mechanics of Innovative Materials in Engineering Applications)

Abstract

:
Static bending, free vibration and buckling of functionally graded porous micro-plates are investigated using a general third order plate theory. In addition, analytical solutions are obtained using the Navier method. The effect of the material length scale factor and the variation of material property through the thickness direction of plates are considered as well as porosity effects. Three different porosity distributions are considered and the effects of porosity variations are examined in the framework of a general third order plate theory. Numerical results show that the effect of each distribution of porosity is distinguished due to coupling between the heterogeneity of the material properties and the variation of porosity.

1. Introduction

Micro electromechanical systems (MEMS) and Nano electromechanical systems (NEMS) are intelligent miniaturized systems used in many electronic devices, information/communication systems, chemical applications, and biomechanics [1]. The most of the structures used in MEMS and NEMS are in the form of beams and plates which can undergo deformation with the action of various loading cases. Studies show that the deformation of these small structures exhibits size dependency [2,3,4,5,6,7,8,9], which cannot be captured using the classical elasticity theory [3]. Hence, researchers have been studying new formulations which take into account the size dependency. The couple stress theory [10,11,12], Erigen’s nonlocal elasticity theory [13] and the strain gradient plasticity theories [2,14] are the examples of the continuum theories which can be used to capture the size dependent behavior of small scaled structures.
Functionally graded materials (FGMs) are composite materials that have a smooth variation of the material properties within a body. These novel materials have several advantages including reduced in-plane and transverse stresses, low residual stresses [15,16], high thermal resistance, low thermal conductivity, and high fracture toughness [17,18]. Also, contrary to classical laminated composites, FGMs do not exhibit interlaminar stresses because of the smooth variation of material composition [18]. In addition, the variation of the material properties can be optimized and adjusted for the requirements. Hence, many researchers have been taken advantages of these novel materials using them in various engineering fields.
Reddy [19] developed a microstructure-dependent nonlinear Euler-Bernoulli and Timoshenko beam models based on the modified couple stress theory [12]. Analytical solutions to bending, natural vibration and the buckling of homogeneous and functionally graded beams were provided. The newly developed model could capture the stiffening effect (size dependency) of small structures in bending and free vibration problems. It was shown that the difference in the deflections and frequencies decrease as the beam thickness increases.
Another study which also accounts for the size effect on functionally graded micro-beams was presented by Ke and Wang [7]. They investigated the influence of length scale parameter, power-law index, and length-to-thickness ratio on the dynamic stability of functionally graded micro-beams with several boundary conditions. Their parametric study showed that the micro-beams whose thicknesses were similar to their material length scale parameter experience a significant size effect in their free vibration, static buckling and dynamic characteristics. Simsek and Reddy analyzed functionally graded micro-beams with a higher order beam theory and modified couple stress theory [18]. Their study also emphasized that size effects became prominent as the thickness of the micro-beam and the material length-scale parameter are similar in magnitude. Li et al. [20] conducted a free vibration analysis of functionally graded beams with size-dependent Timoshenko beam model based on the strain gradient theory. Their study indicated that variation of material properties through the beam thickness had an important effect on the natural frequencies which could be used to control the natural frequencies. Their results also showed that the natural frequencies were in an increasing trend with increasing material length scale parameter.
In addition to FGM beams, mechanical responses of FGM plates have been investigated in the literature. Javaheri and Eslami [21] analyzed the bucking responses of functionally graded simply supported rectangular plate subjected to in-plane loading using the classical plate theory. Their result showed that the functionally graded plates resulted in smaller critical buckling loads than those of homogeneous plates. Hence, they noted that although functionally graded plates had numerous superior properties like heat resistance, their resistance should be checked against buckling failure. Nonlinear bending and the post-buckling behavior of functionally graded circular plates under mechanical and thermal loads were studied by Ma and Wang [22] based on von Kármán nonlinearity. Their study showed higher performance of functionally graded plates to resist thermal load or combination of thermal and mechanical load cases compared to fully metallic plate. They obtained lower temperature and deflection values than in pure metallic plates. Also, the effect of power-law index were discussed. Another study that investigated the thermal buckling of FGM plate was presented by Lanhe [23]. In this study, the equilibrium and stability equations of a simply supported moderately thick rectangular FGM plate were derived based on the first order shear deformation theory. The closed form solutions for the critical buckling temperature were obtained, and effects of plate aspect ratio, the relative thickness, the gradient index and the transverse shear on the buckling temperature were investigated. The results of the thick plates and plates with high aspect ratio showed that the critical buckling temperature was noticeably affected by the transverse shear deformation. In 2012, Reddy and Kim established a general third-order plate theory (GTPT) that accounted for the geometric nonlinearity, functionally graded material property, and microstructure-based size effects [24]. In GTPT, the material property variation was introduced in the plate constitutive relation. Also, the size effect was included by adopting Yang’s modified couple stress theory [12]. Later, Kim and Reddy presented the Navier solutions of the bending, free vibration and buckling problems of functionally graded plates using GTPT [25]. Similarly, Thai and Choi developed size-dependent functionally graded Kirchhoff and Mindlin plate models based on Yang’s modified couple stress theory [12]. The analytical solutions for deflection, buckling and vibration problems were presented in their study. Ansari et al. [26] investigated the size dependent vibration behavior of FG rectangular Mindlin micro-plates taking into account the geometric nonlinearity. Another size-dependent plate theory was developed by Zhang et al. [27]. Their model was based on the strain gradient theory and a newly developed refined shear deformation theory. The analyses of static bending, free vibration and buckling for functionally graded micro-plates resting on elastic foundation were conducted. Similar to the previous studies, the material length scale factor resulted in an increase in the stiffness of the functionally graded micro-plate. Also, the size dependent model was more influenced by the effect of the transverse shear deformation compared to the classical models which did not include material length scale parameter in the formulation.
It is shown that FGMs can be further improved in terms of weight reduction and energy absorption by introducing the porosity while maintaining a significant amount of strength [28,29]. There are various studies in the literature which show that the functionally graded porous materials (FGPMs) have promising results for several engineering applications such as enhanced filtration [30], automotive industry [31], and medical implants [32,33]. Hence, numerous researchers investigated the mechanical responses of the FGPMs, as the new class of composite materials.
Elastic buckling and static bending analysis of functionally graded porous beams were performed by Chen at al. [9] based on Timoshenko beam theory. Their study showed that as the porosity and the slenderness of the beam increased, the critical buckling load decreased and the maximum deflections obtained from bending analysis increased. Moreover, the buckling and the bending analysis were significantly affected by the porosity distribution. Later, Chen et al. [34] conducted free and forced vibration analysis of functionally graded beams with symmetric and asymmetric porosity distributions. They concluded that the symmetric porosity distribution resulted better beam stiffness; hence, higher fundamental frequency and lower dynamic deflection compared to the asymmetric one. Shafiei and Kazemi [35] presented the buckling analysis of two-dimensional functionally graded porous tapered Euler-Bernoulli nano- and micro-beams. The size effect was included using the modified couple stress theory. A modified series solution for free vibration analyses of moderately thick functionally graded porous deep curved and straight beams with general boundary conditions were proposed by Zhao et al. [36]. Their results showed that the highest frequencies were obtained by the FGP beams with symmetric porosity distribution, but the lowest frequencies were resulted by the FGP beams with uniform porosity distribution. In addition to geometric and material parameters, the boundary conditions strongly influenced the vibration characteristics of FGM beams. A finite element analysis was performed on free and forced vibration of functionally graded porous beam and frame structures by Wu et al. [37]. Their approach was tested with several numerical examples and the results were compared with both analytical and experimental results.
Shahverdi and Barati [38] developed a general nonlocal strain gradient elasticity model, and the vibration analysis of functionally graded porous nano-plates were carried out. To obtain more accurate predictions, two scale parameters were included in their proposed model. The derived governing equations of nano-plate were solved for natural frequencies by using the Galerkin method. It was demonstrated that the porosity had an important effect on mechanical behavior. They concluded that both nonlocal and strain gradient parameters were needed to capture both stiffness-softening and stiffness-hardening effects in the modeling of nano-plates. Wang et al. [39] investigated the vibrations of the longitudinally traveling functionally graded plates with both evenly and unevenly distributed porosity. The effects of porosity, damping and constituent fraction on the dynamic response of FGPM plates were examined. An analytical approach for natural frequencies of the functionally graded porous plates were established by Rezaei et al. [40]. Their model was based on the four variable refined plate theory of Thai an Choi [41], which is a modified first order shear theory. The effect of porosity, power-law index, thickness-side ratio, aspect ratio, porosity distribution and the boundary conditions of the rectangular plate were examined. Akbas [42] presented a vibration and the static bending analysis of a simply supported functionally graded porous rectangular plate within the framework of first order shear deformation theory. Wang and Zu [43] investigated the large-amplitude vibration of sigmoid functionally graded thin plates with porosities. The geometric nonlinearities were included by adopting von Kármán nonlinear plate theory. Their results indicated that even and uneven distribution types of porosity had different effects on the resonance domain of sigmoid functionally graded plates. Besides, the effect of the constituent volume distribution was clearly shown. The nonlinear thermomechanical buckling and post buckling response of porous functionally graded plates were studied by Cong et al. [44]. Their formulations were based on Reddy’s higher-order shear deformation theory. Their results demonstrated that better performance of the plate with evenly distributed porosity in terms of buckling and post-buckling analysis. Also, Pasternak shear layer stiffness gave a significant effect over Winkler elastic modulus on the porous functionally graded plate nonlinear stability. The transient response of porous functionally graded nano-plates under various pulse loads were investigated by Mirjavadi et al. [45] in the framework of nonlocal strain gradient theory. The effects of nonlocal coefficient, strain gradient parameter, porosity, type of impulse and the loading time on the transient response were discussed. Zhao et al. [46] established a novel three-dimensional exact solution for vibration analysis of thick FGPM plates with three different porosity distributions. Then, Zhao et al. [47] conducted an extensive study on the dynamic analysis of FGP circular, annular an sector plates with general elastic restraints. Thang et al. [48] performed elastic buckling and free vibration analyses of porous-cellular plates with various porosity distributions. Bending, free vibration and buckling analyses of FGP micro-plates were conducted by Kim et al. [49]. In their study, both the classical and first order shear theory were used with modified couple stress theory. Because they only studied the classical and first order shear deformation theory, their models cannot properly account for the parabolic transverse shear deformation.
Contrary to majority of the previous studies on functionally graded porous plates, this work presents the analysis of FGP plates with a higher order shear deformation theory that provides more accurate transverse shear stresses. Moreover, by using the general third-order plate theory [24], this study accounts for the thickness deformation of plates. In addition, the effect of porosity distribution, the size-dependency and the smooth change in the material property through the thickness direction of the plate are considered. The effect of several types of porosity distributions to the FGM micro-plates is investigated more precisely, and the numerical results of the bending, free vibration and bucking analysis are presented in this work.

2. General Third-Order Plate Theory

The classical plate theory (CPT) is based on the assumption that straight lines perpendicular to the mid-plane of the plate in the undeformed configuration remains straight, inextensible and perpendicular to the mid-plane of the plate after the deformation. The first order shear deformation plate theory +(FSDT) releases the normality condition. However, it requires a shear correction factor because the resulted transverse shear strains are nonzero but independent of the thickness coordinate, z. In addition to the normality condition, the general third order shear deformation plate theory (GTPT) [50] further releases the straight and the inextensibility conditions by expanding the displacement fields as cubic functions of z. Hence, the necessity of the shear correction factor is avoided by achieving the desired quadratic variation of the transverse shear strains and stresses through the plate thickness.

2.1. Displacements and Strains

In the GTPT, displacement field is extended up to third power of thickness coordinate for in plane displacements and second power of thickness coordinate for out of plane displacement [24]. The assumed displacement field is
u 1 ( x , y , z , t ) = u ( x , y , t ) + z θ x ( x , y , t ) + z 2 ϕ x ( x , y , t ) + z 3 ψ x ( x , y , t ) u 2 ( x , y , z , t ) = v ( x , y , t ) + z θ y ( x , y , t ) + z 2 ϕ y ( x , y , t ) + z 3 ψ y ( x , y , t ) u 3 ( x , y , z , t ) = w ( x , y , t ) + z θ z ( x , y , t ) + z 2 ϕ z ( x , y , t )
where u, v, w, θ x , θ y , θ z , ϕ x , ϕ y , ϕ z , ψ x , and ψ y are unknown generalized displacements. The linearized strains of the GTPT are (see Reddy and Kim [24]):
ε x x ε y y ε z z γ x y γ x z γ y z = ε x x ( 0 ) ε y y ( 0 ) ε z z ( 0 ) γ x y ( 0 ) γ x z ( 0 ) γ y z ( 0 ) + z ε x x ( 1 ) ε y y ( 1 ) ε z z ( 1 ) γ x y ( 1 ) γ x z ( 1 ) γ y z ( 1 ) + z 2 ε x x ( 2 ) ε y y ( 2 ) ε z z ( 2 ) γ x y ( 2 ) γ x z ( 2 ) γ y z ( 2 ) + z 3 ε x x ( 3 ) ε y y ( 3 ) ε z z ( 3 ) γ x y ( 3 ) γ x z ( 3 ) γ y z ( 3 )
where the nonzero strains are given as
ε x x ( 0 ) = u x , ε x x ( 1 ) = θ x x , ε x x ( 2 ) = ϕ x x , ε x x ( 3 ) = ψ x x , ε y y ( 0 ) = v y , ε y y ( 1 ) = θ y y , ε y y ( 2 ) = ϕ y y , ε y y ( 3 ) = ψ y y , ε z z ( 0 ) = θ z , ε z z ( 1 ) = 2 ϕ z , γ x y ( 0 ) = u y + v x , γ x y ( 1 ) = θ x y + θ y x , γ x y ( 2 ) = ϕ x y + ϕ y x , γ x y ( 3 ) = ψ x y + ψ y x , γ x z ( 0 ) = θ x + w x , γ x z ( 1 ) = 2 ϕ x + θ z x , γ x z ( 2 ) = 3 ψ x + ϕ z x , γ y z ( 0 ) = θ y + w y , γ y z ( 1 ) = 2 ϕ y + θ z y , γ y z ( 2 ) = 3 ψ y + ϕ z y .

2.2. Modified Couple Stress Model

In the classical couple stress theory [10,11], a couple that leads to the particle to rotate as well as a force to translate are considered as the applied loads on material particle. Equilibrium equations of the forces and the moments of the forces are satisfied as those in the case of classical elasticity. The strain and curvature tensors are the deformation tensors which are conjugated with the stress and the couple stress tensors, respectively. In addition to the Lamé constants, two independent length scale parameters exist associated with the symmetric and the antisymmetric parts of the curvature tensor in the constitutive relation. However, it requires extensive experiments to determine these two length scale parameters. In 2002, Yang et al. proposed the modified couple stress theory [12] and showed that only one single length scale parameter is adequate to include size effect for linear isotropic elastic materials. They achieved this result by enforcing higher order equilibrium equations which restrict the couple stress tensor to be symmetric [12]. Hence, only the symmetric part of the curvature tensor contributes to the strain energy density and the only one length scale parameter associated with the symmetric part of the curvature tensor is required [12,24].
The expression for the virtual strain energy δ U using the modified couple stress theory is [24]
δ U = V δ ε : σ + δ χ : m d v = V σ i j δ ε i j + m i j δ χ i j d v
where summation on repeated indices is implied; here σ i j and m i j denote the Cartesian components of symmetric stress tensor and the deviatoric couple stress tensor, respectively. For the linear elastic solid material, they can be expressed as the following form
σ i j = λ ε i i δ i j + 2 μ ε i j , m i j = 2 μ l 2 χ i j
where λ and μ are the Lamé constants, l is the material length scale parameter that can represents the size dependency in small scale structures, and δ i j is the Kronecker’s delta. The components of the linear strain tensor, i.e., ε i j are given in Equations (2) and (3). The components of the symmetric curvature tensor χ i j are
χ = 1 2 ω + ω T , ω = 1 2 × u
or
χ i j = 1 2 ω i x j + ω j x i , i , j = 1 , 2 , 3
and ω i ( i = 1 , 2 , 3 ) are the components of the rotation vector can be expressed as
ω x = ω 1 = 1 2 u 3 x 2 u 2 x 3 , ω y = ω 2 = 1 2 u 1 x 3 u 3 x 1 , ω z = ω 3 = 1 2 u 2 x 1 u 1 x 2 .
Substituting the rotation vector ω i ( i = 1 , 2 , 3 ) into the above vector equation, the components of the symmetric curvature tensor result
χ x x = χ 11 = 1 2 2 u 3 x 1 x 2 2 u 2 x 1 x 3 , χ y y = χ 22 = 1 2 2 u 1 x 2 x 3 2 u 3 x 1 x 2 , χ z z = χ 33 = 1 2 2 u 2 x 1 x 3 2 u 1 x 2 x 3 , χ x y = 2 χ 12 = 1 2 2 u 3 x 2 2 2 u 2 x 2 x 3 + 2 u 1 x 1 x 3 2 u 3 x 1 2 , χ x z = 2 χ 13 = 1 2 2 u 3 x 2 x 3 2 u 2 x 3 2 + 2 u 2 x 1 2 2 u 1 x 1 x 2 , χ y z = 2 χ 23 = 1 2 2 u 1 x 3 2 2 u 3 x 1 x 3 + 2 u 2 x 1 x 2 2 u 1 x 2 2 .
which can be written as
χ x x χ y y χ z z χ x y χ x z χ y z = χ x x ( 0 ) χ y y ( 0 ) χ z z ( 0 ) χ x y ( 0 ) χ x z ( 0 ) χ y z ( 0 ) + z χ x x ( 1 ) χ y y ( 1 ) χ z z ( 1 ) χ x y ( 1 ) χ x z ( 1 ) χ y z ( 1 ) + z 2 χ x x ( 2 ) χ y y ( 2 ) χ z z ( 2 ) χ x y ( 2 ) χ x z ( 2 ) χ y z ( 2 ) + z 3 χ x x ( 3 ) χ y y ( 3 ) χ z z ( 3 ) χ x y ( 3 ) χ x z ( 3 ) χ y z ( 3 )
where the nonzero components in terms of generalized displacements are obtained as
χ x x ( 0 ) = 1 2 x w y θ y , χ x x ( 1 ) = 1 2 x θ z y 2 ϕ y , χ x x ( 2 ) = 1 2 x ϕ z y 3 ψ y , χ y y ( 0 ) = 1 2 y θ x w x , χ y y ( 1 ) = 1 2 y 2 ϕ x θ z x , χ y y ( 2 ) = 1 2 x 3 ψ x ϕ z x , χ z z ( 0 ) = 1 2 θ y x θ x y , χ z z ( 1 ) = ϕ y x ϕ x y , χ z z ( 2 ) = 3 2 ψ y x ψ x y , χ x y ( 0 ) = 1 2 y w y θ y + x θ x w x , χ x y ( 1 ) = 1 2 y θ z y 2 ϕ y + x 2 ϕ x θ z x , χ x y ( 2 ) = 1 2 y ϕ z y 3 ψ y + x 3 ψ x ϕ z x , χ x z ( 0 ) = 1 2 x v x u y + θ z y 2 ϕ y , χ x z ( 1 ) = 1 2 x θ y x θ x y + 2 ϕ z y 3 ψ y , χ x z ( 2 ) = 1 2 x ϕ z y 3 ψ y + x 3 ψ x ϕ z x , χ x z ( 3 ) = 1 2 x ψ y x ψ x y , χ y z ( 0 ) = 1 2 y v x u y + 2 ϕ x θ z x , χ y z ( 1 ) = 1 2 y θ y x θ x y + 2 3 ψ x ϕ z x , χ y z ( 2 ) = 1 2 y ϕ z y 3 ψ y + x 3 ψ x ϕ z x , χ y z ( 3 ) = 1 2 y ψ y x ψ x y .

2.3. FGPM Plate Constitutive Equations

In this study, isotropic plates with variation of two constituents and distribution of porosity through thickness are assumed. The considered three different porosity distributions [49] are
Type 1 : Φ ( z ) = ϕ cos π z h Type 2 : Φ ( z ) = ϕ cos π 2 z h + 1 2 Type 3 : Φ ( z ) = ϕ cos π 2 z h 1 2
where ϕ is the maximum porosity value achieved by the distribution. A typical material property of the functionally graded porous materials (FGPM) can be considered as in the following power-law relations [49],
P ( z ) = P t P b z h + 1 2 n + P b 1 Φ ( z )
where P t and P b are the typical material properties at the top and the bottom surfaces, respectively and n is the power-law index.
Figure 1 shows the normalized porosity and the material property distributions through the plate thickness. For the illustration purpose, the material property at the top surface is assumed as ten times larger than the material property at the bottom surface of the plate and the porosity value of ϕ = 0.5 is taken. As can be seen in Figure 1a, all three types of porosity distribution result in the same volume of void through the material. However, the Type 1 distribution is symmetric about the mid-plane of the plate. Also, the middle section of the plate is more porous and the porosity of the material decreases through the top and the bottom surfaces in the case of Type 1 distribution. On the other hand, this symmetry is not valid for porosity distribution Type 2 and 3. In the case of Type 2, the porosity increases zero at the top surface to assigned ϕ value at the bottom surface of the plate. Type 3 porosity distribution has the opposite variation to the Type 2.
Figure 1b–d demonstrate the effect of power-law index, n and porosity distribution types on the variation of material properties. As can be seen from Figure 1b–d, the distributions of the normalized material property with Type 1 and Type 2 show similar trend, but Type 3 exhibits quite different variation. This is because the material on top surface is assumed to be stiffer than the material on bottom surface.
Under isothermal conditions, the linear constitutive relation for an isotropic plate is given by
σ x x σ y y σ z z σ x y σ x z σ y z = E ( 1 + ν ) ( 1 2 ν ) 1 ν ν ν 0 0 0 ν 1 ν ν 0 0 0 ν ν 1 ν 0 0 0 0 0 0 1 2 1 2 ν 0 0 0 0 0 0 1 2 1 2 ν 0 0 0 0 0 0 1 2 1 2 ν ε x x ε y y ε z z ε x y ε x z ε y z
In the present study, Young’s modulus, E, varies through thickness direction but Poisson’s ratio, ν , is assumed to be a constant in the constitutive Equation (12). The generalized forces and couples can be expressed as then
M x x i M y y i M z z i = h 2 h 2 σ x x σ y y σ z z z i d z = k = i 3 + i A 11 k A 12 k A 12 k A 12 k A 11 k A 12 k A 12 k A 12 k A 11 k ε x x k i ε y y k i ε z z k i
M x y i M x z i M y z i = h 2 h 2 σ x y σ x z σ y z z i d z = k = i 3 + i B 11 k 0 0 0 B 11 k 0 0 0 B 11 k γ x y k i γ x z k i γ z z k i
M x x i M y y i M z z i = h 2 h 2 m x x m y y m z z z i d z = k = i 2 + i B 11 k 0 0 0 B 11 k 0 0 0 B 11 k χ x x k i χ y y k i χ z z k i
M x y i M x z i M y z i = h 2 h 2 m x y m x z m y z z i d z = k = i 3 + i B 11 k 0 0 0 B 11 k 0 0 0 B 11 k χ x y k i χ x z k i χ y z k i
where m i j and σ i j are expressed in Equation (5), and plate stiffness A 11 , A 12 , B 11 , and B 11 are
A 11 ( k ) = 1 ν 1 + ν 1 2 ν h 2 h 2 z k E z d z , A 12 ( k ) = ν 1 + ν 1 2 ν h 2 h 2 z k E z d z B 11 ( k ) = 1 2 1 + ν h 2 h 2 z k E z d z , B 11 ( k ) = l 2 1 + ν h 2 h 2 z k E z d z .

2.4. Equation of Motion

Reddy and Kim [24] use the principle of virtual displacements to obtain the equations of motion of GTPT. The dynamic case of the principle of virtual displacements is generalized to Hamilton’s principle [51].
0 T δ K δ U δ V d t = 0
where δ K is the virtual kinetic energy, δ U is the virtual strain energy, and δ V is the virtual work done by external forces. Each virtual energy terms and the virtual work are defined as
δ K = Ω h 2 h 2 ρ u i t δ u i t d z d x d y
δ U = Ω h 2 h 2 σ i j δ ε i j + m i j δ χ i j d z d x d y
δ V = Ω h 2 h 2 f ¯ i δ u i + c ¯ i δ ω i d z d x d y + Ω + q i t δ u i + p i t δ ω i d x d y + Ω q i b δ u i + p i b δ ω i d x d y + S t ¯ i δ u i + s ¯ i δ ω i d S
and ρ is mass density, f ¯ i is the body force measured per unit volume, c ¯ i is the body couples measured per unit volume, q i t , p i t , q i b , and p i b are distributed forces and couples on the top surface ( Ω + ) and bottom surface ( Ω ), respectively, t ¯ i and s ¯ i are surface forces and couples measured per unit area acting on lateral surfaces of a plate, respectively, with the subscript i taking the values of x, y, or z. The Ω indicates the mid surface of plates. For a rectangular plate in the Cartesian coordinate system, Ω is define as the x y plane at z = 0 , Ω + is defined at z = + h 2 , and Ω is defined at z = h 2 .
Substituting the generalized forces Equations (13) to (16) into Equation (20) and performing integration by parts, the equation of motion of the GTPT in terms of stress resultants are obtained as [24].
δ u : M x x ( 0 ) x + M x y ( 0 ) y + 1 2 y M x z ( 0 ) x + M y z ( 0 ) y + F x ( 0 ) + 1 2 c z ( 0 ) y = I 0 u ¨ + I 1 θ ¨ x + I 2 ϕ ¨ x + I 3 ψ ¨ x
δ v : M x y ( 0 ) x + M y y ( 0 ) y 1 2 x M x z ( 0 ) x + M y z ( 0 ) y + F y ( 0 ) 1 2 c z ( 0 ) x = I 0 v ¨ + I 1 θ ¨ y + I 2 ϕ ¨ y + I 3 ψ ¨ y
δ w : M x z ( 0 ) x + M y z ( 0 ) y 1 2 y M x x ( 0 ) x + M x y ( 0 ) y + 1 2 x M y y ( 0 ) y + M x y ( 0 ) x + F z ( 0 ) + 1 2 c y ( 0 ) x c x ( 0 ) y = I 0 w ¨ + I 1 θ ¨ z + I 2 ϕ ¨ z
δ θ x : M x x ( 1 ) x + M x y ( 1 ) y M x z ( 0 ) + 1 2 M x y ( 0 ) x + M y y ( 0 ) y M z z ( 0 ) y + 1 2 y M x z ( 1 ) x + M y z ( 1 ) y + F x ( 1 ) + 1 2 c y ( 0 ) + 1 2 c z ( 1 ) y = I 1 u ¨ + I 2 θ ¨ x + I 3 ϕ ¨ x + I 4 ψ ¨ x
δ θ y : M x y ( 1 ) x + M y y ( 1 ) y M y z ( 0 ) 1 2 M x x ( 0 ) x + M x y ( 0 ) y M z z ( 0 ) x 1 2 x M x z ( 1 ) x + M y z ( 1 ) y + F y ( 1 ) 1 2 c x ( 0 ) 1 2 c z ( 1 ) x = I 1 v ¨ + I 2 θ ¨ y + I 3 ϕ ¨ y + I 4 ψ ¨ y
δ θ z : M x z ( 1 ) x + M y z ( 1 ) y M z z ( 0 ) 1 2 y M x x ( 1 ) x + M x y ( 1 ) y + 1 2 x M y y ( 1 ) y + M x y ( 1 ) x + 1 2 M x z ( 0 ) y M y z ( 0 ) x + F z ( 1 ) + 1 2 c y ( 1 ) x c x ( 1 ) y = I 1 w ¨ + I 2 θ ¨ z + I 3 ϕ ¨ z
δ ϕ x : M x x ( 2 ) x + M x y ( 2 ) y 2 M x z ( 1 ) + M x y ( 1 ) x + M y y ( 1 ) y M z z ( 1 ) y M y z ( 0 ) + 1 2 y M x z ( 2 ) x + M y z ( 2 ) y + F x ( 2 ) + 1 2 c y ( 1 ) + 1 2 c z ( 2 ) y = I 2 u ¨ + I 3 θ ¨ x + I 4 ϕ ¨ x + I 5 ψ ¨ x
δ ϕ y : M x y ( 2 ) x + M y y ( 2 ) y 2 M y z ( 1 ) M x x ( 1 ) x + M x y ( 1 ) y M z z ( 1 ) x M x z ( 0 ) 1 2 x M x z ( 2 ) x + M y z ( 2 ) y + F y ( 2 ) 1 2 c x ( 1 ) 1 2 c z ( 2 ) x = I 2 v ¨ + I 3 θ ¨ y + I 4 ϕ ¨ y + I 5 ψ ¨ y
δ ϕ z : M x z ( 2 ) x + M y z ( 2 ) y 2 M z z ( 1 ) 1 2 y M x x ( 2 ) x + M x y ( 2 ) y + 1 2 x M y y ( 2 ) y + M x y ( 2 ) x + M x z ( 1 ) y M y z ( 1 ) x + F z ( 2 ) + 1 2 c y ( 2 ) x c x ( 2 ) y = I 2 w ¨ + I 3 θ ¨ z + I 4 ϕ ¨ z
δ ψ x : M x x ( 3 ) x + M x y ( 3 ) y 3 M x z ( 2 ) + 3 2 M x y ( 2 ) x + M y y ( 2 ) y M z z ( 2 ) y 2 M y z ( 1 ) + 1 2 y M x z ( 3 ) x + M y z ( 3 ) y + F x ( 3 ) + 3 2 c y ( 2 ) + 1 2 c z ( 3 ) y = I 3 u ¨ + I 4 θ ¨ x + I 5 ϕ ¨ x + I 6 ψ ¨ x
δ ψ y : M x y ( 3 ) x + M y y ( 3 ) y 3 M y z ( 2 ) 3 2 M x x ( 2 ) x + M x y ( 2 ) y M z z ( 2 ) x 2 M x z ( 1 ) 1 2 x M x z ( 3 ) x + M y z ( 3 ) y + F y ( 3 ) 3 2 c x ( 2 ) 1 2 c z ( 3 ) x = I 3 v ¨ + I 4 θ ¨ y + I 5 ϕ ¨ y + I 6 ψ ¨ y
where F ξ i = f ¯ ξ i + h 2 i q ξ t + 1 i q ξ b , f ¯ i = h 2 h 2 z i f ξ d z , c ξ ( i ) = h 2 h 2 z i c ¯ ξ d z for i = 0 , 1 , 2 , 3 and ξ = x , y , z . The superposed dot are used to indicate the time derivative, hence u ¨ = 2 u / t 2 . The mass resultants I 0 , I 1 , …, I 6 are defined as
I i = h / 2 h / 2 ρ ( z ) i d z , i = 0 , 1 , 6
The equations of motion for GTPT in terms of generalized displacements can be found in the references [24,25].

3. Solution Procedure

Navier solution technique is applied to obtain analytical solutions to static bending, free vibration and buckling problems for a simply supported rectangular FGPM plate. The boundary conditions of the simply supported plate are demonstrated in Figure 2a. Figure 2b shows the cross-sectional view of typical FGP plate used in this study (e.g., Type 3 is illustrated in Figure 2b). Also, x z cross section is identical to y z cross section of the plates.
To employ the Navier solution technique, the generalized displacements (u, v, w, θ x , θ y , θ z , ϕ x , ϕ y , ϕ z , ψ x , and ψ y ) and applied load ( q z ) are expanded in double trigonometric series which are selected to satisfy the boundary conditions of the simply supported plate problem. The following double trigonometric series are used to expand the generalized displacements
u x , y , z , t = m = 1 n = 1 U m n t cos α x sin β y v x , y , z , t = m = 1 n = 1 V m n t sin α x cos β y w x , y , z , t = m = 1 n = 1 W m n t sin α x sin β y θ x x , y , z , t = m = 1 n = 1 Θ x m n t cos α x sin β y θ y x , y , z , t = m = 1 n = 1 Θ y m n t sin α x cos β y θ z x , y , z , t = m = 1 n = 1 Θ z m n t sin α x sin β y ϕ x x , y , z , t = m = 1 n = 1 Φ x m n t cos α x sin β y ϕ y x , y , z , t = m = 1 n = 1 Φ y m n t sin α x cos β y ϕ z x , y , z , t = m = 1 n = 1 Φ z m n t sin α x sin β y ψ x x , y , z , t = m = 1 n = 1 Ψ x m n t cos α x sin β y ψ y x , y , z , t = m = 1 n = 1 Ψ y m n t sin α x cos β y
where α = m π a , β = n π b . The generalized displacement coefficients U m n , V m n , W m n , Θ x m n , Θ y m n , Θ z m n , Φ x m n , Φ y m n , Φ z m n , Ψ x m n , and Ψ y m n are treated as time independent variables for static bending and buckling problems.
For the static bending problem, a uniformly distributed transverse load q z is defined as acting to the top surface of the plate, which is also expanded in double trigonometric series in Equation (35). Because q z is assumed as an uniformly distributed load, the coefficient for the applied force is simplified as shown in Equation (36).
q z = m = 1 n = 1 Q m n t sin α x sin β y
where
Q m n = 4 a b 0 a 0 b q z sin α x sin β y d x d y = 16 q 0 m n π 2
where q 0 is the magnitude of the applied uniformly distributed load.
After substituting these expanded forms into the equation of motions (22) to (23), the algebraic system of equations are obtained for static bending problem as
C m n Δ m n = F m n
where C is the 11 by 11 coefficient matrix whose expanded form can be found in the Kim and Reddy’s study [25], Δ is the vector of generalized displacements, and F is the force vector given in Equation (38).
Δ m n = U m n V m n W m n Θ x m n Θ y m n Θ z m n Φ x m n Φ y m n Φ z m n Ψ x m n Ψ y m n , F m n = 0 0 Q m n 0 0 h 2 Q m n 0 0 h 2 4 Q m n 0 0
Equation (37) is solved for each m = 1 , 3 , 5 , and n = 1 , 3 , 5 , in order to find generalized displacements Δ . Then, generalized displacements for each m and n are inserted into the double trigonometric form of displacements given in Equation (34) to obtain displacements u, v, w, θ x , θ y , θ z , ϕ x , ϕ y , ϕ z , ψ x , and ψ y .
In the case of free vibration analysis, the time dependent generalized displacements in Equation (34) are assumed as
Δ ( t ) m n = e i ω t Δ m n
where ω is the natural frequency. By substituting Equations (34) and (39) into the equation of motion, we setup the eigenvalue problem to determine eigenfrequencies:
C m n ω m n 2 M m n Δ m n = 0
where C is the coefficient matrix which is the same as the one in the bending analysis, M is the matrix of inertias (please see [25]), and Δ is the vector of generalized displacements. The free vibration analysis of simply supported square FGPM plate is performed by solving Equation (40). The fundamental frequencies are then obtained when m = 1 and n = 1 .
For buckling analysis, the simply supported plate is assumed under in-plane compressive forces acting on x and y edges which are N ^ x x and N ^ y y , respectively. Then, additional N ^ x x 2 w x 2 + N ^ y y 2 w y 2 terms are inserted in right hand side of Equation (24). The following system of equations is obtained to find the buckling loads N 0 ;
C m n N m n Δ m n = 0
where N is the matrix with only nonzero element is N 33 = N 0 ( α 2 + k β 2 ) , k = N ^ y y N ^ x x . After performing a static condensation to Equation (41) and setting m = 1 and n = 1 , the critical buckling load is found as
N 0 = 1 α 2 + k β 2 C 0303 C 0301 C 0302 C 0304 C 0305 C 0310 C 0311 T C ¯ C 0103 C 0203 C 0403 C 0503 C 1003 C 1103 .
where
C ¯ = C 0101 C 0102 C 0104 C 0105 C 0110 C 0111 C 0201 C 0202 C 0204 C 0205 C 0210 C 0211 C 0401 C 0402 C 0404 C 0405 C 0410 C 0411 C 0501 C 0502 C 0504 C 0505 C 0510 C 0511 C 1001 C 1002 C 1004 C 1005 C 1010 C 1011 C 1101 C 1102 C 1104 C 1105 C 1110 C 1111 1 .
In the numerical solutions, it is assumed that the square FGPM plate is under the same amount of compressive force acting in x and y edges (i.e., their ratio, k = 1 ).

4. Numerical Results and Discussions

Numerical examples of the analytical solutions given in the previous section are obtained using the material properties and the dimensions of the functionally graded plate from the study of Kim and Reddy [25]. The dimensions of functionally graded plate in Figure 2 are a = 20 h , b = 20 h , and h = 17.6 × 10 6 m. The moduli and mass densities of two constituents are E t = 14.4 GPa, E b = 1.44 GPa, ρ t = 12.2 × 10 3 kg/m, and ρ b = 1.22 × 10 3 kg/m a with the Poisson ratio ( ν ) of 0.38.

4.1. Verification

To show the consistency of the present study with the literature, the static bending, free vibration and buckling results are compared with the results of the previously published papers. One of the selected paper is based on the general third order plate theory [25], a generalized third order shear deformation theory (TSDT) for plates, for functionally graded plates without considering the porosity effect. The second paper has the same porosity distributions with the current study; however, it is based on classical plate theory (CPT) and first order shear deformation theory (FSDT) [49]. As can be seen from Figure 3, the present study shows good agreement with the previous studies. For non-porous cases, the current results and the Kim’s results are identical because they both uses the GTPT. On the other hand, the comparison with the second paper shows a small deviation resulted from the differences in the applied theories.

4.2. Static Bending

In the static bending problem, the uniformly distributed load of q 0 = 1 N/m 2 is applied to the top surface of the simply supported FGP plate. The generalized displacements and the applied load are expanded in double trigonometric series which are given in Equations (34) and (35) up to m = 31 and n = 31 terms. Then, the system equation given in Equation (37) is solved.
The nondimensional transverse deflections at y = b / 2 along x for porous plates without variation of material constituents are shown in Figure 4a. The length scale parameter effect is also neglected to see the effect of porosity variations clearly. Because the power-law index is set to zero (see Equation (11)), the plate consists of a single material which is the top surface material. As expected, while the porosity of the plate increases, the bending stiffness of the plate decreases. Hence, larger deflection values are obtained for more porous plates. It is important to remind that although distribution types are different, the total porosity volume of the plates result from them are the same. Hence, a reasonable comparison of the type of porosity distribution can be conducted. Another important point is that in the case of Type 1 distribution of homogeneous plate, we have symmetric variation of the porosity distribution about the mid-plane of the plate, which decrease the effect of coupling terms in the coefficient matrix given in Equation (37). On the other hand, this is not the case for Type 2 and Type 3 distributions even for the n = 0 case.
Figure 4 shows the center deflections of simply supported homogeneous plates with various porosity variations and FGM plates with various power-law indices. Figure 4a clearly shows that the porosity distribution types result in different bending rigidity; although they have the same amount of porosity values. The plate with Type 1 distribution are stiffer than the those with Type 2 and 3 with same porosity. The reason is that Type 1 distribution results in high porosity values at the mid-plate and the very low values as one moves through the top and the bottom surfaces (please see Figure 1a). Furthermore, the transverse deflection values for Type 2 and 3 are the same, because the plate is homogeneous and the resultant material property of them are the same. The only difference is that porous volume increases from top surface to the bottom for Type 2, while it increases from bottom to the top surface of the plate for the case of Type 3. It makes no difference in terms of bending rigidity for homogeneous plates.
Figure 4b shows the nondimesional transverse deflection along x while y = b / 2 for functionally graded porous plate without including the length scale parameter. High deflection values are obtained as the power-law index increases for each distribution. The reason is that while the power-law index increases, the material property distributions become highly asymmetric and therefore the effects of coupling terms come into the picture and start to influence the results dramatically. Also for larger values of power-law index, the defections of Type 2 and Type 3 start to deviate from each other. For example, when n = 5 , Type 2 shows the higher transverse deflections than Type 3. Because Type 2 is the most sensitive porosity distribution type for the given material configuration (i.e., E t / E b = 10 ) and the most influenced type by the effect of power-law index. On the other hand, Type 3 is least effected one by the change of n.
The effect of the length scale parameter is presented in Figure 5 for FGPM plate. This figure clearly shows that the length scale parameter has an important contribution to capture the deformation behavior of the micro plates. Indeed, to find the more accurate results from the numerical analysis, the length scale factor for the micro-plates should be determined first. Regardless of the porosity and the variations in the constituents, the plate stiffness increases as the length parameter increases. For a porous plate without material variation ( n = 0 ), the effect of porosity distribution type becomes negligible for length scale parameter values higher than 0.6 h , because the influence of the length scale parameter has a dominant effect in this regime. Similarly, for functionally graded plates ( n 0 ), the central deflection values obtained from Type 1 and 2 become quite similar for large values of length scale parameter (larger than 0.6 h ). On the other hand, the effect of length scale parameter is not much dominant for Type 3 compared to Type 1 and 2. Moreover, for large values of power-law index and the small length parameter, Type 2 results high transverse deflection values as compared to Type 1 and 3, due to differences in the effects of coupling terms.
For an illustration of stress variations of higher order shear deformation theory, the bending stresses and transverse shear stresses of the porosity distribution Type 1 are presented in Figure 6. The black line without a marker shows the variations of stresses in homogeneous materials ( n = 0 ). It clearly shows that the distribution of bending stress of higher order plate theory is similar to ones obtained from the classical and first order shear deformation theory (i.e., linear variation of bending stress) but the distribution of transverse shear stresses take the form of a parabolic variation that cannot be presented by the classical and the first order shear deformation theory. Also, Figure 6 shows effects of porosity distribution and length scale parameters. The stress variations are highly affected by the variations of material properties in Figure 1b.

4.3. Free Vibration

The fundamental frequencies of the functionally graded porous plates are obtained by solving Equation (40) for m = 1 and n = 1 .
Figure 7 shows the nondimensional fundamental frequencies versus porosity distributions for various power-law indices without including length scale parameter. The porous plate without material variation, i.e., n = 0 , the natural frequencies of the plate with Type 1 porosity distribution show higher values compared the ones with Type 2 and 3 porosity distributions. Although the total mass of the plates decreases same amount as the porosity increases regardless of the distribution types, Type 2 and 3 result a decrease in natural frequency values as the porosity increases because the stiffness of the plate is significantly degraded in the case of Type 2 and 3. For the FGPM plates with nonzero power-law index, Type 1 and 3 distribution show increase in natural frequency values as the porosity increases. However, Type 2 exhibits the reverse behavior because of the strong contribution coming from the coupling terms which are resulted from the heterogeneity in both porosity distribution and the enhancement in material variation (the density of the plate porosity is very small above the mid-plane section in the case of Type 2 distribution and E t / E b = 10 , ρ t / ρ b = 10 ). The variation of the nondimensional form of fundamental frequencies with varying power-law index and the length scale parameter is presented in Figure 8. The length scale parameter causes significant change in fundamental frequencies of the plates because the stiffening effect in small scale structures is captured by the length scale parameter. Although Type 1 and 3 result higher frequencies compared to the non-porous plate, the frequency values obtained from Type 2 distribution are lower than both non-porous plate and the plate with porosity distribution of Type 1 and 3. As stated before, the reason is that Type 2 is the most effected distribution types by the heterogeneity (with high Elastic modulus and low porosity density at the top sections and reverse at the bottom sections of the plate) and hence this causes high impact of coupling terms which change the behavior significantly.

4.4. Buckling

The critical buckling load is obtained by solving Equation (42) for the functionally graded porous plate subjected to in-plane compressive forces.
Figure 9 shows the variation of non-dimensional critical buckling load with porosity values for FGPM square plate. It is clearly shown that as the power-law index increases the critical buckling load decreases, which makes the plate susceptible to buckling failures. This is also valid for the porosity factor. As expected, as the porosity of the material increases, plate become softer. Also, Type 1 distribution results in stiffer behavior against buckling regardless of the power-law index.
The results of the investigation of the effect of length scale parameter for FGPM are given from Figure 10. To examine the effect of the different porosity distribution types, the porosity value is kept constant as 0.5. As in the case of bending and the free vibration analysis, the stiffness of the plates increase as the length scale parameter increases. On the other hand, the differences in the porosity variation and the enhanced modulus for the top portion of the plate cause deviations and different trends for each porosity distribution.

5. Conclusions

Static bending, free vibration and buckling analysis of the functionally graded porous micro-plates are conducted using general third order plate theory (GTPT) developed by Reddy and Kim [24] and solutions are obtained using Navier’s solution technique. The effect of the material length scale factor and material property variation through the plate thickness direction are included while using the third order plate theory. Three different porosity distributions are used in the analysis and the effects of porosity variations are investigated. In the numerical problems, the elastic modulus and density of the plate have power-law distribution through the thickness while keeping top material is stiffer and denser ( E t / E b = 10 and ρ t / ρ b = 10 ).
The results of the bending, free vibration and the buckling problems show that each distribution has different responses due to effect of coupling terms resulted from the heterogeneity of the material properties and symmetric or/and asymmetric variation of porosity. In general, because of the assumption that E t / E b = 10 , the plate becomes softer as the power-law index increases. Also, the material length scale parameter has a strong effect which makes the plate stiffer. The numerical solutions of this study is limited to rectangular simply supported plate due to the limitation of Navier solution technique.

Author Contributions

Conceptualization, S.C., J.K., and H.T.; methodology, J.K.; software, S.C. and J.K.; validation, S.C. and J.K.; formal analysis, S.C.; investigation, S.C.; resources, J.K. and H.T.; data curation, S.C.; writing—original draft preparation, S.C.; writing—review and editing, J.K. and H.T.; visualization, S.C.; supervision, J.K. and H.T.; project administration, J.K. and H.T.; funding acquisition, H.T.

Funding

This research was funded by Western Michigan University.

Acknowledgments

The authors gratefully acknowledges the support of Western Michigan University through research start-up fund.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bhushan, B. Nanotribology and nanomechanics of MEMS/NEMS and BioMEMS/BioNEMS materials and devices. In Nanotribology and Nanomechanics: An Introduction, 4th ed.; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  2. Fleck, N.A.; Muller, G.M.; Ashby, M.F.; Hutchinson, J.W. Strain gradient plasticity: Theory and experiment. Acta Metall. Mater. 1994, 42, 475–487. [Google Scholar] [CrossRef]
  3. Lam, D.C.; Yang, F.; Chong, A.C.; Wang, J.; Tong, P. Experiments and theory in strain gradient elasticity. J. Mech. Phys. Solids 2003, 51, 1477–1508. [Google Scholar] [CrossRef]
  4. Mcfarland, A.W.; Colton, J.S. Role of material microstructure in plate stiffness with relevance to microcantilever sensors. J. Micromech. Microeng. 2005, 15, 1060–1067. [Google Scholar] [CrossRef]
  5. Motz, C.; Schöberl, T.; Pippan, R. Mechanical properties of micro-sized copper bending beams machined by the focused ion beam technique. Acta Mater. 2005, 53, 4269–4279. [Google Scholar] [CrossRef]
  6. Demir, E.; Raabe, D.; Roters, F. The mechanical size effect as a mean-field breakdown phenomenon: Example of microscale single crystal beam bending. Acta Mater. 2010, 58, 1876–1886. [Google Scholar] [CrossRef]
  7. Ke, L.L.; Wang, Y.S. Size effect on dynamic stability of functionally graded microbeams based on a modified couple stress theory. Compos. Struct. 2011, 93, 342–350. [Google Scholar] [CrossRef]
  8. Liu, D.; He, Y.; Dunstan, D.J.; Zhang, B.; Gan, Z.; Hu, P.; Ding, H. Toward a further understanding of size effects in the torsion of thin metal wires: An experimental and theoretical assessment. Int. J. Plast. 2013, 41, 30–52. [Google Scholar] [CrossRef]
  9. Chen, D.; Yang, J.; Kitipornchai, S. Elastic buckling and static bending of shear deformable functionally graded porous beam. Compos. Struct. 2015, 133, 54–61. [Google Scholar] [CrossRef] [Green Version]
  10. Toupin, R.A. Elastic materials with couple-stresses. Arch. Ration. Mech. Anal. 1962, 11, 385–414. [Google Scholar] [CrossRef]
  11. Mindlin, R.D.; Tiersten, H.F. Effects of couple-stresses in linear elasticity. Arch. Ration. Mech. Anal. 1962, 11, 415–448. [Google Scholar] [CrossRef]
  12. Yang, F.; Chong, A.C.; Lam, D.C.; Tong, P. Couple stress based strain gradient theory for elasticity. Int. J. Solids Struct. 2002, 39, 2731–2743. [Google Scholar] [CrossRef]
  13. Eringen, A.C. Nonlocal polar elastic continua. Int. J. Eng. Sci. 1972, 10, 1–16. [Google Scholar] [CrossRef]
  14. Fleck, N.A.; Hutchinson, J.W. A reformulation of strain gradient plasticity. J. Mech. Phys. Solids 2001, 49, 2245–2271. [Google Scholar] [CrossRef] [Green Version]
  15. Yang, Y.Y. Stress analysis in a joint with a functionally graded material under a thermal loading by using the Mellin transform method. Int. J. Solids Struct. 1998, 35, 1261–1287. [Google Scholar] [CrossRef]
  16. Neubrand, A.; Chung, T.J.; Rödel, J.; Steffler, E.D.; Fett, T. Residual stresses in functionally graded plates. J. Mater. Res. 2002, 17, 2912–2920. [Google Scholar] [CrossRef]
  17. Birman, V.; Byrd, L.W. Modeling and Analysis of Functionally Graded Materials and Structures. Appl. Mech. Rev. 2007, 60, 195. [Google Scholar] [CrossRef]
  18. Şimşek, M.; Reddy, J.N. A unified higher order beam theory for buckling of a functionally graded microbeam embedded in elastic medium using modified couple stress theory. Compos. Struct. 2013, 101, 47–58. [Google Scholar] [CrossRef]
  19. Reddy, J.N. Microstructure-dependent couple stress theories of functionally graded beams. J. Mech. Phys. Solids 2011, 59, 2382–2399. [Google Scholar] [CrossRef]
  20. Li, L.; Hu, Y. Nonlinear bending and free vibration analyses of nonlocal strain gradient beams made of functionally graded material. Int. J. Eng. Sci. 2016, 107, 77–97. [Google Scholar] [CrossRef]
  21. Javaheri, R.; Eslami, M.R. Buckling of functionally graded plates under in-plane compressive loading. J. Appl. Math. Mech. 2002, 82, 277–283. [Google Scholar] [CrossRef]
  22. Ma, L.S.; Wang, T.J. Nonlinear bending and post-buckling of a functionally graded circular plate under mechanical and thermal loadings. Int. J. Solids Struct. 2003, 40, 3311–3330. [Google Scholar] [CrossRef]
  23. Lanhe, W. Thermal buckling of a simply supported moderately thick rectangular FGM plate. Compos. Struct. 2004, 64, 211–218. [Google Scholar] [CrossRef]
  24. Reddy, J.; Kim, J. A nonlinear modified couple stress-based third-order theory of functionally graded plates. Compos. Struct. 2012, 94, 1128–1143. [Google Scholar] [CrossRef]
  25. Kim, J.; Reddy, J.N. Analytical solutions for bending, vibration, and buckling of FGM plates using a couple stress-based third-order theory. Compos. Struct. 2013, 103, 86–98. [Google Scholar] [CrossRef]
  26. Ansari, R.; Faghih Shojaei, M.; Mohammadi, V.; Gholami, R.; Darabi, M.A. Nonlinear vibrations of functionally graded Mindlin microplates based on the modified couple stress theory. Compos. Struct. 2014, 114, 124–134. [Google Scholar] [CrossRef]
  27. Zhang, B.; He, Y.; Liu, D.; Shen, L.; Lei, J. An efficient size-dependent plate theory for bending, buckling and free vibration analyses of functionally graded microplates resting on elastic foundation. Appl. Math. Model. 2015, 39, 3814–3845. [Google Scholar] [CrossRef]
  28. Al Jahwari, F.; Anwer, A.A.; Naguib, H.E. Fabrication and microstructural characterization of functionally graded porous acrylonitrile butadiene styrene and the effect of cellular morphology on creep behavior. J. Polym. Sci. Part B Polym. Phys. 2015, 53, 795–803. [Google Scholar] [CrossRef]
  29. Al Jahwari, F.; Huang, Y.; Naguib, H.E.; Lo, J. Relation of impact strength to the microstructure of functionally graded porous structures of acrylonitrile butadiene styrene (ABS) foamed by thermally activated microspheres. Polymer 2016, 98, 270–281. [Google Scholar] [CrossRef]
  30. Yadroitsev, I.; Shishkovsky, I.; Bertrand, P.; Smurov, I. Manufacturing of fine-structured 3D porous filter elements by selective laser melting. Appl. Surf. Sci. 2009, 255, 5523–5527. [Google Scholar] [CrossRef]
  31. Koohbor, B.; Kidane, A. Design optimization of continuously and discretely graded foam materials for efficient energy absorption. Mater. Des. 2016, 102, 151–161. [Google Scholar] [CrossRef]
  32. Giannitelli, S.M.; Basoli, F.; Mozetic, P.; Piva, P.; Bartuli, F.N.; Luciani, F.; Arcuri, C.; Trombetta, M.; Rainer, A.; Licoccia, S. Graded porous polyurethane foam: A potential scaffold for oro-maxillary bone regeneration. Mater. Sci. Eng. C 2015, 51, 329–335. [Google Scholar] [CrossRef] [PubMed]
  33. Han, C.; Li, Y.; Wang, Q.; Wen, S.; Wei, Q.; Yan, C.; Hao, L.; Liu, J.; Shi, Y. Continuous functionally graded porous titanium scaffolds manufactured by selective laser melting for bone implants. J. Mech. Behav. Biomed. Mater. 2018, 80, 119–127. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, D.; Yang, J.; Kitipornchai, S. Free and forced vibrations of shear deformable functionally graded porous beams. Int. J. Mech. Sci. 2016, 108, 14–22. [Google Scholar] [CrossRef] [Green Version]
  35. Shafiei, N.; Kazemi, M. Buckling analysis on the bi-dimensional functionally graded porous tapered nano-/micro-scale beams. Aerosp. Sci. Technol. 2017, 66, 1–11. [Google Scholar] [CrossRef]
  36. Zhao, J.; Wang, Q.; Deng, X.; Choe, K.; Xie, F.; Shuai, C. A modified series solution for free vibration analyses of moderately thick functionally graded porous (FGP) deep curved and straight beams. Compos. Part B Eng. 2018, 165, 155–166. [Google Scholar] [CrossRef]
  37. Wu, D.; Liu, A.; Huang, Y.; Huang, Y.; Pi, Y.; Gao, W. Dynamic analysis of functionally graded porous structures through finite element analysis. Eng. Struct. 2018, 165, 287–301. [Google Scholar] [CrossRef]
  38. Shahverdi, H.; Barati, M.R. Vibration analysis of porous functionally graded nanoplates. Int. J. Eng. Sci. 2017, 120, 82–99. [Google Scholar] [CrossRef]
  39. Wang, Y.Q.; Wan, Y.H.; Zhang, Y.F. Vibrations of longitudinally traveling functionally graded material plates with porosities. Eur. J. Mech. A Solids 2017, 66, 55–68. [Google Scholar] [CrossRef]
  40. Rezaei, A.S.; Saidi, A.R.; Abrishamdari, M.; Mohammadi, M.H. Natural frequencies of functionally graded plates with porosities via a simple four variable plate theory: An analytical approach. Thin-Walled Struct. 2017, 120, 366–377. [Google Scholar] [CrossRef]
  41. Thai, H.T.; Choi, D.H. A simple first-order shear deformation theory for the bending and free vibration analysis of functionally graded plates. Compos. Struct. 2013, 101, 332–340. [Google Scholar] [CrossRef]
  42. Akbaş, S.D. Vibration and static analysis of functionally graded porous plates. J. Appl. Comput. Mech. 2017, 3, 199–207. [Google Scholar]
  43. Wang, Y.Q.; Zu, J.W. Large-amplitude vibration of sigmoid functionally graded thin plates with porosities. Thin-Walled Struct. 2017, 119, 911–924. [Google Scholar] [CrossRef]
  44. Cong, P.H.; Chien, T.M.; Khoa, N.D.; Duc, N.D. Nonlinear thermomechanical buckling and post-buckling response of porous FGM plates using Reddy’s HSDT. Aerosp. Sci. Technol. 2018, 77, 419–428. [Google Scholar] [CrossRef]
  45. Mirjavadi, S.S.; Afshari, B.M.; Barati, M.R.; Hamouda, A. Transient response of porous FG nanoplates subjected to various pulse loads based on nonlocal stress-strain gradient theory. Eur. J. Mech. A Solids 2018, 74, 210–220. [Google Scholar] [CrossRef]
  46. Zhao, J.; Choe, K.; Xie, F.; Wang, A.; Shuai, C.; Wang, Q. Three-dimensional exact solution for vibration analysis of thick functionally graded porous (FGP) rectangular plates with arbitrary boundary conditions. Compos. Part B Eng. 2018, 155, 369–381. [Google Scholar] [CrossRef]
  47. Zhao, J.; Xie, F.; Wang, A.; Shuai, C.; Tang, J.; Wang, Q. Dynamics analysis of functionally graded porous (FGP) circular, annular and sector plates with general elastic restraints. Compos. Part B Eng. 2019, 159, 20–43. [Google Scholar] [CrossRef]
  48. Thang, P.T.; Nguyen-Thoi, T.; Lee, D.; Kang, J.; Lee, J. Elastic buckling and free vibration analyses of porous-cellular plates with uniform and non-uniform porosity distributions. Aerosp. Sci. Technol. 2018, 79, 278–287. [Google Scholar] [CrossRef]
  49. Kim, J.; Żur, K.K.; Reddy, J.N. Bending, free vibration, and buckling of modified couples stress-based functionally graded porous micro-plates. Compos. Struct. 2019, 209, 879–888. [Google Scholar] [CrossRef]
  50. Reddy, J.N. Theory and Analysis of Elastic Plates and Shells, 2nd ed.; CRC Press Taylor & Francis Group: Boca Raton, FL, USA, 2006. [Google Scholar]
  51. Reddy, J.N. Energy Principles and Variational Methods in Applied Mechanics; John Wiley & Sons: Hoboken, NJ, USA, 2017. [Google Scholar]
Figure 1. The normalized porosity distribution through the plate thickness (a) and normalized material property resulted from Type 1 distribution (b); Type 2 distribution (c); and Type 3 distribution (d).
Figure 1. The normalized porosity distribution through the plate thickness (a) and normalized material property resulted from Type 1 distribution (b); Type 2 distribution (c); and Type 3 distribution (d).
Jcs 03 00015 g001
Figure 2. Boundary conditions of the simply supported plate [50] (where i = 0, 1, 2, 3 and j = 0, 1, 2) (a) and the cross-section of the FGPM plate (Type 3) (b).
Figure 2. Boundary conditions of the simply supported plate [50] (where i = 0, 1, 2, 3 and j = 0, 1, 2) (a) and the cross-section of the FGPM plate (Type 3) (b).
Jcs 03 00015 g002
Figure 3. Comparisons of the non-dimensional deflection, w ¯ x , b 2 , 0 (a), the non-dimensional fundamental frequency, ω ¯ (b), and the non-dimensional critical buckling load, N 0 ¯ (c).
Figure 3. Comparisons of the non-dimensional deflection, w ¯ x , b 2 , 0 (a), the non-dimensional fundamental frequency, ω ¯ (b), and the non-dimensional critical buckling load, N 0 ¯ (c).
Jcs 03 00015 g003
Figure 4. Non-dimensional deflection w ¯ x , b 2 , 0 of homogeneous plate with different porosity distributions (a) and FGM plate with various power-law index, n (b).
Figure 4. Non-dimensional deflection w ¯ x , b 2 , 0 of homogeneous plate with different porosity distributions (a) and FGM plate with various power-law index, n (b).
Jcs 03 00015 g004
Figure 5. Non-dimensional central deflection w ¯ a 2 , b 2 , 0 versus length scale parameter l / h along a FGM simply supported plate with constant porosity, ϕ = 0.5 and various power law index n.
Figure 5. Non-dimensional central deflection w ¯ a 2 , b 2 , 0 versus length scale parameter l / h along a FGM simply supported plate with constant porosity, ϕ = 0.5 and various power law index n.
Jcs 03 00015 g005
Figure 6. Bending stresses ( σ ¯ y y ) distribution (a) and transverse shear stresses ( σ ¯ x z ) distribution (b) through plate thickness.
Figure 6. Bending stresses ( σ ¯ y y ) distribution (a) and transverse shear stresses ( σ ¯ x z ) distribution (b) through plate thickness.
Jcs 03 00015 g006
Figure 7. Non-dimensional fundamental frequency, ω ¯ versus porosity values, ϕ with zero length scale parameter l / h = 0 .
Figure 7. Non-dimensional fundamental frequency, ω ¯ versus porosity values, ϕ with zero length scale parameter l / h = 0 .
Jcs 03 00015 g007
Figure 8. Non-dimensional fundamental frequency, ω ¯ versus power-law index, n.
Figure 8. Non-dimensional fundamental frequency, ω ¯ versus power-law index, n.
Jcs 03 00015 g008
Figure 9. Non-dimensional critical buckling load, N 0 ¯ versus porosity values, ϕ with zero length scale parameter l / h = 0 .
Figure 9. Non-dimensional critical buckling load, N 0 ¯ versus porosity values, ϕ with zero length scale parameter l / h = 0 .
Jcs 03 00015 g009
Figure 10. Non-dimensional critical buckling load, N 0 ¯ versus power-law index, n.
Figure 10. Non-dimensional critical buckling load, N 0 ¯ versus power-law index, n.
Jcs 03 00015 g010

Share and Cite

MDPI and ACS Style

Coskun, S.; Kim, J.; Toutanji, H. Bending, Free Vibration, and Buckling Analysis of Functionally Graded Porous Micro-Plates Using a General Third-Order Plate Theory. J. Compos. Sci. 2019, 3, 15. https://doi.org/10.3390/jcs3010015

AMA Style

Coskun S, Kim J, Toutanji H. Bending, Free Vibration, and Buckling Analysis of Functionally Graded Porous Micro-Plates Using a General Third-Order Plate Theory. Journal of Composites Science. 2019; 3(1):15. https://doi.org/10.3390/jcs3010015

Chicago/Turabian Style

Coskun, Semsi, Jinseok Kim, and Houssam Toutanji. 2019. "Bending, Free Vibration, and Buckling Analysis of Functionally Graded Porous Micro-Plates Using a General Third-Order Plate Theory" Journal of Composites Science 3, no. 1: 15. https://doi.org/10.3390/jcs3010015

APA Style

Coskun, S., Kim, J., & Toutanji, H. (2019). Bending, Free Vibration, and Buckling Analysis of Functionally Graded Porous Micro-Plates Using a General Third-Order Plate Theory. Journal of Composites Science, 3(1), 15. https://doi.org/10.3390/jcs3010015

Article Metrics

Back to TopTop