Section 1.

1.3 Compatibility of Strain

As seen in the previous section, the displacements can be determined from the strains
through integration, to within a rigid body motion. In the two-dimensional case, there are
three strain-displacement relations but only two displacement components. This implies
that the strains are not independent but are related in some way. The relations between
the strains are called compatibility conditions.

1.3.1 The Compatibility Relations

Differentiating the first of 1.2.5 twice with respect to y , the second twice with respect to
x and the third once each with respect to x and y yields

∂ 2 ε xx ∂ 3u x ∂ 2 ε yy ∂ 3u y ∂ 2 ε xy
1 ⎛⎜ ∂ 3 u x ∂ 3u y ⎞
= , = , = + ⎟
∂y 2 ∂x∂y 2 ∂x 2 ∂x 2 ∂y ∂x∂y 2 ⎜⎝ ∂x∂y 2 ∂x 2 ∂y ⎟⎠

It follows that

∂ 2 ε xx ∂ ε yy ∂ 2 ε xy
2

+ = 2 2-D Compatibility Equation (1.3.1)

∂y 2 ∂x 2 ∂x∂y

This compatibility condition is an equation which must be satisfied by the strains at all
material particles.

Physical Meaning of the Compatibility Condition

When all material particles in a component deform, translate and rotate, they need to meet
up again very much like the pieces of a jigsaw puzzle must fit together. Fig. 1.3.1
illustrates possible deformations and rigid body motions for three line elements in a
material. Compatibility ensures that they stay together after the deformation.

deformed
- compatibility ensured

deformed
undeformed - compatibility not satisfied

Figure 1.3.1: Deformation and Compatibility

Solid Mechanics Part II 17 Kelly

 Section 1.3

The Three Dimensional Case

There are six compatibility relations to be satisfied in the three dimensional case :

∂ 2 ε yy ∂ 2 ε zz ∂ 2 ε yz ∂ 2 ε xx ∂ ⎛ ∂ε yz ∂ε zx ∂ε xy ⎞
+ = 2 , = ⎜⎜ − + + ⎟
∂z 2 ∂y 2
∂y∂z ∂y∂z ∂x ⎝ ∂x ∂y ∂z ⎟⎠
∂ 2 ε zz ∂ 2 ε xx ∂ 2 ε zx ∂ 2 ε yy ∂ ⎛ ∂ε yz ∂ε zx ∂ε xy ⎞
+ = 2 , = ⎜⎜ + − + ⎟ (1.3.2)
∂x 2 ∂z 2
∂z∂x ∂z∂x ∂y ⎝ ∂x ∂y ∂z ⎟⎠
∂ 2 ε xx ∂ 2 ε yy ∂ 2 ε xy ∂ 2 ε zz ∂ ⎛ ∂ε yz ∂ε zx ∂ε xy ⎞
+ = 2 , = ⎜+ + − ⎟
∂y 2 ∂x 2 ∂x∂y ∂x∂y ∂z ⎜⎝ ∂x ∂y ∂z ⎟⎠

By inspection, it will be seen that these are satisfied by Eqns. 1.2.19.

1.3.2 Problems
1. The displacement field in a material is given by
u x = Axy, u y = Ay 2 ,
where A is a small constant. Determine
(a) the components of small strain
(b) the rotation
(c) the principal strains
(d) whether the compatibility condition is satisfied

Solid Mechanics Part II 18 Kelly

 Section 3.2

3.2 The Stress Function Method

An effective way of dealing with many two dimensional problems is to introduce a new
“unknown”, the Airy stress function  , an idea brought to us by George Airy in 1862.
The stresses are written in terms of this new function and a new differential equation is
obtained, one which can be solved more easily than Navier’s equations.

3.2.1 The Airy Stress Function

The stress components are written in the form

 2
 xx 
y 2
 2
 yy  (3.2.1)
x 2
 2
 xy 
xy

Note that, unlike stress and displacement, the Airy stress function has no obvious physical
meaning.

The reason for writing the stresses in the form 3.2.1 is that, provided the body forces are
zero, the equilibrium equations are automatically satisfied, which can be seen by
substituting Eqns. 3.2.1 into Eqns. 2.2.3 {▲Problem 1}. On this point, the body forces,
for example gravitational forces, are generally very small compared to the effect of
typical surface forces in elastic materials and may be safely ignored (see Problem 2 of
§2.1). When body forces are significant, Eqns. 3.2.1 can be amended and a solution
obtained using the Airy stress function, but this approach will not be followed here. A
number of examples including non-zero body forces are examined later on, using a
different solution method.

