-
Analytical solution
nonrectangular plate with in-plane
Variable stiffness
Tian-
chong YU
,
Guo-jun
NIE
,
Zheng
ZHONG
,
Fu-yun CHU
(School of Aerospace Engineering and
Applied Mechanics, Tongji
University,Shanghai 200092, P.R. China)
Abstract:
The
bending
problem
of
a
thin
rectangular
plate
with
in-
plane
variable
stiffness is
studied. The basic equation is formulated for the
two-opposite-edge simply
supported
rectangular plate under the distributed loads. The
formulation is based on
the assumption
that the flexural
rigidity of the plate
varies in
the plane
following
a
power
form
,
and
Poisson
’
s
ratio
is
constant
.
A
fourth-order
partial
differential
equation with
variable coefficients is derived by assuming a
Levy-type form for the
transverse
displacement. The governing equation can be
transformed into a Whittaker
equation
,
and an
analytical solution is obtained for a thin
rectangular plate subjected
to the
distributed loads
.
The
validity of the present solution is shown by
comparing
the
present
results
with
those
of
the
classical
solution
.
The
influence
of
in-
plane
variable
stiffness
on
the
deflection
and
bending
moment
is
studied
by
numerical
examples
.
The
analytical solution presented here is useful in
the design of rectangular
plates with
in-plane variable stiffness.
Keywords:
in-
plane variable stiffness ,power
form
,
Levy-type
solution
,
rectangular
plate
Chinese
Library Classification 0343
2010
Mathematics Subject Classification 74B05
第
1
页
1 Introduction
The
term” variable
stiffness” implies
that the
stiffness parameters vary spatially
throughout.
The
s
tructure
[1]
.
Funct
ionally
graded
materials
(FGMs)
are
inhomogeneous
composites
,
in which. the
mechanical properties vary smoothly with
the
position
to
meet
the
predetermined
functional.
performance
.
The
structures
composed of the
FGMs are of variable
stiffness
.
There
are
extensive
literatures
on
the
bending
,
vibration
,
and
fracture
of
the
FGM structures
[2-9]
.
The deformation of a functionally
graded beam was studied by
the direct
approach
[10]
.
An efficient and simply refined theory
was presented for the
buckling analysis
of functionally graded plates by Thai and
Choi
[11]
.
Jodaei
et a1
[12]
dealt
with the three
—
dimensional
analysis of functionally graded annular plates
using
the
state
—
space
based
differential
quadrature
method(SSDQM)
.
Wen
and
A1iabadi
[13]<
/p>
investigated functionally graded plates
under static and dynamic loads by
the
local integral equation
method(LIEM)
.
There are also
some works on the FGM
shells and
cylinders[14-17]
.
However
,
most of
the studies on the FGMs deal with material
stiffness varying
along
the
thickness
direction
.
The
studies
on
the
plates
with
in
—
plane
variable
stiffness
are
quite
few
.
Shang
[18]
studied
the
rectangular
plates
with
bidirectional
linear
stiffness
with
two
opposite
edges
simply
supported
and
the
other
two
edges
arbitrarily supported under the
distributed
loads
.
Yang
[19]
investigated the structural
analysis of the plates with
unidirectionally varying flexural
rigidity by the Galerkin
line
method
.
Liu
et
a1.
[20]
analyzed
the
free
vibration
of
a
Functionally
graded
isotropic
rectangular
plate
with
in
—
plane
material
in
homogeneity
using
the
Levy-type
so1ution
.
Uymaz et
a1
.
[21]
Considered the functionally graded plates with
properties
Varying
in
an
in
—
plane
direction
based
on
a
five
< br>—
degree
—
of-
freedom
shear deformable plate theory
with different boundary
conditions
.
In
this paper
,
the Levy-type
solution[22-23]is presented for the bending of a
thin
rectangular plate with
in
—
plane variable stiffness
under the distributed load
.
2 Basic equations
Consider a
thin rectangular plate of length A and with B with
in-plane Variable
stiffness
,
as
shown in
Fig.1. Introduce a Cartesian coordinate system
0
—
XYZ such
that
0
≤
X
≤
p>
A
,
0
≤
Y
≤
B
We
assume. That the flexural rigidity of the plate
D=D(X,Y)is a function of X
and Y. The
governing equation of the plate with
in
—
plane Variable stiffness
can be
obtained as
Where W is
the transverse
displacement
,
V is
Poisso
n
′
s
ratio
,
and Q is
the
normal pressure on The
plate
.
第
2
页
It
is
assumed
that
the
flexural
rigidity
of
the
plate
varies
only
along
the
Y-direction according to
the following power form:
Where Y and P are two material
parameters describing the in homogeneity of D,
Do is the flexural rigidity
at
,
Y=0, and
D
b
is the flexural rigidity
at Y=b. In
this ‘case
,
Eq. (1) can be reduced to
3
Solution
The
rectangular
plate
is
assumed
to
be
simply
supported
along
two
opposite
edges parallel to the Y-direction. To
solve the governing equation with the prescribed
boundary conditions, a generalized
Levy-type approach is employed as
where
Y
m
(y) is an unknown function
to be determined.
Substituting Eq. (4)
into Eq. (3) yields the following differential
equation:
第
3
页
The solution of the above
equation consists of two parts, i.e,the general
solution
Y
m
o
of
the homogeneous differential equation and the
particular solution
Y
m
o
of the
no homogeneous differential equation
.
Thus
,
the
solution of Eq.(5) can be expressed
as
.
3.1 General
solution of governing equation (5)
To
find
the
solution
of
Eq.
(5),
we
first
consider
the
following
fourth-order
homogeneous
differential equation with variable coefficients:
The
above
equation
can
be
transformed
into
a
Whittaker
equation,
and
the
solution can be expressed as
with
for p
≠
1 and
p
≠
1/v
for
p=1
,
and
第
4
页
for p=1/v
.
W
k,g
(t) is the
Whittaker function
,
is the exponential integral
function,
I(t) and Kg(t) are the gth-
order modified Bessel
functions
of
the
first
and
second
kinds,
respectively,
and
c1,
c2,
c3,
and
c4
are
the
constants to be
determined.
3.2 Particular solution of
governing equation (5)
From the
solution of the homogeneous differential equation
(7), the solution of
the nonhomogeneous
differential equation (5) can be expressed as
Where
c
l
(t),
c
2
(t),
c
3
(t), and
c
4
(t) satisfy the following
equations:
from which we obtain
第
5
页
where
Then. we have
where
C
1
,
C
2
,C
3
,
and C
4
are the
constants to be determined.
By
substituting Eq. (15) into Eq. (12), the
particular solution of Eq. (5) can be
written as
where
From Eqs.(8) and
(16)
,
the solution of Eq.(5)
can be expressed as
The
coefficients dl, d2, d3, and d4 are determined
from the boundary conditions.
In view
of Eqs.(4) and (17)
,
the
transverse displacements of the rectangular plate
can be expressed as
3.3 Boundary conditions
It
is assumed that the two opposite edges parallel to
the y-direction are simply
supported,and
the
other
two
edges
have
arbitrary
boundary
conditions
such
as
free,
第
6
页
simply supported, or clamped
conditions.
The boundary conditions for
the remaining two edges(Y=0 and Y=b) are given
as follows:
The simply
supported boundary conditions (S) are
The clamped boundary conditions (C) are
The free boundary condition
8(F) are
We consider the
simply supported conditions in Eq. (19) as an
example
.
Thus,
we
have
From
Eq.
(22),
dl,
d2,
d3,
and
d4
can
be
determined,
and
the
results
are
given
in
Appendix
A.
4 Results and discussion
In this section, we make numerical
studies on the static behavior of a rectangular
plate
with
in-
plane
variable
stiffness
under
a
uniform
pressure
(a=b=1
m,v=0.3,
Do=100 Nm,Db=500
kNm, and q=100 kN/m2). The variation of the
bending stiffness
is
described
by
Eq.
(2).
For
different
parameters
p
in
Eq.
(2),
we
calculate
the
deflection w and the
bending Moments M
X
and
M
Y
along the y-direction (at
x=a/2) of
a
four-edge
simply
supported
rectangular
plate
and
a
two-
opposite-edge
simply
supported plate with the other two
edges clamped, respectively. The results are shown
in Figs.2-7, and the following
observations can be made:
(i) The
present results for a homogeneous plate (p=0) are
exactly coincident with
those
obtained
from
the
classical
plate
theory,
which
verifies
the
correctness
of
the
present solution.
(ii) If the flexural rigidities Do and
Db satisfy Do
第
7
页
with the stiffness parameter p >0 is
smaller than that of the plate with P <0 (see
Figs.
2 and 5). This is because the
average bending rigidity of the plate with p >0 is
larger
than that of the plate with p
<0. It can also be seen that the maximum
deflection of the
plate is no longer at
the center of the plate when
p
≠
0.
(iii)
It
can
be
observed
from
Figs.2
and
5
that
the
deflections
of
the
homogeneous plates with
the stiffness Do and Db are the upper and lower
limits of
the deflections of a plate
with variable stiffness from Do to Db.
(iv) For different stiffness parameters
p, the variations of the bending moments
M
X
and
M
Y
with the coordinate y (at
x=a/2) shown in Figs.3, 4, 6, and 7 are similar.
However,
the
maximum
magnitudes
of
the
bending
moments
M
and
My
decrease
with
the increase of P(when p >0), and they increase
with the increase of the absolute
value
of p (when p <0). The maximum bending moments of
the plate with different
stiffness
parameters p do not occur at the same point.
第
8
页
5 Conclusions
An
exact solution of the bending problem is obtained
for a thin rectangular plate
subjected
to the distributed loads by assuming that the
flexural rigidity is of the power
form,
and
Poisson’s ratio
is a
constant. The validity of the present results is
shown by
comparing the present results
with those of the classical solution. The
influence of the
variable parameters on
the deflection and bending moments is studied by
numerical
examples. The obtained
solution
can be
used to
assess the validity
and
accuracy of
various approximate theoretical and
numerical models of plates with in-plane variable
stiffness.
References
[1] Tatting, B.
F. and Girdle, Z. Analysis and Design of Variable
Stiffness Composite Cylinders,
PhD.
dissertation, Virginia Polytechnic Institute and
State University, 1-200 (1998)
[2]
Cheng, Z. (a. and Batra, R. C. Three-dimensional
thermo elastic deformations of a functionally
graded elliptic plate. Composites: Part
B, 31, 97-106 (2000)
[3] Chakrabortya,
A., Gopalakrishnana, S
,
and
Reddy, J. N. A new beam finite element for the
analysis
of
functionally
graded
materials.
International
Journal
of
Mechanical
Sciences,
45,
519-539 (2003)
[4] Dai, K.
Y., Liu, G. R., Lim, K. M., Han, X., and Du, S.
Y
. A mesh free radial point
interpolation
method for analysis of
functionally graded material (FGM) plates.
Computational Mechanics,
34,213-223
(2004)
[5] Ferreira, A. J. M., Batra,
R. C., Roque, C. M. C., Qian, L. F., and Martins,
P. A. L. S. Static
analysis of
functionally graded plates using third-order shear
deformation theory and a meshes
method.
Composite Structures, 69, 449-457 (2005)
[6]
Yu,
T.
and
Zhong,
Z.
Vibration
of
a
simply
supported
functionally
graded
piezoelectric
rectangular
plate. Smart Materials and Strictures, 15,
1404-1412 (2006)
[7]
Shang,
E.
T.
and
Zhong,
Z.
Closed-
form
solutions
of
three-dimensional
functionally
graded
plates.
Mechanics of Advanced Materials and Structures,
15, 355-363 (2008)
[8]
Zhong,
Z.
and
Cheng,
Z.
(a.
Fracture
analysis
of
a
functionally
graded
strip
with
arbitrary
第
9
页
distributed material properties.
International Journal of Solids and Structures,
45, 3711-3725
(2008)
[9]
Chen,
S.
P.
and
Zhong,
Z.
Three-dimensional
elastic
solution
of
a
power
form
functionally
graded
rectangular plate. Journal of Mecha
epics’
and MEMS, 1(2),
349-358 (2009)
[10]
Birsan,
M
,
Altenbach,
H.,
Sadowski,
T.,
Eremeyev,
V
.
A.,
and
Pietas’
,
D.
Deformation
analysis
of
functionally
graded
beams
场
the
direct
approach.
Composites:
Part
B,
43,
1315-1328
(2012)
[11] Thai, H. T. and Choi, D. H.
An efficient and simple refined theory for
buckling analysis of
functionally
graded plates. Applied Mathematical Modeling, 36,
1008-1022 (2012)
[12] Jodaei, A.,
Jalal, M., and Yes, M. H. Free vibration analysis
of functionally graded annular
plates
be
state-space
based
differential
quadrature
method
and
comparative
modeling
be
ANN.
Composites: Part B, 43, 340-353 (2012)
[13] Wen, P. H. and Aliabad, M. H.
Analysis of functionally graded
plates
勿
mesh less method: a
purely analytical formulation.
Evgiv,eering Artalysi
。
with
Boundary Elemi arts, 36, 639-650
(2012)
[14]
Liew,
K.
M.,
Zhao,
X.,
and
Lee,
Y
.
Y.
Post
buckling
responses
of
functionally
graded
cylindrical
shells
under
axial
compression
and
thermal
loads.
Composites:
Part
B,
43,
1621-1630
(2012)
[15]
Malekzadeh,
P.,
Fiouz,
A.
R.,
and
Sahraouian,
M.
Three-dimensional
free
vibration
of
functionally graded truncated conical
shells subjected to thermal environment.
Journal of Pressure V
essels
aced Piping, 89, 210-221 (2012)
[16]
Sadeghi,
H.,
Baghani,
M.,
and
Naghdabadi,
R.
Strain
gradient
elasticity
solution
for
functionally graded
micro-cylinders. International Journal of
Engineering Science, 50, 22-30
(2012)
[17]
Zenkour,
A.
M.
Dynamical
bending
analysis
of
functionally
graded
infinite
cylinder
with
rigid core. Applied Mathematics and
Computation, 218, 8997-9006 (2012)
[18]
Shang, X. C. Exact solution on a problem of
bending of double-direction rectangular elastic
plates
of
variable
rigidity
(in
Chinese).
Journal
of
Lanzhou
University
(Natural
Sciences),
27(2), 24-32(1991)
[19]
Yang,
J.
The
structural
analysis
of
plates
with
unidirectional
varying
rigidity
on
Gale
kin
line method (in Chinese).
Journal of Wuhan, Institute of Chemical Techs
logy, 18(1), 57-60
(1996)
[20] Liu, D. Y., Wang, C.
Y
., and Chen, W. Q. Free vibration of
FGM plates with in-plane material
in
homogeneity. Composite Struct,92, 1047-1051 (2010)
[21]
Uymaz,
B.,
Aided,
M.,
and
File,
S.
Vibration
analyses
of
FGM
plates
with
in-plane
material in homogeneity by
Ritz method. Composite Structures, 94, 1398-1405
(2012)
[22] Badeshi, M. and Saudi, A.
R. Levy-type solution for buckling analysis of
thick functionally
graded rectangular plates
based on the higher-order shear deformation plate
theory. Applied
Mathematical Modeling,
34, 365
于
3673 (2010)
[23]
Thai,
H.
T.
and
Kim,
S.
E.
Levy-type
solution
for
buckling
analysis
of
orthotropic
plates
based on two variable
refined plate theories. Crocoites Structures, 93,
1738-1746 (2011)
第
10
页
-
-
-
-
-
-
-
-
-
上一篇:罗夏墨迹心理测验,你看懂了吗?一个投射法人格测验
下一篇:呼吸机相关词汇