Three-dimensional elasticity problems are difficult to solve.

Thus we first
develop governing equations for two-dimensional problems, and explore
four different theories:
- Plane Strain
- Plane Stress
- Generalized Plane Stress
- Anti-Plane Strain

Since all real elastic structures are three-dimensional, theories set forth
here will be approximate models. The nature and accuracy of the
approximation will depend on problem and loading geometry.
The basic theories of plane strain and plane stress represent the
fundamental plane problem in elasticity. While these two theories apply
to significantly different types of two-dimensional bodies, their formulations
yield very similar field equations.

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Two vs Three Dimensional Problems
Three-Dimensional Two-Dimensional

x x

y y

z z

z
Spherical Cavity

y
x

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Plane Strain
Consider an infinitely long cylindrical (prismatic) body as shown. If the body
forces and tractions on lateral boundaries are independent of the z-coordinate
and have no z-component, then the deformation field can be taken in the
reduced form
u u ( x, y ) , v v ( x, y ) , w 0

R
z

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Plane Strain Field Equations
u v 1 u v
Strains ex , ey , exy , ez exz e yz 0
x y 2 y x

x (ex e y ) 2ex , y (ex e y ) 2e y

Stresses z (e x e y ) ( x y )
xy 2exy , xz yz 0

Equilibrium Equations Strain Compatibility

x xy 2 ex e y
2
2 exy
Fx 0 2
x y y 2 x 2 xy
xy y
Fy 0
x y Beltrami-Michell Equation
Navier Equations 1 Fx Fy
2 ( x y )
2 u ( )
u v
Fx 0 1 x y
x x y
u v
2 v ( ) Fy 0
y x y

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Examples of Plane Strain Problems
y

x
x

Long Cylinders Semi-Infinite Regions

Under Uniform Loading Under Uniform Loadings

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Plane Stress
Consider the domain bounded two stress free planes z = h, where h is
small in comparison to other dimensions in the problem. Since the region is
thin in the z-direction, there can be little variation in the stress components
z , xz , yz through the thickness, and thus they will be approximately zero
throughout the entire domain. Finally since the region is thin in the z-
direction it can be argued that the other non-zero stresses will have little
variation with z. Under these assumptions, the stress field can be taken as
y

x x ( x, y ) 2h
y y ( x, y )
xy xy ( x, y )
z xz yz 0
R
z
x

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Plane Stress Field Equations
Strains Strain Displacement Relations
1 1
ex ( x y ) , e y ( y x ) u v w 1 u v
E E ex , ey , ez , exy
x y z 2 y x
e z ( x y ) (e x e y )
E 1 1 v w 1 u w
e yz 0 , exz 0
1 2 z y 2 z x
exy xy , exz e yz 0
E
Equilibrium Equations
Strain Compatibility
x xy
Fx 0 2 ex e y
2
2 exy
x y 2
xy y
y 2 x 2 xy
Fy 0
x y
Beltrami-Michell Equation
Navier Equations
F Fy
E u v 2 ( x y ) (1 ) x
2 u Fx 0 x y
2(1 ) x x y
E u v
2 v Fy 0
2(1 ) y x y

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Examples of Plane Stress Problems

Thin Plate With

Central Hole

Circular Plate Under

Edge Loadings
Elasticity Theory, Applications and Numerics
M.H. Sadd , University of Rhode Island
 Plane Elasticity Boundary Value Problem
Displacement Boundary Conditions
u ub ( x, y ) , v vb ( x, y ) on S Si
So

Stress/Traction Boundary Conditions

Txn Tx( b ) ( x, y ) (xb ) n x (xyb ) n y
y R
T T
y
n
y
(b)
( x, y ) n x n y
(b)
xy
(b)
y on S S = Si + So

Plane Strain Problem - Determine in- Plane Stress Problem - Determine in-
plane displacements, strains and stresses plane displacements, strains and stresses
{u, v, ex , ey , exy , x , y , xy} in R. Out-of- {u, v, ex , ey , exy , x , y , xy} in R. Out-of-
plane stress z can be determined from plane strain ez can be determined from
in-plane stresses via relation (7.1.3)3. in-plane strains via relation (7.2.2)3.

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Correspondence Between Plane Formulations
Plane Strain Plane Stress
u v E u v
2 u ( ) Fx 0 2 u Fx 0
x x y 2(1 ) x x y
u v E u v
2 v ( ) Fy 0 2 v Fy 0
y x y 2(1 ) y x y

x xy x xy
Fx 0 Fx 0
x y x y
xy y xy y
Fy 0 Fy 0
x y x y