3.2.2 The Biharmonic Equation

The Compatability Condition and Stress-Strain Law

In the previous section, it was shown how one needs to solve the equilibrium equations,
the stress-strain constitutive law, and the strain-displacement relations, resulting in the
differential equation for displacements, Eqn. 3.1.4. An alternative approach is to ignore
the displacements and attempt to solve for the stresses and strains only. In other words,
the strain-displacement equations 3.1.2 are ignored. However, if one is solving for the
strains but not the displacements, one must ensure that the compatibility equation 1.3.1 is
satisfied.

Eqns. 3.2.1 already ensures that the equilibrium equations are satisfied, so combine now
the two dimensional compatibility relation and the stress-strain relations 3.1.1 to get
{▲Problem 2}

Solid Mechanics Part II 46 Kelly

 Section 3.2

  4  4  4 
plane stress :  4  2 2 2  4   0
 x x y y 
(3.2.2)
  4  4  4 
plane strain :  4  2 2 2  4 1     0
 x x y y 

Thus one has what is known as the biharmonic equation:

 4  4  4
 2 2 2  4  0 The biharmonic equation (3.2.3)
x 4 x y y

The biharmonic equation is often written using the short-hand notation  4  0 .

By using the Airy stress function representation, the problem of determining the stresses
in an elastic body is reduced to that of finding a solution to the biharmonic partial
differential equation 3.2.3 whose derivatives satisfy certain boundary conditions.

Note that the biharmonic equation is independent of elastic constants, Young’s modulus E
and Poisson’s ratio . Thus for bodies in a state of plane stress or plane strain, the stress
field is independent of the material properties, provided the boundary conditions are
expressed in terms of tractions (stress) 1; boundary conditions on displacement will bring
the elastic constants in through the stress-strain law. Further, the plane stress and plane
strain stress fields are identical.

3.2.3 Some Simple Solutions

Clearly, any polynomial of degree 3 or less will satisfy the biharmonic equation. Here
follow some elementary examples.

(i)   Ay 2
 2
one has  xx   2 A,  yy   xy  0 , a state if uniaxial tension
y 2

(ii)   Bxy
here,  xx   yy  0,  xy   B , a state of pure shear

(iii)   Ay 2  Bxy
here,  xx  2 A,  yy  0,  xy   B , a superposition of (i) and (ii)

1
technically speaking, this is true only in simply connected bodies, i.e. ones without any “holes”, since
problems involving bodies with holes have an implied displacement condition (see, for example, Barber
(1992), §2.2).

Solid Mechanics Part II 47 Kelly

 Section 3.2

3.2.4 Pure Bending of a Beam

Consider the bending of a rectangular beam by a moment M 0 , as shown in Fig. 3.2.1.

The elementary beam theory predicts that the stress  xx varies linearly with y, Fig. 3.2.1,
with the y  0 axis along the beam-centre, so a good place to start would be to choose, or
guess, as a stress function   Cy 3 , where C is some constant to be determined. Then
 xx  6Cy,  yy  0,  xy  0 , and the boundary conditions along the top and bottom of
the beam are clearly satisfied.

b
M0 M0

Figure 3.2.1: a beam in pure bending

The moment and stress distribution are related through

b b
M 0    xx ydy  6C  y 2 dy  4Cb 3 (3.2.4)
b b

and so C  M 0 / 4b 3 and  xx  3M 0 y / 2b 3 . The fact that this last expression agrees with
the elementary beam theory (    My / I with I  2b 3 h / 3 , where h is the depth “into
the page”) shows that that the beam theory is exact in this simple loading case.

Assume now plane strain conditions. In that case, there is another non-zero stress
component, acting “perpendicular to the page”,  zz   ( xx   yy )  3My / 2b 3 . Using
Eqns. 3.1.1b,

1   1   2 3M 
 xx   
(1   ) xx   yy  
2
 y  y, say
3 
E  E b  (3.2.5)
1   v(1   ) 3M 
 yy   
  xx  (1   ) yy   
E 2b 3 
 y   y, say
E 

and the other four strains are zero.

As in §1.2.4, once the strains have been found, the displacements can be found by
integrating the strain-displacement relations. Thus

Solid Mechanics Part II 48 Kelly

 Section 3.2

u x
 xx   y
x
 u x  xy  f ( y )
u y
 yy   y
y (3.2.6)
 u y   y  g ( x)
1
2
2

