Soil Mechanics II
Soil Mechanics II
Soil Mechanics II
Faculty of Technology
CHAPTER ONE
STRESS DISTRIBUTION IN SOIL
1.1
Introduction
Bossinesq (1885) has given the solution for the stresses caused by the application of a
point load at the surface of a homogenous, elastic, isotropic and semi-infinite medium.
Q
x
r
R
z zx
z
A (x,y,z)
zy
y
xy
yz
yx
r = x2 + y2
R = r 2 + z2 = x2 + y2 + z2
r
r rz
xz
x
r
z
; cos =
R
R
v = poisson' s ratio
sin =
Z
Faculty of Technology
In rectangular coordinates
3Q z 3
z =
2 R 5
3Q x 2 z 1 2v
(2 R + z ) x 2
z
1
x =
5 +
2 R
3 R ( R + z ) R 3 ( R + z ) 2 R 3
(1.1)
(1.2)
3Q y 2 z 1 2v
(2 R + z ) y 2
z
1
y =
3
5 +
2 R
3 R( R + z ) R ( R + z ) 2 R 3
xy =
(1.3)
3Q xyz 1 2v (2 R + z ) xy
2 R 5
3 R 3 (R + z) 2
(1.4)
3Q xz 2
2 R 5
3Q yz 2
yz =
2 R 5
In cylindrical coordinates
3Q z 3
z =
2 R 5
Q 3zr 2
1 2v
5
r =
2 R
R( R + z )
xz =
(1.5)
(1.6)
(1.7)
(1.8)
Q
1
z
(1 2v)
3
2
R( R + z ) R
(1.9)
3Q rz 2
(1.10)
2 R 5
In most foundation problems it is very necessary to be acquainted with the increase in
vertical stresses (for settlements) and the increase in shear stresses (for shear strength
analysis).
The equation for the vertical stress can be rewritten as follows
rz =
3Q z 3 3Q
3Q
z3
z =
5 =
2
=
2 52
2 R
2 (r + z )
2 z 2
= Iz
Q
z2
2
1 + (r z )
52
(1.11)
52
3
1
Iz =
2
2 1 + (r z )
Values of Iz for different values of r/z ratio can be tabulated (Table 1.1) or plotted as Iz
versus r/z ratio (Fig 1.2) and hence can be used for routine stress calculation. Note that
the maximum vertical stress is observed directly below the load (r=0).
Table 1.1
Iz
r/z
0 0.4775
0.2 0.4329
0.4 0.3295
0.6 0.2214
0.8 0.1386
1.0 0.0844
1.2 0.0513
1.4 0.0317
1.6 0.0200
1.8 0.0129
2.0 0.0085
2.2 0.0058
2.4 0.004
2.6 0.0028
2.8 0.0021
3.0 0.0015
Iz
Faculty of Technology
0.5
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0
0.5
1.5
2.5
3.5
r/z
Fig 1.2 Non dimensional influence factor
(1.11)
Pressure distributions
Graphical vertical stress distribution on a horizontal plane at any depth z below the ground
surface can be drawn as shown here under.
The vertical stress on a horizontal plane at a depth z is given by
Q
z = Iz
2
For several assumed values of r, r/z is calculated and Iz is found for each, the value of z
is then computed. As an example consider the following cases
Case 1: Let z =c (where c= constant number) and r be varied as 0, 0.25c, 0.5c,
0.75c, c, 1.25c, 1.5c, 1.75c, 2c, etc
Case 2: Let z=2c and r be varied as case 1
Case 3: Let z=4c and r be varied as case 1
Now the stresses are calculated (tabulated) and plotted together below.
2
zc /Q
r
0
0.25c
0.5c
0.75c
c
1.25c
1.5c
1.75c
2c
Faculty of Technology
r
0
0.25c
0.5c
0.75c
c
1.25c
1.5c
1.75c
2c
Case 2 (z=2c)
(zc2)/Q
Iz
r/z
0
0.4775 0.1194
0.125 0.4593 0.1148
0.250 0.4103 0.1026
0.375 0.3436 0.0859
0.500 0.2733 0.0683
0.625 0.2094 0.0523
0.750 0.1565 0.0391
0.875 0.1153 0.0288
1.000 0.0844 0.0211
r
0
0.25c
0.5c
0.75c
c
1.25c
1.5c
1.75c
2c
Case 3 (z=4c)
(zc2)/Q
Iz
r/z
0
0.4775
0.0298
0.0625 0.4728
0.0296
0.125 0.4593
0.0287
0.1875 0.4380
0.0274
0.25
0.4103
0.0256
0.3125 0.3782
0.0236
0.375 0.3436
0.0215
0.4375 0.3082
0.0193
0.5
0.2733
0.0171
0.50
0.45
0.40
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
z=c
z=2c
z=4c
-3
-2
-1
r/c
From the above curves one can see that the stress diminishes as we move down from the
ground surface and also as we move away from the point of load application.
Similarly vertical stress distribution on a vertical plane at any radial distance r from the load
can be drawn as shown here under.
For several assumed values of z, r/z is calculated and Iz is found for each, the value of z
is then computed. As an example consider the following cases
Case 1: Let r=a (where a= constant number) and z be varied as 0, 0.5a, a, 2a,
5a, 10a, etc
Case 2: Let r=2a and z be varied as case 1
Case 3: Let r=4a and z be varied as case 1
Now the stresses are calculated and plotted below.
z
0
0.5a
a
2a
5a
10a
Case 1 (r=a)
Iz
r/z
2.00 0.0085
1.00 0.0844
0.50 0.2733
0.20 0.4329
0.10 0.4657
(za2)/Q
indet.
0.0342
0.0844
0.0683
0.0173
0.0047
z
0
0.5a
a
2a
5a
10a
Case 2 (r=2a)
Iz
r/z
4.0 0.0004
2.0 0.0085
1.0 0.0844
0.4 0.3295
0.2 0.4329
(za2)/Q
indet.
0.0016
0.0085
0.0211
0.0132
0.0043
z
0
0.5a
a
2a
5a
10a
Case 3 (r=4a)
Iz
r/z
8
0.0000
4
0.0004
2
0.0085
0.8
0.1386
0.4
0.3295
(za2)/Q
indet.
0.0001
0.0004
0.0021
0.0055
0.0033
2
za /Q
0.00
0
0.02
0.04
0.06
2
4
z/a
0.08
0.10
r=a
r=2a
r=4a
6
8
10
12
Faculty of Technology
Isobars
Z
Fig 1.3 Isobar diagram
The procedure for plotting an isobar is as follows. Let it be required to plot an isobar for
which z=0.2Q per unit area (20% isobar)
From = I
z
z
Q
z
Iz =
Iz =
0.2 Q z
2
= 0.2 z 2
Assuming various values for z, the corresponding Iz- values are computed; for these
values of Iz, the corresponding r/z-values are obtained; and, for the assumed values of z,
r-values are got. It is obvious that, for the same value of r on any side of the z-axis, or
line of action of the point load, the value of z is the same; hence the isobar is
symmetrical with respect to this axis.
When r=0, Iz=0.4775; the isobar crosses the line of action of the load at a depth
of: I z = 0.2 z 2 z 2 = I z / 0.2 z = I z / 0.2 = 0.4775 / 0.2 = 1.545 units
The calculations are best performed in the form of a table as given below. The plot is
shown too.
z
0.5
1.0
1.5
1.5452
Iz
0.0500
0.2000
0.4500
0.4775
r/z
1.211
0.645
0.155
0.0
r (unit)
0.605
0.645
0.232
0.000
0.2Q
0.2Q
0.2Q
0.2Q
-1.0
0
-0.5
0.0
0.5
1.0
0.4
0.8
1.2
1.6
2
Westergaard (1938) has obtained an elastic solution for stress distribution in soil under a
point load. He has assumed the soil to be laterally reinforced by numerous, closely spaced
horizontal sheets of negligible thickness but of infinite rigidity, which prevent it from
5
Faculty of Technology
undergoing lateral strain. The vertical stress z caused by a point load, as obtained by
Westergaard, is given by (all notations as above):
HL\CH1.doc
z =
1 1 2v
2 2 2v
32
2
2
z
1 2v r
+
2 2v z
(1.12)
q/unit length
x
z
A
x
Z
z =
2q z 3
( z 2 + x 2 )2
(1.13)
x =
2q x 2 z
( z 2 + x2 )2
(1.14)
xz =
2q z 2 x
( z 2 + x2 )2
(1.15)
3
z =
B
q/unit length
o
x
y
-z
(1.13a)
y
1
y
3q z
22+ x 2 )2 2
3
2
2
2
2 2
2(
(z
x +y +z
x +y +z
z
x
A
z =
or
q
2 z m 2 + 1
q
I
z =
z
where,
I=
3n
n
2
2
2
2
n + m +1 n + m +1
(1.13b)
(1.13c)
1
2 m 2 + 1
3n
n
2
2
2
2
n + m +1 n + m +1
(1.13d)
q/unit area
x
B=2b
z =
x =
[ + sin cos( + 2 )]
(1.16)
[ sin cos( + 2 )]
(1.17)
q
z
xz = sin sin( + 2 )
(1.18)
x A
Fig 1.5: Uniform load
The vertical stresses at different depths
below the center of a uniform load of intensity q and width B are as follows:
Depth z
0.1B
0.2B
0.5B
B
2B
5B
10B
z
0.997q 0.977q 0.818q 0.550q 0.306q 0.126q 0.064q
Faculty of Technology
A few typical pressure bulbs for this case of strip loading are shown in Fig 1.6.
B=2b
q/unit area
z/q=1/8
z/q=1/2
z/q=1/4
qx
1
sin 2
B
2
q x
z R
1
x = ln 1 + sin 2
B R 2 2
B
2
z =
q/unit area
R1
R2
z
xz =
q
2
(1.19)
(1.20)
2z
1 + cos 2 B
(1.21)
z = q 1
2
1 + ( R / z )
q/unit area
2R
q/unit area
2R
(a)
q/unit area
R
(1.22)
z
A
(b)
(1.23)
Faculty of Technology
R
2R
0.15q
3R
4R
0.05q
0.1q
Faculty of Technology
1.7
Stress Due to Uniformly Loaded Rectangular Area
The more common shape of a loaded area in foundation engineering practice is a
rectangle, especially in the case of buildings. Based on Bossinesqs theory, Newmark
(1935) has given an expression for the vertical stress at a point below the corner of a
uniformly loaded rectangular area (Fig 1.11) as
rea
a
t
ni
q/u
m=B/z
n=L/z
A
Fig 1.11: Vertical stress at the corner of a uniformly loaded rectangular area
q
z =
4
2
2
2
2
2
2
(2mn m + n + 1) (m + n + 2) + tan 1 (2mn m + n + 1)
2
m 2 + n 2 + 1 + m 2n 2 (m 2 + n 2 + 1)
m + n 2 + 1 m 2n 2
(1.24)
Since this equation is symmetrical in m and n, the values of m and n are interchangeable.
The above equation can be rewritten in the form:
z = qI
(1.25)
1
4
2
2
2
2
2
2
(2mn m + n + 1) (m + n + 2) + tan 1 (2mn m + n + 1)
2
m 2 + n 2 + 1 + m 2n 2 (m 2 + n 2 + 1)
m + n 2 + 1 m 2n 2
Based on this equation, Fadum (1941) has prepared a chart for the influence values for
sets of values for m and n, as shown in Fig 1.12.
Faculty of Technology
Fig 1.12 Values of influence factor I for calculating vertical stress at the corner
By superposition we can determine stress at any point within or out side the loaded area.
For instance the stresses at point A and point B are
U
T
I
II
IV
B
V
III
P
PWBT:I
SVBT:II
QWBU:III
RVBU:IV
At po int A : z = q I
+I
+I
+I
II
III
IV
I
At po int B : z = q I
I
I
+I
II
III
IV
I
10
Faculty of Technology
shown to be
1 + ( R / z ) 2
z = q 1
3/ 2
z
q
2 3
1
(1.26)
If a series of values is assigned for the ratio z/q, such as 0, 0.1, 0.2, ., 0.9 and 1.0, a
corresponding set of values for the relative radii, R/z, may be obtained. If a particular
depth is specified, then a series of concentric circles can be drawn. Since the first has zero
radius and the 11th has infinite radius, in practice only 10 circles are drawn. Each ring or
annular space causes a stress of q/10 at a point beneath the center at the specified depth z,
since the number of annular spaces (m) is 10.
The relative radii can be tabulated as shown below:
S. No of circle z/q Relative Radii (R/z)
1
2
3
4
5
6
7
8
9
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.000
0.270
0.400
0.518
0.637
0.766
0.918
1.110
1.387
Now lets assume that a set of equally spaced rays, say n in number, is drawn emanating
from the center of the circles, thus dividing each annular area in to n sectors, and the total
area in to m*n sectors. If the usual value of 20 is adopted for n, the total number of
sectors in this case will be 10x20=200. Each sector will cause a vertical stress of 1/200th
of the total value at the center at the specified depth and is referred to as a mesh or an
influence unit. The value 1/200 or 0.005 is said to be the influence value (or influence
factor) for the chart.
Construction of the Newmarks chart
1. For the specified depth (say 20m), calculate the radii of the circles as:
R/z 0 0.270 0.4 0.518
R (m) 0 5.40 8.00 10.36
0.637
12.74
0.766
15.32
1.908
38.16
Faculty of Technology
A
20cm
Influence value =0.005 B
Q1
Q4
R3
R4 R2
B
R1
A
Fig 1.14: Equivalent Point Load Method
Referring to Fig 1.14, if the influence values are I1, I2, I3, . for the point loads Q1,
Q2, Q3,, then z is found as:
z = (Q1I1+ Q2I2+ Q3I3+..)
(1.27)
This method gives a good result if z/B> 3.
12
Faculty of Technology
Faculty of Technology
and base embedment depth D. The influence factor Is (see Figure 3.1 for identification of
terms) can be computed using equations given by Steinbrenner (1934) as follows:
1 2v
I S = I1 +
I2
with I1 and I 2 as follows:
1 v
2
2
(1 + M 2 + 1)( M 2 + N 2 )
1
+ ln (M + M + 1)( 1 + N )
I1 = M ln
M (1 + M 2 + N 2 + 1)
M + M 2 + N 2 + 1
N
M
I2 =
tan 1
tan 1 in radians
2
2
2
N M + N +1
where; M =
B' =
L'
B'
and
N=
H
B'
B
for center and B' = B for corner I i ;
2
L' =
L
for center and L' = L for corner I i
2
IF = influence factor from the Fox (1948b) equations, which suggest that the settlement
is reduced when it is placed at some depth in the ground, depending on Poissons ratio
and L/B. Figure 3.1 can be used to approximate IF.
Note: if your base is "rigid" you should reduce the Is factor by about 7 percent (that is, Is,
rigid = O.931 Is, flexible)
Figure 1.17: Influence factor IF for footing at a depth D. Use actual footing width and
depth dimension for this D/B ratio.
14
Faculty of Technology
CHAPTER TWO
SHEAR STRENGTH OF SOILS
2.1 Introduction
HL\CH2.doc
The shear strength of a soil is its maximum resistance to shear stresses and a consequent
tendency for shear deformation. Its value determines such factors as the stability of slopes
(cuts, earth dams), the allowable bearing capacity of a foundation and the thrust of a soil
against retaining walls. Knowledge of the shear strength is thus an essential prerequisite
to any analysis concerned with the stability of a soil mass.
Basically, a soil drives its shearing strength from the following:
1) Resistance due to the interlocking of particles.
2) Frictional resistance between individual soil grains, which may be sliding
friction, rolling friction, or both.
3) Adhesion between soil particles or cohesion
Granular soils or sands may derive their shear strength from the first of the two sources,
while cohesive soils or clays may derive their shear strength from the second and third
source. Highly plastic clays, however, may exhibit the third source alone for their
shearing strength. Most natural soil deposits are partly cohesive and partly granular and
as such, may fail in to the second of the three categories just mentioned, from the point of
view of shearing strength.
2.2 Basic Concepts of Shearing Resistance
a) Friction between Solid Bodies
Consider a prismatic block resting on the surface MN (Fig 2.1). The block is subjected to
two forces: force Pn acting at right angle to surface MN and force Fa acting tangentially to
the plane MN.
As Fa increases, Fr also increases such that Fa=Fr. The block will start sliding
when angle of obliquity reaches a maximum value, m.
o If the block and surface MN are the same material, m= (=angle of
internal friction) and thus tan = ( =coefficient of friction)
o If the block and surface MN are different materials, then m= angle of skin
(wall) friction.
Faculty of Technology
b) Sign convention
Now lets consider an element of soil whose sides are chosen as the principal planes
(planes which do not have shear stresses), the major and the minor, as shown in Fig 2.3
(a). Let it be required to determine the stress conditions on a plane inclined at angle
from the major principal plane measured CCW (Fig 2.3(b)). Considering the equilibrium
of this element redrawn as Fig 2.3 (c), we have
= 1 cos2 + 3 sin 2 = 3 + (1 3 ) cos2
( ) ( )
= 1 3 + 1 3 cos 2
2
( )
1 3 sin 2
2
(2.1)
(2.2)
16
Faculty of Technology
2
+3
+ 3 1 3
3
2 1
cos 2 + 2 = 1
2 1
2
2
+ 3
3
+ 2 = 1
1
2
(2.3)
1
Coulombs Envelope
f = c + f tan
(2.4)
f = c + f tan
Where: is the total stress normal to the failure plane
is the angle of internal friction
17
Faculty of Technology
The use of Coulombs equation didnt always result in successful design of soil
structures. This is due to the fact that water cannot sustain shear stress; the shear
resistance of a soil must result solely from the frictional resistance arising at the particle
contact points, the magnitude of which depends solely on the magnitude of the effective
stress carried by the soil skeleton. Hence, expressing Coulombs equation in terms of
effective stress:
f = c' + 'f tan'
(2.5)
in which the parameters c' and ' are properties of the soil skeleton, referred to
respectively as the effective cohesion and the effective angle of friction.
b) Mohr-Coulomb Theory of Shear strength
The Mohr- Coulomb theory of shearing strength of soil, first postulated by Coulomb
(1976) and later generalized by Mohr, is the most commonly used concept. Mohr stated
that when shear stress on the failure plane reaches a unique function of the normal stress
on the failure plane, failure takes place. Thus if we conduct several tests and obtain
principal stresses at failure, the Mohrs circles at failure can be constructed as in Fig 2.5.
A line drawn tangent to the failure circles is called an envelope (Mohr failure envelope).
Therefore, it can be said that the Mohrs circle of stress relating to a given stress
condition would represent, incipient failure condition if it just touches or is tangent to the
strength or failure envelope (circles C, D, E); otherwise, it would wholly lie below the
envelope as shown in circle A. circles lying above the Mohr failure envelope (circle B)
cannot exist.
Faculty of Technology
can be placed if required to allow the specimen to drain. For undrained tests, metal plates
and solid metal grilles may be used.
And the following plots can be produced (Example Fig 2.7 for overconsolidated clay)
2) Shear strain vs. shear stress (h vs. )
1) Normal stress vs. shear stress ( vs. )
19
Faculty of Technology
20
Faculty of Technology
21
Faculty of Technology
V
V 1
V
o
V V
o
A= o
=A
=
o
h h
h
o
h 1
o h
o
V=A (ho - h) = Vo - V
A
o .
1
a
V
o
1
a
1
1
factor. Once the corrected area is determined, the additional axial stress or the deviator
stress, , is obtained as
= =
1 3
reading)
22
Faculty of Technology
( ) / 2
( )
1 3
1 3
=
c cot + ( + ) / 2 2c cot + ( + )
1 3
1 3
(1 3 ) = 2c cos + (1 + 3 ) sin
1 (1 sin ) = 3 (1 + sin ) + 2c cos
sin =
1 =
3 (1 + sin ) 2c cos
+
(1 sin )
(1 sin )
1 = 3 N + 2c N
(2.6)
These state of stress at failure is defined as plastic equilibrium condition, when failure is
imminent.
From one test, a set of 1 and 3 is known; but it can be seen from equation 2.6 that at
least two such sets are necessary to evaluate the parameters c and . Conventionally, three
or more such sets are used from a corresponding number of tests. However, the usual
procedure is to plot the Mohrs circles for a number of tests and take the best common
tangent to the circles as the strength envelope. A small curvature occurs in the strength
envelope of most soils, but since the effect is slight, the envelope for all practical
purposes, may be taken as a straight line. The intercept of the strength envelope on the axis gives cohesion and the angle of the slope of this line with the -axis gives the angle
of internal friction.
In derivation of equation 2.6, we have ( 1 3 ) = 2c cos + ( 1 + 3 ) sin
( 1 3 ) / 2 = c cos + (( 1 + 3 ) / 2) sin
q = c cos + p sin
which plots as a straight line as q against p (Fig 2.9), the slope of the line defining sin
and the intercept ccos, and thus provide an alternate method for determining shear
strength parameters from the results of triaxial tests.
Faculty of Technology
This test is applicable to the short term stability analysis of works formed on or in clay
deposits where it is considered that insufficient time has elapsed for any dissipation of
excess porewater pressure to have occurred by the end of construction. Such works
generally include small embankments, cuttings, retaining walls and foundation of
buildings.
Test on Saturated Clay
If the clay sample is saturated, the increase in cell pressure is carried entirely by the
porewater as an excess porewater pressure with no change in effective stress in the
sample and therefore no change in shear strength. Thus, the deviator stress required to
fail the sample is independent of the cell pressure at which the test is run. Fig 2.10 shows
the corresponding Mohr stress circles, the common tangent to which defines the failure
envelope for the soil, which in this case is horizontal giving u=0; the intercept on the
vertical shear stress axis defines the undrained cohesion Cu. The undrained shear
strength, u, in terms of the total stress is thus u= Cu + tanu, with u =0,
fu = Cu.
(2.7)
The undrained cohesion thus defines the undrained shear strength. And Cu is more
generally referred to as the undrained shear strength.
24
Faculty of Technology
correlationHL\CH2.doc
Fig 2.11 Failure envelopes for undrained tests on partially saturated clay w.r.t (a) total stress (b) effective stress
25
Faculty of Technology
The greater the pressure to which a sample is consolidated, the greater the deviator stress
required to cause failure. Fig 2.13(a) shows typical Mohr circles of total stress. The
intercept and slope of the failure envelope define the total stress shear strength parameters
for the soil, which for a consolidated-undrained test are denoted by Ccu and cu
respectively. If the porewater pressure is measured during the test, as is the usual practice,
then the corresponding Mohr circles of effective stress may be drawn (Fig 2.13(b)), the
failure envelope now defining the effective stress shear strength parameters C and '.
The effective stress circles may lie either to the left or right of their respective total stress
circles, depending on whether the porewater pressure is positive or negative.
Fig 2.12: Failure envelopes for consolidated-undrained tests on clay w. r. t (a) total stress, (b) effective stress.
26
Faculty of Technology
For normally consolidated clays, the failure envelopes with respect to total and effective
stress pass through the origin giving Ccu and C' equal to zero; for overconsolidated clays
Ccu and C' usually range between 5 and 30 KPa. The value of ' is not influenced
appreciably by overconsolidation and ranges from about 30o to 20o, decreasing with
increasing plasticity index. The value of cu varies similarly and its relation to ' is
determined by the magnitude of the porewater pressure at failure.
In the standard triaxial test the sample is failed in a conventional manner by holding 3
constant and increasing 1. In a field situation the stress changes leading to a potential
failure of an element are more complex than this and in general range from 3 constant
and 1 increasing to 1 constant and 3 decreasing. It is important to consider the possible
influence of this on the measured values of the shear strength parameters. This is
illustrated in Fig 2.14 for the case of normally consolidated clay. In test 1 a soil sample is
fully consolidated under a cell pressure 3.
Faculty of Technology
the shearing stage. The form of the stress-strain-volume change relationship for sands and
clays are similar to those obtained from direct shear tests and are given in Fig 2.15.
(a) Sand
(b) Clay
Fig 2.15: Stress-strain-volume change relationships for drained tests (a) Sands (b) clay
Since u=o through out the shearing process, ' 3=3, ' 1=1 and the Mohr circles of
effective stress and total stress coincide (Fig 2.16). The failure envelope defines the
effective stress parameters C' and '. These are often denoted Cd and d respectively.
Generally there is little difference between the effective stress parameters obtained from
drained and consolidated-undrained tests. However, for sands and heavily
overconsolidated clays d is slightly higher than ', because of the work done by the
sample as it expands against the confining pressure during shear.
28
Faculty of Technology
To introduce the idea of pore pressure parameters, consider the soil element shown in Fig
2.17 wherein the element is subjected to a triaxial loading in which 2=3. We may
consider this stress system to be composed of an isotropic change of stress 3 plus a
uniaxial change of deviator stress 1-3. This, of course, is the stress system imposed
in the triaxial test and it is clear that the excess porewater pressure generated within the
element will result firstly from the change in the all-round stress and secondly from the
change in the deviator stress.
Af
1.5-2.5
0.7-1.3
0.3-0.7
-0.5-0
29
Faculty of Technology
1 +
2
q =
Actually the p and q values represent the coordinates of a point on Mohrs circle, with p
representing the center of the circle (always located on the abscissa, or axis) and q
representing the maximum shear stress, equal to the radius of the circle. The locus of p-q
points for a test series is known as a stress path. Such a graphical representation is known
as the p-q diagram. Fig 2.18 shows the Mohr envelope ( line) developed from tangents
to the circles at points 1, 2, and 3 and line Kf (Kf line), which passes through p-q points
A, B, and C- points of maximum shear for the respective circles. Thus, the Kf line
represents a limiting state of stress at impending failure. The following relationship
between the line and the Kf line exist as illustrated in Fig 2.9.
sin = tan
(2.11)
d = c cos ,
where d is the intercept on the q-axis and is the angle made by the Kf line with the paxis.
Remark
0.5
increses
0
Triaxial compression
0
increases Triaxial extension
decreases increases
Point A:
30
Faculty of Technology
ecr depends on 3 in the triaxial test and v in the direct shear test. The larger
3 and v, the smaller ecr.
Saturated Cohesionless Fine-Grain Soils
The angle of internal friction of saturated sands and some inorganic silts is only slightly
less than that of the soil in a dry state and of the same relative density (as the drainage is
instantaneous). However the shear strength might be altered significantly by a change in
the pore pressures (equation 2.12). Quite apparently, when the porewater pressure
approaches , the shear strength approaches zero. When that happens, we may approach
impending instability or perhaps motion (e.g., slope failures, boiling). Fluctuation in the
water table is a common cause of significant variations in the pore stress and, thereby, in
the shear strength of the soil.
31
Faculty of Technology
If a saturated and/or inorganic silt is totally saturated and under hydrostatic neutral stress
such that it is not subjected to any effective stress, the mass is in a state of liquefaction.
Under such circumstances, the pore water pressure u equals the total normal stress ,
there by reducing the shear strength to zero, and then the soil is in a quick or boiling or
flowing condition (conditions that result in impending upward movement of soil and
water).
If submerged fine sand undergoes a sudden decrease in the void ratio, an increase in the
porewater pressure, u, may result such that the pore pressure may equal or exceed the
value for . For example, pile driving, earthquakes, blasts, or other forms of vibration or
shock may cause a sudden decrease in the volume, thereby increasing the pore pressure u
as a result of a surge in hydrostatic excess pressure. Should the value of u reach sufficient
magnitude, say u > , the shear strength of the soil may be totally lost, resulting in what
is known as spontaneous liquefaction. Loose, fine silty sands are most vulnerable to such
effects from shock or dynamic loads or sudden fluctuations in the water table.
Compacting loose sand stratum is frequently a viable option to decrease the possibility of
liquefaction.
HL\CH2.doc
Shear strength behavior of clays is influenced by the fact whether the clay is normally
consolidated or overconsolidated, by the fact whether it is undisturbed or remolded, by
the drainage conditions during testing, consistency of the clay, by certain structural
effects, by the type of test and by the type and rate of strain.
Normally Consolidated Clays
When the sample is extracted from the ground, the overburden pressure is removed and
the pore pressure altered significantly, that is, negative pore pressures are developed. In
order to simulate a somewhat realistic in-situ state of stress, the characteristics of
saturated, normally consolidated clays extracted from a given stratum are commonly
investigated via a CU triaxial test. A confining pressure, 3, and a deviator stress, , are
applied for undrained conditions. A confining pressure, say o, of the in-situ value may
be estimated as the overburden pressure for the depth of the stratum from which the
sample was extracted. If several such tests are run for varying confining cell pressures, a
Mohr envelope may be obtained as indicated in Fig 2.20 (a). If the confining pressure is
less than the in-situ value o, the Mohr envelope depicts a range of preconsolidation of
the soil; that is, relative to the confining pressures, the soil specimen appears
overconsolidated. The shear strength of the clay specimen tested in this range is higher
than that indicated by a straight line through the origin. The relationship between the
shear strength and the normal stress in this range is designated by line portion ab, which
is slightly curved, but frequently interpreted as a straight line. On the other hand, if the
tests are run under confining cell pressures larger than o, the envelope to the rupture
circles is approximately a straight line, represented by segment bc. Although the effect of
porewater pressure is present for all ranges of normal stress within the Mohr envelope,
32
Faculty of Technology
the pore pressure is larger for the case where the confining pressures are smaller than o.
Generally, there is an increased drainage of water with increasing confining pressure.
However, the effective stress is more significant in the region where the confining
pressures are larger than o.
Fig 2.20 relationships from CU and drained tests (a) CU test (b) CD test
Figure 2.20 (b) shows the relationship between the shear strength and the effective stress
for a case of CD tests. These tests are appreciably more time consuming than the CU test,
and are therefore not as common during investigations. The specimen is subjected to a
confining pressure, and then the deviator stress is applied at a very slow rate and drainage
is permitted at each end of the sample. A series of such tests run under varying confining
cell pressures provides a strength envelope as shown in Fig 2.20 (b). The strength
envelope for such test is again somewhat curved for the range of confining pressures less
than o; for confining pressures greater than o, the Mohr envelope is approximately a
straight line. Briefly, a comparison of the two cases shows the following
1. The corresponding slopes of straight-line segments of the Mohr envelopes (Fig 2.20)
are significantly different.
2. For the drained condition the neutral (pore-water) stress is virtually negligible.
3. The effective friction is significantly larger for the drained case than for the undrained
case.
The relationship between angles and cu may be illustrated by means of Fig 2.20 (c).
The shear strengths for the CD and CU tests are and 1 as shown. The actual strength
may lie somewhere between that for a CD and that for an undrained condition since total
drainage is unlikely; that, the actual strength of the sample may range between values of
HL\CH2.doc
and 1.
Overconsolidated Clays
HL\0.doc
A soil being evaluated or tested may have been subjected to a great deal of
precompression (i.e. overconsolidation pressure, Pc) induced by loads, which since then
may have disappeared. Fig 2.21 shows the Mohr strength envelope for an
overconsolidated clay. One notes that for < Pc, the Mohr envelope deviates from a
straight line, that is, the shear strength is larger than that given by the dotted straight line.
Generally preconsolidation causes or results in a smaller void ratio at failure than would
otherwise exist, even though the specimen tends to expand as a result of extraction from
an in-situ condition. Cohesion and general capillary forces tend to resist the volume
increase, thereby resulting in somewhat greater shear strength, as indicated by the curved
portion of the envelope. Beyond the preconsolidation pressure Pc, the effective normal
stress and shear strength relationship is given by a reasonably straight line.
33
Faculty of Technology
Fig 2.21 Shear strength vs. effective normal stress for OCC
The shear strength of clays that have fine discontinuities, hairline cracks, and the like,
generally referred to as fissures, may be appreciably different from that of the typical
overconsolidated clay described above. Depending on the magnitude and orientation of
these fissures, test results may be particularly misleading in the overall evaluation. For
example, the results from a direct shear test on a sample where fissures are parallel to the
shearing force may be appreciably smaller than for one where the orientation of fissures
is 90o to the shear force. On the other hand, the triaxial test yields somewhat more reliable
results, improved perhaps by the lateral restraint of the confining pressure.
Sensitivity of Clay
If the strength of undisturbed sample of clay is measured and its strength is again
measured after remolding at the same water content to the same dry density, a reduction
in strength is often observed. This is an important phenomenon which is quantitatively
characterized by sensitivity, defined as follows:
Unconfined compression strength, undisturbed
Sensitivity S f =
Unconfined compression strength, remolded
A comparison of stress-strain curves for sensitive clay in the undisturbed and remolded
sates is shown in the figure below. Sensitivity classification is given in the table below.
Table: Sensitivity classification of clays
Sensitivity, Sf
1
1-2
2-4
4-8
8-16
>16
Classification
Insensitive
Low
Medium
Sensitive
Extra-sensitive
Quick (Sf can be up to 150)
Overconsolidated clays are rarely sensitive, although some quick clays have been found
to be overconsolidated.
3) In-Situ Evaluation of Shear Strength
As mentioned previously, the extraction of a soil sample and the subsequent changes
induced by the extraction process, the handling, and the testing procedures may greatly
alter the characteristics of the specimen and, therefore, the test results. Furthermore, it is
not always feasible and practical to duplicate the in-situ conditions. Frequently, it is both
practical and desirable to test the soil in the in-situ condition.
A number of field tests are used to estimate the shear strength of the soil. Some of
these give results that fit in to the theoretically based expressions used to designate the
shear strength of the soil. Others are empirical in nature and greatly dependent on the
engineers judgment and experience.
34
Faculty of Technology
Is used for the determination of shear strength of soft clays (clays which may be disturbed
during the extraction and testing process). The test is performed at any given depth by
first augering to the prescribed depth, cleaning the bottom of the borings, and then
carefully pushing the vane instrument (Fig 2.22) in to the stratum to be tested. The torque
is then applied gradually and the peak value noted. The shear strength of the soil can then
be estimated by using the formulae derived below.
28 r
26 r
The number of blows required to drive the split-spoon sampler is still another means used
to estimate the soils resistance to shear (shear strength of the soil). Here a split-spoon
sampler (Fig 2.23 (a)) is lowered to the bottom of the bore hole by attaching it to the drill
rod and then driven by forcing it in to the soil by blows from a hammer (64Kg) falling
from a height of 76cm. The sampler is initially driven 15cm below the bottom of the bore
hole to exclude the disturbed soil while boring. It is then further driven 30cm in two
stages (each 15cm). The number of blows required to penetrate the last 30cm is termed as
the standard penetration value, N. Corrections have to be made for overburden pressure,
dilatancy etc. Then corrected N value is used for correlation. HL\6.doc
35
Faculty of Technology
Consistency
Very soft
Soft
Medium
Stiff
Very stiff
Hard
N
0-2
2-4
4-8
8-15
15-30
>30
qu(KPa)
<25
25-50
50-100
100-200
200-400
>400
The correlation for clays is unreliable. Hence, vane shear test is recommended for more
reliable information.
Usually SPT is conducted at every 2m depth or at the change of stratum. If refusal is
noticed at any stage, it should be recorded.
Other Penetrometers
Penetrometers are some times used to test the shear strength of the soil at the surfacelateral or vertical. Their use is primarily applicable to fine-grained soil; coarse and
gravelly strata tend to give erroneous results. The procedure for using the penetrometer
consist of first cleaning the surface of any loose material, pushing the pentrometer in to
the stratum to the calibration mark on the head of the penetrometer, and recording the
maximum reading on the penetrometer scale. This reading represents the pressure in force
per unit area necessary to push the penetrometer to the designated mark. The reliability of
the results must be interpreted with a view to the condition present at the time of testing
(like water content).
36
Faculty of Technology
CHAPTER THREE
LATERAL EARTH PRESSURE AND RETAINING WALLS
1.
Introduction
Soil is neither a solid nor a liquid, but it exhibits some of the characteristics of both. One
of the characteristics similar to that of a liquid is that of its tendency to exert a lateral
pressure (earth pressure) against any object in contact. Thus structures which retain or
support soil like retaining walls, abutments, sheet pile walls, basement walls and under
ground conduits need estimation of the lateral pressure for their design.
Various types of retaining walls are shown in Fig 3.1 and are widely employed in civil
engineering works ranging from their use in road and rail construction to support cuts and
fills where space is limited to prevent the formation of appropriate side slopes, to the
construction of marine structures such as docks, harbours and jetties.
Faculty of Technology
provides the resistance for maintaining stability. The pressure or resistance which soil
develops in response to movement of the structure toward it is called the passive earth
pressure which may be very much greater than the active earth pressure. The surface over
which the sheared-off soil wedge tends to slide is referred to as the surface of sliding or
rupture.
Fig 3.2 depicts the relationship between the earth pressure and the wall movement. Po
represents the magnitude of pressure when no movement of the retaining wall takes place;
it is commonly referred to as earth pressure at rest. As the wall moves toward the backfill,
the pressure increases, reaching a maximum value of Pp (passive earth pressure) at point
C. On the other hand, if the wall moves away from the backfill, the force decreases,
reaching a minimum value of Pa (active earth pressure) at point B.
Fig 3.3 Relationship between earth pressure and wall movement for cohesionless soil
The magnitude of the active and passive forces Pa and Pp could be derived from the basic
condition of static equilibrium as follows.
38
Faculty of Technology
P a = Wtan ( - )
but
W =
1
1
H ( H cot ) = H 2 cot
2
2
1
P a = H 2 cot tan ( - )
2
The maximum value of Pa may be obtained when cPa /c=0. Thus,
Pa 1 2
= H cot sec 2 ( ) + tan ( ) cosec 2 = 0
= cr = 45 +
Pa = H 2 tan 2 45
2
2
1 2
Pa = H K a
where : K a = tan 2 45
2
2
but
W =
(3.1)
1
1
H ( H cot ) = H 2 cot
2
2
1
P p = H 2 cot tan ( + )
2
1
H 2 cot sec 2 ( + ) + tan ( + ) cos ec 2 = 0
= cr = 45
1
Pa = H 2 tan 2 45 +
2
2
Pp =
1
H 2 K p
2
(3.2)
where : K p = tan 2 45 +
2
Ka and Kp are generally referred to as coefficients for active and passive pressure,
respectively. They are constants for any given soil where = constant. It is clear that the
value of Kp is significantly larger than Ka.
3.
Faculty of Technology
The soil deforms vertically under its self weight but is prevented from deforming laterally
because of an infinite extent in all lateral directions. Let E and be the modulus of
elasticity and Poissons ratio of the soil respectively.
Lateral strain, h = h v v + h = 0
E
E
E
h
v
=
= Ko
v 1 v
Ko is the ratio of the intensity of the lateral earth pressure at rest to the vertical stress at a
specified depth.
As the vertical pressure at any depth is v=z, then h=Ko v = Ko z. The distribution of
the earth pressure at rest with depth is linear (hydrostatic nature) for constant properties
such as E,, and , as shown in Fig 3.4 (b)
If a structure such as a retaining wall of height H is interposed from the surface and
imagined to be held without yield, the total thrust on the wall per unit length Po is given
by
Po =
H
h
dz =
. z . dz =
1
K
2
. . H
Faculty of Technology
principal stress. The corresponding Mohr circle for this case is depicted by circle I in Fig
3.5 (b). On the other hand, if the wall were to push against the backfill, a case of passive
pressure would be developed. The vertical stress then becomes the minor principal stress,
and the lateral stress would thus become the major principal stress. The Mohr circle for
this condition is depicted by circle II in Fig 3.5 (b).
Fig 3.5 The orientation of slip planes in granular soil mass with a level surface
Considering circle I (active case), we have
sin =
(
(
1
1
) / 2
(
=
) / 2
(
3
3
1
1
)
)
3
3
ht an
(45
- /2
) (3.3)
( ' 3
'
3
)/ 2
1 ) / 2
(
(
'
3
'
3
)
1 )
1 + sin
= 1 t an
1 sin
(45
+ /2
)=
ht an
(45
+ /2
(3.4)
We note that tan2 (45- /2) and tan2 (45+ /2) in Eqs. 3.3 and 3.4 are constants for
constant values of . Hence, the corresponding pressures against the retaining wall vary
linearly with depth, as indicated by Fig 3.6. The corresponding resultant pressures, active
and passive, can be calculated for a unit length of retaining wall as
1
H 2 tan
2
1
=
H 2 tan
2
Pa =
Pp
1
H 2K
2
(45 + /2 ) = 1 H 2 K
2
(45
/2
)=
( 3 .5 )
( 3 .6 )
41
Faculty of Technology
Fig 3.7: Lateral pressure and slip planes in granular sloped backfill (active case)
Fig 3.7(c) symbolizes an active state of stress. The magnitude of the vertical stress is
depicted by distance OC; the lateral stress, acting parallel to the sloped surface, is
represented by the distance OA. Hence, from Fig 3.7 (c) we have
DE
OD
( + 3 )
OD = 1
2
OB = ODcos
DE = OD sin = AD
sin =
AB =
( 1 3 )
2
BD = ODsin
AD =
(AD) 2 ( BD ) 2 =
[
= OD [cos +
]
]
Ka =
Ka =
[
[
[
[
Thus
] [
] [
]
]
sin sin
2
]
]
( 3 .7 )
h = v K a = z cos K a
For a given slopped surface and uniform soil properties Ka becomes a constant. Thus, the
intensity of load, or stress, varies linearly with depth. Hence, as before the total resultant
active force may be given by
1
Pa = H 2 cos K a where Ka is given in equation 3.7. The direction of the
2
resultant is parallel to the sloped surface.
Now suppose there is inclined backfill and the wall face is also
inclined. To calculate the total lateral force acting on the wall,
combine the lateral force acting on an imaginary vertical surface
passing through the heal of the wall with the weight of the wedge
of soil between the imaginary surface and back of wall to get the
resultant.
[
[
cos +
h
=
v
cos
]
]
cos 2 cos 2
cos 2 cos 2
(3 .8 )
1
H 2 cos K p where Kp is given in equation 3.8. The direction of the
2
resultant is parallel to the sloped surface.
For the case of level surface, equations 3.7 and 3.8 reduce to equations 3.5 and 3.6
respectively.
Pp =
42
Faculty of Technology
Coulombs Theory
In 1776 Coulomb introduced an expression for determining the thrust on retaining walls.
Basic assumptions were
1. The soil is isotropic and homogenous and posses both cohesion and internal friction
2. The failure surface is plane. The back fill surface is planner.
3. Frictional forces are distributed uniformly along the plane failure surface
4. The failure wedge is a rigid body undergoing translation
5. There is a wall friction which develops as the sliding wedge moves along the back of
the wall.
For the active case derivation was made considering Fig 3.8 as follows.
H
H
sin[180 ( + )] =
sin( + )
sin
sin
H sin( + )
AE
AB
using sin law :
=
AE =
sin( + ) sin( )
sin sin( )
1
1 H sin( + ) H
sin( + )
1 H2
*
Area (ABE) = AE * BD =
sin( + )
sin( + ) =
2
2 sin sin( ) sin
2
sin( )
2 sin
2
sin( + )
1H
Weight W = * area (ABE) =
sin( + )
2
sin( )
2 sin
From Fig 3.8 (b) :
P
W
a
=
sin( ) sin[180 - - + + )
1 H2
Wsin( )
sin( + )
sin( )
=
sin( + )
Pa =
sin[180 - - + + )
2
sin( ) sin[180 - - + + )
sin
AB = H/sin ;
BD = ABsin[180 ( + )] =
(3 . 9 )
To determine the orientation of the failure plane ( cr) that produces a maximum Pa, set
Pa
= 0 . After determining cr and substituting it in equation 3.9 we have,
Pa =
H2
2
sin 2 ( + )
(sin )sin( ) 1 +
where K a =
sin( + )sin( - )
sin( )sin( + )
sin 2 ( + )
(sin )sin( ) 1 +
H2
2
Ka
sin( + )sin( - )
sin( )sin( + )
(3 .10 )
The corresponding passive trust can similarly be derived and is given as follows
Pp =
H2
2
sin 2 ( )
(sin )sin( + ) 1
sin( + )sin( + )
sin( + )sin( + )
H2
2
Kp
(3 .11)
43
where K p =
Faculty of Technology
sin ( )
2
(sin 2 )sin( + ) 1
sin( + )sin( + )
sin( + )sin( + )
For a vertical smooth back face and horizontal back fill surface (that is, =90, =0 and
=0) equations 3.10 and 3.11 reduce to equations derived based on the Rankine theory
(3.7 &3.8).
Coulomb arbitrarily placed the resultant thrust Pa (Pp) at the third point from the bottom.
Correspondingly, he assumed the pressure distribution to vary linearly with depth.
Although this assumption appears to give results acceptable for very rigid walls and
granular backfill, it is not valid for relatively flexible bulk heads, for cohesive backfills,
or where the retaining wall rotates about points not close to the bottom.
5. Lateral Earth Pressures In cohesive (C-) soils
The Mohrs circle may be used to determine the lateral thrust on retaining walls
supporting cohesive soils.
Lets assume a vertical and smooth retaining wall and a C- soil backfill with a horizontal
surface. At a given depth z on the element shown in Fig 3.9 (a), the vertical stress 1= z.
The lateral stress is 3. Both of these are principal stresses. The corresponding Mohrs
circle for this case is shown in Fig 3.9 (b).
Sin =
1 sin
cos
2c
3 = 1
+
1
sin
1 + sin
1 - sin 2
(1 - sin ) (1 + sin )
cos
=
=
=
But
1 + sin
1 + sin
(1 + sin ) (1 + sin )
1 - sin
1 + sin
1 sin
1 - sin
2c
3 = 1
+
+ sin
1
sin
1 - sin
= tan 2 45 - . Hence,
But = h and
1
1 + sin
2
3 = h tan 2 45 2ctan 45 = h K 2c K
a
a
2
2
where K
= tan 2 45 -
2
Faculty of Technology
The maximum horizontal stress (or pressure) for the active case occurs when h=H. The
pressure distribution is shown in Fig 3.10(b). The corresponding resultant Pa is
1
P = H 2 K 2cH K
a 2
a
a
(3 .12 )
For the passive state, from Fig 3.9 (c) we note that
( 3 1 ) /2
( 3 1 )
=
(1 + 3 )/2 + c cot
(1 + 3 ) + 2 c cot
Rearrangin g terms,
Sin =
1 + sin
cos
3 = 1
+ 2c 1 sin
1
sin
2
1
sin
(1
- sin ) (1 + sin )
cos
=
=
=
But
1 sin
1 sin
(1 sin ) (1 sin )
1 + sin
1 sin
1 + sin
1 + sin
+ 2c
3 = 1
1 sin
1 sin
1 + sin
= tan 2 45 + . Hence,
But = h and
1
1 sin
2
3 = h tan 2 45 + + 2ctan 45 + = h K + 2c
p
2
2
where K
= tan 2 45 +
2
The maximum horizontal stress (or pressure) for the passive case occurs when h=H. The
pressure distribution is shown in Fig 3.10(c). The corresponding resultant Pp is
1
P = H 2 K + 2cH K
p 2
p
p
6.
(3 .13 )
3 = 2c K
(tension)
This implies the formation of a crack as depicted in Fig 3.11(a). The corresponding
pressure distribution based on eqn. 3.12 is shown in Fig 3.11(b).
Fig 3.11:Tension crack in c- soil and corresponding pressure distribution (active state)
The theoretical depth of crack ht can be determined by recognizing that, at the bottom of
the crack, 3=0. Thus from equation 3.12,
0 = h K 2c K
a
a
2c
h =
t K
a
(3 .14 )
2c K
= h K 2c K
a
a
4c
H =
c K
= 2h
t
(3.15)
45
Faculty of Technology
Though eqn 3.15 provides a theoretical depth to which an excavation may be made with
out lateral support, it should be used cautiously.
HL\30-1.doc30-4
7. Effect of Uniform Surcharge
The extra loading carried by a retaining structure is known as surcharge. If the uniform
surcharge is of intensity q per unit area, then the vertical stress at every elevation in the
backfill is considered to increase by q. As such, the lateral pressure has to increase by
Ka. q. Thus, at any depth z, h= h Ka + Ka. q.
Fig 3.12 (b) and (c) show two different ways in which the pressure distribution may be
shown. In Fig 3.12 (c), the uniform surcharge is also considered to have been converted
into an equivalent height He, of backfill, which is easily established, as shown.
46
Faculty of Technology
Fig 3.15: Rebhanns condition for Coulombs wedge theory- Location of failure plane for
the active case
Let BD be a line inclined at to the horizontal through B, the heel of the wall, D being
the intersection of this -line with the surface of the backfill. The value of Pa depends
upon the angle relating to the location of the failure plane. Pa will be zero when = (as
no wall required to retain the soil mass at angle ), and increases with an increase in up
to a limit, beyond which it decreases and reaches zero again when =180-(here the
failure wedge has no mass). Thus, the failure plane will lie between the -line and the
back of the wall.
Let AE be drawn at an angle (+) to the wall face AB to meet the -line in E. Let CG be
drawn parallel to AE to meet the -line in G. Let the distances be denoted as follows: AE
= a, BG = c, CG = x, BD = b, BE = d.
Now it is required to determine the criterion for which Pa is the maximum, which is supposed to give the correct
location of the failure surface.
Weight of the soil in the sliding wedge is
47
Faculty of Technology
This equation signifies that for BC to be the failure plane the requirement is that the area
of the failure wedge ABC be equal to the area of the triangle BCG (Rebhanns condition).
The triangles ABC and BCG have a common base BC, so altitudes on BC should be
equal. Thus, AJ [sin (<AJB)] = CG [sin (<BCG)]. But sin<AJB = sin<BCG as CG is
parallel to AJ. This leads to CG=AJ=x; and JE= a- x.
Triangles DAE and DCG are similar (AA similarity). Hence b d . x = a
(i)
Also triangles BCG and BJE are similar and consequently
bc
d
.x = a - x
c
(ii)
c 2 = bd
or
c=
bd
Thus if c is known the position of G and hence that of the most dangerous rupture plane
(surface BC) can be determined and the weight of the sliding wedge, W, and the active
thrust, Pa, can be calculated.
Summary
ab
The value of x can now be determined as c x = b(a x)
x=
b+c
=-
Substituting c = b (a x) in the equation of Pa, one gets
c = bd
x
x=bd/ (b+c)
1
(3.16)
Pa = x 2 sin
2
This is a semi-analytical procedure for the computation of the active thrust by Coulombs
wedge theory.
However, elegant graphical methods have been devised and preferred to the analytical
approach, in view of their versatility, coupled with simplicity.
A graphical method is given by Poncelet and sometimes known as the Rebhanns
graphical method, since it is based on Rebhanns condition. The steps involved are as
follows with reference to Fig 3.16.
48
Faculty of Technology
Faculty of Technology
1. With a suitable scale draw the ground line and the retaining wall. This should include
height and slope of the retaining wall, surface configuration of the backfill, location
and magnitude of concentrated (line) surcharge loads, uniformly distributed
surcharge, and so on.
2. From point B draw line BD making an angle with the horizontal (-line)
3. Draw line BK at angle (=) with line (-line)
4. Draw rays BC1, BC2, BC3, and so on, that is , assumed failure surfaces
5. Determine the weight of each wedge, accounting for variations in densities if the
backfill is layered system, for variable moisture content, and so on.
6. To a convenient scale, plot these weights along line BD (-line). For example, the
distance from B to W1 along line BD equals W1; similarly, the distance from W1 to
W2 along line BD equals W2, and so on.
7. From each of the points located on line BD, draw lines parallel to line BK to intersect
the corresponding assumed failure surfaces; that is, the line from W1 will intersect
line BC1 at F1; the one from W2 will intersect line BC2 at F2, and so on.
8. Connect these points of intersection (F1, F2, F3) with a smooth line, Culmanns
curve.
9. Parallel to BD, draw a tangent line to Culmanns curve. In Fig 3.17 point F represents
such a tangent point. More than one tangent is possible if the Culmann line is
irregular.
10. From the point of tangency, draw line FW parallel to line BK. The magnitude of FW,
based on the used scale, represents the active pressure Pa. If several tangents to the
curve are possible, the largest of such values becomes the value of Pa. The failure
surface passes through F and B, as shown in the figure.
b) Passive Earth Resistance -Cohesionless soil
Fig 3.18 illustrates the procedure for determining the passive resistance via Culmanns
method. The approach is similar to that for the active pressure, with some notable
differences: (1) line BD makes an angle of below rather than above the horizontal; (2)
the reference line makes an angle of with line BD, with measured as indicated. For
the assumed sliding wedges, the weights W1, W2, and so on, are plotted along line BD.
From these points, lines are drawn parallel to the -line to intersect the corresponding
rays. The Culmann line represents a smooth curve connecting such points of intersection.
A tangent to the Culmann curve parallel to line BD is drawn, with the resultant earth
pressure being the scaled value of line FW.
Fig 3.18: Culmanns graphical method for passive earth pressure (cohesionless soil)
50
Faculty of Technology
Faculty of Technology
4. Coulombs theory is more versatile than Rankines in that it take into account any
shape of the back fill surface, break in the wall face or in the surface of the fill, effect
of stratification of the backfill, effect of various kinds of surcharge on earth pressure,
and the effects of cohesion, adhesion and wall friction. It leads itself to elegant
graphical solutions and gives more reliable results, especially in the determination of
the passive earth resistance; this in spite of the fact that static equilibrium condition
does not appear to be satisfied in the analysis.
5. Rankines theory is relatively simple and hence is more commonly used, while
Coulombs theory is more rational and versatile although cumbersome at times;
therefore, the use of the latter is called for in important situations or problems.
52
Faculty of Technology
CHAPTER FOUR
STABILITY OF SLOPES
4.1 Introduction
The term slope as used here refers to any earth mass, natural or man-made, whose
surface forms an angle with the horizontal. Hills and mountains, river banks and coastal
formations are common examples of sloped earth masses formed by nature perhaps as a
result of glacial movements, weathering, erosion, deposit buildup and sedimentation, or
other factors. Examples of man-made slopes may include fills such as embankments,
earth dams, levees; examples associated with cuts may include highway and rail way cuts,
canal banks, foundation excavations, and trenches.
Virtually every slope experiences gravitational forces and subsequent changes.
Indeed, the combining effects of gravity and water are primary and direct influences on
changes for most of the slopes (failure of an earth slope). Some slopes may possibly also
be influenced by such natural phenomena as chemical actions, earth quakes, glacial
forces, or wind. In virtually all instances, the effect is a general flattening of the slope,
either suddenly or slowly and cumulatively. Predicting the change with any degree of
accuracy may be a difficult task; preventing such change may be an even greater task for
the soil engineer.
Determination of the potential failure surface and the forces tending to cause slip
and those tending to restore or stabilize the mass of the earth are essential steps in the
stability analysis (functional) of earth slopes and available margin of safety. The soil
mass is assumed to be homogenous. It is also assumed that it is possible to compute the
seepage forces from the flow net and the shearing strength of the soil from the MohrCoulomb theory. The same principles of limiting equilibrium mechanics are used to
analyse the stability of slopes (unretained soil masses) as for retained soil masses. For
purposes of analysis, slopes may be classified either as an infinite or finite. In practice, a
slope may be considered to be infinite where the soil properties at corresponding depths
are the same and where the depth D to a hard stratum is constant and small compared
with the overall length of the slope, as shown, for example, in Fig 4.1(a). For this
geometry, any mass movement or failure of the slope usually involves the displacement
of a soil mass having a length many times greater than its depth, as indicated by the
potential slip path [Fig 4.1(a)]. Neglecting end effects, a typical vertical slice through the
soil of width b can be taken as representative of the whole failure mass, in which case a
consideration of the simple statics governing equilibrium of the typical slice can be used
to derive stability relationship for the slope.
Faculty of Technology
complex than that for infinite slopes. This category includes cuttings for roads, railways,
canals, etc., and embankments for roads and earth dams.
4.2 Infinite slopes
Lets consider an infinite slope inclined at angle to the horizontal in a general C- soil
as shown in Fig 4.2.
Normal load
Wcos Wcos W
N
=
=
=
=
cos 2 = Hcos 2
Area
l*l
l*l
b/cos
b
(4.2)
Shear force
Wsin Wsin W
T
=
=
=
=
sin cos = Hsin cos
Area
l*l
l*l
b/cos
b
The shear strength of the soil is, = c' + tan' = c' + ( Hcos 2 - u)tan '
Factor of safety against slippage, F =
(4.1)
(4.3)
(4.4)
Shear strength
c' + ( Hcos - u)tan '
=
shear stress
Hsin cos
2
At
( 4 .5 )
the verge of failure the shear strength is equal to the acting shear stress, i.e. F=1. If we
make F=1 in equation 4.5, we can find that the slope will be stable to a maximum depth
Hc, called the critical depth of the stratum. Thus
Hc =
( 4 .6 )
From equations (4.5) and (4.6) different cases can be seen as:
Case (i) Cohesionless soil - No seepage
c' + ( Hcos 2 - u)tan ' tan '
=
Hsin cos
tan
F=
Hc =
=0
sin cos - cos 2 tan' = 0 tan = tan ' cr = ' is called the critical slope
Hc =
54
Faculty of Technology
Hc =
sec 2
t an - tan '
= ABcos = Hcos2
Thus, u = h = Hcos2
w w
w
( 4 .7 )
( 4 .8 )
The critical depth Hc obtained from eqn. 4.8 when F=1 is,
Hc =
c'
( 4 .9 )
c'
Example: A relatively cohesive soil at a constant infinite slope has a negligible seepage
and negligible pore water pressure. If =18 KN/m3, c' = 36Kpa, =140, H=3m, and
=220, find the factor of safety against sliding.
2
2
0
14 0
Solution: From eqn 4.7 F = c' + H( - w )cos tan ' = 36 + 3 * (18 0)(cos0 22 )tan
= 2.54
0
Hsin cos
18 * 3 * sin 22 cos 22
Faculty of Technology
Investigation of the stability of finite slopes involves the following steps according to the
commonly adopted procedure:
a) Assuming a possible slip surface,
b) Studying the equilibrium of forces acting on this surface (activating & resisting)
c) Repeating the process until the worst slip surface, that is, the one with minimum
margin of safety is found.
Methods of analysis can be by:
1) Considering the whole sliding mass as a free body which include,
Culmann method
The Circular arc analysis
2) Dividing the sliding mass in to a number of vertical slices and considering the
equilibrium of each slice
Swedish method of slices- (Fellenius and Peterson)
Bishop Method
Taylors method
4.3.1 Culmann Method
This method of analysis is based on the assumption that failure occurs on a slip plane
through or above the toe of the slope. Fig 4.4 shows the forces acting on the sliding mass.
1
hL * 1
2
Hsin( )
H
h
=
h =
sin
sin( )
sin
sin( )
1
W = HL
2
sin
The resisting force to sliding along AC
F R = cL + Wcos tan
AB =
c+
sin( )
1
H
cos tan
2
sin
sin( )
1
H
sin
2
sin
A value of FS ranging from 1.25 to 1.5 is considered to be satisfactory for an earth slope. For
economic reasons, a value greater than 1.5 is not desired. Hence, FS=1.5 may be considered
to be necessary as well as sufficient.
Referring to the force polygon shown in Fig 4.4(c)
1/2 LHsin( - )
1 sin( - )sin( - )
c
W
cL
=
=
=
( 4 .11) w
sin( - )
cos
sin cos
sin cos
here the quantity c/H is known as the stability number. The condition for impending slip
occurs when c/H is a maximum. Thus, from eqn 4.11,
56
Faculty of Technology
(c / H )
1
= 0 cr = ( + )
( 4 .12 )
The method gives good results for very steep or vertical slopes; it does not provide
satisfactory results for relatively flat slopes.
Example: Find the factor of safety of the slope if =17KN/m3, c=26KPa, =180, =480
and H=15m.
Solution: The critical failure plane orientation from eqn 4.12 is cr = 1 ( + ) = 33 0 and
2
factor
of
safety
from
eqn
4.10
is
then
sin( )
1
sin(48 3 3 )
1
cos tan
c + H
26 + * 15 * 17
cos 3 3tan 1 8
2
sin
2
sin48
=
FS =
= 1 . 58
sin( )
sin(48 3 3 )
1
1
H
sin
* 15 * 17
sin 3 3
2
sin
2
sin48
r c r + tan ' i
i =1
FS =
We
(4.13)
If there is a tension crack, the cracked part will not have shear resistance and hence the
total cohesion resistance will be c(AB'). And AB' =r ' where ' is in radian
If the backfill is purely cohesive =0 [or if analysis is on total stresses (=0 analysis)],
57
Faculty of Technology
c r 2
We
In eqn 4.13 the term 'i varies over the arc and is not symmetrical. Its value is influenced
mainly by the weight of the soil above a point i. This approach leads to the method of
slices discussed below.
4.3.3 The Swedish Method of Slices
In this method of analysis the area within the slip circle is divided in to a number of
vertical slices generally of equal width [Fig 4.6 (a)]. The forces on each slice are
evaluated from the limiting equilibrium of slices. These forces are shown in Fig 4.6 (b).
The equilibrium of the entire mass is determined by summation of the forces on each slice
as forces activating failure and forces resisting failure.
Here
Normal stress, i =
T i = W i sin i
and
Ni
W i cos i
=
Li *1
Li
Shear stress, i =
and
Ti
W i sin i
=
Li *1
Li
Wi
cos i tan
Li
L
i =1
cL
i =1
+ W i cos i tan
i =1
W sin
i =1
cL
i =1
+ W i cos i tan
i =1
W sin
i =1
( 4 . 14 )
activating moments. The shearing resistance with in the tension cracked zone is zero.
Ai (m2)
Wi (KN)
i (deg)
Li (m)
CLi (KN)
Wi cosi (KN)
Wi sini(KN)
1
2
etc
Sum
In this method a number of slip circles are analyzed and the minimum factor of safety is
finally obtained. The circle with this minimum factor of safety is the failure circle.
58
Faculty of Technology
W i cos i
ui
Li
n
FS =
c' L + (W cos
i =1
i =
and
i =1
W i sin i
Li
u i L i )(tan ' )
( 4 . 15 )
W sin
i
i =1
The effect of pore water pressure is to reduce the effective stress and thereby reduce the
stability because the shear strength mobilized would be decreased.
The factor of safety is calculated from,
n
FS =
i =1
(4.16)
W sin
i =1
In the absence of a flow net, the approximate value of FS may be obtained from:
n
FS =
i =1
c' L i + (W i cos i
i =1
n
i =1
W i sin
)(tan ' )
(4.17)
Here the normal component Ni (=Wi cosi) of the weights of the slices have to obtained
using effective or buoyant unit weight sub (b or ) and the tangential components Ti
(=Wisini) using the saturated unit weight.
59
Faculty of Technology
FS =
c' L + (W cos
i =1
i =1
u i L i )(tan ' )
W sin
i =1
FS =
c' L
i =1
(4.18)
i =1
i =1
W i sin i
In the absence of a flow net, the approximate value of FS may be obtained from eqn 4.17
given previously.
60
Faculty of Technology
Earth quake
In the case of an earth quake, a horizontal seismic force is applied at the centroid of the
slice as shown in Fig 4.9 and the factor of safety is estimated as usual.
n
FS =
W sin + C W a
i =1
(4.19)
i =1
i =1
i i
N i = Wi cos i u i li
and
(4.20)
(4.21)
Substituting eqn 4.20 to eqn 4.21 and dividing by cos i and rearranging, we have
c ' li
sin i + (F1i F2i )
FS
Ni =
sin i tan '
cos i +
FS
c' l +
i
FS =
c' li
sin i + (F1i F2i )
FS
tan '
sin i tan '
cos i +
FS
Wisini
Wi ui li cos i
61
FS =
Calculatio
Faculty of Technology
c' l i cos
[W
u i l i cos
+ (F1i F 2i
)]tan '
W i sin i cos i +
FS
l be
c' l i cos
[W
W i sin i cos
u i l i cos
]tan '
( 4 . 22 )
FS
The solution of this equation requires an iterative analysis since FS appears on both sides
of the equation. Trial and error approach consists of assuming an FS value on the right
hand side of the equation, then solving for FS on the left. If the difference between the
assumed FS and the computed FS is significant, a new FS is assumed and the procedure
repeated until a satisfactory FS is determined. Computer programs are generally preferred
in this method. If '=0, the eqn 4.22 may be solved directly for FS.
The following table summarizes different situations and their corresponding preferred
methods.
S. No.
Situation
Preferred Method
1
End of construction with partly Total stress analysis (=0 analysis)by
saturated soil-construction period is using Cu and u from UU tests
short compared to consolidation time
2
Steady Seepage
Effective stress analysis by using
C', ' with the pore water pressure
obtained from flow nets
3
Rapid draw down
Effective stress analysis with C', '
from CU tests. Pore water pressure
estimated from flow nets or other
methods
4
Earth quake
Total or Effective stress analysis
depending on the data available
4.3.5 Friction Circle Method
The reaction R against the failure surface makes an angle to the normal to the failure
surface. So, all reactions will be tangent to the circle with O as center and rsin as a
radius (Fig 4.11). This circle is called the friction circle.
Faculty of Technology
In the friction circle method with known , the following quantities are known:
1. Magnitude of weight sliding wedge W and its line of action
2. Direction of R (at angle to the normal to the failure surface). Magnitude of R and its
distribution along the failure surface not known
3. The direction of total cohesion Cm=cm L along the rupture surface parallel to chord AB
Lets assume a uniform distribution of cohesion along the failure surface. If we resolve
the cohesion force parallel and perpendicular to the chord AB, the components
perpendicular to chord are oppositely directed, they do not contribute any moments. If cm
is mobilized cohesion, then
c m L r = c m L' r c
rc =
Lr
L'
The line of action of W and Cm are located first. Then the triangle of forces may be
completed as shown in Fig 4.11(b) and the values of R and Cm can then be determined
and the factor of safety is calculated.
The factor of safety with respect to cohesion, assuming the friction is mobilized in full, is
given by:
FS
c
cm
(4.23)
The factor of safety with respect to friction, assuming the cohesion is mobilized in full, is
given by:
FS
tan
tan m
(4.24)
Where and m are the total (effective) angle and mobilized friction angle. If the factor of
safety with respect to the total shear strength FS is required, m is to be chosen such that
FSc and FS are equal. Here FS is:
FS =
shear strength ( = )
mobilized shear resistance ( =
(4.25)
Where: = c + 'tan
and
m = cm + ' tanm
4.3.6 Taylors Method
Taylor (1948) prepared two charts relating the stability number (N=c/H) to the angle of
slope, , based on the friction circle method and analytical approach. The first is for the
general case of a c- soil with the angle of slope less than 530 (Fig 4.12). The second is
for a soil with =0 and a layer of rock or stiff material at a depth DH below the top of the
embankment, as shown in Fig 4.13. Here, D is known as the depth factor; depending
upon its value, the slip circle will pass through the toe or will emerge at a distance nH in
front of the toe (the value of n may be obtained from the curves).
Fig 4.12: Taylors charts for slope stability (After Taylor, 1948)
(for =00 and <530, use Fig 4.13)
63
Faculty of Technology
Fig 4.13: Taylors charts for slope stability (After Taylor, 1948)
Taylors charts are based on the assumption of full mobilization of friction, that is, these
give the factor of safety with respect to cohesion. This is all right for a purely cohesive
soil; but, in the case of c- soil, where the factor of safety FS with respect to shearing
strength is desired, m should be used for :
tan m =
(Also
tan
FS
FS
(4.26)
)
The charts are not applicable for a purely frictional soil (c = 0). The stability then depends
only upon the slope angle, irrespective of the height of the slope.
64
Faculty of Technology
CHAPTER FIVE
BEARING CAPACITY OF SOILS
1. Introduction
All civil engineering structures impose a loading on the underlying soil or rock. The part
of the structure, usually lying below ground level which transmits the load to the
supporting strata is referred to as the foundation. To ensure stability, foundations must
provide an adequate factor of safety against shear or bearing failure of the underlying
soil and the structure must be capable of withstanding the settlements that will result, in
particular the differential settlements. In the vast majority of cases the choice of the
foundation is governed by the second of theses two factors, thus in general the need to
limit settlements of the structure will control the design of its foundation.
Thus the criteria for the determination of the bearing capacity of a foundation are based
on the requirements for the stability of the foundation. The design value of the safe
bearing capacity would be the smaller of the two values, obtained from the two criteria:
1.
Shear strength criterion
HL\51.doc
2. Settlement criterion
The soils limiting shear resistance is referred to as the ultimate bearing capacity, qu, of
the soil. For design, one uses an allowable bearing capacity, qa, obtained by dividing the
ultimate bearing capacity by a suitable safety factor (i.e. qa=qu/FS). HL\52.doc
The bearing capacity is affected by factors like
a) Nature of the soil and its physical and engineering properties
b) Size, shape, depth, rigidity and roughness of the foundation
c) Water table conditions and initial stresses in the foundation soil
d) Total and differential settlements that the structure can withstand without
functional failure
In general, foundations constructed at a depth below the ground level which is
approximately less than or equal to 2 (i.e. Df <2B) like footings and mat foundations are
designated as shallow foundations. In contrast, piles, piers and caissons (Df >2B) are
designated as deep foundations.
Numerous proposals have been advanced regarding considerations, criteria, and
procedures for the evaluation of the bearing capacity of soils. Some analytical methods
of estimating bearing capacity are given below.
Faculty of Technology
directions. The P associated with wedge I represents the active pressure resultant,
whereas the P for wedge II is the passive thrust.
Pa =
(a)
Pp =
(b)
The two resultants at the interface are assumed to have same magnitude. Hence equating
eqn (a) and (b), we have
1
K a + qu K
K a H 2 2 cH
2
solving for q u , we have
qu =
1
2c
(K p K a ) +
K
K
K
a
a
H =
- /2)
2tan(45
+ 2cH
)+
DfK
+ D fK pH
(c )
B
2
1
B K
4
Let
Thus
1
K
2
H =
1
H
2
But K
Hence,
) + 2c (K + K
1
(K K )
=
2
)
= 2 (K
+ K
= K
DfK
(d)
and
)+
1 + sin
= tan
1 sin
(45
+ /2
= cN
1
BN
2
+ DfN
( 5.1)
For a purely cohesive soil (=0): Kp =1, N =0, Nc =4, Nq =1 and hence the ultimate
bearing capacity from eqn 5.1 is qu = 4c + Df . So in purely cohesive soils increase or
decrease in the width, B, of the footing has no effect on qu. But increasing Df increases qu
by a little factor.
The above derivation is base on less-than-accurate assumptions: (1) the shear at the
interface of the two wedges was neglected and (2) the failure surfaces are not straight
lines as was indicated for the two wedges.
3. Prandtls Theory for Ultimate Bearing Capacity
Prandtls theory of plastic equilibrium reflects on the penetration (deformation) effects of
hard objects in to much softer material. The assumptions made were: (1) Soil is isotropic
and homogenous (2) Soil is weightless (3) Footing is long and has a smooth base and
66
Faculty of Technology
placed at the surface of the soil. Fig 5.2 shows three zones developed with in the soil
resulting from bearing failure.
Zone 1: Active zone
Zone 2: zone of radial shear or
plastic zone
Zone 3: passive zone
(5.2)
(5.3)
Fy = 0
q u B = 2 P p + 2(bd)csin
But bd = (B/2)/cos
. Thus,
67
Faculty of Technology
q u B = 2 Pp + Bctan
(5.4)
Terzaghi represented the value of Pp as the vector sum of the three components: (1) that
from cohesion, (2) that from surcharge, and (3) that resulting from the weight of the soil
(bdef in Fig 5.3 (b)). With the addition for the shape factors in the cohesion and base
terms, Terzaghi obtained expressions for the ultimate bearing capacity for general shear
conditions as:
Long
footings
Square
footings
Circular
where
:
:
footings
:
:
=
= cN
= 1 . 3 cN
= 1 . 3 cN
2 cos
a 2
(45 + / 2
+ D
+ D
+ D
;
K p
1
tan
1
2
2
cos
+ 33
2
45 +
= 3tan
1
B N
2
+ 0 . 4 B N
( 5.6)
+ 0 . 3 B N
( 5.7)
q
f
N
c
[N
( 5.5)
1 cot
with
a = e
(After
3
4
2
tan
S. Husain)
Table 5.1 below gives the values for the various bearing capacity factors recommended
for the above equations.
Table 5.1: Terzaghis N-factors
0
Nc
Nq
2
6.3
1.22
0.18
4
6.97
1.49
0.38
6
7.73
1.81
0.62
8
8.6
2.21
0.91
10
9.61
2.69
1.25
12
10.76
3.29
1.7
14
12.11
4.02
2.23
16
13.68
4.92
2.94
26
28
0
Nc 27.09 31.61
Nq 14.21 17.81
N 11.35 15.15
30
37.16
22.46
19.73
32
44.04
28.52
27.49
34
52.64
36.51
36.96
36
63.53
47.16
51.7
38
77.5
61.55
73.47
40
95.67
81.27
100.39
42
44
46
48
50
119.67 151.95 196.2 258.29 347.52
108.75 147.74 204.2 287.86 415.16
165.69 248.29 427 742.61 1153.2
0
5.7
1
0
18
15.52
6.04
3.87
20
17.69
7.44
4.97
22
20.27
9.19
6.61
24
23.36
11.4
8.58
The results obtained here are quite within acceptable limits for shallow footings (e.g.
Df/B<1) subjected to only vertical loads. But they are limited to concentrically loaded
horizontal footings; they are not suitable for footings that support eccentrically-loaded
columns or to tilted footings. Furthermore, they are regarded as somewhat overly
conservative.
Terzaghi developed his bearing-capacity equations assuming a general shear failure in a
dense soil and a local shear failure for a loose soil. For the local shear failure he proposed
reducing the cohesion and as:
c ''=
2
c
3
tan
1
InclinedLoad: q u = cNc s c d c ic + Df N q s q d q iq + B N s d i
2
2
tan
where: N q = e
tan (45 + / 2)
N c = (N q 1)cot
(5.8)
N = (N q 1) tan(1.4 )
Faculty of Technology
Nc
Nq
N
0
5.1
1
0
2
5.63
1.2
0.01
4
6.19
1.43
0.04
0
26
28
30
6
6.81
1.72
0.11
8
7.53
2.06
0.21
10
8.34
2.47
0.37
12
9.28
2.97
0.6
14
10.37
3.59
0.92
16
11.63
4.34
1.37
18
13.1
5.26
2
20
14.83
6.4
2.87
32
35.49
23.18
22.02
34
42.16
29.44
31.15
36
50.59
37.75
44.43
38
61.35
48.93
64.08
40
75.32
64.2
93.69
42
44
46
48
50
93.71 118.37 152.1 199.27 266.89
85.38 115.31 158.51 222.31 319.07
139.32 211.41 328.74 526.47 873.89
D
B
Any
d c = 1 + 0.2 K p
s c = 1 + 0.2K p
B
L
0
d q = d = 1.0
s q = s = 1.0
For =0
D
For >100 s = s = 1 + 0.1K B
d q = d = 1 + 0.1 K p
q
p
B
L
22
16.88
7.82
4.07
24
19.32
9.6
5.72
Inclination
ic = iq = 1 0
90
i = 1 .0
i = 1 -
V Q
H
B
When triaxial tr is used for plain strain, adjust tr to obtain ps = 1.1 0.1 tr
L
Meyerhof suggested that footing dimensions B'=B-2ey and L'= L-2ex be used in
determining the total allowable load eccentrically applied in the x and y directions,
respectively (i.e., Qu=qu B' L'), and in the corresponding terms in the ultimate bearing
capacity equations and in the various correction factors for shape and inclination.
6. Hansens Bearing Capacity Equation
Hansen proposed the general bearing capacity equation (eqn 5.9) which includes ground
factors and base factors to include conditions for a footing on a slope.
1
q u = cN c s c d c i c b c g c + D f N q s q d q i q b q g q + B N s d i b g
(5.9)
2
where : N q = e tan tan 2 (45 + / 2 ); N c = (N q 1)cot ; N = 1 .5 (N q 1) tan E
xpressions for inclination, shape, depth, base, and ground inclination expressions
proposed by Hanson are given in Table 5.3 (b).
Table 5.3 (a): Hansens N- factors
69
0
5.1
1
0
Faculty of Technology
2
5.63
1.2
0.01
4
6.19
1.43
0.05
6
6.81
1.72
0.11
8
7.53
2.06
0.22
10
8.34
2.47
0.39
12
9.28
2.97
0.63
14
10.37
3.59
0.97
16
11.63
4.34
1.43
26
28
0
Nc 22.25 25.8
Nq 11.85 14.72
N 7.94 10.94
30
30.14
18.4
15.07
32
35.49
23.18
20.79
34
42.16
29.44
28.77
36
50.59
37.75
40.05
38
61.35
48.93
56.18
40
75.32
64.2
79.54
42
44
46
48
50
93.71 118.37 152.1 199.27 266.89
85.38 115.31 158.51 222.31 319.07
113.96 165.58 244.65 368.68 568.59
>0)
B'
sin
L'
B'
s = 1.0 - 0.4
0 .6
L'
s q = 1 .0 +
ic =
d c = 1 . 0 + 0 .4 k (for >0)
d q = 1 + 2 tan (1 sin ) k
2
d = 1.0
D
1
B
D
k = tan 1
if
B
k in radians
k =
D
B
20
14.83
6.4
2.95
22
16.88
7.82
4.13
24
19.32
9.6
5.75
Inclination factors
L'
sc = 1 +
18
13.1
5.26
2.08
if
D
> 1
B
1 1
2 2
1 iq
ic = iq
Nq 1
Hi
(for =0)
AfCa
(for >0)
0 .5 H i
i q = 1
V + A f C a cot
2 1 5
0 .7
H i
0
450
i = 1
V + A f C a cot
22 5
g c = 1 .0
(for >0)
147 0
5
g q = g = (1 0 . 5 t an )
bc = 1
0
147
(for >0)
b q = e 2 tan
b = e 2 .7 tan
in radians
in radians
Notes
9
9
Failure can take place either along the long side or along the short side and thus
shape , depth and inclination factors shall be calculated in both sides
Use Hi as either HB or HL for inclination factors
7. Vesics Bearing Capacity Equation
70
Faculty of Technology
The Vesic procedure is essentially the same as the method of Hansen with select changes.
The Nc and Nq terms are those of Hansen but N is slightly different as is given by:
N = 2(N q + 1) tan (also see Table 5.4 (a))
There are also differences in the ii, bi and gi, terms (Table 5.4 (a)).
Table 5.4 (a): Vesics N - factors
0
5
10
15
20
25
26
0
0
0.4
1.2
2.6
5.4
10.9
12.5
N
30
32
34
36
38
40
45
0
41
56.2
77.9 109.3 271.3
N 22.4 30.2
sc
B
s q = 1 . 0 + tan
L
B
s = 1.0 - 0.4
0 .6
L
d c = 1 . 0 + 0 .4 k (for >0)
d q = 1 + 2 tan (1 sin ) 2 k
28
16.7
50
761.3
Inclination factors
mH i
(for =0)
ic = 1
AfCaNc
Where:
D
B
D
1
B
D
k = tan 1
if
B
k in radians
if
D
> 1
B
2 + B/L
1 + B/L
2 + L/B
=
1 + L/B
m = mB =
m = mL
d = 1.0
k =
HL\53.doc
1 iq
ic = iq
Nq 1
(for >0)
Hi
i q = 1
V + A f C a cot
Hi
i = 1
V + A f C a cot
m +1
gc =
in rad (for =0)
5 . 14
g c = iq
0
1 iq
(for >0)
5 . 14 tan
g q = g = (1 . 0 t an )
2
5 . 14 tan
bq = b = (1 tan )
(for >0)
2
in radians
Notes
9
Compute m=mB when Hi=HB (H parallel to B) and m=mL when Hi=HL (H // L). If
you have both HB and HL use m = mB2 + mL2 . Note use of B and L, not B',L'.
Always iq, i > 0. For Vesic use B' in the N term even when Hi=HL
71
Faculty of Technology
Example: Compute the allowable bearing capacity via Terzaghis, Meyerhofs, and
Hansens equations for the footing and soil parameters shown in Fig (Ex1). Use a safety
factor of 2.5 to obtain qall.
Solution:
1. Terzaghi:
q
= cN
+ D
F or = 24 , N
all
= cN
=
1
B N
2
a 2
=
= 11 . 4 ;
2
2 cos (45 + / 2 )
N
K
1
tan
2
cos
= 3tan
+ D
p
2
1 = 8 .9
with
[N
1 ] cot = 23 . 4 ;
a = e
tan
2
4
= 2 . 60
+ 33
= 34 . 19
2
45 +
1
B N
2
= 662.2KPa
q u
= 264 . 88 KPa
FS
2. Meyerhofs equation:
1
B N s d i
2
tan 2 (45 + / 2 ) = 9 .60 ,
q u = cN c s c d c ic + D f N q s q d q iq +
For = 24 : N q = e tan
N = (N q 1) tan( 1 . 4 ) = 5 .72
K p = tan 2 45 + = 2 .37
2
B
Hence, s c = 1 + 0.2K p = 1.02
L
D
d c = 1 + 0.2 K p = 1.31
B
i c = i q = i = 1.0
B
s q = s = 1 + 0.1K p = 1.01
L
D
d q = d = 1 + 0.1 K p = 1.15
B
Thus : q u = cN c s c d c ic + D f N q s q d q iq +
q all =
N c = (N q 1) cot = 19 .32 ,
1
B N s d i = 678.81Kpa
2
qu
= 271 .52 KPa
FS
3. Hansens equation
72
Faculty of Technology
1
B N s d i b g
2
tan 2 (45 + / 2 ) = 9.60 ; N c = (N q 1)cot = 19 .32 ;
q u = cN c s c d c i c b c g c + D f N q s q d q i q b q g q +
For = 24 : N q = e tan
B
B
B
= 1.02 s q = 1 + sin = 1.02 s = 1 0.4 = 0.98
L
L
L
D
as D/B < 1, k = D/B d c = 1 + 0.4 k = 1.40 d q = 1 + 2 tan (1 sin ) 2 = 1.31
B
d = 1.0; i c = i q = i = 1.0; b c = b q = b = 1.0; g c = g q = g = 1.0
Nq
Hence, s c = 1 +
Nc
q u = cN c s c d c i c b c g c + D f N q s q d q i q b q g q +
q all =
1
B N s d i b g = 729 .1Kpa
2
qu
= 291 .64 KPa
FS
73
Faculty of Technology
dw
H
'
(H d w )2
2
H
Where: H=0.5B tan (45+/2); dw =depth of water table below base of footing
and (=-w) are wet and submerged unit weight of the soil respectively
4. The Water Table Below the wedge zone
When the water table is below the wedge zone [depth approximately H=0.5Btan(45+/2)
from base of footing], the water table effects can be ignored for computing the bearing
capacity.
9. Gross and Net Bearing Capacities
The bearing capacity obtained in the above equations is the gross bearing capacity. The
net bearing capacity (qnet) is obtained by deducting the original overburden pressure or
surcharge pressure from the gross bearing capacity (qult).
q
= cN
+ D fN
q net = q u - D f
1
B N
2
q net = cN c + D f N q + B N D f
2
1
q net = cN c + D f (N q 1) + B N
2
Floating foundation: If a foundation on the soil exerts a pressure of intensity less or equal
to the overburden pressure removed during excavation, there will be no need of checking
bearing capacity and such a foundation is called floating (compensating) foundation. For
example in a clay soil (=20KN/m3), a building producing a pressure of 80Kg/m2 with a
basement placed at 5m depth has a floating foundation. (H=100 KPa >80KPa).
10. Bearing Capacity Based on Tolerable Settlement
As discussed previously, the bearing capacity of a foundation is based on two criteria-the
pressure that might cause shear failure of the foundation soil and the maximum allowable
pressure such that the settlements produced are not more than the tolerable values. The
first criterion has already been discussed in detail. For the second criterion, the tolerable
values of the total and differential settlements which a particular structure, on a particular
74
Faculty of Technology
type of foundation in a given soil, can undergo without sustaining any harmful effects are
to be decided up on. These values have already been specified, basing on experience and
judgement.HL\54.doc Once the limiting values of settlement are fixed, the procedure
involves determining that pressure which causes settlements just equal to the limiting
value. This is allowable bearing capacity on the basis of the settlement criterion. It is to
be noted that there is no need to apply a further factor of safety to this pressure, since it
would have been applied even at the stage of fixing up tolerable settlement values.
The smaller pressure of the values obtained from the two criteria is termed the allowable
bearing pressure, which is used for design of the foundation.
The bearing capacity based on settlement criterion may be determined from the field load
tests or plate load tests, standard penetration tests or from the charts like those prepared
by Terzaghi and Peck, based on extensive investigation.
i. Bearing Capacity From SPT
The SPT is widely used to obtain the bearing capacity of soils directly.
9 Terzaghi and Peck (1967) made the first suggestion
9 Meyerhof come up (1974) with modification
9 Bowels (1995) made another modification which he claims to yield an average of
50% increase in allowable bearing capacity. This equation is given as:
q all (KPa)
= 1 2.5 N
q all (KPa)
where
ii.
9
55
B + 0.3
55
B 1.2m
K d;
Df
1 . 33
3B
= corrected N for E
B 1.2m
=1+
value)
55
all
S(mm)
25mm
(q all )25mm
qc
;
30
q all (KPa)
q c B + 0 .3
;
50
B
B 1.2m
for
(a)
for
B > 1.2m
(b)
iii.
q c in KPa
75
Faculty of Technology
Obviously the most reliable method of obtaining the ultimate bearing capacity and the
settlement characteristics at a site is to perform a load test. The test is also used in the
design of highways and runways. The probable settlement of the soil for a given loading
and at a given depth can also be determined.
Round plate with standard diameter (30cm and 70cm) or square plate of side (0.3m x
0.3m and 60cmx60cm) is loaded in a pit excavated in the ground, at a depth equal to the
roughly estimated depth of the foundation for which the bearing capacity is to be
estimated. The procedure is:
9 Excavate a pit to the depth the test is to be performed. The test pit should be at least 4B
(or 4R) wide as the plate to the depth the foundation is to be placed.
9 A load is applied on the plate by increments (P=qult, estimated/5), and settlements are
recorded from dial gauges (at least 3 in no.) for each load increment. Then plots of
settlement vs. time and settlement vs. applied stress are made.
9 The test is continued until a total settlement reaches 25mm or until the capacity of the
testing apparatus is reached or until the soils fails by shear(plate starts to sink rapidly)
Fig 5.5 presents the essential features of the test and typical plots obtained from the plate
load test.
b (B + 0.3)
SP = SF P
for sands and
B(b p + 0.3)
SP =
bP
SF for clays
B
where: B= width of footing (least dimension) and bp= width (diameter) of plate
The permissible settlement value, such as 25mm, should be substituted for SF in the
above equations and the SP value will be calculated. Then from the load-settlement curve,
the pressure corresponding to the computed settlement SP, is the required value of the
ultimate bearing capacity, qult, for the footing.
The coefficient of subgrade reaction, ks, can also be estimated as:
76
ks =
Faculty of Technology
0.4 max
=
KN/m3
S
S
Compactness**
or
Consistency***
Presumed Design
Bearing
Resistance (KPa)
5600
2800
2800
Soft
1400
Soft shale
Soft
850
Description
Massively crystalline igneous and
metamorphic rock (granite, basalt,
gneiss)
Foliated metamorphic rock (slate,
schist)
Sedimentary rock (hard shale,
siltstone, sandstone, limestone)
Remarks
77
Faculty of Technology
Silt
Cohesive soils
Clay
Dense
560
Medium dense
420
Loose
280
Dense
420
Medium dense
280
Loose
140
Hard
280
Stiff
Medium stiff
Soft
200
140
70
Hard
Stiff
Medium stiff
Soft
Very soft
420
280
140
70
Not Applicable
Width of foundation
(B) not less than 1 m
Ground water level
assumed to be depth
not less than B below
the base of the
foundation
* The given design bearing values do not include the effect of the depth of embedment of the
foundation
** Compactness:
*** Consistency:
78