Effects of elastic foundation and the Poisson’s ratio v = v z ( ) on the nonlinear buckling and postbuckling behaviors of imperfect FGM plates subjected to mechanical loads.. Pham [r]
Trang 1226
Effects of elastic foundation and the Poisson’s ratio v = v z ( )
on the nonlinear buckling and postbuckling behaviors
of imperfect FGM plates subjected to mechanical loads
Pham Hong Cong*
VNU University of Engineering and Technology
Received 11 May 2012, received in revised form 20 May 2012
Abstract: This paper presents an analytical approach to investigate effects of elastic foundation
and the Poisson’s ratio v=v z( ) on the nonlinear buckling behavior of imperfect FGM plates, subjected to mechanical loads Material properties are assumed to be temperature independent, and graded in the thickness direction according to a power law distribution in terms of volume fractions of constituents Equilibrium and compatibility equations are derived by using classical plate theory taking into account geometrical nonlinearity, initial geometrical imperfection and elastic foundation with Pasternak model Galerkin method is used to determine explicit expressions of buckling loads and postbuckling paths Analysis is carried out to assess the effects
of material, geometrical, elastic foundation parameters on the stability of FGM plates
Keywords: Buckling and postbuckling, Functionally graded material, Plate, Elastic foundations, Poisson’s ratio v=v z( )
1 Introduction∗∗∗∗
Due to high performance is heat resistance capacity and excellent characteristics in comparison with conventional composites, Functionally Graded Materials (FGMs) which are microscopically composites and composed from mixture of metal and ceramic constituents have attracted considerable attention recent years By continuously and gradually varying the volume fraction of constituent materials through a specific direction, FGMs are capable of withstanding ultrahigh temperature environments and extremely large thermal gradients Therefore, these novel materials are chosen to use in structure components of aircraft, aerospace vehicles, nuclear plants as well as various temperature shielding structures widely used in industries Buckling and postbuckling behaviors of FGM structures under different types of loading are important for practical applications and have received considerable interest Eslami and his co-workers used analytical approach, classical and higher order plate theories in conjunction with adjacent equilibrium criterion to investigate the _
∗
Tel.:+84 1649 589 562
Email: congph_54@vnu.edu.vn
Trang 2buckling of FGM plates with and without imperfection under mechanical and thermal loads [3,8,9] According to this direction, Lanhe [10] also employed first order shear deformation theory to obtain closed-form relations of critical buckling temperatures for simply supported FGM plates Zhao et al [14] analyzed the mechanical and thermal buckling of FGM plates using element-free Ritz method Liew et al [15,16] used the higher order shear deformation theory in conjunction with differential quadrature method to investigate the postbuckling of pure and hybrid FGM plates with and without imperfection on the point of view that buckling only occurs for fully clamped FGM plates The postbuckling behavior of pure and hybrid FGM plates under the combination of various loads were also treated by Shen [17,18] using two-step perturbation technique taking temperature dependence of material properties into consideration Recently, Lee et al [19] made of use element-free Ritz method
to analyze the postbuckling of FGM plates subjected to compressive and thermal loads
The components of structures widely used in aircraft, reusable space transportation vehicles and civil engineering are usually supported by an elastic foundation Therefore, it is necessary to account for effects of elastic foundation for a better understanding of the postbuckling behavior of plates and shells Librescu and Lin have extended previous works [20] to consider the postbuckling behavior of flat and curved laminated composite panels resting on Winkler elastic foundations [20] In spite of practical importance and increasing use of FGM structures, investigation on FGM plates and shells supported by elastic media are limited in number The bending behavior of FGM plates resting on Pasternak type foundation has been studied by Huang et al [21] using state space method, Zenkour [22] using analytical method and by Shen and Wang [23] making use of asymptotic perturbation technique To the best of authors’ knowledge, there is no analytical studies have been reported in the literature on the postbuckling of thick FGM plates resting on elastic foundations.In [11], the authors Dao Van Dung and Nguyen Thi Nga have studied the stability of the composite FGM plate when ( )
v=v z and E=E z( )(without elastic foundation) In [12], the author Do Nam has studied the stability of the FGM plate on the elastic foundation with classical plate theory, in [6] the authors Nguyen Dinh Duc and Hoang Van Tung have studied postbuckling of the high order shear deformable FGM plates on elastic foundation, but these studies assume E=E z( ) and v=const
The aim of the paper proposed is of studying the nonlinear stability of FGM plate on the elastic foundation under the effect of the load in the case both elastic modules are variable v=v z( ) and ( )
E=E z , the study of the effect of the initial imperfect shape, proportion metal-ceramic, and the elastic foundation parameters and the geometric parameters on the nonlinear stability of the FGM plate
2 FGM plates on elastic foundations
Consider a rectangular functionally graded plate of length a, width band thickness h An orthogonal coordinate system Oxyz is choose so that the plane coincides Oxy with the middle surface of the plate and the axis Oz is in the thickness direction (−h/ 2≤ ≤z h/ 2)as shown in Fig 1
Trang 3y
z
h a
b
shear layer
Fig.1.Geometry and coordinatr system of an FGM plate on elastic foundation
The volume fractions of constituents are assumed to vary through the thickness according to the following power law distribution
2
N
h
+
where N is volume fraction index (0 ≤ N < ∞) Effective properties Preff of FGM panel are determined by linear rule of mixture as
Pr ( )eff z = PrmV zm( ) + PrcV zc( ) (2)
where Pr denotes a temperature independent material property, and subscripts m and c stand for the metal and ceramic constituents, respectively
Specific expressions of modulus of elasticity E, Poisson ratio ν the coefficient of thermal expansion α are obtained by substituting Eq (1) into Eq (2) as [11]
2
2
N
N
z h
h
+
1
1 2
2
N
N
z h
h
Where
2
2
+
It is evident from Eqs (3), (4) that the upper surface of the plate (z = − h / 2) is ceramic-rich, while the lower surface (z = h / 2) is metal-rich, and the percentage of ceramic constituent in the panel is enhanced when N increases
The plate – foundation interaction is represented by Pasternak model as [7]
qe = k w k1 − ∇2 2w (5)
Trang 4where ∇ = ∂2 2/ ∂ x2+ ∂2/ ∂ y2, w is the deflection of the plate, k1 is Winkler foundation modulus and k2 is the shear layer foundation stiffness of Pasternak model
2.1 Theoretical formulation
For imperfect plates, the strain components on the reference surface with the geometrical nonlinearity in von Karman sense, are [1]
0 , , , , , , , ,
u
v
ε
ε
γ
(6)
Where u=u x y v( , ), =v x y( , ) and w=w x y( , )are displacements along x y, and z
respectively; w*=w x y*( , ) denotes an initial imperfections of plate The quantity w* is assumed small
The strains across the plate thickness at a distance z from the mid-plane are [1]
Hooke law for an FGM plate is defined as
2
1
E
v E
v
=
+
(8)
where we assume that the plate is subjected to a uniform temperature rise i.e ∆T is a constant So, the force and moment resultants are expressed a
/ 2 , ,
/ 2 / 2
/ 2
h
h h
h
σ σ σ
σ σ σ
−
−
=
=
∫
∫
(9)
Substituting relations (3), (6) ÷ (8) into (9), after series of calculations, we have
10 20 0 20 10 0 11 21 21 11 1
30 0 31
11 21 0 21 11 0 12 22 22 12 2
2
2
γ
γ
(10)
Trang 5ij( 1, 2, 3; 0,1, 2), 1, 2
defined as follows
/ 2
/ 2
/ 2
/ 2
/ 2
/ 2
/ 2
1
/ 2
/ 2
2
/ 2
( )
;
( ) ( )
;
( ) ( )
;
( ) ( )
h
j j
h
h
j j
h
h
j
h
h
h
h
h
E z
v z
E z v z
v z
E z
v z
dz
v z
zdz
v z
α φ
α φ
−
−
−
−
−
=
−
=
−
+
= −
−
= −
−
∫
∫
∫
∫
∫
(11)
The equilibrium equations of a imperfect plate on elastic foundations are in the form [2÷4]:
, ,
*
0 0
We introduce Airy’s stress function f = f x y( , ) so that
N = f N = f N = −f
(13)
It is easy to see that the first two equations in (12) are automatically satisfied
Substituting relations (13) in Eqs (10), we obtain
0 1 10 , 20 , 2 , 3 , 4 1
0 1 10 , 20 , 2 , 3 , 4 1
γ
(13) where
2 2
1 10 20 2 10 11 20 21
3 10 21 20 11 4 10 20
,
= − = − (15)
Substituting once again the expressions of Eq (14) into the relations of internal moments Mijin (10) we obtain
1 , ( 11 20 21 10 ) 1 , ( 11 10 21 20 ) 1 w , ( 11 2 21 3 ) 1 w , ( 11 3 21 2 ) 12 w , 22 w
M =I f −J J +J J +I f J J −J J +I J I +J I +I J I +J I −J −J
Trang 61 , ( 21 20 11 10 ) 1 , ( 21 10 11 20 ) 1 w , ( 21 2 11 3 ) 1 w , ( 21 3 11 2 ) 22 w , 12 w
M =I f −J J +J J +I f J J −J J +I J I +J I +I J I +J I −J −J
31
31 , , 32 , 30
J
J
(16) The substituting (16) into the third equation of (12) we have
1 2 w ,yy(w,xx w ) 2,xx ,xy(w,xy w ),xy ,xx(w,yy w,yy) 1w+ 2 w 0
(17) where
4
4 2 2 2 4
1 1 3; 2 1( 11 2 21 3) 12
A =I I A =I J I +J I −J
The equation (17) includes two unknowns functions w and f , so it is necessary to find a second equation relating to these two unknowns functions by using the compatible equation:
2
0 ,x yy 0 ,y xx 0xy xy, (w,xy) w w,xx ,yy
(18) The substituting the above expressions of εij in (14) into Eqs (18), we obtain
3 w 4 w,xy w w,xx ,yy 2 ,xyw,xy w w,xx yy w,yyw,xx 0
(19)
in wich
3
10 1 10
1
;
I
The couple of Eqs (17) and (19) are the governing equations used to investigate the nonlinear stability of imperfect FGM plates with the Poisson’s Ratio ν ν= ( )z resting on elastic foundation
In the case w* =0, from (17) and (19) we obtain the basic stability equations for perfect FGM plates
2.2 Boundary conditions and the solution of the problem
Suppose that three cases boundary conditions for a rectangular plate will be considered follow as [5, 6]:
Case 1 The edges of plate are simply supported and freely movable (FM) The associated
boundary conditions are:
x=0,x=a; w=M x =N xy =0;N x =N x0
(20)
Case 2 The edges of plate are simply supported and immovable (IM) The associated boundary
conditions are:
Trang 7
0
x=0,x=a; w =M x =u=0;N x =N x
(21)
y=0,y=b; w=M y = =v 0;N y =N y0
Case 3 The edges of plate are simply supported Uniaxial edge loads are applied in the direction of
the x -coordinate The edges x=0, x = a are considered freely movable, the remaining two edges being unloaded and immovable The boundary conditions, for this case, are
x=0,x=a; w=M x=N xy =0;N x =N x0
(22)
y=0,y=b; w=M y = =v 0;N y =N y0
where u, v are the displacement components x, y directions, respectively, N xy,M x,M y are force and moment resultants Moreover,N , Nx0 y0 are prebuckling force resultants in the direction x and
y respectively, for Case 1 and the first of Case 3 or are fictitious compressive edge loads at immovable edges (Case 2 and the second of Case 3)
Approximate solutions of basic equations (17) and (19) are assumed as [5, 6]:
*
=
=
(23)
which fulfill exactly the out-of-plane boundary conditions and satisfy in-plane boundary conditions in
an average sense Moreover, λm =mπ/ ;a δn =nπ/ ; ,b m n=1, 2, Wis amplitude of deflection and
µ is imperfection parameter (0≤µ≤1)
By substituting expression (23) into Eq (19), we obtain
2 4
2 4
4
32
32 W 0
n m
m n
C
C
C
λ
δ
+
=
+
=
= −
= (24)
Introducing Eqs (24) and (23) into the left side of Eq (17) and the applying Galerkin method we obtain equation for determining buckling loads and postbuckling curves of rectangular FGM plates subjected to mechanical, thermal and combined loads
Trang 8
2
32
3
16
m n
A
ab
mn
λ δ
π
+
where m n, are odd numbers This equation will be used to analyze the buckling and postbuckling
behaviors of thick FGM plates under mechanical, thermal and thermomechanical loads on elastic
foundation
3 Nonlinear stability analysis FGM plate on elastic foundation
3.1 Mechanical stability analysis
Consider a rectangular imperfect FGM plate being simply supported at its edges and freely
movable (Case 1) and subjected to the in-plane compressive only loads P x uniformly distributed
along the edges x=0,a In this case, the prebuckling force resultants are given [7]
Substituting this expression (26) into Eq (25) we receive :
2 2 2 2 2 2 2 2 2
2 4 4 4
4
2 2 2
1
8
x
a
a h
P
B B m
+
(27) where :
3
3
h
a
A
k a k a D
(28)
for perfect FGM plate we have
0
2 2 2 2 2 2 2 2 2
2
2 4 4 4
4
2 2 2
1
8
x
a
a h
P
B B m
π
π π
+ +
(29)
Trang 9If v = c ons t we have A1 =0, A3 =0
2 2 2 2 2 ' 2 ' 2 2 2 2 4 4 4
8 W
x
P
+
(30)
If ν = const ; µ = 0, for perfect FGM plate we receive :
2
8
x
P
π
Where
1
1 2; 1 D3
h
;
Looking at the expression (27) and (29÷31) we see the beneficial effects of elastic foundation to rainbow load capacity and after rainbow
From the equation (29) we can see that for the perfect plate (µ =0) the function Px(W) will reach a minimum at W = 0 and P x(0)is the lowest point of the deflection-load graph
3.2 Numerical results and discussions
The purpose of this section is to explore the dependence of the critical force on the coefficientsK1, K2 of the elastic foundation in some cases when the plate is perfect and imperfect
To illustrate the present approach for buckling and postbuckling analysis of thick FGM plates resting on elastic foundations, consider a square ceramic–metal plate consisting of aluminum and alumina with the following properties [3,8,9,10]
In figures, W/h denotes the dimension-less maximum deflection and the FGM plate foundation interaction is ignored, unless otherwise stated
Trang 10Fig.2 The influence of imperfections on the stability of FGM plates under compression
Fig.3 The effects of the area ratio coefficient N on the stability of FGM plates under compression