1  u u y  1
 xy   x    x  f ( y )  g ( x)   0
2  y x  2
 g ( x)  x   f ( y )

Therefore f ( y ) must be some constant,  C say, so f ( y )  Cy  A , and

g ( x)  Cx  12 x 2  B . Finally,

u x  xy  Cy  A
(3.2.7)
u y   12 x 2  12 y 2  Cx  B

which are of the form 1.2.17. For the case when the mid-point of the beam is fixed, so
has no translation, u x (0,0)  u y (0,0)  0 , and if it has no rotation there,  z (0,0)  0 , then
the three arbitrary constants are zero, and

u x  xy
(3.2.8)
u y   12 x 2  12  y 2

3.2.5 A Cantilevered Beam

Consider now the cantilevered beam shown in Fig. 3.2.2. The beam is subjected to a
uniform shear stress  xy   over its free end, Fig. 3.2.2a. The boundary conditions are

 xx (0, y )  0,  xy (0, y )   ,  yy ( x,b)   xy ( x,b)  0 (3.2.9)

It is difficult, if not impossible, to obtain concise expressions for stress and strain for
problems even as simple as this2. However, a concise solution can be obtained by
relaxing one of the above conditions. To this end, consider the similar problem of Fig.
3.2.2b – this beam is subjected to a shear force F, the resultant of the shear stresses. The
applied force of Fig. 3.2.2b is equivalent to that in Fig. 3.2.2a if

b


b
xy (0, y )dy  F (3.2.10)

2
an exact solution will usually require an infinite series of terms for the stress and strain

Solid Mechanics Part II 49 Kelly

 Section 3.2

This is known as a weak boundary condition, since the stress is not specified in a point-
wise sense along the boundary – only the resultant is. However, from Saint-Venant’s
principle (Part I, §3.3.2), the stress field in both beams will be the same except for in a
region close to the applied load.

y y

 b x x
F

(a ) ( b)

Figure 3.2.2: A cantilevered beam subjected to; (a) a uniform distribution of shear
stresses along its free end, (b) a shear force along its free end

The elementary beam theory predicts a stress  xx   My / I  Fxy / I . Thus a good place
to start is to choose the stress function   xy 3 , where  is a constant to be determined.
The stresses are then

 xx  6xy,  yy  0,  xy  3y 2 (3.2.11)

However, it can be seen that  xy ( x,b)  3b 2  0 . To offset this, one can superimpose
a constant shear stress 3b 2 , in other words amend the stress function to

  xy 3  3b 2 xy (3.2.12)

The boundary conditions are now satisfied and, from Eqn. 3.2.10,

F
 (3.2.13)
4b 3

and so

 xx 
3F
3
3F

xy,  yy  0,  xy  3 b 2  y 2  (3.2.14)
2b 4b

3.2.6 Problems

1. Verify that the relations 3.2.1 satisfy the equilibrium equations 2.2.3.

2. Derive Eqn. 3.2.2.

3. A large thin plate is subjected to certain boundary conditions on its thin edges (with
its large faces free of stress), leading to the stress function

Solid Mechanics Part II 50 Kelly

 Section 3.2

  Ax 3 y 2  Bx 5
(i) use the biharmonic equation to express A in terms of B
(ii) calculate all stress components
(iii) calculate all strain components (in terms of B, E,  )
(iv) derive an expression for the volumetric strain, in terms of B, E,  , x and y.
(v) check that the compatibility equation is satisfied
(vi) check that the equilibrium equations are satisfied

4. A very thick component has the same boundary conditions on any given cross-
section, leading to the following stress function:
  x 4 y  4x 2 y 3  y 5
(i) is this a valid stress function, i.e. does it satisfy the biharmonic equation?
(ii) calculate all stress components (with   1 / 4 )
(iii) calculate all strain components
(iv) find the displacements
(v) specify any three displacement components which will render the arbitrary
constant displacements of (iv) zero

5. For the cantilevered beam discussed in §3.2.5, evaluate the resultant shear force and
moment on an arbitrary cross-section x  x . Are they as you expect? (You will find
that the beam is in equilibrium, as expected, since the equilibrium equations have
been satisfied.)

6. For the cantilevered beam discussed in §3.2.5, evaluate the strains and displacements,
assuming plane stress conditions.
Note: to evaluate the three arbitrary constants of integration, one would be tempted to
apply the obvious u x  u y  0 all along the built-in end. However, since only weak
boundary conditions were imposed, one cannot enforce these strong conditions (try
it). Instead, apply the following weaker conditions: (i) the displacement at the built-in
end at y  0 is zero ( u x  u y  0 ), (ii) the slope there, u y / x , is zero.

