Geotech Lecture Notes & Tut Combine
Geotech Lecture Notes & Tut Combine
1
https://www.polyu.edu.hk/cee/people/academic-staff/dr-zhen-yu-yin/
2
Self introduction – education & experiences
Education NTNU
97, BEng in ZJU NGI
03, Msc in ECN (Fr) HUT
06, PhD in ECN (Fr) Umass TUChalmers
(~9 months) Glasgow
SJTU, TongJi
Research fellow (~2 years)
(~4 years)
06-08, StrathU PLAXIS
08-09, Umass ZJU
09-10, ECN U-Stuttgart HK
Engineering
applications
Deformation & strength
degradation
Microstructure &
micromechanics
Contact information
Dr Zhenyu YIN
Department of Civil and Environmental Engineering
The Hong Kong Polytechnic University
Tel: 3400 8470; Office: ZS915
Email: zhenyu.yin@polyu.edu.hk
4
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
Introduction of
Geotechnical Engineering
by Dr. Zhenyu YIN
Office: ZS915, Tel: 3400-8470
Email: zhenyu.yin@polyu.edu.hk
5
❑ Underground geotechnics
Hong Kong
Underground space
(infrastructures)
❑ Challenges
Shanghai
Salisbury Road, HK
Leakage of tunnel Fukuoka Nov. 8 2016 June 3 2007
❖ Long-term subsidence ❖ Internal erosion induced damage
6
❑ Underground geotechnics
Example: tunneling induced problems
Subsidence
Face stability
7
❑ Underground geotechnics
8
❑ Foundation engineering
Tower of Pisa
❖ Long-term settlement Long-term
(time effects of clayey soils) settlement
(Inclination)
Ma On Shan, HK 2002
9
❑ Foundation engineering
10
❑ Foundation engineering
Example: long-term settlement of Pisa tower
11
❑ Transportation geotechnics
Settlement of
embankment
Instability of
embankment
Problem of ballast
Differential settlement
and embankment
(Train in and out of the station)
13
❑ Transportation geotechnics
Example: development of high-speed railway in the world
14
❑ Transportation geotechnics
The differential
settlement has to
be strictly
controlled!
15
❑ Transportation geotechnics
Example: Mud pumping in railway engineering
16
❑ Transportation geotechnics
Breakage mechanism?
Mode?
Crushed particles?
17
❑ Environmental geotechnics
Anti-seepage
18
❑ Environmental geotechnics
19
❑ Environmental geotechnics
20
❑ Energy geotechnics
21
❑ Energy geotechnics
22
❑ Geotechnical hazards
Unstable?
Displacement
Ya An debris flow,
Sichuan 2013 Stable?
Time
24
❑ Geotechnical hazards
Example:
massive landslides
after a rainstorm
25
❑ Geotechnical hazards
26
❑ Offshore / submarine geotechnics
❑ Challenges
❖ Long-term instability
(clay)
❖ Grain breakage induced damage
(carbonate sand)
❖ Erosion induced damage
(water flow coupling)
27
❑ Offshore / submarine geotechnics
28
❑ Offshore / submarine geotechnics
29
❑ Offshore / submarine geotechnics
Failure of anchor
Pumping
Water
Suction pressure
Uplift Seepage
of clay liquefaction
Torpedo anchor penetration
Soft clay Sand
30
❑ Offshore / submarine geotechnics
32
❑ Earth works (or hydraulic works)
West Kowloon Cultural Hong Kong airport
District, 40 hectares (reclamation)
(9.46h-PolyU)
33
❑ Earth works (or hydraulic works)
Lantau tomorrow (https://www.lantau.gov.hk/en/sustainable-lantau-office/index.html)
Ma Liu Shui
Lung Kwu Tan
Sunny Bay
(1700 hectares)
East Lantau Metropolis
(1700 hectares)
Siu Ho Wan
34
❑ Earth works (or hydraulic works)
Fill height
Rupture of dam (rockfill)
Kansai airport, Japan
Creeping
Data of
Kansai airport Sinkhole on
(Settlement up to 9 m) Failure of dike top of dam
35
❑ Earth works (or hydraulic works)
36
❑ Earth works (or hydraulic works)
37
Geotechnical Design
Topic 1: Site investigation
Topic 2: Slope Stability
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
Textbook:
❑ B. M. Das. (2011). Principles of Foundation
Engineering, 7th Edition.
❑ B. M. Das & K. Sobhan (2016). Principles of
Geotechnical Engineering, Ninth Edition.
❑ Craig, R.F. (2004). Soil Mechanics, 7th
edition (6thor 5th edition), Spon Press,
London and New York. 38
Break
39
Additional books/references for this course
40
Additional books/references for this course
41
Additional references for this course
Code of Practice for Foundation, Buildings Department, 2017.
Hong Kong Foundation Handbook, HKIE, 2015/2017
Guide to Retaining Wall Design, GEO Guide 1, GEO
Guide to soil and rock description, GEO Guide 3, GEO
Foundation Design and Construction, GEO, 2006.
Prescriptive measures for man-made slopes and retaining walls, GEO 2009
Guide to Site investigation, GEO Guide 2, GEO
Review of design methods for excavation, GCO, 1990
Hong Kong Geology Guidebook, GEO, 2007
Engineering Geology Practice in HK, GEO, 2007
Practice note from Buildings Department
Advanced Soil Mechanics, BM Das, 4ed, CRC Press, 2014.
Shallow Foundations Bearing Capacity and Settlement, 3rd edition, B.M. Das, CRC Press, 2017
Soil Slope stability analysis and stabilization – new methods and insights, 2nd edition, Cheng
Y.M., Spon Press 2014
Geotechnical Correlations for Soils and Rocks, Verbrugge, Wiley, 2018.
In Situ Tests in Geotechnical Engineering, Monnet, Wiley, 2015.
42
The design of geotechnical structures
- generally requires a knowledge of such factors as
(a) the load that will be transmitted by the superstructure to the foundation
system, (e.g. stress in soil)
(b) the requirements of the local engineering code, (e.g. factor of safety)
(c) the behavior and stress-related deformability of soils that will support the
foundation system, (e.g. mechanical properties of soils)
(d) the geological conditions of the soil under consideration. (e.g. type of soils
and depth)
For example:
when determining which
ICC (484 m) IFC (415 m)
foundation is the most
economical, the engineer
must consider the
superstructure load, the
subsoil conditions, and the
desired tolerable settlement.
43
Importance of Geotechnical Engineering in HK
Special features of Hong Kong:
(a) Scarce of land (not exactly true), high population
(b) Reclamation and construction adjacent to slopes
(c) Deep basement for various commercial uses
(d) Large amount of materials are involved
Outcome:
(a) Lots of slope failures – average 300 per year
(b) Large diameter deep piles
(c) Deep excavation (underground space)
(d) High ground tables, loose soil creates various problems
(e) Loose backfill create settlement and failure problems
44
Nature of Geotechnical Engineering
45
Analysis of Geotechnical Engineering
(a) Rule of thumb: rely mainly on past experience, useful, but not
applicable to new structure or situations
(b) Statistics: Interpretation of field test results, settlement, soil
movement etc., useful but may not be representative
(c) Analytical solutions: Only for some simple problems. Use of design
graphs and tables are also common to replace the use of long
equations. Many classical methods rely on assumptions which may
be approximations.
(d) Computational method: Rely on use of computer programs with
less assumptions, needs knowledge and judgment in using these
programs. otherwise, just Rubbish in -> rubbish out.
47
Break
49
TOPICS & SYLLABUS:
I. Geotechnical Properties of Soil
Topic 1: Site investigation
II. Natural Soil Deposits and Subsoil
Topic 2: Slope Stability Exploration
50
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
Site investigation
Part I- Geotechnical Properties of Soil
51
Main content:
Recall of Soil
Mechanics
I. Physical properties
❖ Size Limits for Soils and Grain-Size Distribution
❖ Weight–Volume Relationships
❖ Relative Density, Atterberg Limits and Activity
❖ Soil Classification Systems
52
1 Size Limits for Soils and Grain Size Distribution
Ballast
Sand
East Rail line of HK Golden Beach in HK 53
1 Size Limits for Soils and Grain Size Distribution
54
1 Size Limits for Soils and Grain Size Distribution
Sieve Analysis
A sieve analysis is conducted by taking a measured amount of dry, well-pulverized
soil and passing it through a stack of progressively finer sieves with a pan at the
bottom. The amount of soil retained on each sieve is measured, and the cumulative
percentage of soil passing through each is determined. This percentage is generally
referred to as percent finer.
D60
Uniformity coefficient Cu =
D10
Coefficient of gradation
(or coefficient of curvature)
( D30 )
2
Cc =
( D60 )( D10 )
56
1 Size Limits for Soils and Grain Size Distribution
Example: Completely Decomposed Granite (CDG), a typical granular soil in Hong Kong
Uniformity coefficient Cu? Coefficient of gradation Cc?
100
80
Percent finer (%)
60
40
D60 = 0.15
20 D30 = 0.0129
D10 = 0.0015
0
0.001 0.01 0.1 1 10 100
Particle size (mm)
( D30 ) = 0.74
2
D
Cu = 60 = 100 Cc =
D10 ( D60 )( D10 )
57
1 Size Limits for Soils and Grain Size Distribution
Hydrometer Analysis
In case of fine grained soil, sieve analysis test does not
give reliable test result. This because a fine grained soil
consist of different sizes of particles starting from 0.075
mm to 0.0002 mm. and it is not practicable to design
sieve having so smaller screen size. Also, there is a George Gabriel Stokes
chance of lost of sample during sieving. Therefore, (1819-1903)
hydrometer analysis is done for grain size analysis of
fine grained soils.
Hydrometer analysis is based on Stokes’ law. According 1 ( Gs − 1) w 2
v= D
to this law, the velocity at which grains settles out of 18
suspension, all other factors being equal, is dependent D = diameter of the soil particle
upon the shape, weight and size of the grain. Gs = specific gravity of soil solids
= viscosity of water
In case of soil, it is assumed that the soil particles are w = unit weight of water
spherical and have the same specific gravity. Therefore, v = Terminal velocity of particle
we can say that in a soil water suspension the coarser
particles will settle more quickly than the finer ones.
58
1 Size Limits for Soils and Grain Size Distribution
Sediment distance of the particle
From Stokes’ law to particle size L within the time t
1 ( Gs − 1) w 2 18 18 L
v= D D= v=
18 ( Gs − 1) w ( Gs − 1) w t
90
80
Deep marine
No.2
soils
Percentage finer (%)
70 No.4
20
Top marine soils From 39.5 to 40.5
From 47.2 to 48.2
10
GSD of HKMC 0
From 31.4 to 32.4
Vv Special case
e=
Vs
60
2 Weight–Volume Relationships
61
3 Relative Density, Atterberg Limits and Activity
In granular soils, the degree of compaction in the field can be measured
according to the relative density, defined as
emax − e d − d (min) d (max)
Dr ( % ) = 100 or 100
emax − emin d (max) − d (min) d
emax = void ratio of the soil in the loosest state
emin = void ratio in the densest state
e = in-situ void ratio Attention: physically
loose or dense, not
mechanically!
Denseness of a Granular Soil
62
3 Relative Density, Atterberg Limits and Activity
When a clayey soil is mixed with an excessive amount of water, it may flow like
a semiliquid. If the soil is gradually dried, it will behave like a plastic, semisolid,
or solid material, depending on its moisture content.
The moisture content, in percent, at which the soil changes from a liquid to a
plastic state is defined as the liquid limit (LL or wL). Similarly, the moisture
content, in percent, at which the soil changes from a plastic to a semisolid
state and from a semisolid to a solid state are defined as the plastic limit (PL
or wP) and the shrinkage limit (SL or wS), respectively.
These limits are referred to as Atterberg limits
PI = LL − PL
or PI = wL − wP
(wS) (wP) (wL)
63
3 Relative Density, Atterberg Limits and Activity
Because the plasticity of soil is caused by the adsorbed water that surrounds
the clay particles, we can expect that the type of clay minerals and their
proportional amounts in a soil will affect the liquid and plastic limits.
Skempton (1953) observed that the plasticity index of a soil increases linearly
with the percentage of clay-size fraction (% finer than 2 mm by weight)
present. The correlations of PI with the clay-size fractions for different clays
plot separate lines. This difference is due to the diverse plasticity
characteristics of the various types of clay minerals.
On the basis of these results, Skempton defined a quantity called activity,
which is the slope of the line correlating PI and % finer than 2 m, expressed as
PI
A=
% of clay-size fraction, by weight
65
4 Soil Classification Systems
Soil classification systems divide soils into groups and subgroups based on
common engineering properties such as the grain-size distribution, liquid limit,
and plastic limit. The two major classification systems presently in use are:
(1) the American Association of State Highway and Transportation Officials
(AASHTO) System, mainly for the classification of highway subgrades;
(2) the Unified Soil Classification System (also ASTM), for all geotechnical work
66
4 Soil Classification Systems
The Unified Soil Classification System was originally proposed by A.
Casagrande in 1942 and was later revised and adopted by the United States
Bureau of Reclamation and the U.S. Army Corps of Engineers.
Arthur Casagrande
(1902-1981)
67
Break
68
5 Hydraulic Conductivity of Soil and Steady-State Seepage
The void spaces, or pores, between soil grains allow
water to flow through them. In soil mechanics and
foundation engineering, you must know how much
water is flowing through a soil per unit time. This
knowledge is required to design earth dams,
determine the quantity of seepage under hydraulic Henry Darcy
structures, and dewater foundations before and (1803-1858)
during their construction.
h
i=
v = ki L
−4 e
3
0.6 2.32 w
k = 3.5 10 Cu D10 for sand
1 + e Type of soil k (cm/sec)
en
k = C 1 + e for clay
69
5 Hydraulic Conductivity of Soil and Steady-State Seepage
For most cases of seepage under hydraulic structures, the flow path changes
direction and is not uniform over the entire area. In such cases, one of the ways
of determining the rate of seepage is by a graphical construction referred to as
the flow net, a concept based on Laplace’s theory of continuity.
2h 2h 2h
kx 2 + k y 2 + kz 2 = 0 ❖ A flow line is a path that a
x y z water particle would
follow in traveling from
the upstream side to the
downstream side.
❖ An equipotential line is a
line along which water, in
piezometers, would rise
x to the same elevation.
❖ flow lines perpendicular
to equipotential lines
❖ A flow net is a
combination of numerous
equipotential lines and
flow lines.
z
70
5 Hydraulic Conductivity of Soil and Steady-State Seepage
A number of flow lines and equipotential lines are drawn by trial and error so
that all the flow elements in the net have the same length-to-width ratio (L/B).
All flow lines must intersect all equipotential lines at right angles (perpendicular).
❖ The space between two consecutive flow lines is defined as a flow channel.
❖ The space between two consecutive equipotential lines is called a drop.
Seepage: q = kh N f n
max
Nd
Nf = number of flow channels
Nd = number of drops
n = B/L
Flow line
(assumed isotropic)
Equipotential line B for flow channel, L for drop
71
6 Effective Stress
The total stress at a given point in a soil mass can be expressed as
=+u
Pore water pressure
Total stress Effective stress
(used to calculate settlement
Example 1 and soil strength)
Example 2
= w h1 + sat h2 u = w ( h1 + h2 + h )
= ( h1 + sat h2 ) − w h2
h
= h1 + h2 = − u = − w h2 = ( − i w ) h2
h2 72
7 Consolidation
In the field, when the stress on a saturated clay layer is increased, for example
by the construction of a foundation, the pore water pressure in the clay will
increase. Because the hydraulic conductivity of clays is very small, some time
will be required for the excess pore water pressure to dissipate and the
increase in stress to be transferred to the soil skeleton.
This gradual increase in the effective stress in the clay layer will cause
settlement over a period of time and is referred to as consolidation.
u = w hi
73
7 Consolidation
Oedometer test:
❖ Specimens are placed inside a ring, with
one porous stone at the top and one at
the bottom of the specimen.
❖ A load on the specimen is then applied
so that the total vertical stress is equal
to Settlement readings for the
Preconsolidation Pressure specimen are taken periodically for 24
hours.
❖ The load on the specimen is doubled
and more settlement readings are
taken.
Compression Index
e
Swelling Index Cc or Cs =
log
74
7 Consolidation
log ( z )
(Compression index)
e c z0
z =
1 + e0 z
Cs (1 + e0 ) (ii)
(iii) From overconsolidated to normally consolidated
with z0 c & z0 + z c : Cc (1 + e0 )
C C + z
Sc ( p ) = s log c + c log z 0 H0 e
1 + e0
z0 1 + e0
c z =
1 + e0 77
7 Consolidation
78
7 Consolidation
79
7 Consolidation
ue (0, t ) = 0
❖ Boundary conditions:
ue (2d , t ) = 0, t 0
Depth i
cv t
❖ We can define an operator b: b=
( z )
2
❖ Exact solution for excess pore water pressure and degree of consolidation:
m =
2u i Mz
ue =
m=0 M
sin
d
exp(− M Tv )
2
Time factor:
M= (2m + 1)
2 ct
m =
Tv = v2 (unitless)
2
U =1−
m=0 M 2
exp( − M 2
Tv ) d
ue S
(Average degree of consolidation U = 1 − = c (t) )
ui Sc (max)
❖ Approximate solutions:
2
U U 0.6
Tv = 4
−0.933log(1 − U ) − 0.085 U 0.6
84
8 Shear Strength
The shear strength of a soil, defined in terms of effective stress, is
(Mohr–Coulomb failure criterion)
N
n =
A f = c + n tan
R
=
A
Charles-Augustin
de Coulomb
(1736-1806)
85
8 Shear Strength
(a) Direct shear n
,
(, )
86
8 Shear Strength
z (z)
(d) Biaxial shear
h
=0 (~h)
r
r
Hollow cylinder pv T
(f) Others
torsional shear
True triaxial
1
z
1 a
po pi z q
q
3 3
2 r=2
87
8 Shear Strength
Triaxial compression tests can be conducted on
sands and clays to measure shear strength
(CD)
(UU) 88
8 Shear Strength
Triaxial compression tests can be conducted on
sands and clays to measure shear strength
89
9 Unconfined Compression Test and Sensitivity
For many naturally deposited clay soils, the unconfined compression strength is
much less when the soils are tested after remolding without any change in the
moisture content. This property of clay soil is called sensitivity. The degree of
sensitivity is the ratio of the unconfined compression strength in an
undisturbed state to that in a remolded state: qu ( undisturbed )
St =
qu ( remolded )
90
9 Unconfined Compression Test and Sensitivity
The sensitivity ratio of most clays ranges from about 1 to 8; however, highly
flocculent marine clay deposits may have sensitivity ratios ranging from about
10 to 80. Some clays turn to viscous liquids upon remolding, and these clays are
referred to as “quick” clays. The loss of strength of clay soils from remolding is
caused primarily by the destruction of the clay particle structure that was
developed during the original process of sedimentation.
Intact and remolded
Leda clay from Ottawa
91
Summary of Part I: Geotechnical Properties of Soil
I. Physical properties
❖ Size Limits for Soils and Grain-Size Distribution
❖ Weight–Volume Relationships
❖ Effective Stress
❖ Consolidation
❖ Shear Strength
92
Summary of HKMC
94
Summary of HKMC (after Yeung & So 2001)
6sin 3M
M= = sin −1
3 − sin 6+M
M
95
Summary of HKMC (after Yeung & So 2001)
96
Summary of HKMC (after Yeung & So 2001)
97
Exercise 1
The moist weight of 3 x 10-3 m3 of soil is 50 N. If the moisture content is 20%, and the
specific gravity of soil solids is 2.7, find the following:
a. Moist unit weight (kN/m3)
b. Dry unit weight (kN/m3)
c. Void ratio
d. Porosity
e. Degree of saturation (%)
f. Volume occupied by water (m3)
98
Exercise 2
The maximum and minimum dry unit weights of a sand are 17.1 and 14.2 respectively. The
sand in the field has a relative density of 70% with a moisture content of 12%. Determine
the moist unit weight of the sand in the field.
100
Exercise 3
A laboratory consolidation test on a normally consolidated clay showed the following
results:
The specimen tested was 25.4 mm in thickness and drained on both sides. The time required
for the specimen to reach 50% consolidation was 5 min.
A similar clay layer in the field 2.8 m thick and drained on both sides, is subjected to a similar
increase in average effective pressure (i.e., ’0 = 150 kN/m2 and ’0 + ’ = 220 kN/m2).
Please determine:
a. the expected maximum primary consolidation settlement in the field.
b. the length of time required for the total settlement in the field to reach 40 mm. (Assume
a uniform initial increase in excess pore water pressure with depth.)
102
References relating to HKMC
[1] S.-F. Zou, X.-Y. Xie, J.-Z. Li, Z.-J. Wang, H.-Y. Wang, Rheological characteristics and one-dimensional isotache modelling of marine soft clays,
Marine Georesources & Geotechnology, 37 (2019) 660-670.
[2] J. Yin, W. Feng, Validation of a new simplified hypothesis B method for calculating consolidation settlement of clayey soils exhibiting creep,
Geotechnical Engineering, (2018).
[3] W.-Q. Feng, B. Lalit, Z.-Y. Yin, J.-H. Yin, Long-term non-linear creep and swelling behavior of Hong Kong marine deposits in oedometer condition,
Computers and Geotechnics, 84 (2017) 1-15.
[4] F. Tong, J.-H. Yin, Experimental and constitutive modeling of relaxation behaviors of three clayey soils, Journal of Geotechnical and Geoenvironmental
Engineering, 139 (2013) 1973-1981.
[5] F. Tong, J.-H. Yin, H. Pei, Experimental study on complete consolidation behavior of Hong Kong marine deposits, Marine Georesources &
Geotechnology, 30 (2012) 291-304.
[6] F. Tong, J.-H. Yin, Nonlinear creep and swelling behavior of bentonite mixed with different sand contents under oedometric condition, Marine
Georesources & Geotechnology, 29 (2011) 346-363.
[7] Z. Fang, J.-H. Yin, Physical modelling of consolidation of Hong Kong marine clay with prefabricated vertical drains, Canadian Geotechnical Journal, 43
(2006) 638-652.
[8] J.-H. Yin, C.-M. Cheng, Comparison of strain-rate dependent stress-strain behavior from K o-consolidated compression and extension tests on natural
Hong Kong marine deposits, Mar. Georesour. Geotechnol. , 24 (2006) 119-147.
[9] C.-M. Cheng, J.-H. Yin, Strain-Rate Dependent Stress--Strain Behavior of Undisturbed Hong Kong Marine Deposits under Oedometric and Triaxial
Stress States, Mar. Georesour. Geotechnol. , 23 (2005) 61-92.
[10] J.-H. Yin, Stress-strain strength characteristics of a marine soil with different clay contents, Geotech. Test. J. , 25 (2002) 459-462.
[11] J.-H. Yin, J.-G. Zhu, J. Graham, A new elastic viscoplastic model for time-dependent behaviour of normally and overconsolidated clays: theory and
verification, Canadian Geotechnical Journal, 39 (2002) 157-173.
[12] J. Zhu, J. Yin, Drained creep behaviour of soft Hong Kong marine deposits, Géotechnique, 51 (2001) 471-474.
[13] J.-H. Yin, Time and strain-rate effects on the stress-strain behaviour of soft Hong Kong marine deposits, in: Proceedings of the international
conference on soil mechanics and geotechnical engineering, aa balkema publishers, 2001, pp. 333-336.
[14] P. Lee, Y. Lam, W. Wong, Evaluation of a rational multistage triaxial test method on the strength parameters of Hong Kong marine clay, in: Soft Soil
Engineering, Routledge, 2017, pp. 683-688.
[15] J.-G. Zhu, J.-H. Yin, Strain-rate-dependent stress-strain behavior of overconsolidated Hong Kong marine clay, Canadian Geotechnical Journal, 37
(2000) 1272-1282.
[16] J.-H. Yin, J.-G. Zhu, Measured and predicted time-dependent stress-strain behaviour of Hong Kong marine deposits, Canadian Geotechnical Journal,
36 (1999) 760-766.
[17] J.-G. Zhu, J.-H. Yin, S.-T. Luk, Time-dependent stress-strain behavior of soft Hong Kong marine deposits, Geotech. Test. J. , 22 (1999) 118-126.
[18] J.-H. Yin, Properties and behaviour of Hong Kong marine deposits with different clay contents, Canadian Geotechnical Journal, 36 (1999) 1085-1095.
[19] X. LEE PKK, A laboratory investigation of the permeability characteristics of Hong Kong marine clay, in: Recent Advances in Soft Soil Engineering:
Proceedings of the International Conference on Soft Soil Engrg (ICSSE). Guangzhou: Science Press, 1993, pp. 2-127.
105
References relating to CDG
[20] R.P. Chapuis, M. Aubertin, Discussion of "Influence of relative compaction on the hydraulic conductivity of completely decomposed granite in Hong Kong", Canadian Geotechnical
Journal, 47 (2010) 704-707.
[21] H. Chen, C.F. Lee, K.T. Law, Causative mechanisms of rainfall-induced fill slope failures, Journal of Geotechnical and Geoenvironmental Engineering, 130 (2004) 593-602.
[22] H. Chen, C.F. Lee, J.M. Shen, Mechanisms of rainfall-induced landslides in Hong Kong, 2000.
[23] L.M. Chu, J.H. Yin, A laboratory device to test the pull-out behavior of soil nails, Geotech. Test. J. , 28 (2005) 499-513.
[24] S.M. Dasaka, L.M. Zhang, Spatial variability of in situ weathered soil, Géotechnique, 62 (2012) 375-384.
[25] J.K.M. Gan, D.G. Fredlund, Shear strength characteristics of two saprolitic soils, Canadian Geotechnical Journal, 33 (1996) 595-609.
[26] M.A. Hossain, J.H. Yin, Behavior of a Compacted Completely Decomposed Granite Soil from Suction Controlled Direct Shear Tests, Journal of Geotechnical and Geoenvironmental
Engineering, 136 (2010) 189-198.
[27] M. Kumruzzaman, J.H. Yin, Influence of the intermediate principal stress on the stress-strain-strength behaviour of a completely decomposed granite soil, Géotechnique, 62 (2012) 275-
280.
[28] G.L.M. Leung, W.G. Wong, Fundamental findings of dynamic behaviour of pavement subgrade material in hong kong, 2010.
[29] A.G. Li, L.G. Tham, Z.Q. Yue, C.F. Lee, K.T. Law, Comparison of field and laboratory soil-water characteristic curves, Journal of Geotechnical and Geoenvironmental Engineering, 131
(2005) 1176-1180.
[30] P. Liu, X.W. Zhou, Y.B. He, Bond Yield Characteristics of Undisturbed Completely Decomposed Granite, Advances in Materials Science and Engineering, (2015).
[31] K.K. Lu, J.H. Yin, S.C. Lo, Modeling Small-Strain Behavior of Hong Kong CDG and Its Application to Finite-Element Study of Basement-Raft Footing, Int. J. Geomech. , 18 (2018).
[32] B.N. Madhusudhan, B.A. Baudet, Influence of reconstitution method on the behaviour of completely decomposed granite, Géotechnique, 64 (2014) 540-550.
[33] A.W. Malone, C.W.W. Ng, J.W. Pappin, I. Publications Comm Of Xiv, Panel discussion: The prediction and the control of displacements around deep excavations in completely
decomposed granite, 1999.
[34] C.W.W. Ng, X.S. Li, P.A. Van Laak, D.Y.J. Hou, Centrifuge modeling of loose fill embankment subjected to uni-axial and bi-axial earthquakes, Soil Dynamics and Earthquake
Engineering, 24 (2004) 305-318.
[35] K.T.W. Ng, I.M.C. Lo, Fines migration from soil daily covers in Hong Kong landfills, Waste Manage. (Oxford) 30 (2010) 2047-2057.
[36] Y.J. Shang, S.J. Wang, Z.Q. Yue, Variation of shearing; strength parameters of undisturbed CDG in Hong Kong, 2003.
[37] S.Y. Wang, D. Chan, K.C. Lam, Experimental study of the effect of fines content on dynamic compaction grouting in completely decomposed granite of Hong Kong, Construction and
Building Materials, 23 (2009) 1249-1264.
[38] S.Y. Wang, D. Chan, K.C. Lam, S.K. Au, L.G. Tham, Soil improvement using compaction grouting - a laboratory investigation on the confining pressure and injection rate in completely
decomposed granite, 2007.
[39] S.Y. Wang, D.H. Chan, K.C. Lam, S.K.A. Au, Numerical and experimental studies of pressure-controlled cavity expansion in completely decomposed granite soils of Hong Kong,
Computers and Geotechnics, 37 (2010) 977-990.
[40] S.Y. Wang, D.H. Chan, K.C. Lam, S.K.A. Au, Effect Of Lateral Earth Pressure Coefficient On Pressure Controlled Compaction Grouting In Triaxial Condition, Soils and Foundations,
50 (2010) 441-445.
[41] W.M. Yan, X.S. Li, Mechanical response of a medium-fine-grained decomposed granite in Hong Kong, Eng. Geol. , 129 (2012) 1-8.
[42] J.H. Yin, Influence of relative compaction on the hydraulic conductivity of completely decomposed granite in Hong Kong, Canadian Geotechnical Journal, 46 (2009) 1229-1235.
[43] J.H. Yin, L.M. Chu, L.J. Su, Millpress, Study on the interface shear strength of soil nails in a completely decomposed granite soil by laboratory pullout tests and large-size direct shear
box tests, 2005.
[44] J.H. Yin, W.H. Zhou, M. Kumruzzaman, C.M. Cheng, A Rigid-Flexible Boundary True Triaxial Apparatus for Testing Soils in a Three-Dimensional Stress State, Geotech. Test. J. , 34
(2011) 265-272.
[45] Z.Q. Yue, Y.J. Shang, R.L. Hu, X.B. Tu, Five test methods for porosity of completely decomposed granite in Hong Kong, Int. J. Rock Mech. Min. Sci. , 41 (2004) 393-394.
[46] T.L.T. Zhan, G.W. Jia, Y.M. Chen, D.G. Fredlund, H. Li, An analytical solution for rainfall infiltration into an unsaturated infinite slope and its application to slope stability analysis, Int.
J. Numer. Anal. Methods Geomech. , 37 (2013) 1737-1760.
[47] L.M. Zhang, S.M. Dasaka, Uncertainties in Geologic Profiles versus Variability in Pile Founding Depth, Journal of Geotechnical and Geoenvironmental Engineering, 136 (2010) 1475-
1488.
106
TOPICS & SYLLABUS:
I. Introduction and Geotechnical
Topic 1: Site Investigation Properties of Soil
Topic 2: Slope Stability II. Natural Soil Deposits and Subsoil
Exploration
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
1
Civil Engineering 2019/2020
Subject: Geotechnical Design (CSE40403)
Site investigation
Part II- Natural Soil Deposits and Subsoil
Exploration
by Dr. Zhenyu YIN
Office: ZS915, Tel: 3400-8470
Email: zhenyu.yin@polyu.edu.hk
2
Main content:
3
I. Natural Soil Deposits Atmosphere:
Air, wind, rain/sunny
life Rocks
Minerals
Weathering
Hydrosphere
Soils
Transported soils
Residual soils
Peats/organic soils
4
I.1 Soil origin
Most of the soils that cover the ❖ Changes in temperature result in expansion
earth are formed by the weathering and contraction of rock due to gain and loss
of various rocks. of heat.
Two general types of weathering: ❖ Continuous expansion and contraction will
result in the development of cracks in rocks.
(1) Physical weathering, is a process
by which rocks are broken down ❖ Flakes and large fragments of rocks are split.
into smaller and smaller pieces by
❖ Frost action is another source of physical
physical forces without any change
weathering of rocks.
in the chemical composition.
❖ Water can enter the pores, cracks, and other
(2) Chemical weathering, is a
openings in the rock. Continuous freezing
process of decomposition or
and thawing will result in the breakup of a
mineral alteration in which the
rock mass.
original minerals are changed into
something entirely different. ❖ Exfoliation by which rock plates are peeled
off from large rocks by physical forces.
❖ Action of running water, glaciers, wind,
ocean waves, etc.
5
I.1 Soil origin
Most of the soils that cover the
earth are formed by the weathering
of various rocks.
Two general types of weathering:
(1) Physical weathering, is a process
by which rocks are broken down
into smaller and smaller pieces by
physical forces without any change
in the chemical composition.
(2) Chemical weathering, is a
process of decomposition or
mineral alteration in which the
Common minerals in igneous rocks are quartz,
original minerals are changed into
feldspars, and ferromagnesian minerals. The
something entirely different.
decomposed products of these minerals due to
chemical weathering listed in the above Table.
6
I.2 Transported soils
7
I.2 Transported soils
(1) Gravity transported soil
Most rock weathering is a combination of (2) Lacustrine (lake) deposits
physical and chemical weathering. (3) Alluvial or fluvial soil deposited
by running water
❖ Soil produced by the weathering of rocks
can be transported by physical processes (4) Glacial deposited by glaciers
to other places. The resulting soil deposits (5) Aeolian deposited by the wind
are called transported soils.
❖ In contrast, some soils stay where they
were formed and cover the rock surface
from which they derive. These soils are
referred to as residual soils.
❖ There are peats and organic soils, which
derive from the decomposition of organic
materials.
8
I.2 Transported soils
Landslides in HK
after a rainstorm
9
I.2 Transported soils
(1) Gravity transported soil
Most rock weathering is a combination of (2) Lacustrine (lake) deposits
physical and chemical weathering. (3) Alluvial or fluvial soil deposited
by running water
❖ Soil produced by the weathering of rocks
can be transported by physical processes (4) Glacial deposited by glaciers
to other places. The resulting soil deposits (5) Aeolian deposited by the wind
are called transported soils.
❖ In contrast, some soils stay where they
were formed and cover the rock surface
from which they derive. These soils are
referred to as residual soils.
❖ There are peats and organic soils, which
derive from the decomposition of organic
materials.
10
I.2 Transported soils
(1) Gravity transported soil
Most rock weathering is a combination of (2) Lacustrine (lake) deposits
physical and chemical weathering. (3) Alluvial or fluvial soil deposited
by running water
❖ Soil produced by the weathering of rocks
can be transported by physical processes (4) Glacial deposited by glaciers
to other places. The resulting soil deposits (5) Aeolian deposited by the wind
are called transported soils.
❖ In contrast, some soils stay where they
were formed and cover the rock surface
from which they derive. These soils are
referred to as residual soils.
❖ There are peats and organic soils, which
derive from the decomposition of organic
materials.
11
I.2 Transported soils
Cross section of a braided-stream deposit
Fine sand
Gravel
Silt
Coarse sand
12
I.2 Transported soils
(1) Gravity transported soil
Most rock weathering is a combination of (2) Lacustrine (lake) deposits
physical and chemical weathering. (3) Alluvial or fluvial soil deposited
by running water
❖ Soil produced by the weathering of rocks
can be transported by physical processes (4) Glacial deposited by glaciers
to other places. The resulting soil deposits (5) Aeolian deposited by the wind
are called transported soils.
❖ In contrast, some soils stay where they
were formed and cover the rock surface
from which they derive. These soils are
referred to as residual soils.
❖ There are peats and organic soils, which
derive from the decomposition of organic
materials.
13
I.2 Transported soils
(1) Gravity transported soil
Most rock weathering is a combination of (2) Lacustrine (lake) deposits
physical and chemical weathering. (3) Alluvial or fluvial soil deposited
by running water
❖ Soil produced by the weathering of rocks
can be transported by physical processes (4) Glacial deposited by glaciers
to other places. The resulting soil deposits (5) Aeolian deposited by the wind
are called transported soils.
❖ In contrast, some soils stay where they
were formed and cover the rock surface
from which they derive. These soils are
referred to as residual soils.
❖ There are peats and organic soils, which
derive from the decomposition of organic
materials.
14
I.3 Residual soils
15
I.4 Peats and organic soils
16
In summary
Rocks
Minerals
19
II.1 Subsurface exploration program
The purpose of subsurface exploration is to obtain information that will aid the
geotechnical engineer in
❖ Selecting the type and depth of foundation suitable for a given structure.
❖ Evaluating the load-bearing capacity of the foundation.
❖ Estimating the probable settlement of a structure.
❖ Determining potential foundation problems (e.g., expansive soil, collapsible soil,
sanitary landfill, etc.).
❖ Determining the location of the water table.
❖ Predicting the lateral earth pressure for structures such as retaining walls, sheet pile
bulkheads, and braced cuts.
❖ Establishing construction methods for changing subsoil conditions
20
II.1 Subsurface exploration program
Subsurface exploration comprises several steps:
❖ Collection of Preliminary Information ✓ Information regarding the type of
structure to be built and its general use;
❖ Reconnaissance
✓ A general idea of the topography;
❖ Site Investigation ✓ Type of soil to be encountered near and
around the proposed site.
21
II.1 Subsurface exploration program
Subsurface exploration comprises several steps:
❖ Collection of Preliminary Information
❖ Reconnaissance
❖ Site Investigation
24
In summary
25
Break
26
II.2 Boring and Sampling
Exploratory Borings in the Field
❖ Auger boring is the simplest method of making exploratory boreholes.
❖ Two types of hand auger (< 3 to 5 m): the posthole auger and the helical auger.
❖ Electrical continuous-flight augers are probably the most common method used
for advancing a borehole (up to 60~70 m):
- Rotary drilling, for soils not too hard;
- Percussion drilling, for hard soil and rock.
Posthole Helical
auger auger
27
II.2 Boring and Sampling
Sampling of soils
Two types of soil samples during subsurface exploration: disturbed and undisturbed.
28
II.2 Boring and Sampling
Sampling of soils
29
II.2 Boring and Sampling
AR ?
[1] Split-Spoon Sampling
Split-spoon samplers can
be used in the field to
obtain soil samples that
are generally disturbed,
but still representative.
31
II.2 Boring and Sampling
[3] Sampling with a Thin-Walled Tube
❖ Thin-walled tubes are made of seamless steel and are frequently used to obtain
undisturbed clayey soils.
❖ The most common dimension: outside diameter of 50.8 mm or 76.2 mm.
❖ The bottom end of the tube is sharpened.
Exercise
A thin-walled tube with 50.8 mm outside diameter has a wall thickness of 1.63 mm.
How is the degree of disturbance?
32
II.2 Boring and Sampling
[4] Sampling with a Piston Sampler
When undisturbed soil samples are very soft
or larger than 76.2 mm in diameter, they tend
to fall out of the sampler. Piston samplers are
particularly useful under such conditions.
❖ Initially, the piston closes the end of the
tube. The sampler is lowered to the
bottom of the borehole, and the tube is
pushed into the soil hydraulically, past the
piston.
❖ Then the pressure is released through a
hole in the piston rod.
To a large extent, the presence of the piston
prevents distortion in the sample by not
letting the soil squeeze into the sampling
tube very fast and by not admitting excess
soil. Consequently, samples obtained in this
manner are less disturbed than those 33
obtained by Shelby tubes.
II.2 Boring and Sampling
Observation of Water Tables
❖ If water is encountered in a borehole during
a field exploration, that fact should be
recorded.
❖ In soils with high hydraulic conductivity (e.g.
in sand), the level of water in a borehole will
stabilize about 24 hours after completion of
the boring.
❖ In highly impermeable layers (e.g. in clay),
the water level in a borehole may not
stabilize for several weeks. A piezometer
should be used with periodic checking until
the water level stabilizes.
34
In summary
Exploratory Borings in the Field
✓ Hand augers (< 3 to 5 m): the posthole auger and the helical auger.
✓ Electrical augers (up to 60~70 m): rotary drilling (soft), percussion drilling (hard)
Sampling of soils
✓ Disturbed samples, used for tests of physical properties
✓ Undisturbed samples, used for hydro-mechanical properties
39
II.3 Field Mechanical Testing
0.689
Overconsolidation ratio, OCR of clay N 60
OCR = 0.193
s
0
Effective vertical stress
Unconfined compression strength (UCS) of clay (in MPa)
w −w
CI = L
wL − wp
40
II.3 Field Mechanical Testing Correction coefficient. For example,
Liao and Whitman’s relationship (1986)
Some correlations by N60
pa atmospheric
= CN N 60 CN = pressure
Correction for N60 in Granular Soil N 60 s 0 (101.325 kN/m2)
Effective vertical
stress
Relative Density of Granular Soil uniformity coefficient
‘
(Cubrinovski and Ishihara 1999)
41
II.3 Field Mechanical Testing Correction coefficient. For example,
Liao and Whitman’s relationship (1986)
Some correlations by N60
pa
= CN N 60 CN =
Correction for N60 in Granular Soil N 60 s 0
42
II.3 Field Mechanical Testing
Approximate relation
between N and the
relative density of soil
(GeoGuide 3, HK)
43
https://en.wikipedia.org/wiki/Shear_vane_test
II.3 Field Mechanical Testing
[2] Vane Shear Test (VST)
The VST is used during the drilling operation
to determine the in-situ undrained shear
strength of clay soils—particularly soft clays.
Testing procedure:
(1) The vanes of the apparatus are pushed
into the soil at the bottom of a borehole
without disturbing the soil appreciably.
(2) Torque is applied at the top of the rod to
rotate the vanes at a standard rate (0.1⁰/sec).
This rotation will induce failure in a soil of
cylindrical shape surrounding the vanes.
(3) The maximum torque, T, applied to cause
failure is measured. The undrained shear
strength is estimated by
T D2 H D
cu = K= 1 + (why?) 44
K 2 3H
II.3 Field Mechanical Testing
T
[2] Vane Shear Test (VST)
T D2 H D (why?)
cu = K= 1 +
K 2 3H
Three failure surfaces: side of cylinder, top and bottom
T2 = 2 2 r dr cu r
0
dr
( D 2)
3
D3
D 2
= 4 cu r 2 dr = 4 cu = cu r
0
3 6
In total:
D D3 T D2 H D
T = T1 + T2 = DHcu + cu cu = with K = 1 +
2 6 K 2 3H 45
II.3 Field Mechanical Testing
[2] Vane Shear Test (VST)
The cu is overestimated by VST. (Morris and Williams 1994)
In real design cu needs to be
corrected: cu = cu (VST)
(Bjerrum 1972)
(Hansbo 1957)
(Larsson 1980) 46
https://en.wikipedia.org/wiki/Cone_penetration_test
II.3 Field Mechanical Testing
[3] Cone Penetration Test (CPT)
❖ The test is also called the static penetration test,
and no boreholes are necessary to perform it.
❖ In the original version, a cone with a base area of
10 cm2 was pushed into the ground at a steady rate
(20 mm/sec).
❖ Two measures by the CPT can be obtained:
(a) the cone resistance (qc) to penetration developed
by the cone, which is equal to the vertical force
applied to the cone, divided by its horizontally
projected area;
(b) the frictional resistance (fc) which is equal to the
vertical force applied to the sleeve, divided by its
surface area—actually, the sum of friction and
adhesion.
47
II.3 Field Mechanical Testing
[3] Cone Penetration Test (CPT)
Friction ratio (Fr) is defined as
Frictional resistance f c
Fr = =
Cone resistance qc
Empirical formulations by
Anagnostopoulos et al. (2003):
49
II.3 Field Mechanical Testing
[3] Cone Penetration Test (CPT)
50
II.3 Field Mechanical Testing
[3] Cone Penetration Test (CPT)
(Robertson and Campanella 1983)
Undrained Shear Strength
Preconsolidation pressure
Overconsolidation ratio
51
II.3 Field Mechanical Testing
[3] Cone Penetration Test (CPT) (Anagnostopoulos et al. 2003)
Correlation between qc and N60
52
https://en.wikipedia.org/wiki/Pressuremeter_test
II.3 Field Mechanical Testing
[4] Pressuremeter Test (PMT)
❖ The PMT is an in-situ test conducted in a borehole. It
was originally developed by Menard (1956) to measure
the strength and deformability of soil.
❖ The PMT consists of a probe with three cells (top and
bottom guard cells, the middle measuring cell).
❖ The common dimension of probe has a diameter of 58
mm and a length of 420 mm.
❖ The probe cells can be expanded by liquid/gas. The
guard cells are expanded to reduce the end-condition
effect on the measuring cell (V0 = 535 cm3).
❖ During the expansion, the soil around the borehole is
first pushed back into the initial state (with an
incremental volume v0), and then continues to expand.
❖ The soil is considered to have failed when the total
volume of the expanded cavity is about twice the
volume of the original cavity (V = 2*(V0 + v0)). 53
II.3 Field Mechanical Testing
[4] Pressuremeter Test (PMT)
Zone I represents the reloading
portion during which the soil
around the borehole is pushed back
into the initial state (i.e., the state it
was in before drilling). The pressure
p0 represents the in-situ total
horizontal stress.
Zone II represents a pseudoelastic
zone in which the cell volume
versus cell pressure is practically
linear (with a pressuremeter
modulus Ep). The pressure pf
represents the yield pressure.
Zone III is the plastic zone. The
pressure pl represents the limit
pressure.
54
II.3 Field Mechanical Testing
[4] Pressuremeter Test (PMT)
Pressuremeter modulus
Dp
E p = 2 (1 + m s )(Vo + vm )
Dv
ms = Poisson’s ratio (=0.33)
Preconsolidation pressure
(Kulhawy and Mayne 1990) s c = 0.45 pl
56
II.3 Field Mechanical Testing
[5] Dilatometer Test (DMT)
Three important parameters:
57
In summary
Field testing methods Measures For clayey soils For sandy soils Applied by
Standard Penetration Test (SPT) N60 cu, OCR, UCS Dr, f', Es Blow
Cone Penetration Test (CPT) qc, fc cu, s'c, OCR Dr, f' penetrate
Pressuremeter Test (PMT) p0, pf, pl, Dv cu, s'c, K0, Ep K0, Ep Expand
Dilatometer Test (DMT) p0, p1 cu, OCR, K0, Es K0, Es, f' Expand
58
Break
59
II.4 Geophysical Exploration
Three typical geophysical exploration techniques:
✓ Seismic refraction survey
✓ Cross-hole seismic survey
✓ Resistivity survey
Advantages: rapid evaluation, large areas, less expensive than field testing.
Drawbacks: definitive interpretation is difficult, thus only for preliminary work.
60
II.4 Geophysical Exploration
[1] Seismic Refraction Survey
Purpose: to obtain the thickness of the soil layer and the depth to rock/hard soil at a site.
Principle and procedure:
❖ To create a point Source (P-waves) by a hammer blow or by a small explosive charge;
❖ To record with time the first arrival of P-waves at several other points (e.g., 1, 2, …… );
❖ To estimate the thickness of soil layer using velocities and times.
v=
Es (1 − ms )
( g ) (1 − 2ms )(1 + ms )
61
II.4 Geophysical Exploration
To determine the velocity v of P waves in various layers and the thicknesses
of those layers:
62
II.4 Geophysical Exploration
[2] Cross-Hole Seismic Survey
Purpose: to obtain the shear modulus of
the soil layer at a specified depth.
Principle and procedure:
❖ Two holes drilled into the ground with
a distance L.
❖ A vertical impulse is created at the
bottom of one borehole by means of
an impulse rod, to generate S-waves.
❖ The velocity of S-waves can be
estimated by
L
vs =
t
Gs vs2
vs = Gs =
( g ) g (Stokoe and Woods, 1972)
63
II.4 Geophysical Exploration
[3] Electrical Resistivity Survey
The electrical resistivity of any conducting
material having a length L and an area of
cross section A, defined as Electrical
AR
= resistance
L
Purpose: to obtain the thickness of the soil
layer at a site.
(a) Wenner method
Principle and procedure:
❖ To make four electrodes driven into the
ground and spaced equally along a straight
line (referred as Wenner method).
2 dV Voltage drop
=
I Electrical current
65
II. Subsoil Exploration Report
General requirements of a Report:
1) A description of the scope of the investigation
2) A description of the proposed structure for which the subsoil exploration has been
conducted
3) A description of the location of the site, including any structures nearby, drainage
conditions, the nature of vegetation on the site and surrounding it, and any other
features unique to the site
4) A description of the geological setting of the site
5) Details of the field exploration—that is, number of borings, depths of borings, types of
borings involved, etc.
6) A general description of the subsoil conditions, as determined from soil specimens and
from related laboratory tests, standard penetration resistance and cone penetration
resistance, etc.
7) A description of the water-table conditions
8) Recommendations regarding the foundation, including the type of foundation
recommended, the allowable bearing pressure, and any special construction procedure
that may be needed; alternative foundation design procedures should also be
discussed in this portion of the report
66
9) Conclusions and limitations of the investigations
II. Subsoil Exploration Report
Necessary attachments to the report:
1) A site location map
2) A plan view of the location of the borings with
respect to the proposed structures and those
nearby
3) Boring logs
4) Laboratory test results
5) Other special graphical presentations
67
Exercise 1
A load of 100 kN is carried on a foundation 2 m square at a shallow depth in a soil
mass. The foundation soil has a water table at the top of ground, and a saturated
unit weight of 17 kN/m3.
Please use the rule of ASCE to determine the approximate minimum depth of
boring. (Note: The calculation of stress in soil is simplified by the equation of point
load.) 52
Q 3 1
s z = 2 I P with I P =
z 2 1 + ( r z )2 Q
q = Q/(B*L)
68
Exercise 2
For a Shelby tube, given: outside diameter = 76.2 mm and inside diameter = 73 mm.
What is the area ratio of the tube?
Do2 − Di2
Ar =
Di2
70
Exercise 3
A soil profile is shown in Figure below along with the standard penetration numbers in
the clay layer. Use the following equations to determine the variation of cu and OCR
with depth. What is the average value of cu and OCR?
0.689
N 60
cu = 0.29 N 600.72 pa OCR = 0.193
s
0
= 100 kPa in MPa
Atmospheric pressure
Dp
E p = 2 (1 + m s )(Vo + vm )
Dv
76
Exercise 6
The results of a refraction survey at a site are given in the following table. Please
determine the P-wave velocities and the thickness of the material encountered.
Distance, x (m)
78
TOPICS & SYLLABUS:
Topic 1: Site Investigation
Topic 2: Slope Stability
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
a
1
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
Slope Stability
Sliding distance
Sliding
surface
Time
4
1 Introduction
Landslide in La Conchita,
California, 1995
Champagne, France
Sliding distance
Time
5
1 Introduction
Results of slope instability:
landslides
Slope
Soil Rock
6
1 Introduction
21/5/2016
slope failure at
Sai Kung,
about 2000m3
failure mass
1 Introduction
Rainfall is a very critical
problem in HK. Average
rainfall 1961-1990 is
2214mm per year. From
1961-2007 is 2304mm.
Average 300 slope failures
in the past._by CEDD
Debris flow at Tsing Shan 1990 Debris flood at Route Twisk, Tsuen Wan,
Flume test of landslide in PolyU (using HK-CDG, simulating rainfall, reduced scale)
1 Introduction
Before After
In summary
Stability analysis
✓ Factor of safety
✓ Failure surface
✓ Shear strength
✓ Equilibrium
13
Break
14
2 Factor of safety
t
❖ Total stress approach t f = cu
15
2 Factor of safety
Three FOS can be defined:
1. Moment equilibrium: generally used for the analysis of the rotational landslides.
Considering a circular slip surface, the factor of safety is defined:
Fsm = Mr/Md
3. Shear strength: the actual shear stress mobilized along the failure surface against
the ultimate shear strength as:
tm =tf/Fst with tf = c’+σ’n tanϕ’ or tf = cu
16
2 Factor of safety
Mr = sum of resisting moments
Md = sum of driving moments
❖ Based on moment: Fsm = Mr/Md tf = ultimate shear strength
❖ Based on force: Fsf = Fr/Fd tm = mobilized shear stress
Fr = sum of the resisting forces
❖ Based on shear strength: Fst = tf/tm Fd = sum of the driving forces
Attention:
✓ Without external forces, these definitions are same.
✓ With external forces, they are different.
Mr M + Me (If Mr = 3, Md = 2, Me = 1, Fs = ?)
− r
Md − Me Md
Active support is more commonly adopted,
=
(Me + Mr − Md ) Me > or < 0??? China code for hydraulic structures – adopt passive
(Md − Me ) Md support concept.
17
2 Factor of safety Exercise
Force against sliding
Example 1: Estimate the factor of safety (1) Fsf =
Force of sliding
using three definitions of FOS?
O Rcos(a+b)
Rotation center
R htanb
htan(a+b)
W
tf h
b a
Shear Strength
(2) Fsτ =
Shear stress
Unit weight: g
Undrained shear strength: cu
Moment against rotation
(3) Fsm =
Moment of overturning
18
2 Factor of safety
Force against sliding
Example 1: Estimate the factor of safety (1) Fsf =
Force of sliding
using three definitions of FOS?
t fl
Fs =
Rcos(a+b) W cos (a + b )
O
Rotation center cu h cos (a + b )
=
R htanb 1 2
g h tan (a + b ) − tan ( b ) cos (a + b )
htan(a+b) 2
2c 1
= u
W g h tan (a + b ) − tan ( b ) cos 2 (a + b )
tf h
b a
Shear Strength
(2) Fsτ =
Shear stress
tf t fl
Fs = =
Unit weight: g W cos (a + b ) l W cos (a + b )
Undrained shear strength: cu
Moment against rotation
1 (3) Fsm =
W= h h tan (a + b ) − h tan ( b ) g
2 Moment of overturning
1
= g h 2 tan (a + b ) − tan ( b ) t f lR t fl
Fs = =
2 WR cos (a + b ) W cos (a + b )
19
2 Factor of safety
This is the simplest of the circular arc methods of analysis. Failure is assumed to occur by
rotation of a rigid block of soil on a cylindrical failure surface along which the undrained
shear strength of the soil cu is mobilized. [known as Fellenius Method]
20
2 Factor of safety
x
1) Length of arc: L = Rq
O
q R 2) The tm is average mobilized shear stress along
C B L, so T = tmL = (cu/Fs)Rq
R
W 3) Restoring moment: TR
h
tf 4) Overturning moment: Wx
Alternatively,
tfLR = cuLR = cuR2q
Moment 𝒂𝒈𝒂𝒊𝒏𝒔𝒕 𝒓𝒐𝒕𝒂𝒕𝒊𝒐𝒏
𝑭𝒔 =
Moment of overturning
Wx
21
2 Factor of safety
22
In summary
23
Break
24
3 Methods of slope stability analysis
❖ Limit equilibrium method (LEM) – most common
❖ Limit analysis method (LAM) – energy balance
❖ Finite element/difference method (FEM/FDM)
❖ Distinct element/rigid element method – uncommon
❖ Slip-line method – uncommon
❖ Variational principle – uncommon
25
3 Methods of slope stability analysis Exercise
htan(b)
❑ Example 3a: Estimate FOS and the critical height
of a vertical cut by LEM (total stress approach) W
tf h
Relating to high & steep slopes b
26
3 Methods of slope stability analysis
htan(b)
Wsin(b) 1
W = g h 2 tan ( b ) W
W
2 tf h
Wcos(b) b
Solution:
Step 1: for an arbitrary slip surface (b ), the general FOS is
t l = W cos ( b )
tf h 1
= g h 2 tan ( b ) cos ( b )
Fs cos ( b ) 2 4cu hcr
Fs = =
Fs =
2cu b = 4
⎯⎯⎯ → gh h
g h tan ( b ) cos 2 ( b )
h −2cu cos ( 2 b )
= = 0 cos ( 2 b ) = 0 b =
b g sin ( b ) cos ( b ) 2
4
corresponds to a minimum FOS for general case
27
3 Methods of slope stability analysis
Unsafe theorem
An unsafe solution to the true collapse loads (for the ideal plastic material) can be
found by selecting any kinematically possible failure mechanism and performing
an appropriate work (or energy) calculation. The loads so determined are either
on the unsafe side or equal to the true collapse loads.
Referred as “Upper bound” theorem. The accuracy of the solution depends on
how close the assumed failure mechanism is to the real one.
28
3 Methods of slope stability analysis
htan(b)
❑ Example 3b: Estimate FOS and the critical height
W
of a vertical cut by Upper bound analysis
tf h
b
Solution:
Step 1: for an arbitrary slip surface, under a relative
displacement between the two rigid blocks u, Wsin(b)
1
the dissipation energy is W = g h 2 tan ( b )
2 W
c h Wcos(b)
E p = t lu = u u
Fs cos ( b )
work done by external (body) forces is
u
1 ucos(b)
Ew = Wu cos ( b ) = g h 2 tan ( b ) u cos ( b )
2 usin(b)
Step 2: equilibrium of work gives:
E p = Ew
cu h 1
u = g h 2 tan ( b ) u cos ( b )
Fs cos ( b ) 2
2cu
Fs =
g h tan ( b ) cos 2 ( b ) The rest is same as the example 3a.
29
3 Methods of slope stability analysis
Safe theorem
If a statically admissible stress field covering the whole soil mass can be found,
which nowhere violates the yield condition, then the loads in equilibrium with the
stress field are on the safe side or equal to the true collapse loads.
Referred as “Lower bound” theorem. The accuracy of the solution depends on
how close the assumed stress field is to the real one.
30
3 Methods of slope stability analysis
❑ Example 3c: Estimate FOS and the critical height
of a vertical cut by Lower bound analysis
h
Solution:
In stress field, to find an equilibrium distribution of
sz = g *z
stress that balances the body stress (geostatic stress):
cu sh =
sz = gh = 2
Fs K0*s’z+u
2cu
Fs =
gh t
sz = gh
cu
Fs =1 2cu cu/Fs
⎯⎯⎯
→ hcr =
g sz sh sz s’n
31
3 Methods of slope stability analysis
z
Unit weight: g
shear strength: c’, f’
(water table below
the bedrock)
t f t f t fl
Fs = = = = ......
t T l T
36
Exercise 2: The figure below shows a selected slice of soil block Tutorials
above an assumed circular slip surface.
O d
F A’ B’
R
y
b
a
x
Unit weight: g A B
shear strength: c’, f’
a
b
(a) Derive formulations for the driving moment Md and (x1, y’1) (x2, y’2)
resisting moment Mr of the selected slice of soil
(xm, y’m)
block, assuming AB a straight-line.
Given coordinates: O(11,20), F(14,14),and soil properties:
c’ =10 kPa, f’=35⁰, g = 11 kN/m3, z
the given slice
(b) For the slice from AA’(x=22) to BB’(x=24), calculate
the Md and Mr?
(c) For the slice from AA’(x=6) to BB’(x=8), calculate the (x2, y2)
Md and Mr? (x1, y1) (xm, ym)
37
Solution of Exercise 2(a) Tutorials
(a) Derive formulations for the driving O d
moment Md and resisting moment Mr of F A’ B’
the selected slice of soil block, assuming R
AB a straight-line. y
b
a
x
A B
a
b
W = g zb
(x1, y’1) (x2, y’2)
N = W cos a = g zb cos a T
T = W sin a = g zb sin a (xm, y’m)
N W
y1 = 20 − 22.82542 − ( 22 − 11) = 0
2
(xm, y’m)
b = x2 − x1 = 2
xm − xc −1 23 − 11
a = sin −1 = sin = 0.55 (= 31.72 )
o
W = g zb a b
N = W cos a = g zb cos a T
(x2, y’2)
T = W sin a = g zb sin a
W N
(xm, y’m)
(x1, y’1)
Tf
z
N g zb cos b T = t f l = ( c + s n tan f ) b cos a
s n = = = g z cos 2 a f
l b cos b = c b cos a + g zb cos a tan f (x1, y1)
(xm, ym)
(x2, y2)
M r = T f R = ( c b cos a + g zb cos a tan f ) R
M d = Wd = g zbR sin a ( or M d = TR = g zb sin a R ) 40
Solution of Exercise 2(c) Tutorials
d O
Given coordinates: O(11,20), F(14,14),and soil R
properties: c’ =10 kPa, f’=35⁰, g = 11 kN/m3, F
(c) For the slice from AA’(x=6) to BB’(x=8),
calculate the Md and Mr? B’
A’
y
M r = ( c b cos a + g zb cos a tan f ) R a
x
M d = g zbR sin a
A B
R = x + y = 11 + 20 = 22.8254
2 2 2 2
c c
a b
yi = yc − R 2 − ( xi − xc )
2
T (x2, y’2)
y1 = 20 − 22.82542 − ( 6 − 11) = −2.27
2
W N (xm, y’m)
y2 = 20 − 22.8254 − ( 8 − 11) = −2.63
2 2
(x1, y’1)
At the slope: y = x tan ( 4 ) = x y1 = 6, y2 = 8
Tf
z = ym − ym = ( 6 + 8 ) 2 − ( −2.27 − 2.63) 2 = 9.45 z
b = x2 − x1 = 2 (x1, y1)
(xm, ym)
x −x −1 7 − 11
a = sin m c
−1
= sin = −0.18 (= −10.09 )
o
❖ Known equations: 4n
No. of unknowns > No. of equations
❖ Unknowns: 6n-2 or 5n-2
For example:
add n-1 equations (n-1 interslices),
hence there will 1 equation more than unknowns.
Assumptions cannot be proved, still problem!
45
4 LEM-based method of slices
LEMs according to equilibrium conditions
O
X1 E2
E1
X2
46
4.1 Translational sliding failure
e.g., a long gentle soil
slope on hard soil/bedrock.
l=b/cosb
Infinite Slope
Consider the infinite slope as shown, the tangential force T down the slope is:
T = W sin b = g zb sin b (1)
The tangential stress τ down the slope is: T zb sin
= = = z sin cos (2)
l b cos
The normal total force N on the slip surface is:
N = W cos b = g zb cos b
(3)
The normal total stress σ on the segment is : N g zb cos b
sn = = = g z cos 2 b
l b cos b (4)
The normal effective stress σ’ on the segment is then:
s 'n = s − u = g z cos 2 b − ( z − hw ) g w cos 2 b = (g z − g w z + g w hw ) cos 2 b (5)
The shearing resistance τf at the base of the segment will be:
t f = c '+ s n tan f ' = c '+ ( g z − g w z + g w hw ) cos 2 b tan f ' (6)
47
4.1 Translational sliding failure Exercise
(b) Wet cohesionless slope c’= 0, hw = 0, this is the most critical condition, when the water
table is at ground level and seepage is occurring downslope. How is Fs?
48
4.1 Translational sliding failure
The factor of safety is given by:
Fs = t f t (7)
For the general case:
c '+ ( g z − g w z + g w hw ) tan f cos 2 b
Fs = (8)
g z sin b cos b
(b) Wet cohesionless slope c’= 0, hw = 0, this is the most critical condition, when the water
table is at ground level and seepage is occurring downslope. From eq.(8):
0 + ( g z − g w z + g w 0 ) tan f cos 2 b g tan f 'cos b g tan f '
Fs = = = (10)
g z sin b cos b g sin b g tan b
b ≈ ϕ’/2 for fully saturated slope (when Fs ≈ 1), a major reduction in slope angle.
49
Exercise
Example 5
The infinitely long slope is shown in the figure. The thickness of the soil is 4 m, the bedrock is
below it, and the dip Angle of the slope is β =17o. The unit weight of soil is 20 kN/m3.
Assuming that the sliding surface in the figure is parallel to the slope.
(1) For a dry cohesionless slope, please calculate the minimum factor of safety for the slope.
(soil parameter: f’=35o, c’=0)
(2) For a wet cohesionless slope, please calculate the minimum factor of safety for the slope.
(soil parameter: f’=35o, c’=0)
(3) For a cohesionless slope, if the water table is 1m, calculate the minimum safety factor for
slope.
(4) For s cohesive slope (soil parameter: f’=25o, c’=20kPa), calculate the minimum safety
factor when the water table is 1m below.
l = b/cosb 50
Solution of Example 5
(1) For a dry cohesionless slope, please calculate the tan f ' tan ( 35 )
Fs = = =2.29
minimum factor of safety for the slope. (soil tan b tan (17 )
parameter: f’=35o, c’=0)
(2) For a wet cohesionless slope, please calculate the g tan f ' ( 20 − 9.81) tan ( 35 )
Fs = = =1.17
minimum factor of safety for the slope. (soil g tan b 20 tan (17 )
parameter: f’=35o, c’=0)
(3) For a cohesionless slope, if the water table is 1m,
calculate the minimum safety factor for slope. β =17o
(4) For s cohesive slope (soil parameter: f’=25o,
c’=20kPa), calculate the minimum safety factor
when the water table is 1m below.
g = 20 kN/m3
l = b/cosb
tan f ( g z − g w z + g w hw ) tan ( 35 ) ( 20 4 − 9.81 4 + 9.81 1)
Fs = = = 1.45
g z tan b tan (17 ) 20 4
c '+ tan f cos 2 b ( g z − g w z + g w hw ) 20+ tan ( 25 ) cos 2 (17 ) ( 20 4 − 9.81 4 + 9.81 1)
Fs = = = 1.86
g z sin b cos b sin (17 ) cos (17 ) 20 4 51
4.1 Translational sliding failure
http://geoinvention.com/teaching.html
52
Exercise
Example 6
Parallel seepage flow,
U g wb ( z − hw ) cos b
u= = = g w ( z − hw ) cos 2 b
l l
Attention:
✓ For parallel seepage flow, the hydraulic gradient height is “(z-hw)cosb ”
perpendicular to the slip surface;
✓ For horizontal seepage flow, the hydraulic gradient height is directly “(z-hw)”,
same value at any direction.
53
Example 6
Parallel seepage flow,
U g wb ( z − hw ) cos b
u= = = g w ( z − hw ) cos 2 b
l l
l=b/cosb 55
Break
56
4.2 Swedish method (or Ordinary method)
X1 E2
E1
X2
57
4.2 Swedish method (or Ordinary method)
O
hi
Wi
A b
(0,0) li
a+ - ai
Sm (xi+1, y’i+1)
(xmi, ymi)
59
4.2 Swedish method
Wi
A
(0,0) li
a+ - ai
Sm
(xmi,ymi) N
Fs =
( cb sec a + b h g cos a
i i i i i tan f )
b h g sin a i i i For clayey soils (undrained),
possibly below the toe! 60
Solution of Example 7
The shear stress at base of slice = cu. The shear resistance moment is equal to the force
along the length of the circular LAB are multiplied by the radius of the circle, R:
Shear resistance moment = cuLABR
The overturning moment is Wd. W is the total weight of the segment determine from γA, γ
is the bulk unit weight of the soil and A is the area of the segment. d is the horizontal
distance from the circle center to the centroid of the failure mass. The factor of safety is
then:
cu LAB R
Fs =
Wd
The determination of W and d can entail lengthy calculation, so it will be more convenient
to split the segment into a number of slices. Since di for i segment is given by Rsinαi, and
the arc length of i segment is approximated by chord length as bi/cosαi = bisecαi, the factor
of safety is then given by
Fs =
Rcu li = Rcubi sec ai = cubi sec ai
Wdi Wi R sin a i bi hig sin ai
Alternatively,
Fs =
( l c + b h g cos a tan f ) f
i i i i l c = c b sec a
= 0 , c = cu
⎯⎯⎯⎯→ F = i u u i i
As the condition of limiting equilibrium develops with the factor of safety close to 1, a tension
crack may form near the top of the slope through which no shear strength can be developed
and if it fills with water. A horizontal hydrostatic force Pw will increase the disturbing moment
Pwxw. The factor of safety will be further reduced because of the shorter length of circular arc
along which shearing resistance can be mobilized. The depth of a tension crack can be taken
as:
(xo, yo) Pw
xw
O
(xc, yc)
62
4.2 Swedish method
Fs =
c b sec a
u i i
W sin a
i i
Solution of Example 8
s z = g zc = 2cu
1
✓ the hydrostatic force is Pw = g w zc2 (Depth of water is zc)
2 zc
(In clayey soils, dissipation of c
2
So, good covers (also should be light) are necessary, like by grass/plants! 63
4.2 Swedish method Exercise
64
4.2 Swedish method
Solution of Example 9
For moment equilibrium, ignoring inter-slice forces, the overturning moment produced by the
forces S must be balanced by the resisting moment of the mobilized shear strength forces T
(=τl= τfl/Fs): Si R = Ti R
t f li
Wi sin a i R = R
Fs
Assume Fs is the same for all slices:
Fs =
t l f i
W sin a
i i
The shear strength of the soil in terms of effective stress is given by Mohr-
Coulomb equation as (in force, not stress)
(s n = g h cos2 a , l = b cos a )
T = t f l = ( c + s n tan f ) l = cl + g hb cos a tan f (like parallel seepage flow)
Then,
Fs =
( l c + b h g cos a
i i i i tan f )
b h g sin a
i i i
and the subscript denote that the factor of safety is determined from a moment equilibrium
65
4.2 Swedish method (Drained analysis in saturated slope)
N ' = W cos a − ul
( s ' = s − u, multiply l to both sides)
or N=Wcosα expressed in total force
c ' LAB + tan f ' (Wi cos a i − uli ) l c + ( b h g cos a − g h b sec a ) tan f
Fm = Fs =
i i i i w i i i
W sin a
i i b h g sin a
i i i
Values of W, α and u can be determined for each slice and presented in a table, the
summations are obtained by addition in the columns. This solution underestimates the
factor of safety compared with more accurate methods of analysis with an error of up
to 20%, and will therefore, be conservative. 66
di
LEM – Swedish Method O
(xc,yc) rR
bi B (xf,yf)
Moment against 𝒓𝒐𝒕𝒂𝒕𝒊𝒐𝒏 Rr
𝑭𝒔 =
Moment of overturning (xmi,y’mi)
hi
Wi
A
(0,0) li
a+ - ai
Sm
(xmi,ymi) N
(or P)
General formulation for non-circular surface The moment due to N is zero because the
direction of N towards O in circular
surface, but P not in non-circular surface!
67
4.3 Method of Bishop Widely accepted by many codes in different countries!
(Bishop, 1955)
In Swedish method, T is the maximum force corresponding to (X1=X2, E1≠E2)
shear strength. So, T is not real applied value when Fsm > 1.
According to Fsf = Tmax/T, the real applied value should be
T = Tmax/Fsf ≈ Tmax/Fsm.
T = cl + N tan f Fsm (Mohr-Coulomb)
According to the moment definition for Fs, we have [try the derivation by yourself]
http://geoinvention.com/teaching.html
70
In summary O
Fs =
( l c + b h g cos a tan f )
i i i 1 i Special cases:
b h g sin a
i i 2 i ✓ Undrained condition c = cu , f = 0, g 1 = g 2 = g sat
1.47
Principle:
The circular surface, which gives the lowest factor of safety, is the critical circle with FOS.
General method:
A number of trial circles must then be analyzed in the same way but with difference circle
centers and different points where the circle cuts the slope. Use each grid point, try 10-20
failure surfaces, the smallest factor of safety is assigned as the factor of safety for that grid
point, then form the contour of factors of safety, and determine the minimum FOS.
73
5 Design of slope stability
Local and global minimum in general
There are several minimum values C, D, E, B, and E is the global minimum which is the
smallest within the solution domain. Points C and D possess the property of gradient =0 and
are called local minimum.
Using classical mathematics, unless the initial solution is close to E, points C or D may be
obtained. For multi-variables problem, search for global minimum will be a difficult
problem. If B is the smallest value, the traditional gradient method will fail to obtain it.
In slope stability analysis, we need to find the global minimum/critical solution for design.
74
5 Design of slope stability
A good way:
To use the artificial intelligence based methods in slope stability,
to overcome the difficulty of nearly infinite possibilities by
Heuristic optimization method.
75
5 Design of slope stability
Min
76
5 Design of slope stability
Competition to
next population
Start
(can be >1)
77
5 Design of slope stability
78
5 Design of slope stability
❖ Stochastic OA: Genetic algorithm (GA)
79
5 Design of slope stability
All these possible failures should be considered in the slope stability analysis.
o
80
5 Design of slope stability
http://geoinvention.com/down_Slop.html
Development of a design tool for slope stability analysis using MATLAB
81
5 Design of slope stability
Stability chart which gives the minimum factor
of safety of a simple slope without the use of a
program, enable simple and quick stability
assessment.
❖ For the sake of simplicity, most charts have
been developed for homogeneous slope and
very simple geometry.
❖ In the case of non-homogeneity of the soil
strata, average parameters should be
evaluated.
83
5 Design of slope stability
SLOPE/W, a Slope Stability Analysis Software
(https://www.geoslope.com/support/support-resources/tutorial-
videos/geostudio-2018)
84
5 Design of slope stability
For more complicated cases, the FEM is recommended for design.
85
5 Design of slope stability
Recommended factors of safety (by GEO, Hong Kong, 1984)
The definitions of factor of safety, the method of analysis and the accuracy
in the optimization search are never specified in either the Hong Kong code,
Euro code or any other engineering codes.
For very important slopes, the monitoring is needed.
86
5 Design of slope stability
Problems with Slope stability analysis: inaccurate
1. Simplified methods satisfy either force (more common) or moment equilibrium
but not both.
2. Fail to converge in the solution of the factor of safety is sometimes obtained from
commercial programs for ‘rigorous’ methods which satisfy both force and moment
equilibrium.
3. If a solution converged in the analysis, there is still a possibility that the answer is
not acceptable (e.g. global minimum search).
4. All failure mass is 3D, but 2D analysis is performed at present due to various
constraints. The difficulties associated with true 3D analysis are : (1) sliding
direction, (2) satisfaction of all 3D force and moment equilibrium, (3) great amount
of computational geometry calculations is required.
5. Soil strength affected by heavy rainfall with a long recurrency period, additional
stabilization from the presence of vegetation and soil suction, etc. is usually not
considered.
For very important slopes, the monitoring is needed.
87
5 Design of slope stability
(Passive support)
88
5 Design of slope stability
Basic Elements of a
Soil-nailed System
90
5 Design of slope stability
For a typical slope design works, the following items are required to be considered:
❖ Detailed geological conditions, topography, surrounding buildings, drainage, loadings,
vegetation.
❖ Soil and rock conditions and strength parameters, water table.
❖ For sandy soil slopes in HK, non-circular failures are measured from past failures which
should be adopted in the analysis.
❖ Design of slope stabilization measures (usually in form of soil nails in HK) or even
trimming of slope to a gentler profile.
❖ Permanent drainage, surfacing (vegetation), net or other accessories.
91
In summary
1
Types of foundation Foundations
2
Design Philosophy
Limit State
Defined as conditions under which a structure or its component members no
longer perform their intended functions.
Whenever a structure or part of a structure fails to satisfy one of its intended
performance criteria, it is said to have reached a limit state
✓ The foundation was a reinforced concrete raft 2 ft thick and founded at a depth of 12 ft
below the ground surface. (1 ft=0.3048 m)
✓ The weight of the silo was 20,000 tons, which was 42.5 percent of the total weight, when it
was filled.
✓ Filling the silo with grain started in September 1913, and in October when the silo contained
875,000 bushels, and the pressure on the ground was 94 percent of the design pressure, a
vertical settlement of 1 ft was noticed.
✓ The structure began to tilt to the west and within twenty-four hours was at an angle of
26.9° from the vertical, the west side being 24 ft below and the east side 5 ft above the
4
original level.
Design Philosophy
Limit State
Defined as conditions under which a structure or its component members no
longer perform their intended functions.
Whenever a structure or part of a structure fails to satisfy one of its intended
performance criteria, it is said to have reached a limit state
Pisa tower
✓Construction in 1173
Gallileo
✓Inclination: 5.5o
✓Settlement: 5m
Clay
6
Design Philosophy
Limit State
Defined as conditions under which a structure or its component members no
longer perform their intended functions.
Whenever a structure or part of a structure fails to satisfy one of its intended
performance criteria, it is said to have reached a limit state
8
Ma On Shan, HK 2002
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
Soil h
l Footing
plate
(dense soil)
(dense soil)
10
I.1 General Concept
MPM-DEM modeling of footing
(Liang & Zhao, 2019, DOI:
10.1002/nag.2921)
11
I.1 General Concept pressure qu
What is Ultimate Bearing Capacity qu?
load Qult
The load per unit area of the foundation at which shear failure in soil occurs.
qu Failure mode 1:
General shear failure
dense sand or
stiff cohesive soil (sudden
failure up to ground)
Bigger Dr , more
failure lines!
2BL
B* = B B = width;
B+L L = length
Square foundations? B*=B=L
Circular foundations? B*=B=L=diameter
(Note: L is always greater than/equal to B) 13
I.1 General Concept
What is Ultimate Bearing Capacity qu?
The load per unit area of the foundation at which shear failure in soil occurs.
S
More settlement for
punching failure!
S (Bigger Df)
Less settlement for general failure!
14
In summary
Design Philosophy
- Some notations: Limit State?
Ultimate Limit State (ULS)? Serviceability Limit State (SLS)?
- Limit State Design (factors)
15
Break 16
I.2 Terzaghi’s Bearing Capacity Theory
(for general failure case, strip or continuous footing)
Equation of ultimate bearing capacity qu is derived
considering force equilibrium.
Pp Pp
Arc of a logarithmic spiral Strong assumption:
r = r0 exp ( tan ) =
17
I.2 Terzaghi’s Bearing Capacity Theory
(for general failure case, strip or continuous footing)
Equation of ultimate bearing capacity qu is derived considering force equilibrium.
(v)
Pp b
W = b tan
2 qu = + c tan − tan
b 2 18
I.2 Terzaghi’s Bearing Capacity Theory
(for general failure case, strip or continuous footing)
Equation of ultimate bearing capacity qu is derived considering force equilibrium.
1
qu = cN c + qN q + BN
2
exp ( 2(3 4 − / 2) tan )
N c = tan ( K c + 1) = cot
− 1 = cot ( N q − 1)
2 cos ( 4 + 2 )
2
1
exp ( 2(3 4 − / 2) tan ) cot x =
N q = K q tan = tan x
2 cos 2 ( 4 + 2 )
K p = 3 tan 2 ( 4 + ( 2 + 33 /180 ) 2 )
20
【Improved accuracy by numerical solution compared to graphical method】
I.2 Terzaghi’s Bearing Capacity Theory 1
(for general failure case, strip or continuous footing) qu = cN c + qN q + BN
2
Terzaghi’s Bearing Capacity factors
Matlab code
clc; % clean display
close all; % close figures
clear; % clean values of variables data1993=importdata('data.dat');
%% phi1993=data1993(:,1);
phi0=0:1:50; Nc1993=data1993(:,2);
phi=phi0/180*pi; Nq1993=data1993(:,3);
Nr1993=data1993(:,4);
Nq=exp(2*(3/4*pi-phi/2).*tan(phi))./(2*(cos(pi/4+phi/2)).^2);
Nc=(Nq-1)./tan(phi); figure(1)
h1=plot(phi0, Nq,'r-');hold on
Kpr=3*(tan(pi/4+(phi+33/180*pi)/2)).^2; h2=plot(phi0, Nc,'b-');hold on
Nr=1/2*(Kpr./(cos(phi)).^2-1).*tan(phi); legend([h1,h2],'Nq','Nc','best');
%%
figure(2)
h1=plot(phi0, Nr,'b-');hold on
h2=plot(phi1993, Nr1993,'r-');hold on
legend([h1,h2],'Nr','Nr-1993','best'); 21
I.2 Terzaghi’s Bearing Capacity Theory 1
(for general failure case, strip or continuous footing) qu = cN c + qN q + BN
2
Terzaghi’s Bearing Capacity factors
exp ( 2(3 4 − / 2) tan )
N c = tan ( K c + 1) = cot − 1 = cot ( N q − 1)
2 cos 2
( 4 + 2 )
(for ’ = 0? Just put 0.0001)
Extending for local shear failure - using reduced (2/3) of tan(friction angle):
2 1
❑ Strip foundation: qu = c N c + qN q + BN
2/3*’
3 2
❑ Square foundation: qu = 0.867cN c + qN q + 0.4 BN
(2/3*1.3
❑ Circular foundation: = 0.867)
qu = 0.867cN c + qN q + 0.3 BN
Alternatively, we can
still use the old table by
the = tan −1 ( 23 tan )
24
I.2 Terzaghi’s Bearing Capacity Theory 2 1
(for local shear failure case, strip/continuous foundation) qu = cN c + qN q + BN
3 2
Terzaghi’s Modified Bearing Capacity factors
Matlab code
clc; % clean display
close all; % close figures
clear; % clean values of variables
%%
phi0=0:1:50;
phi=atan(tan(phi0/180*pi)*2/3);
Nq=exp(2*(3/4*pi-phi/2).*tan(phi))./(2*(cos(pi/4+phi/2)).^2);
Nc=(Nq-1)./tan(phi);
Kpr=3*(tan(pi/4+(phi+33/180*pi)/2)).^2;
Nr=1/2*(Kpr./(cos(phi)).^2-1).*tan(phi);
%%
25
1
I.3 Effect of Groundwater Table qu = cN c + qN q + * BN
e.g. 2
* replaced by = ( sat − w )
q = ( D f − D ) + D
* replaced by = ( sat − w )
q = Df
26
1
I.3 Effect of Groundwater Table qu = cN c + qN q + * BN
e.g. 2
q = Df
* replaced
1
D + ( B − D ) ( for D B )
B
by av =
B-D B
or by av = ( for D B )
27
I.4 Factor of Safety qu − q0 q 1
qall − q0 = qall = u + 1 − q0
FS FS FS
Factor of Safety (FS) is defined in two ways:
❑ Using gross ultimate bearing capacity qu ❑ Using net ultimate bearing
depending on shear strength of soils: capacity qnet(u)=qu-q0
q q qu ( net ) qu − q0
qall = u or FS = u qall ( net ) = =
FS qall FS FS
The qall is gross allowable bearing capacity WF + WS
q0 = D f
A
FS is in the range 3 ~ 6
Dead and live loads above the footing
according to the design of the superstructure
Qult qu
FS = = =4
Qall qall
Df = 1 m
’ = 20°
𝑐 ′ = 15.2 kN/m2
𝛾 = 17.8 kN/m3
B = L = 1.5 m
29
Example 1 A square foundation is 1.5 m×1.5 m in plan. The soil supporting the foundation
has a friction angle of ’ = 20°, and 𝑐 ′ = 15.2 kN/m2. The unit weight of soil, 𝛾, is
17.8 kN/m3. Determine the allowable gross load on the foundation with a factor
of safety (FS) of 4. Assume that the depth of the foundation (𝐷𝑓 ) is 1 m and that
general shear failure occurs in the soil.
(a) How is total allowable gross load?
Solution:
(deg)
30
Example 1 A square foundation is 1.5 m×1.5 m in plan. The soil supporting the foundation
has a friction angle of ’ = 20°, and 𝑐 ′ = 15.2 kN/m2. The unit weight of soil, 𝛾, is
17.8 kN/m3. Determine the allowable gross load on the foundation with a factor
of safety (FS) of 4. Assume that the depth of the foundation (𝐷𝑓 ) is 1 m and that
general shear failure occurs in the soil. q = 130 kN/m 2 , Q = 292.5 kN
all all
(b) For the same cost of concrete, how will be the Qall for a circular footing?
(c) To achieve the same Qall, how should be the diameter for a circular footing?
(deg)
Terzaghi’s Bearing Capacity equations
1
qu = cN c + qN q + BN
2
qu = 1.3cN c + qN q + 0.4 BN
qu = 1.3cN c + qN q + 0.3 BN
Qult qu
FS = = =4
Qall qall
Df = 1 m
’ = 20°
𝑐 ′ = 15.2 kN/m2
𝛾 = 17.8 kN/m3
B = L = 1.5 m 31
Example 1 A square foundation is 1.5 m×1.5 m in plan. The soil supporting the foundation
has a friction angle of ’ = 20°, and 𝑐 ′ = 15.2 kN/m2. The unit weight of soil, 𝛾, is
17.8 kN/m3. Determine the allowable gross load on the foundation with a factor
of safety (FS) of 4. Assume that the depth of the foundation (𝐷𝑓 ) is 1 m and that
general shear failure occurs in the soil. q = 130 kN/m 2 , Q = 292.5 kN
all all
(b) For the same cost of concrete, how will be the Qall for a circular footing?
(c) To achieve the same Qall, how should be the diameter for a circular footing?
(d) Assuming that the ground water table is located 1 m below the bottom of the footing,
determine the allowable gross load per unit area with Fs = 3 and sat = 19 kN/m3.
(deg)
Terzaghi’s Bearing Capacity equations ( w = 9.81 kN/m3)
1
qu = cN c + qN q + BN
2
qu = 1.3cN c + qN q + 0.4 BN
qu = 1.3cN c + qN q + 0.3 BN
Qult qu
FS = = =3
Df = 1 m
Qall qall
’ = 20°
𝑐 ′ = 15.2 kN/m2
𝛾 = 17.8 kN/m3
D=1m
𝛾𝑠𝑎𝑡 = 19 kN/m3
33
B = L = 1.5 m
Example 1 A square foundation is 1.5 m×1.5 m in plan. The soil supporting the foundation
has a friction angle of ’ = 20°, and 𝑐 ′ = 15.2 kN/m2. The unit weight of soil, 𝛾, is
17.8 kN/m3. Determine the allowable gross load on the foundation with a factor
of safety (FS) of 4. Assume that the depth of the foundation (𝐷𝑓 ) is 1 m and that
general shear failure occurs in the soil. q = 130 kN/m 2 , Q = 292.5 kN
all all
(d) Assuming that the ground water table is located 1 m below the bottom of the footing,
determine the allowable gross load per unit area with Fs = 3 and sat = 19 kN/m3.
1
* replaced by av = D + ( B − D ) ( for D B )
B
34
Example 1 A square foundation is 1.5 m×1.5 m in plan. The soil supporting the foundation
has a friction angle of ’ = 20°, and 𝑐 ′ = 15.2 kN/m2. The unit weight of soil, 𝛾, is
17.8 kN/m3. Determine the allowable gross load on the foundation with a factor
of safety (FS) of 4. Assume that the depth of the foundation (𝐷𝑓 ) is 1 m and that
general shear failure occurs in the soil. q = 130 kN/m 2 , Q = 292.5 kN
all all
(e) Assuming that local shear failure occurs in the soil supporting the foundations.
(deg)
Terzaghi’s Bearing Capacity equations
1
qu = cN c + qN q + BN
2
qu = 1.3cN c + qN q + 0.4 BN
qu = 1.3cN c + qN q + 0.3 BN
Qult qu
FS = = =4
Qall qall
Df = 1 m
’ = 20°
𝑐 ′ = 15.2 kN/m2
𝛾 = 17.8 kN/m3
B = L = 1.5 m 35
Example 1 A square foundation is 1.5 m×1.5 m in plan. The soil supporting the foundation
has a friction angle of ’ = 20°, and 𝑐 ′ = 15.2 kN/m2. The unit weight of soil, 𝛾, is
17.8 kN/m3. Determine the allowable gross load on the foundation with a factor
of safety (FS) of 4. Assume that the depth of the foundation (𝐷𝑓 ) is 1 m and that
general shear failure occurs in the soil. q = 130 kN/m 2 , Q = 292.5 kN
all all
(e) Assuming that local shear failure occurs in the soil supporting the foundations.
(deg)
Solution: qu = 0.867cN c + qN q + 0.4 BN
(’ should be reduced to tan-1(2/3*tan’) = 13.64° 2/3*’=13.3°
(Wrong!)
N c = 11.86, N q = 3.88, N = 1.18
qu = 0.867 15.2 11.86 + (1 17.8) 3.88 + 0.4 17.8 1.5 1.18
= 156.3 + 69.1 + 12.6 = 237.96 kN/m 2
Thus:
q 237.96
qall = u = = 59.5 kN/m 2
FS 4
The total allowable gross load:
Qall = qall A(area) = 59.5 ( B B) = 59.5 (1.5 1.5) = 133.9 kN
36
In summary
40
41
Presumed allowable vatical bearing pressure under foundations on horizontal ground
(The Code of Practice for Foundation 2017, Hong Kong)
42
I.5 General Bearing Capacity Equation
N q = tan 2 ( + )e tan
4 2
New Bearing Capacity factors N c = ( N q − 1) cot
N = 2( N q + 1) tan
43
I.5 General Bearing Capacity Equation
N q = tan 2 ( + )e tan
4 2
New Bearing Capacity factors N c = ( N q − 1) cot
N = 2( N q + 1) tan
Matlab code
clc; % clean display
close all; % close figures
clear; % clean values of variables
%%
phi0=0:1:50;
phi=phi0/180*pi;
Nq=exp(2*(3/4*pi-phi/2).*tan(phi))./(2*(cos(pi/4+phi/2)).^2);
Nc=(Nq-1)./tan(phi);
Kpr=3*(tan(pi/4+(phi+33/180*pi)/2)).^2;
Nr=1/2*(Kpr./(cos(phi)).^2-1).*tan(phi);
%%
44
I.5 General Bearing Capacity Equation
1
qu = cN c + qN q + BN
2
B≠L
c’ = effective cohesion
= ’ (effective), (average) or total unit weight g of the soil
B = width of footing (= diameter of a circular footing)
Nc, Nq, N = bearing capacity factors
Fcs , Fqs , F s = shape factors Note:
Fcd , Fqd , F d = depth factors Equations are derived for
plane strain case. One
Fci , Fqi , F i = load inclination factors may need corrections for
from triaxial data.
Attention:
If L B 2 use tr
✓ The Equation is proposed for effective stress
parameters using ’ and c’; If L B 2 use ps = 1.5tr − 17
✓ The Equation can also be adopted for total If tr 34 use tr = ps
stress parameters using and c.
Keep consistent ! 46
Direct shear sn
t,
Triaxial shear
sa
sr
sr
Note:
Equations are derived for If L B 2 use tr
plane strain case. One If L B 2 use ps = 1.5tr − 17
may need corrections for If tr 34 use tr = ps
from triaxial data.
47
N q = tan 2 ( + )e tan
I.5 General Bearing Capacity Equation 4 2
N c = ( N q − 1) cot
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i N = 2( N q + 1) tan
❑ Shape Factors:
(empirical relations based on extensive laboratory tests by De Beer 1970)
B Nq
Fcs = 1 +
L Nc
B
Fqs = 1 + tan
L
B
F s = 1 − 0.4 Put B/L
L
L ( footing length ) B ( width )
(1) For = 0?
Nq = 1
Put equations needed
N c = 5.14
N = 0
(2) For square footing? (3) For continuous footing? 48
N q = tan 2 ( + )e tan
I.5 General Bearing Capacity Equation 4 2
N c = ( N q − 1) cot
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i N = 2( N q + 1) tan
❑ Depth Factors :
(empirical relations based on extensive laboratory tests by Hansen 1970)
If D f B
If D f B
(1) Df = 0?
(2) Df = B?
49
N q = tan 2 ( + )e tan
I.5 General Bearing Capacity Equation 4 2
N c = ( N q − 1) cot
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i N = 2( N q + 1) tan
❑ Inclination Factors:
(empirical relations based on extensive laboratory tests by Meyerhof 1963)
2
o
Fci = Fqi = 1 − o
90
2
F i = 1 − ( = 0 N = 0 )
= load inclination angle from vertical
Attention: qu is vertical!!!
Df
qu BL
Qult (i) cos = Qult Qult (i) =
cos
50
Example 2 Qall(i), not Qall
A square column foundation has to carry a gross allowable total load Q of 150 kN.
The depth of the foundation is 0.7 m. The load is inclined at an angle of 20°to the
vertical. Determine the width of the foundation, B. (factor of safety is 3)
General bearing capacity equation:
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i
51
Example 2
A square column foundation has to carry a gross allowable total load Q of 150 kN.
The depth of the foundation is 0.7 m. The load is inclined at an angle of 20°to the
vertical. Determine the width of the foundation, B. (factor of safety is 3)
Solution:
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i
No mention of water table → no water
c = 0 → qu = 0 + qN q Fqs Fqd Fqi + 12 BN F s F d F i
q = D f = 0.7 18 = 12.6 kN/m 2
N q = tan 2 ( + )e tan
4 2
N c = ( N q − 1) cot
for = 30o → N q = 18.4, N = 22.4 N = 2( N q + 1) tan
B B B
Fqs = 1 +tan = 1 + tan 30o = 1 + 0.577 = 1.577 F s = 1 − 0.4 = 0.6
L B L
Assuming D f B 1:
Df 0.7 0.202
Fqd = 1 + 2tan (1 − sin ) 2
= 1 + 2tan 30o (1 − sin 30o ) 2 =1+ F d = 1
B B B
2 2
o 20 20
2 2
Solution:
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i
qu 44.68 14.89
qall = = 221.2 + + 13.3B 3 = 73.73 + + 4.43B
FS B B
Q 150 150 14.89
qall = all = 2 2 = 73.73 + + 4.43B
Area B B B
(Wrong!)
By trial & error, find B = 1.3 m
Check D f / B = 0.7 /1.3 1.0 → under the assumption
Attention:
❖ Previous water table effect is applicable.
❖ The ultimate or allowable Q with inclination needs to be modified.
54
Break 55
I.6 Eccentrically Loaded Foundations (one-way case)
Calculation of foundation base pressure for e ≤ B/6:
(moment = force*distance)
1 BL 2 B
M = q(ii ) 2
2 2 3 2
6M
q(ii ) = 2
B L
M (i ) (ii )
q = q(i ) q(ii ) =
Q 6M
2 e=
BL B L Q
Q 6M Q 6e Q Q
= 1 = 1 M
BL BQ BL B e
= +
q(ii )
q(i ) q(ii )
Q
6𝑒 q(i ) =
(1 − >0) BL 56
𝐵
I.6 Eccentrically Loaded Foundations (one-way case)
Calculation of foundation base pressure for e ≤ B/6:
(moment = force*distance)
1 BL 2 B
M= q(ii ) 2
2 2 3 2
6M
q(ii ) = 2
B L
M (i ) (ii )
q = q(i ) q(ii ) =
Q 6M
2 e=
BL B L Q
Q 6M Q 6e Q Q
= 1 = 1 M
BL BQ BL B e
𝐵 = +
(𝑒 = )
6 q(ii )
q(i ) q(ii )
Q
6𝑒 q(i ) =
(1 − =0) BL 57
𝐵
I.6 Eccentrically Loaded Foundations (one-way case)
M
Calculation of foundation base pressure for e > B/6: e=
Q
B 1
Vertical force Fv (total area of pressure): e= − B
2 3
1
Fv = BLqmax
2
2Q
Force equilibrium: Fv = Q qmax = 1
BL Fv B
Moment equilibrium: 3
→ Fv has same location as Q B
B 1 3
e= − B B = B − 3e
2 3 2 Tension will cause large penetration
2Q 2Q 4Q (settlements). So, for design the
qmax = = =
BL ( 3 B − 3e) L 3L( B − 2e) case e > B/6 has to be avoided!!!
2 58
I.6 Eccentrically Loaded Foundations (one-way case)
M
e=
M Q
e= (eccentricity)
Q
For e B 6 :
Q 6M Q 6e
qmax = + 2 = 1 +
BL B L BL B
Q 6M Q 6e
qmin = − 2 = 1 −
BL B L BL B
For e B 6 :
4Q
qmax =
3L( B − 2e)
qmin = 0
Qult = qu BL with B = B − 2e
(Empirical formula, not theoretical) 59
I.6 Eccentrically Loaded Foundations (one-way case)
[1] Effective Area Method (Meyerhoff, 1953)
5. To double check the qmax due to applied eccentric load smaller than q’u/Fs?
60
I.6 Eccentrically Loaded Foundations (one-way case)
[2] Prakash and Saran Theory (1971)
qu = cN c (e) Fcs (e) + qN q (e) Fqs (e) + 12 BN (e) F s (e)
1. The influence of effective width length was included in the modified bearing
capacity factors under eccentric loading.
2. To calculate ultimate load, use full area (not effective area) Qult = qu A
3. The bearing capacity factors and factors of shape effect was modified as
L
Fcs (e) = 1.2 − 0.025 and Fcs (e) 1.0; Fqs (e) = 1
B
2
2e B 3 e B
F s (e) = 1 + − 0.68 + 0.43 −
B L 2 B L
62
I.6 Eccentrically Loaded Foundations (one-way case)
[3] Reduction Factor Method (For Granular Soil) (Purkayastha and Char 1977)
(Based on the stability analysis of eccentrically loaded continuous footings
supported by a layer of sand using the method of slices)
qu (eccentric) qu
Rk = 1 − =1 −
qu (centric) qu e
k
qu = qu 1 − a
e
k
B
Rk = a
B
considering Shape, inclination, depth!
* To calculate ultimate load, use full area (not effective area)
e=
M Qult = qu A
Q
63
Example 3
A square foundation is shown in the figure below. Assume that the one-way load
eccentricity e =0.15 m. Determine the ultimate load, Qult. (check e/B < 1/6?)
Bearing capacity equation according to effective area method
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i
64
Solution 1
qu = qN q Fqs Fqd Fqi + BN F s F d F i
With 𝑐 ′ = 0, the q’u becomes 2
with q = D f = ( 0.7 )(18 ) = 12.6 kN m 2
For 𝜙 ′ = 30°, from the Table, 𝑁𝑞 = 18.4 and 𝑁𝛾 = 22.4. (eccentricity effect ignored)
B = B − 2e = 1.5 − 2 ( 0.15 ) = 1.2 m
L = L = 1.5 m
B 1.2
Fqs = 1 + tan = 1 + tan30 = 1.462
L 1.5
Shape factors using B’ and L’,
B 1.2
F s = 1 − 0.4 = 1 − 0.4 = 0.68 Depth factors using B and L.
L 1.5
So 1
qu = qN q Fqs Fqd Fqi + BN F s F d F i
2
1
= (12.6 )(18.4 )(1.462 )(1.135 ) + (18 )(1.2 )( 22.4 )( 0.68 )(1)
2
= 384.7 + 164.50 = 549.2 kN m 2
Hence
Qult = BL ( qu ) = (1.2 )(1.5 )( 549.2 ) 988 kN 65
Example 4
A continuous foundation is shown in the figure below. If the load eccentricity is 0.2
m, determine the ultimate load, Qult, per unit length of the foundation. (a) Use
Meyerhof’s effective area method.
Bearing capacity equation (effective area method)
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i
66
Example 4
A continuous foundation is shown in the figure below. If the load eccentricity is 0.2
m, determine the ultimate load, Qult, per unit length of the foundation. (a) Use
Meyerhof’s effective area method.
(check e/B < 1/6?)
Solution:
For 𝑐 ′ = 0, it gives
1
qu = qN q Fqs Fqd Fqi + BN F s F d F i
2
Where q = (16.5) (1.5) = 24.75 kN/m2.
For 𝜙 ′ = 40°, from the Table, 𝑁q = 64.2 and 𝑁𝛾 =109.41. Also,
B = 2 − ( 2 )( 0.2 ) = 1.6 m
Because the foundation in question is a continuous foundation, 𝐵′ Τ𝐿′ is zero. Hence,
𝐹𝑞𝑠 = 1, 𝑁𝛾𝑠 = 1. From the Table, (For continuous foundation, shape factors are 1)
Fqi = F i = 1
Df 1.5
Fqd = 1 + 2tan (1 − sin ) = 1 + 0.214 = 1.16
2
B 2
F d = 1
1
qu = ( 24.75 )( 64.2 )(1)(1.16 )(1) + (16.5 )(1.6 )(109.41)(1)(1)(1) = 3287.39 kN m 2
2
Consequently Qult = ( B )(1) ( qu ) = (1.6 )(1)( 3287.39 ) 5260 kN 67
Example 4
A continuous foundation is shown in the figure below. If the load eccentricity is 0.2
m, determine the ultimate load, Qult, per unit length of the foundation. (b) use the
theory of Prakash and Saran (1971).
qu = cN c (e) Fcs (e) + qN q (e) Fqs (e) + 12 BN (e) F s (e)
L
Fcs (e) = 1.2 − 0.025 and Fcs (e) 1.0; Fqs (e) = 1
B
2
2e B 3 e B
F s (e) = 1 + − 0.68 + 0.43 −
B L 2 B L
Since 𝑐 ′ = 0
1
Qult = B qN q( e ) + BN ( e )
2
e 0.2
= = 0.1
B 2
Since 𝜙 ′ = 40°and e/B = 0.1, Figures give 𝑁q(𝑒) = 56.09 and 𝑁𝛾(𝑒) ≈ 71.8.
Hence,
1
Qult = 2 ( 24.75 )( 56.09 ) + (16.5 )( 2 )( 71.8 ) = 5146 kN
2
(attention, lack of correction of depth!) 69
Example 4
A continuous foundation is shown in Figure 3.18. If the load eccentricity is 0.2 m,
determine the ultimate load, Qult, per unit length of the foundation. (c) Use the
Reduction Factor Method.
qu = qu (1 − Rk )
k
e
Rk = a
B
70
Example 4
A continuous foundation is shown in Figure 3.18. If the load eccentricity is 0.2 m,
determine the ultimate load, Qult, per unit length of the foundation. (c) Use the
Reduction Factor Method.
Case I Case II
74
(not required)
I.7 Eccentrically Loaded Foundations (two-way case)
Extension2: Continuous Foundation Subjected to Eccentric Inclined Loading
(Saran and Agarwal 1991)
75
Example 5 (not required)
A continuous foundation is shown in the Figure. Estimate the ultimate load, per
unit length of the foundation.
76
Example 5 (not required)
A continuous foundation is shown in the Figure. Estimate the ultimate load, per
unit length of the foundation.
12.5
15
77
In summary
78
TOPICS & SYLLABUS:
Topic 1: Site Investigation
Topic 2: Slope Stability
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
1
Topic 3: Shallow Foundations
Key points of a foundation design:
(a) shall be safe against overall shear failure (Ultimate Limit State);
(b) cannot undergo excessive displacement (Serviceability Limit State).
I. Ultimate Bearing Capacity (ULS) II. Settlement (SLS)
❖ General Concept ❖ Types of Settlement
❖ Terzaghi’s Bearing Capacity Theory ❖ Elastic Settlement
❖ Effect of Groundwater Table ❖ Consolidation Settlement
❖ Factor of Safety
❖ General Bearing Capacity Equation
❖ Eccentrically Loaded Foundations
❖ Special Cases
2
Ma On Shan, HK 2002
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
Sc(p)
Sc(s)
sz=gh s h = K 0s z
sh=? (s h = K 0s z ? ) (x)
(s h = s h + u ) (√)
Example 1:
Please plot the distribution of vertical stress (total/effective) and horizontal
stress (total/effective) along the depth for four cases.
Cas-1 Cas-2 Cas-3 Cas-4
H g, c,
Cas-2
H gsat, c,
6
II.1 Stresses in ground
Example 1:
Please plot the distribution of vertical stress (total/effective) and horizontal stress
(total/effective) along the depth for four cases.
Cas-1
H g, c,
s'z=sz=gH s'h=sh=K0gH
Cas-2
Attention: s h = K 0s z (No)
H gsat, c, s h = s h + u = K 0s z + u K 0s z (Yes )
sz=gsatH sh=s'h+u=K0(gsat-gw)H+gwH
s'z=(gsat-gw)H s'h=K0s'v=K0(gsat-gw)H
7
II.1 Stresses in ground
Example 1:
Please plot the distribution of vertical stress (total/effective) and horizontal stress
(total/effective) along the depth for four cases.
Cas-3
g, c,
H
hw
(assuming g ≈ gsat)
Cas-4
g1, c1, 1
H
hw
h2 g2, c2, 2
(assuming g1 ≈ g1sat)
8
II.1 Stresses in ground
Example 1:
Please plot the distribution of vertical stress (total/effective) and horizontal stress
(total/effective) along the depth for four cases.
Cas-3
g, c, s'z=sz=g(H-hw) s'h=sh=K0g(H-hw)
H
hw
sz=gsatH sh=s'h+u=s'h+gwhw
(assuming g ≈ gsat) s'z=g(H-hw)+(gsat-gw) hw s'h=K0g(H-hw)+K0(gsat-gw) hw
Cas-4
s z1 = s z1 = g 1 ( H − hw ) , s h1 = s h1 = K 0s z1
g1, c1, 1 s'z1(sz1)
H s z 2 = s z1 + ( g 1 − g w )( hw − h2 ) , s h 2 = K 0s z 2
s'z2(sz2)
hw s z 2 = s z1 + g 1 ( hw − h2 ) , s h 2 = K 0s z2 + g w ( hw − h2 )
h2 g2, c2, 2
s'z3(sz3) s z3 = s z 2 + ( g 2 − g w ) h2 , s h 3 = K 0s z3
(assuming g1 ≈ g1sat)
s z 3 = s z 2 + g 2 h2 , s h 3 = K 0s z3 + g w h2 9
II.1 Stresses in ground (induced by surface loading of different geometries)
Forms of external surface load
10
II.1 Stresses in ground (induced by surface loading of different geometries)
(Google map:
https://www.google.com/maps/place/Phase+8+(Block+Z),+The+Hong+Kong+Polytechnic+University/@22
.3066481,114.1783823,18z/data=!3m1!4b1!4m5!3m4!1s0x340400e80cc5a513:0xf260aecd5ec53016!8m
2!3d22.3066456!4d114.1794793?hl=en&authuser=0)
11
II.1 Stresses in ground (induced by surface loading of different geometries)
Recall
❑ Stress due to a Concentrated Load P
Stresses at any point inside homogeneous, x
elastic, and isotropic mediums due to a r
concentrated point load located at the surface:
3Pz 3 3P
s =
z 2 R 5 = y R z
2 52
2 z 1+ ( r z )
2
3P zx 2 1 − 2 R 2 − Rz − z 2 x ( 2 R + z )
2
z
A (x,y,z)
s x = 5 + 3 − 3
2 R 3 R ( R + z ) R ( R + z )2 sz
3P zy 2 1 − 2 R 2 − Rz − z 2 x ( 2 R + z )
2
s y = 5 + 3 − 3 2
2 R 3 R ( R + z ) R ( R + z )
3P xyz (1 − 2 ) xy ( 2 R + z )
xy = = 5 − 2
( )
yx
2 R 3 R 3
R + z
2
= = − 3Pyz
yz zy
2 R 5
2
= = − 3Pxz J. V. Boussinesq
xz zx
2 R 5 (1842-1929, France) 12
II.1 Stresses in ground (induced by surface loading of different geometries)
3Pz 3 3P
Recall s z = 2 R 5 = 2 52
2 z 1+ ( r z )
2
❑ Stress due to a concentrated load
3P zx 1 − 2 R − Rz − z 2 x ( 2 R + z )
2 2 2
s x = 5 + 3 − 3
Matlab code and plot 2 R 3 R ( R + z ) R ( R + z )2
3P zy 1 − 2 R − Rz − z
2 2 2
x ( 2 R + z )
2
clc s y = 5 + 3 − 3
2 R 3 R ( R + z ) R ( R + z )2
close all
3P xyz (1 − 2 ) xy ( 2 R + z )
clear xy = yx = − 2
2 R 5 3 R 3 ( R + z )
%% 2
= = − 3Pyz
p=1; yz zy
2 R 5
FontSize=16; 2
x=0.1:0.1:10; set(1,'color','w');
z=(0.1:0.1:10)'; set(1,'Position',[500 200 400 700]);
% [x,z]=meshgrid(X,Z); axes1 = axes('Parent',figure1);
r=x; h1=contour(z,x,sigz,'LevelStep',0.01,'ShowText','on');hold on;
% R=sqrt(r.^2+z.^2); set(gca,'Ydir','reverse')
%% set(axes1,'BoxStyle','full','Layer','top','XAxisLocation','top');
for i=1:length(x) axis([0 5 0 10]);
for j=1:length(z) set(gca,'Ytick',0:2:10);
r=x(i); legend('{\it\sigma}_{z}')
R=sqrt(r^2+z(i)^2); shading interp ;
r_z(j,i)=r/z(j); shading flat;
I1(j,i)=3/(2*pi)/((1+(r/z(j))^2)^2.5); %r/z(j) xlabel('\itx','FontWeight','bold');
sigz(j,i)=p/z(j)^2*I1(j,i); ylabel('{\itz}','FontWeight','bold');
end % colorbar('peer',axes1);
end set(gca,'Fontsize',FontSize,'FontName','Times new Roman');
set(get(gca,'XLabel'),'FontSize',FontSize,'FontName','Times new Roman');
set(get(gca,'YLabel'),'FontSize',FontSize,'FontName','Times new Roman');13
II.1 Stresses in ground (induced by surface loading of different geometries)
Recall
❑ Stress due to a uniform circular load
3P 3 ( q0 rd dr )
s z = 52
ds z = 52
2 z 2 1+ ( r z ) 2 z 2 1+ ( r z )
2 2
clc 2 B 2
2 −3 2
close all s =
ds z = q0 1 − 1 + ( B 2 z )
clear
z
=0 r =0
%%
p=1; Matlab code and plot?
nu=0.3;
x=0.1:0.1:10;
z=(0.1:0.1:10)';
% [x,z]=meshgrid(X,Z);
r=x; sz
% R=sqrt(r.^2+z.^2);
%%
for i=1:length(x)
for j=1:length(z)
Much simplified from
r=x(i); this program!
R=sqrt(r^2+z(i)^2);
r_z(j,i)=r/z(j);
I1(j,i)=3/(2*pi)/((1+(r/z(j))^2)^2.5); %r/z(j)
sigz(j,i)=p/z(j)^2*I1(j,i);
end
end 14
II.1 Stresses in ground (induced by surface loading of different geometries)
Recall
❑ Stress due to a uniform circular load
3P 3 ( q0 rd dr )
s z = 52
ds z = 52
2 z 2 1+ ( r z ) 2 z 2 1+ ( r z )
2 2
2 B 2 −3 2
s z = ds z = q0 1 − 1 + ( B 2 z )
2
=0 r =0
15
II.1 Stresses in ground (induced by surface loading of different geometries)
Recall
❑ Stress due to a uniform rectangular load
3P 3 ( q0 dxdy )
s z = 52
ds z = 52
2 z 2 1+ ( r z ) 2 z 2 1+ ( r z )
2 2
L B
s z = ds
y =0 x =0
z = q0 I
1 2mn(m 2 + n 2 + 2) −1 2mn m + n + 1
2 2
I= + tan
4 (1 + n 2 )(1 + m 2 ) m 2 + n 2 + 1 m 2 + n 2 + 1 − m 2 n 2
with m = B z , n = L z sz
General case
General case
s z = q0 ( I1 + I 2 + I 3 + I 4 )
17
z m=B z
n=L z
18
s z = q0 I c (A special case: O in the center) m1 = L B , n1 = 2 z B
19
❑ Stress due to a uniform rectangular load dQ p
Example 2: applying the principle of superposition to calculate the stress at M y
A D
Case-1 x l
M b
z
B C
M
A D
Case-2
z
B C
M
A F
Case-3
M
B E
20
C D
❑ Stress due to a uniform rectangular load dQ p
Example 2: applying the principle of superposition to calculate the stress at M y
A G D
Case-1 x l
M
E F b
AEMG+EBHM+MHCF+GMFD z
B H C
M
A G D
Case-2 z
AEMG-BEMH+GMFD-HMFC
B H C
E F
M
A F
Case-3
M
B E ABCM+MDFE
D 21
C
❑ Stress due to a uniform rectangular load dQ p
Example 2: applying the principle of superposition to calculate the stress at M y
A E
x l
Case-4 b
M D z
M
B C z
K J
Case-5
A M I H
B C F G
D E
22
❑ Stress due to a uniform rectangular load dQ p
Example 2: applying the principle of superposition to calculate the stress at M y
A E
AFME+FBGM+MGCD x l
Case-4 b
M D
F z
M
B G C z
K J
Case-5
A M I H
B C F G
ABCM+KMIJ+MDEI+MCGH-MCIF
D E
23
II.1 Stresses in ground (induced by surface loading of different geometries)
q0 BL
s z =
( B + z )( L + z )
z/2 z/2
24
Example 3 A flexible rectangular area measures 2.5 m×5 m in plan. It supports
a load of 150 kN/m2.
Determine the vertical stress increase due to the load at a depth of
6.25 m below the center of the rectangular area.
(1) Using general rectangular index method;
(2) Using simplified center point index method;
(3) using 2:1 method
25
Example 3 A flexible rectangular area measures 2.5 m×5 m in plan. It supports
a load of 150 kN/m2.
Determine the vertical stress increase due to the load at a depth of
6.25 m below the center of the rectangular area.
Solution 2: s = q I m = L B , n = 2 z B
z 0 c 1 1
Solution 1: s z = q0 I m = B z; n = L z
I a = f ( m2 , n2 )
m2 = B H
n2 = L H
27
II.1 Stresses in ground (induced by surface loading of different geometries)
28
Example 4
Refer to the Figure below. Determine the average stress
increase below the center of the loaded area between z
= 3 m to z = 5 m (that is, between points A and A’).
B1 + B2
1 = tan −1 − 2
z
sz
(1 & 2 in radians) = tan −1 B1
2
z
30
Example 5
An embankment is shown in the figure.
Determine the stress increase under the
embankment at points A1 and A2.
31
Example 5
An embankment is shown in the figure.
Determine the stress increase under the
embankment at points A1 and A2.
B2 (B1=0)
32
In summary
33
Break 34
II.2 Settlement based on the theory of elasticity
Settlement can be estimated from the strain (perfectly flexible) by Bowles (1987):
1 H
( s z − s y − s x ) dz
H
Se = z dz =
0 E 0
Es = pa’N60
pa = 100 kN/m2
’= 10
Given: μs = 0.3 and qo = 167.7 kN/m2. Estimate the elastic settlement
( = s ) below the center of the foundation. Use Eqation
( = s )
Solution: B = L = 2.44 m
H = 10.98 m < 5B
40
Example 7
A rigid shallow foundation 1 m × 2 m is shown in
Figure below. Calculate the elastic settlement at
the center of the foundation.
( = s )
Es =
Es( i ) z
z
=
(10 ,000 )( 2 ) + ( 8,000 )(1) + (12 ,000 )( 2 )
= 10 ,400 kN m 2
5
For the center of the foundation,
=4
L 2
m = = =2
B 1
and
H 5
n = = = 10
B 1
2 2
42
From Table, F1=0.641 and F2=0.031. then
1 − 2s
I s = F1 + F2
1 − s
1 − 2 0.3
= 0.641 + 0.031 = 0.659
1 − 0.3
Again, Df/B = 1/1 = 1, L/B = 2, and μs = 0.3. From the table, If = 0.709.
Hence,
1 − s2
Se( flexible ) = qo ( B ) Is I f
Es
1 1 − 0.3
2
= (150 ) 4 ( 0.659 )( 0.709 ) = 0.0122 m = 12.2 mm
2 10 ,400
43
II.2 Settlement based on the theory of elasticity
I F = + 4.6 + 10 f
4
E0 + k Be 2 Be
0.8 −1
1 + 0.6 ( 0 ( e ) ) flexible
E kB
IG = −1
1.27 + 0.75 ( E ( kB ) )0.8 rigid
0 e
Alternatively for IG
−1
Ef 2t
3
I F = + 4.6 + 10
4 E0 + k Be 2 Be
45
II.2 Settlement based on the theory of elasticity
1− 2 1 − 2
Se = qo ( B) Is I f Se = qo Be I G I F I E
E E (Key parameters of soil for design)
✓ Suggestion of E E = N 60 pa
by Kulhawy and
Mayne (1990)
✓ Suggestion of E by E = cu
Duncan and
Buchignani (1976)
46
For a shallow foundation supported by a silty
Example 8 clay,
Length: L = 3.05 m
Width: B = 1.52 m
Depth of foundation: Df = 1.52 m
Thickness of foundation: t =0.305 m
Load per unit area: qo = 239.6 kN/m2
Ef = 15.87 × 106 kN/m2
The silty clay soil has the following properties:
H = 3.66 m
μs = 0.3
Solution: Eo = 9660 kN/m2
k = 565.6 kN/m2/m
Estimate the elastic settlement of the
foundation below center.
I G = 0.69
I F = 0.785
I E = 0.908
47
In summary
48
Break 49
e
II.3 Consolidation settlement Vertical strain z =
1 + e0
For homogenuous field:
❑ Recall of one-dimensional consolidation settlement Hc
Sc = z dz = z H c
0
e
Solid
Solid
Not constant!
50
II.3 Consolidation settlement
s c = s z0 log (s z )
(Compression index)
e
z =
1 + e0 s z
Cs (1 + e0 ) (ii)
(iii) From overconsolidated to normally consolidated
with s z0 s c & s z0 + s z s c : Cc (1 + e0 )
C s C s + s z
Sc ( p ) = s log c + c log z 0 H0 e
1 + e0 s
z0 1 + e0 s
c z =
1 + e0 51
II.3 Consolidation settlement
Cc = 0.1576e0 wL + 0.193
53
❑ Settlement during consolidation
Terzaghi’s theory of 1-D consolidation
Time factor:
M= (2m + 1)
2 ct
m =
Tv = v2 (unitless)
2
U =1−
m=0 M 2
exp( − M 2
Tv ) d
ue S
(Average degree of consolidation U = 1 − = c (t) )
ui Sc (max)
❖ Approximate solutions:
2
U U 0.6 Sc (t ) = USc ( p )
Tv = 4
−0.933log(1 − U ) − 0.085 U 0.6
55
II.3 Consolidation settlement
❑ Secondary consolidation settlement
At the end of primary consolidation (i.e., after the complete dissipation of excess
pore water pressure) some settlement is observed that is due to the plastic
adjustment of soil fabrics. This stage of consolidation is called secondary
consolidation.
Attention to definitions: C = e log t
C e = e ln t
C = log t
v
C t
Sc (s) = H c log 2
1 + ep t1
e2
Zhao (2017)
57
II.3 Consolidation settlement
Under
effective
stress
Creep movement !
58
II.3 Consolidation settlement Creep always exists under the effective
❖
stresses (loading), independent of the
excess pore water pressure.
❑ Secondary consolidation settlement
❖ Therefore, creep has nothing to do with
the “primary” consolidation.
❖ Creep exists during and after “primary”
consolidation.
❖ Creep rate depends on stress states:
▪ Creep rate is large in a normally
consolidated state.
▪ Creep rate is small in a over-
consolidated state.
❑ For C/Cc
❖ Inorganic clays and silts = 0.04±0.01
Mersri and Godlewski (1977)
❖ Organic clays and silts = 0.05±0.01
❖ Peats = 0.075±0.01 59
II.3 Consolidation settlement
C C
(g) l=0.05 and 0.1 (h) l=0.05 and 0.1
R2=0.895
Training Testing R2=0.857
(120 measurements) 61
Example 9
A plan of a foundation 1 m × 2 m. is shown in Figure below. Estimate the
consolidation settlement of the foundation (using the 2:1 method for s)
16.5
eo
Cc
s l =
(150 )(1)( 2 ) = 25kN m 2
(1 + 2 )( 2 + 2 )
63
Similarly
s m =
(150 )(1)( 2 )
= 13.45kN m2
(1 + 3.25)( 2 + 3.25)
and s b =
(150 )(1)( 2 ) = 8.39 kN m2
(1 + 4.5)( 2 + 4.5 )
1
Thus, = 25 + 4 (13.45 ) + 8.39 = 14.53kN m 2
s av
6
So Sc =
( 0.32 )( 2.5) log 52.84 + 14.53 = 0.0469 m = 46.90 mm
1 + 0.8 52.84
❖ “Elastic settlement” is required for both sand and clay layers. For clay, it refers to
immediate settlement without sense of time.
❖ “Primary settlement” and “Secondary settlement” are only required for clay layer.
65
Break 66
II.4 Field load test
Plate load test
in the field
67
II.4 Field load test
68
II.4 Field load test
Definition of parameters
for differential settlement
(tilt)
Polshin and Tokar (1957) suggested the following allowable deflection ratios
for buildings as a function of the ratio of the length to the height of a building:
71
II.5 Settlement restriction
72
Laurits Bjerrum (1918-1973)
II.5 Settlement restriction
73
In summary
BF
Field load test SF = SP
qu ( F ) = qu ( P ) BP (size effect)
❖ For tests in clayey soils
❖ For tests in sandy soils BF 2
qu ( F ) = qu ( P ) 2B F
BP SF = SP
F
B + B P
(tilt)
1
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
Q2 L3
X=
Q1 + Q2
B=
A L = 2 ( L2 + X )
L
(L = L1+L2+L3)
3
III.1 From combined footings to Mat Foundation
4
III.1 From combined footings to Mat Foundation
❑ Cantilever footing
To strengthen the integrity of the foundation, adjust the uneven settlement between the
foundations, and eliminate or reduce the sensitivity of the superstructure to settlement.
5
III.1 From combined footings to Mat Foundation
is a combined footing that may cover the entire area under a structure
supporting several columns and walls.
is a combined footing that may cover the entire area under a structure
supporting several columns and walls.
8
Attention:
III.2 Bearing capacity of mat foundations In footing: Qall ( net ) = Qlive + dead
In mat foundation: Qall ( net ) = Qlive + dead − D f
Gross ultimate bearing capacity pressure
qu = cN c Fcs Fcd Fci + qN q Fqs Fqd Fqi + 12 BN F s F d F i
Special case:
for saturated clays (f’=0) and a vertical loading condition, qu(net)=?
N q = 1; N c = 5.14; N = 0
B Nq B 1 0.195B
f = 0 Fcs = 1 + = 1+ = 1+
L Nc L 5.14 L
D
Fcd = 1 + 0.4 f
B
1
B
12
Example 2
Following the results of a standard penetration test in Depth (m) Field value of N60
the field (sandy soil) as the table at right:
1.5 9
Estimate the net allowable bearing capacity of a mat
foundation 6.5 m × 5 m in plan. Here, Df = 1.5 m and 3.0 12
allowable settlement = 50 mm. Assume that the unit 4.5 11
weight of soil, γ = 16.5 kN/m3.
6.0 7
7.5 13
N 60 D f Se Se 9.0 11
qall( net ) = 1+0.33 16.63 N 60
0.08 B 25 25 10.5 13
The bottom of the foundation = 1.5 m,
2B below the bottom = 1.5+2x5=11.5 m
Solution:
9 + 12 + 11 + 7 + 13 + 11 + 13
N 60 = 11
7
N D S S
qall( net ) = 60 1+0.33 f e 16.63 N 60 e
0.08 B 25 25
11 1.5 50
= 1 + 0 .3 3 = 302.23kN m 2
365.86 kN m 2
0.08 5 25
13
III.3 Settlement of mat foundations
qu ( net ) L B
FS = =
Q
Depth of mat for a fully compensated foundation qall ( net )
A
− Df
Q Q
qall ( net ) = − D f =0 D f =
A A
Also, stress increase is zero since q=0,
no settlement. So, absolutely safe!
qu ( net ) qu ( net ) Q
FS = = f
D
qall ( net ) Q
− Df A
A
15
Example 3
A mat foundation on a saturated clay soil has dimensions of 20 m × 20 m. Given: dead
and live load =48 MN, cu = 30 kN/m2, and γclay = 18.5 kN/m3.
a. Find the depth, Df, of the mat for a fully compensated foundation.
b. What will be the depth of the mat (Df) for a factor of safety of 3 against bearing
capacity failure?
Q Q
Compensated: qall ( net ) = − D f =0 D f =
A A
0.195 B Df
5.14cu 1 + +
B
1 0 .4
L
FS =
Q
− Df
A
16
Example 3
A mat foundation on a saturated clay soil has dimensions of 20 m × 20 m. Given: dead
and live load =48 MN, cu = 30 kN/m2, and γclay = 18.5 kN/m3.
a. Find the depth, Df, of the mat for a fully compensated foundation.
b. What will be the depth of the mat (Df) for a factor of safety of 3 against bearing
capacity failure?
19
Example 4
Solution:
20
Summary
cu = 120 kN m 2 , f = 0, B = 8 m, L = 18 m, D f = 3 m
0.195 B Df
qu( net ) = 5.14cu 1 + 1 + 0.4
L B
22
Solution for exercise 1
Determine the net ultimate bearing capacity of mat foundations with the following
characteristics:
cu = 120 kN m 2 , f = 0, B = 8 m, L = 18 m, D f = 3 m
Solution:
0.195 B Df
qu ( net ) = 5.14cu 1 + 1 + 0.4
L B
0.195 8 3
= 5.14 120 1+
1+0.4
18 8
= 770.8 kN m 2
23
Exercise 2
Following are the results of a standard penetration test in the field (sandy soil):
Estimate the net allowable bearing capacity of a mat foundation 6.5 m × 5 m in plan. Here,
Df = 1.5 m and allowable settlement = 30 mm. Assume that the unit weight of soil, γ =
16.5 kN/m3.
N 60 D f Se Se
qall( net ) = 1+0.33 16.63 N 60
0.08 B 25 25
24
Solution for exercise 2
Estimate the net allowable bearing capacity of a mat foundation 6.5 m × 5 m in plan. Here,
Df = 1.5 m and allowable settlement = 30 mm. Assume that the unit weight of soil, γ =
16.5 kN/m3.
Solution:
The bottom of the foundation = 1.5 m,
2B below the bottom = 1.5+2x5=11.5 m
9 + 12 + 11 + 7 + 13 + 11 + 13
N 60 = = 10.86
7
N 60 D f Se Se
qall( net ) = 1+0.33 16.63 N 60
0.08 B 25 25
10.86 1.5 30
= 1 + 0 .3 3 = 179 k N m 2
216.7 kN m 2
0.08
5 25
25
Exercise 3
A mat foundation is shown in Figure. The design considerations are L = 12 m, B = 10 m,
Df = 2.2 m, Q = 30 MN, x1 = 2 m, x2 = 2 m, x3 = 5.2 m, and preconsolidation pressure σ’c
≈ 105 kN/m2. Calculate the consolidation settlement under the center of the mat.
Q Sc(p) = 378 mm
q= − Df
A
H 2 I a( H 2 ) − H1I a( H1 )
(H
av H1 )
= qo
2
H 2 − H 1
CH + av
Sc( p ) = c c log o
( For NC Clay )
1 + eo o
26
Solution for exercise 3
A mat foundation is shown in Figure. The design considerations are L = 12 m, B = 10 m,
Df = 2.2 m, Q = 30 MN, x1 = 2 m, x2 = 2 m, x3 = 5.2 m, and preconsolidation pressure σ’c
≈ 105 kN/m2. Calculate the primary consolidation settlement under the center of the mat.
Solution: o = ( D f + x1 ) (16 ) + x2 (18 − 9.81) + x3 (17.5 − 9.81)
2
5.2
= ( 4.2 )(16 ) + ( 2 )( 8.19 ) + ( 7.69 ) = 103.574 kN m
2
2
o c Normally consolidated clay.
Cc H c + av
Sc( p ) = log o
1 + eo
o
H c = x3 = 5.2 m; Cc = 0.38; eo = 0.88.
Q 30 103
q = − Df = − 16 2.2 = 214.8 kN m 2
A 12 10
H 2 I a( H 2 ) − H1I a( H1 ) ( 2 + 2 + 5.2 ) I a( H 2 ) − ( 2 + 2 ) I a( H1 )
(H
av H1 )
= qo = 214.8
2
H 2 − H 1 5.2 27
Solution for exercise 3
A mat foundation is shown in Figure. The design considerations are L = 12 m, B = 10 m,
Df = 2.2 m, Q = 30 MN, x1 = 2 m, x2 = 2 m, x3 = 5.2 m, and preconsolidation pressure σ’c
≈ 105 kN/m2. Calculate the primary consolidation settlement under the center of the mat.
For I a( H 2 ) ,
B2 10 2
m2 = = = 0.54
H 2 2 + 2 + 5.2
I a( H 2 ) =0.19
L2 12 2
n2 = = = 0.65
H 2 2 + 2 + 5.2
Again, for I a( H1 )
B 2 10 2
m2 = = = 1.25
H1 2 + 2
I a( H1 ) = 0.235
L 2 12 2
n2 = = = 1.5
H1 2 + 2
( 2 + 2 + 5.2 )( 0.19 ) − ( 2 + 2 )( 0.235 )
( H H ) = 214.8
av = 33.38 kN m 2
5.2
2 1
So, the stress increase below the center of the 12 m 10 m area is 4 33.38 = 133.52 kN m 2 .
0.38 5200 103.574 + 133.52
Sc( p ) = log = 378 mm
1 + 0.88 103 .574 28
Exercise 4
For the mat foundation in Exercise 3, estimate the primary consolidation settlement under
the corner of the mat.
Sc(p) = 168 mm
Q
q= − Df
A
H 2 I a( H 2 ) − H1I a( H1 )
(H
av H1 )
= qo
2
H 2 − H 1
CH + av
Sc( p ) = c c log o
( For NC Clay )
1 + eo o
29
Solution for exercise 4
For the mat foundation in Exercise 3, estimate the primary consolidation settlement under
the corner of the mat.
5.2
= ( 4.2 )(16 ) + ( 2 )( 8.19 ) + ( 7.69 ) = 103.574kN m
2
2
o c Normally consolidated clay.
Cc H c + av
Sc( p ) = log o
1 + eo
o
H c = x3 = 5.2 m; Cc = 0.38; eo = 0.88.
Q 30 103
q = − Df = − 16 2.2 = 214.8kN m 2
A 12 10
H 2 I a( H 2 ) − H1I a( H1 ) ( 2 + 2 + 5.2 ) I a( H 2 ) − ( 2 + 2 ) I a( H1 )
(H
av H1 )
= qo = 214.8
2
H 2 − H 1 5.2
30
Solution for exercise 4
For the mat foundation in Exercise 3, estimate the primary consolidation settlement under
the corner of the mat.
For I a( H 2 ) ,
B 10
m2 = = = 1.09
H 2 2 + 2 + 5.2
I a( H 2 ) =0.23
L 12
n2 = = = 1.3
H 2 2 + 2 + 5.2
Again, for I a( H1 )
B 10
m2 = = = 2.5
H1 4
I a( H1 ) = 0.25
L 12
n2 = = =3
H1 4
( 2 + 2 + 5.2 )( 0.23) − ( 2 + 2 )( 0.25 )
av ( H H ) = 214.8 = 46 .1 kN m 2
5.2
2 1
Se =
qo Be I G I F I E
Eo
(1 − s2 ) with Be =
4 BL
1
IF = + 3
4 Ef 2t
4.6 + 10
Eo + Be k 2 Be
1
IE = 1−
B
3.5 exp (1.22 s − 0.4 ) e + 1.6
Df 32
Solution for Exercise 5
4 BL 4 3 3
Be = = = 3.385 m
Eo 16000
= = = 11.82
kBe 400 3.385
H 20
= = 5.91
Be 3.385
1
IF = +
4 3
Ef 2t
4.6 + 10
Be Be
Eo + k
2
1
= + = 0.815
4
15 10 6
2 0.25 3
4.6 + 10
16000 +
3.385
400 3.385
2
33
Solution for Exercise 5
1
IE = 1−
B
3.5 exp (1.22 s − 0.4 ) e + 1.6
Df
1
=1 − = 0.923
3.385
3.5 exp (1.22 0.3 − 0.4 ) + 1.6
1.5
Se =
qo Be I G I F I E
Eo
(1 − s2 )
34
TOPICS & SYLLABUS:
Topic 1: Site Investigation
Topic 2: Slope Stability
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
1
Topic 5: Earth Retaining Structures
Key points of a retaining structures design:
(a) Lateral earth pressure theories;
(b) Factor of safety in terms of overturning, sliding and bearing capacity.
I. Retaining Walls II. Sheet Pile Walls
❖ Lateral earth pressure theories ❖ Cantilever sheet pile walls
❖ Stability of retaining walls ❖ Anchored sheet-pile walls
❖ Gravity and cantilever walls ❖ Penetration into soils
❖ Mechanically stabilized retaining III. Braced Cuts
walls ❖ Pressure envelope
❖ Design of components
❖ Stability and deformation
2
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
3
I.1 Types of retaining wall
✓ Economical to a
height of about 8 m
5
I.3 Application of lateral earth pressure theories to design
K p = tan 2 ( 4 + 2 )
p = 0 K p + 2c K p K 0 (1 − sin ) OCR sin
K a = tan 2 ( 4 − 2 )
a = 0 K a − 2c K a
6
I.3 Application of lateral earth pressure theories to design
K p = tan 2 ( 4 + 2 )
p = 0 K p + 2c K p
7
I.3 Application of lateral earth pressure theories to design
❑ Rankine’s active earth pressure
Cohesionless soils
8
I.3 Application of lateral earth pressure theories to design
❑ Rankine’s passive earth pressure
Kp=?
Cohesionless soils
9
I.3 Application of lateral earth pressure theories to design
❑ Rankine’s active earth pressure
Ka=?
Cohesive soils
10
I.3 Application of lateral earth pressure theories to design
❑ Rankine’s active earth pressure
Ka=?
11
I.3 Application of lateral earth pressure theories to design
❑ Rankine’s passive earth pressure
Kp=?
Cohesive soils
12
I.3 Application of lateral earth pressure theories to design
❑ Rankine’s passive earth pressure
Kp=?
13
I.3 Application of lateral earth pressure theories to design
15
I.3 Application of lateral earth pressure theories to design
(active)
16
I.3 Application of lateral earth pressure theories to design
❑ Coulomb’s earth pressure theory (active)
17
I.3 Application of lateral earth pressure theories to design
❑ Coulomb’s earth pressure theory (active)
18
I.3 Application of lateral earth pressure theories to design
❑ Coulomb’s earth pressure theory (active)
For a special case: = ; = 0; = 0; = + KA=?
2 4 2
sin 2 ( + )
KA = 2
sin ( + ) sin ( − )
sin 2 ( ) sin ( − ) 1 +
B sin ( − ) sin ( + )
C
sin 2 +
= 2
2
pa
sin ( ) sin ( )
2
sin sin 1 +
2
2 sin sin
2 2
A
sin 2 +
( ) = tan 2 −
2
2 = cos
=
1 + sin ( ) 1 + sin ( ) 4 2
2
Same as Rankine’s Ka
19
I.3 Application of lateral earth pressure theories to design
(passive)
20
I.3 Application of lateral earth pressure theories to design
❑ Coulomb’s earth pressure theory (passive)
pp R
21
I.3 Application of lateral earth pressure theories to design
❑ Coulomb’s earth pressure theory (passive)
22
I.3 Application of lateral earth pressure theories to design
❑ Coulomb’s earth pressure theory (passive)
For a special case: = ; = 0; = 0; = − Kp=?
2 4 2
sin 2 ( − )
KP = 2
sin ( + ) sin ( + )
sin ( ) sin ( + ) 1 −
2
sin ( + ) sin ( + )
sin 2 −
= 2
2
pp
sin ( ) sin ( )
2
sin sin 1 −
2
45 − 2 2 sin sin
2 2
sin 2 −
( )
2
2 = cos 2
= = tan +
1 − sin ( )
2
1 − sin ( ) 4 2
Same as Rankine’s Kp
23
In summary
24
Break
25
I.4 Stability of retaining walls
❑ Failure modes
Sliding
Overturning
bearing capacity
failure
deep-seated
shear failure
(due to weak
soil layers)
26
I.4 Stability of retaining walls
(Attention: inclined)
Overturning point
Attention:
when estimating overturning,
suppose no supporting pressure.
27
I.4 Stability of retaining walls
Attention:
For more safe design, Mpp is not
accounted in resisting moment.
28
I.4 Stability of retaining walls
❑ Check for sliding along the base (design requirement normally FS > 1.5)
29
I.4 Stability of retaining walls
(2)
(1)
30
I.4 Stability of retaining walls
V 6e
qmax = 1 +
B B
Attention:
For more safe design, e = B − CE = B − M net
Pp is not accounted for 2 2 V
inclination.
31
Example 1
The cross section of a cantilever retaining wall is shown in the figure below.
Calculate the factors of safety with respect to overturning, sliding, and
bearing capacity.
32
Solution:
33
Solution:
34
35
Solution:
36
Example 2
A gravity retaining wall is shown in the following figure. Use ’ = 2/3’1 and
Coulomb’s active earth pressure theory. Determine:
(a) The factor of safety against overturning.
(b) The factor of safety against sliding.
(c) The pressure on the soil at the toe and heel.
37
38
Solution:
39
Solution:
Attention:
❖ Using Coulomb’s
theory, no soil above
the back face;
❖ Using Rankine’s theory,
account the soil above.
40
Solution:
41
Solution:
42
I.5 Additional consideration of design for construction
A retaining wall may be constructed with one or more of the following joints:
(1) Construction joints are vertical and horizontal joints that are placed between two successive
pours of concrete.
(2) Contraction joints are vertical joints (grooves) placed in the face of a wall (from the top of
the base slab to the top of the wall) that allow the concrete to shrink without noticeable harm.
(3) Expansion joints allow for the expansion of concrete caused by temperature changes; These
joints may be filled with flexible joint fillers.
43
I.5 Additional consideration of design for construction
44
In summary
45
Break
46
I.5 Mechanically stabilized retaining walls - Soil reinforcement
Geogrids
Metal Strips
47
I.5 Mechanically stabilized retaining walls - General design considerations
48
I.6 Retaining walls with metallic strip reinforcement
❑ Main components
✓ 1. Backfill, which is granular soil
✓ 2. Reinforcing strips, which are thin, wide strips placed at regular intervals
✓ 3. A cover or skin, on the front face of the wall
49
I.6 Retaining walls with metallic strip reinforcement
2q
a (2) = M ( − sin cos 2 )
a (1) = 0(1)
K a − 2c K a ⎯⎯⎯⎯ c '= 0
→ a (1) = 1 zK a
= z
0 (1) 1
0.4b
K a = tan 2 ( 4 − 2 ) M = 1.4 − 1
0.14 H
50
I.6 Retaining walls with metallic strip reinforcement
L = lr + le
lr =
( H − z)
tan ( 45 + 1 2 )
L=
( H − z)
+
FS(P) a SV S H The dimension of L should be
based on the depth of the layer
tan ( 4 + 1 2 ) 2 w 0 tan u next to the top of backfill.
52
I.6 Retaining walls with metallic strip reinforcement
(trial-error)
(depth-dependent)
53
I.6 Retaining walls with metallic strip reinforcement
54
I.6 Retaining walls with metallic strip reinforcement
56
Solution:
- Check for internal stability
58
I.7 Retaining walls with geotextile reinforcement
59
I.7 Retaining walls with geotextile reinforcement
60
I.7 Retaining walls with geotextile reinforcement
62
Example 4
A geotextile-reinforced retaining wall 5 m high is shown in the figure below. For the granular
backfill, 1 = 15.7 kN/m3, ’1 = 36o. For the geotextile, Tult = 52.5 kN/m.
For the design of the wall, determine SV , L and ll. (use RFid = 1.2, RFcr = 2.5, and RFcbd = 1.25)
For the retaining wall, calculate the factor of safety against overturning, sliding, and bearing
capacity failure.
63
Solution:
- Check for internal stability
(z = 4.5 m should be
considered or recheck)
65
I.8 Retaining walls with geogrid reinforcement
66
I.8 Retaining walls with geogrid reinforcement
67
I.8 Retaining walls with geogrid reinforcement
68
Example 5
Consider a geogrid-reinforced retaining wall. Referring to the figure below, given: H = 6 m,
1 = 16.5 kN/m3, ’1 = 35o, Tall = 45 kN/m, FS(B) = 1.5, FS(P) = 1.5, Cr = 0.8, and Ci = 0.75. For the
design of the wall, determine SV and L.
70
TOPICS & SYLLABUS:
Topic 1: Site Investigation
Topic 2: Slope Stability
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
1
Topic 5: Earth Retaining Structures
Key points of a retaining structures design:
(a) Lateral earth pressure theories;
(b) Factor of safety in terms of overturning, sliding and bearing capacity.
I. Retaining Walls II. Sheet Pile Walls
❖ Lateral earth pressure theories ❖ Cantilever sheet pile walls
❖ Stability of retaining walls ❖ Anchored sheet-pile walls
❖ Gravity and cantilever walls ❖ Penetration into soils
❖ Mechanically stabilized retaining III. Braced Cuts
walls ❖ Pressure envelope
❖ Design of components
❖ Stability and deformation
2
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
3
II.1 Types of sheet pile
❑ Wooden sheet piles: only for temporary, light
structures that are above the water table
(cross section 50 mm x 300 mm)
5
II.3 Cantilever sheet piling penetrating sandy soils (c’ = 0, f ’ ≠ 0)
P
1 1
No fills FH = 0 p D − a D − P = 0
Rankine's theory: 2 2
p = 0 K p + 2c K p 2P
L D2 K p − D2 Ka = 2P D =
( K p − Ka )
a = 0 K a − 2c K a ?
1 D 1 D
M = 0 p D − a D − P ( L + D ) = 0
2 3 2 3
D=? z (for Mmax) D 3 K p − D 3 K a − 6 PD − 6 PL = 0
Passive Active = ( K p − K a ) D 3 − 6 PD − 6 PL = 0 D = ?
DK p DK a
p − a = D ( K p − K a )
O [Note: More simple if only one side of
pressure by “passive-active”]
M = P ( L + z ) + a z − p z = P ( L + z ) + ( K a − K p ) z 3
1 z 1 z 1
2 3 2 3 6
Strong assumption
(free earth support) M = 0 P + 1 ( K a − K p ) z 2 = 0 z = 2P
M max = ?
z 2 ( K p − Ka )
Note: 2 equations,
Or zero shear stress: P − ( K p − K a ) z = 0 z =
z 2P
( K p − Ka )
but 1 unknown (D). 2 6
II.3 Cantilever sheet piling penetrating sandy soils (c’ = 0, f ’ ≠ 0)
P P ❖ The distribution is only
No fills due to Rankine’s earth
pressures
L L ❖ In one single soil layer,
the pressure distribution
should be continuous
L5=?
D ( K p − Ka ) D ( K p − Ka )
Note:
Linear distribution is idealized, the
reality should be nonlinear curved.
L4 L4 L5 L5 L L3
M = 0 − ( + ) − P1 D + −
2 P D − =0
L5 = ?
3 3 4
2 3 2 3 3 3
L L
F = 0 − ( + ) − P = 0 L =
4 5
3 L4 − 2 P D=?
H 3
2
3 4
2
5
( 3 + 4 )
P − ( K p − Ka ) z
z 2P
=0 z =
( K p − Ka )
By zero shear stress: 2
9
II.3 Cantilever sheet piling penetrating sandy soils (c’ = 0, f ’ ≠ 0)
D = L3 + L4
4 = 5 + L4 ( K p − K a )
⎯⎯⎯⎯
→
with 5 = L3 ( K p − K a ) + ( L1 + L2 ) K p
L1 L L
P = 1 + ( 1 + 2 ) 2 + 2 3
2 2 2
2 = ( L1 + L2 ) K a
L3 = 2 ( K p − K a )
L1 L1 L2
1 2 3 + L2 + L3 + 1L2 2 + L3
1
z=
P L2 L2 L3 2 L3
+ ( 2 − 1 ) + L3 + 2 L44 + A1 L34 − A2 L24 − A3 L4 − A4 = 0
2 3 2 3
5 = L3 ( K p − K a ) + ( L1 + L2 ) K p
L4 = ...
D = L3 + L4 12
II.3 Cantilever sheet piling penetrating sandy soils (c’ = 0, f ’ ≠ 0)
E
D=?
O L4 = D − L3 z (for Mmax)
L5=?
3 4
14
Solution
Part a
Using the Figure of the pressure
distribution diagram, one can now
prepare the following table for a
step-by-step calculation.
15
Solution
Part b
The total length of the sheet pile is
Part c
Finally, we have the following table:
16
II.4 Cantilever sheet piling penetrating clayey soils (c ≠ 0, f = 0)
P Example 1: Please demonstrate the whole
No fills procedure with formulas for
(a) the theoretical depth of penetration of the
L
clay layer by the sheet pile D;
(b) the position and the value of maximum
bending moment of the sheet pile.
Rankine's theory: K a = tan 2 ( 4 − f 2 ) = 1
p = 0 K p + 2c K p
K p = tan ( 4 + f 2 ) = 1
2
D=?
a = 0 K a − 2c K a
O
Strong assumption
17
II.4 Cantilever sheet piling penetrating clayey soils (c ≠ 0, f = 0)
P Example 1: Please demonstrate the whole
No fills procedure with formulas for
Rankine's theory:
p = 0 K p + 2c K p
(a) the theoretical depth of penetration of the
L
clay layer by the sheet pile D;
a = 0 K a − 2c K a (b) the position and the value of maximum
bending moment of the sheet pile.
2c
K a = tan 2 ( 4 − f 2 ) = 1
z (for Mmax)
K p = tan ( 4 + f 2 ) = 1
2
19
II.4 Cantilever sheet piling penetrating clayey soils (c ≠ 0, f = 0)
P P K a = tan 2 ( 4 − f 2 ) = 1
No fills K p = tan ( 4 + f 2 ) = 1
2
Rankine's theory: L L
Solution:
p = 0 K p + 2c K p
a = 0 K a − 2c K a
p − a = 4c p − a = 4c
More general assumption
D L4 L4 P ( P + 12cL )
M = 0 4cD − 8c − P ( L + D ) = 0
4cD 2
− 2 PD − =0 D =?
2 2 3 2c
L
F = 0 4cD − 8c 4 − P = 0 L = 4 cD − P → L = 4cD − P
H 2
4
4c 4
4c
2
P z P P
By zero shear stress: P − 4cz = 0 z = M max = P ( L + z ) − 4cz = P L + − 2c
4c 2 4c 4c
20
II.4 Cantilever sheet piling penetrating clayey soils (c ≠ 0, f = 0)
Example 3: Please demonstrate the
Granular fills
whole procedure with formulas for
(, f, c = 0)
L (a) the theoretical depth of penetration
P1 = 2
L 1 2
= L Ka
of the clay layer by the sheet pile D;
2 2 (b) the position and the value of
LK a = 2 maximum bending moment of the
sheet pile.
z (for Mmax) Rankine's theory:
D=?
p = 0 K p + 2c K p
L4=? O
a = 0 K a − 2c K a
4c 4c
For saturated clay, f = 0
K a = tan 2 ( 4 − f 2 ) = 1
K p = tan ( 4 + f 2 ) = 1
2
21
II.4 Cantilever sheet piling penetrating clayey soils (c ≠ 0, f = 0)
23
II.4 Cantilever sheet piling penetrating clayey soils (c ≠ 0, f = 0)
E
z ( L1 + L2 ) K a = 2 FH = 0
1
z (for Mmax) P −6D + ( 6 + 7 ) L4 = 0
D=? 2
L4=? O 4c − ( L1 + L2 ) D − P
L4 =
4c
6 = 4c − ( L1 + L2 ) 4c+ ( L1 + L2 ) = 7 M = 0
( p −a ) ( p −a )
6D
D 1 L
− M P − ( 6 + 7 ) L4 4 = 0
2 2 3
M at the point E : Pz = M E z = M E P
P ( P + 12cz )
L1 L1 L2 L2 L2 ( 4c − L1 − L2 ) D 2 − 2 PD − =0
with M E = 1 + L2 + 1L2 + ( 2 − 1 ) L1 + L2 + 2c
23 2 2 3
D=?
By zero shear stress:
P 6 M max
P −6z = 0 z = M max = p ( z + z ) − z2 S=
6 2 all
L1 L1 L2 L2 L2 6 2
= 1 + L2 + z + 1L2 + z + ( 2 − 1 ) + z − z
23 2 2 3 2 24
Tutorials
Exercise 2
In the figure below, for the sheet pile wall penetrating into saturated clay, determine
(a) The theoretical and actual depth of penetration, using Dactual = 1.5Dtheory.
(b) The minimum size of sheet pile section necessary, using all = 172.5 MN/m2.
25
Solution
26
Solution
27
II.5 Anchored sheet-pile walls
❑ When the height of the
backfill material behind a
cantilever sheet-pile wall
exceeds about 6 m, tying the
wall near the top to anchor
plates, anchor walls, or
anchor piles becomes more
economical.
D free D fixed free earth support method
❑ Anchors minimize the depth
of penetration required by the
sheet piles and also reduce
the cross-sectional area and
weight of the sheet piles
needed for construction.
❑ Two design methods
depending on assumptions
(free earth support method;
fixed earth support method)
fixed earth support method
28
II.6 Free Earth Support Method for Penetration of Sandy Soil (c’ = 0, f ’ ≠ 0)
L1 K a = 1
Active
L1 L L
Passive P = 1 + ( 1 + 2 ) 2 + 2 3
2 2 2
ME
M E = Pz z =
8 = p − a = L4 ( K p − K a ) P
L1 L1 L2 L2 L2 L3 2 L3
M E = 1 + L2 + L3 + 1L2 + L3 + ( 2 − 1 ) + L3 + 2
23 2 2 3 2 3
29
II.6 Free Earth Support Method for Penetration of Sandy Soil (c’ = 0, f ’ ≠ 0)
Note:
❖ Without anchor, we have 2 equations
(sum(M)=0, sum(FH)=0) but 1 unknown (D).
❖ With anchor, we have 2 unknowns (D, Fanchor).
L4 2 L4
P ( L1 + L2 + L3 ) − ( z + l1 ) − 8 l 2 + L2 + L3 + =0
2 3
L4 = ? (for length of sheet pile)
F H =0 (for anchor)
L4 L
P − 8 − F = 0 F = P − 8 4 = ?
2 2
( z − L1 ) = 0 z = ?
2
L
1 1 − F + 1 ( z − L1 ) + K a
2 2 (Zero shear stress/maximum
L L x 2
x x moment occurs within L )
M max = 1 1 1 + x − F ( l2 + x ) + 1x + K a 2
23 2 2 3 30
Exercise 3
Let L1 = 3.05 m, L2 = 6.1 m, l1 = 1.53 m, l2 = 1.52 m, c’ = 0, f’ = 30o, = 16 kN/m3, sat = 19.5
kN/m3, and E = 207,103 MN/m2 in the figure below.
(a) Determine the theoretical and actual depths of penetration. (Note: Dactual = 1.3Dtheory)
(b) Find the anchor force per unit length of the wall.
(c) Determine the maximum moment, Mmax.
31
Part a
Solution We use the following table.
32
Part b
Solution The anchor force per unit length of the wall is
Part c
33
II.6 Free Earth Support Method for Penetration of clayey Soil (c’ ≠ 0, f ’= 0)
L1 K a = 1
Active
L1 L
Passive P = 1 + ( 1 + 2 ) 2
2 2
P
MD
M D = Pz z =
P
6 = p − a = 4c − ( L1 + L2 )
L1 L1 L2 L2 L2
M D = 1 + L2 +
1 L2 + (
2 −
1 )
23 2 2 3 34
II.6 Free Earth Support Method for Penetration of clayey Soil (c’ ≠ 0, f ’= 0)
35
II.6 Free Earth Support Method for Penetration of clayey Soil (c’ ≠ 0, f ’= 0)
Solution:
x = z − L1
M O =0
D
P ( L1 + L2 ) − ( z + l1 ) − 6 D l2 + L2 + = 0
2
P
D=?
F H =0
P − 6D − F = 0 F = P − 6D = ?
( z − L1 ) = 0 z = ?
2
L
1 1 − F + 1 ( z − L1 ) + K a
2 2
(Zero shear stress/maximum
L1 L1 x x2 x
M max = 1 + x − F ( l2 + x ) + 1x + K a moment occurs within L2)
23 2 2 3 36
Exercise 4
In the figure below, let L1 = 3 m, L2 = 6 m, and l1 = 1.5 m. Also let, = 17 kN/m3, sat = 20 kN/m3,
f’ = 35o, and c = 41 kN/m2.
(a) Determine the theoretical depth of embedment of the sheet-pile wall.
(b) Calculate the anchor force per unit length of the wall.
37
Solution
Part a
6 = 4c − ( L1 + L2 )
1.53
Part b
38
Design of anchors will be studied in
II.7 Anchors “Advanced Geotechnical Design”
Various types of anchoring for sheet-pile walls
Tieback
39
Summary
f f ME
K a = tan 2 45 − ; K p = tan 2 45 + ; z =
2 2 P
5 = ( L1 + L2 ) K p + L3 ( K p − K a )
5 8P
A1 = ; A2 =
( K p − Ka ) ( K p − Ka )
6 P 2 z ( K p − K a ) + 5 P ( 6z 5 + 4P )
A3 = ; A4 =
2 ( K p − K a ) 2 ( K p − K a )
2 2
Dtheory = L3 + L4
2P
z =
( K p − Ka )
1 z
M max = P ( z + z ) − z2 ( K p − K a )
2 3
Solution for Exercise 1
Part a:
f 32
K a = tan 2 45 − = tan 2 45 − = 0.307
2 2
f 32
K p = tan 2 45 + = tan 2 45 + = 3.255
2 2
2 40.38
L3 = = = 1.633 m
( K p − Ka ) (18.2 − 9.81) ( 3.255 − 0.307 )
1 1 1
P = 1L1 + 1L2 + ( 2 − 1 ) L2 + 2 L3
2 2 2
1 1 1
= 19.77 4+19.77 8+ ( 40.38 − 19.77 ) 8 + 40.38 1.633
2 2 2
= 39.54 + 158.16 + 82.44 + 32.97 = 313.11kN m
= 16.1 4 + (18.2 − 9.81) 8 3.255 + (18.2 − 9.81) 1.633 ( 3.255 − 0.307 )
= 2283.16
Thus
Dtheory = L3 + L4 = 1.633 + 11.683 = 13.316 m
Part b:
The total length of the sheet piles is
Part c:
2P 2 313.11
z = = =5.03 m
( K p − Ka ) (18.2 − 9.81) ( 3.255 − 0.307 )
1 z
M max = P ( z + z ) − z2 ( K p − K a )
2 3
1 5.03
= 313.11 ( 5.48 + 5.03) − (18.2 − 9.81) 5.032 ( 3.255 − 0.307 )
2 3
= 3290.79 − 524.62 = 2766.17 kN m m
Tutorial Refer to Figure below, for which L1 = 2.4 m, L2 = 4.6 m, = 15.7 kN/m3, sat = 17.3
Exercise 2 kN/m3 and f’= 30°for sand, and c = 29 kN/m2 for clay.
a. What is the theoretical depth of embedment, D?
b. Increase D by 40%. What length of sheet piles is needed?
c. Determine the theoretical maximum moment in the sheet pile.
f ME
K a = tan 2 45 − ; z =
2 P1
P1 ( P1 + 12cz1 )
D 2 4c − ( L1 + L2 ) − 2 DP1 − =0
( L1 + L2 ) + 2c
P1
z =
6
D 4c − ( L1 + L2 ) − P1
L4 =
4c
6 z 2
M max = P1 ( z + z1 ) −
2
Solution for Exercise 2
Part a:
f 30
K a = tan 2 45 − = tan 2 45 − = 0.333
2 2
1 1
P1 = 1L1 + 1L2 + ( 2 − 1 ) L2
2 2
1 1
= 12.56 2.4+12.56 4.6+ ( 24.04 − 12.56 ) 4.6
2 2
= 15.072 + 57.776 + 26.415 = 99.26 kN m
= 2.567 m
P1 ( P1 + 12cz1 )
D 2 4c − ( L1 + L2 ) − 2 DP1 − =0
( L1 + L2 ) + 2c
D 2 4 29 − 15.7 2.4 + (17.3 − 9.81) 4.6 − 2 D 99.26
Part c:
P1 99.26
z = = = 2.263 m
6 43.866
6 z 2
M max = P1 ( z + z1 ) −
2
43.866 2.2632
= 99.26 ( 2.567 + 2.263) − = 367.085kN m m
2
TOPICS & SYLLABUS:
Topic 1: Site Investigation
Topic 2: Slope Stability
Topic 3: Shallow Foundations
Topic 4: Mat (Raft) Foundations
Topic 5: Earth Retaining Structures
Topic 6: Pile Foundations
1
Question:
Skyscrapers are one of the
favorite topics of the public on
civil engineering projects.
On what are the skyscrapers
founded?
2
Ranking: Height to pinnacle (highest point)
https://en.wikipedia.org/wiki/List_of_tallest_buildings 3
The Burj Khalifa set several world records:
Jan.2004 – Jan.2010
❖ Tallest existing structure: 829.8 m
❖ Tallest structure ever built: 829.8 m
❖ Tallest freestanding structure: 829.8 m
❖ Tallest skyscraper (to top of spire): 828 m
❖ Tallest skyscraper to top of antenna: 829.8 m
❖ Building with most floors: 163
❖ World's highest elevator installation
❖ World's longest travel distance elevators: 504 m
❖ Highest vertical concrete pumping for a building:
606 m
❖ World's tallest structure that includes residential
space
❖ World's highest installation of an aluminium and
glass façade: 512 m
❖ World's highest nightclub: 144th floor
❖ World's highest restaurant (Atmosphere): 122nd
floor at 442 m
❖ World's highest New Year display of fireworks
❖ World's largest light and sound show staged on a
single building
https://en.wikipedia.org/wiki/Burj_Khalifa 4
Boreholes for site investigation
In 2003, 33 test holes were drilled to study the strength of
the bedrock underlying the structure. "Weak to very weak
sandstone and siltstone" was found, just meters below the
surface. Samples were taken from test holes drilled to a
depth of 140 meters, finding weak to very weak rock all the
way.
6
Topic 6: Pile Foundations
Key points of foundation design: (Pile: bottom, side; Shallow: bottom)
(a) shall be safe against overall shear failure (Ultimate Limit State);
(b) cannot undergo excessive displacement (Serviceability Limit State).
I. Single pile II. Pile group
❖ General information ❖ Configuration & Group Efficiency
❖ Pile loading capacity ❖ Ultimate Capacity of Group Piles
❖ Point load capacity ❖ Settlement of Group Piles
❖ Frictional resistance
III. Bored pile
❖ Negative Skin Friction
❖ Pile driving ❖ Construction Procedures
❖ Pile load test ❖ Load-Bearing Capacity
❖ Settlement of piles ❖ Elastic settlement
7
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
❖ General information
❖ Pile loading capacity
❖ Point load capacity
❖ Frictional resistance
❖ Negative Skin Friction
❖ Pile driving
❖ Pile load test
❖ Settlement of piles
9
I.1 General information (a) When one or more upper soil layers are
highly compressible and too weak to
support the load transmitted by the
Why use piles? superstructure.
(earth-retaining structures,
foundations of tall structures) (b) When bedrock is not encountered at a
reasonable depth below the ground
surface, the pile will be long with frictional
resistance at the soil–pile interface.
(c) When subjected to horizontal forces, pile
foundations resist by bending, while still
supporting the vertical load.
(d) When expansive and collapsible soils may
be present at the site of a proposed
structure. Piles may be used extended into
stable soil.
(e) The foundations of some structures
subjected to uplifting forces. Piles are
sometimes used for these foundations to
resist the uplifting force.
(f) Bridge abutments and piers are usually
(transmission towers,
constructed over pile foundations to avoid
offshore platforms, and the loss of bearing capacity (erosion).
basement mats below
the water table)
10
I.1 General information
11
I.1 General information
13
I.1 General information
14
I.1 General information
Qall = Ap f w
15
I.1 General information
GFRP tube
Fam et al (2003)
FRP rebar
16
I.1 General information
17
I.1 General information
Table. Typical design load ranges for different pile types and sizes
18
Break 19
I.2 Pile loading capacity – load transfer mechanism
Pile load test
(Mansur and Hunter, 1970)
= Q1 + Q2 Q
Q2
Q1 = Q – Q2
During pile load test:
s = 2.5 mm: Q = 700 kN, Q1 = 600 kN, Q2 = 100 kN Q Qu = 45%, Q1 Qs = 55%, Q2 Q p = 22%
s = 5 mm: Q = 1100 kN, Q1 = 800 kN, Q2 = 300 kN Q Qu = 71%, Q1 Qs = 73%, Q2 Q p = 67%
s = 11 mm: Q = 1550 kN, Q = 1100 kN, Q = 450 kN
u s p
Q1 is developed at a much smaller pile
Ultimate Skin friction Point displacement (or much faster) than Q2
20
I.2 Pile loading capacity – load transfer mechanism
The load on the pile is gradually increased from zero to Q(z=0) at the ground surface.
dQ( z ) = f( z ) pdz
Q(z) dQ( z )
f( z ) =
pdz
dQ(z) = f(z).p.dz
Q(z+dz)
21
I.2 Pile loading capacity – load transfer mechanism
The load on the pile Q(z=0) is increased up to yield value Qu
mostly
adopted
failure mode
❖ Pile foundations are deep foundations and
that the soil fails mostly in a punching mode
(a triangular zone I developed at the pile tip).
❖ In dense sands and stiff clayey soils, a radial
shear zone, II, may partially develop.
22
I.2 Pile loading capacity – pile length
The ultimate load-carrying capacity of a pile is given by
where Qu = Q p + Qs
Qp = load carried at the pile point
Qs = load carried by skin friction developed at the side of the pile (caused by shearing
resistance between the soil and the pile)
Extreme case:
The lengths of point
bearing piles up to Extreme case:
the rock-bed (or The lengths of friction
more embedded in piles depend on the
rock-bed). shear strength of the
soil, the applied load,
and the pile size.
Figure. (a-b) Point bearing piles; (c) friction piles
23
I.2 Pile loading capacity – general equation
Qp = load-capacity of the pile point
Qu = Q p + Qs Qs = frictional resistance (skin friction)
derived from the soil–pile interface.
Qs = pLf
Q p = Ap q p
(Point Bearing Capacity)
qu = q p = cN c + qN q + DN
25
I.3 Point bearing capacity – Derivation of Qp
Because the width D of a pile is relatively small, the term γDNγ* may be dropped from
the right side of the preceding equation without introducing a serious error; thus, we
have
q p = cN c + qN q + DN
Note that the term q for pile is the effective vertical stress at the pile tip, better to
expressed in q’. Thus, the point bearing of piles is
Q p = Ap q p = Ap ( cN c + qN q )
Where
Ap = area of pile tip
qp = unit point resistance
c’ = cohesion of the soil supporting the pile tip
q’ = effective vertical stress at the level of the pile tip
Nc*, Nq* = the bearing capacity factors
26
I.3 Point bearing capacity – Development of Qp
Ultimate Point Resistance, qp of Test Pile in sand at the Ogeechee River Site (Vesic 1970)
Q p = Ap q p = Ap ( cN c + qN q )
Variation of qp &
N60 with depth
27
I.3 Point bearing capacity – Development of Qp (Meyerhof’s Method)
Q p = Ap q p = Ap ( cN c + qN q )
Q p = Ap q p = Ap qN q
Q p = Ap qN q Ap ql
Q p N ccu Ap = 9cu Ap
where
1 1 + 2K0
m = ( v + 2 h ) = q ( 0 = v = q = z , K 0 = 1 − sin )
3 3
N = f ( , I rr )
I
where Irr = reduced rigidity index for the soil I rr = r
1 + Ir
Es Gs
Ir = rigidity index = =
2 (1 + s ) q tan q tan
Table. Bearing Capacity Factors Nσ* Based on the Theory of Cavity Expansion
30
I.3 Point bearing capacity – Development of Qp
Vesic’s Method for clay
In saturated clay (’=0), the net ultimate point bearing capacity of a pile can be
approximated as
Q p = Ap q p = Ap cu N c
4
N c = ( ln I rr + 1) + + 1
3 2
Ir
where Irr = reduced rigidity index for the soil I rr = (=0 for saturated clay)
1 + Ir
Es c
( 347 u − 33 300) Es s =0.5 Es
Ir = rigidity index = ⎯⎯⎯
2 (1 + ) q tan f =cu → =
3cu pa 3cu
s
= c + tan
Es = modulus of elasticity of soil f n
cu = undrained shear strength of saturated clay
pa = atmospheric pressure (=101.325 kPa)
31
I.3 Point bearing capacity – Development of Qp
Coyle and Castello’s Method for sand
Q p = Ap q p = Ap qN q
32
Example 1 Solution
Consider a 15-m long concrete pile
with a cross section of 0.45 m x 0.45
m fully embedded in sand. For the
sand, given: unit weight = 17 kN/m3;
(Figure left-down)
and soil friction angle ’ = 35o.
Estimate the ultimate point Qp with
each of the following:
a. Meyerhof’s method
b. Vesic’s method (m = 250)
c. The method of Coyle and Castello
d. based on the results of parts a, b, ( better to use m to replace 0 )
and c, adopt a value for Qp
143
33
Example 1 Solution
Consider a 15-m long concrete pile
with a cross section of 0.45 m x 0.45
m fully embedded in sand. For the
sand, given: unit weight = 17 kN/m3;
and soil friction angle ’ = 35o.
Estimate the ultimate point Qp with
each of the following:
a. Meyerhof’s method
b. Vesic’s method (m = 250)
c. The method of Coyle and Castello
d. based on the results of parts a, b,
and c, adopt a value for Qp
48
34
Example 2
Consider a pipe pile (flat driving point)
having an outside diameter of 406 mm. Solution
The embedded length of the pile in
layered saturated clay is 30 m. The
following are the details of the subsoil.
The groundwater table is located at a
depth of 5 m from the ground surface.
Estimate Qp by using
a. Meyerhof’s method
b. Vesic’s method
4
N c = ( ln I rr + 1) + + 1 = 11.57
3 2
35
I.3 Point bearing capacity – Development of Qp
Q p = Ap q p
Correlations for Calculating with SPT and CPT Results
On the basis of field observations, Meyerhof (1976) also suggested that the ultimate
point resistance in a homogeneous granular soil may be obtained from standard
penetration numbers as
L
q p = 0.4 pa N 60 4 pa N 60
D
N60 = the average value of the standard penetration number near the pile point (about
10D above and 4D below the pile point)
pa = atmospheric pressure ( = 101.325 kPa)
According to Briaud et al. (1985), in granular soil with the standard penetration
resistance N60
q p = 19.7 pa ( N 60 )
0.36
According to Meyerhof (1956), in granular soil with the cone penetration resistance qc
q p = qc
36
Example 3
Consider a concrete pile that is 0.305 m x 0.305 m in cross section in sand. The pile is 15.2 m
long. The following are the variations of N60 with depth.
a. Estimate Qp using q p = 0.4 pa N 60 L 4 pa N 60
D
b. Estimate Qp using
q p = 19.7 pa ( N 60 )
0.36
Solution
37
I.3 Point bearing capacity – Point Bearing Capacity of Piles Resting on Rock
Sometimes piles are driven to an underlying layer of rock. In such cases, the engineer
must evaluate the bearing capacity of the rock. The ultimate unit point resistance in rock
(Goodman, 1980) is approximately
Where
q p = qu N +1 ( )
Nф = tan2(45+ ф’/2)
qu = unconfined compression strength of rock
ф’ =drained angle of friction qu( design ) = qu( lab ) 5
Unconfined Compressive Strength of Rocks qu(lab) Typical Values of Angle of Friction ’ of Rocks
38
Break 39
I.4 Frictional Resistance – Derivation of Qs
Frictional Resistance, Qs
The frictional, or skin, resistance of a pile may be written as
Qs = pLf
Where
p = perimeter of the pile section
ΔL = incremental pile length over which p and f are taken to be constant
f = unit friction resistance at any depth z
p = D
L
p = 2(d1+d2)
40
I.4 Frictional Resistance – Frictional (Skin) Resistance in Sand
(L’≈15*D)
f = K o tan if 0 z L
f = f z = L if z L
41
I.4 Frictional Resistance – Frictional (Skin) Resistance in Sand
Effective Earth Pressure coefficient
Based on load test results in the field, Mansur and Hunter (1970) reported the following
average values of K.
H-piles . . . . . . K = 1.65
Steel pipe piles . . . . . . K = 1.26 Qs = pLf
f = K o tan
Precast concrete piles . . . . . . K = 1.5
f = f z = L
Replace the summation by average
f 0.93
1702
z
43
For 14/4
I.4 Frictional Resistance – Frictional (Skin) Resistance in Sand
Qs = pf av L
Correlation with Standard Penetration Test Results
Meyerhof (1976) indicated that the average unit frictional resistance, fav, for high
displacement driven piles may be obtained from average standard penetration
resistance values as
f av = 0.02 pa ( N 60 )
where
( N60 ) =
average value of standard penetration resistance.
Pa = atmospheric pressure (= 101.325 kPa or kN/m2)
For low-displacement driven piles
f av = 0.01 pa ( N 60 )
f av = 0.224 pa ( N 60 )
0.29
44
I.4 Frictional Resistance – Frictional (Skin) Resistance in Sand
Qs = pLf
Correlation with Cone Penetration Test Results
Nottingham and Schmertmann (1975) and Schmertmann (1978) provided correlations
for estimating Qs using the frictional resistance (fc) obtained during cone penetration
tests. According to this method
f = fc
The variations of α’ with z/D for electric cone and mechanical cone penetrometers are
shown in Figure below, respectively. We have
Qs = p ( L ) f = p ( L ) f c
Solution
If asked to determine Qall based on Meyerhof’s method and Briaud’s method (FS=3)?
46
Example 6
Consider an 18-m long concrete pile (cross section: 0.305 m x 0.305 m) fully embedded in a
sand layer. For the sand layer, the following is an approximation of the cone penetration
resistance qc (mechanical cone) and the frictional resistance with depth.
Estimate the allowable load that the pile can carry. Use FS=3.
Solution
0.44
47
I.4 Frictional Resistance – Frictional (Skin) Resistance in Clay
Qs = pf av L
λ Method
This method, proposed by Vijayvergiya and Focht (1972), Variation of with pile
is based on the assumption that the displacement of soil embedment length, L
caused by pile driving results in a passive lateral pressure
at any depth and that the average unit skin resistance is
f av = ( o + 2cu )
where
o = mean effective vertical stress for
the entire embedment length
cu = mean undrained shear strength (ф = 0)
In general case
o = ( A1 + A2 + A3 + ... ) L
cu = ( cu(1)L1 + cu(2)L2 + cu(3)L3 + ... ) L
48
I.4 Frictional Resistance – Frictional (Skin) Resistance in Clay
Qs = pLf
α Method
According to the method, the unit skin resistance in clayey soils can be represented by
the equation
f = cu
→α1 Variation of
(interpolated values
→α3 based on Terzaghi, Peck
and Mesri, 1996)
→α3
49
I.4 Frictional Resistance – Frictional (Skin) Resistance in Clay
Qs = pLf
β Method
When piles are driven into saturated clays, the pore water pressure in the soil around
the piles increases. The excess pore water pressure in normally consolidated clays may
be four to six times cu. However, within a month or so, this pressure gradually dissipates.
Hence, the unit frictional resistance for the pile can be determined on the basis of the
effective stress parameters of the clay in a remolded state (c’=0). Thus, at any depth,
f = o
Where
σ’o = verticaleffective stress
β = Ktanф’R
ф'R = drained friction angle of remolded clay
K = pressure coefficient
Qs = pLf = pLK tan R 0
50
Example 7
Refer to the pile in saturated clay shown in the figure below. For the pile, calculate the skin
resistance by (1) the -method, (2) the -method, and (3) the -method. For the -method,
use ’R = 30o for all clay layers. The top 10 m of clay is normally consolidated. The bottom clay
layer has an OCR = 2. (Note: diameter of pile = 406 mm)
51
Solution
52
Example 8
A concrete pile 305 mm x 305 mm in cross section is driven to a depth of 20 m below the
ground surface in a saturated clay soil. A summary of the variation of frictional resistance fc
obtained from a cone penetration test is as follows.
Estimate the frictional resistance Qs for the pile.
Solution
0.84
0.71
0.63
53
Break 54
I.5 Negative Skin Friction (NSF)
by integrating with depth in clay fill, the total downward drag force on a pile is
1
Qn =
Hf
f n pdz =
Hf
K f z tan pdz = K f tan pH 2f
0 0 2
where
Hf = Height of the fill (attention: if the fill is above water table, the ’f should be
replaced by the moist unit weight)
p = perimeter of the pile
56
I.5 Negative Skin Friction (NSF)
fn of sand fill is
(b) Granular soil fill over clay layer
not considered
the negative (downward) skin stress on the pile is
f n = K o tan
where
K’ = K0 = 1-sin’ (of clay layer)
’0 = ’fHf + ’z ( ’ = effective unit weight of clay layer)
’ = 0.5~0.7’ (of clay layer)
Solution
59
Solution (Attention: for end-bearing piles:
Example 9b L1 = L-Hf)
In the figure below, let Hf = 2 m
and L = 20 m. The pile is circular
in cross section with a diameter
of 0.305 m. For the fill that is
above the water table, and f =
16.5 kN/m3, ’clay = 34o and
sat(clay) = 17.2 kN/m3 Determine
the total drag force. ( ’ = 0.6’)
60
Break 61
I.6 Pile driving
Piles are driven into the ground by means of
❖ Hammers (for big piles)
❖ vibratory drivers (for small piles)
❖ jetting or partial augering (for small piles in soft soils)
❖ drilling (Bored Pile, for big piles)
62
I.6 Pile driving
Piles may be divided into two categories based on the nature of their placement:
❖ displacement piles
✓ Driven piles are displacement piles, because
❖ non-displacement piles they move some soil laterally; hence, there is
a tendency for densification of soil
surrounding them.
✓ Concrete piles and closed-ended pipe piles
are high-displacement piles.
✓ Steel H-piles displace less soil laterally during
driving, so they are low-displacement piles.
63
I.6 Pile driving
For pile-driving by hammers
❖ In the driving operation, a cap is attached to the top of the pile.
❖ A cushion may be used between the pile and the cap.
The cushion has the effect of reducing the impact force and
spreading it over a longer time; however, the use of the
cushion is optional.
❖ A hammer cushion is placed on the pile cap.
❖ The hammer drops on the cushion
Types of hammer
(a) drop hammer
(b) single-acting air of steam hammer
64
I.6 Pile driving
Types of hammer
(c) double-acting and differential air or steam hammer
(d) diesel hammer
(e) vibratory pile driver
65
I.6 Pile driving
Table. Examples of Commercially Available Pile-Driving Hammers
66
I.6 Pile driving
Pile-Driving Formulas: (a) EN formula
Dynamic equations are widely used in the field to determine whether a pile has reached a
satisfactory bearing value at the predetermined depth. The principle is (so-called
“(Engineering News Record (EN) Formula)”):
Energy imparted by the hammer per blow = (pile resistance)(penetration per hammer blow)
WR h
Qu =
S +C C = 25.4 mm (=1 inch) for drop hammers
WR = weight of the ram C = 2.54 mm for steam hammers
h = height of fall of the ram
S = penetration of pile per hammer blow the average value obtained from
C = a constant the last few driving blows.
For single- and double-acting hammers, the term WRH can be replaced by EHE
where
EH E
E is the efficiency of the hammer Qu =
HE is the rated energy of the hammer S +C
67
I.6 Pile driving
Pile-Driving Formulas: (b) Modified EN Formula
EWR h WR + n W p
2
Qu =
S + C WR + W p
E = efficiency of hammer
C = 2.54 mm if the units of S and h are in mm (whatever drop or steam)
Wp = weight of the pile with the cap
n = coefficient of restitution between the ram and the pile cap
Typical values for E
Single- and double-acting hammers 0.7–0.85
Diesel hammers 0.8–0.9
Drop hammers 0.7–0.9
Qu =
25.4
+ 2.54 22.24 + 40.08
N
9639.08
=
25.4
+ 2.54
N
Normal value of N:
4-5 for wooden piles
6-8 for concrete piles
12-14 for steel piles
69
I.6 Pile driving
Pile-Driving Formulas: (c) Danish formula EH E
Qu =
(Olson and Flaate, 1967) EH E L
S+
2 Ap E p
E = efficiency of hammer
HE = rated hammer energy
Ep = modulus of elasticity of the pile material
L = length of the pile
Ap = cross-sectional area of the pile
Wp EH E L
with Cd = 0.75 + 0.14 and = 2
R
W A E S
p p
70
Example 10
( = WR h )
W
EWR h R + n 2
Wp
Qu =
S + C WR + W p
EH E
Qu =
EH E L
S+
2 Ap E p
71
Solution
72
Break 73
I.7 Pile Load Test
❖ The primary reason is the unreliability of prediction methods.
❖ Force control: Step loads are applied (each step load = 1/4 of working load).
❖ Or displacement control: a constant rate of penetration 0.25-2.5 mm/min.
❖ The load test should be carried out to at least a total load of two times proposed
working load.
❖ After the desired pile load is reached, the pile is gradually unloaded.
74
I.7 Pile Load Test
snet = st − se
Curve (1): Pile settlement may increase with load to a certain point, beyond which
the load–settlement curve becomes vertical. The load corresponding to the point
where the curve of Q versus becomes vertical is the ultimate load, for the pile.
Curve (2): In many cases, the latter stage of the load–settlement curve is almost
linear, showing a large degree of settlement for a small increment of load. The
ultimate load, for such a case is determined from the point of the curve of Q versus
where this steep linear portion starts.
But not always clear? 75
I.7 Pile Load Test Total settlement
One of the methods to obtain the ultimate load Qu from the load-settlement plot is that
proposed by Davisson (1973). Davisson’s method is used more often in the field and is
described here. Referring to Figure below, the ultimate load occurs at a settlement level
(su) of
D Qu L
su ( mm ) = 0.012 Dref + 0.1 +
D A E
Where ref p p
Qu is in kN
D = pile diameter or width (in mm)
Dref = reference pile diameter
or width (= 300 mm)
L = pile length (mm)
Ap = pile cross section (mm2)
Ep = Young’s modulus of pile material
(kN/mm2) Slope su/Qu
= L/(ApEp)
76
Example 11
The figure below shows the load test results of a 20-m
long concrete pile (406 mm × 406 mm) embedded in
sand. Using Davisson’s method, determine the ultimate
load Qu. Given: Ep = 30x106 kN/m2.
Qu = 1460 kN
s = 3.735 mm
Solution
D Qu L
su = 0.012 Dref + 0.1 +
D
ref Ap E p Slope = 0.004
406 Qu ( 20000 )
= 0.012 300 + 0.1 +
300 406 406 30
= 3.735 + 0.004Qu
77
I.7 Pile Load Test
Soils surrounding the pile would be disturbed during pile driving, but they will be recovered
to some degree after certain time. So the test should be conducted after certain time.
Se(1) =
( Q
wp + Qws ) L
Ap E p
Qwp = load carried at the pile point under working load condition
Qws = load carried by frictional (skin) resistance under working load condition
Ap = area of cross section of pile
L = length of pile
Ep = modulus of elasticity of the pile material
The magnitude of ξ varies between 0.5 and 0.67 depending on the nature of the
distribution of the unit friction (skin) resistance f along the pile shaft.
81
I.8 Settlement of Piles
The total settlement of a pile under a vertical working load Qw is given by
Se = Se(1) + Se( 2) + Se( 3)
(1 − ) I
qwp D
Se( 2) = 2
s wp
Es
D = width or diameter of pile
qwp = point load per unit area at the pile point = Qwp/Ap
Es = modulus of elasticity of soil at or below the pile point
μs = Poisson’s ratio of soil
Iwp = influence factor ≈ 0.85
82
I.8 Settlement of Piles
The total settlement of a pile under a vertical working load Qw is given by
Se = Se(1) + Se( 2) + Se( 3)
Q D
Se( 3) = ws (1 − s2 ) I ws
pL Es
Qws = load carried by frictional (skin) resistance under working load condition
p = perimeter of the pile
L = embedded length of pile
D = width or diameter of pile
Es = modulus of elasticity of soil at or below the pile point
μs = Poisson’s ratio of soil
L
Iws = influence factor = 2 + 0.35
D
Qws Cs
Se( 3) = L
Lq p Cs = an empirical constant = 0.93 + 0.16 C p
D
by Vesic (1977), just for reference
83
I.8 Settlement of Piles
Example 12: The allowable working load on a prestressed concrete pile 21-m long that
has been driven into sand is 502 kN. The pile is octagonal in shape with D = 356 mm (see
Table below). Skin resistance carries 350 kN of the allowable load, and point bearing
carries the rest. Use Ep = 21 X 106 kN/m2, Es = 25 X 103 kN/m2, μs = 0.35 and ξ = 0.62.
Determine the settlement of the pile.
84
I.8 Settlement of Piles
From the table, for D = 356 mm, the area of pile cross section. Ap = 1045 cm2, Also,
perimeter p = 1.168m. Given: Qws = 350 kN, so
Se(1) =
(Qwp + Qws ) L 152 + 0.62 ( 350 ) ( 21)
= = 0.00353 m = 3.35 mm
Ap E p ( 0.1045 m 2
)( 21 10 6
)
152 0.356
(1 − ) I 3 (
1 − 0.352 ) ( 0.85 ) = 0.0155 m = 15.5 mm
qwp D
Se( 2) = 2
=
0.1045 25 10
s wp
Es
Q D
Se( 3) = ws (1 − s2 ) I ws I ws = 2 + 0.35
L 21
= 2 + 0.35 = 4.69
pL s
E D 0.356
0.356
3 (
1 − 0.352 ) ( 4.69 ) = 0.00084 m = 0.84 mm
350
Se( 3) =
(1.168 )( 21) 25 10
❖ Configuration
❖ Group Efficiency
❖ Ultimate Capacity of Group Piles in Saturated Clay
❖ Elastic Settlement of Group Piles
❖ Consolidation Settlement of Group Piles
88
II.1 Configuration
Interference between two piles
89
II.2 Group Efficiency
The efficiency of the load-bearing capacity of a group pile may be defined as
Qg ( u )
=
Q u
Where
ŋ = group efficiency
Qg(u) = ultimate load-bearing capacity of the group pile
Qu = ultimate load-bearing capacity of each pile without the group effect
(for pure friction pile)
Qg ( u ) f av pg L ( )
pg = 2 Lg + Bg f av 2 ( n1 + n2 − 2 ) d + 4 D L 2 ( n1 + n2 − 2 ) d + 4 D
= = ⎯⎯⎯⎯⎯⎯→ = =
Qu f av pL number of pile=n1n2 n1n 2 pLf av pn1n2
Attention: ≤ 1 !!!
90
II.2 Group Efficiency
Equation for Group Efficiency of Friction Piles
Converse-Labarre equation
( n1 − 1) n2 + ( n2 − 1) n1
= 1−
90 n n
1 2
where θ (deg) = tan-1(D/d)
Los Angeles Group Action equation
D
= 1− n1 ( n2 − 1) + n2 ( n1 − 1) + 2 ( n1 − 1)( n2 − 1)
dn1n2
11d n + n − 2 0.3
= 1 −
1 2
+
7 ( d − 1) n1 + n2 − 1 n1 + n2
2
where d is in ft
91
II.3 Ultimate Capacity of Group Piles in Saturated Clay
93
II.3 Ultimate Capacity of Group Piles in Saturated Clay
Qg (u ) = Lg Bg cu( p ) N c + 2 ( Lg + Bg ) cu L
94
II.3 Ultimate Capacity of Group Piles in Saturated Clay
Step 3. Compare the two values obtained from two methods. The lower value is taken
as the final Qg(u) for design.
Step 4. If step 1 is adopted, need to reduce it by .
Qg (u ) = Lg Bg cu( p ) N c + 2 ( Lg + Bg ) cu L
96
Solution of Example 1
97
Solution of Example 1
Remark: No good, in
practice should be Qu*
98
Break 99
II.4 Elastic Settlement of Group Piles
In general, the settlement of a group pile under a similar working load per pile increases
with the width of the group and the center-to-center spacing of the piles (d). Several
investigations relating to the settlement of group piles have been reported in the
literature, with widely varying results. The simplest relation for the settlement of group
piles was given by Vesic (1969), namely,
100
II.4 Elastic Settlement of Group Piles
For group piles in sand and gravel using standard penetration test results, for elastic
settlement, Meyerhof (1976) suggested the empirical relation
0.96q Bg
sg ( e ) (mm) = I with q =
Qg (u )
(in kN/m 2 )
N 60 Lg Bg
Similarly, the group pile settlement is related to the cone penetration resistance by the
formula
qBg
sg ( e ) = I Average cone penetration resistance
2qc within the seat of settlement.
102
II.5 Consolidation Settlement
of Group Piles
The consolidation settlement of a group
pile in clay can be estimated by using
the 2:1 stress distribution method. The
calculation involves 5 steps:
C C + av
Sc ( p ) = s log c + c log 0 Hc
1 + e0 0 1 + e0 c
103
Example 3
A group pile in clay is shown in the figure below. Determine the consolidation
settlement of the piles. All clays are normally consolidated.
104
Solution
105
Solution
106
Break 107
Civil Engineering 2020/2021
Subject: Geotechnical Design (CSE40403)
❖ Basic information
❖ Construction Procedures
❖ Estimation of Load-Bearing Capacity
❖ Elastic settlement of Drilled Shafts at Working Load
109
III.1 Basic information
Definition:
We use the term drilled shaft (or Bored Pile) for a hole drilled or excavated to the bottom
of a structure’s foundation and then filled with concrete. Depending on the soil conditions,
casings may be used to prevent the soil around the hole from caving in during construction.
The diameter of the shaft is usually large enough for a person to enter for inspection.
Advantages:
1. A single drilled shaft may be used instead of a group of piles and the pile cap.
2. Constructing drilled shafts in deposits of dense sand and gravel is easier than others.
3. Environment friendly (avoid impact, vibration, noise, displacements during construction).
4. Because the base of pile can be enlarged, it provides great resistance to the uplifting load.
5. The surface of the base of the drilled shaft can be visually inspected.
Drawbacks:
1. The concreting operation may be delayed by bad weather.
2. As in the case of braced cuts, deep excavations for drilled shafts may induce substantial
ground loss and damage to nearby structures.
110
III.1 Basic information
Types of drilled shaft
(a) straight shaft
(b) and (c) belled shaft
(d) straight shaft socketed into rock
111
III.2 Construction Procedures
112
III.2 Construction Procedures
113
III.2 Construction Procedures
115
III.2 Construction Procedures
Other Considerations
For an ordinary drilled shafts without casings, a minimum amount of vertical steel
reinforcement is always desirable. Minimum reinforcement is 1% of the gross cross-
sectional area of the shaft. For drilled shafts with nominal reinforcement, most building
codes suggest using a design concrete strength, fc, on the order of f ’c/4. Thus, the
minimum shaft diameter becomes (implied FS = 4 for pile structure)
f c = 0.25 f c
Qw D =
Ds = diameter of the shaft
Qw Qw Q
fc = = = 2.257 w f’c = 28-day concrete strength
Ags D 2
s
f c
( 0.25) fc Qw = working load of the drilled shaft
4
s
4 Ags = gross cross-sectional area of the shaft
116
III.3 Estimation of Load-Bearing Capacity
(1) General equation of the ultimate load-bearing capacity of a drilled shaft:
Qu = Q p + Qs (B is small compared to the depth)
1
q p = cN c Fcs Fcd Fcc + qN q Fqs Fqd Fqc + BN F s F d F c
2 Qp = load-capacity at the base
Qu = Ap ( cN c Fcs Fcd Fcc + qN q Fqs Fqd Fqc ) + Qs
Qs = frictional (skin) resistance
118
III.3 Estimation of Load-Bearing Capacity
(2) Ultimate load-bearing capacity of a drilled shaft in Granular Soil
❖ Estimation of Qp
Method of Chen and Kulhawy (1994)
119
III.3 Estimation of Load-Bearing Capacity
(2) Ultimate load-bearing capacity of a drilled shaft in Granular Soil
❖ Estimation of Qp
Method of Chen and Kulhawy (1994)
120
III.3 Estimation of Load-Bearing Capacity
(degree)
(2) Ultimate load-bearing capacity of a drilled shaft in Granular Soil
❖ Estimation of Qp N q* = 0.21e0.17
Method of Berezantzev et al. (1961)
Q p( net ) = Ap q ( N q* − 1)
121
III.3 Estimation of Load-Bearing Capacity
(2) Ultimate load-bearing capacity of a drilled shaft in Granular Soil
❖ Estimation of Qs p = shaft perimeter = πDs
f = unit frictional (or skin) resistance = Kσ’otanδ’
Qs = pfdz K = earth pressure coefficient ≈ Ko = 1-sin’
L1
0
σ’0 = effective vertical stress at any depth z
0 0
Q p ( net ) + Qs
Qall(net) =
FS
122
III.3 Estimation of Load-Bearing Capacity
(2) Ultimate load-bearing capacity of a drilled shaft in Granular Soil
Example 1a: A soil profile is shown in the figure
below. A point bearing drilled shaft with a bell is
placed in a layer of dense sand and gravel. Determine
the allowable load the drilled shaft could carry. Use
the Method of Chen and Kulhawy (1994) and a factor
of safety of 4. Take and For the dense sand layer,
Ignore the frictional resistance of the shaft.
Solution:
123
124
III.3 Estimation of Load-Bearing Capacity
(2) Ultimate load-bearing capacity of a drilled shaft in Granular Soil
Example 1b: A soil profile is shown in the figure below.
A point bearing drilled shaft with a bell is placed in a
layer of dense sand and gravel. Determine the
allowable load the drilled shaft could carry. Use the
Method of Berezantzev et al. (1961) and a factor of
safety of 4. Take and For the dense sand layer, Ignore
the frictional resistance of the shaft.
Solution:
125
Break 126
III.3 Estimation of Load-Bearing Capacity
(3) Ultimate load-bearing capacity of a drilled shaft in Clay For simplicity
❖ Estimation of Qp cu
1 then N c* = 9
(
Q p( net ) Ap cN c Fcs Fcd Fcc + q ( N q − 1) Fqs Fqd Fqc ) pa
c=c , =0 N =1
⎯⎯⎯⎯⎯⎯
u q
→ Q p( net ) Ap cu N c Fcs Fcd Fcc
Approximation to:
Q p( net ) = Ap cu N c*
Conservatively, * = 0.4
127
III.3 Estimation of Load-Bearing Capacity
(3) Ultimate load-bearing capacity of a drilled shaft in Clay
( *
= 0.4 )
128
L1 = 8.23 m
III.3 Estimation of Load-Bearing Capacity L2 = 2.59 m
Ds = 1.0 m
(3) Ultimate load-bearing capacity of a drilled shaft in Clay cu (1) = 50 kN/m 2
cu (2) = 108.75 kN/m 2
Solution:
129
III.3 Estimation of Load-Bearing Capacity
(4) Drilled Shafts extending into Rock
Ultimate load-bearing capacity of a drilled shaft in rock:
Q p ( MN ) = 4.83 ( qu MN/m 2 ) Ap ( m 2 )
0.51
Solution:
131
III.4 Elastic settlement of Drilled Shafts at Working Load
The total settlement of a pile under a vertical working load Qw is given by
Se = Se(1) + Se( 2) + Se( 3)
Se(1) =
( Q
wp + Qws ) L
Ap E p
(1 − s2 ) I wp or Se( 2) =
qwp D Qwp C p
Se( 2) =
Es Dq p
Qw = 1005 kN
Example 4: The figure below shows a drilled
shaft without a bell. Estimate the elastic
settlement at the working load. (taking =
L1 = 8 m
0.65, Ep = 21x106 kN/m2, Es = 14,000 kN/m2,
s = 0.3, Qwp = 250 kN) L2 = 3 m
Ds = 1.5 m
Se(1) =
(Q wp + Qws ) L cu (1) = 50 kN/m 2
Ap E p cu (2) = 105 kN/m 2
Qwp C p
Se( 2) =
Dq p
Q D
Se( 3) = ws (1 − s2 ) I ws
pL Es
133
Qw = 1005 kN
Solution:
L1 = 8 m
L2 = 3 m
Ds = 1.5 m
cu (1) = 50 kN/m 2
cu (2) = 105 kN/m 2
= 0.65
Ep = 21x106 kN/m2
Es = 14,000 kN/m2
s = 0.3
Qwp = 250 kN
134
Break 135
Tutorial A driven closed-ended pile, circular in cross section, is shown in the Figure bleow.
exercise 1 Calculate the following.
a. The ultimate point load using Meyerhof’s procedure.
b. The ultimate point load using Vesic’s procedure. Take Irr = 50.
c. An approximate ultimate point load on the basis of parts (a) and (b).
d. The ultimate frictional resistance Qs. [Use Eqs. (1) through (2), and take K = 1.4
and δ’ = 0.6’.]
e. The allowable load of the pile (use FS = 4).
L = 15 D (1)
f = f z = L ( 2)
f
=381 mm 136
Example 1
137
Solution for Example 1
Qs = pLf =
( f z =0 + f z =3.05 ) pL + ( f z =3.05 + f z =5.715 ) p
1 ( L − L1 ) +f z =5.715 p ( L − L) z
2 2
0 + 23.38 23.38 + 34.33
= ( 0.381 ) 3.05 + ( 0.381 ) ( 5.715 − 3.05 )
2 2
+ 34.33 ( 0.381 ) (15.24 + 3.05 + 3.05 − 5.715 ) (Part of constant f, smaller ’ used
for safe design, but not realistic)
= 42.68 + 92.04 + 642.05 = 776.77 kN 139
Solution for Example 1
47.946 23.38
z = L’ 70.41 34.33
43.89
1
Solution for Exercise 1
According to the rule of ASCE:
(1) Determine the depth z1 at which D’= 0.1*q (q= net stress on the foundation).
52
Q 3 1
D z = 2 I P with I P =
z 2 1 + ( r z )2
3 Q
D z = (because of r = 0)
2 z 2
Q Q
q= =
BL 4
3 Q Q z independent of
D z = 0.1q = 0.1 z = 4.37 m
2 z 2 4 applied load!
0 = ( sat − w ) z = 7 z
3 100
D z = 0.05 0 = 0.05 7 z z = 5.15 m
2 z 2
3
Solution for Exercise 2
For a Shelby tube, given: outside diameter = 76.2 mm and inside diameter = 73 mm.
What is the area ratio of the tube?
Do = 76.2 mm
Di = 73 mm
( 76.2 ) − ( 73)
2 2
D −D
2 2
Ar = o
=
i
= 8.96 %
( 73)
2 2
Di
Change the “inside diameter = 73 mm” by “the wall thickness of the tube = 1.63 mm”?
Ar = = = 9.14 %
( 72.94 )
2
Di2
4
Tutorials
Exercise 3
A soil profile is shown in Figure below along with the standard penetration numbers in
the clay layer. Use the following equations to determine the variation of cu and OCR
with depth. What is the average value of cu and OCR?
0.689
N 60
cu = 0.29 N 600.72 pa OCR = 0.193
0
= 100 kPa in MPa
Atmospheric pressure
0.689
N 60
OCR = 0.193
0
When z = 3 m
0 = z1 + ( sat − w ) z2 = 16.5 1.5 + (19 − 9.81) 1.5 = 38.535 kPa
cu = 0.29 ( N 60 ) pa = 0.29 50.72 100 = 92.4 kPa
0.72
0.689
N
0.689
5
OCR = 0.193 60 = 0.193 = 5.5
0 0.038535
Dp
E p = 2 (1 + m s )(Vo + vm )
Dv
9
Solution for Exercise 5
Dp
E p = 2 (1 + m s )(Vo + vm )
Dv
vo + v f
vm = = 113 10−6 m3
2
Dp = p f − po = 284.1 kN/m 2
Dv = v f − vo = 134 10−6 m3
Dp −6 284.1
E p = 2 (1 + m s )(Vo + vm )
Dv
= 2 (1 + 0.5 ) ( 535 10 −6
+ 113 10 ) −6
134 10
= 4121 .57 kN/m 2
10
Tutorials
Exercise 6
The results of a refraction survey at a site are given in the following table. Please
determine the P-wave velocities and the thickness of the material encountered.
Distance, x (m)
11
Solution for Exercise 6
Plot the figure:
Chart Title
80
Time of first arrival, t (sec x 10-3)
70
60
Ti2
50
40
30
20
10
0
xc
0 10 20 30 40 50
Distance, x (m)
12
Solution for Exercise 6
Velocity
In the Figure, the times of first arrival of the P waves are plotted against the
distance of the geophone from the source of disturbance. The plot has three
straight-line segments. The velocity of the top three layers can now be calculated as
follows:
13
Can’t
Solution for Exercise 6
Thickness of Layers
14
Exercise 1: Estimate the factor of safety of the translational slope Tutorials
below according to an arbitrary slice of soil block.
Force against sliding Shear Strength
(1) Fsf = (2) Fsτ =
Force of sliding Shear stress
z
Unit weight: g
shear strength: c’, f’
(water table below
the bedrock)
f f fl
Fs = = = = ......
T l T
2
Exercise 2: The figure below shows a selected slice of soil block Tutorials
above an assumed circular slip surface.
O d
F A’ B’
R
y
b
a
x
Unit weight: g A B
shear strength: c’, f’
a
b
(a) Derive formulations for the driving moment Md and (x1, y’1) (x2, y’2)
resisting moment Mr of the selected slice of soil
(xm, y’m)
block, assuming AB a straight-line.
Given coordinates: O(11,20), F(14,14),and soil properties:
c’ =10 kPa, f’=35⁰, g = 11 kN/m3, z
the given slice
(b) For the slice from AA’(x=22) to BB’(x=24), calculate
the Md and Mr?
(c) For the slice from AA’(x=6) to BB’(x=8), calculate the (x2, y2)
Md and Mr? (x1, y1) (xm, ym)
3
Solution of Exercise 2(a) Tutorials
(a) Derive formulations for the driving O d
moment Md and resisting moment Mr of F A’ B’
the selected slice of soil block, assuming R
AB a straight-line. y
b
a
x
A B
a
b
W = g zb
(x1, y’1) (x2, y’2)
N = W cos a = g zb cos a T
T = W sin a = g zb sin a (xm, y’m)
N W
y1 = 20 − 22.82542 − ( 22 − 11) = 0
2
(xm, y’m)
b = x2 − x1 = 2
xm − xc −1 23 − 11
a = sin −1 = sin = 0.55 (= 31.72 )
o
W = g zb a b
N = W cos a = g zb cos a T
(x2, y’2)
T = W sin a = g zb sin a
W N
(xm, y’m)
(x1, y’1)
Tf
z
N g zb cos b T = f l = ( c + n tan f ) b cos a
n = = = g z cos 2 b f
l b cos b = c b cos a + g zb cos a tan f (x1, y1)
(xm, ym)
(x2, y2)
M r = T f R = ( c b cos a + g zb cos a tan f ) R
M d = Wd = g zbR sin a ( or M d = TR = g zb sin a R ) 6
Solution of Exercise 2(c) Tutorials
d O
Given coordinates: O(11,20), F(14,14),and soil R
properties: c’ =10 kPa, f’=35⁰, g = 11 kN/m3, F
(c) For the slice from AA’(x=6) to BB’(x=8),
calculate the Md and Mr? B’
A’
y
M r = ( c b cos a + g zb cos a tan f ) R a
x
M d = g zbR sin a
A B
R = x + y = 11 + 20 = 22.8254
2 2 2 2
c c
a b
yi = yc − R 2 − ( xi − xc )
2
T (x2, y’2)
y1 = 20 − 22.82542 − ( 6 − 11) = −2.27
2
W N (xm, y’m)
y2 = 20 − 22.8254 − ( 8 − 11) = −2.63
2 2
(x1, y’1)
At the slope: y = x tan ( 4 ) = x y1 = 6, y2 = 8
Tf
z = ym − ym = ( 6 + 8 ) 2 − ( −2.27 − 2.63) 2 = 9.45 z
b = x2 − x1 = 2 (x1, y1)
(xm, ym)
x −x −1 7 − 11
a = sin m c
−1
= sin = −0.18 (= −10.09 )
o
c’=10 kPa
f’=35°
gd=11 kN/m3
x = 22m, 24m
(2) Mr=? Md=?
(1) x = 6m, 8m
Please calculate the resisting moment Mr Mr=? Md=?
and driving moment Md for two given slices.
− ymi
hi = ymi
(R = x +y
2
c
2
c ) bi = xi +1 − xi
yi = yc − R 2 − ( xi − xc ) xmi − xc
2
a i = sin −1
R M d i = bi hig sin a i R
yi = xi tan b Average of “i” & “i+1”
li = bi cos a i M r i = ( li c + bi hig cos a i tan f ) R
93
Tutorials
LEM – Method of slices (Swedish Method)
(to establish an excel file)
c’=10 kPa
f’=35°
gd=11 kN/m3
x = 22m, 24m
x = 6m, 8m Mr=? Md=?
(a) The reservoir is empty and soil of fill is dry? Mr=? Md=?
94
Tutorials
LEM – Method of slices (Swedish Method)
96
Tutorials
LEM – Method of slices (Swedish Method)
b h g sin ai i i (xmi,y’mi)
hi
Wi
A
(0,0) li
a+ - ai
Sm
(xmi,ymi) N
Fm ,bishop =
( cl + N tan f ) R = cl cos a + (W − u l cos a ) tan f
i i i i i i i i ma
W sin a Ri i W sin a
i i
bi = li cos a i
tan f sin a i tan f tan a i
with ma = cos a i + = cos a i 1 + Wi = bi hig
Fm F m
ui li = g w hi bi sec a i
98
Tutorials
LEM – Method of slices (from Swedish to Bishop)
99
Tutorials
LEM – Method of slices (Swedish Method)
(e) For the last case of rapid emptying, please try to find the minimum value of Fs of
Bishop for the rotation center within the domain of (-5, 15), (5, 15), (-5, 25), (5, 25)
using the grid-point method. (online teaching tool is allowed)
(f) If the FOS is smaller than 1.1, please state possible design modifications for the slope.
Assignment 2:
A final report needed for all questions from (a) to (f).
Deadline for submitting report: 12th March 2021.
100
Exercise 1 For the following cases, determine the allowable gross vertical load-bearing
capacity of the foundation. Use Terzaghi’s equation and assume general shear
failure in soil. Use FS = 4 (based on gross..).
Part B Df 𝜙′ 𝒄′ γ Foundation type
a. 1.22 m 0.91 m 25° 28.75 kN/m2 17.29 kN/m3 Continuous
b. 2m 1m 30° 0 17 kN/m3 Continuous
c. 3m 2m 30° 0 16.5 kN/m3 Square
Part a: qall = 252.6 kN/m2
Part b: qall = 176.8 kN/m2
Part c: qall = 280.0 kN/m2
Qall = ?
Part c
From the equation :
qu = 1.3cN c + qN q + 0.4 BN
From the Table, when = 30 ,
N c = 37.16
N q = 22.46
N = 19.13
Thus :
qu = 1.3 0 37.16 + ( 2 16.5 ) 22.46 + 0.4 16.5 3 19.13
= 1119.95 kN m 2
q 1119.95
qall = u = = 280.0 kN m 2
FS 4
Qall = qallB2 =
Exercise 2 A square column foundation has to carry a gross allowable
′
load of 1805 kN (FS = 3,
′
based on gross..). Given: Df = 1.5 m, γ = 15.9 kN/m3, 𝜙 = 34°, and 𝑐 = 0. Use
Terzaghi’s equation to determine the size of the foundation (B). Assume general
shear failure. B = 2.0 m
B B BB B2
From the equation :
2
qu = 1.3cN c + qN q + 0.4 BN
3
2
= 1.3 0 18.92 + ( 2 16.5 ) 8.26 + 0.4 16.5 B 4.31
3
= 272.58 + 28.446 B
Combine the two equations :
13320
qu (net) = qu − q = 272.58 + 28.446 B − 2 16.5
B2
B = 5.75 m
Exercise 4 A′column foundation (see Figure) is 3 m × 2 m in plan. Given: Df = 1.5 m, 𝜙 ′ = 25°,
𝑐 = 70 kN/m2. Using Terzaghi′s bearing capacity equations assuming general
shear failure and FS = 3 (based on net..), determine the net allowable load the
foundation could carry. (w = 9.81 kN/m3) Q = 5240.87 kN
net
2 BL 2 2 3
B = = = 2.4 m
B+L 2+3
q = D1 + D2 ( sat − w ) = 1 17 + (1.5 − 1) (19.5 − 9.81) = 21.845kN m 2
From the equation :
qu = 1.3cN c + qN q + 0.4 B* N
From the Table, when = 25 ,
N c = 25.13
N q = 12.72
N = 8.34
Thus :
qu = 1.3 70 25.13 + 21.845 12.72 + 0.4 (19.5 − 9.81) 2.4 8.34
= 2642.28kN m 2
qu − q ( 2642.28 − 21.845 )
qall ( net ) = = = 873.48 kN m 2
FS 3
Qnet = qall ( net ) BL = 880.76 2 3 = 5240.87 kN
Exercise 5 An eccentrically loaded foundation is shown in the Figure. Use FS of 4 (based on
gross..) and determine the maximum allowable load that the foundation can carry.
Use Meyerhof’s effective area method.
F d = 1
Thus :
qu = (17 0.8 ) 23.18 1.542 1.147 1
+ 0.5 17 1.3 30.22 0.653 1
= 775.63kN m 2
Qult = qu BL = 775.63 1.3 1.5 1512 kN
Qult 1512
Qall = = = 378 kN
FS 4
Exercise 1
Refer to the Figure below. Using the equations outlined below, determine the average
stress increase in the clay layer below the center of the foundation due to the net
foundation load of 900 kN Δσ’ =68.60 kN/m2 av
H 2 I a( H 2 ) − H1I a( H1 ) (Attention: this equation
( H H ) = qo
av is based on corner…)
2 1
H 2 − H1
B L
m2 = ; n2 =
H H
general, can be H1 or H2
Clay
1
Solution for Exercise 1
P 900
qo = = = 268.74 kN m 2
A 1.83 1.83
H 2 I a( H 2 ) − H1I a( H1 ) 4.27 I a( H 2 ) − 1.22 I a( H1 )
(H
For 1 / 4 of footing : av H1 )
= qo = 268.74
2
H 2 − H 1 3 . 05
For I a( H1 ) and I a( H 2 ) :
B1 1.83 2
m2(1) = n2(1) = = = 0.75
H1 1.22
B1 1.83 2
m2( 2) = n2( 2 ) = = = 0.214
H2 4.27
(H
For whole footing : 4 av H1 )
= 4 17.15 = 68.60 kN m 2
2
2
Exercise 1
Refer to the Figure below. Using the equations outlined below, determine the average
stress increase in the clay layer below the center of the foundation due to the net
foundation load of 900 kN Δσ =68.60 kN/m2
av
H 2 I a( H 2 ) − H1I a( H1 )
( H H ) = qo
av
2 1
H 2 − H 1
B L
m2 = ; n2 =
H H
4
Solution for Exercise 2
qo = H = 17 10 = 170kN m 2
1) Stress increase at A
B1 = 3 m; B2 = 10 2 = 20 m.
B1 3 B2 20
= = 0.6 and = =4
z 5 z 5
I = 0.472
= (1) + (2)
(
= qo I (left side ) + I (right side ) )
= 2 170 0.472 = 160.48 kN m 2
5
Solution for Exercise 2
2) Stress increase at B
At left side, B1 = 0 m; B2 = 20 m.
B1 0 B2 20
= =0 and = =4
z 5 z 5
I (left side ) = 0.415
At right side, B1 = 6 m; B2 = 10 2 = 20 m.
B1 6 B2 20
= = 1.2 and = =4
z 5 z 5
I (right side ) = 0.488
3) Stress increase at C
At middle section, B1 = 26 m; B2 = 20 m.
B1 20 + 6 B2 20
= = 5.2 and = =4
z 5 z 5
I = 0.5
B1 0 B2 20
= =0 and = =4
z 5 z 5
I (extra ) = 0.415
9
Exercise 3
10
Solution for Exercise 3
1.91
B = = 0.955 m; = 4;
2
L 3.05
m = = = 1.60
B 1.91
H 9.76
n = = = 10.22
B 1.91
2 2
F1 = 0.597; F2 = 0.025;
1 − 2s 1 − 2 0.3
I s = F1 + F2 = 0.597 + 0.025 = 0.611
1 − s 1 − 0.3
Df 0.76 B 1.91
= = 0.4; = = 0.626
B 1.91 L 3.05
I f is located in 0.81 ~ 0.86
0.625 − 0.5
If = ( 0.81 − 0.86 ) + 0.86 = 0.8475
1 − 0.5
1 − s2 1 − 0.32
Se( flexible,center ) = qo ( B ) I s I f = 144 ( 4 0.955 ) 0.611 0.8475
Es 22080
= 0.0117 m = 11.7 mm
Se( rigid ) 0.93Se( flexible,center ) = 0.93 11.7 = 10.9 mm 11
Exercise 4
A square column foundation is shown in following figure. Determine the average increase
of pressure in the clay layer below the center of the foundation:
(a) by using elastic theory (for ’av, elastic simplified method)
(b) by using the 2:1 method (for ’av, 2:1 simplified method)
(c) Estimate the 1-D consolidation settlement of the clay layer using the results (a)
(d) Assume that the primary consolidation settlement is completed in 5 years. Also let
Cα=0.003. Estimate the secondary consolidation settlement at the end of 12 years.
= 4qo I or =qo I c ( elastic theory )
(same as ex.1)
qo B L
= (2 : 1 method)
( B + z )( L + z )
= ( t + 4 m + b ) 6
av
Cs H c C H + av
Sc( p ) = log c + c c log o
1 + eo o 1 + eo c
Sc + av
e = (1 + eo ) ; e = Cs log c + Cc log o
Hc o c
C H c t
Sc( s ) = log 2
1 + ep t1 12
Exercise 4
13
Exercise 4
14
Exercise 4
15
Solution for Exercise 4
16
Solution for Exercise 4
Cs H c C H + av
Sc( p ) = log c + c c log o
1 + eo o 1 + eo c
0.06 3050 100 0.25 3050 110.02
= log + log
1 + 0.68 49.75 1 + 0.68 100
= 52.17 mm
17
Solution for Exercise 4
(d) Cα=0.003
The change of void ratio of clay during primary consolidation, Δe
Sc 52.17
e = (1 + eo ) = (1 + 0.68 ) = 0.0287
Hc 3050
+ av
or for this case, e = Cs log c + Cc log o
o c
C H c t 0.003 3050 12
Sc( sec ondary ) = log 2 = log = 2.11 mm
1 + ep t1 1 + 0.6513 5
18
Exercise 1
The figure below shows a cantilever sheet pile wall penetrating a granular soil. Here, L1 = 2 m,
L2 = 3 m, = 15.9 kN/m3, sat = 19.33 kN/m3, and f’ = 32°.
(a) What is the theoretical depth of embedment, D?
(b) For a 30% increase in D, what should be the total length of the sheet piles?
(c) What should be the minimum section modulus of the sheet piles? Use all = 172 MN/m2.
14
Solution
Part a
Using the Figure of the pressure
distribution diagram, one can now
prepare the following table for a
step-by-step calculation.
15
Solution
Part b
The total length of the sheet pile is
Part c
Finally, we have the following table:
16
Tutorial Figure below shows a cantilever sheet pile wall penetrating a granular soil. Here,
Exercise 1 L1 = 4 m, L2 = 8 m, = 16.1 kN/m3, sat = 18.2 kN/m3, and f’= 32°
a. What is the theoretical depth of embedment, D?
b. For a 30% increase in D, what should be the total length of the sheet piles?
c. Determine the theoretical maximum moment of the sheet pile.
f f ME
K a = tan 2 45 − ; K p = tan 2 45 + ; z =
2 2 P
5 = ( L1 + L2 ) K p + L3 ( K p − K a )
5 8P
A1 = ; A2 =
( K p − Ka ) ( K p − Ka )
6 P 2 z ( K p − K a ) + 5 P ( 6z 5 + 4P )
A3 = ; A4 =
2 ( K p − K a ) 2 ( K p − K a )
2 2
Dtheory = L3 + L4
2P
z =
( K p − Ka )
1 z
M max = P ( z + z ) − z2 ( K p − K a )
2 3
Tutorials
Exercise 2
In the figure below, for the sheet pile wall penetrating into saturated clay, determine
(a) The theoretical and actual depth of penetration, using Dactual = 1.5Dtheory.
(b) The minimum size of sheet pile section necessary, using all = 172.5 MN/m2.
25
Solution
26
Solution
27
Tutorial Refer to Figure below, for which L1 = 2.4 m, L2 = 4.6 m, = 15.7 kN/m3, sat = 17.3
Exercise 2 kN/m3 and f’= 30°for sand, and c = 29 kN/m2 for clay.
a. What is the theoretical depth of embedment, D?
b. Increase D by 40%. What length of sheet piles is needed?
c. Determine the theoretical maximum moment in the sheet pile.
f ME
K a = tan 2 45 − ; z =
2 P
P1 ( P1 + 12cz1 )
D 2 4c − ( L1 + L2 ) − 2 DP1 − =0
( L1 + L2 ) + 2c
P1
z =
6
D 4c − ( L1 + L2 ) − P1
L4 =
4c
6 z 2
M max = P1 ( z + z1 ) −
2
Exercise 1 Refer to the braced cut shown in the figure below. Given: g = 16 kN/m3, f’=
38°and c’ = 0. The struts are located at 3.5 m center-to-center in the plan. Draw
the earth-pressure envelope and determine the strut loads at levels A, B, and C.
f
K a = tan 2 45 −
2
a = 0.65g HK a
M B1 = 0; M B2 =0
Strut load = reactions horizontal spacing
Solution for Exercise 1
f 38
K a = tan 2 45 − = tan 2 45 − = 0.238
2 2
a = 0.65g HK a = 0.65 16 (1 + 2.5 + 3 + 1.5 ) 0.238 = 19.79 kN m 2
M B1 =0
1 + 2.5
19.79 (1 + 2.5 )
A= 2
= 48.49 kN m
2.5
B1 = 19.79 (1 + 2.5 ) − 48.49 = 20.78kN m
M B2 =0
3 + 1.5
19.79 ( 3 + 1.5 )
C= 2
= 66.80 kN m
3
B2 = 19.79 ( 3 + 1.5 ) − 66.80 = 22.27 kN m
The strut loads are
At A, 48.49 ( spacing ) = 48.49 3.5 = 169.71 kN
At B, ( B1 + B2 ) ( spacing ) = ( 20.78 + 22.27 ) 3.5 = 150.67 kN
At C , 66.80 ( spacing ) = 66.80 3.5 = 233.79 kN
Exercise 2 Refer to the figure below. For the braced cut, given: H = 6 m, Hs = 2.5 m; gs = 16.5
kN/m3; angle of friction of sand f’s = 35°; Hc = 3.5 m, gc = 19.2 kN/m3 and
unconfined compression strength of clay layer, qu = 62 kN/m2.
a. Estimate the average cohesion (cav) and average unit weight (gav) for the
construction of the earth-pressure envelope.
b. Plot the earth-pressure envelope.
1
cav = g s K s H s2 tanfs + ( H − H s ) nqu
2H
K s : a lateral earth pressure coefficient for the sand layer ( 1)
ranging from 0.5 to 1.0,
n: a coefficient of progressive failure
average value 0.75
1
g av = g s H s + ( H − H s ) g c
H
gH 4c
4: a = max g H 1 − , 0 .3g H
c g H
gH
4: a = 0.3g H
c
Solution for Exercise 2
Part a:
1
cav = g s K s H s2 tanfs + ( H − H s ) nqu
2H
16.5 1 2.52 tan ( 35 ) + ( 6 − 2.5 ) 0.75 62
1
=
26
= 19.58kN m3
1 1
g av = g s H s + ( H − H s ) g c = 16.5 2.5 + ( 6 − 2.5 ) 19.2
H 6
= 18.075kN m3
Part b:
g av H 18.075 6
= = 5.54 4
cav 19.58
4c
a = max g H 1 − , 0 .3g H
gH
4c 4 19.58
a = g H 1 − = 18.075 6 1 − = 30 .13kN m 2
g H 18.075 6
a = 0.3g H = 0.3 18.075 6 = 32.54 kN m3
Thus,
a = 32.54 kN m3
Exercise 3 Refer to the figure below. Given: g = 17.5 kN/m3, c = 30 kN/m2. and center-to-
center spacing of struts in the plan = 5 m. Draw the earth-pressure envelope and
determine the strut loads at levels A, B, and C.
gH 4c
4: a = max g H 1 − , 0 .3g H
c gH
gH
4: a = 0.3g H
c
M B1 = 0; M B2 = 0
Strut load = reactions horizontal spacing
gH gH
4 4
c c
Solution for Exercise 3
1m
=1.75 m
A
Part a:
3m
H = 1+ 3 + 2 +1 = 7 m
g H 17.5 7 B
= = 4.08 4 =5.25 m
c 30
2m
4c
a = max g H 1 − , 0 .3g H C
g H
1m
4c 4 30
a = g H 1 − = 17.5 7 1− = 2.5kN m
2
g H 17.5 7
a = 0.3g H = 0.3 17.5 7 = 36.75kN m 3
Thus
a = 36.75kN m3
Part b:
To calculate the strut loads, taking the moment about B1 , we have M B1 = 0, and
1 0.25 H 1 + 3 − 0.25 H
A ( 3) − 36.75 0.25 H 1 + 3 − 0.25 H + − (1 + 3 − 0 .25 H ) 36 .75
2 3 2
1 2 4 − 0.25 7
3 A = 36.75 0.25 7 4 − 0.25 7 + ( 4 − 0.25 7 ) 36.75
2 3 2
A = 61.38kN m
Solution for Exercise 3
1m
=1.75 m
A
Also, horizontal forces = 0. Thus,
3m
1
36.75 0.25 H + 36.75 (1 + 3 − 0.25 H ) = A + B1
2 B
1 =5.25 m
36.75 0.25 7 + 36.75 (1 + 3 − 0.25 7 ) = 61.38 + B1
2m
2
B1 = 53.46 kN m C
Taking the moment about B2 , we have M B2 = 0, and
1m
1+ 2
C ( 2 ) − (1 + 2 ) 36.75 2C = 3 36.75 1.5
2
C = 82.69 kN m
Also, horizontal forces = 0. Thus,
36.75 (1 + 2 ) = C + B2
36.75 (1 + 2 ) = 82.69 + B2
B2 = 27.56 kN m
Thus
At A, PA = 61.38 ( spacing ) = 61.38 5 = 306.9 kN
At B, PB = ( B1 + B2 ) ( spacing ) = ( 53.46 + 27.56 ) 5 = 405.1 kN
At C , PC = 82.69 ( spacing ) = 82.69 5 = 413.45 kN
Exercise 4 Determine the sheet-pile section modulus for the braced cut described in Exercise
3. Use σall = 170 MN/m2
M max
S=
all
gH gH
4 4
c c
Solution for Exercise 4
1m
1m
=1.75 m
Part a A
For the maximum moment, the shear
3m
force should be zero. The location of
point E and F can be given as B
=5.25 m
2m
x
53.46 C
B1 − a x = 0 x = = 1.45 m;
1m
36.75
27.56
B2 − a y = 0 y = = 0.75 m
36.75 y
Also,
1 36.75 1
Magnitude of moment at A = 1 1
2 1.75 3
= 3.5kN m meter of wall
1m
and
x 1.45
Magnitude of moment at E = B1 x − a x = 53.46 1.45 − 36.75 1.45
2 2
= 38.88kN m meter of wall
1
Magnitude of moment at C = 36.75 1 = 18.375kN m meter of wall
2
y 0.75
Magnitude of moment at F = B2 y − a y = 27.56 0.75 − 36.75 0.75
2 2
= 10.33kN m meter of wall
Solution for Exercise 4
Hence, the maximum moment is 38.88 kN·m/meter of wall.
The section modulus of the sheet piles is thus
M max 38.88
S= = = 2.287 10−4 m3 m of the wall
all 170,000
Exercise 5 Determine the factor of safety against bottom heave for the braced cut described
in Exercise 3. And assume the length of the cut, L = 18 m.
0.2 B cH
5.14c 1 + +
FS = L B
gH +q
B = B 2
B = 2 B
Solution for Exercise 5
From the Equation:
0.2 B cH
5.14c 1 + +
FS = L B
gH +q
B = B 2 =6 2 = 4.243
B = 2 B = 6
0.2 6 30 7
5.14 30 1 + +
FS = 18 4.243 = 1.75
17.5 7 + 0
Tutorial A driven closed-ended pile, circular in cross section, is shown in the Figure bleow.
exercise 1 Calculate the following.
a. The ultimate point load using Meyerhof’s procedure.
b. The ultimate point load using Vesic’s procedure. Take Irr = 50.
c. An approximate ultimate point load on the basis of parts (a) and (b).
d. The ultimate frictional resistance Qs. [Use Eqs. (1) through (2), and take K = 1.4
and δ’ = 0.6’.]
e. The allowable load of the pile (use FS = 4).
L = 15 D (1)
f = f z = L ( 2)
f
=381 mm 136
Example 1
137
Solution for Example 1
Qs = pLf =
( f z =0 + f z =3.05 ) pL + ( f z =3.05 + f z =5.715 ) p
1 ( L − L1 ) +f z =5.715 p ( L − L) z
2 2
0 + 23.38 23.38 + 34.33
= ( 0.381 ) 3.05 + ( 0.381 ) ( 5.715 − 3.05 )
2 2
+ 34.33 ( 0.381 ) (15.24 + 3.05 + 3.05 − 5.715 ) (Part of constant f, smaller ’ used
for safe design, but not realistic)
= 42.68 + 92.04 + 642.05 = 776.77 kN 139
Solution for Example 1
47.946 23.38
z = L’ 70.41 34.33
43.89
(Recommended!)
z z
(higher ’ gives higher Qs, more realistic!)
f = f z = L = K o tan = 1.4 ( 70.41) tan ( 0.6 40 ) = 43.89
f z =5.715 p ( L − L ) = 34.33 ( 0.381 ) ( 3.05 + 3.05 − 5.715 ) + 43.89 ( 0.381 ) 15.24 = 816.96 kN
Qs = 42.68 + 92.04 + 816.96 = 951.68 kN
Qall = 609.8 kN 140
Answer for Assignment 1
1(a)
2 3H 2 3 127
T 0.051
cu ( field ) = = = 54.342 kPa
K 9.385 10−4
PI (%) = LL − PL = 46 − 21 = 25
1(b)
cu ( field ) 54.342
OCR = = 4.693 = 4.283
0 59.55
2
pa = 100 kN/m3
Using ,
Variation of Dr: The Dr increases with the depth following a power function (Attention: x is Dr, y
is depth).
Dr (%)
0 10 20 30 40 50 60 70
0
6
y = 6E-07x4.0813
R² = 0.9961
8
10
3(a)
u0 = w hw = 9.8 (8 − 3) = 49 kN/m2
p0 − u0 280 − 49
KD = = = 2.432
0 95
0.47 0.47
K 2.432
K0 = D − 0.6 = − 0.6 = 0.655
1.5 1.5
3(b)
3(c)