Ce656 1
Ce656 1
Ce656 1
PLATES
• Flat plates are extensively used in many
engineering applications like roof and floor of
buildings, deck slab of bridges, foundation
footings, water tanks, bulk heads, turbine disks
etc….
• Plates are normally subjected to lateral loads
causing bending of the plate.
• The geometry of the plate is normally defined
by the middle plane which is a plane equidistant
from top and bottom faces of the plate.
• The thickness of the plate (h) is measured in
direction normal to the middle plane of the
plate.
• The flexural properties of plate depends on the
thickness.
• In general, plate problems can be classified in
to three major categories.
• i) Thin plate ii) Moderately thick plate and
• Iii) Thick plate depending upon the thickness of
the plate.
• If the thickness of the plate is very small when
compared to other dimensions, such a plate is called
thin plate. As a thumb rule,
• a/h > 20 ------- Thin plate
• 5<a/h< 20 ------ Moderately thick plate
• a/h <5 ------ Thick plate
• Where a is the least lateral dimension of the
plate.
• The theories associated with thin plate and
moderately thick plates are two dimensional
theories.
• Very thick plates have to be analysed using 3D
theories.
• Analysis of a plate by 3D theory is the most
accurate.
• Reduction of the above theory to 2D theory is
based on certain assumptions.
• 2D theories are classified as follows.
• i) Classical thin plate theory (CTPT)
• Ii) First order shear deformation theory (FSDT)
• Iii) Higher order shear deformation theory
(HSDT)
• Classical theory is used for thin plates.
• This is based on Love-Kirchhoff hypothesis
which is similar to Bernoulli’s assumptions for
thin beams.
• This is basically plane cross sections remain
plane after bending.
• (i.e.) Effect of transverse shear stresses are
neglected.
• Though for a majority of practical applications,
the CTPT yields sufficiently accurate results,
the accuracy decreases with increasing plate
thickness , with rapidly increasing loads and
also in problems of stress concentration in
plates.
• The shortcomings of thin plate theory can be
partly overcome by refined theories like shear
deformation theories.
• Basically these theories are based on
assumed displacement fields.
• The accuracy of results depend up on the
consistency of the assumed displacement field.
• It has been found that shear deformation
theories predict the accurate results only for
moderately thick plates.
• In the case of very thick plates, only three
dimensional analysis based on theory of
elasticity gives accurate results.
• For many practical problems, solutions based
on three dimensional theory of elasticity would
be very difficult to obtain.
• Only simple problems have been solved using
this theory.
• 3D solutions are used to assess the accuracy of
CTPT and shear deformation thoeries.
• Hence accurate results can be obtained using
CTPT for thin plates.
• In the case of moderately thick plates, effect of
shear deformation has to be considered.
• That is why we use FSDT and HSDT for
moderately thick plates.
• FSDT is also called Mindlin’s theory
• The presence of shear deformation is to
increase the deflection, making the plate more
flexible.
• Very thick plates have to be analysed using 3D
theory.
• Other aspects may also come in to picture
when we discuss plates.
• One aspect is the material of the plate.
• We may have isotropic plates, orthotropic
plates or anisotropic plates.
• Second aspect is the type of analysis we do.
• We may do linear analysis or nonlinear
analysis.
• Also another type of analysis is the inelastic
analysis or plastic analysis.
• Still other types analyses include dynamic analysis,
buckling analysis and post buckling analysis.
• Hence there are different types of plate analysis.
• The following schematic diagram can be used
to illustrate the different types of plate analysis.
•
Analysis of Plates
Isotropic Anisotropic
First Higher
Order Order
1
1 1
11
1 1
• Here 1 means the following part of flow chart
•
Isotropic(Small Anisotropic(Small
Anisotropic(Small
or Large or Large
or Large Deflection
Deflection) Deflection)
• In the present subject, we will linear elastic
analysis.
• Also the type of materials will be isotropic as
well as orthotropic.
• Theories adopted will be mostly Classical thin
plate theory (CTPT). Later First Order Shear
Deformation Theory will be introduced.
Basic Equations of Thin Plate
Theory
• The following fundamental assumptions are
made in the small deflection theory of
homogeneous thin plates.
• Assumptions
• i) The deflection of the plate is small compared
to thickness.
• Ii) The middle plane of the plate remains
unstrained after bending and hence will be a
neutral plane
• Iii) Plane sections initially normal to middle
plane remain plane and normal after bending.
• Iv) Transverse normal stress and strain
(through the thickness of the plate) is neglected
• The material of the plate is assumed elastic and
isotropic in the beginning(Later orthotropic and
anisotropic plates may be discussed.)
• Deflection small means around 1/10 of
thickness of the plate.
• If it is 1/5,1/2 then we have to use nonlinear
theory.(Geometric nonlinearity.)
• In this case middle plane does not remain
neutral. Stretching of midplane occurs.
Slope and Curvatures of a Bent
Plate.
• Consider a rectangular plate with x and y axes
lying in the middle plane and z axis
perpendicular downward.
• u, v and w are the displacements of a typical
point in the plate after bending.
• w is called lateral deflection and is deflection of
the middle plane(neutral plane) u=v=0 for this
plane.
• Now the plate deflection is w. This is a function
of x and y.
𝜕𝑤
• Now represents the slope of the plate in x
𝜕𝑥
direction (Gradient in x direction)
𝜕𝑤
• represents the slope of the plate in y
𝜕𝑦
direction (Gradient in y direction).
• If we consider two very near points a and b on the
plate in any direction, then we have change in
deflection between these points given by
𝜕𝑤 𝜕𝑤
• 𝑑𝑤 = 𝑑𝑥 + 𝑑y.
𝜕𝑥 𝜕𝑦
• If the corresponding direction is n, then the slope in
that direction is given by
𝜕𝑤 𝜕𝑤 𝑑𝑥 𝜕𝑤 𝑑𝑦
• 𝜕𝑛
=
𝜕𝑥 𝑑𝑛
+ 𝜕𝑦 𝑑𝑛
x
dx
𝜕𝑤
w w+ 𝜕𝑥 . 𝑑𝑥 dx
a
α
dy
dn b
n
z
z
y
• The above diagram shows the bending moments.
• Just as in the case of beams, all bending
moments cause bending normal stresses and
vice versa for plates.
• Now for the above quantities, we have
ℎ/2
• 𝑀𝑥𝑥 = −ℎ/2 σ𝑥𝑥 𝑧𝑑𝑧. (Moment/unit length)
ℎ/2
• And 𝑀𝑦𝑦 = −ℎ/2 σ𝑧𝑧 𝑧𝑑𝑧.
• These can be verified from the figure below.
• (For 𝑀𝑥𝑥 )
Neutral axis h/2 z
σ𝑥𝑥
h/2
z
• Now to get expressions for moments, we
should first obtain expressions for stresses from
strains.
• This is from stress- strain law for plane stress
case(Assumption iv)
𝐸
• We know that σ𝑥𝑥 = (ε𝑥𝑥 + υε𝑦𝑦 )
1−υ2
𝐸
• and σ𝑦𝑦 = (ε𝑦𝑦 + υε𝑥𝑥 )
1−υ2
• We assume that we have isotropic material.
• Substituting the expressions for strains, we
have
𝜕2 𝑤 𝜕2 𝑤
• σ𝑥𝑥 =
𝐸
1−υ2
𝑧
𝑟𝑥
+ υ
𝑧
𝑟𝑦
=
−𝐸𝑧
1−υ2 𝜕𝑥 2
+ υ 2
𝜕𝑦
• Similarly, we have
𝐸 𝑧 𝑧 −𝐸𝑧 𝜕2 𝑤 𝜕2 𝑤
• σ𝑦𝑦 = + υ = + υ 2
1−υ2 𝑟𝑦 𝑟𝑥 1−υ2 𝜕𝑦 2 𝜕𝑥
A
C B
𝑄𝑥
y
A
C B
y 𝑀𝑥𝑦
𝜕𝑀𝑥𝑦
𝑀𝑥𝑦 + dy
𝜕𝑦
x
• V) Beam Supported edge
• If the edge of a rectangular plate is
monolithically connected to a supporting beam,
then the deflection and rotation at the edge
would not be zero but would be equal to the
deflection and rotation of the beam.
𝑉𝑥
𝑀𝑥𝑥
x
𝑉𝑥
𝑀𝑥𝑥
𝑀2 𝑀1 x
y
• For the given problem, we have
𝑎 𝜕2 𝑤 𝜕2 𝑤
• At 𝑥 = ± , 𝑀𝑥𝑥 = −𝐷 + υ 2 = 𝑀1
2 𝜕𝑥 2 𝜕𝑦
𝑏 𝜕2 𝑤 𝜕2 𝑤
• At y= ± , 𝑀𝑦𝑦 = −𝐷 +υ 2 = 𝑀2
2 𝜕𝑦 2 𝜕𝑥
• Since it is case of pure bending, B. M at any will
be 𝑀1 and 𝑀2 causing bending in xz and yz
planes
• Hence in general we may say that for any point
in the plate,
𝜕2 𝑤 𝜕2 𝑤
• −𝐷
𝜕𝑥 2
+ υ 2
𝜕𝑦
= 𝑀1
𝜕2 𝑤 𝜕2 𝑤
• And −𝐷 + υ 2 = 𝑀2
𝜕𝑦 2 𝜕𝑥
𝜕2 𝑤 −(𝑀1 −υ𝑀2 )
• Solving, we get =
𝜕𝑥 2 𝐷(1−υ2 )
𝜕2 𝑤 −(𝑀2 −υ𝑀1 )
• And =
𝜕𝑦 2 𝐷(1−υ2 )
• Integrating the above by parts, we get
𝑀1 −𝜐𝑀2 2 𝑀2 −𝜐𝑀1
• 𝑤=− 𝑥 − 𝑦2
2𝐷 1−υ2 2𝐷 1−υ2
• Now when 𝑀1 = 𝑀2 = 𝑀, then
𝑀
• 𝑤= − (𝑥 2 +𝑦 2 )
2𝐷 1+υ
• This curved surface is a paraboloid of
revolution.
• Here both the moments are sagging(+ve)
• Hence curvatures in both directions are of the
same sign.
• The surface so generated are called synclastic
surface. (𝑀1 , 𝑀2 >0).
• When 𝑀1 = 𝑀 and 𝑀2 = −𝑀, we get
𝑀
• w= − (𝑥 2 −𝑦 2 )
2𝐷 1+υ
• Such surfaces are called anticlastic surfaces.
• Here the curvatures in both directions are
opposite in sign.
Cylindrical Bending of Plates
• Consider a long rectangular plate which is
subjected to a lateral load q(x) which does not
vary in y direction (as shown in the figure)
a
x
Unit width
y
SS
SS
a
x
b SS
SS
y
• For all sides simply supported rectangular
plates, there is a method based on double
Fourier Series.
• This is called Navier’s method.
• Here we expand the loading as well as the
deflection in a double Fourier sine series.
• The basis of this is that the deflection satisfies
the boundary condition that it is =0 on the
boundary. Also the bending moment is zero on
the boundary. Hence SS condition is satisfied.
Fourier Series-Basics
• Any integrable periodic function can be
expanded in terms of a Fourier sine series or
cosine series.
• We will be seeing Fourier sine series
• Consider a function f(x, y).
• In our case, this may represent deflection or
applied load.
• This can be expanded in terms of a Fourier
Sine series as follows
∞ ∞ 𝑚π𝑥 𝑛π𝑦
• 𝑓 𝑥, 𝑦 = σ𝑚=1 σ𝑛=1 𝑓𝑚𝑛 𝑠𝑖𝑛 𝑠𝑖𝑛
𝑎 𝑏
• Here a and b are the sides of the rectangular
plate we analyse.
• 𝑓𝑚𝑛 are called the Fourier coefficients of the
function f(x, y)
• Expression for 𝑓𝑚𝑛 are derived from the
orthogonal properties of sine functions.
• Multiplying the above equation by
𝑚π𝑥 𝑛π𝑦
s𝑖𝑛 𝑠𝑖𝑛 and integrating , we get
𝑎 𝑏
𝑎 2 𝑚π𝑥 𝑏 2 𝑛π𝑦
• 𝑓𝑚𝑛 =𝑥0 𝑠𝑖𝑛 dx. =𝑦0 𝑠𝑖𝑛 dy
𝑎 𝑏
𝑎 𝑏 𝑛π𝑥 𝑛π𝑦
• = =𝑥0 =𝑦0 f(x, y)𝑠𝑖𝑛 𝑠𝑖𝑛 dx dy
𝑎 𝑏
• (i.e.)
𝑎 𝑏 𝑎 𝑏 𝑛π𝑥 𝑛π𝑦
• 𝑓𝑚𝑛 .
2 2
= =𝑥0 =𝑦0 f(x, y)𝑠𝑖𝑛
𝑎
𝑠𝑖𝑛
𝑏
dx dy
• (Or) we get the Fourier coefficient of the
function as
4 𝑎 𝑏 𝑛π𝑥 𝑛π𝑦
• 𝑓𝑚𝑛 = 𝑥=0 𝑦=0
f(x, y)𝑠𝑖𝑛 𝑠𝑖𝑛 dx dy
𝑎𝑏 𝑎 𝑏
• Hence the Fourier coefficient 𝑓𝑚𝑛 depends on
the function f(x,y).
• This function can be anything like
• c ,x, y, xy, 𝑥 2 etc…. Where c is a constant.
• Hence 𝑤 𝑥, 𝑦 = σ∞ σ∞
𝑚=1 𝑛=1 𝐴𝑚𝑛 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
𝑚π 𝑛π
• Where α𝑚 = and β𝑛 =
𝑎 𝑏
• And 𝐴𝑚𝑛 is a constant(Fourier coefficient)
• This form of deflection satisfies the simply
supported boundary conditions.
• The lateral load q(x,y) is also expressed in
terms of a double Fourier series as
• q 𝑥, 𝑦 = σ∞ σ ∞
𝑚=1 𝑛=1 𝑞𝑚𝑛 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
• Here 𝑞𝑚𝑛 is the Fourier coefficient of the load.
• Substituting the above expressions in the thin
plate equation equation (5), we get
𝑞𝑚𝑛
• σ∞ σ∞
[𝐴 (α
𝑚=1 𝑛=1 𝑚𝑛 𝑚
2 + β2𝑛 )2 − ]𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦 = 0.
𝐷
• The above is valid for all values of x and y.
• Hence we can see that
2 2 𝑞𝑚𝑛
• [𝐴𝑚𝑛 (α2𝑚 + β𝑛 ) − ] =0
𝐷
𝑞𝑚𝑛
• Thus 𝐴𝑚𝑛 =
𝐷(α2𝑚 +β2𝑛 )2
• Hence the expression for deflection becomes
• 𝑤 𝑥, 𝑦 = σ∞ σ ∞
𝑚=1 𝑛=1 𝐴𝑚𝑛 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
∞ ∞ 𝑞𝑚𝑛
• = σ𝑚=1 σ𝑛=1 2 2 2 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
𝐷(α +β ) 𝑚 𝑛
b/2 x
y
a b/2
a
b
η x
ξ
y x
𝑞0 𝑥
• Here load intensity is given by 𝑞 = .
𝑎
• The Fourier coefficient of the load is given
below.
4
• 𝑞𝑚𝑛 = 𝑛𝑖𝑠 𝑞α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦 dx dy
𝑎𝑏
4 𝑞0 𝑥
• = 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦 dx dy
𝑎𝑏 𝑎
4𝑞0 𝑎 𝑏
• = 2 =𝑥0 𝑥𝑠𝑖𝑛α𝑚 𝑥 𝑑𝑥 =𝑦0 𝑠𝑖𝑛β𝑛 𝑦𝑑𝑦
𝑎 𝑏
4𝑞0 (1−cosnπ) 𝒂
• = 2 𝒙=𝟎
𝒙𝒔𝒊𝒏𝜶𝒎 𝒙 𝒅𝒙
𝑎 𝑏 β𝑛
4𝑞0 (1−cosnπ) −𝑎 𝑐𝑜𝑠𝑚π
• = 2 ( ) Integration
𝑎 𝑏 β𝑛 α𝑚 by parts
8𝑞0 𝑚+1
• = 2 (−1) , m=1,2,3,..
𝑚𝑛π
• n=1,3,5…..
• Hence the expression for deflection is given as
∞ ∞ 𝑞𝑚𝑛
• 𝑤 𝑥, 𝑦 = σ𝑚=1 σ𝑛=1 2 2 2 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
𝐷(α +β )𝑚 𝑛
8𝑞0 𝑚+1
(−1)
∞ 𝑚𝑛π2
• = σ∞
𝑚=1 σ 𝑛=1 𝐷(α2 +β2 )2 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
𝑚 𝑛
8𝑞0 ∞ (−1)𝑚+1
• =
π 2
∞
σ𝑚=1,2,3 σ𝑛=1,3,5 2 2
𝐷𝑚𝑛(α𝑚 +β𝑛 ) 2 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
Influence Function for SS Plates
• Influence function is a more general form of
influence lines.
• For a simply supported plate, if we apply a unit
load at some point and measure the deflection
at some other point ,that can be taken as
influence function.
• This we have already discussed
• Earlier case, we put P=1 to get the following
expression for influence function.
4 ∞ 𝑠𝑖𝑛α𝑚 ξ 𝑠𝑖𝑛β𝑛 η
• G 𝑥, 𝑦, ξ, η = σ∞ σ 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
𝑎𝑏𝐷 𝑚=1 𝑛=1 (α2𝑚 +β2𝑛 )2
• =
𝟒𝑴 ∞
𝒂𝒃𝑫
σ𝒎=𝟏 σ∞
𝒏=𝟏
𝜶𝒎 𝒄𝒐𝒔𝜶𝒎 𝝃 𝒔𝒊𝒏𝜷𝒏 𝜼
𝟐 𝟐 𝟐 𝒔𝒊𝒏𝜶𝒎 𝒙 𝒔𝒊𝒏𝜷𝒏 𝒚
𝜶 +𝜷𝒎 𝒏
y
• Now the deflection surface is taken as
π𝑥 π𝑦
• w 𝑥, 𝑦 = 𝑤0 𝑠𝑖𝑛 𝑠𝑖𝑛
𝑎 𝑏
• This satisfies the SS boundary conditions.
• Substituting in the plate equation, we get
π4 π4 π4 π𝑥 π𝑦 π𝑥 π𝑦
• 𝐷 + 2 2 2 + 𝑤0 𝑠𝑖𝑛 𝑠𝑖𝑛 = 𝑞0 𝑠𝑖𝑛 𝑠𝑖𝑛
𝑎4 𝑎 𝑏 𝑏4 𝑎 𝑏 𝑎 𝑏
𝑞0 𝑞0
• Hence 𝑤0 = π4 π4 π4
= π2 π2 2
𝐷 +2 2 2 + 4 𝐷( 2 + 2 )
𝑎4 𝑎 𝑏 𝑏 𝑎 𝑏
𝜕2 𝑤 𝜕2 𝑤 𝑞0 π2 π2 π𝑥 π𝑦
• 𝑀𝑦𝑦 = −𝐷 + υ 2 = 2 (υ 2 + )𝑠𝑖𝑛 𝑠𝑖𝑛 .
𝜕𝑦 2 𝜕𝑥 π 2 π 2 𝑎 𝑏2 𝑎 𝑏
2 + 2
𝑎 𝑏
• And
𝜕2 𝑤
• 𝑀𝑥𝑦 = −D(1 − υ)
𝜕𝑥𝜕𝑦
𝑞0 1−𝜐 π2 π𝑥 π𝑦
• = 2 𝑐𝑜𝑠 𝑐𝑜𝑠
π2 π2 𝑎𝑏 𝑎 𝑏
𝑎2+ 2
𝑏
𝜕(𝛻2 𝑤)
• Now shear force 𝑄𝑥 = −𝐷
𝜕𝑥
𝜕3 𝑤 𝜕3 𝑤
• = −D( 3 + 2 )
𝜕𝑥 𝜕𝑥𝜕𝑦
π3 π3
𝑞0 ( 3 + 2 ) π𝑥 π𝑦 𝑞0 π π𝑥 π𝑦
• = 𝑎
π2 π2 2
𝑎𝑏
𝑐𝑜𝑠 𝑠𝑖𝑛
𝑎 𝑏
= π2 π2
𝑐𝑜𝑠 𝑠𝑖𝑛
𝑎 𝑏
( 2+ 2) 𝑎( + )
𝑎 𝑏 𝑎2 𝑏2
• For a square plate a=b, the bending moments
• and twisting moment at centre will be
(x=a/2,y=a/2)
𝑞0 π2 π2 π π
• 𝑀𝑥𝑥 = 2 2 2( 2 + υ 2 )𝑠𝑖𝑛 𝑠𝑖𝑛 .
π π 𝑎 𝑎 2 2
2 + 2
𝑎 𝑎
𝑞0 (1+υ)𝑎2
• = = 0.1317𝑞0 𝑎2 (assuming υ=0.3)
π2
• By symmetry 𝑀𝑥𝑥 = 𝑀𝑦𝑦 at the centre.
• Also by symmetry, 𝑀𝑥𝑦 = 0 at the centre.
• Now let’s find the support reactions for the
square plate.
• Support reactions consists of side reactions in
the form of distributed forces and corner
reactions in the form of concentrated forces.
• The side reactions are nothing but the Kirchhoff
shear 𝑉𝑥 and 𝑉𝑦 .
• Now we have already got expressions for
Kirchhoff reactions as
𝜕 𝜕2 𝑤 𝜕2 𝑤
• 𝑉𝑥 = −𝐷 [ 2 + (2 − υ) 2 ]
𝜕𝑥 𝜕𝑥 𝜕𝑦
π3 π3
𝑞0 ( 3 +(2−υ) 2 ) π𝑥 π𝑦
• = 𝑎
π2 π2 2
𝑎𝑏
𝑐𝑜𝑠 𝑠𝑖𝑛 (for a=b)
( 2+ 2) 𝑎 𝑏
𝑎 𝑏
π3 π3
𝑞0 ( 3 +(2−υ) 3 ) π𝑥 π𝑦
𝑎 𝑎
• 𝑉𝑥 = π2 π2 2
𝑐𝑜𝑠 𝑠𝑖𝑛
( 2+ 2) 𝑎 𝑎
𝑎 𝑎
𝑞0 𝑎(3−υ) π𝑥 π𝑦 π𝑥 π𝑦
• = 𝑐𝑜𝑠 𝑠𝑖𝑛 = 0.2148𝑞0 𝑎𝑐𝑜𝑠 𝑠𝑖𝑛
4π 𝑎 𝑎 𝑎 𝑎
π𝑦 π𝑥
• Because of symmetry 𝑉𝑦 = 0.2148𝑞0 𝑎𝑐𝑜𝑠 𝑠𝑖𝑛
𝑎 𝑎
• Coming to the corner reaction, all the corner
reactions will be equal because of symmetry.
• All of them R = 2𝑀𝑥𝑦 (at corner say x=0,y=0)
2𝑞0 1−𝜐 π2 π𝑥 π𝑦
• = 2 𝑐𝑜𝑠 𝑐𝑜𝑠 (a=b)
π2 π2 𝑎𝑏 𝑎 𝑏
𝑎2+ 2 𝑏
𝑞0 𝑎2 1−𝜐
• = = 0.0355𝑞0
2π2
• 𝑉𝑥 and 𝑉𝑦 will be upward and R will be
downward (as shown below).
π𝑥 π𝑦
𝑞 = 𝑞0 𝑠𝑖𝑛 𝑠𝑖𝑛
𝑎 𝑏
x
R
𝑉𝑥
y R 𝑉𝑦 R
SS b/2
a
y b/2 x
𝜕2 𝑤
• The simply supported B.C are 𝑤 = =0
𝜕𝑥 2
• The solution to equation a can be taken in the
following form
• 𝑤0 = σ∞ 𝑚=1 𝑓𝑚 𝑦 𝑠𝑖𝑛α𝑚 𝑥
𝑚π
• Where 𝑓𝑚 𝑦 is a function of y only. α𝑚 =
𝑎
• The assumed solution satisfies the simply
supported boundary condition.
• Substituting the above in equation (a), we get
the following ordinary differential equation.
𝑑 2𝑓 𝑑 4 𝑓𝑚
• σ∞ [α4
𝑚=1 𝑚 𝑓𝑚 −
2
2α𝑚 2 𝑚
+ 4 ]𝑠𝑖𝑛α𝑚 𝑥 =0
𝑑𝑦 𝑑𝑦
• This should be satisfied by all x.
𝑑 2 𝑓𝑚 𝑑 4 𝑓𝑚
• Hence we have α4𝑚 𝑓𝑚 − 2
2α𝑚 2 + =0
𝑑𝑦 𝑑𝑦 4
• The solution to the ODE can be taken as
• 𝑓𝑚 (𝑦) = 𝐴𝑒 λ𝑚𝑦
• Substituting this, we get the following algebraic
equation.
• α4𝑚 − 2α2𝑚 λ2𝑚 +λ4𝑚 = 0
• (Or) (α2𝑚 − λ2𝑚 )2 = 0
• Here we have real and repeated roots.
• The roots are λ𝑚 = ±α𝑚 occurring twice.
• Since we have repeated roots, the solution will
be 𝑒 ±α𝑚𝑦 and y𝑒 ±α𝑚𝑦 . Since exponentials can
be written using hyperbolic sine and cosine
forms, we have the following solution for 𝑓𝑚 (y).
• 𝑓𝑚 𝑦
• = (𝐴𝑚 𝑐𝑜𝑠ℎα𝑚 𝑦 + 𝐵𝑚 𝑠𝑖𝑛ℎα𝑚 𝑦 +
• 𝐶𝑚 α𝑚 𝑦𝑠𝑖𝑛ℎα𝑚 𝑦 + 𝐷𝑚 α𝑚 𝑦 𝑐𝑜𝑠ℎα𝑚 𝑦)
• Hence we get the homogeneous solution as
• 𝑤0 =
σ∞
𝑚=1(𝐴𝑚 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + 𝐵𝑚 𝑠𝑖𝑛ℎ𝛼𝑚 𝑦 + 𝐶𝑚 𝛼𝑚 𝑦𝑠𝑖𝑛ℎ𝛼𝑚 𝑦
+ 𝐷𝑚 𝛼𝑚 𝑦 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦)𝑠𝑖𝑛𝛼𝑚 𝑥
• Next we have to get the particular solution.
• This depends on the type of loading present.
• The loading is a function of x and y. This is also expressed in
terms of a single Fourier series.
σ ∞
• 𝑞 𝑥, 𝑦 = 𝑚=1 𝑞𝑚 𝑦 𝑠𝑖𝑛α𝑚 𝑥
2 𝑎
• Where 𝑞𝑚 𝑦 = =𝑥0 𝑞 𝑥, 𝑦 𝑠𝑖𝑛(α𝑚 𝑥)𝑑𝑥
𝑎
• Substituting in the governing D.E for plates, we
get the following equation.
𝑑 2𝑝 𝑑 4 𝑝𝑚
• σ∞ [α4
𝑚=1 𝑚 𝑝𝑚 −
2
2α𝑚 2 𝑚
+ 4 ]𝑠𝑖𝑛α𝑚 𝑥
𝑑𝑦 𝑑𝑦
σ∞
𝑚=1 𝑞𝑚 𝑦 𝑠𝑖𝑛α𝑚 𝑥
• =
𝐷
• Since the above equation has to be satisfied for
all values of x, we have
𝑑 2𝑝 𝑑 4 𝑝𝑚 𝑞𝑚
• α4𝑚 𝑝𝑚 − 2
2α𝑚 2 𝑚
+ = .
𝑑𝑦 𝑑𝑦 4 𝐷
• By determining 𝑝𝑚 ,satisfying the above ODE,
the particular integral 𝑤1 can be determined.
• Thus we will be able to get the complete solution.
• Sometimes the Levy’s solution can be used to
consider the variation in loading in x direction
only.
• In that case, the particular integral can be
𝑞𝑚
written as 𝑝𝑚 = 4 where 𝑞𝑚 is the Fourier
α𝑚 𝐷
coefficient of the load.
• The total solution for Levy’s method is given by
the following equation.
• 𝑤=
(𝐴𝑚 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + 𝐵𝑚 𝑠𝑖𝑛ℎ𝛼𝑚 𝑦 + 𝐶𝑚 𝛼𝑚 𝑦𝑠𝑖𝑛ℎ𝛼𝑚 𝑦
σ∞ 𝑞 𝑠𝑖𝑛𝛼𝑚 𝑥 .
𝑚=1 +𝐷𝑚 𝛼𝑚 𝑦 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + 4𝑚 )
α𝑚 𝐷
x
ξ
𝟐 ξ+𝒄 4𝑞
• 𝒒𝒎 = ξ−𝒄
𝑞 𝑠𝑖𝑛α𝑚 𝑥 𝑑𝑥 = 𝑠𝑖𝑛α𝑚 ξ 𝑠𝑖𝑛α𝑚 𝑐
𝒂 𝑚π
• m=1,2,3….
• d) Line load p per unit length
p
ξ
x
4𝑞
• 𝑞𝑚 = 𝐿𝑖𝑚𝑖𝑡𝑐 0 𝑚π 𝑠𝑖𝑛α𝑚 ξ𝑠𝑖𝑛α𝑚 𝑐
4𝑞
• = 𝑠𝑖𝑛α𝑚 ξ. α𝑚 𝑐
𝑚π
4𝑞 𝑚π 2𝑝
• = 𝑠𝑖𝑛α𝑚 ξ. 𝑐 = 𝑠𝑖𝑛α𝑚 ξ,
𝑚π 𝑎 𝑎
• m=1,2,3..
• (Here we have put 𝑝 = 2𝑞𝑐)
• Now let’s take some examples.
• 1) As a first example, let’s take an all edges simply
supported plate subjected to udl. This has been
discussed in Navier’s method.
• Now the starting solution for Levy’s method is given as
follows
• 𝑤=
(𝐴𝑚 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + 𝐵𝑚 𝑠𝑖𝑛ℎ𝛼𝑚 𝑦 +
σ∞𝑚=1 𝐶 𝛼 𝑦𝑠𝑖𝑛ℎ𝛼 𝑦 + 𝐷 𝛼 𝑦 𝑐𝑜𝑠ℎ𝛼 𝑦 + 𝑞𝑚 )𝑠𝑖𝑛𝛼 𝑥
𝑚 𝑚 𝑚 𝑚 𝑚 𝑚 4 𝑚
α𝑚 𝐷
4𝑞
• Here for a udl, we have 𝑞𝑚 = , m=1,3,5..
𝑚π
𝑚π
• Also α𝑚 =
𝑎
• Hence 𝑤=
(𝐴𝑚 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + 𝐵𝑚 𝑠𝑖𝑛ℎ𝛼𝑚 𝑦
σ∞
𝑚=1,3,5 4𝑞𝑎4
+ 𝐶𝑚 𝛼𝑚 𝑦𝑠𝑖𝑛ℎ𝛼𝑚 𝑦 + 𝐷𝑚 𝛼𝑚 𝑦 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + )𝑠𝑖𝑛𝛼𝑚 𝑥.
𝑚 5 π5 𝐷
• Now considering symmetry in y direction,
• We can see that odd terms of the deflection function w.r.t y,
should vanish.
• Hence 𝐵𝑚 = 𝐷𝑚 = 0 in the above expression for deflection.
• Hence the expression for deflection for this case would be
• 𝑤(𝑥, 𝑦) =
4𝑞𝑎4
σ∞
𝑚=1,3,5(𝐴𝑚 𝑐𝑜𝑠ℎ𝛼𝑚 𝑦 + 𝐶𝑚 𝛼𝑚 𝑦𝑠𝑖𝑛ℎ𝛼𝑚 𝑦 + 5 5 )𝑠𝑖𝑛𝛼𝑚 𝑥
𝑚 π 𝐷
• The above deflection function satisfy the following
equations
• a) Plate differential equation
• b) The boundary condition at x=0 and x=a.
• Now it has to satisfy the B.C on other two edges
𝑏
• Namely y= ± .
2
• For the present case, these edges are simple
supports.
• Hence the conditions at these supports is
𝑏 𝜕2 𝑤
• Along y= ± , 𝑤= =0
2 𝜕𝑦 2
• Using these conditions, we get the following two
equations
α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 4𝑞𝑎4
• 𝐴𝑚 𝑐𝑜𝑠ℎ + 𝐶𝑚 sinh + 5 5 =0 and
2 2 2 𝑚 π 𝐷
α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏
• 𝐴𝑚 𝑐𝑜𝑠ℎ + 𝐶𝑚 (2𝑐𝑜𝑠ℎ + 𝑠𝑖𝑛ℎ ) =0
2 2 2 2
Some Relations for Hyperbolic
𝑦 −𝑦
Functions 𝑦 −𝑦
𝑒 +𝑒 𝑒 −𝑒
• 𝐶𝑜𝑠ℎ 𝑦 = , Sinℎ 𝑦 =
2 2
𝑑 𝑒 𝑦 −𝑒 −𝑦
• 𝐶𝑜𝑠ℎ𝑦 = = 𝑆𝑖𝑛ℎ𝑦
𝑑𝑦 2
𝑑 𝑒 𝑦 +𝑒 −𝑦
• 𝑆𝑖𝑛ℎ𝑦 = = 𝐶𝑜𝑠ℎ𝑦
𝑑𝑦 2
• 2𝑆𝑖𝑛ℎ𝑦. 𝐶𝑜𝑠ℎ𝑦 = 𝑆𝑖𝑛ℎ2𝑦
• 𝐶𝑜𝑠ℎ2 𝑦 − 𝑆𝑖𝑛ℎ2 𝑦 = 1
• 𝐶𝑜𝑠ℎ2 𝑦 + 𝑆𝑖𝑛ℎ2 𝑦 = 𝐶𝑜𝑠ℎ2𝑦
• Solving the two equations, we get
𝛼 𝑏
𝑞𝑎4 4+𝛼𝑚 𝑏𝑡𝑎𝑛ℎ 𝑚
• 𝐴𝑚 = − [ 2
𝛼𝑚 𝑏 ]
𝐷 5 5
𝑚 𝜋 𝑐𝑜𝑠ℎ
2
𝑞𝑎4 2
• And 𝐶𝑚 = [ 5 5 𝛼𝑚 𝑏 ]
𝐷 𝑚 𝜋 𝑐𝑜𝑠ℎ
2
Clamped
a SS
b/2
b/2
α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 4𝑞𝑎4
• 𝐴𝑚 𝑐𝑜𝑠ℎ + 𝐶𝑚 sinh + 5 5 =0
2 2 2 𝑚 π 𝐷
α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏
• 𝐴𝑚 𝑠𝑖𝑛ℎ + 𝐶𝑚 sinh + cosh =0
2 2 2 2
• Solving the above equations, we get
2 𝛼𝑚 𝑏
4𝑞𝑎4 1 𝛼𝑚 𝑏 𝑠𝑖𝑛ℎ
• 𝐴𝑚 = − 𝛼𝑚 𝑏 [ 2
+ 1]
𝐷 5 5
𝑚 𝜋 𝑐𝑜𝑠ℎ 𝛼𝑚 𝑏+𝑠𝑖𝑛ℎα𝑚 𝑏
2
𝛼𝑚 𝑏
4𝑞𝑎4 2𝑠𝑖𝑛ℎ
• And 𝐶𝑚 = [ 2
]
𝐷𝑚5 𝜋5 𝛼𝑚 𝑏+𝑠𝑖𝑛ℎα𝑚 𝑏
• Hence the expression for deflection is obtained
in the following form.
𝛼 𝑏
4𝑞𝑎4 2𝑠𝑖𝑛ℎ 𝑚
• 𝑤= σ∞
𝑚=1,3,5{
2
[α𝑚 𝑦𝑠𝑖𝑛ℎα𝑚 𝑦
𝐷𝜋5 𝛼𝑚 𝑏+𝑠𝑖𝑛ℎα𝑚 𝑏
𝛼𝑚 𝑏 𝛼𝑚 𝑏 cosh α𝑚𝑦
• − tanh cosh α𝑚 𝑦] − 𝛼𝑚 𝑏 + 1}
2 2 cosh 2
𝑠𝑖𝑛α𝑚 𝑥
•
𝑚5
• The maximum deflection occur at the centre of the
plate (x=a/2, y=0).
• The expression for maximum deflection can be written
as
5𝑞𝑎4
• 𝑤𝑚𝑎𝑥 = −
384𝐷
𝑚−1 𝛼 𝑏
4𝑞𝑎4 −1 2 2𝑠𝑖𝑛ℎ 𝑚 1
σ∞
𝑚=1,3,5 [ 2
α𝑚 𝑏𝑡𝑎𝑛ℎα𝑚 𝑏 + ]
𝐷𝜋5 𝑚5 𝛼𝑚 𝑏+𝑠𝑖𝑛ℎα𝑚 𝑏 𝑐𝑜𝑠ℎα𝑚 𝑏
𝑞𝑎4 4
• 𝑤𝑚𝑎𝑥= 0.01302 − 5 0.8499 − 0.0004217
𝐷 π
𝑞𝑎4
• = 0.001917 (Exact answer for a square
𝐷
plate)
• As far as moments are concerned, max +ve
moment occur at the centre of the
plate(x=a/2,y=0) and maximum –ve moments
occur at the centre of the fixed
edge(x=a/2,y=b/2).
• Now let’s discuss the moments for this problem.
• The moment expressions can be shown to be
as below.
α 𝑏
4𝑞𝑎2 2𝑠𝑖𝑛 𝑚
• 𝑀𝑥𝑥 = σ∞
𝑚=1,3,5[(
2
){(1 − ν)α𝑚 𝑦
π3 α𝑚 𝑏+𝑠𝑖𝑛ℎα𝑚 𝑏
α𝑚 𝑏 α𝑚 𝑏
• 𝑠𝑖𝑛ℎα𝑚 𝑦 − 𝑐𝑜𝑠ℎα𝑚 𝑦(2ν + 1 − ν 𝑡𝑎𝑛ℎ )} −
2 2
𝑐𝑜𝑠ℎα𝑚 𝑦 𝑠𝑖𝑛α𝑚 𝑥
(1 − ν) α 𝑏 + 1]
𝑐𝑜𝑠ℎ 𝑚 𝑚3
2
α 𝑏
4𝑞𝑎2 2𝑠𝑖𝑛 𝑚
• 𝑀𝑦𝑦 = σ∞
𝑚=1,3,5[(
2
){(ν − 1)α𝑚 𝑦
π3 α𝑚 𝑏+𝑠𝑖𝑛ℎα𝑚 𝑏
α𝑚 𝑏 α𝑚 𝑏
• 𝑠𝑖𝑛ℎα𝑚 𝑦 − 𝑐𝑜𝑠ℎα𝑚 𝑦(2 − 1 − ν 𝑡𝑎𝑛ℎ )} + (1 −
2 2
𝑐𝑜𝑠ℎα𝑚 𝑦 𝑠𝑖𝑛α𝑚 𝑥
ν) α 𝑏 + ν]
𝑐𝑜𝑠ℎ 𝑚 𝑚3
2
a/2
SS b/2
SS
y
x
𝑀1 = 𝑓1 (𝑥)
𝑀2 = 𝑓2 (𝑥)
y
• The boundary conditions for the problem can be
written as
𝑎 𝜕2 𝑤
• a) Along 𝑥 =± , 𝑤= 2 =0
2 𝜕𝑥
𝑏
• b) Along y = ± , w=0
2
𝑏 𝜕2 𝑤
• c) Along y = − , −𝐷 2 = 𝑀1
2 𝜕𝑦
𝑏 𝜕2 𝑤
• d) Along y = , −𝐷 2 = 𝑀2
2 𝜕𝑦
• Where 𝑀1 and 𝑀2 are the applied moments
along the edges
• Now In general 𝑀1 ≠ 𝑀2
• But they can be split in to i) moments symmetric
about x axis and ii) moments antisymmetric
about x axis.
• Now since the applied udl=0, the general
solution for the homogeneous plate equation is
the solution for deflection.
σ ∞
• 𝑤 = 𝑚=1,3,5(𝐴𝑚 cosh α𝑚 𝑦 + 𝐵𝑚 sinh α𝑚 𝑦
• +𝐶𝑚 α𝑚 𝑦 sinh α𝑚 𝑦 + 𝐷𝑚 α𝑚 𝑦 𝑐𝑜𝑠 α𝑚 𝑦)𝑐𝑜𝑠α𝑚 𝑥
• If a plate is subjected symmetric edge
moments is considered,𝐵𝑚 and 𝐷𝑚 would be
zero.
• If it is anti-symmetric moment, then 𝐴𝑚 and 𝐶𝑚
would be zero.
(𝑀1 +𝑀2 )
• Consider a symmetric moment 𝑀0 =
2
𝑏
acting along the edges 𝑦 = ± .
2
• Then 𝑀0 can be expanded in Fourier cosine
series as
• 𝑀0 = σ∞
𝑚=1,3,5 𝑓𝑚 cos α𝑚 𝑥
𝑎 𝑚−1
2 4𝑀0
• Now 𝑓𝑚 = 2
𝑎 𝑀0 cos α𝑚 𝑥 𝑑𝑥 = (−1) 2
𝑎 − 𝑚π
2
• m=1,3,5…
• The solution for the present problem is
σ ∞
• 𝑤 = 𝑚=1,3,5(𝐴𝑚 cosh α𝑚 𝑦 +𝐶𝑚 α𝑚 𝑦 sinh α𝑚 𝑦)
• 𝑐𝑜𝑠α𝑚 𝑥
• The BCs for the present case are
𝑏 𝜕2 𝑤
• Along 𝑦 = ± w=0 and −𝐷 2 = 𝑀0 =
2 𝜕𝑦
σ∞
𝑚=1,3,5 𝑓𝑚 cos α𝑚 𝑥
• From the expression for w as a series shown
above, we can put in the BCs and get two
equations to solve for 𝐴𝑚 and 𝐶𝑚 .
• They are as follows.
α𝑚 𝑏 α 𝑏
𝑓𝑚 𝑡𝑎𝑛ℎ 𝑚 𝑓𝑚
2 2
• 𝐴𝑚 = 2 α𝑚 𝑏 and 𝐶𝑚 = − 2 α𝑚 𝑏
2𝐷α𝑚 𝑐𝑜𝑠ℎ 2𝐷α𝑚 𝑐𝑜𝑠ℎ
2 2
• −α𝑚 𝑦sinhα𝑚 𝑦]
𝑚−1
4𝑀0
• We now can substitute 𝑓𝑚 = (−1) 2
𝑚π
• Hence we get the deflection as follows.
• 𝑤=
𝑚−1
2𝑀0 𝑎2 ∞ (−1) 2 α𝑚 𝑏 α𝑚 𝑏
3
σ𝑚=1,3,5 α𝑚 𝑏 [ 𝑡𝑎𝑛ℎ 𝑐𝑜𝑠ℎα𝑚 𝑦 −
π 𝐷 cosh 2 2
2
• α𝑚 𝑦sinhα𝑚 𝑦] 𝑐𝑜𝑠α𝑚 𝑥
• The above result is useful in solving the plate
problem using method of superposition.
Method of Superposition
• This is a method to solve plates with arbitrary
boundary conditions.
• Consider a rectangular plate with one pair of
Opposite edges SS and the other pair fixed.
𝑎
• Let the plate be SS on the edges 𝑥 = ± and
2
𝑏
let it be fixed along the edges 𝑦 = ± where a
2
and b are the sides and the origin of
coordinates be at the centre of the plate.
• Let the plate be subjected to a udl of intensity q.
• The solution to this problem can be obtained
from two parts.
• In the first part, we consider the the plate to be
SS and subjected to udl.
• In the second, part we consider an all sided SS
plate subjected to bending moment applied
𝑏
along 𝑦 = ± , which when superposed on the
2
first part eliminates the rotation of the edges 𝑦 =
𝑏
± due to lateral loads .
2
• This is explained below.
• The solution for the first part(SS plate subjected
to udl) has already been done earlier and the
result is given by the following expression.
𝑞𝑎4 ∞ 1 α𝑚 𝑏
• 𝑤1 = σ [ {− 4 + α𝑚 𝑏𝑡𝑎𝑛ℎ cos α𝑚 ℎ𝑦 +
π5 𝐷 𝑚=1,3,5 𝑚5 coshα𝑚 𝑏 2
2
𝑚−1
4
• α𝑚 𝑦𝑠𝑖𝑛α𝑚 ℎ𝑦} + ](−1) 2 𝑐𝑜𝑠α𝑚 𝑥.
𝑚5
𝑏
• The slope at the edges 𝑦 = ± is obtained as
2
𝑚−1
𝜕𝑤1 2𝑞𝑎3 (−1) 2 𝑐𝑜𝑠α𝑚 𝑥 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏
• = σ∞
𝑚=1,3,5 [ − 𝑡𝑎𝑛ℎ (1+ tan )]
𝜕𝑦 𝐷π4 𝑚4 2 2 2 2
• Next we consider the solution for the SS plate subjected to end
𝑏
moments along 𝑦 = ±
2
• The solution has been done earlier for this and the result is
given below.
𝑎2 ∞ 𝑓𝑚 𝑐𝑜𝑠α𝑚 𝑥 α𝑚 𝑏 α𝑚 𝑏
• 𝑤2 = σ [ 𝑡𝑎𝑛ℎ 𝑐𝑜𝑠ℎα𝑚 𝑦
2π2 𝐷 𝑚=1,3,5 𝑚2 𝑐𝑜𝑠ℎα𝑚 𝑏 2 2
2
• −α𝑚 𝑦sinhα𝑚 𝑦]
𝑏
• Now, the slope along 𝑦 = ± is given by
2
𝜕𝑤2 −𝑎 𝑓𝑚 𝑐𝑜𝑠α𝑚 𝑥 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏
• = σ∞ [ + 𝑡𝑎𝑛ℎ (1 − 𝑡𝑎𝑛ℎ )]
𝜕𝑦 2π𝐷 𝑚=1,3,5 𝑚 2 2 2 2
𝜕𝑤1 𝜕𝑤2
• Now, we have the condition that + =0
𝜕𝑦 𝜕𝑦
𝑏
along 𝑦 = ± .
2
• From this condition, upon substituting the above
expressions, we get the expression for 𝑓𝑚 as
follows.
𝑚−1 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏
4𝑞𝑎 2 −𝑡𝑎𝑛ℎ (1+ 𝑡𝑎𝑛ℎ )
2 2 2 2
• 𝑓𝑚 = (−1) 2 [ ]
𝑚3 π3 α𝑚 𝑏+𝑡𝑎𝑛ℎα𝑚 𝑏(1−α𝑚 𝑏𝑡𝑎𝑛ℎα𝑚 𝑏)
2 2 2 2
𝑚−1 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏 α𝑚 𝑏
1 −𝑡𝑎𝑛ℎ (1+ 𝑡𝑎𝑛ℎ )
2 2 2 2
• −α𝑚 𝑦sinhα𝑚 𝑦] (−1) 2 [ ]
𝑚5 α𝑚 𝑏+𝑡𝑎𝑛ℎ α𝑚 𝑏(1−α𝑚 𝑏𝑡𝑎𝑛ℎ α𝑚 𝑏)
2 2 2 2
q
y
1 2
y y y
q
• The bending of each span can be analysed
using Levy’s method by superposing laterally
loaded plates with plates subjected to end
moments.
• Since the structure and loading are symmetric,
it is only necessary to consider only one half the
plate(Hence only one intermediate support)
• Consider support 22.
• Now plate 2 will be subjected to lateral load q
and unknown bending moment along edge 22.
• Plate 3 will be subjected to only unknown
bending moment along edge 22.
• All edges are simply supported.
• Free Body Diagrams of Plates 2 and 3
Plate 2 Plate 3
b x
x
𝑀𝑥
y y
a a
• The deflection of the plate 2 considering the
lateral load q on the plate can be written as
follows.
𝑚−1
2𝑞𝑏4 ∞ (−1) 2 𝑐𝑜𝑠β𝑚 𝑦 1
• 𝑤2 = 5
σ𝑚=1,3 [ β𝑚 𝑎
𝐷π 𝑚5 𝑐𝑜𝑠ℎ
2
β𝑚 𝑎 β𝑚 𝑎
• β𝑚 𝑥𝑠𝑖𝑛ℎβ𝑚 𝑥 − 2 + 𝑡𝑎𝑛ℎ 𝑐𝑜𝑠β𝑚 𝑥} + 2]
2 2
• The slope along the edge x=a/2 can now be
determined from the above expression as
follows.
𝑚−1
𝜕𝑤2 2𝑞𝑏3 (−1) 2 𝑐𝑜𝑠β𝑚 𝑦 β𝑚 𝑎 β𝑚 𝑎
• (at x=a/2)= σ∞
𝑚=1,3 [ 2 β𝑚 𝑎
− 𝑡𝑎𝑛ℎ ]
𝜕𝑥 𝐷π4 𝑚4 2𝑐𝑜𝑠ℎ 2
2
𝑎
• The variation of B.M 𝑀𝑥 along 𝑥 = ± can be
2
represented as follows.
𝑚−1
𝑎
• 𝑀𝑥 along 𝑥 = ± = σ∞
𝑚=1,3(−1) 2 𝐸𝑚 𝑐𝑜𝑠β𝑚 𝑦
2
• The deflection function for a plate subjected to edge
𝑎
moments along 𝑥 = ± can be written as
2
𝑚−1
𝑏2 ∞ (−1) 2 𝐸𝑚 𝑐𝑜𝑠β𝑚 𝑦 β𝑚 𝑎 β𝑚 𝑎
• 𝑤′2 = σ [ 𝑡𝑎𝑛ℎ
2π2 𝐷 𝑚=1,3,5 β 𝑎
𝑚2 𝑐𝑜𝑠ℎ 𝑚 2 2
2
• 𝑐𝑜𝑠ℎβ𝑚 𝑥 − β𝑚 𝑥sinhβ𝑚 𝑥]
• The slope along the edge x=a/2 due to edge moment
can be determined as follows.
𝑚−1
𝜕𝑤′2 𝑏 (−1) 2 𝐸𝑚 𝑐𝑜𝑠β𝑚 𝑦
• (along x=a/2)= − σ∞
𝜕𝑥 2π𝐷 𝑚=1,3 𝑚
β𝑚 𝑎 β𝑚 𝑎
• [𝑡𝑎𝑛ℎ + β 𝑎 ]
2 2𝑐𝑜𝑠ℎ2 𝑚
2
β 𝑎
•
𝑏 𝐸
= 4π𝐷 ( 𝑚𝑚 )[𝑡𝑎𝑛 𝑚
2
+ 𝑐𝑜𝑡ℎ β𝑚2 𝑎 + β𝑚 𝑎
2 β𝑚 𝑎 − β𝑚 𝑎
2 β𝑚 𝑎]
2𝑐𝑜𝑠ℎ 2 2𝑠𝑖𝑛ℎ 2
ℎ
C 3
y h
𝑎 3
• The side of the plate a. Then altitude h=
2
• Now applied load is in the form of edge moment
𝑀𝑛 on the edges.
• We can see that at the support,
• 𝑀𝑥𝑥 + 𝑀𝑦𝑦 = 𝑀𝑛 + 𝑀𝑡 (due to invariance
• property of curvature)
• Since 𝑀𝑡 = 0, we see that 𝑀𝑥𝑥 + 𝑀𝑦𝑦 = 𝑀𝑛 on the
boundary.
• Since applied distributed load q=0, we get the
equation 2 as 𝛻 2 𝑀 = 0
• This 2 is satisfied by taking 𝑀𝑥𝑥 + 𝑀𝑦𝑦 = 𝑀𝑛
throughout the plate.
𝑀𝑛
• Thus we get 𝑀 = throughout the plate.
(1+ν)
• Thus equation 1 is written as follows.
2 𝑀𝑛
• 𝛻 𝑤= −
(1+ν)𝐷
• Thus we have to find deflection w satisfying the
above equation.
• We know that w=0 on the edges of the plate.
• Thus we can write the equations of the edges.
𝑥 𝑥
• Edge OA-- Equation is 𝑦 = (Or) (𝑦 − ) =0
3 3
𝑥 𝑥
• Edge OB-- Equation is = − (Or) (𝑦 + ) =0
3 3
• Edge AB--- Equation is x = h (Or) (x − ℎ) = 0
𝑥 𝑥
• Let 𝑤 = 𝐶(𝑦 − ) (𝑦 + )(𝑥 − ℎ) where C is a
3 3
constant.
• The above function satisfies the support
conditions for w.(=0 on the edges)
2 𝑥2
• Hence w= 𝐶(𝑦 − )(𝑥 − ℎ)
3
2 2 𝑥3 𝑥2ℎ
• = 𝐶(𝑦 𝑥 − 𝑦 ℎ − + )
3 3
• Now substituting this in the equation
2 𝜕2 𝑤 𝜕2 𝑤 𝑀𝑛
• 𝛻 𝑤= + 2 = − , we get
𝜕𝑥 2 𝜕𝑦 (1+ν)𝐷
2ℎ 𝑀𝑛
• 𝐶(−2𝑥 + + 2𝑥 − 2ℎ) = −
3 (1+ν)𝐷
3𝑀𝑛
• Thus we get 𝐶 =
4(1+ν)𝐷ℎ
𝑥3 𝑥2ℎ
• Thus w = 𝐶(𝑦 2 𝑥 − 𝑦 2 ℎ − + )
3 3
3𝑀𝑛 𝑥3 𝑥2ℎ
• = (𝑦 2 𝑥 −𝑦 ℎ 2
− + )
4(1+ν)𝐷ℎ 3 3
• Now deflection on x axis(y=0) will be
3𝑀𝑛 𝑥3 𝑥2ℎ
• w= (− + )
4(1+ν)𝐷ℎ 3 3
𝑑𝑤
• Maximum occurs when =0
𝑑𝑥
2𝑥ℎ 2ℎ
• This gives −𝑥 2 + = 0 Or x(𝑥 − ) =0
3 3
𝑴𝒏 𝒉𝟐
• Taking x=2h/3, we get 𝑤𝑚𝑎𝑥 = (At centroid)
𝟐𝟕 𝟏+𝝂 𝑫
Analysis of a Clamped Elliptical
Plate
• Consider a clamped elliptical plate as shown in
the figure subjected to udl.
Clamped
a
x
y
• The semi major axis =a
• Semi minor axis =b
• The BCs here are
• Deflection, w=0 and normal slope =0
𝑥2 𝑦2
• Now let w(x,y)= + − 1)2 where C is a
𝐶( 2
𝑎 𝑏2
constant.
• We can see that both BCs will be satisfied
because equation of the boundary is
𝑥2 𝑦2
• + −1 =0
𝑎2 𝑏2
• Now the deflection should satisfy the plate
differential equation 𝐷𝛻 4 𝑤 = 𝑞.
𝑥4 𝑦4 𝑥2𝑦2 𝑥2 𝑦2
• Now w(x,y)= 𝐶( 4 + +1+ 2 2 2 −2 2 − 2 2)
𝑎 𝑏4 𝑎 𝑏 𝑎 𝑏
• Substituting in the differential equation, we get
24 24 16
• CD + + 2 2 =𝑞
𝑎4 𝑏 4 𝑎 𝑏
𝑞
• (Or), constant C= 24 24 16
+ + 𝐷
𝑎4 𝑏4 𝑎2 𝑏2
z y
• The plate is subjected to a concentrated force
P.
• Here we replace the triangular plate by an
equivalent square plate up on which two loads
+P and –P are acting.
• The point of application of the fictitious load –P
is the mirror image point of the actual load with
respect to the diagonal.
• We can see that the deflection of the original
and equivalent plate are the same.
• Hence the solution for the given problem is
given as below.
• Deflection of the plate, w =𝑤1 + 𝑤2
• The deflection for downward load P is given by
4𝑃
𝑠𝑖𝑛α𝑚 ξ 𝑠𝑖𝑛β𝑛 η
• 𝑤1 = σ∞
𝑚=1 σ
∞ 𝑎𝑏
𝑛=1 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
𝐷(α2𝑚 +β2𝑛 )2
• Here it is a square plate . Hence a=b
• Hence
4𝑃𝑎2 ∞ ∞ 𝑠𝑖𝑛α𝑚 ξ 𝑠𝑖𝑛β𝑛 η
• 𝑤1 = σ σ 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
π4 𝐷 𝑚=1 𝑛=1 𝐷(𝑚2 +𝑛2 )2
• Due to upward load P (image load whose location is
given by the coordinates (𝑎 − ξ, 𝑎 − η)), we have
• 𝑤2 =
4𝑃𝑎2 ∞ ∞ 𝑠𝑖𝑛α𝑚 (𝑎−ξ) 𝑠𝑖𝑛β𝑛 (𝑎−η)
− 4
σ𝑚=1 σ𝑛=1 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
π 𝐷 𝐷(𝑚2 +𝑛2 )2
• Hence the deflection w is given by
4𝑃𝑎2 ∞ ∞ 𝑠𝑖𝑛α𝑚 ξ 𝑠𝑖𝑛β𝑛 η
• 𝑤= 4
σ𝑚=1 σ𝑛=1 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
π 𝐷 𝐷(𝑚2 +𝑛2 )2
4𝑃𝑎2 ∞ ∞ 𝑠𝑖𝑛α𝑚 (𝑎−ξ) 𝑠𝑖𝑛β𝑛 (𝑎−η)
• − 4 σ𝑚=1 σ𝑛=1 𝑠𝑖𝑛α𝑚 𝑥 𝑠𝑖𝑛β𝑛 𝑦
π 𝐷 𝐷(𝑚2 +𝑛2 )2
• Similar method may also be adopted for triangular
plates subjected to uniformly distributed load.