7. Show that the stress function


p
20h 3
 
 20 y 3 L2  x 2   4 y 5  15h 2 x 2 y  2h 2 y 3  5h 3 x 2
satisfies the boundary conditions for the simply supported beam subjected to a
uniform pressure p shown below. Check the boundary conditions in the weak (Saint-
Venant) sense on the shorter left and right hand sides (for both normal and shear
stress). Since the normal stress  xx is not zero at the ends, but only its resultant,
check also that the moment is zero at each end.

y p

Lp x Lp h

L L

Solid Mechanics Part II 51 Kelly

 Section 3.2

Note that the elementary beam theory predicts an approximate flexural stress but an
exact shear stress:
6 p  h2 
 xx   3 y L2  x 2 ,  xy  3 x  y 2 
6p
h h  4 

8. Consider the dam shown in the figure below. Assume first a general cubic stress
function
1 1 1 1
  C1 x 3  C 2 x 2 y  C3 xy 2  C 4 y 3
6 2 2 6
Apply the boundary conditions to determine the constants and hence the stresses in
the dam, in terms of  , the density of water. (Use the stress transformation equations
for the sloped boundary and ignore the weight of the dam.)
[Just consider the effect of the water; to these must be added the stresses resulting
from the weight of the dam itself, which are given by
 1 
 xx  0,  yy   s g  x  y ,  xy  0
 tan  
where  s is the density of the dam material.]




Solid Mechanics Part II 52 Kelly

Alexandria University ‫جامعة االسكندرية‬

Faculty of Engineering ‫كلية الهندسة‬

Mechanical Engineering Department ‫قسم هندسة القوى الميكانيكية‬

Aug. 2020 2020 ‫أغسطس‬

THEORY OF ELASTICITY (MEC 0708747)

Take Home Exam

Answer All Questions:

Question One:

a) In the following phrase, add the missing words:

“At the mating surface between a disc shrunk on a shaft, the radial stress is ………… and
when rotation starts there is a superposition of ……………stress. It is therefore necessary to
ensure that the shrink-fit stress is always ………. the ……….. stress at the mating surface so
that the disc does not become …….. on the shaft.

b) Discuss briefly the difference between the plane stress case and the plane strain case.

d) Discuss the difference between isotropic and anisotropic materials.

Question Two:

The following stress distribution has been determined for a machine component:

σxx = 3x3 – 3y2 – z , σyy = 3y2 , σzz = 3x + y – z + 5/4

τxy = z- 6xy – 3/4 , τyx = 0 , τzx = x + y – 3/2

Determine if equilibrium is satisfied in the absence of body forces.

____________________________________________________________________________-

Question Three:

For a prismatical bar stretched by its own weight, discuss the following knowing that the bar
area is A and its length is L:

a- The boundary conditions of the bar.

b- The compatibility equations.
c- The displacements.
d- The vertical displacement of the lower end of the bar.
Question Four:

a- Explain what is meant by speed limit of a rotating disk shrunk onto a shaft.
b- A steel bush is to be shrunk onto a steel shaft so that the internal diameter
is extended by 0.152 mm above its original size. The inside and outside
diameters of the bush are 203 mm and 305 mm respectively. Find:
i- The normal pressure intensity between the bush and the shaft.
ii- The hoop stress at the inner and outer surfaces of the bush.

Take E = 207 GN/ m² & ʋ = 0.28

Question Five:

Determine whether the following strain fields are compatible:

Field 1 Field 2

ϵxx = 2x2 + 3y2 + z + 1 ϵxx = 3y2 + xy

ϵyy = 2y2 + x2 + 3z + 2 ϵyy = 2y + 4z + 3

ϵzz = 3x + 2y + z2 + 1 ϵzz = 3zx + 2xy + 3yz + 2

ϒxy = 8xy ϒxy = 6xy

ϒyz = 0 ϒyz = 2x

ϒzx = 0 ϒzx = 2y

Good Luck

Dr. Aly Helmy

 Alexandria University
Faculty of Engineering
Mechanical Dept.
Theory of Elasticity

Theory of Elasticity
Report

Student name: Karim Wagdy Mohamed Mahmoud

Submitted to: Prof.Dr / Aly Helmy

c)

f)
Question 5
Question 6

Question 7
Question 8
Question 9

Question 10
Question 11