Agricultural Engineering Soil Mechanics
Agricultural Engineering Soil Mechanics
6. Housing of Animals
by A. Maton, J. Daelemans and J. Lambrecht
1985 xii + 4 5 8 pp.
Agricultural
Engineering
Soil Mechanics
E. M c K Y E S
ELSEVIER
Amsterdam — Oxford — New York — Tokyo 1989
ELSEVIER SCIENCE PUBLISHERS Β. V.
Sara Burgerhartstraat 25
P.O. Box 211, 1000 AE Amsterdam, The Netherlands
ISBN 0-444-88080-1
All rights reserved. No part of this publication may be reproduced, stored in a retrieval system or
transmitted in any form or by any means, electronic, mechanical, photocopying, recording or
otherwise, without the prior written permission of the publisher, Elsevier Science Publishers B.V./
Physical Sciences & Engineering Division, P.O. Box 330, 1000 AH Amsterdam, The Netherlands.
Special regulations for readers in the USA - This publication has been registered with the Copyright
Clearance Center Inc. (CCC), Salem, Massachusetts. Information can be obtained from the CCC
about conditions under which photocopies of parts of this publication may be made in the USA. All
other copyright questions, including photocopying outside of the USA, should be referred to the
publisher.
No responsibility is assumed by the Publisher for any injury and/or damage to persons or property
as a matter of products liability, negligence or otherwise, or from any use or operation of any meth
ods, products, instructions or ideas contained in the material herein.
FOREWORD
LIST OF SYMBOLS
δ- A n g l e of s u r f a c e f r i c t i o n - P e r m e a b i l i t y , m/s
€- A n g l e of a plane - T h e r m a l c o n d u c t i v i t y , J/(ms°)
- Engineering strain Κ - E r o d i b i l t y f a c t o r , ton/acre
7]- Porosity Ka - A c t i v e pressure c o e f f i c i e n t
- N e g a t i v e shear slip line K -
c Cohesion coefficient
β- A n g l e of a plane KQ - A t rest c o e f f i c i e n t
μ - 7Γ/4 - φ/2 Kp - P a s s i v e pressure c o e f f i c i e n t
£- P o s i t i v e shear slip line L - Length, m
p- M a s s density, t / m 3
Ls - L a t e n t heat of f u s i o n , J/kg
- Angle LL - Liquid limit, %
(7- Stress level, kPa m - Slope stability number
(7 - D Blade stress, k P a m - v C o e f f i c i e n t of c o m p r e s s i o n , k P a
tfjw' Ice-water energy, J / m 2
Μ - Moment, kNm
σ - η N o r m a l stress, k P a η - Exponent
(7'- Effective stress, kPa Πςΐ - N u m b e r of potential drops
7 - Shear s t r e s s , k P a nf - N u m b e r of flow paths
φ- c cot0 Ν - Normal force, k N
<jf>- A n g l e of internal f r i c t i o n - N u m b e r of passes
a - Length, m Ny- Frictional cutting factor
A - Area, m 2
N - c Cohesional cutting factor
- Soil loss, t/ha N - c a Adhesion cutting factor
b - Width, m Nq - Overburden cutting factor
Β - Foundation width, m ρ - Pressure, kPa
c - Cohesion, kPa Ρ - Wall or blade f o r c e , k N
c -
a Adhesion, kPa P -
!
a E f f e c t i v e wall f o r c e , k N
C - Undrained strength, kPa Pa - A c t i v e wall f o r c e , k N
- Cropping factor PI - P l a s t i c i t y index, %
C -
c C o e f f i c i e n t of c o m p r e s i o n PL - Plastic limit, %
C -
v C o e f f i c i e n t of Q - Bearing load, k N
consolidation, m / y 2
Quit- U l t i m a t e load, k N
d - Depth, m q - B e a r i n g pressure, k P a
D - Grain size, m m - D r a i n a g e rate, mm/day
- O v e r b u r d e n depth, m - H e a t flow, J / ( m s ) 2
s - F i e l d slope, %
- Side w i d t h , m
S - D e g r e e of s a t u r a t i o n
- Wheelslip, %
SF - Safety factor
t - Time, s
Τ - Torque, k N m
- Tangential force, k N
- Tensile s t r e n g t h , k N / m
- Temperature, ° C
T -
v Time factor
u - P o r e w a t e r pressure, k P a
U - Uniformity coefficient
- D e g r e e of consolidation
- Water force, kN/m
ν - W a t e r v e l o c i t y , m/s
- Specific volume, m / t 3
V - Volume, m 3
V -
s Solids v o l u m e , m*3
V -
v V o i d s v o l u m e , m.3
Vw - Water volume, m 3
1.1. O R I G I N S O F SOILS
S o i l s c a n be defined as m a s s e s of m i n e r a l p a r t i c l e s m i x e d with v a r y i n g
p r o p o r t i o n s of w a t e r , gases and s a l t s , and o f t e n o r g a n i c m a t t e r . S o i l s are
s o m e t i m e s f o r m e d physically and c h e m i c a l l y in place f r o m parent r o c k
m a t e r i a l , a n d c a n o c c u r f r o m this o r i g i n in l a y e r s up to 25 m deep.
A l t e r n a t i v e l y , they c a n be r e m o v e d f r o m r o c k m a t e r i a l by physical and
chemical w e a t h e r i n g and subsequently be t r a n s p o r t e d and s e d i m e n t e d to a new
location.
S a n d a n d g r a v e l , and m a n y silt m a t e r i a l s (individual p a r t i c l e s larger than
0 . 0 0 2 m m d i a m e t e r ) , are generally of the s a m e m i n e r a l o g i c a l c o m p o s i t i o n as
t h e i r parent r o c k s s u c h as q u a r t z , f e l d s p a r s , c a l c i t e , d o l o m i t e , m i c a s and
o t h e r s . T h e y are r e m o v e d f r o m m a s s i v e r o c k f o r m a t i o n s usually by physical
agents including w i n d , w a t e r and f r e e z i n g . T r a n s p o r t a t i o n of these granular
m a t e r i a l s by w i n d , w a t e r , ice or falling d o w n slopes c a n alter the shape and
size of individual p a r t i c l e s by abrasion and i m p a c t , and c a n also s o r t different
g r a i n s i z e s one f r o m another. T h a t is the r e a s o n w h y , for e x a m p l e , m a n y
shoreline beaches or desert sand dunes are found to have r e m a r k a b l y u n i f o r m
s i z e s of sand p a r t i c l e s . Y e t , s o m e g r a n u l a r m a t e r i a l s s u c h as g l a c i a l till
c o m p r i s e a wide m i x t u r e of g r a i n s i z e s , f r o m silt to boulders larger than
300 m m in d i a m e t e r , because of their mode of transport in ice.
Exceptions to the above case do o c c u r , s u c h as the c r y p t o c r y s t a l l i n e quartz
s h o w n in F i g . 1 . 1 . T h e s e silt s i z e d q u a r t z p a r t i c l e s have been r e f o r m e d f r o m
dissolved c o n s t i t u e n t s in a hot spring in A r k a n s a s , but they have f o r m e d into
the c h a r a c t e r i s t i c rounded shape of q u a r t z c r y s t a l s found in m o s t sands.
C l a y m i n e r a l s , on the other h a n d , are f o r m e d m o r e usually by c h e m i c a l
changes to produce new m i n e r a l s s u c h as k a o l i n i t e , illite, c h l o r i t e , v e r m i c u l i t e ,
m o n t m o r i l l o n i t e and others. T h e s e new m i n e r a l s tend to break easily into s m a l l
( l e s s t h e n 0.002 m m ) s i z e d p a r t i c l e s and have rather c h e m i c a l l y a c t i v e
s u r f a c e s , as well as a large s u r f a c e area per unit m a s s . F i g . 1.2 s h o w s a
kaolinite m a t e r i a l w h i c h has been f o r m e d as a residual soil on parent r o c k in
G e o r g i a , U . S . A . T h e large b l o c k s of c l a y platelets have not yet been broken
up into individual p a r t i c l e s due to the l a c k of s u f f i c i e n t physical and c h e m i c a l
a c t i o n at this point in their h i s t o r y . F i g . 1.3 s h o w s p a r t i c l e s of the s a m e
m i n e r a l in a finer s t r u c t u r e . F i g . 1.4 illustrates a m i x e d clay soil f r o m the
C h a m p l a i n S e a deposits in Q u e b e c and O n t a r i o . It c o m p r i s e s illite, chlorite
and m u s c o v i t e c l a y s i z e d p a r t i c l e s , as well as silt s i z e d r o c k g r a i n s . F i g . 1.5
2 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
is a m i c r o p h o t o g r a p h of a m o n t m o r i l l o n i t e clay f r o m the N o r t h A m e r i c a n
p r a i r i e s . T h i s is a m i n e r a l w h i c h c a n absorb considerable a m o u n t s of water
and swell to s e v e r a l t i m e s its dry v o l u m e .
Gravel Sand
AASHO Boulders Silt Clay
coarse medium fine coarse fine
Sand
ASTM Gravel Silt Clay
coarse fine
Sand
FA A Gravel Silt Clay
coarse fine
Gravel Sand
USDA Cobbles Silt Clay
e r
c o a r s e mediurrι fine Y *
coarse fine coTrse fine
F i g , 1.6. C l a s s i f i c a t i o n s y s t e m s o f soils by g r a i n s i z e s .
ORIGINS AND CLASSIFICATION OF SOILS
5
6 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
1.2. C L A S S I F I C A T I O N O F SOILS
100 90 80 70 60 50 40 30 20 10 0
Percent S a n d
P r o b l e m 1.2. A m e c h a n i c a l a n a l y s i s is p e r f o r m e d on a p p r o x i m a t e l y 1000 g of
a pulverized soil s a m p l e , w i t h the f o l l o w i n g results. W h a t is the textural class
of this s o i l ?
r e s p e c t over the y e a r s .
W h e t h e r a soil is well g r a d e d or poorly graded depends on the v a r i e t y of
d i f f e r e n t g r a i n s i z e s c o n t a i n e d in the m a t e r i a l . O n e m e a s u r e of g r a d a t i o n is
called the ' u n i f o r m i t y c o e f f i c i e n t , U , w h e r e ;
1
U = D 6 0 /D 1 Q (1.1)
I m p o r t a n t properties
N a m e of group Symbols Compacted Compacted Compressi- Worka-
permeability strength bility bility
Well graded
gravels. GW pervious excellent negligible exec.
P o o r l y graded
gravels. GP high good negligible good
Silty gravels. GM f a i r to poor good negligible good
Clayey gravels. GC impervious low v e r y low good
Well graded sands. SW pervious excellent negligible exec.
P o o r l y graded sands SP pervious good v e r y low fair
Silty sands. SM f a i r to poor good low fair
C l a y e y sands. SC impervious fair-good low good
Inorganic silts. ML f a i r to poor fair medium fair
I n o r g a n i c c l a y s of fair to
low p l a s t i c i t y . CL impervious fair medium good
Organic silts. OL f a i r to poor poor medium fair
Inorganic silts,
elastic. MH f a i r to poor fair-poor high poor
Plastic clays. CH impervious poor high poor
Organic clays. OH impervious poor high poor
PI = L L - PL (1.2)
1.3. S O I L PHASES
V o i d ratio e = V v /V
' s (1.3)
Porosity η = Vv (1.4)
D e g r e e of s a t u r a t i o n S
( W x 100% (1.5)
Weiqht of w a t e r
y W
w
=V V (1.6)
S p e c i f i c g r a v i t y of solids G g = (1.7)
Weight of solids W g = 7s s = 7w s s
V G V (1.8)
1.4. P R O B L E M S
P l a s t i c l i m i t test N o . 1 2 3 4 5
M a s s of can plus wet soil, g 30.96 29.56 30.46 29.19 28.22
M a s s of c a n plus dry s o i l , g 28.84 27.76 28.37 27.42 26.66
M a s s of drying c a n , g 19.94 20.22 19.78 20.06 20.15
L i q u i d limit test N o . 1 2 3 4 5
N u m b e r of blows 21 23 26 28 30
M a s s of c a n plus wet soil, g 85.67 96.79 93.18 83.32 92.27
M a s s of c a n plus dry s o i l , g 58.30 64.42 61.67 56.20 60.13
M a s s of drying c a n , g 20.34 20.80 19.93 21.02 19.85
A n s w e r : S = 7 9 % , J= 17.0 k N / m ,
3
γ. = 13.1 k N / m . 3
12 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
s = c + a tan</>
n (2.1)
where:
s = soil shear s t r e n g t h ( m a x i m u m shear stress)
c = c o h e s i o n , the part of s t r e n g t h independent of n o r m a l pressure
(Tn = n o r m a l stress on the sliding plane
φ = angle of internal f r i c t i o n of the m a t e r i a l .
Normal stress, σ π
I V
\ * ^
^.Developed
shear strength, s
Soil'
internal sliding plane
U n i t m a s s , t/m?
Moist 1.5 1.8 1.9 2.1 2.1
Submerged 0.9 0.95 1.05 1.2 1.2
1
1 n c r e a s e f i v e degrees for soils c o n t a i n i n g less than 5 % fines.
U n d r a i n e d shear
strength, kPa 0-25 25-50 50-100 100-200 200-400 400
S t a n d a r d penetration
r e s i s t a n c e , blows/ft 0 2 4 8 16 32
2.7).
W i t h the knowledge of the cohesional and f r i c t i o n a l s t r e n g t h p a r a m e t e r s
of a soil, E q n . 2.1 c a n be used directly to find the u l t i m a t e f o r c e on a plane
of a g i v e n area in a soil, if it is k n o w n that that plane is the failure s u r f a c e .
O n e s u c h simple example is the failure s u r f a c e beneath a d r i v e n lugged tire
o f a v e h i c l e . A t the m a x i m u m tire t r a c t i o n f o r c e , slip will o c c u r in the soil
o n a h o r i z o n t a l plane beneath the tire. If the tire lugs have penetrated the
soil surface, then the failure plane will be essentially a soil to soil i n t e r f a c e ,
between the m a t e r i a l trapped in the s p a c e s b e t w e e n lugs and the soil beneath.
Η = sA = cA + σ A tan<f> = c A + W tan φ
m r
n
r
n
η r ^ r r ^
I n t h i s c a s e , the cohesion and angle of internal f r i c t i o n of the soil are
k n o w n , t h e c o n t a c t area of the rear tires is the loaded w i d t h of e a c h
SOIL S H E A R S T R E N G T H 17
H m = 30 kPaxO.432 m x 0 . 6 6 m x 2 + 0.75x4x9.8 t a n 2 5 ° k N
s = c + σ tan δ (2.2)
α Π
c i r c l e of stresses.
P a r a l l e l to the plane:
T dx/cos0 + a d x t a n 0 c o s 0
3 - ajdxs^ = 0
Τ = ( σ - a^sin0cos0
1 = | 1~ 3 |
σ a
s i n 2 g (2Λ)
Problem 2.2. A n example is given in Fig. 2.5 wherein normal stresses o and zz
particular Mohr's circle of stresses for this point can be constructed from the
stress combinations on these two mutually perpendicular planes as shown in
Fig. 2.5. The centre of the circle can be calculated as;
V{ xx (a
* ζζ
σ ) / 2
} 2 +
r
xl = ^ { ( 6 0 - 30)/2J + 1 0 2 2
kPa = 18.03 kPa
Then the principal stresses are equal to the normal stress at the centre
of Mohr's circle plus or minus the radius as follows.
In addition, the clockwise angle of the major principal plane, 0 1 ? from the
ζ plane can be found from the geometry of Mohr's circle as:
20 2 = 90° - 2 0 l 5 0 2 = 28.15°
Note in Fig. 2.5 that the S point in Mohr's circle is 90° from the σ-j point,
and therefore the S plane in the material body is inclined at 45° relative to
the major principal (σ-j) plane in the same direction. Further examples of this
nature are provided for the purpose of practice at the end of this Chapter.
SOIL S H E A R S T R E N G T H 21
Failure strength
criterion
s = c + (7 tan0
n
20 f = (180° - ^ F A O ) = (90° + φ)
0 f = 4 5 ° + φ/2 (2.5)
OO 1
= c/tan</> = c c o t 0 (2.6)
= sin<p = (2.7)
AO' (σ + σ )/2 + c οοίφ
λ 3
-40L
{^} [ ° = (9 +30)721
y{< 90
" 30)/2} + 10 ] k P a =
2 2
|J1.6j k R a
φ = sin- 1
I 91.6 - 28 4) ο
\ 91.6 + 28.4/
0 = {90° + 31.8° - tan- (10/30)}/2 = 51.7°
1
1
^ i ^ - 1
{*$}-*] (2J0)
σ
χχ =
^^Y 1
' sin0cos2e} (211)
Τ χ ζ = { σ ΐ
2 ° 3
}{sin</>sin2e} (212)
2.4. T H E M E T H O D O F S T R E S S CHARACTERISTICS
p r i n c i p a l s t r e s s , θ\ m e a s u r e d c l o c k w i s e f r o m the h o r i z o n t a l or χ a x i s . W h e n
s o i l f a i l u r e h a s been a s s u m e d , these t w o v a r i a b l e s are s u f f i c i e n t to describe
c o m p l e t e l y the stress state at a point, and the principal s t r e s s e s or stresses
on the horizontal and v e r t i c a l planes in the m a t e r i a l c a n be obtained directly
f r o m these variables using the equations below. B y e m p l o y i n g these t w o
v a r i a b l e s , the number o f p a r a m e t e r s needed to describe the state o f s t r e s s e s
a t a point has been reduced f r o m three to t w o , by utilising the a s s u m p t i o n
of the soil failure condition.
σ ζ ζ = σ ( 1 - sin</>cos20) - ψ (2.15)
r
x z = '^ζχ =^ n0sin20 s i
(2.16)
In the χ d i r e c t i o n :
σ dz - q dz - -RV**dxdz + T d x + tf^zXdxdz - r dx = 0
xx xx Ox ζχ σζ ζχ
S i m i l a r l y , in the ζ d i r e c t i o n :
C218)
28 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
In their present form, Eqn. 2.19 and 2.20 are virtually impossible to
integrate over finite χ and ζ displacements because they comprise partial
differentials of the stress variables a and 0. However, Sokolovsky (1956)
demonstrated that the two expressions can be rearranged by multiplying them
by sin(0±/>O and -cos(0±/i), respectively, with the following result.
A l o n g £ slip lines:
£ = tan(0 + μ) (222)
A l o n g η slip lines:
S = tan(0 - μ) (224)
σ η = σ'
η + u (2.26)
or
σ ' η = σ η - u (2.27)
where:
σ η
!
= e f f e c t i v e i n t e r g r a n u l a r s t r e s s (average per unit area),
s
= ' + a 'tan0
c
n
!
(2.28)
Fig. 2.11. M o h r ' s circle of s t r e s s e s for a soil in the failure condition using
total and intergranular (effective) s t r e s s e s .
σ η = 7dry x 2 m +
%at x 3 m =
^ 1 3 x 2 + 1 8 x 3
^ k P a = 8 0 , 0 k P a
u = 9.8 k N / m 3
χ 3 m = 29.4 k P a
F i g . 2.15. R e s u l t s of shear f o r c e v e r s u s h o r i z o n t a l d e f o r m a t i o n m e a s u r e m e n t s
for a sandy soil in a d i r e c t shear box.
f 4
*
7
i
kPa
80 40!
CA φ=0
60
Test points 3$·
Τ, Ν c
3
4
40 20]
c-
20! 10
0 0
0 20 40 60 80 100 120 140 10 20 30 40 50 σ =
Ν, Ν
a -u,kPa
n
Cell
Top cap
Porous stone
Soil sample
Membrane
Cell fluid
Flexible tube
Porous stone
Elastic ring
Stress
model
Ρ 1
= Ρ - σ Α 3 ρ (2.29)
A = A /(i - e0 z z ) (2.30)
where
A = the soil cylinder c r o s s s e c t i o n area at any t i m e ,
A ο = the initial area,
€ = the v e r t i c a l s t r a i n of the soil s a m p l e = Δ ζ / ζ (compression).
z z 0
A = V/z = ( V Q + A V ) / { z ( l - e )\
Q zz (2.31)
where
V = the soil s a m p l e v o l u m e at f a i l u r e ,
ζ = the s a m p l e height at f a i l u r e ,
V 0 = the initial soil cylinder v o l u m e ,
A V = the change in soil v o l u m e f r o m the initial value (expansion),
z = the
Q initial soil cylinder height.
SOIL S H E A R S T R E N G T H 39
1 13 78 3
2 70 157 15
3 120 238 28
4 155 302 35
F i g . 2.20. P r o b l e m 2.5 c o n s t r u c t i o n of M o h r s c i r c l e s of e f f e c t i v e
f
failure
s t r e s s e s to obtain e f f e c t i v e c o h e s i o n and f r i c t i o n angle.
is d r a w n in F i g . 2.20. T h e r e s u l t i n g e f f e c t i v e s t r e n g t h p a r a m e t e r s are c = 1
20 k P a a n d 0 ' = 1 6 ° .
T h e e s t i m a t i o n of s t r e n g t h p a r a m e t e r s f r o m s e t s of M o h r ' s c i r c l e s at soil
f a i l u r e , s u c h as in F i g . 2.20, c a n be a c c o m p l i s h e d by g r a p h i c a l m e a n s on a
l a r g e scale g r a p h , and the values of c o h e s i o n and f r i c t i o n angle m e a s u r e d .
A l t e r n a t i v e l y , E q n . 2.32 below, f r o m E q n . 2.9, c a n be adapted to an a n a l y t i c a l
e s t i m a t i o n of c and φ , or c and 0 ' f o r e f f e c t i v e s t r e s s e s , using the total or
f
e f f e c t i v e principal s t r e s s e s at soil f a i l u r e .
σ
3 = a
l K
a "
2 c
V^ ( 2
· 3 2 )
(233)
Test No. 1 2 3
c = 19.9 k P a , φ = 2 4 . 9 ° , c 1
= 21.4 k P a , φ' = 2 7 . 5 °
E q n . 2.29 has been used to find the net f o r c e , P , applied to the top cap
f
1100 n a s D
established during a f e w trials w i t h i n c r e a s i n g loads, P, the
e e n
0.1
0.2
0.3
Settlement
A h , mm 0.4
l 0.5
0.6
0.7
e f = wG s (2.34)
e = e + * - f h
(235)
f h (l + e )
f f
w %
p
dry = TT^ ( 2
· 3 6 )
Pdry =
vi/S + 1 / G q
( 2 3 7 )
2.7. F I E L D S T R E N G T H TESTS
TP—Ψ
1- S h e a r ring 2. Shear plate
Ν J,
Spring-*
J-Scale
Piston
4. S h e a r vane 5. Penetrometer
6. Pocket
penetrometer
soil s t r e n g t h p a r a m e t e r s by the f o l l o w i n g r e l a t i o n .
T m = ( c A + Ntan</>)r (2.38)
F m = c A + Ntan</> (2.39)
in w h i c h A is the h o r i z o n t a l plate a r e a .
I n p r a c t i c e , the shear ring is m o r e c o n v e n i e n t to operate in the field
because it is both m a n u a l l y or m e c h a n i c a l l y easier to provide a f o r c i n g torque
to the d e v i c e , than to find the a n c h o r a g e and s t r e n g t h needed to apply large
h o r i z o n t a l f o r c e s to the plate. A l t h o u g h neither the shear r i n g nor plate are
readily available for purchase, m a n y m o d e l s have been built over the y e a r s by
r e s e a r c h e r s and soil testing t e a m s . T h e y have ranged f r o m s i m p l e manual
d e v i c e s , i n w h i c h the v e r t i c a l f o r c e is p r o v i d e d by dead w e i g h t s , and the
torque or h o r i z o n t a l f o r c e applied by handles, to a u t o m a t e d and e l e c t r o n i c a l l y
m o n i t o r e d h y d r a u l i c or e l e c t r i c p o w e r e d m a c h i n e s m o u n t e d on v e h i c l e s . T h e
t y p e of d e v i c e required depends on the r e s o u r c e s a v a i l a b l e , and the n u m b e r
of m e a s u r e m e n t s to be m a d e . L a r g e r e s e a r c h institutes w i t h s u f f i c i e n t funding
o f t e n opt f o r an a u t o m a t e d v e h i c l e - m o u n t e d m o d e l , w h i l e typically
p o s t g r a d u a t e r e s e a r c h e r s w i t h l i m i t e d r e s o u r c e s are f o r c e d to build and
m a n u a l l y operate their o w n a p p a r a t u s . In the latter c a s e , the applied torque
i s g e n e r a l l y m e a s u r e d by s o m e s o r t of s p r i n g g a u g e d e v i c e . T h i s does not
imply that the a c c u r a c y of m e a s u r e m e n t s m u s t be c o m p r o m i s e d in either c a s e ,
p r o v i d e d that c a r e is taken to insert and operate the shear testing d e v i c e s
properly.
A = (250 2
- 200 )/4 m m
2
2
= 17,700 m m 2
= 0.0177 m 2
T /r
m = c A + Ntan</>
Τ = /VrdA = [ΐττν ύτ
2
= 1Τ7Γγ 3
(2.40)
F i g . 2.28. R e s u l t s of a s e r i e s of s h e a r g r a p h t e s t s on a clay l o a m s o i l .
54 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
F i e l d d e v i c e n u m b e r 5 in F i g , 2.27 is a p e n e t r o m e t e r . U n l i k e the
a f o r e m e n t i o n e d tests, the i n f o r m a t i o n r e c e i v e d f r o m the p e n e t r o m e t e r does
not relate to the C o u l o m b soil s t r e n g t h p a r a m e t e r s , c and φ . I n s t e a d , a single
quantity is m e a s u r e d , n a m e l y the f o r c e required to m o v e the p e n e t r o m e t e r at
a s e l e c t e d c o n s t a n t speed and a p a r t i c u l a r depth in the soil. T h i s f o r c e is
o f t e n n o r m a l i z e d by dividing by the the base area of the p e n e t r o m e t e r tip,
t o a r r i v e at the a v e r a g e applied v e r t i c a l pressure, k n o w n as the cone index.
While either the penetration f o r c e or the cone index by itself does not allow
t h e e s t i m a t i o n of soil cohesion or f r i c t i o n angle, it has been used o f t e n to
g i v e a n indication of the r e l a t i v e soil s t r e n g t h , density, root p e n e t r a t i o n ,
t r a f f i c a b i l i t y by vehicles or other design a p p l i c a t i o n . T h e s e uses of the
p e n e t r o m e t e r require extensive testing and c a l i b r a t i o n of the cone index
m e a s u r e m e n t s w i t h respect to the application d e s i r e d , but m a n y e x a m p l e s of
this procedure are available in the soil m e c h a n i c s l i t e r a t u r e .
The pocket p e n e t r o m e t e r is designed for an application d i f f e r e n t f r o m the
g e n e r a l p e n e t r o m e t e r test. A c i r c u l a r flat tipped rod is pressed into a soil
s u r f a c e up to an indicated depth. B y the theory of bearing c a p a c i t y of
f o o t i n g s on cohesive soils ( C h a p t e r 3) the required f o r c e is c o n v e r t e d to the
u n d r a i n e d shear s t r e n g t h , C , and this is indicated on the c a l i b r a t e d spring
s c a l e o n the i n s t r u m e n t . T h i s device is useful for a rapid and a p p r o x i m a t e
d e t e r m i n a t i o n of the undrained shear s t r e n g t h and bearing c a p a c i t y of shallow
c o h e s i v e undrained soils.
C o m p a r i n g the applicability of all the tests illustrated in F i g . 2.27, only
t h e shear vane and pocket p e n e t r o m e t e r are a l w a y s r e s t r i c t e d to use in
s a t u r a t e d c o h e s i v e soils. B o t h of these d e v i c e s g i v e an indication of the
u n d r a i n e d shear s t r e n g t h , C , w i t h the a s s u m p t i o n that total s t r e s s angle of
f r i c t i o n is z e r o . T h e shear r i n g , shear plate and s h e a r g r a p h c a n give drained
strength p a r a m e t e r s , provided that the soil in question has adequate hydraulic
c o n d u c t i v i t y to release developed pore w a t e r pressures over the t i m e period
of each test. E v e n for the few c e n t i m e t r e s of d i s t a n c e b e t w e e n the soil under
o n e o f t h e s e shearing d e v i c e s , it c a n take m o r e than one hour for the pore
water to m o v e and reduce the pore pressure, if the soil h y d r a u l i c c o n d u c t i v i t y
is l e s s than 10 ~ c m / s e c . B e c a u s e this is not a p r a c t i c a l t i m e delay for the
4
2.8. PROBLEMS
A n s w e r : c = 8.2 k P a , φ = 2 5 . 8 ° .
Test N o . σ , kPa
3 Ρ, Ν Δζ, m m Δ\/, c m 3
A n s w e r : c = 10.2 k P a , φ= 30.2°.
58 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
Answer: t 1 Q Q = 4 hours.
w, % 5 10 15 20 25 30
T o t a l m a s s , kg 4.84 5.00 5.22 5.41 5.42 5.37
3.1. B E A R I N G CAPACITY
A. Punching B. Rotation
F i g . 3 . 1 s h o w s t w o c o m m o n m o d e s of f o u n d a t i o n failure under a v e r t i c a l
f o r c e , Q. If the foundation is c o n s t r a i n e d to descend in a level attitude, for
e x a m p l e by a r i g i d c o n n e c t i o n to a v e r t i c a l c o l u m n , then the foundation
punches into the soil as s h o w n in c a s e A . A f r e e standing load usually c a u s e s
a foundation to r o t a t e as in case B, s i n c e the soil tends to m o v e to one side
o f t h e f o u n d a t i o n before the other. I n f a c t the a n a l y s i s for the f o r c e w h i c h
c a u s e s incipient failure is p r a c t i c a l l y identical f o r both c a s e s .
smooth c o n t a c t area w i t h the soil, and that the depth of the a c t i v e zone (soil
m o v i n g downwards) is the s a m e as that of the passive zone (soil m o v i n g up).
T o t h e s i d e s of the foundation base and at the s a m e e l e v a t i o n , there is a
vertical pressure a c t i n g in the soil a r i s i n g f r o m a depth D of o v e r b u r d e n . The
m a g n i t u d e of this pressure is:
q =γΌ (3.1)
where y %
is the e f f e c t i v e w e i g h t density of the soil above the f o u n d a t i o n .
<7i = % y
+ z
where Κ = } '
a 1 + sincp
I f it is a s s u m e d that there is p r a c t i c a l l y no s t r e s s c h a n g e b e t w e e n points
A and C , then σ at point A equals
3 at point C , since the major principal
s t r e s s c h a n g e s d i r e c t i o n b e t w e e n the t w o points. A t point C , the v e r t i c a l
s t r e s s is;
# 3 = yz + q = yz + y'D
q = γζ(Κ 2
- 1) + q K 2
+ 2cv1<~(l + Κ )
Ό ' ρ ρ ρ ρ
where Κ = 1/Κ = } l ^ + s n
02)
ρ a 1 - sin0
q n = (K 2
- 1)γΒ/2 + 2c\/K"(l + Κ ) + qK 2
ο ρ ' ρ ρ ρ
= (7Β/2)Νγ+ cN c + qN q (3.3)
Y^ 7T/2/'
= 2 t a n 0 d 0 (3.5)
σ Α = a e*" c
t a n
0 (3.6)
σ
Α " 1 ι sin0 " c a
- e
χ 1 - sin</> / e
r = r e" ^ t a n
0
c
64 A G R I C U L T U R A L ENGINEERING SOIL MECHANICS
_ d _ = _S_J -ftan0l
e
cos/i sin/i\ j
d» = d tan/A e f tan<
A
Vyfp^e^^-cot^ »*φ[{β^^-ΐ] + q
= (γΒ/2)Νγ + c N c + q N q (3.8)
Qo
(a)
Rigorous solution
for c=0; γ , 0 *0
P r o b l e m 3 . 1 . C o n s i d e r a long n a r r o w r i g i d s t r i p f o o t i n g to be installed at a
d e p t h o f 1 m in a clay l o a m soil h a v i n g a m o i s t density of 17.6 k N / m 3 .
c=15 k P a and φ = 3 0 ° . If the strip f o o t i n g is to be designed to c a r r y a load of
500 k N per m e t r e length s a f e l y , w h a t should be the w i d t h of the f o o t i n g ?
= 1500 kN/m/B
Thus 1500 k N / m = ( 1 4 1 B + 7 7 4 B ) k N / m
2
—3
\/ /
V /
'111-
ς = | Β Ν [ ΐ - 0 . 3 5] [l-1.5 Hj2
ο Ύ γ + y D N q [i + 0 . § ][l
2 + 0 .l g][l-1.5 &12)
F o r c o h e s i v e soils ( φ = 0):
β°
F i g . 3.9. B e a r i n g c a p a c i t y f o r m u l a Ν f a c t o r s under inclined loads as
e s t i m a t e d by S o k o l o v s k i (1956) and B r i n c h H a n s e n (1957).
7'ii
1 I l q
7
Total density y
% =7 b u o y (B/2)N<y + c N c + qN q
(3.15)
where Ύ, = Ύ - Ύ (3.16)
'buoy ι /w
and y = the w e i g h t density of w a t e r = 9.8 k N / m 3
q- = 7 » D l + (y . 7 W )D 2
(3.17)
7 b u o y = (18.6-9.8) k N / m 3
= 8.8 k N / m 3
q f
= [ 1 3 . 7 x 0 . 6 + 8.8x0.7 ] k P a = 14.4 k P a
q 0 = [8.8x0.5x22.5(1-0.3x0.5) + 14.4x23(1-0.2x0.5) ] k P a = 3 8 2 k P a
Q u l t = 382 k P a x l m x 2 m = 764 k N
% = c N^ x (3.18)
f* Β
Qo
ν ν
I
1
7 = typi + γ ο 2 2
+
73 3 D ) / ( D
1 + D
2 + D
3> ( 3
· 1 9 )
c
= [ l( l + L ) + c (L
c L
5 2 2 + L ) + c L3]/[L
4 3 1 + L 2 + L3 L + 4 + L ] 5 (3.20)
The same principle could be used for two or four l a y e r s if so desired, and
a t l e a s t a crude e s t i m a t e of u l t i m a t e pressure c a n be obtained using these
a v e r a g e p a r a m e t e r s in the bearing c a p a c i t y e q u a t i o n .
Κ Β'- 5
•ζ Β ^
Effective width B f
= Β - 2e (3.22)
E f f e c t i v e bearing area A !
= BM_ (3.23)
i.Oi
-V \ -
Dhesive s oil
0.8
v^ C <
\ -
^Granu ar soil
0.6
Reduction
factor, R
- -
e
= '/A
A
0.4
- -
0.2
-
0
0.5
Eccentricity ratio, e/B
F i g . 3.15. The bearing c a p a c i t y r e d u c t i o n f a c t o r for e c c e n t r i c loading in the
c a s e s of c o h e s i v e and g r a n u l a r soils ( M e y e r h o f , 1953).
76 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
E f f e c t i v e bearing area A 1
= B'L 1
= (B - 2 e ) ( L - 2e )
b { (3.24)
B f
= Β - 2e = (1.0 - 0.4)m = 0.6 m , and L ' = (2.0 - 0.4)m = 1.6 m
q 0 = 0 . 5 γ Β Ή γ ( 1 - 0.3BVL ) + c N ( l + 1
c 0.2B'/U)
Q
design = q 0
B
' L I
/ S F
= 1 8 8 k P a x
°* 6 m x
l - 6 m / 3 = 60 kN
3.2. F O U N D A T I O N SETTLEMENT
ζ = qBC 2 (3.26)
where ζ = foundation s e t t l e m e n t ,
q = a v e r a g e foundation bearing pressure
^ 2 = a soil s e t t l e m e n t c o n s t a n t (in units of pressure" ) 1
ζ = qC^ (3.27)
z_ = 4 028)
z " (1 + B / B )
0 0
2
SHALLOW FOUNDATIONS 79
q = kz n
(3.29)
ζ = total plate s e t t l e m e n t
η = a c o n s t a n t exponent for a soil
q = (k /B + k ) z
c 0
n
= kz n
(3.30)
k :
c kPa/m" " 1 1
or p s i / i n " " ,
1
k : kPa/m
0
n
or p s i / i n n
f r o m F i g . 3.18.
k c = B B (k
1 2 1 -k )/(B
2 2 - Β ) = -0.323 k P a / m "
χ
n 1
(3.32)
k
0 = 0<2 2 - k j B j V ^
B
- Β ) = 71.1 k P a / m
χ
n
(3.33)
η = (n 2 + n ) / 2 = 0.566
2 (3.34)
i n t e r c e p t , k 0 , a r e m o r e a c c u r a t e l y c o m p u t e d using a c a l c u l a t o r linear
regression.
F a n (1985) d e m o n s t r a t e d that the c o n f i d e n c e in s i n k a g e c o n t e n t s increases
m a r k e d l y if sets of three or four plates h a v i n g d i f f e r e n t s i z e s are used for
s e t t l e m e n t t e s t s . T a b l e 3.1 g i v e s the a v e r a g e c o e f f i c i e n t s o f v a r i a t i o n in the
constants for four series of s i n k a g e t e s t s , using i n c r e a s i n g n u m b e r s of plates.
I t is evident that in the type of soil t e s t e d , i m p r o v e m e n t s in the c o n f i d e n c e
o f the k values c a n be i m p r o v e d by m o r e than five t i m e s if four plates of
d i f f e r e n t w i d t h s are used instead of t w o .
82 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
3601 I 1 1 1 1 1
0 10 20 30 40 50 60
1/B, m" 1
F i g . 3.19. G r a p h i c a l d e t e r m i n a t i o n of k c and k^ p a r a m e t e r s .
T A B L E 3 . 1 . A v e r a g e c o e f f i c i e n t s of v a r i a t i o n of s i n k a g e c o n s t a n t s w h e n
d i f f e r e n t n u m b e r s of plates are used ( M c K y e s and F a n , 1985).
N o . plates A v e r a g e c o e f f i c i e n t of v a r i a t i o n , %
used η
k
c
d^, c m p, k P a z, c m d ,
2 cm p, k P a z, c m d^, c m p, k P a z, c m
d, c m 5 7.5 10
k, k P a / m n
2683 1751 1602
η 0.719 0.694 0.711
η = (n 1 + n 2 + n ) / 3 = 0.71
3
κ = (l)k ^ c+
( C o r r e l a t i o n c o e f f i c i e n t R = 0.98)
3.3. P R E S S U R E D I S T R I B U T I O N U N D E R FOUNDATIONS
V/Q I
1 Shear force
,5
X/B
Raft footing
Bending moment
Q Pressure distribution ^
.5
Κ Β
X — * — * — * — * — τ — *
F i g . 3 . 2 1 . A wide r a f t f o o t i n g .
86 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
bending m o m e n t d i a g r a m s are s i m i l a r t o s u c h a c a s e . I f t h e f o o t i n g is
c o n s i d e r e d to be r i g i d , then considerable bending m o m e n t s c a n be generated
along the centreline due to the r e l a t i v e l y large s i z e of w i d t h , B .
Differential
settlement
h Β >-
q = q /SF = Q/A = ( Q
0 e U a g e + Q s i i Q + Q f o o t i n g )/(jB ) 2
(3.35)
a s s u m e d to be a c o m p l e t e disc.
Q f o o t i n g /A = 23.6 k N / m 3
χ 0.6 m = 14.2 k P a
T h e u l t i m a t e bearing c a p a c i t y is c a l c u l a t e d by E q n . 3.14 f o r a c i r c u l a r
f o o t i n g , using the N f a c t o r f r o m A p p e n d i x 1 .
Q
30 -τ— I ι 1 •— 1 1
\ D
Ti It 0.3 percent
Ti t > 0.3 percent "
ι Δ
\• •
20 δ\ -
Average ID
foundation
\
• A
settlement
V
•
cm
\
10
- • \ \ Δ
γ •
\ •
• V 0
•
, •.• \ ttp
V n π
Ρ • •
π
υ0 1 2 3 4 5 6
Safety factor
2.5 V e r y good
3.5. P I L E FOUNDATIONS
(a) End or point-bearing pile (b) Friction pile (c) Tension pile
Q
ult = 7ΓΓ (0.67ΓΝ 2
Γ + 1.3cN c + yDN )
q (3.36)
Q u l t = 2 7TrrD (3.37)
7 = s = c
a +
tf ^ tan (3.38)
n
94 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
<Tn = σ
ζζ ο
Κ
= 7 z K
o (3.39)
F o r s t r a i g h t s m o o t h pile sides, the c o e f f i c i e n t of earthpressure at rest
suggested by T e r z a g h i (1943), K = 0 . 5 , c a n be used. H i g h e r c o e f f i c i e n t s m a y be
o
allowed for different shapes of pile. Ireland (1957) has s u g g e s t e d that K 0 can
be a s h i g h as 1.75 for step-tapered c o n c r e t e piles, w h i c h b e c o m e s m a l l e r in
diameter w i t h depth. T h e s a m e E q n . 3.37 c a n be used to find the tension load
capacity, T , for u n i f o r m d i a m e t e r tension piles as well, since the c o n f i g u r a t i o n
of shear stress along the pile length is the s a m e , only r e v e r s e d in d i r e c t i o n .
Q g = q u l ( . B L + D ( 2 B + 2L)s (3.40)
F i g . 3.28. Negative skin friction on deep end bearing piles due to soil
settlement.
Q fc = 35,0007Tr 2
kPa = Q + 2 7Tr(c D a n π-γϋ^Κ tan6/2)
3.6. PROBLEMS
Answer: S.F. = 3 . 1 .
A n s w e r : Q = 339 k N .
A n s w e r : Β = L = 2.1 m.
98 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
and η = 0.5. A 1000 Ν load is placed on t w o plates, (a) 3x4 c m and (b)
4x4 c m , on the s u r f a c e of the f i e l d . W h a t are the e x p e c t e d s e t t l e m e n t s
of plates (a) and (b) under this l o a d ? A r e the s e t t l e m e n t s of the plates
proportional to the pressures a c t i n g on t h e m ?
A n s w e r : b = 312 m m .
A n s w e r : Β = 8.8 m.
W A T E R FLOW IN SOILS 99
4.1. C O N S O L I D A T I O N A N D S E T T L E M E N T RATES
Pore pressure u
F i g . 4 . 1 . S c h e m a t i c d i a g r a m of soil c o n s o l i d a t i o n .
V = ki = - k f r ' = (4.1)
dz 7 dz
w
dR = du (4.2)
AT " AT
where m v is the c o e f f i c i e n t of v o l u m e c o m p r e s s i b i l i t y , w h i c h is a s s u m e d
to be p r a c t i c a l l y a constant over a l i m i t e d r a n g e of soil v o l u m e c h a n g e , and
WATER FLOW IN SOILS 101
( ω
>
du = k dfu ( ^ 5 )
dt 7 m 9z
w v
2
Γ kn 7T t"|
2 2
U = i L = £ ^(l-cosn7r)sin[^]el Tw v 4 m H 2
J (40
Τ = = (4.7)
v
7w v m H 2
Time factor, T v
- C - _ E2 - GO (4.8)
c "Alogp log(p /p )
2 0
AV
Ah _ Ae
=
(4*9)
Vhi " 1 + e„
o
0 0 ο
To calculate the s e t t l e m e n t of a soil w h i c h has been n o r m a l l y consolidated
(that is has not s u f f e r e d a larger pressure at s o m e t i m e in the past), the
change in pressure, A p , c a n be used in the f o l l o w i n g e x p r e s s i o n .
(4J0)
104 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
Void
ratio
log ρ
F i g . 4.4. Consolidation c u r v e for a preconsolidated c o m p r e s s i b l e soil, and the
m e t h o d of e s t i m a t i o n of the preconsolidation pressure, p , of Q
C a s a g r a n d e (1936).
exceeded.
A commonly used method for e s t i m a t i n g the preconsolidation pressure f r o m
t h e results of a consolidation test w a s proposed by C a s a g r a n d e (1936) and is
shown in Fig. 4.4. The point Β where the m i n i m u m radius of c u r v a t u r e o c c u r s
is located, and a tangent B D and horizontal line B E d r a w n t h r o u g h this point.
P o i n t F above the e s t i m a t e d preconsolidation pressure, p , is found at the 0
thin, wide soil layer by inserting this p value in E q n . 4.10. Δ ρ in the equation
0
WATER FLOW IN SOILS 105
Δ ρ = P 2 - p Q (4.11)
AV Ah / x
<W>
m / l l o
m
v = N ^ = h ^
h
• ^ H ^ ] • #?}°^] - °·' 20 88 m
C v = 1.02xl0" m kPa/(9.69xl0"
3 2 4
k P a day) = 1.05 m / d a y 2
T v = C t/H
y
2
= 1.05t ( m / d a y ) / 1 0 0 m
2 2
= 0.20
7 / / / / / / / / / / / /
/ / / / / / /
(413)
* P = P I [ B T Z ]
Δ ρ =
PI[(B + zjo. + z)] («A)
108 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
1 4 2 12 42 39 127
2 4 6 16 31 118 41
3 4 10 20 25 196 21
4 8 16 26 19 314 20
5 10 25 35 14 490 12
Total 30 221
4.2. W A T E R F L O W I N S A T U R A T E D SOIL
q = ^ = vb = k A h k (4J5)
dt a
total v o l u m e of water flow per unit t i m e for a unit soil depth is;
% = k(h.-h ) n Πί (416)
at 1 ο Πς)
WATER FLOW IN SOILS 111
P r o b l e m 4 . 3 . In F i g . 4.9 an e x a m p l e is g i v e n of a s m a l l u n i f o r m e a r t h d a m on
a r e l a t i v e l y i m p e r m e a b l e base. (If the subbase is m o r e permeable than the
d a m s o i l m a t e r i a l , then m o r e w a t e r w o u l d f l o w below the d a m than t h r o u g h
it.) The technique w h i c h has been used in F i g . 4.9 is to draw flow lines, s u c h
a s A C and B D , n o r m a l to equipotential lines s u c h as A B . T h e equipotential
lines are d r a w n in also between A B and C D , s u c h that a p p r o x i m a t e squares
are f o r m e d by the i n t e r s e c t i n g lines. N e a r the u p s t r e a m d a m slope, A B , it is
difficult to m a k e e x a c t squares because the f l o w lines c u r v e c o n s i d e r a b l y , but
a n e f f o r t is m a d e to m a k e the distance between equal head lines about the
s a m e as the w i d t h of the f l o w p a t h .
The u p p e r m o s t line, A C , is called the p h r e a t i c line or s u r f a c e , as it is at
zero gauge pressure. T h e potential of this line at any point is therefore only
i t s height above the r e f e r e n c e e l e v a t i o n , line B D . In d r a w i n g equipotential
lines intersecting line A C , the d i f f e r e n c e in height of A C should be the s a m e
b e t w e e n p o t e n t i a l drops. T h u s the phreatic line c u r v e s d o w n w a r d s w i t h
i n c r e a s i n g n e g a t i v e slope as the equipotential lines b e c o m e closer together,
and the flow paths b e c o m e n a r r o w e r near C D .
I n the e x a m p l e s h o w n , the n u m b e r s of flow paths and head drops in the
s t r u c t u r e are n f = 2 and n ^ l l . T h e s e n u m b e r s are then substituted into
E q n . 4.16, t o g e t h e r w i t h the soil p e r m e a b i l i t y , in order to c a l c u l a t e the v o l u m e
of w a t e r l e a k a g e t h r o u g h the d a m per unit t i m e and per unit l e n g t h of the
dam. F o r i n s t a n c e , if in the e x a m p l e η-|=10 m , h = l m and k = 3 m/year, then
0
q = ^ = 3 ( 1 0 m - l m ) - | = 4.9 m / y e a r / m
m 3
length
dt y 11 -—
112 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
u = h - z Q = 0.80^ - h ) - zQ Q (4.17)
k 12.8m -H
Fig. 4.11. Water flow through soil under a c o n c r e t e d a m , and the water
pressure distribution on the d a m .
U
l = YYVJ h
l = 1 2 2
' 5 k N
/ m ( 4 e l 8 )
WATER FLOW IN SOILS 115
Ν» = W - U 2 = 300 k N / m (4.19)
T
max = c
a ( 1 2
- 8 m ) + N
' t a n
& = 178 k N / m (4.21)
D
F o r D > L/4; d = _8D l n ["Dl + ± (4.23)
7TL n
|7TrJ
For D<I_M; d =
" ^ f f f (4.24)
|_7rrJ
Fig. 4.12. The transformed flow net to parallel subsurface drain tiles (Hooghoudt, 1940).
WATER FLOW IN SOILS 117
d =7TL/[81n(L/7Tr)]
q = 8k 7TLh/[8L ln(L/7Tr)]
2
2
+ 4k h /L
x
2 2
= 0.048(m /day)/[L
2
ln(L/.05 m ) ] + 0.676(m /day)/L
2
4.3. P R O B L E M S .
A n s w e r : k = 0.015 m/y = 4 . 9 x l 0 ~ 8
cm/s.
To. 5 m
3m Silt
s / / / / / / / / / / / / / / / / 7 T / / / / / // / //////// ι / ξ
k 3m y 1m*« 3m Rock
F i g . 4.14. C o n c r e t e d a m in P r o b l e m 4.9.
|*-6m -»|
F i g . 4.15. E a r t h d a m in P r o b l e m 4.10.
A n s w e r : F l o w = 7.7 m^/m/y.
4.11. A n agricultural field has a deposit of fine sand 3.3 m deep situated on
r e l a t i v e l y i m p e r m e a b l e r o c k . P e r f o r a t e d tile drains of 10 c m d i a m e t e r
a r e placed in the field at a s p a c i n g of 15 m and an a v e r a g e depth
1.3 m . T h e s a t u r a t e d h y d r a u l i c c o n d u c t i v i t y of the soil has been
measured to be a fairly c o n s i s t e n t 0.05 m/day at all depths down to the
rock. If it is desired to m a i n t a i n the w a t e r table no higher than 30 c m
f r o m t h e s u r f a c e , w h a t p r e c i p i t a t i o n rate c a n this drainage s y s t e m
handle?
A n s w e r : q = 3.2 mm/day.
5.1. S L O P E S T A B I L I T Y DETERMINATION
0 10 20 30 40 50 60 70 80 90
Slope angle β°
AXcg
(5.2)
m
- H(S°.F.)7 " LrH
Η =lOm
7 — 7 — 7 / / / / / / — / / / / /
F r o m E q n . 5.2, m = ο / ( γ Η ) = 0.15
F i g . 5.2 shows that the m a x i m u m slope angle for φ=0 and n = l is 3 8 ° .
d
of w a t e r . T h u s the w a t e r f o r c e U is:
u
= Tw w h L ( 5
' 3 )
3 " ?w (5-4)
U = ( h z )
Ν = N + U = Wcos0
!
(5.5)
Force Τ parallel to the slice base is the actual f o r c e required to m a i n t a i n
e q u i l i b r i u m . I t s m a x i m u m possible v a l u e , T a x * found f r o m C o u l o m b ' s
m
c a n D e
T
max = N l t a n <
£' + c M
- = (Wcos0- U ) t a n φ 1
+ c'L (5.6)
5.2 S T A B I L I T Y O F W A T E R C O U R S E BANKS
The lowering of the bank slope angle will increase the overall stability of
t h e b a n k , but by itself cannot guarantee the p r e v e n t i o n of the local s l a k i n g
and subsequent failure s h o w n in F i g . 5.9. O n e m e t h o d w h i c h m a y be c o s t l y ,
but p r o m i s e s to be e f f e c t i v e in preventing bank failure due to seepage
pressures, is the installation of an interceptor s u b s u r f a c e p e r f o r a t e d drain tile
parallel to the w a t e r c o u r s e , and at a depth equal to or g r e a t e r than the ditch
b o t t o m . A s s h o w n in F i g . 5.10(a), the drain e f f e c t i v e l y l o w e r s the w a t e r table
near the w a t e r c o u r s e bank, and reduces the h y d r a u l i c pressures a c t i n g within
the slope and at its f a c e . T h e drain tile c a n have a s m a l l e r longitudinal slope
than the watercourse bed, and thus have an outlet in the ditch itself at s o m e
point d o w n s t r e a m , as s h o w n in F i g . 5.10(b).
5.3. E R O S I O N O F BANKS
ν = 0.152D 4 / 9
(G S - D 1 / 2
(5.8)
Soil m a t e r i a l M a x i m u m w a t e r v e l o c i t y , m/s
Clear water W a t e r w i t h suspension
5.4. S T A B I L I T Y O F S M A L L EARTHDAMS
5.5. PROBLEMS
Answer: β = 23°.
5.5. A 9 m high slope has a slope angle of 4 0 ° and the soil is a clay l o a m
h a v i n g a density of 17.6 k N / m , cohesion 10 k P a and angle of f r i c t i o n
3
A n s w e r : S.F. = 1.26.
A n s w e r : S.F. = 1.75.
SLOPE STABILITY 135
A n s w e r : S.F. = 1.9.
LATERAL EARTH PRESSURES 137
6.1. A C T I V E A N D P A S S I V E W A L L PRESSURES
Stresses at point Β
F i g . 6 . 1 , 6.2 and 6.3. C o u l o m b noted that to find the soil failure plane, it is
n e c e s s a r y to locate that path in the soil w h i c h will result in the l a r g e s t
p r e s s u r e on the wall during a c t i v e soil f a i l u r e . I n other w o r d s , the w e a k e s t
p l a n e in the soil, taking a c c o u n t of internal f r i c t i o n , cohesion s t r e n g t h and
soil w e i g h t , is the likely rupture s u r f a c e .
However, C o u l o m b stated also that the appropriate t h e o r e t i c a l s t r a i g h t line
soil rupture plane, s u c h as B a in F i g . 6.4(a), will yield a wall retaining f o r c e
w h i c h " d o e s not differ except by a v e r y s m a l l q u a n t i t y " f r o m that of the
c u r v e d f a i l u r e line, B e g . H e proceeded to establish the f o r c e s y s t e m on the
s o i l m a s s bounded by a h o r i z o n t a l s u r f a c e , the s t r a i g h t internal failure plane
and a s m o o t h v e r t i c a l w a l l . F i g . 6.4(b) s h o w s the f o r c e s y s t e m , using present
day t e r m i n o l o g y for the force v e c t o r s .
LATERAL EARTH PRESSURES 141
^ 7 x ( z - xtan</>) - c ( x 2
+ z ) 2
Ρ = - (6.1)
χ + ztan</>
( γ ζ ί 8 η φ + c ) ( z - 2xztan</>- χ )
2 2
d£ = = Ο (6.2)
d x
(χ + ztan</>) 2
χ = z(y}l + tan 0 - t a n $
2
(6.3)
z/x = tan/3 = (V 1
+ tan 0 2
- tan</>) _1
= tan(7T/4 + φ/2) (6.4)
ρ = Ι γ ζ 2Γ]^4ηό1_ J l ^ i O 0 ( 6 . 5 )
K 2
c
= κ
a
= }- Si
1 + sin<p
. %
n
(6.7)
P =
z 0
h
/P d z
= [ > z K
2
a - οζΚ ]" ε = |7(h -z2)K
2
a - c(h-z )K
0 c
2 2 Z
°
= ! γ η
2
Κ - chK + c K c
(6.9)
2 ' a c 2γΚ 3
LATERAL EARTH PRESSURES 143
In the general case, retaining w a l l s are not s m o o t h nor are they necessarily
vertical. T h e y exhibit a c o e f f i c i e n t of f r i c t i o n , tano, in c o n t a c t w i t h the soil,
a s w e l l a possible adhesion s t r e n g t h , c . T h e s e c o m p l e x i t i e s , as well as a
a
Ρ = L J ± Ζ ZJ1 (6.10)
βίηία+β-δ-φ)
144 A G R I C U L T U R A L E N G I N E E R I N G SOIL MECHANICS
P = \yh K 2
a - chK c - c hKa c a + qhK a (6.11)
a " βΙη(α+/5-δ-φ)
ca " β1η(α + | 8 - δ - ψ ) ( }
dK / 3 = ο
a dj
sin(a - δ - φ)
be e m p l o y e d .
z
o = < c K
c + c
a ca " ^ 7 Κ
K
3 (6.16)
Ρ = I K ( 2 - 2) - cK (h - z ) - c K
a h z c 0 a c a (h - z ) 0 + qK (h - z )
a 0 (6.17)
ρ = W 1 + W 2 + Q - U(cot|3) - c h ) ( l + cot
2 2 βοοί[β-φ\)
sino + coso cot(/3- φ)
= Γ 30 + H 0 . 2 5 c o t f l _ 25q cot/3 c o t [ ^ ]
+ k N / m = ρ ρ
|_sino + c o s o c o t ( p - 0 ) s i n o + c o s o cotCjy-φΤJ a c
65 40.73 18.12
66 40.79 18.34
67 40.80 18.49
68 40.75 18.58
ρ (W + Q) - c h ( l + cotff cot[j8-<fl) Q D
P =
s i n £ + οοβδ c o t ( / j j - 0 ) = a " p p
c
7TPq
P r o b l e m 6.3. A n e x a m p l e s i t u a t i o n w h i c h c a n be a p p r o x i m a t e d as a strip
s u r c h a r g e load beside a retaining wall is s h o w n in F i g . 6.12. A 28 t tracked
tractor w i t h the d i m e n s i o n s noted rests on sand at a distance of 1 m from a
r e t a i n i n g w a l l of 3 m height. T h e t r a c t o r t r a c k s are 2.90 m long on the
g r o u n d . W h a t is the pressure distribution w i t h depth on the w a l l ?
E a c h t r a c k will be a p p r o x i m a t e d as a s t r i p l o a d . The a v e r a g e s u r c h a r g e
pressure is the w e i g h t on e a c h t r a c k divided by the t r a c k ground c o n t a c t area
as f o l l o w s .
q = (28tx9.8kN/t)/(2x0.6mx2.9m) = 78.9 k P a
Β /χ
1 1 = 0.60/1.00 = 0.60, B /X
2 2 = 0.60/3.20 = 0.19
152 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
p
ah = 7 a h
K z = 1 4
· 7 χ 0
- 2 4 : 5 ζ
kN/m 3
= 3.57 ζ kN/m 3
z, m z/x 1 P /2q
x
z/x 2
p /2q
2 Pi P 2 Pah
T o t a l p^, k P a
6.2. D E S I G N O F R E T A I N I N G WALLS
Τ = P s i n ( a - o) + U s i n a - c ^ c o t a = Η + U s i n a - c h c o t a a (6.19)
Ν = P c o s ( a - δ) + U c o s a + c h + W = V + U c o s a + c h + W a (6.20)
LATERAL EARTH PRESSURES 153
S.F. = ( N t a n 6 + c B ) / T a (6.21)
q
ult =
^ T 3 X N
7 + c N
c ( -24)
6
V t 2 Ύ 3 Χ Ν
Ύ + c N
c
b.h. - q = Ν / ( 3 χ ) (6.25J
Problem 6.4. F i g . 6.14 s h o w s a c o n c r e t e g r a v i t y wall having a height of 4 m,
r e t a i n i n g dry c o m p a c t sand w i t h the properties noted. W h a t is the s a f e t y
f a c t o r of this wall in sliding and t i p p i n g ?
P a =γΗ Κ /2 2
3 = 43.0 k N / m
Η = Ρ s i n ( a - δ) = 32.5 k N / m
V = Ρ cos(Q! - δ ) = 28.2 k N / m
LATERAL EARTH PRESSURES 155
Ν = W + V = 169.3 k N / m
and the s a f e t y f a c t o r in sliding;
q r N t a n o _ 169.3 t a n 2 7 ° k N / m 9
X = (Vx p - Hz p + Wx )/N c g
q
ult =
^Ύ3ΧΝ7= ^ 1 4 . 7 x 3 x 0 . 6 7 3 x 3 9 k P a = 579 k P a
S , F
* = q
ult « = / c 5 7 9
k P a / ( N / 3 X ) = 579 kPa/83.9 k P a = 6.9
156 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
C A
6.3. P A S S I V E S O I L RESISTANCE
Ρ = i 7 h 2
K + ZchyjK (6.27)
Ρ =iyh K 2
p + chK c + c hK
a c a + qhK q (6.28)
' ah
c
Impending
Mction
is _ ( c O t Q ? + C O t f f ) SJN(/3 + Φ) (fnq\
ρ " 2βίη(α+|8 + δ + φ ) ^
Κ = . ^ . , *Φ ^ 00
JN (630)
cβιηβ s i n f o +β + δ + </>)
ca s i n a S I N ( a +β + £ + φ)
K q = 2K p (6.32)
isinta+δ) 8ίη(δΤ^Γ. c o s ( a + δ + φ)
= * S m a s
' n
s i n ( a + 6 ^ ) ( & 3 3 )
Φ 0 δ° a 0
β° Κ
Ρ
Κ
c K
ca *Kp 1
*Kp 2
1
E x t r a p o l a t e d by the power law: K p = K (K^/K )
0 Q
φ
(6.34)
2
E x t r a p o l a t e d by linear proportion: K p = K Q + (Κφ - K )(8/(f>) Q (6.35)
(W x + W ) + ( C J L J + c L ) ( s i n / 3 + cos/3 cot [ β + φ ] )
2 2 2
P =
c o s ( a +o ) + sin(o: + o ) c o t ( j g + 0 ) P
p + P
c
7 2 . 7 7 ( c o t a + cotjg) + 2 4 ( 1 + cot/3 c o t [ β + φ ] ) /
c o s ( a + o) + s i n ( a + 6)cot( j8 + <p)
k f s J m
β° P , kN/m
p Ρ = P p + P ,c kN/m
24 252.9 341.9
25 252.5 340.3
26 252.6 339.4
27 253.1 339.2
28 254.0 339.6
Material γ, kN/m 3
00 δ°
Metal 2
Wood 3 Concrete
S i l a g e , 1 0 % m.c. 1
2.9 20 20 20 20
5 0 % m.c. 3.9 33 33 31 33
C o r n , 1 0 % m.c. 6.9 — 14 17 22
1 4 % m.c. 7.8 — 19 19 33
O a t s , 1 0 % m.c. 4.9 — 12 14 22
1 6 % m.c. 4.9 — 22 19 25
R i c e , 1 4 % m.c. 5.9 — 23 27 27
S o y b e a n s , 7 % m.c. 12.7 -- 12 17 21
1 2 % m.c. 13.7 — 12 24 29
W h e a t , 1 0 - 1 5 % m.c. 7.8 27 11 22 27
1 6 - 2 0 % m.c. 6.9 27 22 27 27
1
m.c. = m o i s t u r e content on a percent wet basis.
^ G a l v a n i z e d sheet s t e e l .
^ D o u g l a s fir or oak, a v e r a g e d .
4
S t e e l trowel f i n i s h .
LATERAL EARTH PRESSURES 163
p
h-a5-[i--- K h t a n o / R
] «*>
w h e r e P h = lateral bin wall pressure,
7 = m a t e r i a l bulk density,
Κ = ratio of lateral to v e r t i c a l pressure,
h = depth f r o m the top s u r f a c e of the m a t e r i a l .
R = 7 Γ γ / ( 2 7 Γ γ ) = r/2
2
f o r a c i r c l e of radius r, (6.37)
R = B / ( 4 B ) = B/4
2
f o r a square of side B, (6.38)
Κ 1 - sin0cos2C ( 6 Λ 0 )
1 + sin<p c o s 2 €
where 2 6 = sin 1
(sinδ /sin φ) - δ
p
h = γ Κ η ( 1 - Khtan6/D) f o r h/D < l/(2Ktan6) (6.41)
h = D / ( 2 K t a n £ ) = 4 m / ( 2 x 0 . 4 1 1 t a n l 7 ° ) = 15.9 m
h, m P ,kPa
n (Janssen) P , k P a (Lvin)
h
5 11.9 19.1
10 18.2 47.0
15 21.6 91.3
20 23.4 102.3
Κ = 1 - sin0 (6Λ3)
1 + sin<p
Surface
1ST
tano
t i m e s the slice height, dh. A n d the shear f o r c e s on the slice sides are the
horizontal forces t i m e s the c o e f f i c i e n t of f r i c t i o n , tano, between the m a t e r i a l s
inside and outisde the t r e n c h . T h e w e i g h t of soil in the slice dh is ' y B d h , per
u n i t l e n g t h of t r e n c h . Thus the s u m of v e r t i c a l f o r c e s a c t i n g on a slice is
equated to zero for equilibrium as f o l l o w s .
LATERAL EARTH PRESSURES 167
V + dV + 2 K V ( t a n 5 ) d h / B - V - 7 B d h = 0 (6.44)
dV + 2KV(tano)dh/B - 7 B d h = 0 (6.45)
V = A 1 + A e2
C h
(6.46)
V = ΎΒ 2
Γ _ -2K(tan5)h/Bl
e + Qe -2K(tan<5)h/B ( 6 Λ ? )
2Ktan5 L J
M a r s t o n and A n d e r s o n (1913) found that E q n . 6.47 w o r k s quite reliably to
p r e d i c t v e r t i c a l loads on pipes of d i a m e t e r s 457 and 9 1 4 m m buried in
t r e n c h e s of widths 683 and 1266 m m , r e s p e c t i v e l y , at depths of about 5 m
w i t h o u t a s u r c h a r g e load, Q.
W h e n a s u r c h a r g e load is of finite l e n g t h , s u c h as a vehicle wheel,
S p a n g l e r (1973) s u g g e s t e d that a B o u s s i n e s q type of pressure distribution c a n
be u s e d to design c o n s e r v a t i v e l y f o r the load c a r r i e d by a c e r t a i n length of
p i p e o r c o n d u i t , and he r e c o m m e n d e d that a length of 1 m be considered.
F i g . 6 . 2 2 s h o w s his theoretical c a l c u l a t i o n , w h i c h is consistently equal to or
h i g h e r than experimental m e a s u r e m e n t s of pipe loads f r o m s t a t i c wheels on
t h e s o i l s u r f a c e . A l s o in F i g . 6.22 is d r a w n the predicted c u r v e f r o m the
s e c o n d t e r m of E q n . 6.47 above w h i c h involves a line load, Q, using t y p i c a l
soil strength quantities. T h a t part, V q , of the v e r t i c a l load due to a line load
is larger than the B o u s s i n e s q r e c t a n g u l a r load on a pipe, because the latter
s p r e a d s out in t w o dimensions in the soil (ref. F i g . 4.7). T h u s the B o u s s i n e s q
c u r v e appears to be the logical design tool for this s i t u a t i o n .
6.6. P R O B L E M S
Dry sand
7 = 15.7 kN/m 3
0= 3 2 °
6m
120 = 28°
Bulldozer
blade, ^ ~ 1
3m \
Loose sand 30 cm average
7=12.7 kN/m 3
. . . . |
30 cm S a n d
' T= 5- kN/m3
1 7
φ = 35°
JL δ = 23.3°
A n s w e r : A b o u t 8 0 % of the silage w e i g h t .
Answer: F = 479 k N .
Α Δ Δ Δ
3m Gravel 7 = 15.7 kN/m 3
A Δ Δ Δ
J* 2m Tiebacks at 3 m spacing
3m
Water table
3m -2m
q=60kPa
SILT
7 = 15.7 kN/m3
c = 6 kPa
0 = 30°
δ-26°
7.1. W A T E R EROSION
L = (l/22m) x
(7.2)
S = (0.43 + 0.30s + 0 . 0 4 3 s ) / 6 . 5 7 4
2
(7.3)
C r o p practice C r o p s t a g e period
0 1 2 3 4
U p and d o w n 1.0
Parallel ( 1 2 % ) 0.8 - _
Problem 7 . 1 : A n example of the use of the universal soil loss equation f o r the
prediction of an a p p r o x i m a t e rate of soil erosion is as f o l l o w s . C o n s i d e r a
silage corn field near D e t r o i t , M i c h i g a n or W i n d s o r , O n t a r i o , on a sandy l o a m
field h a v i n g a slope of 2 % and a slope length of 100 m . T h e c u l t i v a t i o n
p r a c t i c e h a s been plowing up and d o w n the slope of the f i e l d . T h e r e l e v a n t
p a r a m e t e r s in this case a r e ;
R = 100 ( F i g . 7.1)
L = (100/22) · 0 3
= 1.57 ( E q n . 7.2)
Y = 11.8(Q χ q ) - p
0 5 6
Κ L S C Ρ (7.4)
where
1. Seasonal rainfall e n e r g y , Ε ( J / m ) , 2
Ζ = KCX (7.5)
0 1000 2000
Annual rainfall, mm
300
200
Specific
soil loss,
K, t/ha/y
100
0.8
0.6
Soil
loss
ratio
C
0.4
0.2
2.0|
1.5
05
0.5
" 025~
0 20 40 60 80 100
Slope length, L, m
Practice Factor
Cotton 1.0 62 55 41 24 9
2.5 84 75 56 32 13
4.5 95 84 63 36 14
C o w p e a s hay 1.5 54 62 67 55 35
2.5 65 76 81 67 43
" silage 2.0 26 29 32 29 20
6.0 43 50 54 49 34
10.0 61 70 83 69 48
11
green crop 6.0 5 17 26 30 22
10.0 7 24 35 41 36
Groundnuts 0.8 50 57 49 32 16
1.5 64 73 62 41 20
3.0 71 80 69 45 22
Maize 2.0 24 29 28 21 12
6.0 43 52 51 37 21
10.0 55 68 66 48 27
Rice 1.5 54 56 53 40 24
2.5 69 72 67 51 31
4.5 84 88 82 62 37
Sorghum 1.0 39 40 35 25 12
2.0 47 49 43 30 15
6.0 81 84 73 52 26
Soybeans 0.8 34 38 36 27 14
1.5 51 57 54 39 20
2.5 69 77 73 53 28
4.0 80 90 85 62 32
Sunflowers 0.1 5 8 11 9 5
0.4 11 19 24 21 12
1.3 28 50 62 55 30
Tobacco 1.0 3 6 8 6 3
2.0 11 21 26 22 11
3.0 19 36 46 38 20
4.0 27 52 66 54 28
Grass ley 97 89 68 42
Poor ley 65 55 33 14
Weed fallow 0 2 10 32
Poor fallow 0 1 5 14
SOIL E R O S I O N A N D P R O T E C T I O N 185
F r o m Table 7.4, the erodibility index for the soil type is found to be 4.0.
F r o m Fig. 7.4, the bare soil erosion is e s t i m a t e d to be about 2 1 t/(ha y). T h u s ,
including the cropping and topography f a c t o r s , C and X , the soil loss e s t i m a t e
w i t h o u t ridges would be:
elevation contours
Grass or legume
buffer strips
Channel
£2:1 Cut Fill
_ % slope
6
1
0 15
1> X, m
30
(a) Broadbase
Embankment
2% slo pe Level
—^ X, m
30 60 90
_ Γ 100 L S "Ι 2
H I
00)
™ -
[ 0.75 + 0.53 S 1 + 0.076 S^J
7.6. W I N D EROSION
7.7. PROBLEMS
7.8. A r a i n s t o r m o c c u r s in a 4 k m a g r i c u l t u r a l area of N e w Y o r k S t a t e
2
A n s w e r : Y = 8 t/ha.
8.1. C U T T I N G FORCES
Ρ = (γη Κ 2
ρ + chK c + c hK
a c a + qhK )w
q (8.1)
Η = P s i n ( a + δ) + c d w c o t a (8.2)
a
V = Pcos(a + δ) - c dw + W a (8.3)
/sinfo+ δ ) 8 ΐ η ( δ Τ ^ _ ϋ 0 8 ( α + δ + φ)"
β = cot" 1
I sinasino = 3 1 6 ? o
μ
sinTa + 0+φ)
ρ 2sin(a + β + δ + φ )
Κ = • cos φ . =
c sin^Q s i n ( a +β+ δ + φ)
κ = . -οο (α + 0 + φ )
8 = 0 β 7 3
ca sin α s i n t o + ρ + δ +
Φ)
Ρ = (γη Κ 2
ρ + c hK c + c hKa c a )w = 70.1 k N
Η = P s i n t o + δ ) + c h w c o t a = 72.8 k N
a
V = P c o s ( a + δ ) - c h w + W = 52.0 a kN
P r o b l e m 8.2.
A n o t h e r example is shown in F i g . 8.3 involving a s t r u c t u r e . A n anchor
made of concrete (specific g r a v i t y 2.4) is designed w i t h the d i m e n s i o n s s h o w n
t o hold down a cable with a design tension of 450 k N . T h e soil has been
compacted around the anchor to have properties γ = 1 9 . 6 k N / m , φ = 3 5 ° , δ = 2 7 ° , 3
d e s i g n , t h e m a x i m u m s t r e n g t h of the a n c h o r - s o i l c o m b i n a t i o n m u s t be
e s t i m a t e d for horizontal m o v e m e n t and v e r t i c a l uplift.
SOIL C U T T I N G A N D T I L L A G E 195
450 kN
Ρ = (γη Κ 2
ρ + c hK c + c hK
a c a )w = 817 kN
H m = Psin(a + 6) + C g d w c o t a = 838 k N
v
m = P c o s ( a + δ) - c d w + W = 369 k N
g
Circular
side
crescents
of w i d t h , w.
ρ =
[(i d qd)(cota.cot/3)sinfa 6)
L 2 _ _ _
T
2
+
sinp
+ + cd^?4- c
a
d ^ ^ ^ l w
sing J (8.5)
1 s i n ( a + j 3+ o φ ) +
SOIL C U T T I N G A N D T I L L A G E 199
dP Q =
[d7dr +2
iqr )sin(a
2
+ δ) +
2
W ^ l d p
sinff J ^ ( Q £ )
2 sin(a+|8 + δ + φ )
2 =
l^2 P
C0S
_
\(\ydr
L Ρ
2
+ i q r ) s i n ( a δ) +
£
2
+ ^ r d ^ l s i n
2 sinff J 1
p'
/q -y\
" 8ίη(α+/5 + δ φ ) +
ρ ΓβΎ^Μτ,)
L -2
l
d 3w
+
^ d v
w J βΐηΟα+δ) + c d cs a
i nsp& i + iw) + c ad c o
s if neg ^ l Jw
s
s i n ( a +β+ δ + φ)
= ( 7 ^ Ν / + c d Nc + q d Nq + c ad N cJa w
2
Ύ n (8.8)
k c o t Q f + cotfi) 1 + | ^ \ c o t j 3 + 2 c o t a c o t j 8 s i n ( a + δ )
2
Ny = ± L -Y Ε l\ n =ί (8.9)
ι sinia+β + δ + φ)
I t c a n be seen in E q n . 8.9 that the slenderness r a t i o , d/w, is an i m p l i c i t
component of the Ν γ f a c t o r , and will thus a f f e c t the f o r c e per unit w i d t h on
a cutting blade. A l s o , if the principle used in C h a p t e r 6 is e m p l o y e d , to find
t h e m o s t likely failure plane by m i n i m i z i n g the Ny t e r m , then the failure
a n g l e , β, will be influenced also by the blade depth to width r a t i o . These
e f f e c t s are in a g r e e m e n t w i t h the p h y s i c a l o b s e r v a t i o n s of P a y n e (1956).
200 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
δ°
a 0
d/w β
^cr
0
Νγ N
c N ca O Q
~0 5 10 15 20
w, cm
6.35err"If" ~^L
Ο ΓΊ
30cm a = 60"
pom
20cm
40 cm
Fig. 8.10. Cross-section of (a) one deep pass and (b) the first of two passes
of a set of chisel plows, P r o b l e m 8.3.
T h e c a l c u l a t i o n s proceed as f o l l o w s .
Η = P s i n ( a + 6 ) + c d w c o t a = 71.1 k N
a
SOIL C U T T I N G A N D T I L L A G E 205
where d = 0.20 m
and q = 0.20 m χ 17.6 k N / m 3 = 3.52 k P a
Ρ = 20.43 k N
H = P s i n t a + δ ) + c d w c o t a = 20.40 k N
2 a
A n o t h e r e x a m p l e of the t h e o r e t i c a l w o r k w h i c h c a n be a c c o m p l i s h e d w i t h
t h e m e c h a n i c a l m o d e l o f passive soil failure by n a r r o w blades is the
o p t i m i z a t i o n of tool design for p a r t i c u l a r purposes, as illustrated by the
p r o b l e m s o l v e d below.
F i g . 8 . 1 1 . C h i s e l plow e x a m p l e in P r o b l e m 8.4.
206 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
a0
Ny Ν Η, Ν
c N
ca
15 17.0 9.6 10.0 932
20 12.0 9.7 5.0 892
25 10.5 9.8 2.6 903
30 9.1 10.2 2.1 970
q' = C N - C + P O N , Q = C N , C + Y Z K O ^ (8.10)
Ν
Ό = cot
* [[π^]« 2 ( 7 Γ / 2 +0 ) t a n
^ - l] (822)
Q = [cN' ( c d - d ) + c |K N» (d0 q
2
- d c
2
)] w (8.13)
208 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
For the force Ρ on that part of a blade above the c r i t i c a l depth, the three
dimensional model of S e c t i o n 8.1.2 is used. T h u s the total h o r i z o n t a l draft
r e q u i r e d to m o v e the i m p l e m e n t is the s u m of Q below the c r i t i c a l depth,
and H-^ above it as s h o w n below.
H = Q + H 1 = Q + P s i n ( a + 6) + c ^ w c o t a (8.14)
P r o b l e m 8.5.
Suppose that a steel flat blade w i t h a w i d t h of 5 c m operates at 100 c m
d e p t h in a l o a m soil w i t h properties Y = 1 9 . 6 k N / m 3 , c = 3 0 k P a , c = 5 k P a , 0=35°
a
a° d ,cm
c Q,MM d /w
c Ny N Q NQa H^,kN Η = Η +0,Ι<Ν
χ
100 cm '
(a)
8.2. V O L U M E O F SOIL C U T
A = (w + s)d (8.15)
soil is not lifted appreciably, then d should be used in E q n . 8.15 rather than
c
SOIL C U T T I N G A N D T I L L A G E 213
V o l u m e cut _ AX _A / P 1 ^
E n e r g y input " HX "H
.10r 1 1 1 1 1 r
d/w
h^-Spacing
W i n g e d subsoiler
The measured results are quite close to the theoretical predictions for the
20 blade rake angle, but are generally lower than predictions for the 3 5 ° t o o l .
T h e trends of c h a n g e s in v o l u m e c u t t i n g e f f i c i e n c y w i t h v a r y i n g r a k e angle
a n d slenderness ratio predicted by the m e c h a n i c a l model are c o n f i r m e d by
these e x p e r i m e n t a l field results in any e v e n t .
A n o t h e r possibility for i m p r o v i n g soil c u t t i n g e f f i c i e n c y is to arrange
implements to follow one another at different depths. Spoor and G o d w i n (1978)
t e s t e d s u c h a s y s t e m in a friable clay s o i l . A s p i c t u r e d in F i g . 8.19, the
c o m b i n a t i o n c o m p r i s e d two chisel plow s h a n k s f o l l o w e d by a deeper w i n g e d
subsoiler. T h e subsoiler had a total w i d t h of 30 c m , w i n g s at a r a k e angle of
2 2 ° and a depth of operation of 40 c m . T h e c h i s e l s were 8 c m wide w i t h a tip
r a k e angle of 2 0 ° and were varied in both depth and s p a c i n g as indicated in
T a b l e 8.4 below. There w a s no g r e a t d i f f e r e n c e in the d r a f t r e q u i r e m e n t of
t h e s y s t e m no m a t t e r w h a t the chisel depths or s p a c i n g . H o w e v e r , the area
d i s t u r b e d by the tool c o m b i n a t i o n increased m a r k e d l y w i t h both chisel depth
and s p a c i n g . The best c o n f i g u r a t i o n f r o m the point of view of e f f i c i e n c y w a s
t h a t having the widest chisel spacing and depth, and w a s about 7 7 % m o r e
e f f i c i e n t in the v o l u m e of soil cut per unit d r a f t than the subsoiler shank
operated alone.
Problem 8.7. A flat chisel blade is 63.5 m m wide and operates 200 m m deep
in a clay loam field having γ = 15.7 k N / m , 0 = 3 0 ° , δ = 2 3 . 3 ° , c o h e s i o n 20 k P a
3
s = dVcot j8 + 2cotacotj8
2
c a ca a
a 0
β° Νγ Ν N
ca
A=(w+s)d, m 2
H, k N H/A, kPa
c
30 61 9.21 10.46 2.25 0.0724 2.69 37.2
60 38.5 10.40 21.66 1.95 0.0824 3.34 40.6
8.3. S O I L LOOSENING
Material γ 0 , kN/m 3
Y f , kN/m 3
Swell factor %
8.4. PROBLEMS
8.9. For the same chisel plow and soil as in P r o b l e m 8.8, how deep c a n the
i m p l e m e n t be operated before the c r i t i c a l depth p h e n o m e n o n will be
encountered?
A n s w e r : d = 0.55 m .
c = 5 k P a , φ = 3 0 ° and δ = 2 0 ° ?
a
A n s w e r : (a) a= 2 7 ° , (b) QC = 9 ° .
220 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
A n s w e r : d = 0.65 m .
F i g . 8 . 2 1 . E x c a v a t i n g b u c k e t in P r o b l e m 8 . 1 1 .
A n s w e r : 15 c h i s e l s .
A n s w e r : (a) A = 0 . 1 9 m , H / A = 3 6 . 6 k P a , (b) A = 0 . 2 1 m , H / A = 7 1 . 6
2 2
kPa.
A n s w e r : d = 2 c m ; Η = 6.5 k N .
SOIL C U T T I N G A N D T I L L A G E 221
9.1. C O M P A C T I O N F O R E A R T H W O R K CONSTRUCTION
1. P r o c t o r (1948)
2. T a y l o r (1948)
3. G r a e c a n (1960)
4. C h a n c e l l o r and S c h m i d t (1962)
5. R o w l e s (1948)
6. L a m b e and W h i t m a n (1979)
1.81 1 1 r
I ι 1 1 1 1
0 5 10 15 20 25
Water content %
F i g . 9 . 1 . C o m p a c t i o n c u r v e results for P r o b l e m 9 . 1 .
224 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
Pdry = Po +
A
^ [ M 1
+
S % ) / p ] + Β log(w%)
Q (9.1)
T h e c o n s t a n t Β is a m e a s u r e of the s e n s i t i v i t y of the c o m p a c t i b i l i t y of a
soil to increases in m o i s t u r e content up to the o p t i m u m w a t e r c o n t e n t . B e y o n d
t h e o p t i m u m m o i s t u r e c o n t e n t , Β will have a n e g a t i v e v a l u e . E x a m p l e s of
s o m e of the c o n s t a n t s in E q n . 9.1 are g i v e n in T a b l e 9.2 below.
Soil p , t/m3
Q p , kPa
Q A , t/m3 B, t / m 3 Source
1 . R a g h a v a n et al., (1977a, b)
2. Soehne (1958)
3. H o v a n e s i a n (1958)
Pdry " P o = [ ° - 1 5
l o g ( N x 6 0 / 7 ) + 0.3 l o g ( 2 5 ) ] t / m 3
U
100 50 20 10 5 2 1 0.5 0.2 O.i 0.05 0.02 0.01 0.001
Grain size, mm
ο.
100 10 1 0.1 0.01 0.001
Grain size, mm
Tired roller
Pressure spread
9.2. C O M P A C T I O N O F A G R I C U L T U R A L SOILS
(9.2)
Υ
* ~ Y = C (
Pdry-P*dry) 2
15
Pdry - t/m 3
p d r y = 1.12 t/m 3
9.3. PROBLEMS
Answer: w = 1 1 % . R a n g e = 10 to 1 4 % .
0.15 t / m 3
a n d Β = 0.3 t / m . F o r the p r o d u c t i o n of g r a i n c o r n in an
3
on the two rear wheels together, and each rear tire has a field c o n t a c t
a r e a o f 0.45 m . If the a v e r a g e topsoil m o i s t u r e content is 2 0 % by
2
Chapter 1 0 . Geotextiles
10.1. U S E S O F GEOTEXTILES
1. S e p a r a t i o n of soil m a t e r i a l s ,
2. L a t e r a l retention of m a t e r i a l l a y e r s ,
3. S p r e a d i n g of loads in the subsoil,
4. D r a i n a g e envelopes and
5. Erosion protection.
T h e p r o p e r t y of g e o t e x t i l e s w h i c h a l l o w s t h e m to p e r f o r m the above
f u n c t i o n s is their tensile s t r e n g t h in t w o d i m e n s i o n s , a property w h i c h is not
s h a r e d by g r a n u l a r soil m a t e r i a l s . W h e n , for e x a m p l e , large p a r t i c l e s rest on
a layer of a finer grained s o i l , especially a w e t soil w i t h large pore pressures,
individual g r a v e l or sand grains c a n sink into the finer soil thus d e s t r o y i n g the
original c o n f i g u r a t i o n and f u n c t i o n of the e a r t h w o r k d e s i g n . F i g . 1 0 . 1 , for
i n s t a n c e , s h o w s the installation of a sheet of g e o t e x t i l e between l a y e r s of
g r a v e l a n d w e t f i n e g r a i n e d s o i l . T h e g e o t e x t i l e allows w a t e r to be
t r a n s m i t t e d t h r o u g h the boundary, thus p e r m i t t i n g drainage and the relief of
s o i l pore pressure. H o w e v e r , by virtue of its tensile s t r e n g t h and s i z e of
openings, the geotextile does not p e r m i t the e n t r a n c e of g r a v e l particles down
i n t o t h e s u b s o i l , or pumping of w e t f i n e s upwards into the g r a v e l l a y e r .
T h e r e f o r e , the g r a v e l layer r e m a i n s i n t a c t as a r o a d b a s e , f o u n d a t i o n subbase
or w h a t e v e r f u n c t i o n d e s i r e d . I n this m a n n e r , a g e o t e x t i l e c a n replace soil
filter m a t e r i a l s .
Tension
Pressure distribution
A n o t h e r i m p o r t a n t a p p l i c a t i o n o f g e o t e x t i l e s is in the drainage of
e a r t h w o r k s , especially r o a d s . R a t h e r than using open ditches if space is at a
p r e m i u m , s u b s u r f a c e drains c a n be used t o m a i n t a i n road bases at l o w w a t e r
contents throughout the y e a r , a n d hence at higher s t r e n g t h l e v e l s . E x c a v a t i o n s
a r e m a d e t o the sides o f a r o a d a t the d e s i r e d depth a n d longitudinal slope,
b e t w e e n 30 and 60 c m in w i d t h depending o n the e x p e c t e d v o l u m e of w a t e r
f l o w . A g e o t e x t i l e is then used t o line the e x c a v a t e d d i t c h , a n d permeable
granular fill is placed to the desired depth o n the g e o t e x t i l e . I f f a s t e r w a t e r
f l o w i s needed in a p a r t i c u l a r a p p l i c a t i o n , then a p e r f o r a t e d drain tile c a n
a l s o be installed in the e x c a v a t i o n a n d s u r r o u n d e d by g r a n u l a r m a t e r i a l . T h e
g r a n u l a r m a t e r i a l s e r v e s t o i m p r o v e h y d r a u l i c c o n d u c t i v i t y and the
p e r f o r m a n c e of the d r a i n , thus e n h a n c i n g the net drainage e f f e c t . F i g . 10.4
s h o w s both o f these options. T h e g e o t e x t i l e is then wrapped around the t o p
of the g r a n u l a r fill m a t e r i a l a n d the s u r f a c e topped up w i t h e x i s t i n g s o i l .
Soil surface
Perforated drain
Fig. 10.4. G e o t e x t i l e s used t o w r a p r o c k drains w i t h or w i t h o u t a drain pipe.
T h e u s e o f a geotextile in this d r a i n a g e a p p l i c a t i o n is p r i m a r i l y as a
f i l t e r i n g m a t e r i a l w h i c h p r e v e n t s fine p a r t i c l e s o f the s u r r o u n d i n g soil f r o m
entering a n d blocking the g r a n u l a r drain f i l l . A g a i n i t s tensile s t r e n g t h allows
t h e g e o t e x t i l e continuity t o be m a i n t a i n e d , y e t a high p e r m e a b i l i t y t o w a t e r
flow h a s a m i n i m a l d e t r i m e n t a l e f f e c t o n the r a t e of w a t e r inflow t o drains.
Some g e o t e x t i l e s are used d i r e c t l y o n drain tiles, e s p e c i a l l y in c a s e s where a
large quantity of tile is p l a c e d . T h i s is p a r t i c u l a r l y true in the c a s e of
s u b s u r f a c e drainage of a g r i c u l t u r a l fields w h e r e up t o 1 0 0 0 m o f s u b s u r f a c e
drains m a y be installed per h e c t a r e of s u r f a c e a r e a . I n these applications, the
a d d i t i o n o f g r a v e l e n v e l o p e s around drains is e c o n o m i c a l l y p r o h i b i t i v e .
H o w e v e r , the l a c k o f g r a v e l renders the entry a r e a f o r w a t e r to drains
considerably s m a l l e r , a n d the flow o f w a t e r longitudinally t h r o u g h g e o t e x t i l e s
t o e n t r y holes in drain tiles m o r e c r i t i c a l . R e s e a r c h projects on these
p r o b l e m s are ongoing t o d a y , and the reader is r e f e r r e d t o r e c e n t articles in
a g r i c u l t u r a l engineering r e s e a r c h journals f o r the l a t e s t d e v e l o p m e n t s .
236 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
10.2. D E S I G N O F G E O T E X T I L E APPLICATIONS
Geotextile.
L+0.8t
t is thus;
q
o =
(B + 0.8t)(L + 0.8t) ( m )
cN (B+0.8t)(L+0.8t)
c 6.28c(B+0.8t)(L+0.8t)
Q
design = SjT. = sF. < ·)10 2
U s i n g E q n . 7.2;
50 k N = 6.28 χ 15 k P a (0.3 m + 0 . 8 t ) / 3 2
t = 1.2 m
q
ult = c N
c ( B +
°- 8 t
V B
= 9 k P a χ 5.14 χ (0.15 + 0.08)/0.15
= 70.9 k P a
GEOTEXTILES 239
Penetration Ζ , cm
Granular fill
V
Geotextile Τ 5 o f t s o j|
<*xx = K
a 7 * (10-3)
The total shear f o r c e per unit length of slope is the integral of the above
expression over the height Η of the slope.
Τ = ^ K a 7 H 2
(10.4)
Κ = ] - ? ^ =
s n
1 - ? ^ ° = 0.271
s n
a 1 + sin0 1 + sm35 u
Τ = K j H 2
/ 2 = 0.271 χ 17.6 k N / m 3
x(3 m ) / 2 = 21.5 k N / m
2
10 15
H = 3m
D
Δ Η
m
10.3. PROBLEMS
A n s w e r : Q = 42.5 k N .
A n s w e r : Q = 54.6 k N .
A n s w e r : Τ = 44.4 k N / m .
GEOTEXTILES 243
11.1. F R O S T PENETRATION
— Temperature, °C
(HI)
(1L2)
%t = a
fez (1L3)
Τ = -T f + (T /z )z
f f (11.4)
Unfrozen soil
f r o s t l i n e of 0 ° C m o v e s d o w n w a r d s v e r y s l o w l y in s u b z e r o air conditions.
Eqn. 11.1 g i v e s the rate at w h i c h heat is t r a n s f e r r e d t h r o u g h soil for a g i v e n
temperature g r a d i e n t . The quantity of heat w h i c h m u s t be m o v e d to a d v a n c e
the frost front in depth is p r i m a r i l y the latent heat of fusion, L , of the pore s
T f d Z f
- k
f — = L
s — (11.5)
z f dt
Ζ 2
!• t
f
f Γ ,
— = — -T dt
f (11.7)
2 L ·'
SOIL F R E E Z I N G 247
40
30
Water
content,
% dry
matter
20
10
0 4 6 8 10 12 14
Thermal conductivity, k , kJ/(m-h-°C)
f
0
F i g . 11.3. T h e r m a l c o n d u c t i v i t y , k f , o f f r o z e n sand or silt and clay as a
function of d r y density and w a t e r content (converted f r o m A l d r i c h ,
1956).
T h e i n t e g r a l of s u b z e r o t e m p e r a t u r e w i t h t i m e is c o m m o n l y called the
f r e e z i n g index, F, a n d is usually e x p r e s s e d in d e g r e e - d a y s . M a p s o f a v e r a g e
a n n u a l f r e e z i n g indices a r e available f o r m a n y p a r t s o f the w o r l d , a n d t h e
f r e e z i n g index f o r a p a r t i c u l a r l o c a t i o n c a n be used in the expression below
w h i c h s o l v e s f o r frost depth z^.
z f = (11.8)
F r o m F i g . 11.3, the f r o z e n t h e r m a l c o n d u c t i v i t y , k , c a n be e s t i m a t e d as
f
V 2 75 m =
7 7 M
Z
f " 113,800 k J ( ° C ) m h X
daT" =
11.2. F R O S T H E A V E PRESSURE
4
%
expansion
~0 10 20 30 40
Water content %
s u r f a c e (Penner, 1959).
L
s A
T f
(11.9)
p
" J
K - V
where L s= the latent heat of fusion of w a t e r ,
Τ = ambient absolute t e m p e r a t u r e ,
v =
w the s p e c i f i c v o l u m e of w a t e r and
v, = the s p e c i f i c v o l u m e of ice per unit m a s s .
334.8 M J / t A T f
ρ = - L
= - 13.6 A T f MPa/o (11.10)
2 7 3 ° χ 0.09 m /t
3
Τ Saturated soil
1m
7 ι^ / / / y y
ν y y y y
Limestone
A h = 3 . 5 % χ 1 m = 3.5 c m
11.2.2. I c e l e n s i n g .
F i g . 11.4(b) s h o w s s c h e m a t i c a l l y a c a s e of an ice lens separating itself
f r o m the soil. T h i s phenomenon o c c u r s essentially at the frost front in
s i t u a t i o n s w h e r e there is a plentiful w a t e r supply w i t h i n the soil, and
s u f f i c i e n t t i m e for the w a t e r to m o v e to the ice lens. W h e n a soil is nearly
s a t u r a t e d , w a t e r does tend to m o v e naturally f r o m w a r m e r depths to colder
p l a c e s , because of the higher free energy in w a r m e r w a t e r . P o r e w a t e r will
move faster, of c o u r s e , in a c o a r s e g r a i n e d soil w i t h a high p e r m e a b i l i t y , and
if there is an upwards hydraulic gradient (artesian conditions).
When the f r o s t front is s t a t i o n a r y for s o m e t i m e , the ice lens will continue
to grow upwards f r o m the f r e e z i n g depth and lift the soil a b o v e , provided that
t h e w a t e r supply r e m a i n s in e x i s t e n c e . T h e r e a s o n for w h i c h the ice lens
grows is that portions of the lower s u r f a c e of the ice are at a higher f r e e z i n g
p o i n t than the local t e m p e r a t u r e itself. T h e s e areas are those of n e g a t i v e
c u r v a t u r e , s h o w n as points A in F i g . 11.7, because the i c e - w a t e r s u r f a c e
tension g i v e s rise to a r e d u c t i o n in pressure on the c u r v e d ice s u r f a c e .
On the other hand, ice l o c a t i o n s w i t h a positive c u r v a t u r e (B in F i g . 11.7)
h a v e a r e l a t i v e pressure increase on t h e m due to i c e - w a t e r s u r f a c e tension,
and their f r e e z i n g point is lower than the e x i s t i n g t e m p e r a t u r e . T h e r e f o r e ,
f r e e z i n g will not o c c u r at these latter p l a c e s , and the ice front does not
advance in the s o i l . H o w e v e r , the f r e e z i n g of w a t e r continues at points A as
w a t e r is supplied f r o m below, and the ice lens g r o w s . G i v e n s u f f i c i e n t t i m e
a n d a h i g h soil hydraulic c o n d u c t i v i t y , the g r o w t h of the ice lens can be
considerable, e v e n at rates up to s e v e r a l m m / h .
The potential heaving pressure is not as g r e a t as that of w a t e r at several
degrees below f r e e z i n g , h o w e v e r , since the f r e e z i n g point depression required
to s t o p g r o w t h of the ice lens is s m a l l at points A in F i g . 11.7. In f a c t the
e f f e c t i v e n e g a t i v e pressure, p, a c t i n g on points A c a n be c a l c u l a t e d as:
252 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
ρ = __20iw_
Ρ
where a j = the i c e - w a t e r s u r f a c e tension (energy per unit area), and
w
% = - J=s (11.12)
dT f jT V
ΔΤ = 20jwTvi ( 1 L 1 3 )
* r L
Ρ 3
L i k e w i s e , the f r e e z i n g point depression for positively c u r v e d s e c t i o n s Β in
F i g . 11.7 is;
ΔΤ = - 2 a
»wTvi 01.14)
f r v Ls L
SOIL F R E E Z I N G 253
ΔΤ = ζσίν,τνιΓι. + U ( 1 L 1 5 )
s
L p
r
vJ
r
ρ = 4.84 χ 1 0 ' 4
+ j L J kPa m (11.17)
A /Ac = (1 - η) 213
(11.18)
'V-P[TV! 1 / 3 m
- i9)
lens in a u n i f o r m soil.
p f = 4.84 χ I P ' 4
|\ + (L^JI -r,] 2 / 3
kPa m (11.20)
The number 2 in the bracket above s h o w s that the particle and void radii
are essentially equal in this c a s e , and 0.63 is the r a t i o of soil area c o n t a c t i n g
the ice lens to the total c r o s s sectional area for this void r a t i o .
ρ = -u = 2 ^ § w ( 1 L 2 2 )
11.3. PROBLEMS
Answer: Heave = 11 c m .
Answer: w = 1 5 % .
Answer: p f = 32 k P a .
References
C A T E R P I L L A R C O . (1981). H a n d b o o k of e a r t h m o v i n g . C a t e r p i l l a r T r a c t o r C o . ,
P e o r i a , I L , 24 pp.
C A U G H Y , R.A., T O O L E S , C . W . and S C H E E R , A . C . (1951). L a t e r a l and v e r t i c a l
p r e s s u r e s of granular m a t e r i a l in deep bins. B u l l . 172, I o w a E n g i n e e r i n g
Experiment Station, A m e s I A .
C H A N C E L L O R , W . J . and S C H M I D T , R . H . (1962). A study of soil d e f o r m a t i o n
beneath surface loads. T r a n s . A m e r . S o c . A g r i c . E n g r s . , V o l . 5, pp. 2 4 0 - 2 4 6 ,
249.
C O R P S O F E N G I N E E R S . (1944). C o n s t r u c t i o n and r o u t e s of c o m m u n i c a t i o n .
W a r D e p a r t m e n t F i e l d M a n u a l N o . F M 5-10, U.S. G o v ' t . P r i n t i n g O f f . ,
W a s h i n g t o n D C , 596 pp.
C O U L O M B , C A . (1776). E s s a i sur une application des r e g i e s des m a x i m i s et
m i n i m i s a q u e l q u e s p r o b l e m e s de statique r e l a t i f s a l'architecture.
Academie royale des S c i e n c e s : M e m o i r e s de M a t h e m a t i q u e et de P h y s i q u e ,
presentes a l ' A c a d e m i e royale des S c i e n c e s , par divers s a v a n t s , et lus dans
les A s s e m b l i e s , P a r i s , V o l . 7, pp. 3 4 3 - 3 8 2 + 2 plates.
D A R C Y , H . (1856). L e s fontaines publiques de la ville de D i j o n . D a l m o n t ,
P a r i s , 674 pp.
D E S I R , F . L . (1981). A field e v a l u a t i o n of the wedge a p p r o a c h to the analysis
of soil c u t t i n g by n a r r o w blades. M . S c . T h e s i s , M c G i l l U n i v . , D e p t . A g r i c .
E n g . , M o n t r e a l , Q C . 245 pp.
E L W E L L , H . A . (1978). M o d e l l i n g soil losses in southern A f r i c a . J . A g r i c . E n g .
R e s . , V o l . 2 3 , pp. 1 1 7 - 1 2 7 .
E L W E L L , H . A . (1980). D e s i g n of s a f e r o t a t i o n a l s y s t e m s . D e p t . C o n s , and
E x t . R e p o r t , M i n . A g r . , H a r a r e , Z i m b a b w e , 50 pp.
E L W E L L , H . A . (1984). Sheet erosion f r o m arable lands in Z i m b a b w e : prediction
a n d c o n t r o l . P r o c . S y m p . C h a l l e n g e s in A f r i c a n H y d r o l o g y and W a t e r
R e s o u r c e s , H a r a r e , I A H S P u b l . N o . 144, pp. 4 2 9 - 4 3 8 .
E V E R T S , R., V A N Z A N T E N , D . C . and R I C H A R D S , P . C . (1977). B u n k e r d e s i g n .
Part 4: R e c o m m e n d a t i o n s . T r a n s . A m e r . S o c . M e c h . E n g r s . J . E n g . Industry,
V o l . 9 9 , pp. 8 2 4 - 8 2 7 .
F A N , T a i - L i n (1985). M u l t i - p l a t e p e n e t r a t i o n tests to d e t e r m i n e soil s t i f f n e s s
m o d u l i . M . S c . T h e s i s , M c G i l l U n i v . , D e p t . A g r i c . E n g . , M o n t r e a l , Q C . 99 pp.
F O R T I E R , S. and S C O B E Y , F . C . (1926). P e r m i s s i b l e canal v e l o c i t i e s . T r a n s .
A m . S o c . of C i v i l E n g . , V o l . 8 9 , pp. 9 4 0 - 9 8 4 .
G A M E D A , S. J U T R A S , P . J . and B R O U G H T O N , R . S . (1983). O c h r e in s u b s u r f a c e
drains in a Q u e b e c fine sandy l o a m . C a n . A g r i c . E n g . , V o l . 25, N o . 2, pp.
209-213.
G I L L , W . R . and M c C R E E R Y , W . F . (1960). R e l a t i o n of size of cut to tillage
tool e f f i c i e n c y . A g r i c . E n g . , V o l . 4 1 , pp. 3 7 2 - 3 7 4 , 3 8 1 .
G O D W I N , R . J . (1974). A n i n v e s t i g a t i o n into the m e c h a n i c s of narrow tines in
f r i c t i o n a l soils. P h . D . T h e s i s , U n i v . of R e a d i n g , 248 pp.
G O D W I N , R . J . and S P O O R , G . (1977). Soil failure w i t h n a r r o w tines. J . A g r i c .
E n g . R e s . , V o l . 22, pp. 2 1 3 - 2 2 8 .
REFERENCES 261
G O R I A T C H K I N , V . P . (1937). K o l e s a z h a t v e n n i h m a s h i n . S o b r a n i e s o c h i n e n i i ,
S e l k h o z g i z , V o l . 2 and 4, M o s c o w .
G O U R C , J . P . ( 1 9 8 3 ) E t u d e s de c a s ; c h a u s s e e s p r o v i s o i r e s , r e m b l a i s ,
s o u t e n e m e n t s . D a n s C o u r s intensifs - les G e o t e x t i l e s . f c o l e polytechnique,
C e n t r e de C o u r s intensifs, M o n t r e a l , Q C .
G R A E C E N , E . L . (1960). W a t e r content and soil s t r e n g t h . J . Soil S c i . , V o l . 1 1 ,
pp. 3 1 3 - 3 3 3 .
H A M M A M J I , Y . (1969). S o m e f a c t o r s a f f e c t i n g heave pressures of f r o z e n soils.
M.Sc. Thesis, McGill University, Montreal, Q C .
H A N S E N , J . B R I N C H (1957). F o u n d a t i o n of s t r u c t u r e s - general r e p o r t . P r o c .
4 t h I n t . C o n f . Soil M e c h . and F o u n d . E n g . , L o n d o n , V o l . I I , pp. 4 4 1 - 4 4 7 .
H A R R , M . E . (1966). F o u n d a t i o n s of t h e o r e t i c a l soil m e c h a n i c s . M c G r a w - H i l l ,
N e w Y o r k , N Y . 3 8 1 pp.
H E T T I A R A T C H I , D . R . P . (1969). The c a l c u l a t i o n of passive e a r t h p r e s s u r e . P h . D .
T h e s i s , U n i v . of N e w c a s t l e - u p o n - T y n e .
H E T T I A R A T C H I , D . R . P . and R E E C E , A . R . (1974). T h e c a l c u l a t i o n of passive
soil r e s i s t a n c e . G e o t e c h n i q u e , V o l . 24, N o . 3, pp. 2 8 9 - 3 1 0 .
H O E K S T R A , P., C H A M B E R L A I N , E . and F R A T E , T . (1965). F r o s t h e a v i n g
p r e s s u r e s . C o l d R e g i o n s R e s . E n g . L a b . R e s . R e p o r t 176, H a n n o v e r , N H .
H O O G H O U D T , S.B. (1940). B i j d r a g e tot de kennis van enige natuurkundige
g r o o t h e d e n v a n de g r o n d . V e r s l a g e n v a n L a n d b o u w k u n d i g e O n d e r z o e k i n g e n ,
V o l . 4 6 , N o . 7, pp. 5 1 5 - 7 0 7 .
H O V A N E S I A N , J . D . (1958). D e v e l o p m e n t and use of a v o l u m e t r i c transducer
f o r s t u d y i n g e f f e c t s of soil p a r a m e t e r s on c o m p a c t i o n . P h . D . thesis,
Michingan State Univ., East Lansing, M I .
H U D S O N , N . (1981). Soil c o n s e r v a t i o n . C o r n e l l U n i v . P r e s s , I t h i c a , N Y , 324 pp.
H V O R S L E V , M . J . (1937). U b e r die F e s t i g k e i t s e i g e n s c h a f t e n g e s t o r t e r bindinger
B o d e n . D a n m a r k s N a t u r v i d e n k a b e l i g e Samfund, Ingeniorvidenskabelige
S k r i f t e r , S e r i e s A , N o . 4 5 , C o p e n h a g e n , 159 pp.
I L R I (1964). C o d e of p r a c t i c e for the d e s i g n of open w a t e r c o u r s e s and
a n c i l l a r y s t r u c t u r e s . I n t . I n s t , f o r L a n d R e c l a m a t i o n and I m p r o v e m e n t ,
W a g e n i n g e n , T h e N e t h e r l a n d s , B u l l . N o . 7.
I R E L A N D , H . O . (1957). P u l l i n g tests on piles in s a n d . P r o c . 4 t h I n t . C o n f . on
Soil M e c h . and F o u n d . E n g . , L o n d o n , V o l . I I , pp. 4 3 - 4 5 .
J A K Y , J . ( 1 9 4 8 ) . P r e s s u r e s in s i l o s . P r o c . 2nd I n t . C o n f . Soil M e c h . ,
R o t t e r d a m , V o l . 1, pp. 1 0 3 - 1 0 8 .
J A N I G A , P.V. (1970). S o m e c o n s i d e r a t i o n s in i n - s i t u f r o s t h e a v i n g . M . S c .
Thesis, McGill University, Montreal, Q C .
J A N S S E N , H . A . (1895). V e r s u c h e uber G e t r e i d e d r u c k in S i l o z e l l e n . Z . V e r e i n e s
D e u t s c h e r I n g e n i e u r e , V o l . 39, pp. 1 0 4 5 - 1 0 4 9 .
K E T C H U M , M . S . (1919). D e s i g n of bins, w a l l s and g r a i n e l e v a t o r s . M c G r a w
Hill, New York, N Y .
K O E R N E R , R . M . (1986). D e s i g n i n g w i t h g e o s y n t h e t i c s . P r e n t i c e - H a l l , E n g l e w o o d
C l i f f s , N J . 4 2 4 pp.
262 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 2. D E S I G N R E Q U I R E M E N T S F O R E X T E N D E D S I L O RING
FOUNDATIONS
Φ δ° a 0
Φ°δ° α ° δ°
0
Κ Κ
ah K
av ah av
φο
α° ah av
36 32 90 0 . 2 0 2 0 . 1 2 6 34 30 90 0 . 2 2 1 0 . 1 2 8 32 28 90 0 . 2 4 3 0 . 1 2 9
85 .222 .167 85 .242 .169 85 .264 .171
80 .240 .216 80 .261 .219 80 .283 .221
75 .256 .275 75 .278 .278 75 .301 .280
70 .269 .345 70 .292 .348 70 .315 .350
65 .279 .430 65 .302 .431 65 .326 .433
60 .283 .533 60 .307 .533 60 .333 .533
30 26 90 0 . 2 6 6 0 . 1 3 0 28 25 90 0 . 2 8 9 0 . 1 3 5 26 23 90 0 . 3 1 6 0 . 1 3 4
85 .287 .173 85 .310 .179 85 .337 .179
80 .307 .223 80 .330 .231 80 .357 .232
75 .325 .283 75 .347 .291 75 .375 .293
70 .340 .352 70 .362 .362 70 .390 .364
65 .352 .435 65 .374 .445 65 .403 .447
60 .360 .533 60 .381 .545 60 .411 .545
24 2 1 90 0 . 3 4 5 0 . 1 3 2 22 19 90 0 . 3 7 7 0 . 1 3 0 20 17 90 0 . 4 1 2 0 . 1 2 6
85 .367 .179 85 .399 .178 85 .434 .175
80 .387 .233 80 .419 .232 80 .454 .231
75 .405 .294 75 .437 .295 75 .472 .295
70 .421 .366 70 .453 .367 70 .488 .368
65 .433 .449 65 .466 .450 65 .501 .451
60 .442 .546 60 .476 .547 60 .511 .548
36 24 90 0 . 2 1 5 0 . 0 9 6 34 23 90 0 . 2 3 4 0 . 0 9 9 32 2 1 90 0 . 2 5 7 0 . 0 9 9
85 .238 .132 85 .257 .137 85 .281 .137
80 .260 .175 80 .280 .182 80 .304 .183
75 .281 .228 75 .301 .235 75 .325 .236
70 .300 .290 70 .320 .298 70 .345 .300
65 .316 .364 65 .336 .373 65 .362 .375
60 .330 .454 60 .349 .464 60 .376 .465
272 AGRICULTURAL ENGINEERING SOIL MECHANICS
</>° δ ° « ° K g h K a v φο δ ο α οφ ο δ ο «ο K g h K j
30 20 90 0 . 2 7 9 0 . 1 0 2 28 19 90 0 . 3 0 4 0 . 1 0 5 26 17 90 0 . 3 3 2 0 . 1 0 2
85 .304 .142 85 .328 .146 85 .357 .144
80 .326 .188 80 .350 .194 80 .380 .194
75 .348 .244 75 .372 .251 75 .401 .251
70 .367 .308 70 .391 .316 70 .421 .317
65 .384 .384 65 .408 .394 65 .438 .395
60 .398 .475 60 .422 .485 60 .453 .486
24 16 90 0 . 3 6 1 0 . 1 0 3 22 15 90 0 . 3 9 1 0 . 1 0 5 20 13 90 0 . 4 2 7 0 . 0 9 9
85 .385 .148 85 .415 .151 85 .451 .147
80 .407 .199 80 .437 .204 80 .473 .201
75 .429 .258 75 .458 .264 75 .494 .263
70 .448 .325 70 .476 .333 70 .513 .333
65 .465 .404 65 .493 .413 65 .529 .414
60 .479 .496 60 .506 .506 60 .544 .507
APPENDIX 4 273
APPENDIX 4. C O N S O L I D A T I O N T I M E FACTORS
274 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 5 275
APPENDIX 5. P A S S I V E S O I L C U T T I N G FACTORS
Ρ = (yd Ny
2
+ cdN c + q dN q + c dN
a c a ) w
276 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
0 10 20 30 40 50 60 70 80 90
Rake angle a0
0 10 20 30 40 50 60 70 80 90
Rake angle a°
APPENDIX 5 277
278 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 5 279
280 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
300;
0 10 20 30 40 50 60 70 80 90
Rake angle a°
APPENDIX 5 281
282 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 5 283
284 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 5 285
300
200
0 10 20 30 40 50 60 70 80 90
Rake angle ct°
0 10 20 30 40 50 60 70 80 90
Rake angle, a0
286 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 5 287
3001 1 1 1 1 1 1 1 r
200
0 10 20 30 40 50 60 70 80 90
Rake angle a 0
0 10 20 30 40 50 60 70 80 90
Rake angle, a 0
288 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
APPENDIX 6. S E L E C T E D V A L U E S O F S O I L M E C H A N I C A L PROPERTIES
( f r o m B e k k e r , 1969).
APPENDIX 6 289
Soil m.c. kc k0 η c
% kPa/m n _ 1
kPa/m n
kPa
Dry sand 3
0 1 1527 1.1 1.0 28
Sandy l o a m 3
15 5 1514 0.7 1.7 29
22 3 43 0.2 1.4 38
L e a n clay* 3
22 16 1723 0.20 68.9 20
24 10 1030 0.17 48.2 18
26 7 580 0.17 34.5 15
28 3 372 0.16 27.6 12
30 3 248 0.16 17.2 11
32 2 120 0.15 13.8 11
3
Land Locomotion Laboratory, Detroit M I .
b
Waterways Experimentation Station, Vicksburg MS.
290 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
AUTHOR INDEX
R o w l e s , W., 222,264
R y c r o f t , D.W., 115,117,132,264
S c h m i d t , R . H . , 222,260
S c h w a b , G.O., 117,132,176,186,264
S c o b e y , F . C . , 131,260
S c o t t , R.F., 84,264
Sill, R . C . , 252,264
S i n g h , D., 163,264
S k a p s k i , A . S . , 252,264
S k e m p t o n , A . W . , 72,264
S m e d e m a , L.K., 115,117,132,264
S m i t h , D . D . , 174,175,264,265
Soehne, W., 224,264
S o k o l o v s k i , V.V., 28,64-66,69,
146,159,264
Soiltest Inc., 52,264
Spangler, M . G . , 165-168,264
Spoor, G . , 197,212-216,260,264
T a y l o r , D.W., 7 8 „ 1 2 0 - 1 2 4 , 2 2 2 , 2 6 4
T e n g , W . C . , 15,264
T e r z a g h i , K., 29,61,64,66,70,78,
79,84,94,95,100,102,106,144,
149,254,264,265
T u r i t z i n , A . M . , 164,265
T u r n b u l l , J.E., 87,91,92,265,268
V a n d e n B e r g , G . E . , 224,265
van S c h i l f g a a r d e , J . , 115,265
V o m i c i l , J.A., 230,265
W a g n e r , A . A . , 7,265
Wall, 173,265
W a r k e n t i n , B.P., 245,248,265
Wendlaar, F.E., 177,265
W h i t m a n , R.V., 31,79,221,265
W i l l i a m s , J.R., 176,265
W i s c h m e i e r , W . H . , 172,174,264,265
W i s m e r , R . D . , 202,262
Y o n g , R . N . , 244,245,248,253,265
Z e l e n i n , A . N . , 207,265
292 AGRICULTURAL E N G I N E E R I N G SOIL MECHANICS
SUBJECT INDEX