Stress Analysis, Strain Analysis, and Shearing of Soils

saL00581_ch04.
qxd 9/14/06 4:54 PM Page 113
C H A P T E R 4
Stress Analysis, Strain
Analysis, and Shearing
of Soils
Ut tensio sic vis (strains and stresses are related linearly).

Robert Hooke
So I think we really have to, first, make some new kind of theories in which we take
regard to the fact that there is no linearity condition between stresses and strains
for soils.
J. Brinch Hansen
Foundation design requires that we analyze how structural loads are transferred to
the ground, and whether the soil will be able to support these loads safely and
without excessive deformation. In order to do that, we treat soil deposits as contin-
uous masses and develop the analyses needed for determining the stresses and
strains that appear in the soil mass as a result of application of loads on its bound-
aries. We then examine whether the stresses are such that shearing (large shear
strains) occurs anywhere in the soil. Other soil design problems involving slopes
and retaining structures are also analyzed by computing the stresses or strains
resulting within the soil from actions done on the boundaries of the soil mass. This
chapter covers the basic concepts of the mechanics of soils needed for understand-
ing these analyses. At the soil element level, coverage includes stress analysis,
strain analysis, shearing and the formation of slip surfaces, and laboratory tests
used to study the stress-strain response and shearing of soil elements. At the level
of boundary-value problems, coverage includes a relatively simple problem in soil
plasticity (Rankine states) and stresses generated inside a semi-infinite soil mass
by boundary loads.
113
saL00581_ch04.qxd 9/14/06 4:54 PM Page 114
114 The Engineering of Foundations
4.1 Stress Analysis

Elements (Points) in a Soil Mass
and Boundary-Value Problems
In this chapter, we are concerned with (1) how soil behaves at the element level
and (2) how this behavior, appropriately described by suitable equations and
ascribed to every element (that is, every point) of a soil mass subjected to certain
boundary conditions, in combination with analyses from elasticity or plasticity the-
ory, allows us to determine how boundary loads or imposed displacements result in
stresses and strains throughout the soil mass (and, in extreme cases, in strains so
localized and large that collapse of a part of the soil mass and all that it supports
happens).
At the element level, we must define precisely what the element is and how
large it must be in order to be representative of soil behavior. We must also have
equations that describe the relationship between stresses and strains for the ele-
ment and that describe the combination of stresses that would lead to extremely
large strains1. At the level of the boundary-value problem, we must define the max-
imum size an element may have with respect to characteristic lengths of the prob-
lem (for example, the width of a foundation) and still be treated as a point. This is
important because we use concepts and analyses from continuum mechanics to
solve soil mechanics problems, and so our elements must be points that are part of
a continuous mass. We must then have analyses that take into account how ele-
ments interact with each other and with specified stresses or displacements at the
boundaries of the soil mass in order to produce values of stress and strain every-
where in the soil mass. Further, we must also be able to analyze cases where the
stresses are such that large strains develop in localized zones of the soil mass. Con-
cepts from plasticity theory are used for that.
It is clear from the preceding discussion that the mechanics of the continuum
is a very integral part of soil mechanics and geotechnical and foundation engineer-
ing. We will introduce the concepts that are necessary gradually and naturally and
with a mathematical treatment that is kept as simple as possible.
Stress
A significant amount of the work we do as foundation engineers is based on the
concept of stress. Stress is a concept from the mechanics of continuous bodies.
Because soil is not a continuous medium, it is useful to discuss the meaning of
stress in soils. Consider a small planar area A passing through point P located
within a soil mass (Fig. 4-1). A normal force FN and a tangential force FT are
applied on A. If soil were a continuum, the normal stress s acting normal to A at
1 Extremely large strains, particularly extremely large shear strains, are closely tied to concepts of
rupture, yield, and failure (which term is used depends on the publication and the subject it deals with).
None of these three terms perfectly describes the range of problems we deal with, so we will introduce
appropriate terms throughout the chapter and the remainder of the text.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 115
CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 115
FN Figure 4-1
Definition of stress in soils: As the area A is allowed to
shrink down to a very small value, the ratios FN/A and
FT/A approach values s and t, the normal and shear
stress at P, respectively.
P
FT
A
point P would be defined as the limit of FN/A as A tends to the point P (that is,
tends to zero, centered around P). The shear stress would be defined similarly.
Mathematically,
FN
s limAS0 (4.1)
A
FT
t limAS0 (4.2)
A
Because soil is not a true continuum, we must modify this definition. A point
within a soil mass is defined as a volume V0 that is still very small compared with
the dimensions of the foundations, slopes, or retaining structures we analyze, but it
is sufficiently large to contain a large number of particles and thus be representa-
tive of the soil.2 To this representative volume V0, we associate a representative
area A0 (also very small, of a size related to that of V0). So we modify Eqs. (4.1)
and (4.2) by changing the limit approached by the area A from zero to A0:
FN
s limASA0 (4.3)
A
FT
t limASA0 (4.4)
A
The preceding discussion brings out one difference between soil mechanics
and the mechanics of metals, for example. In metals, the representative elementary
volume (REV) is very small. The REV for a given metal is indeed so small that, in
ordinary practice or introductory courses, we tend to think of it as being a point,
forgetting that metals are also made up of atoms arranged in particular ways, so
that they too have nonzero REVs, although much smaller than those we must use
in soil mechanics. In soil mechanics, our REVs must include enough particles and
the voids between them that, statistically, this group of particles will behave in a
way that is representative of the way larger volumes of the soil would behave.
2Mechanicians like to use the term “representative elementary volume” to describe the smallest volume
of a given material that captures its mechanical properties.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 116
Figure 4-2
Plane p11
(a) Elemental representation of two-dimensional
stress state at a point; (b) illustration of how the s33
element would distort when acted upon by a
positive s13 for the simple case with zero normal s13 tu su
stresses s11 and s33; (c) sectioned triangular
prismatic element, where su and tu depend on s11 s11 s11
u u
the angle u.
x3 Plane p33
s13 s13
s33 s33
x1 s13
(a) (c)
Vertical leg
of “L” parallel to p11
s13
Horizontal leg
of “L” parallel to p33
(b)
Stress Analysis
Stress analysis allows us to obtain the normal and shear stresses in any plane pass-
ing through a point,3 given the normal and shear stresses acting on mutually per-
pendicular planes passing through the point.4 We will see in Section 4.6 some
examples of how these stresses can be calculated at a point inside a soil mass from
a variety of boundary loadings common in geotechnical engineering. Later in the
text, we discuss some applications requiring corrections for three-dimensional (3D)
effects, but for now we will focus on two-dimensional (2D) stress analysis, which
turns out to be applicable to most problems of interest in practice.
Figure 4-2(a) shows a small, prismatic element of soil representing a “point”
within the soil. The faces of the element are aligned with the directions of the ref-
erence axes x1 and x3. The soil element is acted upon by the normal stresses s11
acting in the x1 direction and s33 acting in the x3 direction and the shear stresses
s13 (which, due to the requirement of moment equilibrium, is identical to s31).5
The first subscript of a stress component represents the direction normal to the
plane on which it acts; the second, the direction of the stress component itself. A
3 A point in the soil is indistinguishable, for our purposes, from a representative soil element, which has
a very small volume (and so is a point for practical purposes) but is sufficiently large to be representative
of the soil in its mechanical behavior.
4 A more proper definition of stress analysis for advanced readers would be that stress analysis aims to
allow calculations of the traction (which has normal and shear components) on a plane given the stress
tensor at the point.
5 Note that many engineering texts use the notation t instead of s when i j to represent shear stress.
ij ij
We have retained here the traditional mechanics terminology, in which s is used for both normal and
shear stresses.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 117
stress component with different subscripts is a shear stress; one with identical sub-
scripts, a normal stress. For example, s11 is the stress acting on the plane normal
to x1 in the x1 direction (that is, it is a normal stress), while s13 is the stress that
acts on the plane normal to x1 in the x3 direction (and is thus a shear stress). It is
simpler (although not required) to solve problems if we adopt the practice of
choosing s11 and s33 such that s11 s33. So if the normal stresses are 1,000 and
300 kPa, then s11 1,000 kPa and s33 300 kPa. Likewise, if the normal
stresses are 100 and 500 kPa, then s11 100 kPa and s33 500 kPa.
The plane where s11 acts is denoted by p11; likewise, p33 is the plane where
s33 acts. The stresses s11, s33, and s13 are all represented in their positive direc-
tions in Fig. 4-2(a). This means normal stresses are positive in compression, and
the angle u is positive counterclockwise from p11. With respect to the shear stress
s13, note that the prism shown in the figure has four sides, two representing plane
p11 and two representing plane p33. Looking at the prism from the left side, we
may visualize the plane p11 as the vertical leg of an uppercase letter “L” and p33
as the horizontal leg of the uppercase letter “L.” We then see that a positive shear
stress s13 acts in such a way as to open up (increase) the right angle of the “L”
formed by planes p11 and p33. Figure 4-2(b) shows the deformed shape that would
result for the element assuming zero normal stresses and positive s13.
If we section the element of Fig. 4-2(a) along a plane making an angle u with
the plane p11, as shown in Fig. 4-2(c), a normal stress su and a shear stress tu must
be applied to this plane to account for the effects of the part of the element that is
removed. While the sign of the normal stress su is unambiguous (positive in com-
pression), the shear stress on the sloping plane has two possible directions: up or
down the plane. So we must decide which of these two directions is associated with
a positive shear stress. The positive direction of the shear stress actually follows
from the sign convention already discussed (that s13 0 when its effect would be
to increase the right angle of the uppercase “L” made up by p11 as its vertical and
p33 as its horizontal leg). It turns out the shear stress tu is positive as drawn in the
figure, when it is rotating around the prismatic element in the counterclockwise
direction. We will show why this is so later, when we introduce the Mohr circle.
Our problem now is to determine the normal stress su and the shear stress tu
acting on the plane making an angle u with p11. This can be done by considering
the equilibrium in the vertical and horizontal directions and solving for su and tu.
The following equations result:
su 12 1s11 s33 2 12 1s11 s33 2cos 2u s13 sin 2u (4.5)
tu 12 1s11 s33 2sin 2u s13 cos 2u (4.6)
where the signs of s11, s33, and s13 are as discussed earlier (positive in compres-
sion for the normal stresses and determined by the “L” rule in the case of s13).
Principal Stresses and Principal Planes

Equations (4.5) and (4.6) tell us that su and tu vary with u. That means that the nor-
mal and shear stresses on each plane through a given point are a unique pair. There
saL00581_ch04.qxd 9/14/06 4:54 PM Page 118
will be two planes out of the infinite number of planes through the point for which
the normal stress will be a minimum and a maximum. These are called principal
stresses. They are obtained by maximizing and minimizing su by differentiating
Eq. (4.5) with respect to u and making the resulting expression equal to zero. The
largest principal stress is known as the major principal stress; it is denoted as s1.
The smallest principal stress is the minor principal stress, denoted as s3. The planes
where they act are referred to as the major and minor principal planes, denoted by
p1 and p3, respectively. When we differentiate Eq. (4.5) with respect to u and make
the resulting expression equal to zero, we obtain the same expression we obtain
when we make tu, given by Eq. (4.6), equal to zero. This means that the shear
stresses acting on principal planes are zero.
An easy way to find the angles up1 and up3 that the principal planes p1 and p3
make with p11 (measured counterclockwise from p11) is then to make tu 0 in
Eq. (4.6), which leads to
tan 12u p 2
2s13
(4.7)
s11 s33
When up is substituted for u back into Eq. (4.5), we obtain the principal stresses s1
and s3, which are the two normal stresses acting on planes where tu 0 and are also
the maximum and minimum normal stress for the point under consideration, given by
s1 12 1s11 s33 2 214 1s11 s33 2 2 s 213 (4.8)
s3 12 1s11 s33 2 214 1s11 s33 2 2 s 213 (4.9)
Note that, given the definition of the tangent of an angle, there is an infinite
number of values of up that satisfy Eq. (4.7). Starting with any value of up satisfy-
ing Eq. (4.7), we obtain additional values that are also solutions of (4.7) by repeat-
edly either adding or subtracting 90°. Using a calculator or a computer program,
the value of up calculated from Eq. (4.7) is a number between 90° and 90°. If s11
s33, we expect the major principal stress s1 to be closer in direction to s11 (the
larger stress) than to s33 (the smaller stress); so if the absolute value of the calcu-
lated value of the angle up is less than 45°, up up1; otherwise, up up3. Once the
direction of one of the principle planes is known, the direction of the other plane
can be calculated easily by either adding or subtracting 90° to obtain an angle with
absolute value less than 90°. For example, if up1 is calculated as 25°, then up3 is
equal to 65°. Alternatively, if up1 is found to be 25°, then up3 is calculated as
25° 90° 65°.
Mohr Circle
Moving 12 (s11 s33) to the left-hand side of Eq. (4.5), taking the square of both
sides of the resulting equation, and adding it to Eq. (4.6) (with both sides also
squared), we obtain
3su 12 1s11 s33 2 4 2 t 2u 14 1s11 s33 2 2 s 213 (4.10)

saL00581_ch04.qxd 9/14/06 4:54 PM Page 119
Recalling the equation of a circle in Cartesian
Line parallel to p11

