This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Hindawi Publishing CorporationJournal of CompositesVolume 2013, Article ID 808764, 12 pageshttp://dx.doi.org/10.1155/2013/808764
Research ArticleBuckling Analysis of Functionally Graded Material Plates UsingHigher Order Shear Deformation Theory
B. Sidda Reddy,1 J. Suresh Kumar,2 C. Eswara Reddy,3 and K. Vijaya Kumar Reddy2
1 School of Mechanical Engineering, R.G.M College of Engineering & Technology, Nandyal, Kurnool, Andhra Pradesh 518 501, India2Department of Mechanical Engineering, J.N.T.U.H College of Engineering, J.N.T. University, Hyderabad,Andhra Pradesh 500 085, India
3The School of Engineering & Technology, Sri Padmavathi Mahila Visvavidyalayam, Women’s University, Tirupati, Chittoor,Andhra Pradesh 517 502, India
Correspondence should be addressed to B. Sidda Reddy; [email protected]
Received 29 May 2013; Revised 7 October 2013; Accepted 11 October 2013
The prime aim of the present study is to present analytical formulations and solutions for the buckling analysis of simply supportedfunctionally graded plates (FGPs) using higher order shear deformation theory (HSDT) without enforcing zero transverse shearstresses on the top and bottom surfaces of the plate. It does not require shear correction factors and transverse shear stresses varyparabolically across the thickness. Material properties of the plate are assumed to vary in the thickness direction according to apower law distribution in terms of the volume fractions of the constituents. The equations of motion and boundary conditions arederived using the principle of virtual work. Solutions are obtained for FGPs in closed-form using Navier’s technique. Comparisonstudies are performed to verify the validity of the present results fromwhich it can be concluded that the proposed theory is accurateand efficient in predicting the buckling behavior of functionally graded plates. The effect of side-to-thickness ratio, aspect ratio,modulus ratio, the volume fraction exponent, and the loading conditions on the critical buckling load of FGPs is also investigatedand discussed.
1. Introduction
Functionally graded materials (FGMs) are the new gener-ation of novel composite materials in the family of engi-neering composites, whose properties are varied smoothlyin the spatial direction microscopically to improve theoverall structural performance. These materials offer greatpromise in high temperature environments, for example,wear-resistant linings for handling large heavy abrasive oreparticles, rocket heat shields, heat exchanger tubes, thermo-electric generators, heat engine components, plasma facingsfor fusion reactors, and electrically insulating metal/ceramicjoints and also these are widely used in many structuralapplications such as mechanics, civil engineering, optical,electronic, chemical,mechanical, biomedical, energy sources,nuclear, automotive fields, and ship building industries tominimize thermomechanical mismatch in metal-ceramicbonding. Most structures, irrespective of their use, will be
subjected to dynamic loads during their operational life.Increased use of FGMs in various structural applicationsnecessitates the development of accurate theoretical modelsto predict their response.
In the past, a variety of plate theories have been proposedto study the buckling behavior of FGM plates. The classicalplate theory (CPT) provides acceptable results only for theanalysis of thin plates and neglects the transverse sheareffects. Javaheri and Eslami [1], Abrate [2], Mohammadi etal. [3], Mahdavian [4], Feldman and Aboudi [5], Shariatet al. [6], and Tung and Duc [7] employed this theoryto analyze buckling behavior of FG plates. However, formoderately thick plates CPT underpredicts deflections andoverpredicts buckling loads and natural frequencies. Thefirst-order shear deformation theories (FSDTs) are based onReissner [8] andMindlin [9] accounts for the transverse sheardeformation effect by means of a linear variation of inplanedisplacements and stresses through the thickness of the plate,
2 Journal of Composites
but requires a correction factor to satisfy the free transverseshear stress conditions on the top and bottom surfacesof the plate. Although the FSDT provides a sufficientlyaccurate description of response for thin to moderatelythick plates, it is not convenient to use due to difficultywith determination of the correct value of shear correctionfactor [10]. The authors [11–16] used FSDT to analyze thebuckling of FG plates. In order to overcome the limitationsof FSDT many HSDTs were developed that involve higherorder terms in Taylors expansions of the displacements in thethickness coordinate. Javaheri and Eslami [17], Najafizadehand Heydari [18], Bodaghi and Saidi [19], Bagherizadeh etal. [20], and Mozafari and Ayob [21] used the HSDT toanalyze the buckling behavior of FG plates. Ma and Wang[22] have investigated the axisymmetric large deflectionbending and postbuckling behavior of a functionally gradedcircular plate under mechanical, thermal, and combinedthermal-mechanical load based on classical nonlinear vonKarman plate theory. They observed from their investigationthat the power law index “𝑛” has a significant effect onthe midplane temperature, critical buckling temperature,and on the thermal post-buckling behavior of FGM plate.Hosseini-Hashemi et al. [23] have developed the closed-form solutions in analytical form to study the bucklingbehavior of in-plane loaded isotropic rectangular FG plateswithout any use of approximation for different boundaryconditions using the Mindlin plate theory. Saidi et al. [24]employed the unconstrained third-order shear deformationtheory to analyze the axisymmetric bending and bucklingof FG solid circular plates in which the bending-stretchingcoupling exists. Oyekoya et al. [25] developed Mindlin typeand Ressner type element formodeling of FG composite platesubjected to buckling and free vibration. Further, they studiedthe plate for the effect of different fiber distribution cases andthe effects of fire distribution on buckling, and free vibration.Ghannadpour et al. [26] applied finite stripmethod to analyzethe buckling behavior of rectangular FG plats under thermalload. The solution was obtained by the minimization ofthe total potential energy and solving the correspondingeigenvalue problem. Thai and Choi [27] presented a simplerefined theory to analyze the buckling behavior of FG plateswhich has strong similarity with classical plate theory inmany aspects, accounts for a quadratic variation of thetransverse shear strains across the thickness, and satisfies thezero traction boundary conditions on the top and bottomsurfaces of the plate without using shear correction factors.The governing equations were derived from the principleof minimum total potential energy. The effects of loadingconditions and variations of power of functionally gradedmaterial, modulus ratio, aspect ratio, and thickness ratiowere also investigated by these authors. Thai and Vo [10]have developed a new sinusoidal shear deformation theoryto study the bending, buckling and vibration of FG platesaccounting for sinusoidal distribution of transverse shearstress and satisfies the free transverse shear stress conditionson the top and bottom surfaces of the plate without usingshear correction factor. Uymaz and Aydogdu [28] analyzedthe rectangular FG plates under different axial loadings forbuckling based on small strain elasticity theory with different
boundary conditions.They also investigated the effects of thedifferent material composition and the plate geometry on thecritical buckling loads and mode shapes.
Lal et al. [29] have examined the second order statis-tics of postbuckling responses of FGM plate subjected tomechanical and thermal load with nonuniform temperaturechanges subjected to temperature independent and depen-dent material properties. The effect of random materialproperties with amplitude ratios, volume fraction index, platethickness ratios, aspect ratios, boundary conditions, andtypes of loadings subjected to temperature independent andtemperature dependent material properties were investigatedthrough numerical examples.
This paper aims to develop analytical formulations andsolutions for the buckling analysis of functionally gradedplates (FGPs) using higher order shear deformation theory(HSDT) without enforcing zero transverse shear stress on thetop and bottom surfaces of the plate. This does not requireshear correction factor. The plate material is graded throughthe thickness direction. The plate’s governing equationsand its boundary conditions are derived by employing theprinciple of virtual work. Solutions are obtained for FGPsin closed-form using Navier’s technique and solving theeigenvalue equation. The present results are compared withthe solutions of Thai and Choi [27] to verify the accuracy ofthe proposed theory in predicting the critical buckling loadsof FG plates. The effect of side-to-thickness ratios, aspectratios, and modulus ratios and the volume fraction exponenton the critical buckling loads are studied after establishing theaccuracy of the present results for FG plates.
2. Theoretical Formulation
In formulating the higher order shear deformation theory,a rectangular plate of length 𝑎, width 𝑏, and thickness ℎ isconsidered, which composed of functionally graded materialthrough the thickness. Figure 1 shows the functionally gradedmaterial plate with the rectangular Cartesian coordinatesystem 𝑥, 𝑦, and 𝑧. Thematerial properties are assumed to bevaried in the thickness direction only and the bright and darkareas correspond to ceramic andmetal particles, respectively.On the top surface (𝑧 = +ℎ/2), the plate is composed of fullceramic and graded to the bottom surface (𝑧 = −ℎ/2) whichcomposed of full metal. The reference surface is the middlesurface of the plate (𝑧 = 0). The functionally graded materialplate properties are assumed to be the function of the volumefraction of constituent materials. The functional relationshipbetween the material property and the thickness coordinatesis assumed to be
𝑃 (𝑧) = (𝑃𝑡− 𝑃𝑏) (
𝑧
ℎ
+
1
2
)
𝑛
+ 𝑃𝑏, (1)
where 𝑃 denotes the effective material property, 𝑃𝑡and
𝑃𝑏denote the property on the top and bottom surface
of the plate, respectively, and 𝑛 is the material variationparameter that dictates the material variation profile throughthe thickness. The effective material properties of the plate,including Young’s modulus, 𝐸, density, 𝜌, and shear modulus,
Journal of Composites 3
a
bz = h/2
z = −h/2
Z
X
Y
Figure 1: Functionally graded plate and coordinates.
𝐺, vary according to (1), and poisons ratio (𝜐) is assumed tobe constant.
2.1. Displacement Models. In order to approximate 3Dplate problem to a 2D one, the displacement components𝑢 (𝑥, 𝑦, 𝑧, 𝑡), V (𝑥, 𝑦, 𝑧, 𝑡), and 𝑤 (𝑥, 𝑦, 𝑧, 𝑡) at any point inthe plate are expanded in terms of the thickness coordinate.The elasticity solution indicates that the transverse shearstress varies parabolically through the plate thickness. Thisrequires the use of a displacement field, in which the in-plane displacements are expanded as cubic functions ofthe thickness coordinate. In addition, the transverse normalstrain may vary nonlinearly through the plate thickness. Thedisplacement field which satisfies the above criteria may beassumed in the form:
𝑢 (𝑥, 𝑦, 𝑧) = 𝑢𝑜(𝑥, 𝑦) + 𝑧𝜃
𝑥(𝑥, 𝑦)
+ 𝑧2𝑢∗
𝑜(𝑥, 𝑦) + 𝑧
3𝜃∗
𝑥(𝑥, 𝑦)
V (𝑥, 𝑦, 𝑧) = V𝑜(𝑥, 𝑦) + 𝑧𝜃
𝑦(𝑥, 𝑦)
+ 𝑧2V∗𝑜(𝑥, 𝑦) + 𝑧
3𝜃∗
𝑦(𝑥, 𝑦)
𝑤 (𝑥, 𝑦, 𝑧) = 𝑤𝑜(𝑥, 𝑦) ,
(2)
where𝑢𝑜, V𝑜, and𝑤
𝑜denote the displacements of a point (𝑥, 𝑦)
on the midplane.𝜃𝑥, 𝜃𝑦are rotations of the normal to the mid plane about
𝑦 and 𝑥-axes.𝑢∗
0, V∗0, 𝜃∗𝑥, and 𝜃∗
𝑦are the higher order deformation terms
defined at the mid plane.By substitution of displacement relations from (2) into
the strain displacement equations of the classical theory ofelasticity, the following relations are obtained:
𝜀𝑥= 𝜀𝑥𝑜
+ 𝑧𝑘𝑥+ 𝑧2𝜀∗
𝑥𝑜+ 𝑧3𝑘∗
𝑥,
𝜀𝑦= 𝜀𝑦𝑜
+ 𝑧𝑘𝑦+ 𝑧2𝜀∗
𝑦𝑜+ 𝑧3𝑘∗
𝑦,
𝜀𝑧= 0,
𝛾𝑥𝑦
= 𝜀𝑥𝑦𝑜
+ 𝑧𝑘𝑥𝑦
+ 𝑧2𝜀∗
𝑥𝑦𝑜+ 𝑧3𝑘∗
𝑥𝑦,
𝛾𝑦𝑧
= 𝜑𝑦+ 𝑧𝜀𝑦𝑧𝑜
+ 𝑧2𝜑∗
𝑦,
𝛾𝑥𝑧
= 𝜑𝑥+ 𝑧𝜀𝑥𝑧𝑜
+ 𝑧2𝜑∗
𝑥,
(3)
where
𝜀𝑥𝑜
=
𝜕𝑢𝑜
𝜕𝑥
, 𝜀𝑦𝑜
=
𝜕V𝑜
𝜕𝑦
, 𝜀𝑥𝑦𝑜
=
𝜕𝑢𝑜
𝜕𝑦
+
𝜕V𝑜
𝜕𝑥
𝑘𝑥=
𝜕𝜃𝑥
𝜕𝑥
, 𝑘𝑦=
𝜕𝜃𝑦
𝜕𝑦
, 𝑘𝑥𝑦
=
𝜕𝜃𝑥
𝜕𝑦
+
𝜕𝜃𝑦
𝜕𝑥
𝑘∗
𝑥=
𝜕𝜃∗
𝑥
𝜕𝑥
, 𝑘∗
𝑦=
𝜕𝜃∗
𝑦
𝜕𝑦
, 𝑘∗
𝑥𝑦=
𝜕𝜃∗
𝑥
𝜕𝑦
+
𝜕𝜃∗
𝑦
𝜕𝑥
𝜀∗
𝑥𝑜=
𝜕𝑢∗
𝑜
𝜕𝑥
, 𝜀∗
𝑦𝑜=
𝜕V∗𝑜
𝜕𝑦
, 𝜀∗
𝑥𝑦𝑜=
𝜕𝑢∗
𝑜
𝜕𝑦
+
𝜕V∗𝑜
𝜕𝑥
𝜑𝑦= 𝜃𝑦+
𝜕𝑤𝑜
𝜕𝑦
, 𝜑𝑥= 𝜃𝑥+
𝜕𝑤𝑜
𝜕𝑥
𝜀𝑦𝑧𝑜
= 2V∗𝑜, 𝜀
𝑥𝑧𝑜= 2𝑢∗
𝑜, 𝜑
∗
𝑦= 3𝜃∗
𝑦, 𝜑
∗
𝑥= 3𝜃∗
𝑥.
(4)
2.2. Elastic Stress-Strain Relations. The elastic stress-strainrelations dependonwhich assumption of 𝜀
𝑧= 0. In the case of
functionally graded materials the constitutive equations canbe written as
{{{{{
{{{{{
{
𝜎𝑥
𝜎𝑦
𝜏𝑥𝑦
𝜏𝑦𝑧
𝜏𝑥𝑧
}}}}}
}}}}}
}
=
[
[
[
[
[
[
𝑄11
𝑄12
0 0 0
𝑄12
𝑄22
0 0 0
0 0 𝑄33
0 0
0 0 0 𝑄44
0
0 0 0 0 𝑄55
]
]
]
]
]
]
[
[
[
[
[
[
𝜀𝑥
𝜀𝑦
𝛾𝑥𝑦
𝛾𝑦𝑧
𝛾𝑥𝑧
]
]
]
]
]
]
, (5)
where 𝜎 = (𝜎𝑥, 𝜎𝑦, 𝜏𝑥𝑦, 𝜏𝑦𝑧, 𝜏𝑥𝑧)𝑡 are the stresses, 𝜀 =
(𝜀𝑥, 𝜀𝑦, 𝛾𝑥𝑦, 𝛾𝑦𝑧, 𝛾𝑥𝑧)𝑡 are the strains with respect to the axes,
and𝑄𝑖𝑗’s are the plane stress reduced elastic coefficients in the
plate axes that vary through the plate thickness given by
𝑄11
= 𝑄22
=
𝐸 (𝑍)
1 − 𝜐2
=
(𝐸𝑐− 𝐸𝑚) ((𝑧/ℎ) + (1/2))
𝑛+ 𝐸𝑚
1 − 𝜐2
,
𝑄12
= 𝑄21
= 𝜐𝑄11,
𝑄33
= 𝑄44
= 𝑄55
=
(1 − 𝜐2)
2 (1 + 𝜐)
𝑄11,
(6)
where 𝐸𝑐is the modulus of elasticity of the ceramic material
and 𝐸𝑚is the modulus of elasticity of the metal.
2.3. Governing Equations of Motion. The work done by theactual forces inmoving through virtual displacements, whichare consistent with the geometric constraints of a body, is setto zero to obtain the equation of motion and this is knownas energy principle. It is useful in (a) deriving governingequations and the boundary conditions and (b) obtainingapproximate solutions by virtual methods.
Energy principles provide alternative means to obtain thegoverning equations and their solutions. In the present study,
4 Journal of Composites
the principle of virtual work is used to derive the equations ofmotion of functionally graded plates.
The governing equations of displacement model in (2)will be derived using the dynamic version of the principle ofvirtual displacements; that is,
The virtual strain energy, work done, and kinetic energyis given by
𝛿𝑈 =∫
𝐴
{∫
ℎ/2
−ℎ/2
[𝜎𝑥𝛿∈𝑥+ 𝜎𝑦𝛿∈𝑦+ 𝜏𝑥𝑦𝛿𝛾𝑥𝑦
+𝜏𝑥𝑧𝛿𝛾𝑥𝑧
+ 𝜏𝑦𝑧𝛿𝛾𝑦𝑧] 𝑑𝑧}𝑑𝑥𝑑𝑦,
𝛿𝑉 = − ∫(𝑞𝑤0+ 𝑁𝑥
𝜕𝑤𝑜
𝜕𝑥
𝛿𝜕𝑤𝑜
𝜕𝑥
+ 𝑁𝑥𝑦
𝜕𝑤𝑜
𝜕𝑦
𝛿𝜕𝑤𝑜
𝜕𝑥
+𝑁𝑦𝑥
𝜕𝑤𝑜
𝜕𝑥
𝛿𝜕𝑤𝑜
𝜕𝑦
+ 𝑁𝑦
𝜕𝑤𝑜
𝜕𝑦
𝛿𝜕𝑤𝑜
𝜕𝑦
)𝑑𝑥𝑑𝑦,
𝛿𝐾 =∫
𝐴
{∫
ℎ/2
−ℎ/2
𝜌0[(��0+ 𝑧
𝜃𝑥+ 𝑧2��∗
0+ 𝑧3 𝜃∗
𝑥)
× (𝛿��0+ 𝑧𝛿
𝜃𝑥+ 𝑧2𝛿��
∗
0+ 𝑧3𝛿
𝜃∗
𝑥)
+ (V0+ 𝑧
𝜃𝑦+ 𝑧2V∗0+ 𝑧3 𝜃
∗
𝑦)
× (𝛿V0+ 𝑧𝛿
𝜃𝑦+ 𝑧2𝛿V∗
0+ 𝑧3𝛿
𝜃
∗
𝑦)
+��0𝛿��0] 𝑑𝑧}𝑑𝑥𝑑𝑦,
(8)
where 𝑞 = distributed load over the surface of the plate.𝑁𝑥and 𝑁
𝑦the inplane loads perpendicular to the edges
𝑥 = 0 and 𝑦 = 0, respectively,𝑁𝑥𝑦
𝑁𝑦𝑥
the distributed shearforces parallel to the edges 𝑥 = 0 and 𝑦 = 0, respectively,𝜌0
= density of plate material, ��0= 𝜕𝑢0/𝜕𝑡, V
0= 𝜕V0/𝜕𝑡,
and so forth, indicate the time derivatives.Substituting for 𝛿𝑈, 𝛿𝑉, and 𝛿𝐾 in the virtual work
statement in (7) and integrating through the thickness,integrating by parts, and collecting the coefficients of𝛿𝑢𝑜, 𝛿V𝑜, 𝛿𝑤𝑜, 𝛿𝜃𝑥, 𝛿𝜃𝑦, 𝛿𝑢∗
𝑜, 𝛿V∗𝑜, 𝛿𝜃∗
𝑥, 𝛿𝜃∗
𝑦, the follow-
ing equations of motion are obtained:
𝛿𝑢0:
𝜕𝑁𝑥
𝜕𝑥
+
𝜕𝑁𝑥𝑦
𝜕𝑦
= 𝐼1��0+ 𝐼2
𝜃𝑥+ 𝐼3��∗
0+ 𝐼4
𝜃
∗
𝑥,
𝛿V0:
𝜕𝑁𝑦
𝜕𝑦
+
𝜕𝑁𝑥𝑦
𝜕𝑥
= 𝐼1V0+ 𝐼2
𝜃𝑦+ 𝐼3V∗0+ 𝐼4
𝜃
∗
𝑦,
𝛿𝑤0:
𝜕𝑄𝑥
𝜕𝑥
+
𝜕𝑄𝑦
𝜕𝑦
+ 𝑞 + �� = 𝐼1��0,
(9)
where
�� = 𝑁𝑥
𝜕2𝑤𝑜
𝜕𝑥2
+ 𝑁𝑥𝑦
𝜕2𝑤𝑜
𝜕𝑦𝜕𝑥
+ 𝑁𝑦𝑥
𝜕2𝑤𝑜
𝜕𝑥𝜕𝑦
+ 𝑁𝑦
𝜕2𝑤𝑜
𝜕𝑦2,
𝛿𝜃𝑥:
𝜕𝑀𝑥
𝜕𝑥
+
𝜕𝑀𝑥𝑦
𝜕𝑦
− 𝑄𝑥
= 𝐼2��0+ 𝐼3
𝜃𝑥+ 𝐼4��∗
0+ 𝐼5
𝜃
∗
𝑥,
𝛿𝜃𝑦:
𝜕𝑀𝑦
𝜕𝑦
+
𝜕𝑀𝑥𝑦
𝜕𝑥
− 𝑄𝑦
= 𝐼2V0+ 𝐼3
𝜃𝑦+ 𝐼4V∗0+ 𝐼5
𝜃
∗
𝑦,
𝛿𝑢∗
0:
𝜕𝑁∗
𝑥
𝜕𝑥
+
𝜕𝑁∗
𝑥𝑦
𝜕𝑦
− 2𝑆𝑥
= 𝐼3��0+ 𝐼4
𝜃𝑥+ 𝐼5��∗
0+ 𝐼6
𝜃
∗
𝑥,
𝛿V∗0:
𝜕𝑁∗
𝑦
𝜕𝑦
+
𝜕𝑁∗
𝑥𝑦
𝜕𝑥
− 2𝑆𝑦
= 𝐼3V0+ 𝐼4
𝜃𝑦+ 𝐼5V∗0+ 𝐼6
𝜃
∗
𝑦,
𝛿𝜃∗
𝑥:
𝜕𝑀∗
𝑥
𝜕𝑥
+
𝜕𝑀∗
𝑥𝑦
𝜕𝑦
− 3𝑄∗
𝑥
= 𝐼4��0+ 𝐼5
𝜃𝑥+ 𝐼6��∗
0+ 𝐼7
𝜃
∗
𝑥,
𝛿𝜃∗
𝑦:
𝜕𝑀∗
𝑦
𝜕𝑦
+
𝜕𝑀∗
𝑥𝑦
𝜕𝑥
− 3𝑄𝑦
= 𝐼4V0+ 𝐼5
𝜃𝑦+ 𝐼6V∗0+ 𝐼7
𝜃
∗
𝑦,
(10)
where the force and moment resultants are defined as
{
{
{
𝑁𝑥
| 𝑁∗
𝑥
𝑁𝑦
| 𝑁∗
𝑦
𝑁𝑥𝑦
| 𝑁∗
𝑥𝑦
}
}
}
=
𝑛
∑
𝐿=1
∫
ℎ/2
−ℎ/2
{
{
{
𝜎𝑥
𝜎𝑦
𝜏𝑥𝑦
}
}
}
[1 | 𝑧2] 𝑑𝑧,
{
{
{
𝑀𝑥
| 𝑀∗
𝑥
𝑀𝑦
| 𝑀∗
𝑦
𝑀𝑥𝑦
| 𝑀∗
𝑥𝑦
}
}
}
=
𝑛
∑
𝐿=1
∫
ℎ/2
−ℎ/2
{
{
{
𝜎𝑥
𝜎𝑦
𝜏𝑥𝑦
}
}
}
[𝑧 | 𝑧3] 𝑑𝑧.
(11)
Journal of Composites 5
And the transverse force resultants and the inertias are givenby
{
𝑄𝑥
| 𝑆𝑥
| 𝑄∗
𝑥
𝑄𝑦
| 𝑆𝑦
| 𝑄∗
𝑦
}
=
𝑛
∑
𝐿=1
∫
ℎ/2
−ℎ/2
{
𝜏𝑥𝑧
𝜏𝑦𝑧
} [1 | 𝑧 | 𝑧2] 𝑑𝑧,
(12)
𝐼1, 𝐼2, 𝐼3, 𝐼4, 𝐼5, 𝐼6, 𝐼7
= ∫
ℎ/2
−ℎ/2
[(𝜌𝑐− 𝜌𝑚) (
2𝑧 − ℎ
2ℎ
)
𝑛
+ 𝜌𝑚]
× (1, 𝑧, 𝑧2, 𝑧3, 𝑧4, 𝑧5, 𝑧6) 𝑑𝑧.
(13)
The resultants in (11)-(12) can be related to the total strains in(5) by the following matrix:
:
{{{{{{{{{{
{{{{{{{{{{
{
𝑁
𝑁∗
—𝑀
𝑀∗
—𝑄
𝑄∗
}}}}}}}}}}
}}}}}}}}}}
}
=[
[
𝐴 | 𝐵 | 0
𝐵𝑡| 𝐷𝑏
| 0
0 | 0 | 𝐷𝑠
]
]
{{{{{{{{{{
{{{{{{{{{{
{
𝜀0
𝜀∗
0
—𝐾𝑠
𝐾∗
—𝜑
𝜑∗
}}}}}}}}}}
}}}}}}}}}}
}
, (14)
where
𝑁 = [𝑁𝑥𝑁𝑦
𝑁𝑥𝑦]
𝑡
; 𝑁∗
= [𝑁∗
𝑥𝑁∗
𝑦𝑁∗
𝑥𝑦]
𝑡
. (15)
𝑁, 𝑁∗ are called the inplane force resultants
𝑀 = [𝑀𝑥𝑀𝑦
𝑀𝑥𝑦]
𝑡
; 𝑀∗= [𝑀
∗
𝑥𝑀∗
𝑦𝑀∗
𝑥𝑦]
𝑡
. (16)
𝑀, 𝑀∗ are called as moment resultants
𝑄 = [𝑄𝑥𝑄𝑦]
𝑡
; 𝑄∗= [ 𝑆𝑥
𝑆𝑦
𝑄∗
𝑥𝑄∗
𝑦 ] .
𝑡 (17)
𝑄, 𝑄∗ denote the transverse force result
𝜀0= [𝜀𝑥0
𝜀𝑦0
𝜀𝑥𝑦0]
𝑡
; 𝜀∗
0= [𝜀∗
𝑥0𝜀∗
𝑦0𝜀∗
𝑥𝑦0]
𝑡
𝐾𝑠= [𝐾𝑥
𝐾𝑦
𝐾𝑥𝑦]
𝑡
; 𝐾∗= [𝐾∗
𝑥𝐾∗
𝑦𝐾∗
𝑥𝑦]
𝑡
𝜑 = [𝜑𝑥𝜑𝑦]
𝑡
; 𝜑∗= [𝜀𝑥𝑧0
𝜀𝑦𝑧0
𝜑∗
𝑥𝜑∗
𝑦]
𝑡
.
(18)
The matrices [A], [B], [D], and [𝐷s] are the plate stiffnesswhose elements can be calculated using (5) and (11)-(12).
3. Analytical Solution for the SimplySupported Plate
Let 𝑎 be simply supported rectangular plate with length𝑎 and width 𝑏 which is subjected to in-plane loading intwo directions (𝑁
𝑥= −𝜆
1𝑁cr, 𝑁𝑦 = −𝜆
2𝑁cr, 𝑁𝑥𝑦 =
0). The following expressions of displacements are chosen
based onNavier’s approach to automatically satisfy the simplysupported boundary conditions of the plate:
𝑢0(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑈𝑚𝑛
cos𝛼𝑥 sin𝛽𝑦, (19a)
V0(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑉𝑚𝑛
sin𝛼𝑥 cos𝛽𝑦, (19b)
𝑤0(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑊𝑚𝑛
sin𝛼𝑥 sin𝛽𝑦, (19c)
𝜃𝑥(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑋𝑚𝑛
cos𝛼𝑥 sin𝛽𝑦, (19d)
𝜃𝑦(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑌𝑚𝑛
sin𝛼𝑥 cos𝛽𝑦, (19e)
𝑢∗
0(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑈∗
𝑚𝑛cos𝛼𝑥 sin𝛽𝑦, (19f)
𝑉∗
𝑜(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑉∗
𝑚𝑛sin𝛼𝑥 cos𝛽𝑦, (19g)
𝜃∗
𝑥(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑋∗
𝑚𝑛cos𝛼𝑥 sin𝛽𝑦, (19h)
𝜃∗
𝑦(𝑥, 𝑦, 𝑡) =
∞
∑
𝑚=1
∞
∑
𝑛=1
𝑌∗
𝑚𝑛sin𝛼𝑥 cos𝛽𝑦. (19i)
The eigenproblem related to governing equations isdefined as
([𝑆]9 × 9
− 𝜆[𝜍]9 × 9
)X = 0, (20)
where [𝑆] collects all stiffness terms and [𝜍] collects all termsrelated to the in-plane forces. In (20) X are the modes ofbuckling associated with the buckling loads defined as 𝜆. Foreach value of 𝑚 and 𝑛, there is a unique value of 𝑁cr. Thecritical buckling load is the smallest value of𝑁cr (𝑚, 𝑛).
4. Results and Discussion
4.1. Comparative Studies. To validate the accuracy of thepresent higher order theory in predicting the critical bucklingload of FG plates subjected to different in-plane loadingconditions (uniaxial compression: 𝜆
1= −1, 𝜆
2= 0; biaxial
compression: 𝜆1
= −1, 𝜆2
= −1; biaxial compression andtension 𝜆
1= −1, 𝜆
2= 1), four numerical examples are
presented and discussed. The material properties adoptedhere are as follows.
Aluminium Young’s modulus (E𝑚): 70GPa, density 𝜌
𝑚=
2702 kg/m3, and Poisson’s ratio (𝜐): 0.3.Alumina Young’s modulus (E
𝑐): 380GPa, density 𝜌
𝑐=
3800 kg/m3, and Poisson’s ratio (𝜐): 0.3.
6 Journal of Composites
Table 1: Comparison of nondimensionalized critical buckling load (𝑁) of simply supportedAl/Al2O3 plate subjected to uniaxial compressionalong the 𝑥-axis (𝜆
[27] are presented in Table 1. It can be observed that thepredicted results are almost identical with that of Thai andChoi [27]. It can also be seen that the critical buckling loaddecreases with the increase of power-law index value whileit increases with the increase of side-to-thickness ratio andaspect ratio. Furthermore, increasing of thickness ratio andaspect ratio not only increases the critical buckling load
values, but also causes the changes in critical bucklingmodes.This can be observed when the plate is subjected to inplanecompression along the 𝑥-axis with aspect ratio value of 2.The critical buckling mode varies from (3,1) to (2,1) as theside-to-thickness ratio increases from 5 to 10. The maximumpercentage error between present HSDT and Thai and Choi[27] is 3.9 at 𝑎/ℎ = 5, 𝑎/𝑏 = 2, and 𝑛 = 2.
8 Journal of Composites
Table 3: Comparison of nondimensionalized critical buckling load (𝑁) of simply supported Al/Al2O3 plate subjected to biaxial compressionand tension (𝜆
aMode for plate is (𝑚, 𝑛) = (2, 1).bMode for plate is (𝑚, 𝑛) = (1, 2).
Example 2. The next comparison is performed for the simplysupported FG plate subjected to inplane biaxial compression(𝜆1= −1, 𝜆
2= −1). The results of critical buckling loads
obtained by the present theory and those reported by Thaiand Choi [27] are presented in Table 2 and observed the
close agreement between the results. It can be seen that, inthis loading condition also, the nondimensionalized criticalbuckling load decreases with the increase of power-law index,while it increases with the increase of aspect ratio and side-to-thickness ratio, but only one critical buckling mode exists
Journal of Composites 9
Table 4: Comparison of nondimensionalized critical buckling load (𝑁) of simply supportedAl/Al2O3 plate subjected to uniaxial compression(𝜆1= −1, 𝜆
in any case of the aspect ratio, thickness ratio, modulusratio, and power-law index. The maximum percentage errorbetween present HSDT and Thai and Choi [27] is 2.90 at𝑎/ℎ = 5, 𝑎/𝑏 = 2, and 𝑛 = 2.
Example 3. The third comparison is carried out for the simplysupported FG plates under inplane biaxial compression andtension (𝜆
1= −1, 𝜆
2= 1). The results predicted by present
theory are compared with Thai and Chois [27] results andgood agreement between the results can be observed. Theresults are presented in Table 3. It can also be seen that, underbiaxial compression and tension, the critical buckling loaddecreases with the increase of power-law index value whileit increases with the increase of side-to-thickness ratio andaspect ratio, same as in uniaxial and biaxial compression.Also, increasing of thickness ratio and aspect ratio not onlyincreases the critical buckling load values, but also causeschanges in critical buckling modes. This can be observedwhen the aspect ratio value is 1.5. The critical buckling modevaries from (2,1) to (1,2) for the values and the side-to-thickness ratio increases from 5 to 10 and 20 to 50. Themaximum percentage error between present HSDT andThaiand Choi [27] is 2.90 at 𝑎/ℎ = 5, 𝑎/𝑏 = 2, and 𝑛 = 2.
Example 4. The last comparison is carried out for the simplysupported FG plates under inplane uniaxial compression(𝜆1= −1, 𝜆
2= 0).The results predicted by present theory are
compared with the classical plate theory (CPT) [16] and first-order shear deformation theory (FSDT) [16] are presentedin Table 4. The maximum percentage error between presentHSDT and CPT and HSDT and FSDT for 𝑛 = 0.1, 1, and10, respectively are 5.633954, 0.19667, 5.195466, and 0.42725,8.294579, 4.39547 at 𝑎/ℎ = 10. It is also seen that, the CPT andFSDT overpredicts the critical buckling load for all thicknessratios.
4.2. Parametric Study. The effect of side-to-thickness ratio,aspect ratio, and the modulus ratio on nondimensionalizedcritical buckling load for simply supported FG plate made ofAl/Al2O3with 𝜀
𝑧= 0 is investigated. Figures 2–4 represent
the variation of nondimensionalized critical buckling loadwith side-to-thickness ratio, aspect ratio, and modulus ratio,respectively, under uniaxial compression. From Figures 2–4,it is important to observe that the nondimensionalized criti-cal buckling loads are higher for ceramic rich plates and lowerformetal rich plates.The critical buckling loads of FGMplates
0.02.04.06.08.0
10.012.014.016.018.020.0
0 5 10 15 20 25 30 35 40
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Side-to-thickness ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 2: Effect of side-to-thickness ratios (a/h) on nondimension-alized critical buckling load (𝑁) under uniaxial compression for asimply supported FG plate for variousmaterial variation parameters(n).
0.0
50.0
100.0
150.0
200.0
250.0
300.0
350.0
400.0
450.0
0 0.5 1 1.5 2 2.5 3
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Aspect ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 3: Effect of aspect ratios (a/b) on nondimensionalizedcritical buckling load (𝑁) under uniaxial compression for a simplysupported FG plate for various material variation parameters (n).
10 Journal of Composites
0.0
5.0
10.0
15.0
20.0
25.0
30.0
35.0
0 0.1 0.2 0.3 0.4 0.5
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Modulus ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 4: Effect of modulus ratio (E𝑚/E𝑐) on nondimensionalized
critical buckling load (𝑁) under uniaxial compression for a simplysupported FG plate for various material variation parameters (n).
0.0
1.0
2.0
3.0
4.0
5.0
6.0
7.0
8.0
9.0
10.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Side-to-thickness ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 5: Effect of side to thickness ratios (a/h) on nondimension-alized critical buckling load (𝑁) under biaxial compression for asimply supported FG plate for variousmaterial variation parameters(n).
are intermediate to that of ceramic and metal. From Figure 2it is seen that the effect of shear deformation is significantfor a side-to-thickness ratio less than 10 and diminishes withincrease of side-to-thickness ratio.The increase of aspect ratioincreases the critical buckling load due to increase of stiffnessof the plate as shown in Figure 3.
From Figure 4, it can be seen that the nondimensionalcritical buckling load decreases as the metal-to-ceramicmodulus ratio increases and decreases as the power law indexincreases.This is due to the fact that higher values ofmetal-to-ceramic modulus ratio correspond to high portion of metal.
0.0
5.0
10.0
15.0
20.0
25.0
30.0
35.0
40.0
0 0.5 1 1.5 2 2.5 3
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Aspect ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 6: Effect of aspect ratios (a/b) on nondimensionalizedcritical buckling load (𝑁) under biaxial compression for a simplysupported FG plate for various material variation parameters (n).
0.0
2.0
4.0
6.0
8.0
10.0
12.0
14.0
16.0
18.0
0 0.1 0.2 0.3 0.4 0.5
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Modulus ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 7: Effect of modulus ratios (E𝑚/E𝑐) on nondimensionalized
critical buckling load (𝑁) under biaxial compression for a simplysupported FG plate for various material variation parameters (n).
The effect of side-to-thickness ratio, aspect ratio,modulusratio, and power law index values on nondiensionalizedcritical buckling load for a simply supported FG plate underbiaxial compression is shown in Figures 5, 6, and 7.The samecan be observed as in the case of uniaxial compression. It isimportant to observe that the critical buckling loads are largerin uniaxial compression and smaller in biaxial compression.Figure 8 shows the variation of nondimensionalized criticalbuckling load of different modulus ratios and power lawindex values under inplane compression and tension. It canbe observed that critical buckling load decreases with theincrease of modulus ratio and power-law index values.
Journal of Composites 11
0
10
20
30
40
50
60
70
0 0.1 0.2 0.3 0.4 0.5
Non
dim
ensio
naliz
ed cr
itica
l buc
klin
g lo
ad
Modulus ratio
Ceramic n = 0.5n = 1 n = 4n = 10 Metal
Figure 8: Effect of modulus ratios (E𝑚/E𝑐) on nondimensionalized
critical buckling load (𝑁) under biaxial compression and tensionfor a simply supported FG plate for various material variationparameters (n) (mode for the plate (𝑚, 𝑛) = (1,2)).
5. Conclusions
A higher order shear deformation theory was successfullydeveloped and applied to study the buckling behavior offunctionally graded plates without enforcing zero transverseshear stresses on the top and bottom surfaces of the plate.The present formulation was compared with the refinedtheory developed by Thai and Choi [27] and proved veryaccurate for the buckling problem. This eliminated the needof using shear correction factors. It can be concluded thatthe present theory is accurate and efficient in predicting thecritical buckling load of simply supported FG plates. Hence,the present findings will be useful benchmark for evaluatingother future plate theories and numericalmethods such as thefinite element and meshless methods.
References
[1] R. Javaheri and M. R. Eslami, “Buckling of functionally gradedplates under in-plane compressive loading,” Journal of AppliedMathematics and Mechanics, vol. 82, no. 4, pp. 277–283, 2002.
[2] S. Abrate, “Functionally graded plates behave like homogeneousplates,” Composites B, vol. 39, no. 1, pp. 151–158, 2008.
[3] M.Mohammadi, A. R. Saidi, and E. Jomehzadeh, “Levy solutionfor buckling analysis of functionally graded rectangular plates,”Applied Composite Materials, vol. 17, no. 2, pp. 81–93, 2010.
[4] M. Mahdavian, “Buckling analysis of simply-supported func-tionally graded rectangular plates under non-uniform in-planecompressive loading,” Journal of Solid Mechanics, vol. 1, no. 3,pp. 213–225, 2009.
[5] E. Feldman and J. Aboudi, “Buckling analysis of functionallygraded plates subjected to uniaxial loading,” Composite Struc-tures, vol. 38, no. 1–4, pp. 29–36, 1997.
[6] B. A. S. Shariat, R. Javaheri, and M. R. Eslami, “Buckling ofimperfect functionally graded plates under in-plane compres-sive loading,” Thin-Walled Structures, vol. 43, no. 7, pp. 1020–1036, 2005.
[7] H. V. Tung and N. D. Duc, “Nonlinear analysis of stabilityfor functionally graded plates under mechanical and thermalloads,” Composite Structures, vol. 92, no. 5, pp. 1184–1191, 2010.
[8] E. Reissner, “The effect of transverse shear deformation on thebending of elastic plates,” ASME Journal of Applied Mechanics,vol. 12, no. 2, pp. 69–77, 1945.
[9] R. D. Mindlin, “Influence of rotary inertia and shear on flexuralmotions of isotropic, elastic plates,” ASME Journal of AppliedMechanics, vol. 18, pp. 31–38, 1951.
[10] H. T. Thai and T. P. Vo, “A new sinusoidal shear deformationtheory for bending, buckling, and vibration of functionallygraded plates,” Applied Mathematical Modelling, vol. 37, no. 5,pp. 3269–3281, 2013.
[11] J. Yang, K.M. Liew, and S. Kitipornchai, “Second-order statisticsof the elastic buckling of functionally graded rectangular plates,”Composites Science and Technology, vol. 65, no. 7-8, pp. 1165–1175, 2005.
[12] M. Mohammadi, A. R. Saidi, and E. Jomehzadeh, “A novelanalytical approach for the buckling analysis of moderatelythick functionally graded rectangular plates with two simply-supported opposite edges,” Proceedings of the Institution ofMechanical Engineers C, vol. 224, no. 9, pp. 1831–1841, 2010.
[13] X. Zhao, Y. Y. Lee, and K. M. Liew, “Mechanical and thermalbuckling analysis of functionally graded plates,” CompositeStructures, vol. 90, no. 2, pp. 161–171, 2009.
[14] H. A. Sepiani, A. Rastgoo, F. Ebrahimi, and A. G. Arani,“Vibration and buckling analysis of two-layered functionallygraded cylindrical shell, considering the effects of transverseshear and rotary inertia,” Materials and Design, vol. 31, no. 3,pp. 1063–1069, 2010.
[15] A. Naderi and A. R. Saidi, “On pre-buckling configurationof functionally graded Mindlin rectangular plates,” MechanicsResearch Communications, vol. 37, no. 6, pp. 535–538, 2010.
[16] R. Saha and P. R. Maiti, “Buckling of simply supported FGMplates under uniaxial load,” International Journal of Civil andStructural Engineering, vol. 2, no. 4, pp. 1036–1050, 2012.
[17] R. Javaheri andM. R. Eslami, “Thermal buckling of functionallygraded plates based on higher order theory,” Journal of ThermalStresses, vol. 25, no. 7, pp. 603–625, 2002.
[18] M. M. Najafizadeh and H. R. Heydari, “Thermal buckling offunctionally graded circular plates based on higher order sheardeformation plate theory,”European Journal ofMechanics A, vol.23, no. 6, pp. 1085–1100, 2004.
[19] M. Bodaghi and A. R. Saidi, “Levy-type solution for bucklinganalysis of thick functionally graded rectangular plates basedon the higher-order shear deformation plate theory,” AppliedMathematical Modelling, vol. 34, no. 11, pp. 3659–3673, 2010.
[20] E. Bagherizadeh, Y. Kiani, and M. R. Eslami, “Mechanicalbuckling of functionally graded material cylindrical shells sur-rounded by Pasternak elastic foundation,”Composite Structures,vol. 93, no. 11, pp. 3063–3071, 2011.
[21] H. Mozafari and A. Ayob, “Effect of thickness variation on themechanical buckling load in plates made of functionally gradedmaterials,” Procedia Technology, vol. 1, pp. 496–504, 2012.
[22] L. S. Ma and T. J. Wang, “Nonlinear bending and post-bucklingof a functionally graded circular plate under mechanical andthermal loadings,” International Journal of Solids and Structures,vol. 40, no. 13-14, pp. 3311–3330, 2003.
12 Journal of Composites
[23] S. Hosseini-Hashemi, K. Khorshidi, and M. Amabili, “Exactsolution for linear buckling of rectangular Mindlin plates,”Journal of Sound and Vibration, vol. 315, no. 1-2, pp. 318–342,2008.
[24] A. R. Saidi, A. Rasouli, and S. Sahraee, “Axisymmetric bendingand buckling analysis of thick functionally graded circularplates using unconstrained third-order shear deformation platetheory,” Composite Structures, vol. 89, no. 1, pp. 110–119, 2009.
[25] O. O. Oyekoya, D. U.Mba, and A.M. El-Zafrany, “Buckling andvibration analysis of functionally graded composite structuresusing the finite element method,” Composite Structures, vol. 89,no. 1, pp. 134–142, 2009.
[26] S. A. M. Ghannadpour, H. R. Ovesy, and M. Nassirnia, “Buck-ling analysis of functionally graded plates under thermal load-ings using the finite strip method,” Computers and Structures,vol. 108-109, pp. 93–99, 2012.
[27] H. T. Thai and D. H. Choi, “An efficient and simple refinedtheory for buckling analysis of functionally graded plates,”Applied Mathematical Modelling, vol. 36, no. 3, pp. 1008–1022,2012.
[28] B. Uymaz and M. Aydogdu, “Three dimensional mechanicalbuckling of FG plates with general boundary conditions,”Composite Structures, vol. 96, pp. 174–193, 2013.
[29] A. Lal, K. R. Jagtap, and B. N. Singh, “Post buckling responseof functionally graded materials plate subjected to mechanicaland thermal loadingswith randommaterial properties,”AppliedMathematical Modelling, vol. 37, no. 5, pp. 2900–2920, 2013.