Plates and Shells
All images are from R. Cook, et al. Concepts and Applications of Finite Element Analysis, 1996.
MAE456 Finite Element Analysis
Plate Formulation
• Plates may be considered similar to
beams, however:
– Plates can bend in two directions
– Plates are flat with a thickness (can’t have an
interesting cross-section)
MAE456 Finite Element Analysis 2
1
Thin Plate Formulation
• Consider a thin plate on the xy plane (z = 0),
with thickness t, & neglecting shear strain.
• If we take a differential slice from plate:
MAE456 Finite Element Analysis 3
Thin Plate Formulation
then:
w = w( x , y )
∂w
u = −z
∂x
∂w γyz = γzx = 0
v = −z
∂y
• Assume σz = 0. Therefore:
MAE456 Finite Element Analysis 4
2
Thin Plate Formulation
• These stresses give rise to moments:
• Maximum stresses are therefore given by:
6M x 2z
σ x ,max = 2 since σ x = σ x ,max ,
t t
6M y
σ y ,max = 2 ,
t
6 M xy
τ xy ,max = 2
t
MAE456 Finite Element Analysis 5
Thin Plate Formulation
• This is similar to the beam formula, but
since the plate is very wide we have a
situation similar to plain strain.
• For a unit width beam, flexural rigidity
D=EI=Et3/12.
• For a unit width plate, flexural rigidity
D=EI/(1-ν 2)=Et3/[12(1-ν 2)].
• This thin plate theory is also called the
“Kirchhoff plate theory.”
MAE456 Finite Element Analysis 6
3
Mindlin Plate Theory
• Mindlin Plate Theory assumes that
transverse shear deformation also occurs.
MAE456 Finite Element Analysis 7
Mindlin Plate Theory
• The deformations and strains are therefore
given by:
MAE456 Finite Element Analysis 8
4
Mindlin Plate Theory
• Mindlin plate elements are
more common than
Kirchhoff elements.
• The displacement interpolation is given by:
Ni can be the same shape functions as for Q4 and Q8
quadrilateral elements.
MAE456 Finite Element Analysis 9
Support Conditions
• Support Conditions are similar to those for
beams:
θn, Mn – rotation and moment normal to edge
θs, Ms– rotation and moment perpendicular to edge
For Mindlin plates, do not restrain θn, to avoid
accuracy problems.
MAE456 Finite Element Analysis 10
5
Test Cases
• For plate elements, patch tests and single
element tests should include the cases
shown:
MAE456 Finite Element Analysis 11
Test Cases
• Plate elements must be able to show
constant σx, σy and τxy at each z level to
pass a patch test. They must pass the test
for constant Mx, My and Mxy.
• Many element formulations perform poorly
for these tests.
MAE456 Finite Element Analysis 12
6
Large Displacements and
Membrane Forces
• A beam with fixed supports will exhibit “string
action” axial forces as shown.
• If we consider both string action and bending
stresses, a beam can carry a distributed load of:
MAE456 Finite Element Analysis 13
Large Displacements and
Membrane Forces
• A similar situation arises with plates,
however basic plate elements are not set up
to handle “membrane” forces.
• If w/t is large (e.g., greater than 0.1), a non-
linear analysis must be performed using
shell elements, which do handle membrane
forces.
• In general, tensile membrane forces will
have a stiffening effect and compressive
membrane forces will decrease stiffness.
MAE456 Finite Element Analysis 14
7
Shell Finite Elements
• Shell elements are different from plate
elements in that:
– They carry membrane AND bending forces
– They can be curved
• The most simple shell element combines a
bending element with a membrane element.
– E.g., combines a plate element and a plane stress
element.
– These elements are flat, therefore it is important that
elements are not all coplanar where they meet at a
node.
MAE456 Finite Element Analysis 15
Shell Finite Elements
• Curved shell elements can be derived
using “shell theory.”
• “Isoparametric” shell elements can also be
obtained by starting with a solid element
and reducing degrees of freedom.
• Thin shell behavior varies widely between
formulations and should be tested before
use.
MAE456 Finite Element Analysis 16
8
Shells and Shell Theory
• A thin shell structure can carry high loads if
membrane stresses predominate.
• However, localized bending stresses will
appear near load concentrations or
geometric discontinuities.
MAE456 Finite Element Analysis 17
Shells and Shell Theory
• Localized bending stresses appear in many
different situations:
MAE456 Finite Element Analysis 18
9
Shells and Shell Theory
• A thin-walled cylindrical tank has high
bending (flexural) stresses at the base.
• Use a finer mesh where there are discontinuities
or abrupt changes in the structure.
MAE456 Finite Element Analysis
Shells and Shell Theory
• For a cylindrical shell of radius R and
thickness t, the localized bending dies out
after a distance λ:
• Membrane stresses do not die out.
MAE456 Finite Element Analysis 20
10
Using Shell Elements to Model Beams
• To do a proper FE analysis, the analyst must
understand how the structure is likely to
behave and how elements are able to
behave.
• In some cases it is more appropriate use
shell elements rather than beam elements.
MAE456 Finite Element Analysis 21
Using Shell Elements to Model Beams
• A curved I-beam reacts to moments as shown,
therefore shell elements would be more accurate
than beam elements.
• Pipe bends react to moments as shown. Use
shell elements or specialized beam elements
with correction factors.
MAE456 Finite Element Analysis
11
Using Shell Elements to Model Beams
• If the load is not applied directly below the
“shear center”, the channel will twist. Use
shell elements instead of beam elements.
MAE456 Finite Element Analysis
Using Shell Elements to Model Beams
• If beam flanges are wide, σx = My/I is not
accurate. Beam elements will not give
accurate results.
• In this case, plate/shell elements should be
used.
MAE456 Finite Element Analysis
12