coordinates, (x a)2 (y b)2 R2, where (a, b)
are the coordinates of the center of the circle and R t
is its radius, Eq. (4.10) is clearly the equation of a
circle with center C[(s11 s33)/2, 0] and radius
R 1 14 1s11 s33 2 2 s 213 in s-t space. Each
S3(s33, s13)
P Line parallel to p33
point of this circle, which is referred to as the Mohr
u
circle, is defined by two coordinates: the first is a S(su, tu)
normal stress (s) and the second is a shear stress (t).
2u s
Figure 4-3 shows the Mohr circle for the stress state
(s11, s33, s13) of Fig. 4-2(a). According to Eqs. (4.5) (s3, 0)
and (4.6), the stresses su and tu on the plane making (s1, 0)
an angle u measured counterclockwise from plane S1(s11, s13)
p11 are the coordinates of a point on the Mohr circle
rotated 2u counterclockwise from point S1 (point S1 Figure 4-3
represents the stresses on p11). Note that the central Mohr circle corresponding to state of stress shown in
angle of the Mohr circle separating S1 from S3, Fig. 4-2.
which is 180°, is indeed twice the 90° angle separat-
ing the corresponding planes, p11 and p33.
In Fig. 4-3, the stresses (s11, s33, s13) plot as two points S1(s11, s13) and
S3(s33, s13), diametrically opposed. The figure illustrates the case of s11 s33. To
understand why the shear stress s13 plots as a positive number with s33 and as a
negative number with s11, refer to Fig. 4-2(b). With the positive values of s11, s33,
and s13 shown in Fig. 4-2(a), the prismatic element deforms as shown in Fig.
4-2(b), with the direction of maximum compression associated with the major
principal stress s1. It is clear that the major principal plane, which is normal to the
direction of maximum compression, is obtained by a rotation of some angle up1
90° counterclockwise with respect to p11. This means that, in the Mohr circle, we s
must have a counterclockwise rotation 2up1 from the point S1(s11, s13) associ-
ated with p11 in order to reach the point (s1, 0) of the circle. This implies that S1
must indeed be located as shown in Fig. 4-3, for if we had S1(s11, s13) instead of
S1(s11, s13), a counterclockwise rotation less than 90° would not take us to (s1,
0). So this means s13 is plotted as negative if spinning clockwise around the ele- Positive normal
ment, as it does for plane p11, and as positive if spinning counterclockwise, as it stress
does for plane p33. This is the basis for the convention we will use for plotting
points in the Mohr circle: Shear stresses rotating around the element in a counter- t
clockwise direction are positive; they are negative otherwise (Fig. 4-4). Note that
this sign convention is not independent from but actually follows directly from the
shear stress sign convention we adopted for the stress s13 appearing in Eqs. (4.5)
and (4.6). Positive shear stress
(for plotting in Mohr
diagram)
Pole Method
Figure 4-4
Mohr circles have interesting geometric properties. One very useful property is that Stress sign
the central angle of a circle corresponding to a certain arc is twice as large as an conventions.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 120
Figure 4-5 t
The geometric property of circles that a
central angle 2u produces the same arc as
an inscribed angle u.
u 2u
inscribed angle corresponding to the same arc (Fig. 4-5). Applying this property to
the Mohr circle shown in Fig. 4-3, the angle made by two straight lines drawn
from any point of the circle to point S(su, tu) representing the stresses on the plane
of interest and to point S1(s11, s13) is equal to u. In particular, there is one and
only one point P on the circle with the property that the line joining P to point
S(su, tu) is parallel to the plane on which su and tu act. The point P with this
property is known as the pole of the Mohr circle. Based on the preceding discus-
sion, the pole can be defined as the point such that, if we draw a line through the
pole parallel to the plane where stresses su and tu act, this line intersects the Mohr
circle at a point whose coordinates are su, tu.
In order to determine the pole, we need to know the orientation of at least one
plane where the stresses are known. We can then use the known stresses (s, t) to
find the pole by drawing a line parallel to the plane acted upon by these stresses.
This line intersects the circle at a point; this point is the pole. Once we know the
pole P, we can determine the stresses on any plane by drawing a straight line
through the pole P parallel to the plane where the stresses are desired. This line
intersects the circle at a point whose coordinates are the desired stresses.
We can use Figs. 4-2 and 4-3 to illustrate the concept of the pole. Consider the
element of Fig. 4-2(a). By plotting the Mohr circle for this state of stress, we
obtain the expected diametrically opposed points S1 and S3 shown in Fig. 4-3. If
we look at point S1 on the circle and consider the corresponding stresses shown in
Fig. 4-2(a), we can easily determine the location of the pole. If we construct a line
through S1 that is parallel to the plane p11 where (s11, s13) acts, we will deter-
mine the pole as the point where this line intersects the Mohr circle. In this case,
since we are dealing with u 0°, the line is vertical, and the pole (point P in Fig.
4-3) lies directly above S1 (also shown in Fig. 4-3). Likewise, if we look at point S3
and construct a line through S3 parallel to plane p33 on which (s33, s13) acts, we
can also determine the pole as the point of intersection of this horizontal line with
the Mohr circle. As expected, the pole is found to be directly to the right of S3 and
to coincide with the point determined previously by examining point S1. This
shows clearly that the pole is unique for a given stress state. Figure 4-6 illustrates
for the same case how the principal directions and principal stresses would be
determined once the pole is known.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 121
t
Line parallel to p11 Figure 4-6
Illustration of the relationship of the principal directions
and their representation in a Mohr diagram.
Line parallel to p3
up3
(s3, 0) P
(s1, 0)
s
up1
Line parallel to p1
EXAMPLE 4-1
A state of stress is represented by the block in Fig. 4-7(a). Determine the location of the
pole and give its coordinates in the s-t system.
■ Solution
First, plot the Mohr circle [Fig. 4-7(b)]. Next, determine the location of the pole. Using
point (200, 100), draw a line parallel to the plane on which s11 (200 kPa) acts. In this case,
it is a vertical line since p11 is vertical. Likewise, if the other point is chosen (100, 100),
we draw a line parallel to the plane on which s33 (100 kPa) acts. In this case, it is a horizon-
tal line since p33 is horizontal. Either method produces the location of the pole: (200,
100).
t
100 kPa
100 kPa
(200, 100)
200 kPa 200 kPa
100 200 s
100 kPa
100 kPa Pole
(a) (b)
Figure 4-7
(a) Stress state for Example 4-1 and (b) corresponding Mohr circle.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 122
Solving Stress Analysis Problems

The steps in solving a 2D stress analysis problem can be outlined as follows:
1. Choose the largest normal stress as s11. For example, if the normal stresses
are 300 and 100 kPa, then s11 300 kPa and s33 100 kPa. Likewise, if
the normal stresses are 100 and 300 kPa, then s11 100 kPa and s33
300 kPa. The plane where s11 acts is denoted by p11, while the one where
s33 acts is denoted by p33.
2. If using Eq. (4.5) or (4.6), assign a sign to s13 based on whether it acts to
increase or decrease the angle between p11 and p33 (with p11 being the vertical
leg of the “L”; see Fig. 4-2). The stress s13 is positive if it acts in a way that
would tend to increase the angle.
3. Recognize that the reference plane for angle measurements is p11 and that
angles are positive counterclockwise.
4. Reason physically to help check your answers. For example, since s11 s33,
the direction of s1 will be closer to that of s11 than to that of s33. Whether s1
points up or down with respect to s11 now depends on the sign of s13, for it
tells us about the directions in which the element is compressed and extended.
Naturally, s1 acts in the direction in which the element is compressed the most.
EXAMPLE 4-2
The state of stress at a point is represented in Fig. 4-8. Find (a) the principal planes, (b) the
principal stresses, and (c) the stresses on planes making angles
15° with the horizontal.
Solve both analytically and graphically.
■ Solution
Analytical Solution
Take s11 200 kPa and s33 50 kPa. So p11 makes an angle equal to 30° with the hor-
izontal, and p33 makes an angle equal to 60° with the horizontal. In order to assign a sign
to s13, we need to consider the right angle made by p11 and p33; we must look at this angle
200 kPa
50 kPa
50 kPa
50 kPa
50 kPa
200 kPa
Figure 4-8
State of stress at a point
30
(Example 4-2).
saL00581_ch04.qxd 9/14/06 4:54 PM Page 123
as if it were an uppercase letter “L,” such that p11 is the “vertical” and p33 is the “horizon-
tal” leg of the letter “L.” Physically rotating the page until we see the “L” may be helpful in
this visualization. The effect of the shear stress on the right angle made by p11 and p33
looked at in this manner is to reduce it; accordingly, s13 50 kPa. We are now prepared
to solve the problem.
Principal Planes
Substituting the values of s11, s33, and s13 into Eq. (4.7):
2 1502
tan 12u p 2
2

200 50 3
from which
arctan a b 16.8°
1 2
up
2 3
The absolute value of 16.8° is less than 45°, so

u p1 16.8°
u p3 16.8° 90° 73.2°
Graphically, up1 would show as an angle of 16.8° clockwise from p11 because the calcu-
lated angle is negative. If the answer is desired in terms of the angles that p11 and p3 make
with the horizontal, then we need to add 30° (the angle that p11 makes with the horizontal) to
these two results: p11 is at an angle 13.2°, and p3 is at 103.2° (or 76.8°) to the horizontal.
Principal Stresses
The principal stresses can be calculated using either Eq. (4.5) with u 16.8° and 73.2° or
Eqs. (4.8) and (4.9). Using Eq. (4.5):
s1 12 1200 50 2 12 1200 50 2cos32 116.8°2 4 1502sin 32 116.8°2 4 215.1 kPa
s3 12 1200 50 2 12 1200 50 2cos32 173.2°2 4 1502sin 32 173.2°2 4 34.9 kPa
Now using Eqs. (4.8) and (4.9):
s1 12 1200 502 214 1200 50 2 2 1502 2 215.1 kPa
s3 12 1200 502 214 1200 50 2 2 1502 2 34.9 kPa
Stresses on Planes Making Angles 15° with Horizontal

These planes make angles 15° and 45° with p11, respectively. Plugging these values
(15° and 45°) into Eqs. (4.5) and (4.6):
su 125 75 cos32 115°2 4 50 sin 32 115°2 4 214.9 kPa
tu 75 sin 3 2 115°2 4 50 cos32 115°2 4 5.8 kPa
Not surprisingly, the su and tu values calculated here are very close to the values
for the major principal plane (215.1 and 0). You should verify that su 175 kPa and tu
75 kPa for the other plane (which makes an angle of 15° with the horizontal).
saL00581_ch04.qxd 9/14/06 4:54 PM Page 124
Figure 4-9 t
Mohr circle and graphical solution of
Example 4-2.
100
p3
(200, 50) 15
50
(214, 6)
100 p1
200 s
P
(175, 75)
(50, 50)
15
Graphical Solution Using the Pole Method

The solution can be seen in Fig. 4-9. The normal and shear stresses on plane p11 are 200
and 50 kPa, respectively; in plane p33 they are 50 and 50 kPa.6 These two points are dia-
metrically opposite each other on the Mohr circle. If we plot these two points in s-t space
and join them by a straight line, this line crosses the s-axis at the center of the circle. We
can then easily draw the Mohr circle using a compass. The principal stresses are now easily
read as the abscissas of the two points with t 0.
If we now draw a line parallel to p11 through the point (200, 50), we obtain the pole P
as the intersection of this line with the circle. Note that if we draw a line parallel to p33
through (50, 50), we obtain the same result. The directions of p1 and p3 are obtained by
drawing lines through P to the points (s1, 0) and (s3, 0) of the Mohr circle, respectively.
The stresses at
15° with the horizontal are found by drawing lines through the pole P
making
15° with the horizontal. These lines intersect the circle at two points with coordi-
nates (214, 6) and (175, 75), respectively.
Total and Effective Stresses

When we plot Mohr circles (as in Fig. 4-10), we are representing the state of stress
at a point in the soil mass. If there is a nonzero pore pressure u at this point, it is
the same in every direction, and thus it does not affect the equilibrium of the point.
We should remember that water cannot sustain shear stresses, so the presence of a
pore pressure affects only normal stresses in the soil. It is useful to examine what
would happen if we plotted both a total stress and an effective stress Mohr circle.
Let us assume Eq. (4.10) applies to total stresses. Referring back to Eq. (3.3), if we
substitute it into Eq. (4.10), we can see that the only effect is to have the Mohr cir-
6 Note that, for Mohr circle construction, counterclockwise shear stresses are plotted as positive, while
clockwise shear stresses are plotted as negative.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 125
Pore pressure u
(s u, t) (s, t)
s, s
Effective-stress Total-stress
Mohr circle Mohr circle
Figure 4-10
Principle of effective stresses: illustration of the difference between the total and the effective
stress state as represented using the concept of the Mohr circle.
cle of total stresses displaced along the s-axis with respect to the Mohr circle of
effective stresses by an amount equal to the pore pressure u.
4.2 Strains*
Definitions of Normal and Shear Strains
Analysis of geotechnical problems cannot always be done only in terms of
stresses. These stresses induce deformations, which are represented by strains.
There are two types of strains: normal and shear strains. At a given point, a normal
strain in a given direction quantifies the change in length (contraction or elonga-
tion) of an infinitesimal linear element (a very small straight line) aligned with that
direction. The so-called engineering shear strain g is a measure, at a given point, of
the distortion (change in shape). With respect to the two reference axes x1 and x3,
the engineering shear strain g13 expresses the increase of the initial 90° angle
formed by two perpendicular infinitesimal linear elements aligned with these axes.
A shear strain, as it is defined in mechanics, is one half the value of the engineer-
ing shear strain.7
Both normal and shear strains eij can be expressed through
e ij a b
1 0u i 0u j
(4.11)
2 0xj 0xi
where i,j 1, 2, 3 reference directions, xi coordinate in the i direction, ui

displacement in the xi direction. The subscripts indicate the directions of linear
7 We have avoided duplication of symbols as much as possible, but there is no good alternative to using
the traditional notation for engineering shear strain, which of course is the same as used for unit weight.
The reader should observe the context in which symbols are used to avoid any confusion.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 126
differential elements and the directions in which displacements of the end points of
the differential elements are considered. When i j, the strain is a normal strain; it
is a shear strain otherwise.
The relationship between a shear strain eij and the corresponding engineering
shear strain gij is
gij 2e ij (4.12)
EXAMPLE 4-3
B B* Derive, in a simple way, the expression for the normal strain at a point in the direction of
A
dx1 du1 reference axis x1.
■ Solution
x3
Let’s consider the case of Fig. 4-11. For an underformed soil mass, we have a differential
element dx1 aligned with the x1 direction (Fig. 4-11). We have labeled the initial point of the
segment A and the end point B. In drafting this figure, we have corrected for rigid body
translation in the x1 direction. In other words, we are plotting the deformed element as if the
x1 displacement of A were zero for easier comparison with the original, undeformed element.
This way, every displacement in the figure is relative to the displacement of A. If, after the
Figure 4-11 soil mass is deformed, point B moves more in the positive x1 direction than point A (that is,
Normal strain e11: if the displacement u1 of B is greater than that of A), as shown in Fig. 4-11, then the element
infinitesimal element
has clearly elongated (this elongation is seen in the figure as B* B).
dx1 shown after
Taking some liberty with mathematical notation, the unit elongation of dx1 is the differ-
correction for rigid
body motion both ence in displacement between points A and B (u1B u1A B* B du1 u1) divided
before deformation by the initial length of the element (x1), or u1/x1. It remains to determine whether elonga-
(AB) and after tion is a positive or negative normal strain. To be consistent with the sign convention for
deformation (AB*). stresses, according to which tensile stresses are negative, our normal strain e11 in the x1
direction must be
0u 1
e 11
0x1
Note that this, indeed, is the expression that results directly from Eq. (4.11) when we make
i 1 and j 1.
EXAMPLE 4-4
Derive, in a simple way, the expression of the shear strain at a point in the x1-x3 plane.
■ Solution
Consider two differential linear elements, dx1 and dx3, aligned with the x1- and x3-axis,
respectively, at a point within an undeformed soil mass (Fig. 4-12). Now consider that the
soil mass is deformed, and, as a result, points B and C (the end points of elements dx1 and
dx3, respectively) move as shown (to new positions labeled B* and C*) with respect to point
A (note that, as for Example 4-3, we are not representing rigid-body translation in the x1 and
saL00581_ch04.qxd 9/14/06 4:54 PM Page 127
x3 directions in the figure). We can see that both point B and point C have displace- C*
ments that have components in both the x1 and x3 directions. Here we are interested in C
just the distortion of the square element made up of dx1 and dx3, not in the elongation
∂u1
or shortening of dx1 and dx3 individually. The distortion of the element clearly results
∂x3
from difference in the displacement u3 in the x3 direction between A and B and in the
dx3 ∂u3
displacement u1 in the x1 direction between A and C. ∂x1
Taking again some liberties with mathematical notation, we can state that the dif- B*
ferences in the displacements of points B and A (in the x3 direction) and C and A (in A dx1 B
the x1 direction) can be denoted as u3 and u1, respectively. Since the deformations we x3
are dealing with are small, the angle by which the differential element dx1 rotates coun-
terclockwise is approximately equal to u3 divided by the length of the element itself,
or u3/x1; similarly, the angle by which dx3 rotates clockwise is u1/x3. These two
rotations create a reduction in the angle between the elements dx1 and dx3, which was x1
originally 90°, characterizing a measure of distortion of the element. If we add them
together, we obtain the absolute value of what has become known as the engineering Figure 4-12
shear strain g13; one-half the sum gives us the absolute value of e13. It remains to deter- Shear strain e13:
mine whether this is a positive or negative distortion. Our sign convention for shear infinitesimal square
strains must be consistent with our shear stress sign convention. Recall from our earlier element shown after
discussion in the chapter that a positive shear stress was one that caused the 90° angle correction for rigid body
of our square to open up (to increase), not to decrease. Therefore, we will need a nega- translation both before
tive sign in front of our sum in order to obtain a negative shear strain for the reduction deformation (defined by
of the 90° angle we found to take place for the element in Fig. 4-12: dx1 AB and dx3 AC)
and after deformation
(defined by dx1* AB*
e 13 a b
1 0u 1 0u 3

2 0x3 0x1 and dx3* AC*).
Note that this equation results directly from Eq. (4.11) when we make i 1 and j 3
(or vice versa).
Strains as expressed by Eq. (4.11) are small numbers; there are other ways of
defining strain that are more appropriate when elongations, contractions, or distor-
tions become very large. However, Eq. (4.11) may still be used for increments of
strains, even if the cumulative strains measured from an initial configuration are
very large. It is appropriate in that case to use a d (the symbol for differential)
before the strain symbol (as in de and dg) to indicate that we refer to a strain
increment.
As was true for stresses, there are also principal strains e1 and e3 (and principal
strain increments de1 and de3). These are strains (or strain increments) in the direc-
tions in which there is no distortion, and the shear strains (or shear strain increments)
are equal to zero. Distortion happens any time the shape of an element changes. It is
important to understand that a point in the soil experiences distortion as long as de1
de3. To illustrate this, Fig. 4-13(b) and (c) show two alternative elements repre-
senting a point P. The larger, outer element is aligned with the principal directions,
and the deformed shape of the element [shown in Fig. 4-13(a)] does not immediately
suggest the notion of distortion. However, there is no doubt that distortion has taken
place when we observe the change of shape of the smaller, inner element as we go
from the initial to the final configuration [Fig. 4-13(b) and (c)].
saL00581_ch04.qxd 9/14/06 4:54 PM Page 128
The volumetric strain increment dev, defined as minus the change in volume
divided by the original volume (the negative sign being required to make contrac-
tion positive), can be easily determined in terms of a cubic element with sides with
length initially equal to 1 and aligned with the principal directions (which means
x1, x2, and x3 are principal directions). The element is then allowed to expand as a
result of elongations equal to du1, du2, and du3 in the three reference directions. As
the cube is aligned with the principal directions, there will be no distortion in the
planes x1x2, x1x3, or x2x3. It is apparent from Fig. 4-14 that
dV 11 du 1 2 11 du 2 2 11 du 3 2 1
Referring back to our definition of normal strain and considering that the initial
length of the sides of the cube are of unit length and the initial volume of the cube
is also equal to 1, we can write the following for the volumetric strain increment:
1 11 de 1 2 11 de 2 2 11 de 3 2
dV
de v
1
which, given that the strain increments are very small (and that second- and third-
order terms would be extremely small and thus negligible), reduces to
de v de 1 de 2 de 3 (4.13)
The volume change at a point is clearly independent of the reference system
and of any distortion, so the following equation would also apply even if x1, x2,
and x3 were not the principal directions:
de v de 11 de 22 de 33 (4.14)
where de11, de22, and de33 normal strains in the arbitrary directions x1, x2, and x3.
Initial outer element
de3 1
P P P de1
1
Final outer element Initial configuration Final configuration

