IIT Delhi: 1 Conditions For Applying Mohr's Circle (Start Time: 00:26)
IIT Delhi: 1 Conditions For Applying Mohr's Circle (Start Time: 00:26)
Lecture - 9
Mohr’s Circle
Abstract
In the last two lectures, we have learnt about planes on which normal and shear components
of traction are maximized/minimized. In this lecture, we will discuss about Mohr’s circle, a
graphical technique to find components of traction on arbitrary planes at a point.
i
lh
De
1 Conditions for applying Mohr’s Circle (start time: 00:26)
IIT
Mohr’s circle is a graphical way to find normal and shear components of traction on arbitrary planes.
The only restriction is that the plane normal has to be perpendicular to one of the principal stress
,
ar
directions. Accordingly, let us consider a coordinate system such that the third coordinate axis
um
is along one of the principal stress directions (the rest two coordinate axes need not be along any
principal direction, also see Figure 1) and we are looking at planes whose normals are perpendicular
tK
to the third principal direction. Then, all such plane normals will have n3 = 0. On such planes,
ee
we want to find the normal and shear components of traction (see Figure 2). The stress matrix
Aj
in this coordinate system will be such that its column will be formed by traction on the plane
along third coordinate axis. But that being a principal plane, the third column will not have any
f.
shear component. Hence, the third column will have its last entry nonzero but the other two zero.
o
Pr
Symmetry of stress matrix will further force the first two entries of last row to be zero, i.e.,
@
× × 0
ht
σ = × × 0 . (1)
rig
0 0 ×
py
Co
Also, the normal vectors of the plane on which we want to obtain normal and shear components
will look like
×
n = × . (2)
0
1
i
lh
De
IIT
(principal axis)
,
ar
um
Figure 1: The cuboid with face normals along the three coordinate axes. The third coordinate axis
tK
have normal component of traction denoted as σxx and shear component of traction denoted as τyx
pointing towards y axis. The third component (τzx ) is absent as per the stress representation (1).
o f.
Similarly, on the e2 plane, we have τxy (same as τyx ) and on the e3 plane, we have σzz (equal to λ3 )
Pr
only. On the −e1 plane, the normal component is σxx in −x direction and shear component is τyx
@
There is another simpler way to draw such a state of stress for whom there are no shear components
in the third direction: we can just draw a square instead of a cube where the sides of the square
denote the faces of the cuboid as shown in Figure 2. The right edge of the square represents the
(+e1 ) face. Similarly, the top edge represents +e2 face and the plane of the square itself represents
e3 face. So, the plane containing the square has just one traction component σzz . We use a dot
enclosed by a circle to denote the traction component coming out of the plane as shown in Figure
2. On the edges of the square, we have both the shear and the normal components of traction. On
the right edge, representing e1 face, we have σxx and τyx (which is equal to τxy in magnitude). So,
we will just write τxy to denote both τyx and τxy . On the top edge, representing e2 face, we have
τxy and σyy . Keep in mind that such a reduce to square to represent stress matrix is possible only
when the third coordinate axis lies along one of the principal directions.
2
i
lh
De
, IIT
ar
um
Figure 2: A square with its sides representing the faces of the cuboid. Traction components are
tK
also drawn.
ee
Aj
Now, our goal is to calculate the normal and shear components of traction on planes whose normals
o
are perpendicular to e3 . The blue line in Figure 2 shows a general plane of such kind. The normal
Pr
to this plane is represented by n and assume that it makes an angle α with e1 axis. To get to this
@
arbitrary plane n, we can rotate our e1 plane by α about the e3 axis. The column representation of
ht
cos(α)
Co
n = sin(α) . (3)
0
An important point to note here is that we have assumed the direction of τ on this plane makes
90◦ (anti clockwise) from the normal vector n (also see Figure 2). Now, σ will be given by
σxx τxy 0 cos(α) ! cos(α)
σ= σ n · n = τxy σyy 0 sin(α) · sin(α)
0 0 σzz 0 0
To get τ , we need to first represent the direction along which it is acting which we denote by n⊥ . If
we look at Figure 2, we can find the angle that n⊥ makes with all the three axes. The representation
3
of n⊥ will thus be
−sin(α)
n⊥ = cos(α) . (5)
0
Now, to get τ , we need to take the component of total traction ( σ n ) on this plane along the
n⊥ direction, i.e.,
σxx τxy 0 cos(α) ! −sin(α)
τ = σ n · n⊥ = τxy σyy 0 sin(α) · cos(α)
0 0 σzz 0 0
= −σxx cos(α)sin(α) − τxy sin2 (α) + τxy cos2 (α) + σyy cos(α)sin(α) (6)
i
Upon doing some algebraic manipulation in equation (4), we get
lh
De
σxx + σyy 2 1 2 1
σ= + σxx cos (α) − + σyy sin (α) − + 2τxy sin(α)cos(α). (7)
IIT
2 2 2
Further, using the following trigonometric identities:
,
ar
cos(2α) = cos2 (α) − sin2 (α) = 2cos2 (α) − 1 = 1 − 2sin2 (α),
um (8)
sin(2α) = 2sin(α)cos(α), (9)
tK
we obtain
ee
Aj
Similar rearrangements and simplifications using trigonometric identities (8) and (9) in equation
rig
(6) gives us
σxx − σyy
py
4
σxx − σyy
Figure 3: A right angled triangle with the two arms as and τxy .
2
i
lh
De
Using equation (13) in equations (10) and (11), we get
IIT
σxx + σyy
σ = + R cos(2φ)cos(2α) + sin(2φ)sin(2α)
h 2
,
i
ar
τ = R − cos(2φ)sin(2α) + sin(2φ)sin(2α) (14)
um
Upon further using following trigonometric identities:
tK
we get
o
σxx + σyy
Pr
σ = + Rcos(2φ − 2α)
2
@
(17)
ht
τ = Rsin(2φ − 2α)
rig
These are the formulae obtained for getting σ and τ on a plane making an angle α with e1 axis. For
py
a given stress matrix, we can find out R and φ using equations (12) and (13) respectively. Then,
Co
using equation (17), we can find σ and τ on the plane which is obtained by rotating e1 by angle α
about e3 axis. Let us try to represent the above two formulas graphically.
5
i
lh
De
IIT
Figure 4: σ − τ plane with various parameters plotted on it.
,
ar
um
plot the point (σxx , τxy ) which corresponds to e1 plane. If we join this point with the center, the
tK
line obtained will give us the radius of the circle. This is because the circle has to pass through
the point corresponding to e1 plane: the circle is the locus of all (σ, τ ) when α is varied and the
ee
point corresponding to e1 plane is obtained for α = 0. We can also verify that the radius R that
Aj
we obtained graphically also matches with the formula for R in equation (12).1 .
o f.
Once we have obtained the radius and the center of the circle, we can draw the complete circle
ht
as shown in Figure 5. This circle is called the Mohr’s circle. When we compare the right angled
rig
triangle in Figure 5 with Figure 3, we see that the angle that the line from center to the point
py
corresponding to e1 plane makes with the σ axis would be 2φ. To find σ and τ on any arbitrary
plane (for general α), let us look at equation (17): the angle in the cosine and sine terms there is
Co
(2φ − 2α). So, for a general plane, the argument in the trigonometric function there reduces by
2α. So, the radial line from the center to the point corresponding to the α-plane on Mohr’s circle
should be at an angle of (2φ − 2α) from the x-axis. Thus, we can obtain the point corresponding
to α-plane on Mohr’s circle by going in the clockwise direction by angle 2α from the e1 -plane point
as shown in Figure 5.
We summarize below the steps involved in drawing the Mohr’s circle and finding the point core-
sponding to α-plane:
σxx + σyy
1. Draw the center of the circle at ,0 .
2
1
In the right angled triangle shown in Figure 4, if we apply the Pythagoras theroem to find the hypotenuse, i.e.,
R, we get the same expression as the algebraic expression for R in (12)
6
i
lh
De
,IIT
ar
um
Figure 5: Mohr’s circle plot
tK
ee
3. the line joining the center and the point for e1 plane forms the radius of the circle.
o f.
Pr
5. To find (σ, τ ) for α-plane, rotate the radial line of e1 plane by 2α clockwise.
ht
rig
Notice that the normal to the required plane made an angle α with e1 in the counter clockwise
py
direction (also see Figure 2). But on the Mohr’s circle, we draw that point by rotating by 2α in the
clockwise direction from the point corresponding to the e1 plane. This is because in equation (17),
Co
we have 2α with a minus sign in the trigonometric functions. So, clockwise rotation of the plane
corresponds to counter-clockwise rotation in the Mohr’s circle and vice versa.
7
i
lh
De
, IIT
ar
um
tK
Mohr’s circle? This is because of our convention for the shear component of traction. In Figure 2,
f.
we had defined positive τ when we go 90◦ in the counter clockwise direction from n. So, to get shear
o
Pr
component on e2 plane, we go 90◦ in the anti-clockwise direction from e2 direction and thus, get
@
to the −e1 direction. So, Mohr’s circle is giving us τ on e2 plane in the −e1 direction whereas τxy ,
ht
by definition, is the shear component in the +e1 direction.Therefore, Mohr’s circle gives us −τxy as
rig
the shear traction on e2 plane. We can remember that this negative sign comes because n⊥ for e2
plane is along −e1 direction by convention.
py
Co
8
λ1 and λ2 will be obtained by adding and subtracting R to the σ for center respectively, i.e.,
σxx + σyy
λ1 = +R (18)
2
σxx + σyy
λ2 = −R (19)
2
This allows us to get the value of principal stress components directly from the Mohr’s circle:
otherwise one has to solve an eigenvalue problem. Writing the center of the circle in terms of
principal stress components, we obtain
λ1 + λ2
Center= ,0 (20)
2
i
lh
We can also get the radius of the circle in terms of principal stress components by subtracting
De
equation (19) from equation (18), i.e.,
λ1 − λ2
R= (21)
IIT
2
Notice that in the previous lecture, we had derived the maximum value of shear component of
,
λ1 − λ2 λ1 − λ2
ar
traction to be . And from Mohr’s circle too, we get τmax = R =
um . From the Mohr’s
2 2
cicrle, we can also see that σ corresponding to the point where we have maximum shear will be the
λ1 + λ2
tK
same as σ for the center of the circle, i.e. . This is the same value that we had derived in
2
ee
the previous lecture. We have thus verified our results of the Mohr’s circle.
Aj
o f.
Pr
@
ht
rig
py
Co