1 Fx Fy F Fy
2 ( x y ) 2 ( x y ) (1 ) x
1 x y x y

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Transformation Between
Plane Strain and Plane Stress
Plane strain and plane stress field equations had identical equilibrium equations
and boundary conditions. Naviers equations and compatibility relations were
similar but not identical with differences occurring only in particular coefficients
involving just elastic constants. So perhaps a simple change in elastic moduli
would bring one set of relations into an exact match with the corresponding
result from the other plane theory. This in fact can be done using results in the
following table.
E
E
Plane Stress to Plane Strain
1 2 1
E (1 2)
Plane Strain to Plane Stress
(1 ) 2 1

Therefore the solution to one plane problem also yields the solution to the other
plane problem through this simple transformation scheme.

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Generalized Plane Stress
The plane stress formulation produced some inconsistencies in particular out-
of-plane behavior and resulted in some three-dimensional effects where in-
plane displacements were functions of z. We avoided these issues by simply
neglecting some of the troublesome equations thereby producing an
approximate elasticity formulation. In order to avoid this unpleasant situation,
an alternate approach called Generalized Plane Stress can be constructed
based on averaging the field quantities through the thickness of the domain.

1 h
Using the averaging operator defined by ( x, y )
2h h
( x, y , z )dz

all plane stress equations are satisfied exactly by the averaged stress, strain and
displacements variables; thereby eliminating the inconsistencies found in the
original plane stress formulation. However, this gain in rigor does not generally
contribute much to applications .

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Anti-Plane Strain
An additional plane theory of elasticity called Anti-Plane Strain involves a
formulation based on the existence of only out-of-plane deformation starting
with an assumed displacement field

u v 0 , w w( x, y )

Strains Stresses
e x e y ez e xy 0
x y z xy 0
1 w 1 w
e xz , e yz xz 2exz , yz 2e yz
2 x 2 y

Equilibrium Equations Naviers Equation

xz yz 2 w Fz 0
Fz 0
x y
Fx Fy 0

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Airy Stress Function Method
Numerous solutions to plane strain and plane stress problems can be determined
using an Airy Stress Function technique. The method reduces the general
formulation to a single governing equation in terms of a single unknown. The
resulting equation is then solvable by several methods of applied mathematics,
and thus many analytical solutions to problems of interest can be found.
This scheme is based on the general idea of developing a representation for the
stress field that will automatically satisfy equilibrium by using the relations
2 2 2
x 2 , y 2 , xy
y x xy
where = (x,y) is an arbitrary form called Airys stress function. It is easily
shown that this form satisfies equilibrium (zero body force case) and
substituting it into the compatibility equations gives
4 4 4
2 2 2 4 4 0
x 4
x y y
This relation is called the biharmonic equation and its solutions are known as
biharmonic functions.

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Airy Stress Function Formulation
The plane problem of elasticity can be reduced to a single equation in terms
of the Airy stress function. This function is to be determined in the two-
dimensional region R bounded by the boundary S as shown. Appropriate
boundary conditions over S are necessary to complete the solution.
Traction boundary conditions would involve the specification of second
derivatives of the stress function; however, this condition can be reduced to
specification of first order derivatives.

4 4 4
2 2 2 4 4 0 Si
x 4
x y y So

2 2
T(n)
x nx xy n y 2 nx ny
y xy
x

y R y
2 2 S = Si + So
T (n)
xy nx y n y nx 2 n y xy
xy x
y
x

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Polar Coordinate Formulation
Plane Elasticity Problem
Strain-Displacement
ur
er
r
1 u
e ur
r
1 1 ur u u
er
2 r r r

Hookes Law
Plane Strain Plane Stress
r ( er e ) 2er 1
er ( r )
( er e ) 2e E
1
z ( e r e ) ( r ) e ( r )
E
r 2er , z rz 0

ez ( r ) ( e r e )
E 1
1
er r , ez erz 0
E

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Polar Coordinate Formulation
Naviers Equations Equilibrium Equations
Plane Strain r 1 r ( r )
Fr 0
ur ur 1 u r r r
2 u r ( ) Fr 0 r 1 2 r
r r r r F 0
r r r
1 ur ur 1 u
2 u ( ) F 0
r r r r
Compatibility Equations
Plane Stress
Plane Strain
E ur ur 1 u
2 u r Fr 0 1 Fr Fr 1 F
2(1 ) r r r r 2 ( r )
1 r r r
E 1 u r ur 1 u
2 u F 0 Plane Stress
2(1 ) r r r r
F F 1 F
2 ( r ) (1 ) r r
r r r
2 1 1 2
2
2

r r r r 2 2

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island
 Polar Coordinate Formulation
Airy Stress Function Approach = (r,)
Airy Representation
1 1 2
r Biharmonic Governing Equation
r r r 2 2
2 2 1 1 2 2 1 1 2
2 2
4
2 2 2 2 2 0
r r r r r r r r r
1
r
r r

r S
R Traction Boundary Conditions
y
r
Tr f r ( r, ) , T f ( r, )


r

x

Elasticity Theory, Applications and Numerics

M.H. Sadd , University of Rhode Island