(a) (b) (c)
1
Figure 4-13 de2
Alternative representations for the state of strain at a point: (a) an x3
Deformed cube
element aligned with the principal strain directions (vertical Original cube
contraction and horizontal elongation); (b) the same element before x2
deformation with an element inside it with sides oriented at 45° to x1
the principal strain directions; (c) the same element after
deformation, showing the distortion of the element with sides not Figure 4-14
aligned with the principal strain directions. Calculation of volumetric strain.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 129
Mohr Circle of Strains

The mathematics of stresses and strains is the same: Normal strains play the same
role as normal stresses, and shear strains, the same as shear stresses. So, just as it
was possible to express the stresses at a point under 2D conditions in terms of a
Mohr circle, the same is possible for strains. Many problems in geotechnical engi-
neering can be idealized as plane-strain problems, for which the strain in one direc-
tion is zero. For example, slopes and retaining structures are usually modeled as rel-
atively long in one direction, with the same cross section throughout and with no
loads applied in the direction normal to the cross sections. Except for cross sections
near the ends of these structures, it is reasonable, based on symmetry considera-
tions, to assume zero normal strain perpendicular to the cross section. The same is
valid for strip footings, which are used to support lines of columns or load-bearing
walls. The result is that we can take de2 0 and do our strain analysis in two
dimensions.
In the case of the Mohr circle of strains, incremental strains, not total strains,
are plotted in the horizontal and vertical axes. Figure 4-15 shows a Mohr circle of
strains plotted in normal strain increment de versus shear strain increment 12 dg
space. As was true for stresses, each point of the Mohr circle represents one plane
through the point in the soil mass for which the Mohr circle represents the strain
state. The points of greatest interest in the Mohr circle of strains are
■ The leftmost and rightmost points, (de3, 0) and (de1, 0), corresponding to the
minor and major principal incremental strain directions
■ The highest and lowest points, (12 dev,
12 dgmax), corresponding to the directions
of largest shear strain (Note that 12 dev 12 (de1 de3) is the de coordinate of the
center and 12 dgmax 12 (de1 de3) is the radius of the Mohr circle.)
■ The points where the circle intersects the shear strain axis, (0,
12 dgz)
The two points with zero normal strain increment correspond to the two
directions along which de 0, that is, the directions along which there is neither
extension nor contraction. It is possible to define a separate reference system for
1

2 dg
1 dgz)
A1 (0,
2
1
Pole P O(
2 dev, 0) c
c
(de3, 0) Z (de1, 0) de
1 dgz)
A2 (0,
2
Figure 4-15
Mohr circle of strains.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 130
each of these two directions such that x1 in each system is aligned with the direc-
tion of zero normal strain. To clearly indicate that x1 is a direction of zero normal
strain, we can use a superscript z, as in x z1. This will be useful in our discussion of
the dilatancy angle, which follows.
Dilatancy Angle
The angle c shown in Fig. 4-15, known as the dilatancy angle, is quite useful in
understanding and quantifying soil behavior. There are two ways of expressing the
dilatancy angle based on the geometry of the Mohr circle of strains:
2 1de 1 de 3 2
1
OZ de 1 de 3 de v
0 OA 1 0 0 dgmax 0
sin c 1
2 1de 1 de 3 2
(4.15)
de 1 de 3
1
OZ 2 de v dev
0 ZA1 0 0 dg 0
tan c 1 z z
0 2 dg 0
(4.16)
where gz shear strain in the x z1-x zp z

3 plane (Fig. 4-16), x 1 is a direction of zero nor-
zp
mal strain, and x 3 is the direction normal to x1.z
The dilatancy angle is clearly related to the volumetric strain resulting from a
unit increase in shear strain.8 By definition, the dilatancy angle c is positive when
there is dilation (volume expansion). This is apparent from Eqs. (4.15) and (4.16),
for the dilatancy angle clearly results positive when volume expands, that is, when
dev 0. Note that the denominators of Eqs. (4.15) and (4.16) are always positive
(hence the absolute values taken), for the dilatancy angle is related to the volumet-
ric strain increment resulting from a unit increment in shear strain, regardless of
the orientation of the shear strain. In other words, the dilatancy angle would still be
the same positive value if the element shown in Fig. 4-16 were sheared to the left
and not to the right as shown. The direction of zero normal strain in Fig. 4-16 is
gz
x3zp
Figure 4-16
State of strain visualized for an
element with one side aligned with x1z Direction of zero
the direction of zero normal strain. normal strain
8 Technically, both the shear and volumetric strain increments in the definition of the dilatancy angle are
plastic strain increments, a distinction that for our purposes is not necessary to make.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 131
represented by x z1, and the direction normal to it, by x zp 3 , where the superscript zp
means that the direction x zp
3 is perpendicular to the direction of zero normal strain.
There are in fact two distinct directions of zero normal strain, as will be shown
later.
The Mohr circle of Fig. 4-15 corresponds to a state of dilation (expansion), as
dev 0. Examining the state of deformation in the x z1-x zp 3 reference system, expan-
sion implies that the normal strain in the x zp 3 direction, normal to x z1, is negative
(that is, that elongation takes place in the x zp
3 direction, as is clearly shown in Fig.
4-16). If we locate the pole in Fig. 4-15 and then draw two lines (one through A1
and one through A2), the two lines perpendicular to these two lines are the direc-
tions of zero normal strain, as we will discuss in detail later.
EXAMPLE 4-5
A soil element is subjected to the following incremental strains: de1 0.03%, de3
0.05%. Knowing that plane-strain conditions are in force (that is, the strain in the x2 direc-
tion is zero), calculate the dilatancy angle.
■ Solution
Because we know the principal strain increments, we can immediately calculate the dila-
tancy angle as
dev de1 de3 0.03 0.05
0 dgmax 0
sin c 0.25
de1 de3 0.03 10.052
from which
c 14.5°
In Problem 4-18, you are asked to continue this by plotting the Mohr circle, finding the
pole for the case when the major principal strain increment is vertical and determining the
directions of the potential slip planes through this element (which is the subject of a subse-
quent section).
For a triaxial strain state, in which e2 e3, the dilatancy angle needs to be
redefined because a maximum engineering shear strain increment dgmax is no
longer possible to define with clarity. In place of it, we work with des, defined as
de s de 1 2de 3 (4.17)
The volumetric strain in the triaxial case follows from Eq. (4.13):
de v de 1 2de 3
Thus, the dilatancy angle for triaxial conditions is written as
de 1 2de 3
sin c (4.18)
de 1 2de 3
saL00581_ch04.qxd 9/14/06 4:54 PM Page 132
A single expression for it, which applies to both plane-strain and triaxial
conditions, is
dev
de1 kde3 dev de1
sin c (4.19)
de1 kde3 de1 kde3 dev
2
de1
where
k e
1 for plane-strain conditions
(4.20)
2 for triaxial conditions
4.3 Failure Criteria, Deformations,

and Slip Surfaces
Mohr-Coulomb Failure Criterion
Soils are not elastic. Referring to the top stress-strain plot of Fig. 4-17, if we apply
repeatedly the same increment of shear strain to an element of soil (this is referred
to as strain-controlled loading), the increment of stress that the soil element is able
to sustain decreases continuously (a process that is sometimes referred to as modu-
lus degradation), until a state is reached (represented by point F in the figure) at
which the stress increment will be zero. At this point, if we continue to increase
the strain, the stress will stay the same. The other stress-strain plot shown in the
figure illustrates another possible response, whereby the stress peaks at point F.
This second response (referred to as strain softening) is common in soils. The lim-
iting or peak stress associated with point F in each case is usually what is meant
by the shear strength of the soil, as it is the maximum stress the soil can take.
Alternatively, we could have started loading the soil element by applying
stress increments to the element (which is referred to as stress-controlled loading).
In this case, when the point F at which the peak stress was observed during strain-
controlled loading is reached, the soil element will not be able to take any addi-
tional stress, and what will happen instead is uncontrolled deformation. It is not
F
Shear stress
c
Figure 4-17
b
Nonlinearity of stress-strain relationship for soils and failure (the
onset of very large deformations at some value of stress, represented
by point F). Note that the first strain increment a generates the stress
increment b and that the second strain increment generates a stress
increment c b b. a 2a Shear strain
saL00581_ch04.qxd 9/14/06 4:54 PM Page 133
possible in the stress-control case to plot the stress-strain relationship after point F
is reached.
The Mohr-Coulomb failure (or strength) criterion has been traditionally used
in soil mechanics to represent the shear strength of soil. It expresses the notion that
shear strength of soil increases with increasing normal effective stress applied on
the potential shearing plane. Although in stress analysis we do not have to concern
ourselves with whether we are dealing with effective or total stresses because the
analysis applies to both, we must now make a clear distinction. Soils “feel” only
effective stresses (that is, any deformation of the soil skeleton happens only in
response to effective stress changes); thus, the response of soil to loading and the
shear strength of soil depends only on the effective stresses. Based on this consid-
eration, we represent the Mohr-Coulomb criterion in s-t space as two straight
lines making angles
f with the horizontal and intercepting the t-axis at distances
c and c from the s-axis. Figure 4-18 shows only the line lying above the s-
axis, since the diagram is symmetric about the s-axis. The distance c is usually
referred to as the cohesive intercept. Mathematically, the Mohr-Coulomb criterion
may be represented in a simple way by
s c s¿ tan f (4.21)
where s is the shear strength of the soil and f is its friction angle. For a given nor-
mal effective stress s on a plane, if the shear stress on the plane is t s as given
by Eq. (4.21), shearing or what is commonly referred to in engineering practice as
“failure” of the material occurs. Failure here means the occurrence of very large
strains in the direction of that plane. This means that the soil cannot sustain shear
stresses above the value given by Eq. (4.21). Equation (4.21) is a straight line in
s-t space (as shown in Fig. 4-18) that is referred to as the Mohr-Coulomb
strength envelope. There can be no combination of s and t that would lie above
the Mohr-Coulomb strength envelope.
Figure 4-19 shows the Mohr circle for a soil element (or point) within a soil
mass, at failure, where s and t act on the two planes corresponding to the two
tangency points between the circle and the envelope (again, only the half of the
diagram lying above the s-axis is shown, as the part below the s-axis is symmet-
ric). All other points, representing all the planes where (s, t) do not satisfy Eq.
(4.21), lie below the Mohr-Coulomb envelope. If we know the directions of the
c
Figure 4-18
s The Mohr-Coulomb failure envelope.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 134
Figure 4-19 t
Mohr circle at failure.
ts
f
f c P 90 f
(s3, 0) (s1 , 0) s
45 f/2
1 (s s ) 1 (s s )
c cot f
2 1 3 2 1 3
principal stresses, we can determine the pole for this circle. Assuming a vertical
major principal stress and a horizontal minor principal stress, the pole P lies at
point (s3, 0). By simple geometry, the central angle corresponding to the arc
extending from (s1, 0) to the point of tangency is 90° f. This means the planes
corresponding to the points of tangency lie at
(45° f/2) with the horizontal.
The direction of a real shear (slip) surface passing through the soil element can be
approximated by one of these two possible directions. Note that this refers to the
direction of the slip surface at the point under consideration, and that at other
points of the soil mass the slip surface direction may be different because the prin-
cipal directions at those points may be different. This means that a slip surface
through a soil mass is a surface tangent at every point to the direction estimated as
described earlier, a direction making an angle of
(45° f/2) with the direction
of the major principal plane at the point. In a subsequent subsection, we discuss in
more detail the issue of slip surfaces, their nature, and their geometry.
It is possible to find the relationship between the principal stresses s1 and s3
at failure from the geometry of Fig. 4-19. It is easier to proceed if we define new,
transformed normal stresses s* through
s* s¿ c cot f (4.22)
Taking Eq. (4.22) into Eq. (4.21) leads to
s s* tan f (4.23)
This transformation (sometimes referred to as Caquot’s principle after the per-
son who first made use of it) is represented graphically in Fig. 4-20. A Mohr circle at
failure is also represented. We can write sin f in terms of the ratio of the radius of
the Mohr circle to the distance from the center of the Mohr circle to the origin of the
transformed system of stress coordinates (that is, where the s*- and t*-axes cross):
s*1 s*3
sin f
s*1 s*3
saL00581_ch04.qxd 9/14/06 4:54 PM Page 135
t* t Figure 4-20
Caquot’s principle: transformed
normal stresses.
1 (s* s*)

2 1 3
f
s, s*
c cot f
−12 (s1* s*3)
The principal stress ratio s*1 /s3* easily follows:

s1*
N (4.24)
s3*
where N is known as the flow number, given by
1 sin f
N (4.25)
1 sin f
If the relationship between the principal stresses in its original form is needed,
we just need to use Eq. (4.22) to rewrite s*1 and s3* in Eq. (4.24):
s¿1 c cot f 1 sin f
N
s¿3 c cot f 1 sin f
which can be rewritten as
s¿1 Ns¿3 1N 12c cot f (4.26)
This expression can be further rewritten as
s¿1 Ns¿3 2c2N (4.27)
by recognizing that
1 sin2 f
cot f
B sin2 f
Note that, for c 0, Eq. (4.26) reduces to
s¿1 Ns¿3 (4.28)
Slip Surfaces*
Figure 4-21(a) shows a slip surface and an element of it, which is expanded in Fig.
4-21(b). A slip surface (also referred to in the literature as a failure surface or the
saL00581_ch04.qxd 9/14/06 4:54 PM Page 136
Figure 4-21 Q Rigid blocks

Details of slip surface
element. dx3
1
Slip surface
Zero-extension line dx1
Expanded in (b)
zp
(a) x3
z
x1
(b)
more technical “shear band”) develops when the shear stresses on every point of
it reach the corresponding shear strengths. The slip surface is a surface along
which a soil mass slides with respect to another, usually stationary, soil mass.
These soil masses behave very much like rigid blocks, not undergoing any defor-
mation. Accordingly, an acceptable model for the slip surface element shown in
Fig. 4-21(b) has the element bounded by two rigid blocks. In truth, slip surfaces
are not truly surfaces, but very thin zones or bands of highly concentrated shear
strains, hence the representation of a slip surface element in Fig. 4-21 as having
nonzero thickness.
The slip surface is bounded by rigid blocks that are on the verge of sliding past
one another. On the onset of failure, these rigid blocks are connected to the slip sur-
face and, being rigid, prevent any contraction or elongation in the direction of the
slip surface. This means the normal strain in the direction x z1 of the slip surface is
equal to zero. The slip surface x z1 is, accordingly, a zero-extension line. The state of
strain at a point of the slip surface, referred to the axis x z1 parallel to the slip surface
and another x zp
3 normal to it, is represented by Fig. 4-16. Taking the thickness of the
slip surface as being equal to 1, the state of strain for the slip surface element of
Fig. 4-21 can be expressed as follows:
de z11 0
dx3
33
de zp dx3
1
0 dg z 0
dx1
dx1
1
de v de z11 de zp
33 dx 3
saL00581_ch04.qxd 9/14/06 4:54 PM Page 137
Recalling the definition of the dilatancy angle c given in Eq. (4.16), we can
write
de v dx3
0 dg 0
tan c z (4.29)
dx1
Equation (4.29) shows that the dilatancy angle represents the angle of the motion
at the top of the element with respect to the horizontal. Equation (4.29) is instru-
mental in understanding the geometry of slip surfaces. Take, for example, the soil
slope of Fig. 4-22. It is common to see such slopes failing along curved slip sur-
faces. The shape of a potential slip surface can be determined if we take a center of
rotation such that every point of the surface is at a variable distance r from the cen-
ter of rotation. Once a relationship between r and the angle of rotation u around
the center of rotation is defined, the shape of the slip surface becomes known.
An infinitesimal rotation du with respect to the center of rotation corresponds
to a tangent displacement rdu along the slip surface and to a possible small
increase of the radius by an amount dr. So dr is analogous to the dx3, and rdu is
analogous to the dx1 of Eq. (4.29). Taking these into Eq. (4.29) gives
dr
tan c (4.30)
rdu
With knowledge of one point of the slip surface, defined by r0 and u0, integra-
tion of Eq. (4.30) leads to
r r0e1uu02 tan c (4.31)
where r0 is the radius corresponding to a reference angle u0.
Equation (4.31) is a geometrical shape referred to as a logarithmic spiral (or
log-spiral, for short). The implication is that a homogeneous soil mass, with the
same dilatancy angle throughout, is expected to fail along a curved slip surface with
the shape of a log spiral. We often approximate these curved surfaces by planes or
cylinders (straight lines or circles in cross section). Both are strictly applicable only
if c 0, for then Eq. (4.31) reduces to that of a circle, and a straight line is nothing
more than a circle with infinite radius; but the approximation is sometimes justified
when the maximum value (umax u0) of u u0 in Eq. (4.31) is small. We analyze
Center of rotation
du
dr
Figure 4-22
rdu
Slip surface geometry.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 138
PB1 PA1 soils that shear under undrained conditions using total
l
2 dg PB2 PA2 stresses and c 0. We often do that for clays. The use of
B2
a slip surface with circular cross section then follows
A1
directly from Eq. (4.31).
c Slip Surface Direction*

c
90 c
The Mohr circle can be used to determine the two direc-
de3 45 c/2 de1
tions of the two potential slip surfaces at a point within a
P 45 c/2 de
soil mass. The direction of the slip surface at a point is the
c direction of its tangent at the point. The pole method can
be used to do this graphically.
Assume that the major principal strain increment is
B1 A2 vertical, therefore normal to the horizontal plane. In such
case, the pole P is located at the leftmost point of the cir-
cle, as shown in Fig. 4-23. We start by drawing lines PA1
Figure 4-23
Determination of direction of slip planes for a and PA2 from P through the points A1 and A2 correspon-
soil element. ding to zero normal strain. These lines are parallel to the
planes that are perpendicular to the directions of zero nor-
mal strain (that is, the slip surface directions). So lines PB1 and PB2, normal to PA1
and PA2, respectively, are the directions of the two potential slip surfaces through
the point under consideration.
From the geometry of the Mohr circle of Fig. 4-23, it can be shown that the
slip surfaces make an angle 45° c/2 with the direction of the minor principal
strain (which is horizontal in this case). If c f (a common assumption, even if
tacitly made, in many geotechnical analyses), then the slip surface directions corre-
spond to the directions of the planes where s and t satisfy the Mohr-Coulomb fail-
ure criterion, which are at 45° f/2 to the horizontal. Even when c differs from
f, so long as the difference is not large, the direction of the slip surfaces can still
be approximated as being 45° f/2 with respect to the horizontal. But, in dilative
sands (sands that exist in a dense state and/or in a state of low effective confining
stresses), the deviations from reality from assuming c f in analyses may be
substantial and should not be ignored. This point will be illustrated in Chapter 10
in the context of bearing capacity analyses of footings.
The Hoek-Brown Failure Criterion for Rocks

Although the Mohr-Coulomb failure criterion is also used for rocks, the Hoek and
Brown (1980, 1988) criterion is more practical, particularly for fractured rock, as
its parameters have been related to quantities usually measured in site investiga-
tions in rock. The criterion is written as
s¿3
s¿1 s¿3 qu m s (4.32)
B qu
where qu uniaxial, unconfined compressive strength of the rock; m and s

model parameters that depend on the degree of fracturing of the rock. The parame-
saL00581_ch04.qxd 9/14/06 4:54 PM Page 139
ters m and s can be expressed in terms of the type of rock, degree of

weathering, and frequency of discontinuities. This is covered in
Chapter 7.
4.4 At-Rest and Active and Passive g, c, f

Rankine States z
(a)
At-Rest State
Let us consider a semi-infinite mass of homogeneous, isotropic soil
[Fig. 4-24(a)]. The word homogeneous means the soil has the same Infinite plane
physical properties (for example, unit weight, hydraulic conductivity,
shear modulus, and shear strength) at every point, and the word (b)
isotropic means the directional properties at any point of the soil mass
are the same in every direction. This soil deposit has a free surface
that is level (horizontal); the soil deposit extends to infinity downward
from the free surface and in the directions parallel to the free surface.
Vertical and lateral stresses within this deposit are principal stresses,
as the shear stresses in horizontal and vertical planes are clearly zero. (c)
This soil deposit idealization of Fig. 4-24(a) allows us to obtain Extension
45 f/2
very important results. If we assume that the soil deposit has not
been disturbed after formation, the soil is said to be in a state of rest.
The ratio K0 sh /sv of the lateral to the vertical effective stress at
any point in a soil mass in a state of rest is referred to as the coeffi-
cient of lateral earth pressure at rest. Note that K0 is a ratio of effec- Active
tive, not total stresses. It is a very important quantity that appears (d)
often in geotechnical design. For a purely frictional soil, with Compression
strength parameters c 0 and nonzero f, Jaky (1944) established an 45 f/2
empirical relationship between K0 and f, given as9
K 0 1 sin f (4.33)
An alternative expression for K0 can also be obtained from elas- (e)
Passive
ticity,10 in which case K0 is related to the Poisson’s ratio of the soil
through Figure 4-24
(a) A semi-infinite soil mass; (b) a
n semi-infinite vertical plate separating
K0 (4.34)
1n two halves of the semi-infinite soil
mass; (c) vertical plate with soil
Equation (4.33) applies only to normally consolidated soils (that completely removed from the left side;
is, soils that have never experienced a larger vertical effective stress (d) soil in the active Rankine state, as
than the stress they currently experience). In a normally consolidated illustrated by drawing the two families
of possible slip surfaces at angles of
9 For clays, some prefer to use 0.95 instead of 1 in Eq. (4.33).
(45° f/2) with the horizontal; (e)
10 To derive this equation, we can write the equation for the lateral strain in terms of its soil in the passive Rankine state, as
coaxial normal stress (that is, the lateral effective stress) and the two transverse normal illustrated by drawing the two families
stresses (the vertical and again the lateral effective stress, but now in the other lateral of possible slip surfaces at angles of
direction) and make the lateral strain equal to zero.
(45° f/2) with the horizontal.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 140
soil, if the vertical effective stress is increased by an amount dsv, the lateral effec-
tive stress increases by an amount dsh K0 dsv in order to keep the ratio of sh to
sv constant and equal to K0. To consider the effects of stress history, assume that
the vertical effective stress is then reduced by dsv down to the original value; if that
is done, the lateral effective stress does not go back to its previous value, retaining a
considerable fraction of the increase dsh it experienced when sv was increased. It
follows that soils that have experienced a larger vertical effective stress previously,
referred to as overconsolidated soils, have higher K0 than normally consolidated
soils. Physically, the reason lateral stresses get locked in is the change in the fabric
and density of the soil required to accommodate the increase in vertical effective
stress, which is to some extent inelastic in nature and thus irrecoverable.
Brooker and Ireland (1965) investigated the effects of stress history on K0,
arriving at the following equation:
K 0 K 0, NC 2OCR 11 sin f2 2OCR (4.35)

where OCR overconsolidation ratio, defined as
s¿vp
OCR (4.36)
s¿v
where svp preconsolidation pressure, which is the maximum vertical effective
stress ever experienced by the soil element and sv is simply the current vertical
effective stress.
Rankine States
Level Ground In order to investigate the behavior of soil deposits when sub-
jected to relatively large strains, it will help us now to go through an imaginary
exercise. Let us consider that we could insert a smooth, infinite, infinitesimally
thin plane vertically into the soil deposit without disturbing the soil [Fig. 4-24(b)],
thus keeping vertical shear stresses equal to zero. Let us consider further that we
could remove all the soil from one side of the plane [Fig. 4-24(c)] so that we could
now either push or pull on the plane in the horizontal direction, thereby either
causing the soil on the other side of the plane to compress or extend in the hori-
zontal direction.
When we pull the vertical plane horizontally, it allows the soil to expand in
the horizontal direction, which leads to a drop in the horizontal effective stress sh .
If we continue to pull on the plane, sh continues to drop, while sv remains
unchanged. This process is illustrated by the Mohr circles in Fig. 4-25, all with the
same major principal stress s1 sv but a decreasing minor principal stress s3
sh . This process, by which sh decreases while sv remains unchanged, cannot go
on indefinitely; it in fact comes to an end when sv /sh becomes equal to the flow
number N s1 /s3 (1 sin f)/(1 sin f), at which point the Mohr circle
touches the Mohr-Coulomb failure envelope (assumed here with c 0), the failure
criterion is satisfied, and slip surfaces can potentially form anywhere in the soil
mass. This state in which the whole soil mass is in a state of incipient collapse is
known as an active Rankine state: active because the self-weight of the soil con-
saL00581_ch04.qxd 9/14/06 4:54 PM Page 141
t Mohr circle Figure 4-25

45 f/2 at rest Sequence of Mohr circles illustrating the horizontal
f Mohr circle unloading of a soil mass until the active Rankine state
45 f/2 for passive is reached and the horizontal compression of a soil
Rankine
state
mass until the passive Rankine state is reached.
s3
shA
sh sv s1
shP
PA PP s
Mohr circle
for active
45 f/2
Rankine state
f
45 f/2
tributes or is active in bringing it about; Rankine because it was first identified by

Lord Rankine in the 19th century (Rankine 1857, Cook 1950).
The coefficient KA of active earth pressure is the ratio of the lateral to the ver-
tical effective stress in a soil mass in an active state, given by
s¿hA s¿3 1 1 sin f
KA (4.37)
s¿v s¿1 N 1 sin f
The vertical effective stress is the major principal effective stress in the active
case, which means the leftmost point of the circle is the pole PA (Fig. 4-25). The
direction of potential slip surfaces is the direction of a line through the pole PA
passing through the point of tangency of the Mohr circle with the failure envelope.
As shown in Fig. 4-25, this direction makes an angle of
(45° f/2) with the hor-
izontal. Figure 4-24(d) shows the two families of potential slip surfaces associated
with the active Rankine state.
While the active state provides a lower bound to the value of the lateral earth
pressure coefficient K, the passive Rankine state is on the other extreme, capping
the possible values of K. The passive state is obtained by pushing the vertical plane
of Fig. 4-24(c) toward the soil mass, which increases sh until sh /sv becomes equal
to N. The progression from the at-rest to the passive Rankine state can be visual-
ized through Mohr circles as shown in Fig. 4-25. The coefficient KP of passive
earth pressure is the ratio of the lateral to the vertical effective stress in a soil mass
in a passive state, given by
s¿hP s¿1 1 sin f
KP N (4.38)
s¿v s¿3 1 sin f
Recognizing that the lateral effective stress is the major principal stress, we
can find the pole PP (see Fig. 4-25). Connecting the pole to the points of tangency
of the Mohr circle for the passive state and the strength envelope, we find that the
directions of the potential slip surfaces in the passive Rankine state make angles of
saL00581_ch04.qxd 9/14/06 4:54 PM Page 142
Figure 4-26 t f
Mohr circles representing the Passive
at-rest state and the active and
passive Rankine states. At rest
Active
s
shA KAs
v
sh 0 K0sv
sv gz
f
shP KPsv

(45° f/2) with the horizontal (Fig. 4-25). Figure 4-24(e) shows the directions
of the potential slip surfaces in the soil mass. Figure 4-26 shows the Mohr circles
for the three states we have examined and how the corresponding lateral effective
stresses relate to the vertical effective stress through the corresponding coefficient
of lateral earth pressure. Note that K0, KP, and KA are ratios of effective, not total,
stresses.
Active and passive earth pressure analysis can be easily extended to c-f mate-
rials (a good reference for that is Terzaghi 1943), but its usefulness in practice is
rather limited. We will therefore not discuss this extension in this text.
EXAMPLE 4-6
For a normally consolidated cohesionless soil with f 30°, calculate the at-rest, active, and
passive values of the earth pressure coefficient K.
■ Solution
The coefficient K0 of lateral earth pressure at rest is calculated using Eq. (4.33):
K 0 1 sin f 1 sin 30° 0.5
The active and passive lateral earth pressure coefficients are calculated using Eqs.
(4.37) and (4.38):
1 sin 30° 1
KA
1 sin 30° 3
1 sin 30°
KP 3
1 sin 30°
Note how the ratio of passive to active pressures is 9 for a 30° friction angle.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 143
ag Figure 4-27
1 Active Rankine state in sloping ground: (a) a prism with
unit width with base at depth z in a soil deposit with
sloping surface; (b) Mohr circle corresponding to active
z
state at base of prism.
T gz cos ag s gz cos2 ag
ag
t gz sin ag cos ag
(a)
f
t
gz cos2 ag
F A ag
E
P gz cos ag sin ag
0
C s
D
B ag
G
f
(b)
Sloping Ground* Consider a soil mass sloping at an angle ag with respect to the
horizontal [Fig. 4-27(a)]. For simplicity, consider the soil mass to be free of water.
Focusing on the prism of unit width shown in Fig. 4-27(a), we see that the weight
of the prism per unit length into the plane of the figure is equal to gz. The weight
of the prism is the only reason its base is subjected to a normal and a shear stress.
By projecting the weight normally and tangentially to the base of the prism and
dividing the component forces by the area of the base (which is 1/cos ag per unit
length of prism normal to the plane of the figure), we obtain the following expres-
sions for the normal and shear stresses at the base of the prism:
s¿ gz cos2 a g (4.39)
t gz sin a g cos a g (4.40)
Note that the effective traction T (the resultant of the normal and shear effec-
tive stresses) on the base of the prism of Fig. 4-27(a) is vertical and given by
T¿ gz cos a g (4.41)
It is important to understand that T is vertical but is not sv (that is, it is not
the vertical effective stress on the horizontal plane) and does not act on the hori-
zontal plane.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 144
We can plot the state of stress defined by Eqs, (4.39) and (4.40) in a Mohr dia-
gram as point A, as illustrated in Fig. 4-27(b). Note that the distance OA from the
origin of the s-t space to point A is equal to gz cos ag (that is, T). When OA is
projected onto the s- and t-axes, we get Eqs. (4.39) and (4.40). The Mohr circle
corresponding to the active Rankine state is drawn in the figure going through A.
Because point A defines the stresses on a plane making an angle ag with the hori-
zontal, the pole P lies on the intersection of the Mohr circle with the OA line.
Now let us say we wish to determine the resultant stress on a vertical plane.
Drawing a vertical line through the pole, we obtain point D as the intersection of
this line with the Mohr circle. The resultant stress (traction) TA (which, note, is not
the horizontal effective stress and is not even horizontal) is thus given by the line
OD; that is,
T¿ OA
T¿A OD
Let us now find the magnitude of the ratio of OD to OA and then its normal
and tangential components acting on the vertical plane. We do that by analyzing
the geometry of the Mohr diagram of Fig. 4-27(b). We first note that, because CE
is perpendicular to OA,
EP EA
We can now write
T¿A OD OP OE EP
(4.42)
T¿ OA OA OE EP
But OE is simply
OE OC cos a g (4.43)
To find EP, we note that AEC is a right triangle and that CA CF. We can
now write
EC OC sin a g
CA CF OC sin f
So we have
CA2 CE 2 EA2
which leads to
EA2 CA2 CE 2 OC 2 sin2 f OC 2 sin2 a g
and that, in turn, leads to
EA OC2sin2 f sin2 a g OC2cos2 a g cos2 f EP (4.44)

saL00581_ch04.qxd 9/14/06 4:54 PM Page 145
Taking Eqs. (4.43) and (4.44) into Eq. (4.42), we obtain an expression for the
active Rankine earth-pressure coefficient for sloping ground:
T¿A cos ag 2cos2 ag cos2 f

KA (4.45)
T¿ cos ag 2cos2 ag cos2 f
The horizontal stress is obtained by projecting the resultant stress TA acting on
the vertical plane onto the horizontal direction:
s¿hA K AT¿cos a g K A gz cos2 a g (4.46)
while the shear stress is obtained by projecting TA on the vertical direction:
t¿vA K AT¿ sin a g K A gz sin a g cos a g (4.47)
The directions of slip planes can be determined in much the same way as we
found directions for level ground, with lines drawn from P through F and G yield-
ing these directions.
4.5 Main Types of Soils Laboratory

Tests for Strength and
Stiffness Determination
Role of Stiffness and Shear Strength Determination
In the preceding sections, we examined stresses, strains, and the relationship
between strains and stresses in soils, particularly shearing (or failure). We did so
because they are an integral part of the calculations we do in the analysis of foun-
dations, slope, retaining structures, and other geotechnical systems. In this section,
we examine the important issue of how to measure stress-strain properties properly
and, in particular, the shear strength of soils. This issue is fundamental both for
work done in practice and for research on soil load response.
Stress Paths in s-t Space

It is not practical to trace the loading of a soil element or soil sample in a labora-
tory test by using Mohr circles because an infinite number of Mohr circles would
be required and no discernible representation of the loading process would result.
A diagram that achieves the same purpose is based on plotting only the highest
points of the Mohr circles, producing a continuous line. The highest point of the
Mohr circle has coordinates s and t (as shown in Fig. 4-28) given by
s 12 1s1 s3 2 (4.48)
t 12 1s1 s3 2 (4.49)
where s corresponds to the center of the Mohr circle, and t is equal to the radius of
the Mohr circle and is thus a measure of shear stress.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 146
Figure 4-28 t
Definition of stress variables s (s, t)
and t.
1

2 (s1 s3)
(s3, 0) (s1, 0) s
1 (s s )

2 1 3
A stress path is a plot in s-t space of the progression of (s, t) points represent-
ing the loading process for a soil element or laboratory soil sample. Figure 4-29
shows the four possible general directions a stress path can take from a point (s0,
t0) depending on whether only s1 or only s3 changes, and whether the change is
an increase or a decrease.
Referring to Eqs. (4.48) and (4.49), we can write the differentials of s and t as
ds 12 1ds1 ds3 2 (4.50)
dt 12 1ds1 ds3 2 (4.51)

When only s1 changes, ds3 0 and
ds 12ds1
dt 12ds1
s1 unchanged
s1
s3
s3 unchanged
(s0, t0)
s1 unchanged
s1
s3
s3 unchanged
Figure 4-29
Possible directions for s-t paths from a starting point (s0, t0). s
saL00581_ch04.qxd 9/14/06 4:54 PM Page 147
When only s3 changes, ds1 0 and
ds 12ds3
dt 12ds3
It follows that all possible stress paths when only one of s1 or s3 changes are
straight lines at 45° with the horizontal, with dt/ds 1 when s1 alone changes and
dt/ds 1 when s3 alone changes. Naturally, other directions are possible if s1
and s3 are allowed to change simultaneously.
Stress plots can also be done in terms of effective stresses:
s¿ 12 1s¿1 s¿3 2 s u (4.52)
t¿ 12 1s¿1 s¿3 2 t (4.53)

where u pore pressure (see Chapter 3).
Effective stress paths clearly coincide with total stress paths if loading is
applied under drained conditions. Under undrained conditions (u 0), effective
stress paths do not coincide with total stress paths and are not linear.
Stress Paths in p-q Space

The stress variables s and t are two-dimensional variables that do not capture the
effects of s2. The mechanical response of soil is more rigorously expressed in
terms of p and q, which are defined as follows in terms of the principal stresses:
p 13 1s1 s2 s3 2 (4.54)
2 1s1 s3 2 2 1s1 s2 2 2 1s2 s3 2 2

1
q (4.55)
22
Those readers with a background in mechanics will recognize p and q as being
related to the octahedral normal and shear stresses sm soct and toct through
p sm soct (4.56)
22
q toct (4.57)
3
where the octahedral normal stress is also known as the mean stress sm. Thus, p
and q are normal and shear stresses that are representative of the three-dimensional
stress state at a point.
Effective stress versions p and q of p and q can also be defined:
p¿ p u (4.58)
q¿ q (4.59)
saL00581_ch04.qxd 9/14/06 4:54 PM Page 148
Main Laboratory Tests

The laboratory soil tests of greatest importance for what we cover in this text are
the one-dimensional compression test (more routinely called the consolidation
test11), the direct shear test, the triaxial compression test, and the unconfined com-
pression test. (See Figs. 4-30 to 4-32.) The consolidation test is a one-stage test,
Q
WT
Porous
stones Sample
Figure 4-30
Consolidometer: Consolidation
(a) photo, ring
(b) schematic diagram. (a) (b)
Slip surface N
T0
T0 Soil sample
Before application
of shear force T
T0
T0 Soil sample
Figure 4-31
Direct shear machine: After application
(a) photo, of shear force T
(b) schematic diagram. (a) (b)
11Note that soil consolidation can be accomplished in many different test setups; however, we will retain
the more traditional name for this test in the remainder of the book.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 149
Figure 4-32
Triaxial test system: (a) photo,
(b) schematic diagram.
Axial force
Top platen
Drainage
Soil Cell
Rubber
sample pressure s3
membrane
Bottom platen
Drainage
(a) (b)
while both the direct shear test and the triaxial compression test have a consolida-
tion stage (in which the sample is taken to a desired stress state) and a shearing
stage (in which the shear stress level is increased within the sample, usually until
very large shear strains, of the order of 20%, develop). The unconfined compres-
sion test is a test that can be done only on soil samples that don’t fall apart when
unconfined and don’t have major defects (fissures, seams of different materials cut-
ting across the soil, and the like). Many clays can be tested in this manner, as we
will see in Chapter 6. The sample is not subjected to a confining stress, that is, the
total stress s3 0. It has only a shearing stage, in which an axial load is applied
on the sample until failure.
The consolidation (one-dimensional compression) test is performed by apply-
ing a vertical load on a cylindrical sample of soil restrained laterally by a steel ring
(Fig. 4-30) and measuring the resulting vertical deflection. If we recall our discus-
sion of the at-rest state and of the coefficient of lateral earth pressure at rest as a
ratio of sh to sv for a normally consolidated, level soil deposit, it is easy to under-
stand why the stress path imposed by this test is as shown in Fig. 4-33(b). The
sample is kept under water to ensure full saturation. Although there is more than
one size of consolidometers, the most typical sample size is 70 mm (diameter) by
25 mm (thickness). The consolidation test is particularly useful in the study of
clays, and we will therefore rediscuss it in Chapter 6.
The direct shear test is performed on a soil sample that is also restrained hori-
zontally by a split box constituted of a top and a bottom part that can move hori-
zontally with respect to each other (Fig. 4-31). The sample is first loaded vertically
saL00581_ch04.qxd 9/14/06 4:54 PM Page 150
t
C
1
t 1
C
t K0 line A B s
t
K0 line B D
A s AB Isotropic consolidation
AB 1D consolidation BC Shearing in compression
s s BC Shearing BD Shearing in extension
(a) (b) (c) (d)
Figure 4-33
Stress paths for (a) isotropic compression, (b) one-dimensional consolidation, (c) direct shear test, and (d) triaxial
compression and extension tests.
without the possibility of lateral expansion (much as in a one-dimensional com-

pression test); then, a horizontal force is applied on the upper part of the box, caus-
ing the sample to shear in the horizontal direction. The stress path for this test is
shown in Fig. 4-33(c). The consolidation portion (AB) of the stress path is just like
the stress path of a one-dimensional compression test. The shearing is a vertical
line BC, indicating that the shear stress added to the sample horizontally does not
change s, just t. In truth, s may increase slightly because of one of the shortcom-
ings of the direct shear test, which is that some passive pressures develop as some
of the push imposed in the upper part of the sample ends up in part as a push
against the lower part of the box on the other side of the sample. This also means
that the shear strength of the soil is not the only resistance to failure; this passive
resistance also contributes to the resistance to shearing of the sample. Another
shortcoming of this test is that it forces shearing to develop horizontally instead of
allowing it to develop naturally along the optimal shearing plane. Nonetheless, the
direct shear test is used in practice often and does provide a reasonable estimate of
shear strength.
In a triaxial (TX) test, a cylindrical sample, typically either 70 or 35 mm in
diameter and at least twice that in height, is wrapped in a close-fitting, cylindrical,
impervious membrane and placed between two platens. The bottom platen is
attached to the base of the triaxial chamber and is stationary. The top platen is
attached to a piston to which an axial force may be applied. During the test, the
sample, platens, and piston are located within a plexiglass or aluminum cylinder
that is sealed at the top and bottom. So, in the course of a test, the sample is first
subjected to an all-around stress by increasing the air pressure around the sample
and then to a change in the axial stress. The first stage is called consolidation
(although the soil may or may not undergo volume change); the second, shearing.
The impervious membrane allows transmission of the applied pressure to the sam-
ple as a total stress. Drainage may be either allowed or prevented during either or
saL00581_ch04.qxd 9/14/06 4:54 PM Page 151
both the consolidation and shearing stages. The effective stress state at the end of
the consolidation stage will be different from the initial state only if drainage is
allowed. Tests in which drainage is allowed during the consolidation stage are
referred to as consolidated or C tests; when drainage is prevented, they are known
as unconsolidated or U tests. Tests in which drainage is allowed during the shear-
ing stage are known as drained (D); otherwise, tests are referred to as undrained
(U). Table 4-1 summarizes the drainage conditions for the three main types of tri-
axial tests.
Table 4-2 shows the evolution of the effective stress state acting on a triaxial
test sample tested under drained conditions. Samples are usually saturated with
water. One of the advantages of doing so is that sample volume change can be eas-
ily measured by measuring the volume of water either expelled from within the
sample or sucked in by the sample. When water is sucked in, the void ratio of the
soil increases, weakening the soil. When water is expelled, the void ratio drops,
and the sample gets stronger. Note that there is no generation of pore pressures due
to consolidation or shearing during drained tests, as the drainage lines are open
during both consolidation and shearing. However, we usually apply a pore pressure
ub (called back pressure) within the sample to dissolve air bubbles into the pore
water, ensuring sample saturation. So the total stresses acting on the sample are
larger than the effective stresses of Table 4-2 by a magnitude ub. CD tests are com-
monly used for sands.
In CU tests, it is possible and desirable to measure the pore pressures during
shearing. The effective stresses within the sample can then be calculated through-
out shearing, and effective stress paths in terms of either (s, t) or (p, q) can be
plotted. CU tests are commonly used for clays. Sands are sometimes tested this
way as well either to enhance understanding of the behavior of the soil or to quan-
tify its response to loads that are applied very fast, as during earthquakes, when the
rate of loading is much faster than the rate of dissipation of excess pore pressure
generated by the loading.
Table 4-1 Types of triaxial tests

TX tests
Drainage during CD CU UU
Consolidation Drained Drained Undrained
Shearing Drained Undrained Undrained
Table 4-2 Evolution of applied effective stresses in TX-CD Tests

Stress Initial After consolidation During shearing
s3 Zero sc sc
s1 Zero sc sc sd
sc sc ub; sc cell pressure; ub back pressure; sd axial force applied
by the piston /sample cross-sectional area.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 152
In a triaxial test, s2 s3; it follows that p, q, p, and q take the following def-

initions:
p sm 13 1s1 2s3 2 (4.60)
p¿ s¿m 13 1s¿1 2s¿3 2 (4.61)

q q¿ s1 s3 (4.62)
4.6 Stresses Resulting from the Most

Common Boundary-Value Problems
Elastic Stress-Strain Relationship and Elastic
Boundary-Value Problems
Stresses appear within a soil mass as a result of loads applied on its boundaries or
loads applied within it, including the soil’s own self-weight. Applied loads can be
point loads, line loads, or distributed loads. The stresses decrease in magnitude with
distance from the applied loads. Solutions from elasticity theory are available to cal-
culate the stress components in two perpendicular planes (for two-dimensional
problems) or three mutually perpendicular planes (for three-dimensional problems)
at a point within a soil mass. For level-ground soil deposits, these directions are
usually the horizontal and vertical directions.
Elasticity solutions are applicable so long as the deformations within the soil
mass can be assumed to be small (which means the soil mass is far from a state of
collapse or sliding). They may be used for well-designed foundations and other
geotechnical structures operating under working loads. The linear elastic stress-
strain relationship of materials is determined by a pair of constants, referred to as
the elastic pair. Solutions of interest in practice are usually formulated in terms of
any two of the following four elastic constants: Young’s modulus E, shear modulus
G, Poisson’s ratio n, and bulk modulus K.
Young’s modulus and Poisson’s ratio are best understood by considering a
stress state in which the only nonzero component is a normal stress acting in some
direction (this is the case of an axially loaded bar, for example, which forms the
basis for many common tests, including uniaxial extension tests of steel bars and
the unconfined compression tests of concrete or soils). Young’s modulus E is the
ratio of this normal stress to the collinear normal strain. Poisson’s ratio n is the ratio
of the strain appearing normal to the direction of the nonzero stress to the normal
strain in the direction of the stress. The shear modulus G is the ratio of a shear
stress to the corresponding shear strain. The bulk modulus K is the ratio of a mean
stress to the corresponding volumetric strain. Table 4-3 shows the relationship
between E, n, G, and K.
An approach used to solve elastic boundary-value problems is usually the
determination of strains from stresses, and integration of the strains, subject to com-
patibility, to yield the displacements. There are other approaches based on the con-
saL00581_ch04.qxd 9/14/06 4:54 PM Page 153
Table 4-3 Relationship between the four most common

elastic constants
Elastic pair
Elastic constant E, N K, G
9KG
Young’s modulus E E
3K G
3K 2G
Poisson’s ratio n n
6K 2G
E
Shear modulus G
2 11 n2
G
E
Bulk modulus K
3 11 2n2
K
E 11 n2 3K 4G
Constrained modulus M
11 n2 11 2n2 3
cept of stress or strain potentials. Some of the solutions can be obtained by direct
integration of the displacements from the Boussinesq solution (discussed later in
this section). Whatever the approach, the solutions are typically somewhat involved
and are outside the scope of this text; they can be found in elasticity texts such as
Timoshenko and Goodier (1970), Fung (1993), or Saada (1993). We will present
only the results that will be of use to engineers solving foundation problems.
Elasticity solutions are given in terms of stresses, strains, and displacements.
The displacements are obtained from the strains by integration. A few displace-
ment solutions are useful in foundation engineering because they form the basis
for calculation of the settlement of shallow foundations; they are presented in
Chapter 9. We restrict ourselves in this chapter to stresses created within soil
masses by boundary loads.
A point anywhere in the soil can be expressed with respect to the origin by
using Cartesian coordinates (x1, x2, x3), cylindrical coordinates (radius r, depth z,
and angle u), or spherical coordinates (radius r and angles u and w). The notation
x, y, and z can be used interchangeably with x1, x2, x3, as convenient. The choice of
which coordinate system to use depends, for any problem, on which system leads
to an easier solution of the problem or the most compact equations. We will use sv
for vertical stress whenever possible, instead of the sz that would typically be used
in elasticity theory, in order to maintain consistency with normal soil mechanics
practice and the other chapters. The same is not done for vertical strain, which is
denoted by ez, because ev is reserved for volumetric strain.12
12Notation can be problematic in a discipline that deals with so many concepts and lies at the crossroads
of many different branches of science. We will keep notation as simple, logical, and consistent as
possible, but keeping track of it will still require some effort on the reader’s part.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 154
Loadings can be applied to the surface of an elastic soil mass in two ways: by
application of point or distributed loads or by application (or imposition) of dis-
placements. When loads are applied, that implies ideal flexibility. There is no inter-
nal rigidity to the loading. When displacements are imposed, there is rigidity
implied. For example, we might push a rigid cylinder vertically down into a soil
mass. This means every point of the cylinder in contact with the soil will settle by
the same amount. This is what happens when a rigid footing resting on soil is
loaded. We will now discuss the stresses that are generated inside a soil mass by
various types of loads.
Vertical Point Load on the Boundary of a Semi-Infinite,

Elastic Soil Mass (Boussinesq’s Problem)
Figure 4-34 shows a vertical point load Q applied at a point on the boundary of a
semi-infinite mass with Young’s modulus E and Poisson’s ratio n. The point of
application of Q is taken as the origin of the system of coordinates. The nonvanish-
ing stress components are the normal stresses sr, su, sv sz and the shear stress
srz; the other shear stresses are zero as a result of the symmetry of the problem.
The equations are as follows:
11 2n2 1r 2 z 2 2 1>2
c d
Q 3r 2z
sr
2p 1r 2 z 2 2 1r 2 z 2 2 3>2 1r 2 z 2 2 1>2 z
(4.63)
11 2n2Q 1r2 z2 2 1>2

c d
z
su
2p 1r z 2 1r z 2 1r z2 2 1>2
(4.64)
2 2 2 2 1>2
z 2
sv
Q
E, n
z sz sv sr
Cross srz
section r su
sr
su Expanded view
(normal to plane)
u sz sv (Normal to plane)
Q r sr
Plan
view su
All points of a circle
Figure 4-34 at depth z with radius
r are subject to the
Vertical point load on boundary of soil mass
same stresses (axial
(Boussinesq’s problem). symmetry).
saL00581_ch04.qxd 9/14/06 4:54 PM Page 155
3Qz3
sz sv
2p 1r2 z2 2 5>2
(4.65)
3Qrz2
srz
2p 1r2 z2 2 5>2
(4.66)
Vertical Point Load Within a Semi-Infinite, Elastic Soil

Mass (Kelvin’s Problem)
Figure 4-35 shows a vertical point load Q applied at a point within a semi-infinite
mass with Young’s modulus E and Poisson’s ratio n. The point of application of Q
is taken as the origin of the system of coordinates (radius r, depth z, and angle u in
cylindrical representation or x1, x2, x3 z in Cartesian coordinates). The nonvan-
ishing stress components are the normal stresses sr, su, sv sz and the shear
stress srz; the other shear stresses are zero as a result of the symmetry of the prob-
lem. The equations are as follows:
11 2n2z
c 2 d
Q 3r 2z
sr 2
8p 11 n2 1r z 2 1r z2 2 3>2
2 5>2
(4.67)
11 2n2z
c 2 d
Q
su
8p 11 n2 1r z2 2 3>2
(4.68)
E, v
sz sv
Cross z
section P
r sr
su
u
Q s r , sz
Plan
view su
All the points on a

circle with radius r
are subject to the Figure 4-35
same stresses Vertical point load within soil mass
(axial symmetry). (Kelvin’s problem).
saL00581_ch04.qxd 9/14/06 4:54 PM Page 156
11 2n2z
c 2 d
Q 3z 3
sv
8p 11 n2 1r z 2 2 3>2 1r 2 z 2 2 5>2
(4.69)
11 2n2r
c 2 d
Q 3rz2
srz
8p 11 n2 1r z2 2 3>2 1r2 z2 2 5>2
(4.70)
Uniform Pressure Distributed over a Circular Area on

the Boundary of a Semi-Infinite, Elastic Soil Mass
Consider the circular loaded area of diameter B and radius b of Fig. 4-36. If the
distributed load is expressed as qb and the center of the circle is taken as the origin
of the system of coordinates, the nonvanishing stress components along the axis of
the loaded area are as follows:
qb c 1 d
z3
sv `
1b 2 z 2 2 3>2
(4.71)
r0
2 11 n2z
qb e 1 2n c 2 d f
3
1 z
sr ` su `
1b z 2 2 1>2 1b 2 z 2 2 1>2
(4.72)
r0 r0 2
Uniform Pressure Distributed over a Rectangular Area

on the Boundary of a Semi-Infinite, Elastic Soil Mass
Figure 4-37 shows a rectangular loaded area, with plan dimensions B and L. If the
magnitude of the distributed load is qb, the vertical stress s33 sz at depth x3 z
under any corner of the rectangle is given as
qb 2mn2C 1 1 C 1
c tan1 a bd
2mn2C 1
s33` sv` (4.73)
corner corner 4p C 1 C 2 C1 C1 C2
qb x
qb
B r L B
y y
sv r 0
sv
P
sr r 0 P
z
z
Figure 4-36 Figure 4-37
Circular loaded area on the boundary of a semi-infinite Rectangular loaded area on the boundary of a semi-
soil mass. infinite soil mass.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 157
where
B
m (4.74)
z
L
n (4.75)
z
C1 1 m2 n2 (4.76)
C 2 1mn2 2 (4.77)
The arctangent term in Eq. (4.73) must be a positive angle in radians; thus,
when C2 C1 and a calculator gives a negative angle, we must add p to that angle.
Using Eq. (4.73), we can calculate the vertical stress at any point with depth z
by using the principle of superposition. For example, the vertical stress at a depth z
below the center of a rectangular loaded area would be the sum of the stresses
under the common corner of the four rectangular loaded areas with half the width
and length of the original area that, considered together, constitute the loaded area.
By using a combination of positive and negative loaded areas, it is possible to cal-
culate the stress induced by a rectangle at a point not lying directly below the
loaded area, as shown in Example 4-7.
EXAMPLE 4-7
The load of a heavy building with plan dimensions 10 and 15 m is assumed to be uniformly
distributed and equal to 150 kPa. Calculate the vertical stress induced by this load at a point
P, which is 5 m from the smaller side and 5 m from the larger side of the building, as shown
in Fig. 4-38, at a depth of 5 m.
■ Solution
Figure 4-38 shows that we can do this calculation by considering the uniformly distributed
load qb applied over a large rectangle 1 (ACGP) with the corner lying directly above the
point where the stress is desired. This generates an excess stress that is now compensated
for by assuming qb applied across rectangles 2 (ABHP) and 3 (FDGP). We are not done
yet because we have now removed the load across the intersection of rectangles 2 and 3
twice. We must add this area, rectangle 4 (FEHP), back in by considering a new load qb
applied over it. We can use Eq. (4.73) directly.
For the large rectangle ACGP:
B 15
m 3
z 5
L 20
n 4
z 5
C 1 32 4 2 1 26
C 2 13 4 2 2 144
saL00581_ch04.qxd 9/14/06 4:54 PM Page 158
Figure 4-38 Building Cross section

(a) Rectangular loaded area for Example
4-7; (b) composition of rectangles used in
P 5m
calculation.
10 m Plan
view
P 5m
15 m 5m
(a)
C D G
E
B H
P
A F
(b)
qb 2 3 4226 1 26 2 3 4226
sv ` c tan1 a bd
corner 4p 26 144 26 26 144
150 2 3 4226 1 26 2 3 4 226

sv ` c tan1 a bd
corner 4p 26 144 26 26 144
30.748 tan1 11.0372 4

150

4p
We see that the arctangent argument is a negative number. If we take the arctangent of
1.037 in a calculator we get 0.804, a negative angle in radians, not the positive angle
producing the same arctangent, which is the one we need. We calculate that number by
adding p to 0.804, obtaining 2.338. This gives us
30.748 2.338 4
qb
sv `
corner 4p
0.246qb 0.246 150 36.9 kPa
Calculations for the other three rectangles and solution of the whole problem are
shown in Table 4-4. We find that the vertical stress at point P is only 2 kPa. The small stress
we calculated suggests that vertical stress dissipates very quickly as we move away from the
saL00581_ch04.qxd 9/14/06 4:54 PM Page 159
Table 4-4 Results for Example 4-7

Rectangle 1 (ACGP) 2 (ABHP) 3 (FDGP) 4 (FEHP) Final
m 3 1 1 1 —
n 4 4 3 1 —
C1 26 18 11 3 —
C2 144 16 9 1 —
sv (kPa) 36.9 30.6 30.5 26.3 2.1
building. In a problem at the end of the chapter, the stress at the same depth as in this exam-
ple under the center of the building is calculated.
Vertical Line Load on the Boundary of a Semi-Infinite,

Elastic Soil Mass
Consider a vertical line load Q (per unit length along the x-axis) applied on the
boundary of a semi-infinite soil mass as shown in Fig. 4-39. The problem is a
plane-strain problem (the normal strain in the x direction, the direction of the line,
is zero everywhere in the soil mass). The depth z is related to r and u through
z r cos u (4.78)
The distance y from the line load to the projection on the soil boundary of the
point where the stresses are calculated is
y r sin u (4.79)
It is best to use cylindrical coordinates x, r, and u, with the line load aligned
with the x-axis. The equations for the nonzero, normal stress components that may
Q Q
3D view Cross-sectional view

y y
u u
sr sr
sx su
su
sv
(a) (b)
Figure 4-39
Vertical line load on the boundary of semi-infinite soil mass.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 160
be useful in typical foundation engineering problems, obtained from integration of

the Boussinesq solution for a point load, are
2Q cos u 2Q z
sr (4.80)
p r p r2
2Qn cos u 2Qn z
sx (4.81)
p r p r2
2Q cos3 u 2Q z3
sz (4.82)
p r p r4
Given that the shear stresses sru and srz are zero, sr is a principal stress.
Uniform Pressure Distributed over an Infinitely

Long Strip on the Boundary of a Semi-Infinite,
Elastic Soil Mass
Figure 4-40 shows a strip load qb applied on the boundary of a semi-infinite soil
mass. The problem is a plane-strain problem (the normal strain in the direction of
the strip is zero). This problem can be solved by integrating the solution for the
vertical line load. In this case, the most compact representation of a point P is in
terms of the angles u1 and u2 and the two vertical lines that pass through the edges
of the loaded area. The nonvanishing stress components13 are
32 1u 2 u 1 2 1sin 2u 2 sin 2u 1 2 4
qb
s11 (4.83)
2p
1u 2 u 1 2
2nqb 2nqba
s22 (4.84)
p p
32 1u 2 u 1 2 1sin 2u 2 sin 2u 1 2 4
qb
s33 (4.85)
2p
3 1cos 2u 1 cos 2u 2 2 4
qb
s13 (4.86)
2p
where a is the angle formed by PP1 and PP2 (shown in Fig. 4-40), which is given by
a u2 u1 (4.87)
The strip load is one of the most important boundary loadings in soil mechan-
ics because the analysis of strip foundations serves as a basis for the core of shallow
foundation design. Recall from our earlier discussion that when a load (as opposed
to a displacement) is applied, as in this case, the results are strictly applicable to a
13 A reminder of this type of stress notation, discussed earlier in the chapter: sij is the stress component
in the direction j acting on the plane normal to direction xi; when i j, sij is a normal stress in that
direction; when i j, sij is a shear stress, for it acts in a given plane in a direction tangential to it.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 161
qb Figure 4-40
Uniform pressure over
B infinite strip on boundary of
semi-infinite soil mass.
P1 P2 x1
u2
u1
x2 a
P
x3
flexible foundation. We will explore some additional facts about this loading later in
the chapter and in the Conceptual Problems section for this chapter.
Rigid Strip and Rigid Cylinder on Boundary

of Semi-Infinite Elastic Soil Mass
A rigid strip (Fig. 4-41) with width B, loaded by a central load Q, expressed as
units of load per length settles uniformly when resting on a semi-infinite, elastic
soil mass. The stress distribution qb at the base of the strip is far from uniform, as
larger stresses near the edges of the load are needed to keep the settlement uni-
form, given the proximity to the free boundary. If we orient the x-axis along the
boundary and normal to the strip, with x 0 at the center of the strip, we obtain
the following equation for qb from elasticity theory:
Q
qb (4.88)
p214B2 x2
Note that qb is equal to
2Q
qb (4.89)
pB
at the center (x 0) and to infinity at the edges (x B/2). The average unit load
(stress) qb,avg is obviously equal to Q/B.
The solutions for a rigid cylinder (a cylindrical footing) with diameter B are very
similar to those for the rigid strip (a strip footing). The contact stress (distributed
x or r
Figure 4-41
Cross section of either a rigid strip or a
x for strip footing rigid cylindrical footing on boundary of
r for cylindrical footing w wz semi-infinite, elastic soil mass.
z0
saL00581_ch04.qxd 9/14/06 4:54 PM Page 162
load qb), expressed in terms of cylindrical coordinates with origin at the center of
the circular cross section of the cylinder, is given by
2Q
qb (4.90)
pB2B2 4r 2
Approximate Stress Distribution Based

on 2:1 Vertical Stress Dissipation
The vertical stress sz under a strip or rectangular load can be approximated by
assuming that the stresses dissipate downward following a 2:1 (V:H) slope as
shown in Fig. 4-42. If the strip load is denoted by qb (in units of stress, that is, load
per unit length per unit length) and the resultant load is denoted by Q (in units of
load per unit length), the resulting equation for a strip footings is
Q qb
sv (4.91)
Bz z
1
B
For a rectangular footing, we have a load Q (units of load) distributed over an
area BL, leading to qb (units of stress). The equation for the stress at depth z is
Q qb
sv
1B z2 1L z2
(4.92)
a1 b a1 b
z z
B L
Figure 4-42 qb
2:1 stress dissipation with
depth.
z
2
1
EXAMPLE 4-8
Consider a 2-m-wide strip load with uniform magnitude 150 kPa placed directly on the sur-
face of the elastic half-space. Calculate and plot the distribution of the vertical stress with
depth z right below the centerline (y 0 m) of the strip load using both elasticity theory
and the 2:1 vertical stress dissipation method.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 163
0 Figure 4-43
1 2:1 method Distribution of vertical stress with depth
along the centerline of the strip footing of
2
Example 4-8 from both elasticity theory
3 Elasticity and the 2:1 method.
theory
Depth (m)
4
5
6
7
8
9
10
0 50 100 150
sv (kPa)
■ Solution
The vertical stress from elasticity theory is given by Eq. (4.85), where s33 sv. For points
lying on the plane of symmetry of the load, u2 u1 u, and u is given, as a function of
depth (z x3), by
u tan1 a b
B>2
z
Thus,
14u 2 sin 2u 2
qb
sv
2p
e 4 tan1 a b 2 sin c 2 tan1 a bdf

150 2>2 2>2
sv
6.28 z z
23.9 e 2 tan1 a b sin c 2 tan1 a b d f

1 1
z z
For the 2:1 vertical stress dissipation method, we have
Q qb 150
sv
Bz z z
1 1
B 2
Figure 4-43 shows results from the application of the preceding equations for the 0- to
10-m depth range. We observe that the approximate 2:1 method underestimates the value of
the vertical stress within the depth range between 0 and 1.5B vertically below the load and
overestimates it for larger depths.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 164
Saint-Venant’s Principle
Saint-Venant’s principle is occasionally useful in design calculations. It states that,
if a point where we may want to calculate the stresses (or strains) is sufficiently
removed from the loading causing these stresses (or strains), any statically equiva-
lent load generates approximately the same stresses (or strains). Statically equiva-
lent loads are loads with the same resultant and same moment with respect to an
arbitrary point. What this principle states is that the details (that is, the geometry)
of the loading have a negligible effect on the calculated stresses so long as the dis-
tance from the point where the stresses are desired to the loading is much larger
than the scale of the load.
EXAMPLE 4-9
A uniform, 2-m-wide square 100-kPa load is placed directly on the surface of a semi-infinite,
elastic soil mass. For this load and for a statically equivalent point load, calculate the verti-
cal stress at points A* and B* located 2 m below points A and B, respectively (refer to Fig.
4-44). Point A is located at the center of the square load; point B is located 3 m away from
both sides of the square load (a distance 4.25 m from the closest corner of the footing along
its diagonal). Discuss the accuracy/applicability of the approximation using a statically
equivalent point load.
■ Solution
To calculate the vertical stress at point A*, we discretize the applied load as four squares,
each with sides 1 m long, sharing point A at one corner. This way, the induced stress will be
four times the stress calculated using Eq. (4.73):
qb 2 mn 2C 1 1 C 1
c tan1 a bd
2 mn2C 1
sv `
corner 4p C 1 C 2 C1 C1 C2
1 1 1 1
2 21.5 2 21.5
1 1.5
≥ tan1 ± ≤¥
100 2 2 2 2
sv ` 4sv `
A corner p 1 1.5 1
1.5 1.5
16 16
31.83 30.653 tan1 10.4262 4 31.8330.653 0.403 4 33.6 kPa 34 kPa
Now we calculate the stress increase due to an equivalent load Q, obtained by multi-
plying the stress qb by the area on which it acts:
Q qb area 100 4 400 kN
Equation (4.65) is used next:
3Qz 3 3 400 2 3
sv 47.7 kPa 48 kPa
2p 1r 2 z 2 2 5>2 2p 10 2 2 2 2 5>2
saL00581_ch04.qxd 9/14/06 4:54 PM Page 165
100 kPa Figure 4-44

The square load and the points where the
vertical stresses are to be calculated in
A B Example 4-9.
2m
Cross
section
A* B*
C F D
A 2m
J
H G
Plan 3m
view
2m 3m
B
E I
Note that this stress is much larger than the 34 kPa calculated earlier using the exact solu-
tion, indicating that use of Eq. (4.65) here is not appropriate because the point at which the
stress is to be calculated is too close to the load.
To calculate the stress at point B*, we use the same two approaches. Using Eq. (4.73)
for rectangles BECD, BDFI, BEHJ, and BIGJ, the desired stress is expressed by the super-
position principle of elasticity as:
sv ` sv ` sv ` sv ` sv `
B,EFGJ B,BECD B,BDFI B,BEHJ B,BIGJ
Considering that the effects on B of BDFI and BEHJ are the same:
sv ` sv ` 2sv ` sv `
B,GHCF B,BECD B,BDFI B,BIGJ
with
5 5 5 5
2 213.5 2 213.5
£ tan1 ° ¢ § 24 kPa
100 2 2 14.5 2 2
sv `
B,BECD 4p 52.6 13.5 25.6
3 5 3 5
2 29.5 2 29.5
£ tan1 ° ¢ § 22.7 kPa
100 2 2 10.5 2 2
sv `
B,BDFI 4p 23.6 9.5 4.6
3 3 3 3
2 25.5 2 25.5
£ tan1 ° ¢ § 21.6 kPa
100 2 2 6.5 2 2
sv `
B,BIGJ 4p 10.6 5.5 0.4
saL00581_ch04.qxd 9/14/06 4:54 PM Page 166
This gives
sv ` 24 2 22.7 21.6 0.2 kPa
B,GHCF
or, for practical purposes, no stress increase.

Using a statically equivalent point load to approximate the rectangular load:
3Qz 3 3 400 2 3
sz 0.2 kPa
2p 1r 2 z 2 2 5>2 2p 132 2 2 2 5>2
Clearly, for point B*, both methods produce the same result because point B* is suffi-
ciently removed from the load that its details (such as its geometry, that is, whether it is a
point load or a distributed load with some shape) become unimportant, and any statically
equivalent load will produce acceptable results.
4.7 Total and Effective Stress Analysis

As seen in both Chapter 3 and this chapter, when loads are applied on the bound-
aries of soil masses, total stress increments appear at every point (or element) of
the soil mass. The magnitude and direction of the total stress increments appearing
in an element of the soil mass depend on the magnitude and direction of the loads,
the geometry of the soil mass, and the distances from the element to the applied
loads. Each total stress increment has two components: a pore pressure increment
and an effective stress increment. The magnitudes of these two components depend
on the hydraulic conductivity of the soil: In free-draining soils, pore pressures do
not build up and the total stress increments are instantaneously converted into
effective stress increments. In soils with very small hydraulic conductivity, the
stress increments are initially balanced by equal pore pressure increments; over
time, the pore pressures dissipate as water flows from points with large to points
with small hydraulic heads, and effective stresses increase by the same amounts as
the pore pressures decrease until, after a sufficiently long time, the pore pressures
return to their initial values and the effective stress increments become equal to the
total stress increments.
A soil mechanics problem can be analyzed in terms of either total or effective
stresses. If an effective stress analysis is used, all the quantities must be expressed in
terms of effective stresses. For example, in stability problems, where a slip surface
will form when the shear strength of the soil is reached, the normal stress of interest
on the slip surface is an effective normal stress, and the shear strength parameters (c
and f) must have been determined based on effective stresses. Laboratory tests used
to determine these parameters must then be performed sufficiently slowly for pore
pressures not to appear (not an economical or practical proposition with slow-draining
soils) or, alternatively, must be performed with equipment that can measure pore
pressures that develop within the sample. On the other hand, if total stress analysis is
used, then the total stresses resulting from the loading of the soil mass must be used
in conjunction with a total-stress-based, undrained shear strength.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 167
In sands, effective stress analysis is almost always used because sands are
free-draining, and induced pore pressures are equal to zero (the exception are seis-
mic loads and other types of loads applied very fast). In clays, in contrast, total
stress analysis is much more frequently used. The main reason is that the determi-
nation of pore pressures, except in simple problems, is often difficult, requiring rel-
atively sophisticated analysis. We will discuss effective and total stress analyses in
more detail in the contexts of the load-carrying capacity of shallow and deep foun-
dations (Chapters 10 and 13) and the stability of retaining structures (Chapter 16)
and slopes (Chapter 17).
4.8 Chapter Summary

Main Concepts and Equations
Stress analysis allows determination of the normal stress and shear stress acting on any
plane through a point, as long as the stresses on two different planes are known. The usual
situation is that the stresses s11, s33, and s13 acting on the perpendicular planes p11 and p33
passing through a point in the soil are known. The normal stresses s11 and s33 are positive
in compression. For determining the sign of s13, we look at the planes p11 and p33 as mak-
ing up a letter “L,” with p11 being the vertical leg of the “L,” and p33, the horizontal leg. If
the shear stress s13 acts in directions that would tend to open up the right angle of the “L,”
then s13 is positive. The stresses su and tu acting on a plane making an angle u with p11
(measured counterclockwise from p11) are given by
su 12 1s11 s33 2 12 1s11 s33 2cos 2u s13 sin 2u (4.5)
tu 12 1s11 s33 2 sin 2u s13 cos 2u (4.6)
The principal planes are the planes where the shear stresses are equal to zero. The
angles up that the principal planes make with the reference plane p11 are determined from
the equation:
tan 12u p 2
2s13
(4.7)
s11 s33
The normal stresses on these planes are the principal stresses s1 and s3, given by
s1 12 1s11 s33 2 214 1s11 s33 2 2 s 213 (4.8)
s3 12 1s11 s33 2 214 1s11 s33 2 2 s 213 (4.9)
Stress analysis problems can also be solved graphically using the Mohr circle and the
pole method; the pole is the point such that any line drawn through it parallel to a plane
where we wish to determine the stresses intersects the Mohr circle at a point with coordi-
nates equal to the stresses on that plane. A separate sign convention exists for plotting shear
stresses in a Mohr diagram. Shear stresses that are applied on the face of an element in the
counterclockwise orientation are positive for that purpose.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 168
Soils are said to have attained a plastic state of “failure” (or to have sheared) when the
principal stress ratio s*/s
1 3* satisfies the Mohr-Coulomb failure criterion, expressed mathe-
matically as
s1*
N (4.24)
s3*
where N is known as the flow number, given by
1 sin f
N (4.25)
1 sin f
and
s* s c cot f
Mohr circles can also be plotted for strains. The space used for plotting the Mohr cir-
cle of strains is defined by a horizontal axis where incremental normal strains de are plot-
ted and a vertical axis where incremental shear strains 12 dg are plotted. The pole method
applies to the Mohr circle of strains as well. The dilatancy angle c can be visualized
through the Mohr circle. The sine of this angle expresses the rate of increase in volumetric
strain for a unit increase in maximum shear strain:
de v
0 dgmax 0
sin c (4.15)
The dilatancy angle can be used to show that slip surfaces in soils with nonzero c are
log-spirals. In a soil where c 0, the slip surface is circular. The dilatancy angle is useful
in other ways, as we will see in Chapter 5.
If undisturbed, soil exists in a state of rest, for which the ratio of lateral to vertical
effective stress is the coefficient of earth pressure at rest K0. If allowed to stretch in the
horizontal direction, a level soil mass reaches an active Rankine state, in which the whole
soil is in a state of incipient plastic failure, with potential slip surfaces making angles of

(45° f/2) with the horizontal. The ratio of lateral to vertical effective stress in this case
is KA. If compressed in the horizontal direction, the soil reaches a passive Rankine state, in
which the whole soil is in a state of incipient plastic failure, with potential slip surfaces
making angles of
(45° f/2) with the horizontal. The ratio of lateral to vertical effective
stress in this case is KP. The values of K0, KA, and KP are calculated using
K 0 11 sin f2 2OCR (4.35)
1 sin f
KP N (4.38)
1 sin f
1 1 sin f
KA (4.37)
N 1 sin f
The stress-strain response and shear strength of soils are studied and measured in the
laboratory using mostly the consolidation, direct shear, triaxial, and unconfined com-
pression tests. Stress paths are a useful way of visualizing the loading of the samples in
these tests. Stress paths are the evolution of p and q (total stress paths) or p and q (effec-
tive stress paths) during loading of a soil element or test sample.
For a triaxial test,
p sm 13 1s1 2s3 2 (4.60)
saL00581_ch04.qxd 9/14/06 4:54 PM Page 169
p¿ s¿m 13 1s¿1 2s¿3 2 (4.61)
q q¿ s1 s3 (4.62)
Applications of loads on the boundary of a soil mass generate stress increments at every
point within the soil mass that have magnitude and direction that depend on the magnitudes
and directions of the loads, the geometry of the soil mass and the distance from the point to
the various loads. These stress increments have a pore pressure and an effective stress com-
ponent, the magnitudes of which depend on the draining properties of the soil and drainage
boundary conditions. An analysis in terms of total stresses is called total stress analysis.
Likewise, an analysis in terms of effective stresses is called effective stress analysis.
Symbols and Notation

Symbol Quantity represented US units SI units
B, L Dimensions ft m
c Cohesive intercept lb/ft2 kPa
E Young’s modulus lb/ft2 kPa
G Shear modulus lb/ft2 kPa
K Bulk modulus lb/ft2 kPa
M Constrained modulus lb/ft2 kPa
N Normal force lb kN
N Flow number unitless unitless
Q Vertical load lb or lb/ft kN or kN/m
qb Stress distribution on top of lb/ft2 kPa
half-space (either as an imposed
load or as a reaction to imposed
displacements by a rigid body)
r Radius of slip surface ft m
r, z Cylindrical coordinates ft m
s Shear strength of soil lb/ft2 kPa
s s coordinate of center of Mohr circle lb/ft2 kPa
T Tangential force lb kN
t Radius of Mohr circle unitless unitless
ub Back pressure lb/ft2 kPa
x, y, z Cartesian coordinates ft m
x1, x2, x3 Cartesian coordinates ft m
gij Engineering shear strain in./in. mm/mm
corresponding to directions xi
and xj, with i j
e1 Major principal strain in./in. mm/mm
e3 Minor principal strain in./in. mm/mm
e11 Normal strain in x1 direction in./in. mm/mm
eij (with i j) Normal strain in./in. mm/mm
eij (with i j) Shear strain in./in. mm/mm
ev Volumetric strain in.3/in.3 mm3/mm3
u Cylindrical/spherical coordinate deg or rad deg or rad
n Poisson’s ratio unitless unitless
p11 Plane on which s11 acts unitless unitless
(continued)
saL00581_ch04.qxd 9/14/06 4:54 PM Page 170
Symbol Quantity represented US units SI units

p33 Plane on which s33 acts unitless unitless
r Spherical coordinate ft m
s Normal stress lb/ft2 kPa
s* Transformed normal stress lb/ft2 kPa
s1 Major principal stress lb/ft2 kPa
s3 Minor principal stress lb/ft2 kPa
s11 Normal stress acting on p11 lb/ft2 kPa
s13 Shear stress acting on p33 in x1 lb/ft2 kPa
direction (numerically equal to s31)
s31 Shear stress acting on p11 in x3 lb/ft2 kPa
direction (numerically equal to s13)
s33 Normal stress acting on p33 lb/ft2 kPa
sc Cell pressure lb/ft2 kPa
sd Applied axial stress (deviatoric lb/ft2 kPa
stress)
su Normal stress acting on plane lb/ft2 kPa
making angle u with p11; where
u is measured counterclockwise
from p11
t Shear stress lb/ft2 kPa
tu Shear stress acting on plane lb/ft2 kPa
making angle u with p11; where
u is measured counterclockwise
from p11
Spherical coordinate deg or rad deg or rad
f Friction angle deg deg
c Dilatancy angle deg deg
4.9 References
References Cited
Brooker, E. W., and H. O. Ireland (1965). “Earth Pressures at Rest Related to Stress History.”
Canadian Geotechnical Journal, Vol. 2, No. 1, 115.
Cook, G. (1950). “Rankine and the Theory of Earth Pressure.” Geotechnique, Vol. 2,
271–279.
Fung, Y. C. (1993). A First Course in Continuum Mechanics. Prentice-Hall, Upper Saddle
River, NJ.
Hoek, E., and E. T. Brown (1980). “Empirical Strength Criterion for Rock Masses.” Journal
of the Geotechnical Engineering Division, ASCE, Vol. 106, No. GT9, 1013–1035.
Hoek, E., and E. T. Brown (1988). “The Hoek-Brown Failure Criterion—A 1988 Update.”
Proc. 15th Canadian Rock Mech. Symp. (ed. J.C. Curran), 31–38. University of
Toronto, Dept. Civil Engineering, Toronto.
Jaky, J. (1944). “The Coefficient of Earth Pressure at Rest.” Journal of Society of Hungarian
Architects and Engineers, October, 355–358.
Rankine W. J. M. (1857). “On the Stability of Loose Earth.” Phil. Trans. Roy. Soc. London,
Vol. 147, No. 1, 9–27.
Saada, A. (1993). Elasticity: Theory and Applications. Krieger, Melbourne, FL.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 171
Terzaghi, K. (1943). Theoretical Soil Mechanics. Wiley, New York.

Timoshenko, S. P., and J. N. Goodier (1970). Theory of Elasticity. McGraw-Hill, New York.
Additional References
Duncan, J. M., and C. Y. Chang (1970). “Nonlinear Analysis of Stress-Strain in Soils.”
Journal of Soil Mechanics and Foundation Engineering Division, ASCE, Vol. 96, No.
SM5, 1629–1653.
Fahey, M., and J. P. Carter (1993). “A Finite Element Study of the Pressuremeter Test in
Sand Using a Non-Linear Elastic Plastic Model.” Canadian Geotechnical Journal, Vol.
30, No. 2, 348–361.
Fung, Y. C. (1993). A First Course in Continuum Mechanics. Prentice-Hall, Upper Saddle
River, NJ.
Kondner, R. L. (1963). “Hyperbolic Stress-Strain Response: Cohesive Soil.” Journal of Soil
Mechanics and Foundation Engineering Division, ASCE, Vol. 89, No. SM1, 115–143.
Lee, J. H., and R. Salgado (1999). “Determination of Pile Base Resistance in Sands.” Journal
of Geotechnical and Geoenvironmental Engineering, Vol. 125, No. 8, 673–683.
Lee, J. H., and R. Salgado (2000). “Analysis of Calibration Chamber Plate Load Tests.”
Canadian Geotechnical Journal, Vol. 37, No. 1, 14–25.
Michalowski, R. L. (2005). “Coefficient of Earth Pressure at Rest.” Journal of Geotechnical
and Geoenvironmental Engineering, Vol. 131, No. 11, 1429–1433.
Parry, R. H. G. (1995). Mohr Circles, Stress Paths and Geotechnics. E&F Spon, London.
Poulos, H. G., and E. H. Davis (1974). “Elastic Solutions for Soil and Rock Mechanics.”
Wiley, New York.
4.10 Problems
Conceptual Problems
4-1 Define all the terms in bold contained in the chapter summary.
4-2* The concept of stress is crucial for geotechnical analyses. By considering soils with
different particle sizes, give an example of a problem in which it is perfectly valid to
use stresses in calculations and another in which it is not.
4-3 Consider Fig. 4-13. Demonstrate visually that there is no distortion in the soil if de1
de3 by replotting the figure for this condition.
4-4* What is a state of plastic equilibrium?
4-5* Why does a Poisson’s ratio of 0.5 imply incompressibility?
4-6* Given that foundation loading in saturated clays is applied under undrained
conditions, what can we state about the volumetric change undergone by the soil after
load application? What value of Poisson’s ratio would be used in settlement equations
for saturated clays derived using elasticity assumptions? What does that imply about
the mechanism by which immediate settlements develop in saturated clays?
4-7 When is a soil said to be in a state of rest?
4-8 What is the range of values that the coefficient of lateral earth pressure at rest can
assume for typical sands and clays?
saL00581_ch04.qxd 9/14/06 4:54 PM Page 172
4-9 What is a Rankine state? What is the difference between an active and a passive
Rankine state?
4-10 What is the pole method? If the stress state at a point in the soil mass is known, how
can the pole be determined in the corresponding Mohr circle?
4-11 How can the direction of potential slip surfaces through a point in a soil mass be
determined?
4-12* Consider drained tests (a triaxial compression and a triaxial extension test) on two
soil samples. Is there a difference between the total stress path and effective stress
path? Plot their stress paths in both the s versus t and the p versus q spaces. During
the shearing phase of the test, what is the slope of the stress path in p versus q space?
Quantitative Problems
4-13 Figure 4-45 shows the state of stress at a point within a soil deposit. Using both an
analytical approach and the pole method, determine (a) the direction of the principal
planes, (b) the magnitude of the principal stresses, and (c) the normal and shear stress
in a plane at 45° with the horizontal.
4-14 For the state of stress defined in Fig. 4-46: (a) draw the Mohr circle; (b) find the pole,
identify it in the Mohr circle, and write its coordinates below; (c) find the principal
planes and stresses (both graphically and analytically); (d) find the stresses (s and t)
on the planes making angles of
45° with the horizontal, indicate which is which,
and solve both graphically and analytically.
4-15 Repeat Problem 4-14 for the state of stress of Fig. 4-47.
4-16 Consider the stress state of Problem 4-13. Redo the problem for shear stresses of 0,
0.2, and 0.3 MPa. What effect do you observe on (a) the magnitude and (b) the
direction of the principal stresses as the applied shear stresses increase?
1 MPa
0.1 MPa
0.4 MPa
0.1 MPa 1 MPa
0.1 MPa
0.1 MPa
0.1 MPa 0.1 MPa

1 MPa
1 MPa
0.4 MPa
60 30
Figure 4-45 Figure 4-46

State of stress at a point for Problem 4-13. State of stress for Problem 4-14.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 173
200 kPa 100 kN/m
30
500 kPa
a b 100 kN/m
100 kN/m d c
500 kPa
30
200 kPa 100 kN/m
Figure 4-48
Element corresponding to Problem 4-17.
30
Figure 4-47
State of stress for Problem 4-15.
4-17* The unit block abcd is acted upon by loadings as shown in Fig. 4-48. (a) Draw the
Mohr circle of stress. (b) Locate the points on the Mohr circle indicating the stresses
on faces ad and dc. (c) Find the pole. (d) Find the directions of principal planes and
the magnitudes of the major and minor principal stresses. (e) Draw a smaller element
inside abcd to represent the principal planes and stresses acting on them. (f) Solve
the problem analytically and compare your results.
4-18 For the soil element of Example 4-5, consider the major principal strain increment to
be vertical and (a) find the pole and (b) find the directions of the zero-extension lines
(slip lines) through the element.
4-19* For the state of strain shown in Fig. 4-49, (a) draw the Mohr circle, (b) find the pole,
(c) find the principal directions and the values of the principal strain increments, and
(d) find the directions of the zero-extension lines.
4-20* Solve Problem 4-19 analytically.

4-21 A dry sand deposit has a friction angle of 35° and unit weight of 20 kN/m3. What is
the lateral effective stress and the coefficient of lateral earth pressure for (a) at-rest,
(b) active, and (c) passive conditions at a depth of 10 m?
4-22 Consider the soil mass of Fig. 4-24(c), retained in place by a semi-infinite,
dimensionless plate. Consider the soil to reach a Rankine state with a dilatancy angle
2.5 103
3 103
Figure 4-49
1.5 103 Incremental strains for Problem 4-19 expressed
103 in terms of the original and deformed shape of a
Initial lengths of element sides 1 square element with sides with length equal to 1.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 174
of 8°. If the plate moves in the direction away from the soil, resulting in an extension
corresponding to a lateral strain increment of 0.5%, what is the vertical strain in the
retained soil mass? What is the direction of the slip plane through the soil mass?
4-23 If the soil of Problem 4-22 has a friction angle of 34.5°, what is the error in the
estimate of the slip plane direction based on friction angle?
4-24* Consider the following ways in which you could shear a triaxial sample: (a) reduce
the radial stress while keeping the axial stress unchanged; (b) increase the radial
stress while keeping the axial stress unchanged; (c) increase the axial stress while
keeping the radial stress unchanged; (d) reduce the axial stress while keeping the
radial stress unchanged. Plot (total stress) p-q diagrams for all four cases. Be specific
about the directions of the total stress paths. Calculate also the following ratio for
each case: (s2 s3)/(s1 s3). What do you believe this ratio tells you about what
you are doing to the sample?
4-25 Three vertical point loads (800, 200, and 250 kN) are applied to the surface of a level
soil mass along a straight line. The loads are separated by distances of 2 m. Using the
Boussinesq equation, compute the vertical stress increase due to these loads at a
depth of 2 m along their vertical lines of action.
4-26 A 100-kN point load is applied at a point on the surface of a soil mass. Calculate
the vertical stress increase due to this load at a point with distances 1.5, 2, and 1 m
with respect to the point of application of the load in the x, y, and z directions,
respectively.
4-27 Two buildings, shown in Fig. 4-50, are separated by 7 m. The uniformly distributed
load of building 1, with dimensions 10 15 m, is equal to 300 kPa; that of building
2, with dimensions 7 25 m, is 500 kPa. Compute the vertical stress increase due to
both buildings at point P, 8 m below the center of building 1.
500 kPa
300 kPa
Cross
8m
section
P
4m
Plan P 15 m
view
6m
7m 7m 10 m
Figure 4-50
Two buildings of Problem 4-27.
saL00581_ch04.qxd 9/14/06 4:54 PM Page 175
E
4 5 6 7 8 9 10
A B C
3 20 24 25 26 21 11
4 columns at
15 m 4-m spacing
2 19 23 28 27 22 12
D
1 18 17 16 15 14 13
7 columns at 4-m spacing

27 m
0.5 m Compacted fill (g 21 kN/m3)
WT
2.0 m
3.0 m
2.5 m
11 m NC clay (gsat 17 kN/m3)
2.5 m
3.0 m
Sand (g 22 kN/m3)
(Not to scale)
Figure 4-51
Building and its foundations for analysis in Problem 4-28.
4-28 A building is to be built on top of the soil profile of Fig. 4-51. A solution based on
footings14 bearing on a layer of compacted soil is proposed. The areas of the
footings are not yet available. The average unit load to be carried by the foundations
will be approximately 100 kPa, but it is expected that it will be distributed among
three classes of columns: A, B, and C. Columns C will carry two times the load
carried by columns of class B, which, in turn, will carry 1.25 times as much load as
columns of class A. Estimate the vertical stress increase due to the building loading
for depths ranging from 2 to 13 m: (a) below point D (center of the corner column),
14 Footings are discussed in detail in Chapters 8–11.

saL00581_ch04.qxd 9/14/06 4:54 PM Page 176
(b) below point E, and (c) below the center of the building. Consider first the
average unit load of 100 kPa applied over the entire plan area of the building and
then the statically equivalent column (point) loads consistent with the 2 and 1.25
ratios discussed earlier.
4-29* In Section 4.6, we studied the strip loading (Fig. 4-40). (a) Find the principal stresses
s1 and s3 acting on principal planes p1 and p3 at point P with coordinates u1 and u2.
Express the results for s1 and s3 in terms of a. (b) Show that the points belonging to
a circle passing through P and through the two edges P1 and P2 of the strip load have
the same state of stress as P. Write the equation for the maximum shear stress for
this circle. (c) Sketch the trajectories of s1 and s3. (d) Using the equation for the
maximum shear stress you obtained in (b), find the maximum shear stress for a given
loading and the circle corresponding to it (that is, the value of a corresponding to
maximum shear stress). The equations derived in Section 4.6 for s11, s33, and s13 are
repeated below for convenience:
32 1u2 u1 2 1sin 2u2 sin 2u1 2 4

qb
s33
2p
32 1u2 u1 2 1sin 2u2 sin 2u1 2 4

qb
s11
2p
1u2 u1 2
2nqb 2nqba
s22
p p
3 1cos 2u1 cos 2u2 2 4

qb
s13
2p
a u2 u1

Stress Analysis, Strain Analysis, and Shearing of Soils

Uploaded by

Copyright:

Available Formats

Stress Analysis, Strain Analysis, and Shearing of Soils

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Stress Analysis, Strain Analysis, and Shearing of Soils

Uploaded by

Copyright:

Available Formats

saL00581_ch04.

qxd 9/14/06 4:54 PM Page 113

Ut tensio sic vis (strains and stresses are related linearly).

114 The Engineering of Foundations

4.1 Stress Analysis

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 115

116 The Engineering of Foundations

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 117

su 12 1s11  s33 2  12 1s11  s33 2cos 2u  s13 sin 2u (4.5)

tu 12 1s11  s33 2sin 2u  s13 cos 2u (4.6)

Principal Stresses and Principal Planes

118 The Engineering of Foundations

s1 12 1s11  s33 2  214 1s11  s33 2 2  s 213 (4.8)

s3 12 1s11  s33 2  214 1s11  s33 2 2  s 213 (4.9)

3su  12 1s11  s33 2 4 2  t 2u 14 1s11  s33 2 2  s 213 (4.10)

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 119

Recalling the equation of a circle in Cartesian

Line parallel to p11

120 The Engineering of Foundations

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 121

100 kPa Pole

122 The Engineering of Foundations

Solving Stress Analysis Problems

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 123

The absolute value of 16.8° is less than 45°, so

u p3 16.8°  90° 73.2°

s1 12 1200  50 2  12 1200  50 2cos32 116.8°2 4  1502sin 32 116.8°2 4 215.1 kPa

s3 12 1200  50 2  12 1200  50 2cos32 173.2°2 4  1502sin 32 173.2°2 4 34.9 kPa

Now using Eqs. (4.8) and (4.9):

s1 12 1200  502  214 1200  50 2 2  1502 2 215.1 kPa

s3 12 1200  502  214 1200  50 2 2  1502 2 34.9 kPa

Stresses on Planes Making Angles 15° with Horizontal

su 125  75 cos32 115°2 4  50 sin 32 115°2 4 214.9 kPa

tu 75 sin 3 2 115°2 4  50 cos32 115°2 4 5.8 kPa

124 The Engineering of Foundations

Graphical Solution Using the Pole Method

Total and Effective Stresses

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 125

where i,j 1, 2, 3 reference directions, xi coordinate in the i direction, ui

126 The Engineering of Foundations

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 127

128 The Engineering of Foundations

Initial outer element

Final outer element Initial configuration Final configuration

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 129

Mohr Circle of Strains

130 The Engineering of Foundations

where gz shear strain in the x z1-x zp z

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 131

132 The Engineering of Foundations

4.3 Failure Criteria, Deformations,

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 133

134 The Engineering of Foundations

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 135

−12 (s1*  s*3)

The principal stress ratio s*1 /s3* easily follows:

136 The Engineering of Foundations

Figure 4-21 Q Rigid blocks

CHAPTER 4 Stress Analysis, Strain Analysis, and Shearing of Soils 137

su 12 1s11 s33 2 12 1s11 s33 2cos 2u s13 sin 2u (4.5)

tu 12 1s11 s33 2sin 2u s13 cos 2u (4.6)

s1 12 1s11 s33 2 214 1s11 s33 2 2 s 213 (4.8)

s3 12 1s11 s33 2 214 1s11 s33 2 2 s 213 (4.9)

3su 12 1s11 s33 2 4 2 t 2u 14 1s11 s33 2 2 s 213 (4.10)

The absolute value of 16.8° is less than 45°, so

u p3 16.8° 90° 73.2°

s1 12 1200 50 2 12 1200 50 2cos32 116.8°2 4 1502sin 32 116.8°2 4 215.1 kPa

s3 12 1200 50 2 12 1200 50 2cos32 173.2°2 4 1502sin 32 173.2°2 4 34.9 kPa

s1 12 1200 502 214 1200 50 2 2 1502 2 215.1 kPa

s3 12 1200 502 214 1200 50 2 2 1502 2 34.9 kPa

su 125 75 cos32 115°2 4 50 sin 32 115°2 4 214.9 kPa

tu 75 sin 3 2 115°2 4 50 cos32 115°2 4 5.8 kPa

−12 (s1* s*3)

The principal stress ratio s1 /s3 easily follows:

K 0 K 0, NC 2OCR 11 sin f2 2OCR (4.35)

EA2 CA2 CE 2 OC 2 sin2 f OC 2 sin2 a g

EA OC2sin2 f sin2 a g OC2cos2 a g cos2 f EP (4.44)

T¿A cos ag 2cos2 ag cos2 f

dt 12 1ds1 ds3 2 (4.51)

s¿ 12 1s¿1 s¿3 2 s u (4.52)

t¿ 12 1s¿1 s¿3 2 t (4.53)

2 1s1 s3 2 2 1s1 s2 2 2 1s2 s3 2 2

In a triaxial test, s2 s3; it follows that p, q, p, and q take the following def-

p¿ s¿m 13 1s¿1 2s¿3 2 (4.61)

11 2n2Q 1r2 z2 2 1>2

150 2 3 4226 1 26 2 3 4 226

30.748 tan1 11.0372 4