Site Investigation in Foundation Engineering
Site Investigation in Foundation Engineering
NOTES
Owned By:
LOW KINSON
CIV4249: Foundation Engineering 1.1
Topic 1: Site Investigation and Exploration
1. GENERAL ................................................................................................... 2
1.1. Objectives ............................................................................................ 2
1.2. Reasons for investigating subsurface .................................................. 2
1.3. Sources of Information......................................................................... 3
2. SITE EXPLORATION OBJECTIVES ........................................................... 4
7. REFERENCES........................................................................................... 31
1. GENERAL
1.1. Objectives
The practice of site investigations is undoubtedly one of the most critical aspects of
the geotechnical engineers’ work. The broad objective of site investigation is to
accumulate all data for the site in question so that the proposal being considered
may be satisfactorily carried out. It follows that the site in question must be examined
thoroughly so that subsequent analyses and design procedures are appropriate to
and representative of the proposal being considered. The type of information
required will depend on the proposal as indicated below.
(a) Investigations for new works: access, rights of way, other rights and
covenants, location of services, presence of mine workings likely, effect of
works on adjoining properties, ground conditions, records of performance of
adjacent structures, topographic, hydrologic and climatologic data (dams),
environmental considerations such as land use population etc.
(c) Reports on existing works (perhaps on likely adverse effects from new works
or the possibility of enlargement or modifications): all available records
concerning design and performance of original structure, observation of
movements and/or movement records, ground conditions relevant to the
problem.
(a) Experience and general knowledge of the locality. (Basis for scientific
investigation may give warning that otherwise might be overlooked etc.).
(b) Existing records e.g. geological survey maps and memoirs; ordinance survey;
municipal and public authority plans; meteorological records; borehole records
from adjacent significant areas; aerial photographs.
(d) Site Exploration (with field and laboratory testing). This is the investigation of
the subsurface conditions and forms the main part of this course.
Unfortunately its practice is still very much neglected particularly on small
infrastructure projects. Often a site investigation is not ever undertaken. In
these instances problems only become apparent on construction with the
result that delays for remedial measures and/or redesign operations which
increase costs quickly often overtaking by considerable amounts the capital
which would have been required for the site exploration. In other instances,
the site exploration may have been of poor quality with the result that soil
stratification and parameters are not representative of the subsurface
conditions. This situation can lead to two extremes:
More important it is not just sufficient to undertake any old site exploration as this
can lead to further excessive expenditure whether of the form of remedial measures
for overestimation of site suitability or costly foundations for underestimation of site
suitability. The site exploration should be carefully planned, undertaken and
supervised in relation to the proposed developments and the subsurface conditions
encountered.
The precise objectives of the site exploration will be determined by the nature and
scope of the geotechnical engineers' brief. It should be noted however that any
unusual conditions encountered and not necessarily in the brief should also be
reported if considered relevant.
To determine the nature, thickness, dip and variability of the strata to be affected
by the proposed development;
To determine the physical properties of the strata relevant to the design;
To determine the ground water level (or levels) variation and chemical
composition;
To obtain the above data in the most economic and technically accurate way.
This requires flexibility in planning and good knowledge of all available
techniques.
3.1. General
3.2. Variables
Apart from economics, the major considerations governing the magnitude and
emphasis of the exploration are the types of soil profile and the type of development.
(a) Two soil profile extremes are possible (and of course everything between).
(1) Preliminary Explorations to asses the general suitability of the site for the given
development or to select the most suitable site from alternatives. This phase
and its results will enable a preliminary assessment of soil and ground water
characteristics so that subsequent phases of exploration can be planned -
location of boreholes, depth, and likely size of boreholes, testing techniques
and equipment necessary. In areas, which are well documented, this phase of
exploration is not always necessary.
(2) General Exploration to select areas of a given site with the most suitable
ground conditions or simply to investigate the ground conditions of a fixed
location. To obtain the required samples and conduct the appropriate testing.
Usually this phase is sufficient for general construction.
3.4. Economics
effective use of the data in both design and construction (e.g. tender process for
below ground work will often be drastically reduced if contractor is confident of
subsurface conditions).
Exposure of relevant strata for inspection, sampling or testing may be by trial pits or
shafts, boreholes or audits. On the other hand, geological methods of investigation
may provide considerable data without even disturbing the ground surface. It is
recommended that some form of boring is undertaken for correlation of results, and
positive identification. The choice of methods will depend on:
(a) Geological Map (e.g. Australian Geological Survey) should always be the first
source of information. This provides a preliminary indication of what to look for
and what to expect; geological succession and general description; contact
zones; faults dips and strikes.
(b) Field Traverse: On major projects, geologists will be available. The engineer
must understand enough to request a geologist when necessary and to
appreciate the results. On small jobs he must be his own geologist.
Seismic survey methods are based on placing a string of geophones along the
ground surface, and measuring the time taken for a sonic event (such as a sledge
hammer hitting a metal plate, or a small explosive charge) to be recorded at each of
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 1.7
Topic 1: Site Investigation and Exploration
the geophones. The sonic wave first reaching a geophone may have travelled
through the ground directly below the surface, or may have travelled to a deeper
more competent layer with higher sonic wavespeed and been refracted up to the
geophone. Knowing the distance between each geophone and the event allows
geologists to predict a model of the subsurface profile based on sonic velocities.
Some features of the method are:
Cheaper than seismic survey but less accurate. For location of bedrock or faults and
determination of groundwater salinity. Generally better for qualitative rather than
quantitative information.
4.3. Boreholes
Hand augering is used for very simple investigations (e.g. house foundations) or
where access would otherwise be difficult or impossible. The hand auger is simply
a fence-post auger. Different heads are available for different soil conditions.
These augers are often truck-mounted and penetration to very substantial depths
(many 10’s of metres) is possible by progressive addition of more auger sections.
With full flight augers, a worm flight brings the soil to the surface. With other
types, the auger is withdrawn regularly to inspect soil.
Depth from which soil is taken can be reasonably estimated.
'Undisturbed' sampling possible by withdrawing auger. Disturbance below
bottom of hole is very small.
Advancing in sands can be difficult due to borehole collapse. Advancing in
sands below the water table can be impossible.
75mm to 300mm diameter is common for site investigation but up to 1m
sometimes used for inspection bores.
Centre post of auger is hollow tube with plug, which can be removed for
sampling.
Auger (full flight) need not be withdrawn for sampling. Therefore, no hole
support is required.
Good access provided to the bottom of hole with no disturbed material able
to fall in and low boring disturbance ahead of auger.
Percussion boring is a very rapid method of drilling a hole into the ground. The drill
is progressively advanced into the soil/rock, and not withdrawn until the target depth
has been reached. Very disturbed soil and rock cuttings are typically blasted from
the top of the hole from the air pressure used for the percussion.
Rotary wash boring is possibly the most commonly used method of drilling used for
geotechnical site investigations in Australia. It is also the method which causes least
soil disturbance. Principle features of this method are as follows:
Rotary core drilling can be used in association with either augered or rotary wash
bored holes. When competent rock is encountered, these techniques are not able
to advance any further because of the crowd (downward thrust) and torque
limitations of the drill rig. Undisturbed sampling of the rock then requires diamond
coring techniques. With diamond coring, it is possible to establish not only the rock
characteristics, but also evaluate the defects such as joints and fault zones.
Effectively unlimited penetration through the rock is possible, although penetration
rates are very slow (typically 0.5 to 6m per hour depending on rock hardness).
Essential features of this technique are as follows:
Boreholes will usually stay open in clays which are stiff to hard. However, there is
a likelihood of borehole collapse in soft clays or granular soils, so some means of
maintaining the hole open is required. The methods used are either steel casing or
drilling mud.
[Link]. Casing
Casing requires the use of steel tubes in a borehole to support the sides. It can be
installed by driving but is more commonly advanced by hydraulic pushing -
remember - it has to come out again! Casing should never be driven ahead of the
base of the borehole indeed it should be kept as far behind as possible.
Where deep boreholes are to be installed the telescoping of casing may be required
in order to overcome installation resistances. The use of casing has the
disadvantage that it must be carried to site. It takes time to install and of course
remove and it involves capital and maintenance costs. However for shallow holes it
is quite efficient and effective.
A common alternative to casing is drilling fluid. The original drilling fluid was based
on natural bentonite clay, which can be mixed in powder form to the drilling water to
create a higher density suspension. Also, the use of synthetic polymer (e.g., partially
hydrolysed polyacrylamide) as drilling fluid has gained significant appreciation. The
support fluid keeps the borehole stable by applying a hydrostatic fluid pressure
(greater than the pore water pressure) on the base and the side wall. The use of
drilling fluid must be evaluated in terms of its environmental effect on the
groundwater regime. Also, not all fluid types (bentonite/polymer) are suitable for all
ground types.
The means of advancement of the borehole are very important in ensuring that
samples retrieved or tests performed are representative of insitu conditions.
Another key requirement is that the drilling procedure and result of the sampling and
testing be logged accurately and in detail. All holes should preferably be logged by
a professional engineer. Logging should include:
Drilling contractor;
Client;
Location;
Borehole number;
Type of boring;
Date;
Position and elevation of bore;
Size of casing;
Rate of progress recovery;
Water table/s and variations with time;
Upper boundary of each new stratum (by sampling or examination of wash);
Borehole support used;
Field classification;
Depth and full description of all samples (soil and water)
Results of SPT tests or other insitu tests;
Results of any quick tests;
Type of tool used;
Reasons for delay or changes in technique;
Any unusual occurrences;
The following sections describe the mechanics of the SPT test, and some of the
many methods proposed for correlating the results engineering properties.
The SPT is one of the most common forms of in-situ test, and is regularly used to
determine the properties of granular soils. The technique is very crude but has the
major advantage of being cheap and easy to carry out. Also the results it gives
whilst far from rigorous are generally accepted right around the world.
W a te r p o rt
H ead P in 465m m 60m m
35 mm
50 mm
D r iv in g
D rillin g B a ll S p lit shoe 20m m
ro d v a lv e b a rre l T h re a d s
C o u p lin g
(300 mm in total) are added together to give the "N" value for the SPT. By conducting
a reasonable number of these tests in the granular formation a reasonable average
N can be obtained. If the ground at the site is hard/very dense then SPT refusal
(N>100=R) is often encountered where a solid SPT bit is recommended to avoid
damage to the standard hollow-bit. If a solid SPT is done then the collection of
sample from the test depth cannot be performed.
The measured N-value at a given depth for a given soil can be influenced by several
factors (e.g., hammer efficiency=ratio of actual energy to the sampler to input energy
(Wh), borehole diameter, sampling method, rod length). Amongst all the hammer
efficiency is the most influencial; a lower efficiency hammer would produce higher
N value which will produce unrealistically higher values of relevant soil strength
parameters or vice-versa. To avoid this issue, the measured N blow count is often
standardised against an average hammer energy of 60% and the converted N is
termed as N60. Hammer efficiency values for different types of hammers (e.g., donut,
safety) used in different countries can be found in standard foundation engineering
text books including some other minor corrections stated above. In case of no
information on the hammer efficiency is available then an industry recommeded
value or assumption is to be utilised.
The hammer efficiency is to be corrected for all soil types. The N60 obtained from
this process will further subject to dilatancy and overburden corrections as
discussed below.
[Link].1. Dilantancy correction for fine or silty sands below the water table
greater the effective stress the greater is the strength of the material being tested.
Therefore, as the strength increase has only been caused temporarily because of
dilation effects, it must be accounted for by some means. Hence the application of
the above correction to reduce the measured N60 value. As this correction is often
needs a judgement on soil permeability, there are different opinions amongst
practicing engineers.
Note that for this correction to apply the sand being tested must
When correlations were originally made between SPT N values and various
engineering properties a major oversight (error?) was made. This oversight involved
N values obtained at shallow depths. It was found that predictions of soil strength
were unrealistically low and compressibility high. Therefore in order to correct this,
a correction factor which depended on the effective overburden pressure at the
depth at which the test was carried out is introduced. Many relationships were
developed and can be found in the text books. However, one common form is that
if the SPT is conducted at depths having effective vertical stress < pa=100 kPa then
the N value is to be amplified and a reduction is applied for effective stress > 100
kPa. A typical correction factor (CN) as proposed by Liao and Whitmann (1986) is
given below. One may wish to use some other relationships and obtain an average
value as the final correction factor. The corrected N60 is term as (N1)60=CN*N60.
.
Relative density (Dr%): The packing of granular soils can be described by its
relative density. The higher the relative density of the soil the higher is the density.
It has a scale from 0 to 100 with the lower value represent a loose packing (i.e. 0-
15=Very loose soil) and the higher value refers to very dense packing (i.e.,
>85=Very dense soil). The packing has direct relationships with its strength and
deformation properties (e.g., internal friction angle (’), Youngs modulus (E)).
Therefore, estimation of Dr% in terms of N is essential for granular soils as
undistustubed samples cannot be collected for laboratory tests. Any correlations
between Dr% and N must also include the depth or effective stress to reflect on the
real field condition. Some simplified correlations are given below and the readers
are advised to learn many other complex form of the correclations from the text
books.
.
Meyerhof (1957): 𝐷 Where, Dr not in % and applicable
.
to clean and medium to fine sands.
Friction angle (’): The peak friction angle for granular soils can be approximately
correlated with the N value. Some approximate corrections are given below.
SPT tests are generally performed in sands and gravels. However, in saturated
clayey soils, due to the proewater pressure development in clay during the
hammering and its influence on the effective stress makes the measured N value
less reliable. As such apart from the hammer energy correction, no other corrections
to N60 are recommended. Approximate correlations with undrained shear strength
(cu) have been postulated - typically cu ≈ 6*N60.
Hara et al. (1974) suggested an approximate correlation below based on SPT tests
carried out on normally to lightly overconsolidated clays (OCR=1-3) from Japan.
𝑐 .
0.24 ∗ 𝑁
𝑝
See also for a very comprehensive summary of all the factors that can affect SPT
tests, and recommendations for interpretation based on the work of others.
An in situ test, the electric Cone Penetrometer Test (CPT), which can be carried out
in most soils. In this test, a probe is jacked into the ground at a constant penetration
rate. This section briefly describes the test, its advantages and disadvantages, and
its interpretation.
The type and consistency of the soil is determined primarily by the cone resistance
(the greater the cone resistance, the denser the sand or stiffer the clay), and the
ratio between the friction and cone resistances (called the friction ratio,
(Fr(%)=100*(fs/qc)) In general terms, cohesive soils have a high friction ratio (say
>5%), whereas sands have low friction ratios (<2%). Many researchers have
proposed empirical correlations between cone resistance, friction ratio and soil type.
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 1.17
Topic 1: Site Investigation and Exploration
Prominent among these are Robertson et al., (1986) and Robertson (2010). Their
results are published widely (e.g. in standard text books and Journals).
It is quick and cheap to use. Many tests to quite reasonable depth may be
accomplished in a day.
It provides a continuous profile very quickly and therefore is an excellent
method of obtaining detailed information between boreholes or identifying
problem areas.
INTERPRETED
SOIL PROFILE
testing (qc, fs, Fr(%)) are typically
plotted in real-time on paper tape, and
SLEEVE CONE FRICTION may also be stored digitally. This
DEPTH RELATED TO GROUND LEVEL IN METRES
-5
-10
-15
-20
-25
-30
LEGEND:
CLAY SAND
SILT GRAVEL
Numerous CPT interpretations have been proposed for soil classifications. Although
newer charts (Robertson, 1990) have been developed based on normalized
parameters (corrected for overburden and water table, often done after the test), the
most popular and simple chart (Robertson, 2010) based on qc and FR(%) which is
easy to use and can be used in real-time during the CPT is given below. In the figure,
ISBT is the soil behaviour index and pa=atmospheric pressure=100 kPa.
qc/pa
Fr (%)
It is common to infer a value of undrained shear strength (su) from the qc values
and the total overburden stress, o. The following expression is used:
where Nk = the cone factor (commonly varies between 10 and 30, and 11 to 19 for
NC clay and may approach 25 for OC clay).
The appropriate value for Nk for a particular site should be determined by correlation
with undrained triaxial tests on undisturbed clay samples.
The outcome from above correlations for a specific site must be verified with
appropriate laboratory tests (e.g., Laboratory consolidation test).
Kulhawy & Mayne (1990) proposed several correlations for sands as:
A useful correlation between SPT N60 value and qc has been proposed by Kulahwy
& Mayne (1990) as follows:
.
5.44 ∗ 𝐷 Where, D50=Mean particle size (in mm), pa=atmospheric
pressure=100 kPa
The pressuremeter is a probe which can be inserted in the ground and then inflated.
It was originally developed in 1956 as a relatively sophisticated instrument for
measuring the in-situ strength and elastic modulus of soils and rocks. Although the
equipment is relatively expensive, its use in the field can be very economic with
many tests being made in a working day. There are basically two forms of the
pressuremeter. One is installed in a prebored hole while the other is a self boring
device.
The Menard
pressuremeter (G-AM
II) probe which is a
cylindrical metal body
with an inner rubber
membrane and outer
protective sheath
mounted so as to
form three
independent cells.
The middle cell is
inflated with water
and the two guard
cells with gas. The
probe is insereted into
a pre-drilled borehole
at a desired depth for
measurements as shown here (Rocktest, 2020). The deformations are measured
by the central measuring cell where conditions of plane strain exist due to the
presence of the guard cells.
The water flow to the measuring cell and gas to the guard cells are regulated via the
control unit fitted with regulators and valves. The volume changes are read directly
from the graduated sight-tube on the reservoir. A coaxial tubing is used to connect
the probe to the control. The inner tube carries the water to the measuring cell while
the space between the two tube walls carries the gas for the two guard cells. This
method prevents the expansion of the inner tubing which would lead to erroneous
readings of the amount of water injected.
The self-boring
pressuremeter is in
principle the same as
the Menard device, but
in addition it has a
specially designed
cutting head at the
bottom end of the
device. This consists of
a rotating cutter which
cuts the soil and this
soil is in turn brought to
the surface by means
of circulating water.
The principle advantage of this device is that the soil disturbance is minimized. For
a Menard pressuremeter it is necessary to drill an oversized hole into which the
pressuremeter is placed; between drilling and installation, the borehole may relax
or even collapse. The self-boring pressuremeter avoids these problems.
For the Menard pressuremeter there are three main testing procedures according
to the soil type allowed.
(a) For free standing soils in which the sides of the borehole do not cave in the
borehole is advanced by the appropriate method and the probe is lowered to
the testing location. Usually tests are conducted at one metre intervals.
(b) For soft clays and granular soils the borehole may be advanced and kept open
with drilling fluid. The test is conducted as above. According to Menard (1957)
the effect of the bentonite drilling mud is negligible, however, others may
reasonably consider this to be an optimistic view.
(c) For the same soft or granular soils a special casing may be used. It consists of
conventional casing which has narrow longitudinal slots machined into it. The
casing will support the borehole walls but also allows the expansion of the
pressuremeter. However it is obvious that the expansion characteristics of the
slotted casing must first be measured, and appropriate correction applied to the
results. The use of drilling mud is preferred.
Self-boring pressuremeters advance to the test depth from the base of a borehole
drilled to just above that depth.
Once the pressuremeter has been installed at its test location the zero readings are
taken with the pressuremeter sides just touching the soil. The pressure is then
progressively increased and the volume of the cell is measured at each incremental
pressure until the soil fails or the volume of the cavity has expanded by twice of the
original volume (V1). Once done, the proble is deflated and moved to another depth.
The pressure readings are ploteed against the volume volume to obtain effectively
a stress-strain plot, as shown.
There are three basic stages involved in the deformation of the soil as shown in the
figure below.
pl 1
3 before the borehole was made.
2 Stage 2 (po to pf) consists of the elastic
pf
deformation of the soil adjacent to the
measuring cell. This is characterised by a
p straight line plot of the pressure/volume curve.
The elastic modulus E is usually derived from
po 1 V this section of the curve.
The pressuremeter may be used to evaluate the properties of most soils and soft
rocks. However, because the test is relatively quick, the analysis is usually based
on undrained condition for cohesive soils and drained condition for granular soils.
Obviously for soils containing a significant portion of silt sized particles or for rocks
containing significant porous fabric such as fissures and joints, it may not be clear
what type of parameter is being derived by the test.
The soil modulus depends on several factors (Briaud, 2013) one of which is the
strain level as can be seen in the above PMT pressure-strain plot at the cavity wall.
Based on the theory of elasticity and using the slope of the initial straight line
(dashed red line) of the PMT pressure-strain plot (between stage 1 and 2), the initial
shear modulus (G) can be calculated as
The shear modulus can then be related to Young’s modulus (E) for isotropic material
(often this is not the case for geomaerials) as
𝐸
𝐺
2 1 𝜐
It is important to remember that the elastic properties determined by pressuremeter
tests are lateral values, and that unless the soil is isotropic, the derived values will
be different from the vertical values because of the manner in which soils are laid
down over time. As the elastic values needed for footing settlement calculations are
vertical, these values may be inappropriate. Pressuremeter results are more
appropriate for piles which may be laterally loaded.
The strength of soils may also be estimated from the pressuremeter test result.
However, in simple terms, it is only possible to evaluate cu for undrained cohesive
soils and ’ for drained granular soils. When a soil possesses both c’ and ’ (e.g.
rocks or mixed soils) a variety of complex testing procedures must be undertaken.
For cohesive soils tested under undrained condition, the undrained cohesion (cu)
may be approximated based on the cavity expansion theory (Baguelin et al. 1978)
as:
𝑐 Where, 𝑁 1 𝑙𝑛
The value of Np typically varies in the range 5 to 12, with an average of 8.5. Ohya
et al (1982) correlated Ep with field SPT N60 as follows:
.
For Clay: 𝐸 1930 ∗ 𝑁 Where, Ep in kPa
.
For sand: 𝐸 908 ∗ 𝑁 Where, Ep in kPa
For granular soils, fully drained conditions are assumed and ’ is evaluated.
According to the original Gibson and Anderson (1961) analysis, at failure the ratio
of circumferential stress to radial stress is given by
1 sin
N [N.B. Not the N as applies to SPT tests]
r 1 sin
Therefore, by plotting log (V/V) against log p’, a reasonably straight line is obtained
for the plastic region. The slope of this line is ½(1-N) and therefore ′ may be
calculated.
While the above simple treatment of the pressuremeter analysis has been
presented, there are still a great number of factors which influence the test results
and which are not fully understood yet.
The pressuremeter was initially seen as the answer to all geotechnical problems,
but over time it has taken a more sensible position as one of a number of less than
perfect tools which are available to the geotechnical engineer. Pressuremeter tests
are used quite extensively for larger projects where the cost of this type of test can
be justified.
The dilatometer is a blade which can be hydraulically jacked into the base of a
borehole, or pushed with a CPT rig. The blade consists of a flat blade 220mm
(length) x 95mm (width) x 14mm (thickness). A thin, expandable 60mm diameter
steel membrane is located flush at the centre of one side of the plate as shown
below. The membrane can be expanded laterally by 1.1mm against the soil. The
pressure A at “lift-off” of the membrane and the pressure B at which the membrane
expands by [Link] in the soil are recorded. In general, at 200mm intervals or
greater this test can be performed. The recored pressures are corrected and
[Link]/flat-
[Link]
The vane shear test is a method used predominantly in soft clays to estimate in-situ
shear strength. The cruciform section vane is generally pushed into the soil at the
base of the borehole, however, it can also be pushed from the surface with a
protective sheath. The vanes can be of different widths and lengths to allow a broad
range of soil strengths to be tested. The height of the vane is twice its diameter.
Once in place, the vane is rotated at a constant rotational speed of 6°/minute, and
the maximum torque achieved is recorded.
The applied torque, T, to turn the rectangular vane is related to the undrained
strength of a soft clay, su(VST) by
𝑇 𝑠 ℎ
Where, h is the height of the rectangular vane and d is the diameter. For other types
of vane, such as the tapered vane, a separate form of the above equation is used.
Both the maximum torque and the residual torque (i.e., 5-10 revolutions after then
maximum torque) can be measured in the field and the relevant undrained shear
strength values can be calculated. The caulcuated undrained shear strength value
from the vane is found to be higher and requires to be corrected for actual design
purpose as:
. %
𝜆 7.01𝑒 0.57
4.4.7. Field Permeability Tests
Permeability tests include both in-situ permeability tests and permeability tests on
samples retrieved from the boreholes for laboratory testing. Although the in-situ tests
have advantages in that tests encompass the global response of the soil, with soil
fabric relatively undisturbed, interpretation can be difficult when the soil is distinctly
layered, or boundary conditions do not match the simplified model in some other
way.
Permeability tests are not an issue with foundations unless the foundations require
an excavation below potential water-table level.
The previous sections have described a number of in-situ testing methods that are
available to the geotechnical engineer. This is by no means an exhaustive list of
possible tools, but covers the principle tools used for typical site investigations.
Tabulated below is an extensive range of geotechnical tests and their applicability
to plan for a site investigation work.
4.5. Sampling
Samples fall into two broad categories - disturbed and undisturbed. The key features
of each of these categories will be discussed hereafter.
The amount of disturbance varies. The returned wash from the wash boring
process is extremely disturbed; the sample obtained during an SPT is
significantly disturbed.
Handling is important. When obtained samples should be placed
immediately in airtight jars, tins, plastic bags etc. labelled to identify the job
number, locality, bore no. depth and date of sampling and sent to lab.
Samples in which the insitu structure and moisture content have been as far
as possible completely preserved.
No sample truly "undisturbed" because of stress release on sampling.
In stiff and over-consolidated clays, sampling not likely to significantly affect
the 'design properties' of soil.
In soft and sensitive clays, sampling disturbance considerably reduces
measured strength, modulus etc. and increases compressibility.
In sands, undisturbed samples are rarely attempted;
The simplest method of obtaining undisturbed samples is pressing a thin
walled tube into the soil i.e. an open drive sampler.
Sampler should not be driven into the soil but pressed in, preferably at a
constant rate of penetration.
Open drive samplers commonly around 50 mm and 75 mm inside diameter.
L/D ratios vary from about 6 to 18.
To avoid disturbing soil the borehole should be cleaned out before sampling
(e.g. with auger or by jetting). Sample ahead of earlier boring. Also maintain
ground water balance. In soft and sensitive soils the thin walled piston
sampler gives best samples. Tube forced into ground with piston blocking
off end. At required depth, the piston is released from ground level and
sampler advanced over the piston. Enables samples to considerable depths
in soft clays without casing. In variable or very hard soils this cannot be
used. In order to obtain "undisturbed" samples, construction of sampler
must be appropriate. Critical parameter is the area ratio which is defined as:
𝐴 % x 100
The table (ICE, 2012) below provides some useful guidance on the planning stage
of a site exploration work.
It is vitally important that all significant aspects of the site investigation are covered
in relation to the proposed development so that the next phase of the development
(probably the design work) may be satisfactorily carried out. Remember that the
report is generally a legal document and any gross misinterpretation of data may
lead to legal proceedings!
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 1.31
Topic 1: Site Investigation and Exploration
The report will also be a record for use in the future should any defect or
unfavourable condition develop during the life of the development. The report may
also be used for assessing any future extensions or modifications to the original
development.
Should any possible source of concern appear during the investigations this should
also be reported and any recommendation as to its solution should be included.
It is also important that the report contains a clear identification of the site and the
location of the boreholes test pits etc. All borehole logs and tests results should also
be included preferably in standard tabular form. Cross-sections are always a good
way of illustrating the soil stratification.
A good report should include plenty of diagrams and figures which always speak
louder than words.
Remember that the report is the final and probably the only product of the site
investigation and generally since it represents a considerable amount of time energy
and cost it is worth doing well.
7. REFERENCES
Bowles, J.E. (1997). Foundation analysis and design. 5th ed. International ed.
McGraw-Hill, Singapore.
Briaud, J.L. (2010). Ménard Lecture - The pressuremeter test: Expanding its use.
Proceedings of the 18th International Conference on Soil Mechanics and
Geotechnical Engineering, Paris, September 2-6.
Das, B.M. & Sivakugan, N. (2019). Principles of Foundation Engineering. 9th SI Ed.,
Cengage Boston, USA.
Lam, C. & Jefferis, S. A. (2017). Polymer support fluids in Civil Engineering. ICE
publishing, London, UK.
Robertson, P.K. (2010). Soil behaviour type from the CPT: an update. 2nd
International Symposium on Cone Penetration Testing, Huntington Beach, CA,
USA. Volume 2&3: Technical Papers, Session 2: Interpretation, Paper No. 56.
TABLE OF CONTENTS
1. THE MAJOR ISSUES 2
2. FOOTING DESIGN - THE BIG PICTURE 2
2.1 Progressive loading of soils 2
2.2 Bearing Capacity 3
2.3 Settlement 4
2.3.1 Compatibility of settlement and structure 4
2.3.2 Settlement Rate 8
3. FOUNDATION PRESSURE 8
3.1 Gross Pressure 8
3.2 Net Pressure 8
3.2.1 What is Net Pressure? 8
3.2.2 Net Pressure for Bearing - General Principle 9
3.2.3 Net Pressure for Settlement - General Principle 9
3.3 Ultimate Bearing Pressure 10
3.4 Maximum Safe Bearing Pressure 10
3.5 Allowable Bearing Pressure 10
3.6 Design Pressure 11
3.7 A last note - avoiding confusion 11
4. NET PRESSURE-EXAMPLE CALCULATIONS 11
4.1 Short-term excavation in Cohesive Soils- Example 1 11
4.1.1 Net pressure for settlement 12
4.1.2 Net pressure for bearing 13
4.2 Long Term Excavation in Cohesive Soils (both consolidation, and elastic) and Any Term
Excavation in Granular Soils (elastic). 13
4.2.1 Schematic - Long term excavation (cohesive) or short-term (granular) 13
4.2.2 Net pressure for settlement 13
4.2.3 Net pressure for bearing 15
4.2.4 Summary 15
4.3 Surface footing 15
4.3.1 Surface footing 15
4.3.2 Net pressure for settlement 15
4.3.3 Net pressure for bearing 16
4.4 Wide excavations 16
4.4.1 Wide excavations 16
4.4.2 Wide excavation in cohesive soils-short term (consolidation only) 16
4.4.3 Wide excavation in cohesive soils (long term) and any term in granular soils 17
4.5 Impervious excavation under water 17
4.5.1 Impervious excavation 18
4.5.2 Impervious excavation in cohesive soils- Short term (consolidation only) 18
4.5.3 Impervious excavation in cohesive soils (long term), and any term in granular soils 19
4.6 Pervious excavation under water 19
4.6.1 Pervious excavation under water 19
4.6.2 Pervious excavation under water in cohesive soils (short term, consolidation only) 20
4.6.3 Pervious excavation under water in cohesive soils (Long term), and any term in granular
soils 20
4.7 Net Pressure Cases Summary 20
4.7.1 Summary Table 21
4.7.2 General Patterns 21
5. PRACTICAL CONSIDERATIONS 22
6. REFERENCES 22
This introduction gives a brief overview of the major issues that influence the design
of shallow foundations. A major portion deals with various definitions of pressure,
and in particular, the computation of net pressure for bearing capacity and
settlement calculation. That is, it concentrates on the correct calculation of the
disturbing forces. In subsequent sections, the methods for computing the capacity
of the soil, or the compressibility of the soil will be introduced.
“The foundation is that part of a structure which serves exclusively to transmit the
weight of a structure onto the natural ground" (Terzaghi & Peck, 1967).
It follows, therefore, that the foundation is that part of a structure directly in contact
with the soil or rock. i.e., structure = foundation + superstructure. The foundation
must be designed so that ground movements transmitted to the superstructure are
no more that the superstructure can adequately withstand.
This requires:
1. That there is an adequate factor of safety against shear failure of the soil, i.e. the
structure does not break into the soil.
2. That settlement, both total settlement and differential settlement, does not exceed
tolerable limits.
Soils deform under the application of footing loads. This section will consider the
deformation effects as a foundation element is progressively loaded. This will be
related to two fundamental concepts in footing design: bearing capacity and
settlement.
PHASE 1 Movement of B
(elastic) (plastic)
to the footing at the ground
Heave
Movement of A
When P = P1: plastic zones will begin to be formed where the stress concentrations
are high i.e., at the edges of the base. These initially small zones will have reached
their maximum shear resistance and will deform plastically, leading to an increase
in the rate of settlement at point A. At the same time, point B will begin to experience
some upward movement or heave.
As P increases above P1 and towards P2, the plastic zones will increase in size and
this will be reflected in the increasing rate of settlement of point A, and the increasing
rate of heave of point B, as the footing experiences more plastic deformation.
Eventually when P = P2, all the soil underlying the footing (or at least all the soil on
the potential failure plane) will have passed into a plastic state so that all shear
resistance is mobilized with no more to develop. The settlement of point A will have
increased dramatically and point B will have experienced considerable heave. The
footing will have exceeded the bearing capacity of the soil. This example illustrates
the two main criteria in foundation design i.e., the bearing capacity and settlement.
Bearing capacity is
dependent primarily on
soil strength. As can be
seen in the previous
example, the ultimate
strength of the soil is
realised at P2. Obviously
this point must never be
reached in the foundation.
Similarly just before P2 is reached, deformations are still quite large and may be
harmful to the structure. Therefore, the actual design load on the footing must be
well clear of this region of excessive plastic deformation, preferably in the region
below P1. It has been shown that the point P1 is about 1/3 to 1/2 of the ultimate load.
Therefore, in design, a factor of safety is applied to the ultimate load, to give a safe
load. It follows that this value is between 2 and 3 but depends on many variables to
be discussed later.
2.3 Settlement
Historically, cities have been built around in river delta areas because of obvious
access to shipping transport. Unfortunately (or fortunately for geotechnical
engineers), these river deltas are areas of deep deposits of soft, normally
consolidated and compressible alluvial deposits. Although this may have been no
problem in the old days when structures were mainly simple and light low-rise
structures, for modern high-rise structures, large settlements must be
accommodated or alternative costly solutions (e.g. piling-deep foundation) must be
used. Good examples of cities built in delta areas are Venice and Bangkok.
Anyone visiting the Doge’s Palace in Venice can see the effects that large
settlements have had on this building.
= differential
settlement
The rate at which settlement occurs will depend on the soil in question. In granular
(free draining) soils, this settlement will occur immediately on loading of the
foundation. However, with cohesive soils (very slow draining soil), a finite length of
time will elapse between loading and reaching a final settlement (in very fatty clays
it may take up to tens of years to get there).
3. Foundation Pressure
The gross pressure exerted by a foundation on a soil is simply the total force
(including the weight of the foundation) acting at the base of the foundation divided
by the area of the foundation base. The gross applied pressure is independent of
the depth of embedment of the foundation for a footing of the same dimensions. In
general, the gross pressure is not used in design of footings.
The net pressure exerted by a foundation on a soil can be different from the gross
pressure. Since in soil engineering it is usually the net pressure which is important,
care must be taken in its correct estimation.
In general terms, the net pressure is that component of the gross pressure which
for bearing resistance, requires mobilisation of the soil strength to support the
foundation;
for settlement, contributes to the volume change of the soil due to loading and
hence the movement experienced by the foundation.
The net pressure for bearing and settlement are not necessarily the same. The net
pressure for bearing is influenced primarily by foundation geometry. However, the
net pressure for settlement is influenced by four factors which are the geometry, the
type of soil, the time which elapses between initial excavation and final construction,
and the type of settlement (i.e. either consolidation or elastic).
Determining the net pressure for bearing is very simple. All that is needed is to
understand the mode in which bearing failure occurs. Both from the schematic in
the previous example and from the various photographs already shown, it can be
appreciated that soil failure involves a rupture that extends below the footing and
emerges, as heaved soil to the side of footing.
The net pressure for bearing therefore involves a difference in the effective pressure
beneath the loaded footing and the effective pressure at the same depth beside the
footing. If the footing is buried a significant distance below the surface, there will be
a surcharge effect next to the footing which will act to reduce the tendency for this
rupture to occur.
𝑞 𝑞′ 𝑞′
𝜎′ 𝛾′𝐷
Where ’ is the effective unit weight of the soil and Df is the depth of soil adjacent to
the base of the footing. You are strongly advised to understand this general principle
and to use the following particular applications only to appreciate how the general
principle is applied.
Settlement does not involve a rupture to the side of the footing, but rather a vertical
settlement below the footing. Therefore, the effective pressure beside the footing is
irrelevant to the determination of net pressure for settlement.
Settlement is nothing but a change in the arrangement of soil particles below the
footing into a closer packed arrangement. This closer packing results from the
applied footing load which imposes higher stresses on the soils.
This rearrangement may either be due to the immediate elastic response of the soil
(sand or clay), which is fully reversible on removal of the load, or it may be due to
long-term squeezing out of water in clays below the footing, known as consolidation.
effective pressure under the footing after application of the load. It is this pressure
difference which will be responsible for settlement.
𝑞 𝑞′ 𝑞′
𝜎′ 𝛾′𝐻
Where ’ is the effective unit weight of the soil and Ho is the depth of soil at the
footing prior to the settlement calculation. As for the net bearing pressure, you are
strongly advised to understand this general principle and to use the following
particular applications only to appreciate how the general principle is applied.
The ultimate bearing pressure (qu) is the value of loading intensity, in terms of either
net or gross pressures (qu net or qu gross), which will cause the soil to fail in shear
loading to give a bearing capacity failure.
The maximum safe bearing pressure (qs) is the pressure that a soil will safely
support without the risk of shear failure, irrespective of any settlement that may
result.
Please note that since qu net and qu gross can be different, the respective values of
FoS (Factor of Safety) can differ. In soils the net stress is used, and therefore, the
second of the above two relationships is more usual. However, in any report, it
should be made quite clear which of net or gross stresses is being considered.
The allowable bearing pressure (qa) is the maximum allowable load intensity on the
soil. This is not only from a bearing pressure aspect but also from the estimated
amount and rate of settlement, and the ability of the structure to accommodate it.
i.e. 𝒒𝒂 𝒒𝒔
For example, it has been estimated that for a certain structure, qs = 50 kPa.
However, the foundation is to support very delicate and settlement sensitive
machinery. If qs was applied, the estimated settlement may be 50 mm, which is not
acceptable. Therefore, qa would have to be lowered to satisfy the lower value of the
settlement.
It follows that qa is a function of the soil characteristics, the depth and dimensions of
the foundation and the type of structure. It is the maximum design pressure.
The design pressure (qd) is the value of foundation pressure used in the final design.
i.e. 𝒒𝒅 𝒒𝒔 𝒖𝒔𝒖𝒂𝒍𝒍𝒚 𝒒𝒅 𝒒𝒔
3.7 A last note - avoiding confusion
The soil engineer is generally concerned with net values while the structural
engineer is often only interested in the gross values. It follows that considerable care
must be taken to make sure that there are no misunderstandings when
communicating design information between the two.
The above definitions do vary from place to place and in practice must be clarified
with those with whom you deal.
The following examples cover a wide range of possible footing configurations, and
are provided only to demonstrate the application of the general principles for
determination of net pressures for bearing and settlement covered elsewhere in this
document. Do not attempt to learn each case separately by rote; the point is to
understand the principles, after which you can quickly compute any case by yourself
especially in the exam or in the design!
Consider an excavation of depth H made into a cohesive soil in which the time
between excavation and construction is short. Assume the clay movement is
controlled by consolidation only i.e., the elastic component of clay rebound and
settlement are ignored.
Stage (a): Imagine a small element of soil A, located immediately below the bottom
of a future excavation. This element will have a depth H below the ground level.
Stage (b): Immediately after the excavation is formed, since insufficient time has
elapsed to permit the clay to swell (or experience negative consolidation), element
A will still be at the same depth.
Stage (c): If the foundation is now placed and loaded to H (the pre-existing soil
stress) then the stress at stage (a) are simply being reimposed and the element A
will still be at the same depth. It follows that for these first three stages, the element
will not have moved.
Stage (d): If we now impose the balance of the gross pressure (q1 - H) to fully load
the foundation, the soil beneath the foundation will consolidate, and element A will
move down.
From the above, it is clear that the consolidation settlement of the foundation is
controlled by the pressure component (q1 - H), which is the net pressure, and not
by the gross pressure, qgross. In other words,
𝑞 𝑞 𝛾𝐻 𝑞 𝛾𝐻
The net pressure for consolidation settlement could be termed as the difference
between the eventual final effective pressure (f’) and the effective pressure which
existed immediately before foundation placement (o’).
From the above, since insufficient time has elapsed for the porewater pressure to
stabilize and swelling to take place, the initial effective pressure is given by
𝜎′ 𝛾𝐻
The final effective pressure is given by,
𝜎′ 𝑞
Now the net pressure for bearing resistance will be considered. Whether the soil has
sufficient time to rebound on excavation and then settle on reimposition of the soil
load (Stages (a) to (c)) is unimportant in this case. Any rebound would only change
the depth H by a few millimetres. Since H is in terms of meters, the influence is
negligible. Supposing we have got as far as stage (c) with H applied to the soil.
Since this quantity of pressure is exactly balanced by the soil alongside the
excavation, there is no resultant force trying to displace the foundation. It is only in
stage (d) after the load is increased beyond H towards q1, that strength must be
mobilised. Therefore, the net pressure for bearing is given by,
𝑞 𝑞 𝛾𝐻 𝑞 𝛾𝐻
Which, in this case, is the same as the net pressure for settlement.
Alternatively, the net pressure for bearing could be estimated as the difference
between the eventual final effective pressure immediately beneath the foundation
(f’) and the effective pressure which exists immediately beside the foundation (b’).
Stage (a) represents the at-rest in-situ equilibrium state prior to the excavation.
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2A.14
Topic 2A: Introduction to Shallow Footings
Stage (b): For cohesive soils- there will now be sufficient time after excavation for
the base and hence element A to swell upwards and therefore raise its elevation.
For granular soils, this movement is immediate ie., this swelling would occur
immediately after excavation. Cohesive soils also exhibit some elastic deformation,
and therefore, this would also occur immediately on excavation. Therefore, it is
assumed that time is sufficient to let both elastic and consolidation to fully occur for
all soils.
Stage (c): The foundation is now placed on the base of the excavation at its new
increased elevation and loaded to H. This will cause the element A to return to its
original position.
Stage (d): When the balance of the gross pressure (q1 - H), is applied, further
settlement will occur. With cohesive soils, for the application of both H and (q1 -
H), consolidation settlement will take time whereas elastic settlement for cohesive
and granular soils will be immediate.
The variation of elevation of Point A with time is shown in the figure below.
i.e. 𝑞 𝑞 𝑞
Using the general principle, the effective pressure which existed immediately before
foundation placement (at Stage (c)) was zero since swelling (or elastic rebound) has
occurred. Therefore,
𝜎′ 0
As before, the eventual final effective pressure, f’, is given by,
𝜎′ 𝑞
𝑞 𝜎′ 𝜎 𝑞
As was argued for the short-term excavation, the fact that rebound and settlement
makes negligible difference to the value of H and therefore, the net pressure for
bearing is the same as in for the short-term excavation case.
𝑞 𝜎′ 𝜎 𝑞 𝛾𝐻
4.2.4 Summary
On the basis of the above discussion it would appear that the net pressure for
consolidation settlement for clays depends on the time lapse between excavation
and construction. If the time is short, then qnet for consolidation is less than qgross by
a value equal to H.
i.e., 𝑞 𝑞 𝛾𝐻
If the time between excavation and construction for clays is long enough to allow
swelling, or for elastic settlements for both clays and granular soils, then qnet is the
same as qgross.
i.e., 𝑞 𝑞
For the net pressures for bearing, it does not matter whether this time is long or
short, but is a function of geometry.
For net pressures for settlement, since no excavation takes place, o’ and f’ are
quite obviously the same for both short and long term cases. Therefore,
𝜎′ 0 𝑎𝑛𝑑 𝜎′ 𝑞
𝑖. 𝑒. , 𝑞 𝑞
For the net pressure for bearing, since time does not matter,
𝜎′ 0 𝑎𝑛𝑑 𝜎′ 𝑞
𝑖. 𝑒. , 𝑞 𝑞
4.4 Wide excavations
An isolated footing within a wide excavation is typical of the situation where isolated
column loads bear on the floor of a wide basement excavation. This section deals
with the net pressures for settlement and bearing for this special case.
q3
H
The case of a wide excavation in cohesive soils differs from the narrow excavation
only in terms of net pressure for bearing.
Since no swelling (or negative consolidation) takes place, the net pressure for
settlement is given by,
𝜎′ 𝛾𝐻 𝑎𝑛𝑑 𝜎′ 𝑞
𝑖. 𝑒. , 𝑞 𝑞
[Link] Net pressure for bearing
𝜎′ 0 𝑎𝑛𝑑 𝜎′ 𝑞
𝑖. 𝑒. , 𝑞 𝑞
[Link] Summary
It would appear in this case that even though the excavation is at depth, the term H
does not influence the net pressure for bearing. The reason for this is that the soil
resisting the bearing pressure exists only to the same level as the base of the
foundation. For narrow
Additional
excavations, where the
resistance soil is immediately
adjacent to the footing
right up to the ground
Failure surface, there is an
surface additional resistance (or
With overburden
to resist
No overburden
to resist weight) of H which must
be mobilised.
4.4.3 Wide excavation in cohesive soils (long term) and any term in granular
soils
The cases of long term wide excavations in cohesive soils, or any wide excavation
for granular soils are entirely equivalent with respect to net bearing pressures for
settlement and bearing. In both cases, the net pressure and gross pressure are
equal.
𝜎′ 0 𝑎𝑛𝑑 𝜎′ 𝑞
𝑖. 𝑒. , 𝑞 𝑞
[Link] Net pressure for bearing
𝜎′ 0 𝑎𝑛𝑑 𝜎′ 𝑞
𝑖. 𝑒. , 𝑞 𝑞
4.5 Impervious excavation under water
The situation corresponds to an excavation in which water has not had time to enter
the excavation, or is prevented from entering the excavation by some impervious
lining (which would not be common), or most likely, the excavation is being pumped
dry from a sump.
The important element to remember with these cases is that the excavation does
not contain water, but the base of the footing experiences a buoyant force from the
water table outside the footing excavation.
Since no swelling can take place the initial stress is equal to the in-situ overburden
stress (we are now dealing with a buoyant effective pressure),
𝜎′ 𝐻 𝛾 𝛾 𝐻𝛾′
Where, =bulk unit weight of soil, w=unit weight of water and ’=effective unit weight
of soil.
The final gross pressure is less the buoyant uplift (remember Archimedes?)
𝜎′ 𝑞 𝛾 𝐻
𝑞 𝜎′𝑓 𝜎′ 𝑜 𝑞4 𝛾𝑤 𝐻 𝐻 𝛾 𝛾 𝑞4 𝛾𝐻
This is the same result as for short term excavations in clay where no water was
involved. Note that despite this difference, the total overburden pressure (H)
appears in the above relationship and not the effective overburden pressure.
𝜎′ 𝑞 𝛾 𝐻 and 𝜎′ 𝐻 𝛾 𝛾
𝑞 𝜎′ 𝜎 𝑞 𝛾 𝐻 𝐻 𝛾 𝛾 𝑞 𝛾𝐻
4.5.3 Impervious excavation in cohesive soils (long term), and any term in
granular soils
As for the short-term cases, the important element to remember here is that the
excavation does not contain water, but the base of the footing experiences a
buoyant force from the water table outside the footing excavation. Whether in
cohesive soils where full rebound may take a substantial time, or for sands where
rebound is elastic and instantaneous, the net pressure will be identical.
For settlement, since swelling can take place as well as any elastic rebound, the net
pressure is derived from
𝜎′ 0 𝑎𝑛𝑑 𝜎′ 𝑞 𝛾 𝐻
𝑞 𝜎′𝑓 𝜎′ 𝑜 𝑞4 𝛾𝑤 𝐻
𝜎′ 𝑞 𝛾 𝐻 and 𝜎′ 𝐻 𝛾 𝛾
𝑞 𝜎′ 𝜎 𝑞 𝛾 𝐻 𝐻 𝛾 𝛾 𝑞 𝛾𝐻
The important element to remember with these cases is that the excavation contains
water at the same height as the water table adjacent to the excavation. I can’t think
why one would want to perform construction in this manner except in the case for a
deep pile foundation!
As for all other similar cases treated previously, it is assumed that construction of
the footing is too rapid not to allow any significant heave to have occurred. Any
elastic component of the settlement is ignored.
Note that there is no buoyant force component this time, if it is assumed that the
foundation is of small volume.
𝑞 𝜎′ 𝜎 𝑞 𝐻 𝛾 𝛾 𝑞 𝛾′𝐻
𝑞 𝜎′ 𝜎 𝑞 𝐻 𝛾 𝛾 𝑞 𝛾′𝐻
4.6.3 Pervious excavation under water in cohesive soils (Long term), and
any term in granular soils
Sufficient time is allowed for the soil to fully rebound - a strange requirement for an
excavation in clay under water. I would think substantial softening would occur.
𝑞 𝜎′𝑓 𝜎′ 𝑜 𝑞5
𝑞 𝜎′ 𝜎 𝑞 𝐻 𝛾 𝛾 𝑞 𝛾′𝐻
In the above cases, 5 different conditions were examined and the results are
summarised in the net pressure cases summary table below. It is possible to
determine general patterns for net pressures from this table.
The following table demonstrates that each particular case should be treated on its
own merits and the respective qnet values should be estimated on the basis of the
geometry of the problem, the soil type involved, the times of construction and the
type of settlement involved.
q
1
(i) q1 - H q1 q1 - H
q2
q2 q2 q2
(ii)
q3
(iii) q3 - H q3 q3
q4
q4 - H q4 - w H q4 - H
(iv)
q5
(v) q5 - 'H q5 q5 - 'H
Despite the apparent complexity of the range of results in the net pressure cases
summary table, they can be very much simplified. For example, the condition of
pervious footings under water is most unlikely as foundations are not surrounded by
water, but are generally impervious and extend to the surface in one way or another.
(a) For consolidation settlements with short term excavations (which is the usual
situation with clays):
𝑞 𝑞 𝛾𝐻
(b) For elastic settlements for both clays and sands (and the more unusual
consolidation settlements with long term excavations):
𝑞 𝑞
(c) For bearing pressure when the excavation comes to the side of the foundation:
𝑞 𝑞 𝛾𝐻
(d) For bearing pressure, when the excavation is significantly wider than the
foundation:
𝑞 𝑞
It is suggested therefore that the results are not memorised but the principles
involved are thoroughly understood.
5. Practical Considerations
1. Footings are not weightless - their weight should be taken into account, and it is
unconservative to assume that they simply replace the weight of the removed
soil.
2. Footing excavations are not left open - they are backfilled with either the
excavated material or select backfill. This material adds to the footing pressure.
3. Elastic rebounds of sands or clays occur during excavation and are never
observed. Any rebound is (unconsciously) removed by trimming the excavation
to depth.
4. The elastic settlements due to the footing weight, or any part of the structure
which is cast or erected prior to connection with other structural members will
never be realized. Consolidation settlements that occur prior to connection will
also not be relevant to the structure.
6. References
Terzaghi, K., Peck, R.B. & Mesri, G. (1996). Soil Mechanics in Engineering Practice,
3rd ed. John Wiley & Sons.
TABLE OF CONTENTS
1. INTRODUCTION 2
1.1 Analytical methods 2
1.2 Historical/Experience 2
1.3 Loading tests 3
2. ANALYTICAL METHODS ( = 0 CONDITION) 3
2.1 Lower Bound Solution 3
2.2 Upper Bound Solution 4
2.3 Prandtl Punch 4
2.3.1 Schematic of Prandtl Smooth Punch 5
2.3.2 Analysis of Prandtl Smooth Punch 5
2.3.3 Rough Punch 6
2.3.4 Rough Punch with > 0 6
3. APPROXIMATE ANALYSES (C-- CASES OF GENERAL SOILS) 6
3.1 Terzaghi’s Approximate Analysis for Strip Footing 7
3.1.1 Application of Terzaghi Equation to Other Foundation Shapes 9
3.2 Meyerhof's Approximate Analysis 9
3.3 Analyses by Hansen and Vesic 10
3.4 Other Analyses 11
3.5 So which Equation do you use? 11
3.6 Deep Foundations 12
4. RELEVANT SOIL PARAMETERS 12
4.1 Granular Soils 12
4.1.1 What happens when the footing is loaded? 12
4.1.2 Strength parameters for granular soils 13
4.2 Cohesive Soils 14
4.2.1 What do these undrained and effective parameters mean? 14
4.2.2 Do I use undrained or effective? 14
4.3 Other forms of Loading 15
4.3.1 Eccentric Loading 16
5. ULTIMATE BEARING CAPACITY OF LAYERED SOILS 17
5.1 Stronger soil over weaker soil (c’-' soil) 18
5.2 Weaker soil over stronger soil (c’-' soil) 19
5.3 Rigid base at shallow depth 19
6. ULTIMATE LIMIT STATE DESIGN METHOD (ULS) 20
7. REFERENCES 21
There are many different ways of applying analytical methods but all essentially
involve the use of soil mechanics theory to produce a relationship between the ultimate
failure stress the relevant properties of the soil and the geometry of the foundation.
The very complex real situation can only be resolved if a number of simplifying
assumptions are made. It follows that the relevancy of any analytical method is
governed by the extent to which these assumptions are reasonable for the situation
examined. What assumptions are acceptable in one case may not be acceptable to
another case. Therefore, some care must be taken in the selection of the analytical
method used for a particular situation. Since modern engineering practice requires a
reasonably accurate estimation of foundation types and costs long before construction
takes place, the use of analytical methods is by far the most commonly encountered
approach.
1.2 Historical/Experience
Historically the system used for foundation design was based on experience of similar
conditions. The practice for centuries was "if it worked before it should work again. If
it did fall down, make the next foundation wider or deeper". Unfortunately such an
attitude cannot be used quite so literally in today's economic and political climate. The
size of buildings today are very much bigger and more sophisticated and labour is far
from cheap. It follows therefore that foundation failures must be avoided but at the
same time not over-designed. Despite this, experience does play a very important role
in soil engineering. Most soil engineers on the basis of initial site inspections and the
qualitative results of the first few boreholes can make a reasonable estimate of likely
foundation types. However such estimations must ultimately be substantiated by more
detailed analyses.
This method basically measures the ultimate bearing capacity directly. While such
tests are often performed on difficult soils (e.g. highly variable) they are very costly in
terms of time and effort. This is particularly so with large foundations. Scale models
can be used, but there are a number of scale errors which can be introduced making
the exercise somewhat less attractive. It follows that before such efforts are
undertaken, a reasonable estimation of bearing capacity must be made (so that the
loading system can be designed) and these can be best achieved using analytical
methods. Loading tests are normally used only where overall project capital costs are
relatively high.
This section on analytical solutions is provided for background and to expand your
knowlede rather than to give you practical solutions to everyday foundation design
problems. Analytical solutions are generally valid for special or particular cases (e.g.
= 0 in this case), and do not provide a means for general solutions. Nevertheless,
some of the more general approaches developed later are based on a combination of
analytical solutions for different special cases.
Analytical solutions can be either upper bound or lower bound. If the upper and lower
bounds of a problem can be determined, the "exact" answer must lie in between.
If an equilibrium distribution of stresses can be found which balances the applied load
and nowhere violates the yield criterion the soil mass will not fail or will be just at the
point of failure. Consider an approximate lower
qu = q B
bound solution for a long strip foundation as
1
q q
2
=0
shown in the figure right.
A A s
ie SOIL IS WEIGHTLESS
In zone 2 the possible principal stresses are
ZONE 1 ZONE 2 plotted.
q q 𝑞 2𝑐
0 0 qA
qA qB
2 Cu 2 Cu
In zone 1 which must be compatible along the
qB 0 boundary between the two zones,
qA qA qA qA STRESSES ON
1 2 SOIL ELEMENT
qB 0
OF EACH ZONE
𝑞 𝑞 𝑞 2𝑐 4𝑐
However although the above is elastically in equilibrium the system is obviously not
kinematically correct. This is therefore a lower bound solution.
The soil mass will collapse if there is any compatible pattern of plastic deformation for
which the rate of work of the external loads exceeds the rate of energy dissipation.
qu
This basically means that if a
plausible collapse mechanism is X
Cu
By taking moments about point X
𝑅
𝑞 𝑅 𝜋𝑅𝑐 𝑅
2
𝑞 2𝜋𝑅𝑐 6.28𝑐
Since we have found the upper and lower bounds the correct answer will be
somewhere in between.
𝒊. 𝒆. , 𝟒𝒄𝒖 𝒒𝒖 𝟔. 𝟐𝟖𝒄𝒖
It is generally this upper bound method which is used, for example slope slability
analysis and a factor of safety applied against total collapse. In this way as long as
there is an adequate factor of safety it does not matter too much about local
overstressing.
It should be noted further that the above methods give no indication of foundation
movements. These however can be calculated by other methods. Also, the upper
bound solution will depend on the correct choice of failure shape. The above gives a
high value of qu but could be improved by choosing a more kinematically admissible
failure surface.
The Prandtl smooth punch is an analytical solution which was developed for the
theoretical solution of punching into a weightless material. The analysis is developed
for various conditions of the punch base, and material strength as discussed below.
90 a 𝑞 𝑞 𝜋𝑐
0
qA
(refer to previous example by taking moments
for Zone II)
In Zone I
2𝑐
𝑞 𝑞 𝑞 𝑐 𝑞 𝜋𝑐 𝑐 2 𝜋 𝑐 5.14𝑐
2
(Coincidentally the above lies between 4.0cu and 6.28cu)
However, again, this is by considering a weightless soil and having a smooth base on
the punch which is clearly not the real situation.
3𝜋
𝑞 1 𝑐 5.7𝑐𝑢
2
2.3.4 Rough Punch with > 0
Let us consider the above rough punch case extended to 0, although the soil is still
weightless. By doing this, the solution can be made more general.
qu
The wedge immediately
beneath the smooth footing PASSIVE
ACTIVE
(Zone I) although still an 45 + /2
45+/2 qA
III 45 - /2
45-/2
I
active wedge changes shape qB
because of a finite value of . II
Log spiral
The inclination of the sides of
the wedge increase from 45
to (45 + /2). Similarly the angle formed by the passive wedge (Zone III) decreases
from 45 to (45 - /2). It follows that the interconnecting shear zone (Zone II) can no
longer be a radial fan but must become a log spiral which complicates the analysis
considerably.
When the weight is introduced to the soil the mathematics becomes extremely
complex. Indeed for a soil which possesses cohesion (c), an angle of friction () and
weight (), a rigorous mathematical solution has yet to be established. It follows that
empirical solutions must be used to provide reasonable (in terms of complexity and
ease of use) methods of analysis but still providing sufficiently accurate and realistic
solutions.
Various solutions have been proposed and these will be discussed in this section
below.
Terzaghi (1943) first extended the rigorous but limited analyses of Prandtl to an
approximate but more general case of a non-surface footing (D > 0) on a c, , soil.
B
The frictional nature of qu
the soil is clearly SURCHARGE D
relevant to the strength
I
contribution of each of
III 45 - /2
, D and c whereas
none of these latter II
parameters is Log spiral
significantly interactive.
Consequently
Terzaghi's approximation was based on the superposition of three simplified situations
as follows:
𝛾𝐵
𝑞 𝑁 .𝑓 𝛾 𝑁
2
For a soil with and D but c = = 0:
𝑞 𝑁 .𝑓 𝑐 𝑐𝑁
These situations are specified by the above N parameters with three different failure
surfaces. The superposition is not theoretically valid. However by numerical analyses
Terzaghi was able to provide a combined solution with a little conservatism.
Terzaghi further considered a perfectly rough footing (i.e. allows shear at interface)
and neglected the shear strength contribution from soil above the footing level (i.e. the
soil above the footing level acts only as a surcharge).
Based on Terzaghi's analysis, the bearing capacity for a strip footing on a dry soil is
given by:
𝟏
𝒒𝒖 𝒄𝑵𝒄 𝒑𝒐 𝑵𝒒 𝜸𝑩𝑵𝜸
𝟐
It should be noted that in the figure above,
This equation is for shallow footings and although originally for continuous strip
footings on the basis of field and experimental evidence it has been extended to cover
rectangular, square and circular footings. Generally, the solution becomes more
conservative with increased depth.
The generalised Terzaghi’s equation for strip footing thus becomes as follows to take
into account the water table should it be present.
𝟏
𝒒𝒖 𝒈𝒓𝒐𝒔𝒔 𝒄𝑵𝒄 𝒑 𝒐 𝑵𝒒 𝜸′𝑩𝑵𝜸
𝟐
𝟏
𝒐𝒓, 𝒒𝒖 𝒏𝒆𝒕 𝒒𝒖 𝒈𝒓𝒐𝒔𝒔 𝒑𝒐 𝒄𝑵𝒄 𝒑 𝒐
𝑵𝒒 𝟏 𝜸′𝑩𝑵𝜸
𝟐
Nc, Nq, and N are known as the Terzaghi bearing capacity factors and are a
function of only.
Strength parameters c and refer to the soil below founding level (i.e. in the
failure zone)
𝟏
’ term in ( 𝜸′𝑩𝑵𝜸 is the effective unit weight of soil below founding depth (i.e.
𝟐
in the failure zone). If the water table is below B depth from the footing base
then ’= otherwise a weighted average ’ is to be used if water table is between
the base and B depth below footing. If the water table is above base then always
use ’.
’ term in (𝒑 𝒐 𝜸′𝑫 is the effective weight of soil above the founding depth
(i.e., the surcharge). If the water table is above the founding depth then the
effective stress at founding depth is to be calculated using the effective unit of
the submerged zone and total unit weight of the dry zone. If the water table is
below the base then use ’= for surcharge.
The N factors (refer to standard text books (e.g., Das & Sivakugan, 2019) for
their values against ) are for long strip footings but may be modified for other
shapes.
The net pressure applied to the underlying soil is the pressure in excess of the
overburden pressure (p’o=’D). As the founding depth (D) for a shallow footing is often
less than 2m the difference between the net and gross ultimate bearing capacities is
expected to be minimal. However, for deep excavations, this may not be the case. In
the above net equation, the unit weight of soil and concrete is assumed to be the same
(though this may not be the reality!).
This solution is strictly applicable to shallow foundations where D/B < 1.0 but in
practice also applied to foundations with D/B > 1.0.
Exact solutions for shape corrections have been obtained. Terzaghi modified his
equations as follows to give suitable approximate solutions for footings other than strip
footings.
Because of the additional resistance to failure at the sides of the footings, the bearing
capacity is increased above that of the strip footing.
[Link] Square
𝒒𝒖 𝒏𝒆𝒕 𝟏. 𝟑𝒄𝑵𝒄 𝒑 𝒐
𝑵𝒒 𝟏 𝟎. 𝟒𝜸′𝑩𝑵𝜸
[Link] Circular
𝒒𝒖 𝒏𝒆𝒕 𝟏. 𝟑𝒄𝑵𝒄 𝒑 𝒐
𝑵𝒒 𝟏 𝟎. 𝟑𝜸′𝑩𝑵𝜸
qu IV
III
I
45 + /2
II
L og spiral
Meyerhof's equation is of the same general form as that of Terzaghi, but includes so-
called shape, depth and load inclination factors as multipliers to the bearing capacity
factors Nc, Nq and N. Note, however that the numerical values of these factors are
different from those proposed by Terzaghi. Meyerhof proposed factors based on
theoretical solutions for cohesion only, friction only and surcharge only solutions,
which were then combined to give a general solution.
Meyerhof assumed the slightly different failure pattern which is shown above.
𝑁 𝑁 1 𝑐𝑜𝑡𝜙′
𝑁 𝑁 1 tan 1.4𝜙
Where, the s values are footing shape factors, the d values are footing depth factors
and the i values are load inclination factors. Meyerhof’s recommendations for these
factors can be found in texts such as Das & Sivakugan (2019) and Bowles, J.E. (1996).
Meyerhof's solution differs from Terzaghi's in that the strength of the soil above the
base of the foundation is considered to assist more than just by providing a surcharge.
This gives a more realistic failure pattern, particularly for deeper foundations.
However, while the Meyerhof analysis does generally give far more realistic
predictions of ultimate bearing capacity for D/B > 1 the performance or design of
foundations is governed by both bearing capacity and settlement. Therefore, the
bearing pressures predicted by Meyerhof's analysis, even after the application of a
factor of safety may result in a footing design with excessive settlements.
Hansen and Vesic were both researchers that developed modified bearing capacity
recommendations based on the original method proposal by Meyerhof. The form of
the generalized bearing capacity equations proposed by Hansen and Vesic are
virtually identical to that of Meyerhof, and their Nc and Nq factors are exactly the same.
However, their N factors and their shape, depth and inclination factors differ. Again,
refer to texts such as Das & Sivakugan (2019) and Bowles, J.E. (1996).
In addition to the other factors defined for the Meyerhof equation, the b values are
factors for sloping footing base and the g values are additional factors for sloping
ground.
There are several other analyses in technical papers on bearing capacity. However,
these are often extremely complex and have not been verified sufficiently by
experimental evidence. However until a more adaptable solution is obtained, the
Terzaghi solution will probably be continued for some time as it is applicable to virtually
all soils based on standard accepted testing techniques.
The following facts should be borne in mind when considering the development of
newer solutions.
Following on from above, engineers tend to use methods which are simple to
apply provided the results are acceptable. Any improved method is likely to be
more complex and therefore less easily applied. Also confidence in its
application would be lacking and therefore there would a strong tendency to
use existing methods which are well proven even though they may be
conservative.
That is to assume that there is a correct answer, which there is not. As in all
engineering, our mathematical models are just convenient ways of (hopefully)
estimating the actual behaviour of complex systems. Because of the indeterminate
nature of soils, the actual behaviour is probably more difficult to determine than for
most other (structural) materials.
For what it’s worth, Bowles (1996) has the following opinion:
For D/B > about 5, the foundation is considered as deep and the equations discussed
above are inadequate. This problem will be discussed later in the course when we
examine deep foundations.
Granular soils behave fundamentally differently from cohesive soils because of their
more open and free-draining structure. This section will deal with the type of strength
parameters which are appropriate, and the way in which these parameters are
obtained.
In order to understand what strength parameters to use for a sand, consider what
happens when the footing is loaded. As the footing is loaded, either progressively
during construction or in service as live or wind loads need to be resisted, the stresses
below the footing will increase. This stress increase will tend to compress the sand
structure, moving particles closer together and reducing the voids between particles.
This compression of the structure happens instantaneously. If there is any water in the
voids or pores between the sand particles, the structure is open enough that this can
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2B.13
Topic 2B: Bearing Capacity of Shallow Footings
be squeezed out immediately to accommodate the new denser packing of the grains.
So, no load is transferred to the water and pore pressures therefore no change of pore
pressure (u. All load is carried by the sand particles, and the closer packing means
increased inter-particle stresses (i.e., increased effective stress though =’ here).
Because of the increased inter-particle stresses, the sand becomes harder to shear.
This is because in order to shear, the particles have to ride up and over each other,
which they are being prevented from doing. The shear strength therefore increases
with the applied footing stresses.
This increase of strength, , with increase of normal stress, , is reflected in the
frictional part of the shear strength equation below.
Δ𝜏 Δ𝜎𝑡𝑎𝑛𝜙′
Because of their free draining nature there is normally no excess porewater pressures
present under a foundation placed on granular soils at any time. Therefore, it is the
drained or effective strength properties which are relevant to the calculation of bearing
capacity on sands and gravels. These properties consist of a finite value of ' with c' =
0. (It may be of interest to note than for cuts in granular soils the same effective
strength properties also apply because excess porewater pressures are not involved).
The main exception to the above general rule is when there are major dynamic or
shock loadings (e.g. earthquakes) applied to a granular soil. Load under these
conditions is applied so rapidly that there is insufficient time for drainage to occur and
liquefaction may occur as effective stresses are thrown onto the porewater pressure.
For such a situation it may be appropriate to deal with the problem in terms of
undrained strength parameters. However this is not a problem that will be considered
here.
In order to obtain strength parameters for sands to use in bearing capacity (or any
other) analysis, it is necessary to perform in-situ tests. Can you guess why?
Contrary to the situation which exists with cohesive soils it is very difficult to obtain the
angle of internal friction of granular soils directly for use in bearing capacity analyses.
As will be discussed at length later it is virtually impossible to obtain an undisturbed
sample for laboratory testing. A disturbed sample is of little use as its properties will
not be the same as the soil as it was in its in-situ undisturbed state. Therefore, in-situ
tests have to be relied on. Unfortunately in the great majority of these tests with the
exception of perhaps the pressuremeter test the angle of internal friction cannot be
measured directly but must be inferred from correlations with the measurement made
by the test technique.
Discussions about the Standard Penetration Test (SPT) and Cone Penetration Test
(CPT) can be found in the notes on Site Investigation, as well as numerous other texts.
Essentially sands and clays behave in a very similar manner. They both respond to
applied load by compressing, and adopting a more closely-packed and therefore
stronger structure. The main difference is that because of their open structure, sands
can do this immediately, whereas clays can take a long time to adopt the new more
compact structure.
Let’s assume that the clays have any amount of time to rearrange and adopt the
smaller structure. Remember the term consolidation? Just as for sands, the strength
of the clays is (virtually) dependent only on the applied stresses i.e., it is predominantly
frictional in behaviour, and is proportional to the effective stress. Why is the effective
stress? Because it is the total stress less the soils own buoyancy which determines
the level of stress between particles. So for clays in the long term,
Δ𝜏 Δ𝜎𝑡𝑎𝑛𝜙′
Note: just as for sands, but mind you, quite often there is a small cohesive component
as well - the so-called “true cohesion”, c’.
What does happen in the meantime as this clay is consolidating? Initially as the load
is applied, if the clay is below the water table, there is no way that the clay particles
can get closer together because there is all the pore water in between. As the clay
particles can’t get closer (i.e., no change in effective stress), there can’t be any
increase in strength. But the extra load applied must be taken by something and in this
case, it is the water which takes the extra load, and the increased load in the water
causes the water pressure (or pore pressure) to increase.
So, here we have a situation where no matter how much load we apply to the clay, in
the short-term there can be no increase in clay strength - that is the strength is stress
independent. What is the stress-independent part of the strength equation? - see if
you were right.
Now, assuming that the footing is loaded, and the pore pressure has increased, what
is going to happen? There will be a pore pressure difference below the footing with the
surrounding pore water, which will tend to drive water away from the footing area. As
this water is driven out (over time), the result will be consolidation of the clay and a
strength increase that will depend on the stress applied and the time allowed for
dissipation of pressures. However, in the long-term, as described above, the strength
of the clay will be simply stress-dependent, or frictional.
In general terms when construction takes place on cohesive soils the undrained
parameters (i.e. the undrained cohesion cu and the undrained angle of internal friction
u which is conventionally taken as zero) are usually used. The reason for this is that
the time taken for construction can be very roughly of the same order of for the
When a structure is placed on a clay consolidation settlement can go for several years
and while this settlement will cause a strength gain in the clay it may not be realised
until long after the structure has been put into service. While it is the larger drained or
effective strength parameters which will eventually apply once drainage has been
completed they would not be the more critical. Nevertheless, there are situations
where some allowance can be made for strength gain during construction. However
this would require a great deal of data on rates of porewater pressure dissipation,
consolidation, settlement and strength characteristics which for real soils is difficult to
determine.
It might be interesting to compare this situation with the construction of cuts in clay
materials where because stresses are removed (rather than applied as is the case
with structures) the most critical time is some time after construction when the induced
negative porewater pressures have dissipated to produce a clay of lower strength. It
follows that for this situation the most critical time is a long time after construction.
Therefore, it is the effective (or drained) strength properties (c', ') which should be
used to check ultimate stability rather than the undrained properties.
Unlike granular soils, cohesive soils are relatively easy to sample and test in the
laboratory. Thin-wall tubes (or u-tubes) are used for taking undisturbed samples of
clays or soils with some clay content. Tubes vary from about 50mm to 80mm diameter,
and wall thickness is about 1mm. The tubes are sealed at either end to retain moisture
content during storage.
For cohesive soils because of the many methods available both in the laboratory and
in the field for the determination of properties, it is comparatively easy to obtain
reasonably accurate undrained strength properties for bearing capacity estimations.
Unfortunately, columns do not always sit nicely in the middle of footings (perhaps
because of architectural constraints), and columns or walls may also impose lateral
loads or moments to the foundations from environmental loading such as earth
pressure or wind loading. Here are some guidelines for accounting for these problems.
It is relatively common to find foundations loaded such that the line of action of a
vertical load does not coincide with the centroid of the foundation area. This is
particularly common at edges of buildings near the site boundary or when axial load
and moment work together. A method for analyzing this situation is given in the
following.
When loading is eccentric, a rigid foundation will no longer give rise to a uniform
pressure distribution on the supporting soil. The variation in contact pressures are
given in the figure to the right.
For 𝑒 c
e
P
𝑃 6𝑒 𝑃 6𝑒
𝑞 1 𝑎𝑛𝑑 𝑞 1
𝐵𝐿 𝐵 𝐵𝐿 𝐵
For 𝑒
q min
4𝑃
𝑞 𝑎𝑛𝑑 𝑞 𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑒 q max e<
B
3𝐿 𝐵 2𝑒 6
The method for handling this situation in design is referred to as Meyerhof's (1953,
1963) effective area method. It
e1
consists of adjusting the area of P
the foundation so that the vertical
load becomes a central load. This
is illustrated for one-way L
eccentricity in the diagram.
P
For an eccentricity of e1 in the L B
direction the new length becomes
𝐿 𝐿 2𝑒 2e 1 L/2 L/2
𝑃 𝑃
𝑞
𝐴 𝐵𝐿′ 𝐵 𝐿 2𝑒
For two-way eccentricity (with e2 in the B direction, not shown here) the new area
becomes
𝐴 𝐵𝐿 𝐿 2𝑒 𝐵 2𝑒
Therefore, with two-way eccentricity, the uniform contact stress (q) becomes
𝑃 𝑃
𝑞
𝐴′ 𝐿 2𝑒 𝐵 2𝑒
This method ensures that the design load employed is always developing stress less
than qmax so that the safety is assured.
The above effective footing dimensions can be used to determine the ultimate net
bearing capacity (q’u) of an eccentrically loaded footing based on Meyerhof (1963)
method discussed earlier. The important issue to be considered here is the effective
L’ and B’ when calculating the shape factors (s). However, when calculating the depth
factors (d) the original width (B) is to be used. Finally, the ultimate load can be
calculated as Pu=q’uA’.
The theories put forward so far deal with generally simple and somewhat idealised
problems. In reality soils are not as simple as this and soil properties often change
dramatically with depth. Therefore, depending on the particular problem, modifications
have to be made for its solution. Some very general applications will be studied here
although often a literature search may have to be undertaken should soil conditions
merit it. There are many specific solutions to be found in the literature.
The ultimate bearing capacity for a strip footing placed in the above layered system
(qu) can be approximated as (Meyerhof & Hanna, 1978):
2 𝐶 𝑃 𝑠𝑖𝑛𝛿
𝑞 𝑞 𝛾 𝐻
𝐵
Where, Pp=passive resistance on the punching surface, =inclination of the passive
force with horizontal, Ca=c’aH=adhesive force, c’a=interface adhesion which depends
on c1.
2𝑐′ 𝐻 2𝐷 𝐾 𝑡𝑎𝑛𝜙′
𝑞 𝑞 𝛾 𝐻 1 𝛾 𝐻 𝑞
𝐵 𝐻 𝐵
Where, Ks=punchinh shear coefficient=(q2/q1, ’1), q1 & q2=bearing capacities of a
surface strip footing placed on thick beds of stronger and weaker soils separately,
qt=the ultimate bearing capacity of a strip footing placed in top stronger
soil=cNc+p’oNq+0.5BN.
The variations of Ks and ca with q2/q1 ratios can be found in Das & Sivakugan (2019)
or Meyerhof & Hanna (1978).
The above equation for strip footing can be extended to a rectangular footing by
incorporating the effect of L and B on the punching volume as:
𝐵 2𝑐′ 𝐻 𝐵 2𝐷 𝐾 𝑡𝑎𝑛𝜙′
𝑞 𝑞 1 𝛾 𝐻 1 1 𝛾 𝐻 𝑞
𝐿 𝐵 𝐿 𝐻 𝐵
For a footing is placed on a thin weaker soil underlain by a stronger soil (i.e., q2/q1>>1)
then the top weaker soil can fail by squeezing out. The ultimate bearing capacity in
this case can be estimated using a semi-imperical relationship as proposed by
Meyerhof (1974) below.
𝐻
𝑞 𝑞 𝑞 𝑞 1 𝑞
𝐻
Where, H=thickness of the top weaker layer and Hf=depth of failure surface beneath
footing in a thick bed of top weaker layer (For vertical loading and strip footing, Hf≈B
for loose sand and clay, Hf≈2B for dense sand; for circular footing, Hf≈B for all soils)
(Meyerhof, 1955). The depth of failure surface decreases significantly with increased
load inlcinations.
qt and qb are to be estimated using a thick layer of each soil separately and using the
Meyherhof bearing capacity equation with shape factors ONLY. The founding depth
(Df) remains the same.
For a strip footing, if a rigid rough base is located at a depth below the footing, which
is within the depth of the failure zone then the full development of the failure surface
will be restricted. Mandel & Salencon (1972) developed an ultimate capacity equation
for strip footing in such case as:
𝑞 𝑐𝑁 ∗ 𝑝′ 𝑁 ∗ 0.5𝛾𝐵𝑁 ∗
Where, N*=modified bearing capacity factors which are dependent on friction angles
and H/B ratios (refer to Das & Sivakugan, 2019). For H>B, the effect of the presence
of a rigid base for lower friction angle (<30o) could be ignored and the general bearing
capacity factors (N) could be used.
For the rough rectangular and circular footing with a rough rigid base, the above strip
footing equation was modified to incorporate the shape factors as:
The previous sections deal with the ultimate bearing capacity of a shallow footing
placed on thick or layered soil deposit. However, in real design, one has to calculate
the size of the footing based on the design load and soil properties so that adequate
factor of safety (FoS) against the shear failure is ensured. In practice, a combined FoS
between 2 to 3 is used on the soil ultimate capacity to determine its allowable capacity,
where a lower value of FoS is associated with higher confidence in soil properties. The
allowable footing load carrying capacity is then compared with the actual load applied
on the footing (ie., Fallowable ≥ Fapplied) to ensure safety against soil bearing failure.
However, this does not address the issue with tolerable settlement which has to be
estimated to ensure that the settlement does not exceed the limiting value during the
serviceability period. This will be discussed later.
Recent design codes use the ultimate limit state (ULS) design approach where the
uncertainities associated with estimating structural loads and soil properties can be
handled separately based on risk classifications. This allows more control on the
design side, for example, if more investigation is done to obtain soil properties then
geotechnical risk can be reduced or FoS can be lowered.
The ULS design equation as you have been practicing in many civil engineering units
can be written as:
𝐸 𝜙 𝑅 .
Where, Ed = Design action effects (actual load multiplied by appropriate load factors);
Rd,ug = design ultimate geotechnical strength;g = geotechnical reduction factor which
depends on how much a designer is confident with the design soil properties. Typically
g can vary between 0.35 and 0.65 (low value for high risk design or low confidence in
soil properties or less investigation has been carried out). Field and lab investigation
may improve the confidence level with soil properties thus less reduction of the
geotechnical capacity. You should follow AS5100.3-2004 Tables 10.3.3(A) & (B) for
further details.
The basic ultimate limit states used in checking strength (Ed) are given below while
many other load combinations can be found in the standards (AS1170, Section 4.2.2).
The long-term and combinations factors can be found in Table 4.1 (AS1170).
7. References
Bowles, J.E. (1996). Foundation Analysis and Design. 5th edition. McGraw-Hill.
Das, B.M. & Sivakugan, N. (2019). Principles of Foundation Engineering. 9th edition.
Cengage, Boston, USA.
TABLE OF CONTENTS
1. INTRODUCTION 3
2. STRESSES IN SOIL MASS 3
2.1 Stress from a vertical point load 4
2.2 Stress below a uniformly loaded circular area 4
2.3 Stress below a line load 6
2.4 Stress below a uniformly loaded flexible strip area 6
2.5 Stress below a horizontal strip load 7
2.6 Stress below a uniformly loaded rectangular area 7
2.7 Stress below a uniformly loaded irregular area 9
3. SETTLEMENT MECHANISMS AND COMPUTATIONS 10
3.1 Immediate settlement 10
3.1.1 Essential features of immediate settlement 10
3.1.2 Appropriate parameters for immediate settlement estimation 10
3.1.3 Immediate settlement based on the theory of elasticity 11
3.1.4 Immediate settlement due to a single concentrated load 11
3.1.5 Immediate settlement due to a uniformly loaded circular area 12
3.1.6 Immediate settlement due to a uniformly loaded rectangular area over an infinite layer 13
3.1.7 Immediate settlement due to a uniformly loaded rectangular area over a finite layer 14
3.1.8 Immediate settlement of shallow foundation in layered soil 15
3.1.9 Typical modulus values for preliminary settlement calculations 16
3.1.10 Typical values of Poisson’s ratio for estimating immediate settlement 17
3.2 Consolidation settlement 17
3.2.1 What is primary consolidation? 17
3.2.2 Calculating consolidation settlement - Oedometer method 18
3.2.3 Typical consolidation parameters 19
3.2.4 Skempton and Bjerrum's correction for overconsolidated soils 20
3.2.5 Other methods for computing consolidation settlement 23
3.3 Secondary consolidation 24
3.3.1 What is secondary consolidation? 25
3.3.2 Calculation of secondary consolidation 25
1. Introduction
Apart from bearing pressure, the other major design consideration for shallow
footings is settlement. Excessive settlement (primarily differential settlement) can
cause a number of problems. Obviously structural distress can be caused
sometimes catastrophically with additional loads being transferred onto other loaded
foundations possibly enough to cause bearing failure.
On applying an increased load to the soil through a footing, the soil is subjected to
increased stress levels, which cause the particles to rearrange in a denser structure
that can accommodate these higher stresses. The volume of voids between the soil
particles decreases, and this decrease is reflected in settlement at the soil surface.
Of course, this rearrangement is only caused by changes in the stresses that the
soil skeleton feels - the effective stresses, rather than the change in total stress,
which includes pore water pressure. So, to compute the settlement of a soil, the
effective stresses in the soil before and after application of the new load must be
determined.
For sands, the rearrangement of soil particles is virtually instantaneous, whether the
sand is above or below the water table. If it is below, the sand structure is so open
that water is not impeded from being squeezed out as the sand grains move to their
denser structure. Sands also respond elastically to applied stress (as long as it is
not excessive). If load is taken off, the sand structure will simply revert to its previous
looser state.
For clays, the settlement is not instantaneous, because water travels slowly through
a clay structure - voids are not continuous, and there are electrochemical attractions
etc. There is a small component of elastic rearrangement, however, the bulk of the
volume change occurs over weeks, months or years as the water is displaced from
the more highly stressed areas. This rearrangement of particles is commonly called
consolidation.
Incidentally, with both sands and clays, this rearrangement of particles in a denser
structure is accompanied by an increase of soil strength.
The stresses due to a vertical point load are in themselves quite irrelevant - point
loads are not applied to soils.
Rather, loads are applied through
footings which distribute the load Q
to an allowable level. However,
this solution can be integrated over
finite areas to determine stresses
due to real footings. Boussinesq R Z
(1885) derived the expressions for
normal and shear stress at any
point considered inside a perfectly z
elastic semi-infinite, weightless r
medium. Since when estimating
the settlement it is primarily the
change of vertical stress which is r
of importance, the equation for z
is presented here.
3𝑄 1
∆𝜎 /
2𝜋𝑧 𝑟
1
𝑧
The above equation cab be simplified to
∆𝜎 𝐼 Where, 𝐼 /
Iis known as the stress influence coefficient, the values with (r/z) can be found in
standard text books.
In order to find the stress directly below the centre of a uniformly loaded circular
area, Boussinesq's equation for z must be integrated over the area.
Consider the small surface element below loaded with q/unit area. The load on this
area can be considered as a point load and is given by q r d dr.
The stress increment (d) at the centre of the circular area at a depth z can be
determined using Boussinesq equation as:
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.5
Topic 2C: Settlement of Shallow Footings
3 𝑞𝑟𝑑𝜃𝑑𝑟 1 q a
𝑑𝜎 / d
2𝜋𝑧 𝑟
1
𝑧
r
dr
For the entire loaded area, the increase of
vertical stress (z) at the centre at depth z z
can be found by integrating the above
equation.
z
⎡ ⎤
1
∆𝜎 𝑞 ⎢1 /
⎥ 𝑞𝐼
⎢ 𝑎 ⎥
⎣ 1 ⎦
𝑧
Where I is the stress influence factor for circular loaded area and the values can
be found in standard textbooks (Das & Sivakugan, 2019) in terms of (z/a, r/a).
10
Influe nce fa ctors for the stre ss be low the ce ntre
of a unifor m ly loa de d fle xible circula r a re a .
2𝑞𝑧
∆𝜎 z
𝜋 𝑥 𝑧
A (x,z)
z
The equation for line load can be further extended to determine the change of
vertical stress at a point below the flexible strip load where the load has a finite width
(B) and infinite length (ie., the case of a wall footing).
𝐵
𝑞 𝑧 𝑧 𝐵𝑧 𝑥 𝑧
4
∆𝜎 𝑡𝑎𝑛 𝑡𝑎𝑛 𝑞𝐼
𝜋 𝐵 𝐵 𝐵
𝑥 𝑥 𝑥 𝑧 𝐵 𝑧
2 2 4
For point A (0≤x≥B/2), the magnitude of 𝑡𝑎𝑛 becomes negative and a term
is to be added to get the solution.
The values of stress influence factor (I) for a uniformly loaded flexible strip footing
in terms of (2z/B, 2x/B) can be found in Das & Sivakugan (2019).
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.7
Topic 2C: Settlement of Shallow Footings
z
2𝐵𝑞𝑥𝑧 A (x,z)
∆𝜎 z
𝐵
𝜋 𝑥 𝑧 𝐵 𝑧
2
∆𝜎 𝑞𝐼 z
z
Where,
1 2𝑚𝑛 𝑚 𝑛 1 𝑚 𝑛 2 2𝑚𝑛√𝑚 𝑛 1
𝐼 . 𝑡𝑎𝑛
4𝜋 𝑚 𝑛 𝑚 𝑛 1 𝑚 𝑛 1 𝑚 𝑛 1 𝑚 𝑛
𝑩 𝑳
Where, 𝒎 𝒂𝒏𝒅 𝒏
𝒛 𝒛
In the above influence factor equation, the arctangent must be a positive angle in
radians. If it becomes negative (m2+n2+1<m2n2) then a term should be added to
that angle.
These stress influence values at the corner of a rectangular loaded area are
presented graphically below. The chart presents I as a function of the parameters
m and n. Please note m and n are interchangeable in this chart.
0.26
I m = 3 .0 m = oc
A rea covered m = 2 .5
0.24 m = 2.0
w ith uniform m = 1.8
norm al load, m = 1.6
0.22 q m = 1.4
m = 1 .2
nz
x m = 1 .0
0.20 y
mz
z m = 0.9
z = qI m = 0 .8
0.18 q.I m = 0.7
z
0.16 m = 0.6
0.12 m = 0 .4
0.10
m = 0 .3
0.08
m = 0 .2
0.06
0.04 m = 0 .1
0.02
m = 0.0
0
0.01 2 3 45 0.1 2 3 4 5 1.0 2 3 45 10
V E R T IC AL S T R E S S BE LO W A C O R N E R
O F A U N IF O R M LY LO AD E D F LE X IBLE
R E C TAN G U LAR AR E A.
When dealing with a foundation of irregular shape, the stress equations become too
difficult for use for conventional design purposes. Newmark however developed a
method to evaluate the change of stress. This method effectively uses a graphical
integration.
NEWMARK CHART FOR VERTICAL STRESS
The method
involves a Newmark Influence Value 0.005
chart which is
placed over a scale
drawing of the
loaded area. The
scale of the area is
such that the
reference line on the
chart ab represents
the depth z to the
same scale.
∆𝜎 𝑞∗𝑁∗𝐼
Where, q= applied stress; N=the number of squares or parts thereof and I=the
influence value given on the chart.
This section describes the essential features of immediate settlement, and the
methods of computation.
For relatively impervious cohesive soils the immediate settlement takes place at
constant volume if the soil is saturated (because the vertical settlement is
accompanied by a compensatory horizontal expansion for = 0.5). As there is no
volume change (V=0), this settlement is elastic and fully recoverable. For cohesive
soils, the volume change takes place over an extended period of time and is known
as consolidation settlement.
For granular soils, which are generally fast draining, the same processes take place
i.e. undrained elastic settlement with no volume change and consolidation
settlement. However, with granular soils, the consolidation settlement occurs so
rapidly that it is not possible to distinguish between the elastic and consolidation
components. Therefore, for granular soils, all settlement is lumped into the one
category of immediate settlement.
Immediate settlement for all soil types are estimated on the basis of elasticity
provided the elastic properties of the soil E and can be evaluated from the site
investigation. Usually, this is possible with cohesive soils, which can be sampled
and tested in a laboratory. For granular soils, however, sampling, in particular
undisturbed sampling is much more difficult, if not impossible in most cases. In-situ
testing techniques do not usually provide elastic parameters but require the use of
empirical methods for settlement estimation as will be discussed later. Some other
methods, such as the pressuremeter do yield elastic parameters. Therefore, the
elastic methods discussed in this section can only be used if we have prior
knowledge of the elastic modulus and to a lesser extent the Poisson's ratio. This
latter value can usually be assumed with sufficient accuracy.
As discussed above all settlement with granular soils is immediate by virtue of the
fast drainage of such soils. The parameters used should reflect this characteristic.
Therefore, the drained modulus (E') and the drained Poisson's ratio (') are
appropriate for granular soils.
For cohesive soils, which are relatively slow draining, the immediate settlement
should be estimated using the undrained modulus (Eu), and the undrained Poisson's
ratio (u). This latter value is 0.5 if the soil is fully saturated.
From the theory of elasticity, the vertical strain at a point within a loaded medium is
given by
Δ𝜎 𝜐
Δ𝜀 Δ𝜎 Δ𝜎
𝐸 𝐸
Where, z = change of vertical stress; x and y = changes of stresses in the
two horizontal directions.
𝜌 𝜀 𝑑𝑧
∴𝑧 𝑟𝑐𝑜𝑡𝜓
𝑟𝑑𝜓 𝑄1 𝜐
∴𝜌 𝜀
𝑠𝑖𝑛 𝜓 𝜋𝑟 𝐸
3.1.5 Immediate settlement due to a uniformly loaded circular area
𝑞 𝑟 𝑑𝜃 𝑑𝑟 1 𝜐
∴ 𝑑𝜌 z
𝜋𝑟 𝐸
Settlement due to whole of loaded area can be obtained by
2 a q
(1 2 ) d dr
0 0 E
2q (1 2 )a
E
It can be similarly shown that the settlement under the edge of a circular load area
is given by
4q (1 2 )a
E
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.13
Topic 2C: Settlement of Shallow Footings
Special cases of the above occur in perfectly elastic soils with = 0.5.
15
. qa 3qa
Centre settlement and Edge settlement
E E
Also, for any point at a depth zo below the centre of the circular area, the settlement
from beneath that point is given by
2𝑞𝑎 1 𝜐
𝜌 /
𝐸 𝑧
1
𝑎
Using the above expression, the settlement of a layer of finite thickness may be
calculated.
The settlement under the corner of a flexible rectangle under uniform q was
proposed by Schleicher (1926).
𝑞𝐵
𝜌 1 𝜐 𝐼
𝐸
The principle of superposition should be used for locations other than at the corners.
𝑞𝐵
𝜌 1 𝜐 𝐼
𝐸
Where, Ι 𝐹 𝐹
Values of F1 and F2 were calculated in graphical form as functions of L/B and H/B
as given below in addition to mathematical expressions.
1 1 𝑚 1 𝑚 𝑛 𝑚 𝑚 1 1 𝑛
𝐹 𝑚 𝑙𝑛 𝑙𝑛
𝜋 𝑚 1 𝑚 𝑛 1 𝑚 𝑚 𝑛 1
𝑛 𝑚
𝐹 𝑡𝑎𝑛
2𝜋 𝑛 𝑚 𝑛 1
Values of F1 ( ) and F2 ( )
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
L/B = 1
2
Depth factor d = H/B
L/B = 2 F1
4
L/B = 5 F2
6
L/B = 5
L/B = 10
8
L/B = oo L/B = 2 L/B = 10
L/B = 1 L/B = oo
10
The elastic settlement below the corner of a shallow footing on a layered soil can be
estimated by considering the theory of superposition. The layer immediately below
the base of the footing is assumed to have a rigid base and the settlement at the
corner of a uniformly loaded rectangular area with a finite thickness (H1) and
modulus (E1) can be calculated. For the second layer, it is assumed that layer 2
exists from the base to a depth (H1+H2) with modulus E2 and having a rigid base.
This assumption allows us to calculate the elastic settlement at the corner of a
rectangular loaded area. However, in this calculation, one has to assume the (H1,
E2) soil as an extra layer which in reality does not exist. Thus to estimate the net
effect of (H2, E2) we have to subtract the settlement effect of (H1, E2) from the
(H1+H2, E2). This concept can be extended to subsequent layers if needed. This
method of superposition is shown in figure below for a two layer case.
The modulus of soil widely varies with tests, strain level, age, compositions, soil
disturbance, soil consistency, stress history etc. It is hard to estimate field modulus
from routine laboratory tests. There is no single modulus value for soils rather a
range is more appropriate. As such in settlement calculations, if the modulus is
approximated from correlations then a range of settlement values is to be estimated.
For approximate calculation, modulus values of a select group of soils are tabulated
below and the values need to be cautiously used as the range is often too large
(Bowles, 1996). It is always desirable that one should use a number of possible
values estimated from a range of tests, such as the values obtained from the
laboratory, field and practical experiences.
The following empirical correlations can be used to estimate modulus of soil (Es)
from undrained shear strength (su) of clay.
Typical values of Poisson’s ratio for common soil types are tabulated on the right
(Bowles, 1996). It is to be noted that for saturated undrained clay, the value of
Poisson’s ratio (u)=0.5 i.e., no change of volume.
Consolidation settlement occurs in all soils, but can only be isolated in fine grained
cohesive soils where the necessary volume changes occur over an extended period.
It takes place after immediate settlement and may take weeks or even years to
complete. The volume changes result in plastic deformations to occur, and this type
In non-cohesive (granular) soils, this form of settlement does not take place as a
time dependent mechanism because of its generally free draining nature. However,
it does take place immediately and therefore it is combined with the immediate
settlement.
Two alternative parameters are often derived from the results of an Oedometer test.
Δ𝑒
𝑚
Δ𝜎′ 1 𝑒
Where, e=change of void ratio; ’=increase of effective stress; eav=average void
ratio during consolidation.
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.19
Topic 2C: Settlement of Shallow Footings
𝜌 𝑚 Δ𝜎′Δ𝐻
𝜌 𝑚 Δ𝜎 Δ𝐻
The second is a pair of indices known as the compression index (Cc) and
recompression index (Cr):
𝐶𝐻 𝜎′ 𝐶𝐻 𝜎′ Δ𝜎′
𝜌 𝑙𝑜𝑔 𝑙𝑜𝑔
1 𝑒 𝜎′ 1 𝑒 𝜎′
The recompression index (Cr) is used (1st part of the eqn. only) for final stress levels
≤ preconsolidation pressure ('c). The compression index is used (2nd part of the
eqn. only) for initial effective stress levels (’o) ≥ preconsolidation pressure ('c).
Where the initial stress is less than the preconsolidation pressure, and the final
stress exceeds it, both indices must be used in conjunction to determine the total
consolidation settlement.
Bowles (1996) summarizes many correlations for Cc for different soil types (Table
2-5) based on parameters such as the soil liquid limit (wL). Actual laboratory data
are always preferred, however, in the absence of such data, typical values may be
inferred from these correlations.
Bowles (1996) suggests that the recompression index (Cr) is typically 10~20% of
Cc i.e.
𝐶 ≅ 0.05 𝑡𝑜 0.1 𝐶
Thus the settlement in an overconsolidated soil will be significantly less than those
generated in a normally consolidated soil.
It is often found that with normally consolidated clays the conventional Oedometer
based methods of consolidation settlement estimation provided a reasonable result.
However, with overconsolidated clays, the conventional method seemed to
overestimate consolidation settlement, particularly when the consolidating layer was
thick relative to the width of the applied load.
In the conventional method the amount of consolidation settlement for any given
layer is provided by the relationship
𝜌 𝑚 Δ𝜎 Δ𝐻
Where, H is the thickness of the sub-layer and is the change in vertical effective
stress.
𝜌 𝑚 Δ𝑢Δ𝐻
𝜌 𝑚 Δ𝜎Δ𝐻
Where, the change in total applied stress () is the same as the change of effective
stress (). The two equations can only become the same when = u. This holds
only when the Oedometer test and the consolidating layer in-situ have the same
boundary conditions.
Since the soil sample in an Oedometer is confined in a ring of steel, the sample
cannot experience lateral yield - all void ratio change must result in a vertical
movement. It follows that for the test to be a model of full scale behaviour, the in-
situ soil must similarly not experience lateral yield. This is approximated only when
the loaded area is very wide compared with the thickness of the consolidating layer.
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.21
Topic 2C: Settlement of Shallow Footings
When this is not the case, lateral yield occurs in-situ. Some of the void ratio change
is accommodated in these lateral movements, and the Oedometer test will
overestimate surface settlements.
The change in porewater pressure due to the change in principal total stresses can
be described by the following:
Δ𝑢 𝐵 Δ𝜎 𝐴 Δ𝜎 Δ𝜎
When no lateral yield occurs in-situ: Δ𝜎 Δ𝜎
and for a saturated soil with B = 1: Δ𝑢 Δ𝜎
Therefore, the Oedometer test is a reasonable model of in-situ conditions.
However, when lateral yield occurs: Δ𝜎 Δ𝜎
Therefore, for overconsolidated clays, where A can be significantly less than 1 (see
table below), from the above pore pressure equation
Δ𝑢 Δ𝜎
𝜌 𝑚 Δ𝑢Δ𝐻 𝜌 𝑚 Δ𝜎 Δ𝐻
It would appear that the difference between the prediction of consolidation
settlement on the basis of the Oedometer test and the actual consolidation
settlement for overconsolidated clays is a function of the geometry of the problem
and the relevant value of A.
oed mv H 1
but c mv H u
z
or c mv u dz
0
3
since u 1A (1 A)
1
z 3
then c mv 1 A 1 Adz
0 1
Relating c to the conventional value of oed
c oed
z 3
0 mv . 1A (1 A)dz
1
where z
o mv 1dz
In order to simplify the above it has been assumed that mv and A are constant.
Then A (1 A)
z
.dz
for 0z 3
0 [Link]
Therefore, only depends on the geometry of the problem. Skempton and Bjerrum
calculated values for for various values of z/B for circular and square footings.
It follows from the above that settlement calculated from the traditional method for
an overconsolidated clay is overestimated although for a normally consolidated clay
(A = 1), the values agree fairly well.
Generally, although shape affects the factors, typical values of are as follows.
For further discussion, see Skempton and Bjerrum (1957) and Parry (1968).
The two methods mentioned below both have in common, very detailed laboratory
testing procedures and a fairly sophisticated analysis. It is important to realize that
the benefits offered by more sophisticated approaches can only be realized if soil
sampling and testing are of the highest standard. Any sample disturbance will
negate the benefits of these methods.
This method involves the stressing of a representative laboratory sample along the
same stress path (both vertical and lateral) as the soil would experience in the field.
First, the in-situ stress state must be determined (not trivial) and replicated. The
changes in the stress-state that result from the new loading must then be imposed.
The strains which occur are measured and the settlement calculated.
This method represents the first truly three dimensional method of analysis. The
necessary assumption is that the solid behaves as an ideal two phase material
consisting of a linearly elastic solid skeleton and an incompressible fluid filling the
voids. In this context, the fluid in the voids acts as a damping agent and does not
affect the ultimate behaviour of the elastic skeleton.
The 3-D field stresses are applied to a representative soil sample in a controlled
manner for both drained and undrained conditions. From elastic theory E’ and Eu
are determined from the laboratory test results. Elastic analyses are performed
using both the drained and undrained elastic parameters. The consolidation
settlement is the difference between those two values of settlement which have
been calculated from elasticity (Davis & Poulos, 1963, Davis & Poulos, 1968).
This method of course raises a most useful tool for the estimation of immediate and
total settlement for clays. If the porewater is regarded as simply a damping agent
and we can conduct laboratory (or indeed field) tests to give the undrained and
drained elastic parameters then immediate settlement for say a rectangular
foundation placed on an infinite soil layer is given by
1 𝜈
𝜌 𝑞𝐵 𝐼
𝐸
1 𝜈′
𝜌 𝑞𝐵 𝐼
𝐸′
It would appear then that the consolidation settlement would be given by
𝜌 𝜌 𝜌
Note: This method could only be regarded as approximate because the
consolidation component is assumed to be elastic which is not strictly true.
Therefore, this method is likely to severely underestimate c.
Some researchers contend that the theory that secondary consolidation occurs at
zero excess pore pressure or constant effective stress is incorrect. Present evidence
seems to indicate that the primary/secondary division is purely artificial and arises
largely as a result of traditional laboratory testing techniques. The relationship
between "primary" and "secondary" appears to be very much influenced by the
thickness of the sample tested and attempts have recently been made to present a
unified theory which eliminates the distinction entirely. Their argument is essentially
that since water flows from the soil during secondary consolidation there must be a
small hydraulic gradient, and therefore small excess pore water pressures. This then
eliminates the distinction between secondary and primary consolidation.
Δ𝑒 Δ𝑒
𝐶 𝑡
𝑙𝑜𝑔𝑡 𝑙𝑜𝑔𝑡 log
𝑡
𝐶 𝐻 𝑡
𝜌 log
1 𝑒 𝑡
of primary consolidation, the time t1 must be chosen with due regard for the
differences between laboratory and field rates of consolidation.
Mesri and Godlewski (1977) have compiled the ratio of secondary compression (C)
to compression index (Cc) for a number of soils as:
𝜌 𝜌 𝜌 𝜌
The above components have been computed with the assumptions that the loaded
area is flexible, and imposes a uniform pressure, and that the loaded area is at the
surface of the soil. These conditions are often not met, and corrections must be
applied to account for footing rigidity, and footings being buried.
Settlement calculations are generally based on the assumption that a uniform load
exists over the foundation for all stages of settlement. This may be nearly true for
earth embankments, dams and tanks. However, it is not true for non-flexible bases
(i.e. those which are rigid or semi-rigid), such as rigid concrete footings or in the
semi-rigid class, compacted road bases.
Q Q
q= q=
C O N TAC T R2 2R (R 2 -a 2 )
STRE SS
SET TL EM EN T
PRO FILE
(a) (b)
For a flexible base in clay, the settlement profile would be non-uniform and pressure
would be uniform as shown in (a) above, whereas a rigid base would apply non-
uniform pressure and experience uniform settlement like (b) above.
The settlement profiles and contact pressures for flexible and rigid footings on sand
are quite different from those shown for clay in the above diagrams. This is because
the modulus of sand generally increases with depth, whereas that for clay tends to
be uniform, and because the lack of lateral confinement at the edge of the footing
tends to push the surface sand outward, reducing edge stresses for both types of
footing. Therefore, in the absence of any other information, the assumed rigidity
correction can be applied to footings on sand as well.
The elastic solutions derived previously are invalid if foundation loads are applied at
some depth below the ground surface. Fox (1948) established the depth correction
factor (If) which showed that the calculated elastic settlement should be reduced by
a factor of between 0.5 and 1.0. The depth correction factors (If) for different values
of D/B, L/B and Poisson’s ratios () can be found in Das & Sivakugan (2019) (refer
to Figure 9.5 in 9ed.).
Just as classic one-dimensional consolidation theory was the basis for estimating
consolidation settlements, this theory also provides the basis for estimating the rate
of settlement. In this theory, which you have learned in CIV3247: Geomechanics-II
unit, the essential time relationship is given by the expression below.
𝐶 𝑡
𝑇
𝐻
Where, T = time constant for the particular initial porewater pressure distribution
(dimensionless); Cv = coefficient of consolidation (m2/year); t = real time (yr); and
H = maximum length of the drainage path (m).
In order to calculate the times of consolidation, it is necessary first to find the stress
distribution produced in the soil due to surface loading. With relatively thin
consolidating layers and wide foundations, this distribution can be approximated as
rectangular, whereas with thick soil formations and/or small loaded areas, a
triangular approximation may be more appropriate. Each case must be judged
individually.
Clay
Clay Approx.
triangular
Sand Approx.
rectangular
Sand
Standard time factor tables are available, which list values for T depending on the
drainage conditions and the initial excess pore water pressure distributions. This
diagram shows three such initial pore water pressure distributions. The following
table lists time constants T for values of U from 10 to 90% of consolidation
completed for the three cases. For difficult stress distributions, T may be found by
analytical methods if necessary i.e., finite difference techniques may be used.
BOTTOM
PERMEABLE
The value of cv is determined from the Oedometer test for the appropriate stress
range of the problem.
As shown in the figure below, for a layer of thickness z, which is double drained (a),
H = z/2, but where only a single drainage layer exists (b), H = z.
Sand Sand
Soft Soft
z z
Clay Clay
Knowing T, Cv, H and t1, the percentage consolidation (U), can be computed, and a
time-settlement curve can be derived for the particular consolidation settlement. At
any time t, the settlement of a foundation is given by
𝜌 𝜌 𝑈𝜌
5.2 Problems in Practice
The previous section presents a picture which seems to indicate that the estimation
of rates of settlement is fairly simple. Unfortunately, in practice, the actual situation
is much more complicated.
In situations where the loaded area is relatively wide in comparison with the
consolidating layer the assumption made above concerning one-dimensional
consolidation theory is probably reasonable. However when the width of the loaded
area is about the same or smaller than the consolidating layer, then it is likely that
considerable radial drainage will take place. Clearly the magnitude of the horizontal
permeability will have a significant influence on this form of drainage. However, as
soil layers including the more permeable are laid down in horizontal layers, it is likely
that in many situations the horizontal permeability will be greater than the vertical.
It follows that the influence of radial drainage will be significant in most situations,
with the exception of those involving very wide loaded areas on relatively thin
formations of clay.
In order to make better estimates of time rates of settlement when radial as well as
vertical drainage occurs it is necessary to undertake much more complex analyses.
This may be in the form of three-dimensional theory or by use of an approximate
solution involving the superposition of vertical and radial drainage. Either way these
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.31
Topic 2C: Settlement of Shallow Footings
The other major influence on the actual rate of settlement is the presence of fabric
within a soil mass (Rowe, 1972). Fabric includes all the many forms of non-
homogeneities which are common in most soils: joints, fissures, fractures,
laminations, lenses and root fibres. These features tend to be much more permeable
than the intact soil and their presence tends to accelerate rates of settlement.
The problem arises because the sample of soil contained in an Oedometer test is
only very small. Therefore, even though a given layer of soil may contain much fabric
there is usually little of this fabric in the test sample. It follows that the results
obtained from such tests will indicate coefficients of consolidation considerably
lower than for the soil mass. A typical value for cv for an intact clay is around 1
m2/year. Depending on the fabric
10 6
Piezom eter readings close content for a soil, cv for the soil mass
to highly fissured rock
From full scale piezom eter
can be well in excess of 100 m2/
Coefficient of consolidation , cv - m2 / year
5
10 and settlement readings year, demonstrating the dramatic
From sm all scale influence of fabric.
laboratory testing.
10 4
This figure on the left shows cv
c lo s e ly sp
a c e d fi ss u values determined or inferred for a
re s
10 3 stiff fissured clay, and serves to
illustrate the effect of fabric. The
methods used to overcome this
10 2 problem generally involve attempts
we ll sp ac ed fis su
to obtain samples for testing so that
re s a representative component of fabric
10
is included. An appropriate size is
one where the sample tested is
0
1 100 200 300 400
several times the spacing of the
Effective stress - kPa
fabric. For example, if fissures are
spaced at about 100 mm then the
Department of Civil Engineering, Monash University
Revised: 5/2020
CIV4249: Foundation Engineering 2C.32
Topic 2C: Settlement of Shallow Footings
sample required would need to be in excess of about 300 mm (about the size of a
"Rowe" consolidation cell), perhaps somewhat more. Clearly there is a limit to the
size of sample that can be retrieved for laboratory testing. Alternatively, larger scale
field testing may have to be carried out. Once again for the larger projects where
more precise rate estimations are required, it may be economically feasible to
construct large instrumented trial embankments from which a realistic cv can be
back-calculated.
This is the opposite of settlement and can occur after opening an excavation i.e.,
the removal of overburden. The swelling involved can be calculated by the normal
methods, but with a negative stress or porewater pressure change. Swelling effects
will be restricted to excavations in clay.
The swelling is associated with an increased moisture content and void ratio, and
therefore with a loss of soil strength.
Should the ground water table change by any significant amount during the lifetime
of a building changes of effective stress in the soil will result. These stress changes
may bring about additional consolidation or swelling.
For example, should the ground water table fall the effective stress in the soil will
increase as the soil particles have to support the full weight rather than the buoyant
weight of the overlying soil. This will result in additional consolidation. Pumping
water from aquifers below the consolidating soil, and thereby lowering the porewater
pressures has the same effect. Such effects are most vividly seen in Venice,
Bangkok and Mexico City.
Mexico City has also suffered this settlement, as much as 7m total and 2m
differential in parts of the city. Bangkok is apparently settling at a rate of about
10mm/year (Zeevaert, 1953).
7. Overall summary
When the elastic and consolidation properties of a soil can be derived by appropriate
laboratory and/or in-situ tests, the total amount of settlement is given by
𝜌 𝐼 𝐼 𝜌 𝜇𝜌 𝜌
Note: Some elastic solutions may include the rigidity and depth factors in it and one
should not apply again to estimate the final settlement. If the stress calculation is
based on the assumption of a flexible surface footing then the above corrections are
applicable.
by the application of the appropriate values of E and . It follows that for the
settlement of non-cohesive soils, the elastic methods described previously should
be used if the elastic parameters can be established. It is almost impossible to
retrieve samples of non-cohesive soils for testing in the laboratory and so in-situ
methods must be used to estimate the appropriate elastic settlements. The correct
estimation of these parameters is particularly important as settlement is generally
the controlling factor for foundations on sands.
Meyerhof (1956) presented a correlation between the net bearing pressure and SPT
blow count (N60) for 25 mm of maximum settlement of shallow footing in sand.
Bowles (1977) proposed a modified form of the correlation to incorporate the depth
factor (Fd) and elastic settlement as:
.
𝑞 𝑘𝑃𝑎 𝐹 (for B ≥ 1.22 m)
.
Burland & Burbidge (1985) proposed a method of calculating the elastic settlement
of shallow footing in sandy soils based on the SPT N60. Details of this method can
be found in Das & Sivakugan (2019).
𝐼
𝑆 𝐶𝐶 𝑞 𝑞 Δ𝑧
𝐸
The pressuremeter stress-strain results over the depth of influence zone below the
footing (refer to CPT) can be used to calculate the elastic settlement in sand based
on a method proposed by Briaud (2007). This method will not be explored here.
The major problem with the plate bearing test is that the zone of influence of the
300mm plate is substantially less than that of the full-scale footing - at approximately
4B, the zone of influence is 1.2m. By comparison, for a 3m wide footing, the zone
of influence is about 12m. If soil conditions improve (or deteriorate) with depth, or
the soil is stratified, the results of the plate bearing test can be very misleading.
Despite the difficulties associated with the plate bearing test, for soils with large
boulders and cobbles (even fill materials) this method may be the only reasonable
method of predicting settlements.
Terzaghi and Peck (1967) give the following relationship for the interpretation of the
test.
2𝐵
𝑆 𝑆
𝐵 0.3
Where, S= settlement of foundation; S1=settlement of 0.30 m square plate; B= width
of foundation.
It is noted that this method is not too satisfactory for loose and organic soils. The
plate bearing test is described in more detail in Bowles (1997).
The theory will not be considered now but suffice it to say that where such
occurrences may occur special care should be taken in the design of foundations.
Similarly, individual footings that support a structure may settle by different amounts
for a variety of reasons. These differential settlements can cause distress in the
structure if they are not taken into account in the initial structural design. This section
covers the reasons for differential settlement, and the approach used for its
estimation.
Often if the whole structure settles uniformly the structure itself may suffer no
distress. However, differential settlements can be extremely troublesome.
Differential settlement can occur for a whole variety of reasons. However, these can
be conveniently grouped as those due to discernible differences in boundary
conditions, and inherent differential settlement. These are outlined in the following
sections.
Variations in soil strata such as lenses of sand in clay strata or vice versa.
Irregular bedrock surface.
Variations in foundation loading such as towers with low side wings and
machinery foundations.
The following table shows some typical differential settlement values for shallow
foundations in terms of maximum settlement values (O’Brien, 2012).
/ L = 1/150. Safe limit for flexible brick walls with L/H>4. Danger of
structural damage to most buildings. Cracking of panels and brick walls.
/ L = 1/250. Visible tilting of high rigid buildings.
/ L = 1/300. First cracking of panel walls
/ L = 1/500. Safe limit for buildings where cracking is not permissible.
/ L = 1/600. Danger limit for frames with diagonals.
Obviously all these are dependent on the structure. Therefore, each structure should
be designed based on its own merits and code specified values.
Frost Heave
Hillside creep
Mining subsidence
Swallow holes
Shock loading - e.g. earthquakes nearby collisions!
Shrinkage of soils and swelling due to seasonal moisture and
temperature changes
As for the previous, but due to vegetation etc.
13. References
Bowles, J.E. (1997). Foundation analysis and design. 5th ed. McGraw-Hill.
Briaud, J.L. (2007). Spread footing in sand: Load settlement approach. J. Geotech.
Geoenv. Engg., ASCE, 133(8), 905-920.
Davis, E.H. & Poulos, H.G. (1968). The use of elastic theory for settlement prediction
under three dimensional condition. Geotechnique, Vol. 18, 1968.
Lambe, T.W. (1964). Method of estimating settlement. ASCE. SM5, Sept. 1964.
Mesri, G. & Castro, A. (1986). C/Cc concept and Ko during secondary compression.
J. Geotech Engg. ASCE, Vol. 113, No. 3 : 230 - 247.
Rowe, P.W. (1972). The relevance of soil fabric to site investigation practice.
Geotechnique June 1972
Terzaghi, K., Peck, R.B. & Mesri, G. (1996). Soil Mechanics in Engineering Practice.
John Wiley & Sons, New York.
TABLE OF CONTENTS
2. GENERAL APPROACH 5
2.1 Ultimate base resistance 6
10. REFERENCES 54
There are a number of methods, which can be used to establish the axial load
capacity of single piles. These include static formulae used by the designer in
advance, static load testing which is performed after installation, dynamic
formulae which are applied during installation of driven piles and dynamic testing
which may be performed during or after installation.
This chapter deals exclusively with the use of static formulae, which is the basis
on which design usually rests. Occasionally, and particularly for very large
projects, static formulae will be correlated against static or dynamic load testing
as part of the design process, however, in general, such testing is only conducted
during the contract period. Such confirmation is very important, as soil and rock
conditions on site may vary considerably from those established in the (limited)
site investigation. In addition, static formulae are based on approximate models
of soil and pile behaviour that must be verified by actual load testing. Some form
of load testing is carried out on most piling projects. The various load test methods
will be covered elsewhere.
Pile driving formulae are used to confirm that the capacity of driven piles is in
accordance with the design and/or contractual requirements. These formulae
relate pile capacity to the hammer energy. The most common formula used in
Australia is the Hiley Formula. In the U.S., the Engineering News Formula is used.
High factors of safety should be used because of the following (incorrect)
assumptions common to all methods:
In order to apply the dynamic formulae, workmen have to “take the set” on the
pile.
Static load tests are performed either to establish the ultimate capacity of the pile,
or to ensure that the pile has a sufficient factor of safety (1.5 to 2.0 is common)
and a satisfactory performance at working load.
It should be noted that this technique is often used for construction control, as the
energy transfer efficiency of the pile driving hammer, the maximum compressive
and tensile stresses in the pile, and the location and severity of any structural
damage in the pile can be established.
2. General approach
The resistance of a pile has two components: the side resistance and the base
resistance. It is generally assumed in design that these components can be
computed independently and then summed.
The expression relating pile load at failure (Pu) to ultimate soil resistance is
therefore given by
Pu Qb Qs Wp
Where, Qs= the ultimate shaft resistance; Qb= the ultimate base resistance; and
Wp= the weight of the pile.
2.1 Ultimate base resistance
The ultimate base resistance can be estimated from the basic Terzaghi’s bearing
capacity equation as:
B
Qb Ab c b N c P 'ob (N q 1) N Pob
2
Where, Ab= the area of the base; cb= the representative cohesion at the base of
the pile; Nc, Nq and N = bearing capacity factors; Pob and P'ob=the total and
effective vertical overburden stresses at the base of the pile; = the
representative effective unit weight at the pile base; and B= the width of the pile
base.
If the weight of the pile is included, the expression for the ultimate base resistance
becomes
B
Qb Ab c b N c P 'ob (N q 1) N Pob Wp
2
However, the weight of the pile is similar to the weight of the soil displaced or
replaced by the installation of the pile, i.e.
Wp Ab Pob
and further, for a pile length, L, being very much greater than its width, B, it is
reasonable to assume that
B
Wp Ab N Ab Pob
2
Similarly, for > 0, Nq Nq - 1.
of undrained shear shtrength is to be used. If a pile just sits on the clay strata
then a reduced value of Nc=6 is to be used and Nc=9 at 3D (D=diameter of the
pile) penetration in the bearing layer (Bell & Robinson, 2012).
For a soil possessing both cohesion and friction, there will be two components of
shaft resistance - these are treated separately.
The cohesive component of the total shaft resistance (Qsc) will be a function of
many factors, but for simplicity, it is assumed that it is directly proportional to the
cohesion of the soil, i.e.
Qsc = As c s
Qs = Qsc + Qsf = As [cs + KP'os tan] This equation is valid for c-soil.
Where, As=area of the pile shaft; =adhesion factor; cs=average undrained shear
strength of the clay (cs=su).
Adding together the individual components of side and base resistance yields the
general expression for the total resistance of a pile loaded in axial compression.
However, since the unit shaft resistance along the full length of the pile may vary
significantly because of differences in the types and properties of soils
encountered, the cross-sectional area of the pile and the effective overburden
pressures, it would be more correct to express this relationship as
L
Pu Ab [c b N c P 'ob N q ] C [ c s KP 'os tan ] dL
0
Where C is the length of the pile perimeter [e.g., circular pile of diameter
D=C=πD; square pile of side B=C=4B].
Within the limitations of the assumptions discussed above, the above equations
are reasonably correct for the general case of any pile in any soil for either drained
or undrained conditions as may be appropriate. However, major difficulties arise
when the equations are applied. Some of the reasons for these difficulties are as
follows.
The aim of any site investigation and laboratory test program is to retrieve
samples with a minimum of disturbance, so that the in-situ strength can be
estimated as closely as possible.
However, the installation of a pile causes a great deal of disturbance. Driven piles
are forcibly inserted into the soil, resulting in soil being displaced down and out.
The stress levels in the soil increase, and with time, the strength of the clay could
be expected to increase accordingly. On the other hand, the soil is failed by the
insertion process, and a remoulded strength could be expected to apply.
If a pile is constructed by drilling, the same remoulding of soil at the wall of the
hole could be expected. However, as opposed to a driven pile, removal of the soil
by drilling will cause a reduction of stress state, and possible loss of strength.
In addition, the water used to assist the drilling can soften the walls of the hole
being drilled.
The values of cohesion are usually established on relatively small samples of clay
(say 50 to 100 mm of diameter, and length to diameter ratio of 2). By contrast,
the volume of soil involved with failure of a pile under load is significantly greater.
If the soil is isotropic and homogeneous, this difference will not have any
consequence. However, soils in general contain significant fabric. In the case of
fissured clays, the structure is one of intact blocks surrounded by weaker
discontinuities. In such a case, the result of a test on an intact sample will
overestimate the mass strength of the soil (Rowe, 1972).
This figure shows the results of a number of strength tests conducted on fissured
London clay. The strength was determined by 5.5 mm diameter penetration tests,
38 and 98 mm diameter triaxial tests, and 865 mm diameter plate load tests. As
the figure shows, the strength measured by the small tests is very much greater
than the larger triaxial
Shear strength kN/m2
tests. These triaxial
tests in turn also show 0 100 200 300 400 500 600
greater strength than 865mm [Link] tests
obtained with large scale 38mm diam. triaxial Sand and gravel
98mm diam. triaxial
tests (i.e. the size of
Depth below G.L. ( metres )
Penetration tests
piles). The reason for
5 Soft brown clay
this is that with fissure
spacings in the range 5.5mm diam. penetration
Brown and grey
between about 100 and tests
mottled clay
250 mm, the smaller 10 38mm and 98mm
tests are unlikely to be diam. triaxials
Stiff grey clay
influenced by the weaker
fissures and therefore a
high strength is 865mm diam.
15
recorded. For the triaxial plate tests
tests, the fissures will
have some effect and
therefore the recorded
strength is reduced. For the large plate tests, the full influence of the fissures will
come into effect to produce an even lower strength.
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3A.10
Topic 3A: Axial Capacity of Piles
For this reason, for piles founded in fissured clays, it is common to introduce a
reduction factor for the cohesion acting at the pile base. The cohesion cb is often
replaced by cb, where < 1. This reduction factor allows for the possibility that
bearing failure will be partly along fissures and partly in intact material.
By contrast, the shaft cohesion is not similarly reduced, as failure of the shaft is
locally at the shaft soil interface and does not allow for failure through fissures in
the body of the clay.
There is also a rate effect which causes a difference between the cohesion that
is derived from laboratory or in-situ tests and the cohesion that applies to a loaded
pile.
Soil investigation tests are normally conducted fairly quickly, and it is common
(although not universal) to base pile capacity calculations on undrained strength
parameters. By contrast, the dead loads and permanent live loads which often
account for the majority of load on a pile are gradually applied over several
months, if not years. This gradual loading allows the soil to consolidate and gain
strength under the imposed loading.
strengths determined from triaxial tests may bear little relevance to the strength
of clay soils during pile installation.
It should be noted that with time, the elevated pore pressures generated during
pile driving dissipate, and with this dissipation, consolidation of the soil occurs.
This consolidation is accompanied by strength increases in the clay. In practice
this is observed as an increase in resistance to driving after pile installation. Piles
that may have driven 50mm with each hammer blow during installation may only
move 0.1mm or less hours or days after installation. This effect is commonly
called pile setup or pile freeze (US). Pile set-up has been observed to continue
(at a decreasing rate) over months or years, as demonstrated by the field
research results shown below.
The set-up effect is also the main reason why pile load tests are preferably
conducted a minimum of 2 days (but preferably 7 days) after installation.
300 30
150 15
300 x 150mm tapered timber (Drammen)
100 10
300 x 125mm I-Beam 150mm (8 in) steel tube (San Francisco)
50 5
(Gothenberg)
0
1 5 10 50 100 500 1000
Time after driving in days
To determine the end bearing and shaft resistance components of pile capacity
in frictional soils, it is necessary to determine the soil friction angle (), which is a
measure of its strength. Although this value doe not appear in the above bearing
capacity equations, it is implied in the bearing capacity factors, Nc and Nq, and
in the pile-soil friction angle ().
Determination of the friction angle () for granular soils is very difficult, as it is not
possible to retrieve undisturbed samples of such materials. In general, it is
necessary to perform in-situ tests, and apply some empirical formulae to the tests
to infer a friction angle. This process is of itself fraught with danger and
complications, as boundary conditions in these tests cannot be controlled, and
the parameters which are important are often not well understood - see for
example the many and conflicting correction factors which are proposed for SPT
tests, or proposed material correlations for CPT tests.
Now add to this the effects of pile installation, which could be densification in the
case of driven piles, or stress-relief (loosening) in the case of drilled piles. The
friction angle of the sand is intimately related to the density (or relative density) of
the granular material i.e., how closely the particles are packed. It is evident that
the post-installation friction angles and in-situ friction angles will differ, and that
proper account of these differences must be made in the design process.
The bearing capacity factors, Nc and Nq, have been developed from approximate
theories by a number of different researchers. These theories have also been
modified empirically to make theory and practice match more closely. Depending
on which theory is adopted, different values of these parameters will be obtained.
The bearing capacity factors are principally dependent on the value of adopted
(in itself a difficult choice). For low values of , the bearing capacity values are
generally insensitive, however, for high values of , Nc and Nq vary considerably,
and a correct choice is critical.
One of the most critical aspects of pile design is proper consideration of the
influence of installation method on the strength of the soils that support the pile.
Driving processes which displace soil tend to increase soil strength and pile
capacity, whereas drilling processes which tend to remove soil tend to decrease
soil strength and pile capacity.
It would seem from the above that driving processes would in all cases be the
preferred method, however, choice of pile type is dependent on column loads,
soil type and stratigraphy, equipment availability, noise and vibration concerns,
historical precedent, labour considerations, and most of all price. The choice is
by no means straight forward.
It would be normal for most projects to check the design by means of one of the
other methods mentioned previously (e.g. static or dynamic test loading of
selected piles). This is particularly appropriate for bored piles for which there is
no driving record. In general, the more significant the job, the more sophisticated
and expensive the testing program justified. For small jobs, because of the
expense of such testing, it is more usual to adopt conservative estimates of the
variables used in the static formulae; driven piles may be checked for these jobs
using only driving formulae with high factors of safety.
The conventional approach for piles installed in clay is to assume that the clay is
saturated and undrained conditions apply. This implies that the angle of friction
of the clay () is zero and the relevant cohesion is the undrained cohesion (cu)
i.e., a total stress approach.
It must be stressed that although this is the traditional approach to predicting pile
capacity, there would appear to be little theoretical justification. After all, the
majority of load borne by piles is often long-term load, and in any case, it is the
effective stresses in a soil which dictate its strength. The undrained approach has
probably been adopted because undrained strength tests of clay are quicker and
cheaper, and static load tests are in any case relatively quick, and therefore
represent an undrained condition.
𝑃 𝑄 𝑄 𝐴 𝑓 𝐴 𝑓 𝐴 𝑐 𝑁 𝐴 𝛼𝑐
The values of Ab and As are determined by the pile geometry, however, the other
variables require careful examination. fb and fs are unit base and shaft
resistances of pile in kPa.
On the basis of Skempton's work, it is now common practice to use the value of
9 for all deep foundations in clay. As such foundations usually have a D/B ratio
much grater than 4, this value is nearly always appropriate. However, it is
important to note that where the D/B ratio is less than 4 (Bell & Robinson (2012)
reported this limit as 3), Nc should be reduced appropriately.
The point bearing capacity in the lower layer could be assumed to increase
linearly from the full bearing capacity of the softer overlying soil at the surface of
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3A.15
Topic 3A: Axial Capacity of Piles
the underlying stiffer soil to the full bearing capacity of the stiffer soil at the critical
penetration ratio.
It was noted previously that there were many factors which contributed to the
differences between the cohesion measured during site investigations and the
value which is applicable to the base resistance of a pile in clay (see Difficulties
in application). Nevertheless, in the absence of any better value, the value of
cohesion determined from high quality laboratory or in-situ tests is used for driven
pile design.
The effect of fabric, and in particular the presence of fissures, may be ignored for
driven piles, but is important for bored piles. For driven piles, the driving process
causes substantial disturbance a remoulding of soils, which tends to destroy the
original structure of the clay.
However, for bored piles, where soil is removed, the remoulding of soil around
the pile does not occur. In the best-case scenario, the soil fabric remains intact,
and in the worst-case scenario (equivalent to poor construction techniques), the
fabric opens due to stress relief resulting from excavation. As noted previously
(see Scale effects), it is usual to incorporate a reduction factor to account for the
effect of fissures on the end bearing of bored piles, thus:
Qb Ab c b N c
Where = a reduction factor which relates the ratio of the full scale strength to
the smaller sample strength.
The value of shaft cohesion (cs), which is adopted in driven pile design is the
undrained value measured in conventional small scale laboratory or in-situ tests.
Any variation between this value and the value which is applicable to a full scale
driven pile is accounted for by the adhesion factor (). This adhesion factor is
determined empirically from the results of both full-scale and model pile load
tests, and ideally accounts for installation, disturbance, scale and loading rate
effects. There have been many such empirical correlations proposed but is
usually given as a function of undrained cohesion (cu) and reduces with
increasing cohesion (see Poulos and Davis, 1980).
The adhesion factor () is a purely empirical factor which relates the likely pile
shaft resistance to values of undrained cohesion (cu) of undisturbed clay. There
is no theoretical basis for .
As discussed previously, there are many factors which can cause the cohesion
at the interface of a pile shaft to differ from the cohesion measured during site
investigations. However, it is generally assumed that the adhesion factor ()
takes all these factors into account. The problem therefore is to select an
adhesion factor which is appropriate for the particular soil and method of pile
installation. The following discussion will consider driven and bored piles
separately.
It demonstrates that the adhesion factor can vary a great deal for a particular type
of pile installed under fairly similar conditions. A very courageous person has
indicated a design curve on this figure, but adopting this recommendation would
not necessarily guarantee a safe design. If a conservative approach were
adopted, a lower bound curve to the data set would result in grossly expensive
foundations. It would seem that some form of testing of the installed piles is
necessary to confirm capacity.
For soft insensitive clays, it is normal practice to adopt adhesion factors of about
unity, or even higher. This is because the immediate loss of strength (from
undisturbed to remoulded values) due to driving of the pile is compensated, or
even over compensated by the subsequent consolidation of the soil. However,
for sensitive soils, the immediate loss of strength will be much higher, and may
not be fully compensated by the subsequent consolidation. In such cases, an
adhesion factor of less than unity would be appropriate.
Based on the design experience in off-shore conditions, the adhesion factor can
be estimated in terms of overconsolidation ratio of the clay as a function of (cu/’v).
The variations of adhesion factor is shown in figure (a) as shown below (Semple
& Rigden, 1984; cited from Tomlinson & Woodward, 2008). For the most of the
onshore pile design, however, the effect of pile slenderness effect is minimal. For
L/B (i.e., length to diameter ratio) ratio exceeding 50, the earlier adhesion factor
is to be further reduced by multiplying it with a length factor (F) as per the figure
(b) below. The length factor will allow both for the flutter effects and the residual
or part-residual shear strength conditions at the interface.
. . .
𝑓 0.5 𝑐 𝜎′ 𝑂𝑅 0.5 𝑐 𝜎′
[take the value whichever is larger]
There are many influences on the adhesion factor of bored piles in clay that are
inherent in the method of construction. These are as follows:
stress relief that occurs on the walls of the borehole to produce swelling and
weakening;
migration of groundwater to this area;
opening of fissures;
formation of a slurry against the borehole walls due to water seepage;
formation of a cake of bentonite if this is used in the drilling process;
water migration from the in-situ concrete;
shrinkage of the concrete.
locations about the world, the adhesion factors generally range from 0.3 to 0.6
(Tomlinson & Woodward, 2007).
Perhaps the best known of these is the study by Skempton (1959), in which it
was found that the average adhesion for London clay was about 0.45 for normal
length piles, which were drilled and filled with concrete soon after. For short piles,
the adhesion factor was about 0.3 is suggested- the reason given being that
heavily fissured clay reduced the resistance (perhaps in pullout the failure was
through the fissures rather than at the pile-soil interface). Skempton also found
that the shaft resistance did not exceed about 100 kPa. For a small diameter
bored piles and piles with enlarged base where long delay beween drilling and
concreting can take place, an adhesion afctor of 0.3 is to be used. If the enlarged
base is to be constructd in a fissured clay, Tominson & Wooodward (2007)
suggested ignoring 2B length of shaft friction behind the enlarge base.
The effective strength parameters for cohesive materials suggest that the
strength of the clay is stress-dependent i.e., the more stress applied, the denser,
and therefore stronger will be the clay in the long term. The (stress-independent)
cohesion value we used for the total stress approach can now be assumed to be
zero, and the general pile equation can be simplified to:
𝑃 𝑄 𝑄 𝐴 𝑓 𝐴 𝑓 𝐴 𝑃′ 𝑁 𝐴 𝐾𝑃′ 𝑡𝑎𝑛𝛿
It should be noted that the assumption that c' = 0 may be reasonable for soft clays
but for stiff clays, a small but important cohesion is generally measured. The
assumption may therefore be an oversimplification.
With regard to the side resistance component, the unit shaft resitance (fs) can be
determined based on the effective stress parameters as:
This should be compared with the total stress adhesion factor (), which is
defined by 𝛼 .
For most types of clay, ' varies between approximately 20 and 40, although it
must be stated that values much less than 20 are fairly common especially for
the more active clays. Mitchell (1976) has proposed the following relationship for
the critical void ratio (i.e. residual) friction angle for clays:
sin cv
0.8 0.094 lnIP
There may be some justification for using a residual friction angle in pile design
given the disturbance caused by installation. If the clay is sensitive, and
remoulding of the soil has not occurred (e.g. in bored piles with careful
construction control), higher values can perhaps be considered.
A range of 6 < IP < 100, yields values of ' between 39 and 21.5,
0.25 0.30
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3A.21
Topic 3A: Axial Capacity of Piles
The above range is relatively narrow, indicating that is not highly dependent on
clay type.
Parry and Swain (1977) adopted a similar approach and derived a slightly more
complex equation for , namely
0.30 0.35
This range is a little higher than derived from Burland's equation and does seem
to fit the test data a little better.
However, while there does seem to have been some form of success with soft,
normally consolidated clays, the same cannot be said for stiff, overconsolidated
soils. The reason for this is that the K value can be very variable and its
relationship to Ko is greatly influenced by the method of installation. The value of
Ko can also vary a great deal in overconsolidated clays.
By contrast, for bored piles in stiff London clay, it was found that the effective
stress method provides an upper bound to the shaft friction measured in field pile
load tests. The bored piles presumably demonstrated lower shaft friction than
predicted because K < Ko as a result of stress relief during drilling.
These general trends are supported by Meyerhof (1976), who demonstrated that
for driven piles in stiff clay, K 1.5Ko, while for bored piles, K 0.75Ko.
For the base resistance component, there appears to be very little available in
the literature. However, according to Poulos and Davis (1980), it would appear
reasonable to use the effective overburden pressure at the toe of the pile, P'ob
along with a value of Nq that has been estimated from a knowledge of the
effective angle of friction in the same manner as granular soils. However, while
in principle this approach seems reasonable, it should be handled with some care
as it has not really been proven in engineering practice.
This method is not commonly used for land-based piles, and should only be
applied where an assumption of normal consolidation is appropriate.
As c m (m
2c m )As
m 2
cm
It must be borne in mind that the lambda method was proposed for marine pile
applications, for which an assumption of normal consolidation is usually
appropriate. For normally consolidated clays, it is not unreasonable to find that
the ratio of undrained cohesion to effective vertical stress is of the order of about
one third (Skempton, 1957), i.e.
cm 1
v 3
5
The accuracy of this can be established in relation to the typical
recommendations for shown previously. A stiff clay has an undrained cohesion
of about 100 kPa. It follows from the above and a consideration of likely unit
weights of a normally consolidated clay, that the stiff consistency should be
encountered at a depth of about 60 m (or about 200 ft). From the above figure,
is about 0.1. From the foregoing, it would appear that this is equivalent to an
adhesion factor ( of about 0.5, which is indeed quite reasonable for a stiff clay.
In a similar manner, taking the undrained cohesion of a soft clay as about 25 kPa,
such a strength should be encountered at about 15 m (50 ft.). At this depth, is
about 0.2 and this corresponds to an value of about 1. Once again, the two
techniques seem to agree fairly well.
For granular soils, drainage is generally so fast that the distinction between
drained and undrained or effective and total stress analysis never arises. As with
piles in clays, the general pile capacity equation may be relatively simple but the
specific values used in the equation depend on a great number of factors and
are, as usual, difficult to evaluate.
Driven piles are a common solution in granular soils. The driving process
densifies and strengthens a granular soil with the result that good pile
performance can be expected.
Bored piles, on the other hand, cannot be constructed in granular soils without
some form of temporary or permanent support. Without such support, the sides
can collapse and the base may heave. Support can be provided by temporary or
permanent casing, which can be driven or screwed into the ground prior to drilling,
or by use of drilling muds (e.g. bentonite/polymer fluid) which apply a stabilizing
hydrostatic stress to the walls. The former solution can be expensive, the second
solution is cheaper, but requires considerable expertise, and is messy. It should
also be borne in mind that the remnant bentonite/polymer may reduce the
available shaft resistance for the pile.
For granular soils, it can be assumed that the cohesion is zero and therefore, the
generalized pile capacity equation reduces to (similar to the ESA for clay):
𝑃 𝑄 𝑄 𝐴 𝑓 𝐴 𝑓 𝐴 𝑃′ 𝑁 𝐴 𝐾𝑃′ 𝑡𝑎𝑛𝛿
For shallow foundations (D/B < 1), it is normally assumed that Nq is a function of
’ only (Terzaghi and Peck, 1967). However, as D/B increases above 1, it is
generally accepted that the Terzaghi and Peck bearing capacity factors become
increasingly conservative (c.f. the effect of D/B on Nc for piles in clay). There
have been a number of suggestions for the Nq proposed over the years. Some
of these, and the methods for their implementation, will be discussed hereafter.
Meyerhof (1951, 1976) proposed bearing capacity factors, which can be found in
most texts [Bowles (1997), Das & Sivakugan (2019)]. These bearing capacity
factors depend on the pile having reached a penetration ratio (D/B) greater than
a critical value, (D/Bcrit) which is a function of the friction angle . Graphs of the
critical penetration ratio will be found together with Meyerhof’s Nq values. Typical
values for a range of friction angles are:
It is important to note that these values of are the in-situ undisturbed values.
Some of the older methods of design have suggested that significant reductions
should be made to the measured ' value when designing bored piles to account
for loosening effects. Although Meyerhof’s method does not explicitly allow for
differences between driven and bored piles, consideration should be given to
reducing the value of ' to between 0.5 and 0.7 if loosening of the soil is
considered likely.
If the critical embedment ratio is not attained, these values must be reduced
accordingly.
For layered soils, where the pile penetrates from a loose soil into a dense soil,
Meyerhof suggests that for all practical purposes, the full bearing capacity of the
denser soil can be assumed if the penetration exceeds 10 pile widths (D/B > 10).
Presuming that the total pile embedment still satisfies the critical ratios mentioned
before, Meyerhof recommends the following reduced unit point bearing pressure:
q b q l (l )
ql (d ) ql (l ) Dd ql (d )
10B
Where, qb is the unit base resistance;
ql(l) is the limiting unit base resistance in the loose sand based on the
appropriate values of Nq and for the loose sand;
ql(d) is the limiting unit base resistance in the dense sand based on the
appropriate values of Nq and for the dense sand;
Dd is the pile penetration into the dense sand, and
B is the pile width.
A number of methods have been proposed, but one of the simplest, and indeed
one that is accepted generally, has been proposed by Kishida (1967). This
method suggests that the value of ’ to be used for the determination of Nq for
driven piles is:
' 0 .5 ( ' 1 4 0 )
Where ’1 is the representative value of ' at the pile toe before pile installation.
Similarly, Kishida proposes that for bored piles, the value of ' to be used for the
determination of Nq is
1 3
N.B. Berezantev’s method is used in combination with limiting depth theory, and
effective overburden pressures reduced from their theoretical values. (See under
The value of P'ob).
The value of effective overburden pressure at the pile base would appear to be a
simple enough calculation, however, there is great debate in the geotechnical
community about the appropriate value. This section will describe some of the
alternative postulations.
In many of the older methods of pile design, it has often been assumed that the
vertical effective stress that is used in pile design in a granular soil increases
linearly with depth just as the vertical effective overburden pressure (Meyerhof,
1951 and 1976).
W.T.
zc
vc
Most recently, some researchers (Kulhawy, 1984) have returned to the idea of a
vertical effective stress that is not limited by some arbitrary critical depth. Others,
(Randolph, 1993) have postulated that a bulb of high stress follows the pile tip
down during driving, and that the shaft resistance behind the tip is progressively
fatigued. The jury is still out! This highlights that in geomechanics we are dealing
with the complex interaction of man-made objects and nature - this is ultimately
what makes geomechanics continually interesting and challenging!
Meyerhof (1951) proposed an alternative set of bearing capacity factors, N'c and
N'q for deep foundations. These are given graphically in Bowles (1988) or Das &
Sivakugan (2019). The method of Meyerhof allows full calculation of effective
overburden pressure at the base, in combination with the Meyerhof bearing
capacity factors for deep foundations (see Meyerhof’s method).
Note, however, that Meyerhof also proposed limiting values of end bearing which
could not be exceeded by the product of Nq and P'ob, i.e.
Qb Ab 50N q tan
Once the value of zc/B is known, the value of P'ob can be calculated. If the toe
of the pile is above this critical depth, then the value of P'ob is simply the vertical
effective overburden pressure. However, if the toe of the pile below this critical
depth, P'ob must be limited to the effective vertical overburden pressure at the
critical depth.
z
Pob for z z c
i.e.
P 'ob = z c for z z c
Vesic (1977) proposed a method for estimating the pile end bearing capacity
based on the theory of expansion of cavity (refer Das & Sivakugan, 2019 for
details). As per this method, the end bearing resistance in sand (fb) is calculated
as:
𝑓 𝜎′ 𝑁 ∗
Refer to Das & Sivakugan (2019) for 𝑁 ∗ values with friction angle (’) and
reduced rigidity index (Irr).
As noted previously, Kulhawy (1984) suggested that the limiting depth approach
to determination of pile capacity in sand was a fallacy which resulted from the
effects of overconsolidation on the capacity of piles of typical (<20m) length. What
is more, the limiting depth theory was based on the results of model pile tests
alone. He proposed that the full overburden pressure should be used, and the
unit tip resistance be computed as:
For D/B > 4~5, the first term becomes less than 10% of the second term, and can
be ignored. The dimensionless portions of both terms are given by Kulhawy in
graphical form against friction angle, . To use these graphs, it is necessary to
estimate the so-called rigidity index, Ir, which he approximates as follows:
30
in loose sand, I r
2
where q is in ton / ft
1/ 2 o o
q tan(28 to 26 )
110
in dense sand, I r
2
where q is in ton / ft
1/ 2 o o
q tan(40 to 37 )
As for the determination of effective overburden at the pile base P'ob, there is
dispute as to whether the full or limited overburden pressure should be used for
computing the effective overburden pressure for shaft resistance, P'os. The
following summarizes the different methods
Vesic (1967) use full calculated effective overburden pressure above critical
depth, and limited effective overburden pressure below the critical depth. (See
also previous).
In this method, the following post-installation friction angle, ', should be used
to compute the critical depth ratio, zc / d:
The parameters K and are sometimes treated separately, but in other design
methods are combined for convenience. This section will initially look at these
parameters individually, and then discuss the methods which use the combined
parameters K tan or K tan .
The value of K must be related to the initial horizontal earth pressure coefficient
(Ko = 1 - sin ), however, because of the very variable effects of pile installation,
the relationship is far from clear.
Tomlinson & Woodward (2008) recommends the following typical values for K /
Ko:
Pile Type K / Ko
The angle of friction between the pile shaft and the soil () has an upper limiting
value of '. Clearly, this upper limit would apply when the side of the pile is so
rough that the shearing occurs entirely within the soil. However, for pile shafts
which are smoother, shearing will occur at the interface with a friction angle <
'. It follows that is a function of the surface characteristics of the pile shaft. There
have been many suggestion as to the values of which should be used for
different surface finishes.
Some design methods (Meyerhof and Vesic among them) consider the two
factors together as the function (K tan but sometimes denoted K tan ). These
methods will be outlined in the following:
Meyerhof also proposed some simple guidelines for quick estimation of shaft
resistance for driven piles based on SPT N values :
K s tan a
on steel piles, and it is likely that the
values given in this figure are 2.0
conservative for other (rougher) surface
finishes.
1.5
The values of K tan given for Vesic’s
method should be compared with those
derived from Meyerhof's (1976) results 1.0
shown previously. It will be seen that the 28 33 38 43
K tan values for Vesic’s method are
about twice those of Meyerhof. Note,
however, that different overburden pressures apply. Poulos and Davis suggest
that Vesic’s method is preferred for driven piles, but recommend using Meyerhof’s
method for bored piles.
Bell & Robinson (2012) summaried the -values of different researchers for
driven and bored piles as given below:
N.B. For long piles in coarse-grained soils, average shaft friction is to be limited
110 kPa unless otherwise proven by suitable pile load testing or experience in
similar ground conditions.
The previous design methods apply to granular soils comprising silica sands -
these are the predominant on-shore sands. In the offshore environment,
calcareous sands predominate, and these behave entirely differently. Although
the angle of friction of calcareous sands prior to pile installation may be in excess
of 35 or even 45, the driving process tends to crush the sand particles
significantly. As a result, normal stresses at the interface are extremely low, with
the result that shaft resistances can be as low as 5 kPa.
There are many forms of in-situ tests, which can be used as part of the overall
design process for any particular type of foundation. Many of these tests are used
to provide conventional soil parameters either directly (e.g. the pressuremeter) or
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3A.35
Topic 3A: Axial Capacity of Piles
indirectly (e.g. SPT/CPT) so that the design methods discussed above can be
used.
Poulos (1989) proposed a direct relationship between the unit ultimate shaft
resistance of a pile with the SPT blow count as:
𝑓 𝛼 𝛽𝑁 𝑖𝑛 𝑘𝑃𝑎
Where, fs = unit ultimate shaft friction; N = SPT value, and α and β are constants
depending upon soil and pile type (refer to the table below).
Poulos (1989) also proposed a direct relationship between SPT N value and unit
ultimate end bearing as:
𝑓 𝐾𝑁 𝑖𝑛 𝑀𝑃𝑎
Where, K=constant (refer to table below), values depend on pile type and soil.
Meyerhof (1976) proposed direct correlations between SPT N-value and unit
shaft/base resitances of a pile as:
The use of CPT for direct correlations with shaft/base resistances of driven piles
is discussed below. It is to be noted that the use of CPT data for bored piles is
highly unreliable unless the CPT results have been adjusted to account for the
changed ground condition affected by the pile installation methods. As such direct
correlations between CPT and bored piles are avoided here.
It would appear that the base resistance of a driven pile in sand is directly related
to the cone resistance, qc.
i.e. Qb Ab qc
For a uniform sand, this would be a relatively simple method of design. The main
problem arises when dealing with 2
natural sands, which can show Cone resistance, kg/cm
100
substantial variations in qc over fairly
short distances. It therefore becomes
largely a scaling problem so that the
representative value of qc can be
Readings at 20cm
extrapolated from the size of the intervals
Idealised curve of
cone resistance
Depth , m
There is another important difference between a CPT and a full-size pile which
concerns the depth which a pile must be driven into a bearing layer to achieve
the resistance predicted from the penetrometer results.
Suppose, as is shown above, the penetrometer tip first enters a bearing point
layer at point A. Because of the scale of the device, it may only require 350 mm
(10 diameters) of further penetration before the resistance of the tip fully reflects
the increased bearing strength.
In the case of a pile of 350mm, the pile penetration necessary to develop the
layer's full capacity may still be 10 diameters, but this would now be 3.5m as
shown in the figure. Therefore, if the required pile resistance was 100 kg/cm2 (9.8
MPa), then the pile must be driven to the depth of point D (which is approaching
a further 3 m of penetration) rather than to some point just.
This procedure was first proposed by De Beer (1963) but seems to give some
very large penetration which according to Vesic (1970) are uneconomical. Vesic
suggests that driving the pile a further 5 diameters after the depth at which the
peak cone tip resistance was recorded is sufficient. The difference is probably
due to the greater densification caused by the full sized pile.
White & Bolton (2005) suggested a direct correlation between ultimate unit base
(fbu) resistance and qc as:
𝑓 0.9𝑞
Where, 𝑞 =weighted average qc value 8D above and 4D below pile tip (Dutch
method).
Once again, there are a great number of different methods for calculating the
shaft resistance of piles directly from CPT results; some based on sleeve friction
(fs) but most based on tip resistance (qc). High degree of caution is required
when using fs as the difference could be significant due to the soil consolidation
effect. An example of the former approach is provided by the following equation:
Qs = S1Asfs
Type of Pile S1
Timber 1.2
Parallel-sided concrete or steel;
flat end 0.6
pointed end 1.1
Driven cast-in-situ* 1.6
Open steel tube and H-pile 0.7
Begemann (1977) has suggested that the factor S1 should be taken as 0.7 for all
types of driven pile. It is also normal to have an upper limit for shaft resistance
which can vary, but according to Dutch practice (Meigh, 1987), is 0.12 MPa or
120 kPa.
As noted above, there are a number of methods for the calculation of shaft
resistance from cone tip resistance. Another method used is:
Qs = S 2 As qc
Type of Pile S2
Timber 0.012
Precast concrete 0.012
Precast concrete, enlarged base* 0.018
Steel displacement 0.012
Open-ended steel tube 0.008
It is important to check that the calculated allowable load does not exceed the
structural capacity of the pile shaft.
If driven or bored piles are driven or placed in rocks where a practical refusal is
expected then the structural capacity of the pile governs the capacity. In the U.K.,
a maximum allowable limit of 0.25f’c for concrete piles is considered when
calculating the pile capacity. This ensures pile materials are not overstressed
from driving. If steel pile is used then the driving stress needs to be limited to
values, which will not cause damage to the rock below the tip.
Goodman recommends for piles driven onto rock that the following equation
applies:
𝑄 𝑞 𝑁 1 𝐴
Nq = tan2 (45 + / 2)
is the drained angle of friction of the intact rock
Zhang (2010) introduced a method for determining the end bearing resistance
based on the unconfined strength of rock mass as:
.
𝑓 6.39 𝜎
where, 𝜎 is the rock mass unconfined compressive strength =()0.7c, c=UCS
of intact rock, E=0.0231RQD-1.32 ≥ 0.15; E = reduction factor to relate intact
rock UCS with rock mass UCS.
Zhang also summarised in the paper a number of methods available to relate the
UCS of intact rock with the UCS of rock mass based on various rock mass
classification systems (i.e., RMR, GSI, Q).
Where piles are driven or more likely drilled some distance into rock, an allowance
can be made for the resistance along the shaft. Socket resistance depends on
many factors such as the strength of the rock, socket roughness, interface bond,
socket contamination, weak rock smear presence etc. The following methods
relate to estimation of shaft resistance.
𝑓 0.05𝑞
This criterion should be used with caution for highly fractured rocks, for which
Poulos recommends a maximum allowable shaft resistance of 75 to 150 kPa.
Poulos also suggests a reduction in the recommendation of 30% for piles in uplift.
Bell and Robinson (2012) provided a table as given below, which summarises a
number of methods for the estimation of ultimate unit shaft resistance (fsu) of
bored pile based on average uniaxial comptessive strength of rock (quc) in MPa.
Williams & Pells (1981) method, widely used in Australia, requires values of rock
socket reduction factor () and rock mass correction factor () can be found using
the figures below if the uniaxial compressive strength of rock (quc) and mass factor
(j) are known from the laboratory and field investigation. Rock socket reduction
factor () for Horvath (1978) & Rosenberg & Journeaux (1976) methods can also
be found in this figure which are lower than the William & Pells values who
included the rock mass correction factor in their method.
Pells (1999) presented an elastic design method for the rock socketed pile based
on the intial work of Pells & Turner (1979). In this method, the compatibility of
deformation between the shaft and base was considered to develop a base load-
sharing chart incorporating the pile-rock stiffness ratio and the socket geometry.
The ultimate limit state and serviceability loads (AS1170) are required in addition
to the settlement limit for this method. A step by step design procedure for full
sockets in uniform rock was provided in Pells (1999). If non-uniform rocks are
encountered and pile shaft slip is considered then Pells suggested using of Rowe
& Armitage (1984) as illustrated in Pells (1999). Readers are advised to consult
Pells (1999), Pells & Turner (1979) and Rowe & Armitage (1986) for an in-depth
understanding of this elastic design method.
The settlement of the pile at the socket level can be estimated using the elastic
equation as:
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3A.45
Topic 3A: Axial Capacity of Piles
The Monash method has been developed by Seidel (1993). The method is
incorporated into the Windows package “Rocket”, which are used in Melbourne
for designing rock socket piles in mudstone rock.
Notes: Most of the methods presented above assumes the socket is clean and
the base is free from debris. If construction methods are poor and base cleanness
can be ensured then the pile capacity is to be determined based on the socket
resistance only with appropriate factor of safety-often this is expensive. Sockets
are also roughened to ensure a ductile pile-rock interface failure. Different
roughness categories used in construction can be found in Pells (1999). In large
projects, where the accuracy of pile capacity estimation is important, a finite
element modelling should be undertaken.
When piles are placed in a settling ground (soft underconsolidated soil/loose fill),
the soil around the pile moves down relative to the pile thus will induce a
downward axial force (known as negative skin friction or drag force) on the pile.
Poulos (2008) pointed out that this additional downward axial force will have no
impact on the ultimate limit state (ULS) analysis as the plunging failure of pile
requires the pile to move down relative to the suppounding ground. Therefore,
the negative skin friction will not generally reduce the ultimate geotechnical
capacity of the pile. However, Tomlinson & Woodward (2008) suggest to ignore
the contribution to the support provided by friction over the length of a pile in soft
clay.
The effect of the negative skin friction on the pile will cause additional settlement
and this has to be accounted for the serviceability limit state (SLS) design.
Sw = S1 + S2 + S3
Where, S1=elastic compression of the portion of the pile shaft in the settling zone
due to the applied load on the pile head (Pw)
S2=elastic compression of the portion of the pile shaft in the settling zone due to
the induced drag forces (Pneg) that can be calculated based on the piles in drained
soil concept (i.e., shaft resistabce piles in granular soil).
S1 and S2 can be computed from simple column compression theory and S3 can
be computed adequately using the elastic theory for the length of the pile
embedded in the stable zone due to applied pile head load (Pw) plus the maximum
drag force in the pile (PNmax).
where Rug2=ultimate geotechnical capacity of the pile in the stable zone below the
depth of soil settlement; FS2= Factor of safety for the portion of the pile in the
stable zone ≥ 1.25. Readers are advised to refer to Poulos (2008) for further
details.
In the previous sections, the axial capacity of a single pile was discussed. Single
piles can, for example, be used to support isolated column loads or under legs of
transmission structures. However, where column loads are particularly high,
groups of single piles may be needed. Structures with raft foundations may also
be supported on arrays of multiple piles, or pile groups.
In addition to the effect of stress overlap, the capacity of a pile may be influenced
by the installation of neighbouring piles. Installation of displacement piles will tend
to increase soil density and lateral stresses locally around a pile. Installation of a
non-displacement pile may have the opposite effect. However, such broad
generalizations cannot be relied upon. For instance, driving piles into sedimentary
rocks can sometimes lead to significant loss of end-bearing known as relaxation.
In any event, the pile group capacity may be less than the sum of the capacities
of the individual piles which make up the group
Pug n Pup
It is common to not allow for any increase in capacity due to densification effects.
However, pile group capacity losses are an effect which engineers must be
careful to account for. Pile group capacity loss is by convention calculated using
a pile group efficiency factor, .
Pug n Pup
The efficiency can vary considerably. Major factors affecting efficiency are soil
type, pile group size and the ratio of pile diameter to spacing.
For pile groups in clay, it appears both from the results of laboratory and full scale
testing that the end bearing component of resistance is largely unaffected by
group action with the result that the efficiency factor for this component is
effectively one, i.e.
end ,clay 10
.
In practice, since piles in clay usually derive most of their capacity through shaft
resistance, it is normal to apply the efficiency factor to the total pile resistance
including both shaft and base resistances, despite the foregoing efficiency factor
of 1.0 for base resistance.
A free standing pile-group is one in which the piles either extend above the soil
surface to a pile cap clear of the founding layer or one in which the pile cap rests
on ground that cannot be relied on to support any load. There are a number of
methods by which the capacity of such a group can be estimated.
The Feld rule is a simple empirical method for estimating an efficiency factor. The
calculated capacity of each individual pile is reduced by one sixteenth for each
neighbouring pile. For example, consider a 3 x 4 pile group as shown below.
A B B A
B C C B
A B B A
In this group, there are three types of pile. Type A (4 No.) is surrounded by 3
neighbours, type B (6 No.) by 5 neighbours, and type C (2 No.) by 8 neighbours.
It follows that the efficiency of the pile group is given by:
12 16 16 16
1 4 13 6 11 2 8 0.70
It should be noted that this method does not take into account the spacing of the
piles. This is clearly a deficiency of this method.
(n 1)m (m 1)n
1
90 mn
Using the example given above for s/d ratios from 2 to 10, the efficiency factor
varies as follows:
s/d d/s
2 0.50 0.58
3 0.33 0.71
4 0.25 0.78
5 0.20 0.82
6 0.17 0.85
8 0.13 0.89
10 0.10 0.91
The pile group is modelled as a shallow footing capacity, however, the footing
capacity, PBL, includes the cohesion at the side, which is not normally taken into
account for shallow footings, i.e.
Model pile tests conducted by Whitaker (1957) appear to confirm that as the pile
spacing decreases, the failure mechanism changes from one of individual pile
failure with slip lines developing around individual piles to one of block failure with
the slip lines around the perimeter of the pile group. However, this method
suggests that there is a critical spacing at which the failure criterion changes
abruptly from one to the other. As this does not appear to be the case it is
suggested that an empirical relationship can be used to give a more realistic
transition from one failure mode to the other. This relationship is given by
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3A.50
Topic 3A: Axial Capacity of Piles
1 1 1
2
2 2
2
Pug n Pup PBL
1 n 2Pup
2
1
2 2
PBL
Poulos and Davis (1980) provided an example of the application of this approach,
where the ultimate load of a pile group of specific dimensions is related to the
number of piles in the group. This is shown below. It can be seen that there is a
general transition between single
pile behaviour and block action. 5m
Furthermore as would be expected,
once block action begins to 5m
dominate, there is little advantage in
increasing the number of piles in the
group. In this particular case,
increasing the number of piles Clay
L = 20m d = 0.3m
beyond about 60 does little to c u = 50 kPa
30
Isolated
single pile
The general applicability of this failure
empirical transition approach was
examined by Poulos and Davis in 20
connection with the model pile tests
conducted by Whitaker. Some of
the comparisons are shown in the 10
following two figures for variations in
pile numbers and pile length to
diameter ratios. These figures
suggest that higher efficiencies can 0 5 10 15
be achieved for greater pile No. of piles
spacings, a smaller number of piles
in a group and for lower length to
diameter ratios.
1.0
0.6
1.0 1.0 1.0 1.0
0.2 52 group 72 group 92 group
1 2 3 4 1 2 3 4 1 2 3 4
s s s
d d d
0.5 0.5 0.5 0.5
0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8
s s s s
d d d d
L L L L
(a) d = 12 (b) d = 24 (c) d = 36 (d) d = 48
This type of pile group has a rigid pile cap which bears directly on soil which
provides reliable support. In this case, all piles must operate as a group because
of the rigid boundary condition at the pile heads.
Pcap cN c B c Lc BL
L
Pcap cN c Bc Lc nAp
Clearly, when the piles are relatively close together, the first mode will govern the
capacity of the group, but as the spacings of individual piles increases, the second
mechanism will predominate. This general behaviour for capped pile groups is also
confirmed by the model pile test results of Whitaker (1960) as shown in the following
figure.
of increase of efficiency
0.7
decreases. For the capped 5 2,7 2 ,and 92
groups, the rate of increase of
0.6
efficiency continues as before 2
5
because of the larger pile cap Tests on freestanding
0.5 groups
influence. 72 Tests on piled
foundations
0.4 Calculated for piled
2
foundations, assuming
9 block failure
0.3
1 2 3 4
Spacing factor, s/d
There is very little detailed information in the literature concerning the capacity of
pile groups in sand. For driven piles, it seems fairly certain that for sand that is less
than dense, the efficiency of a group will generally be greater than unity. It would
appear from various model tests (e.g. Vesic, 1969), that the maximum efficiency
ranges from about 1.3 to 2 at spacings in the general range from 2 to 3 pile widths.
For dense and very dense sand, the driving action could loosen the sand, but as it
is virtually impossible to drive piles through dense sand, this situation is unlikely to
have any practical significance. In the case of bored piles in granular soils, the shaft
capacity does not generally dominate design. Therefore, when considering this
relatively small component of total resistance, the use of an efficiency factor of one
should not cause any major overestimation of capacity.
For end bearing in piles in granular soils, it seems again that an efficiency factor of
unity should be used for both driven and bored piles. However, it must be stressed
in the case of bored piles, that some care should be taken in assessing individual
pile capacity so that any loosening during pile formation is allowed for.
It may be worth pointing out that with piles which are primarily end bearing (which
generally means piles founded in sand, gravel, rock, etc. although it is possible for
this to occur with piles in clay), it is important to check that there is a sufficient
thickness of the bearing stratum under the pile bases so that the bearing capacity is
not reduced significantly.
The underlying softer layer in the following schematic would not affect the capacity
of a single pile, but would be within the zone of influence of the pile group indicated.
Bearing Layer
It is important to check that the bearing capacity of the underlying layer is not
exceeded so that the group does not punch through the bearing layer as illustrated
below. Shallow footing analyses would be applied.
Test pile
Fill or
weak soil
Com pact
stratum
11. References
Bell, A. & Robinson, C. (2012). Chapter 54-Single piles. ICE manual of Geotechnical
Engineering. ICE publishing, London, doi: 10.1680/moge.57098.0803
Das, B.M. & Sivakugan, N. (2019). Principles of Foundation Engineering. 9th ed.
Cengage, Boston, USA.
Haberfield, C.M. & Lochaden, A.L.E. (2019). Analysis and design of axially loaded
piles in rock. Journal of Rock Mechanics and Geotechnical Engineering, 11(3), 535-
548.
Kulhawy, F. (1984). Limiting tip and side resistance: fact or fallacy? Analysis and
design of pile foundations. Proc. ASCE Geotech. [Link]. Symposium. Oct., 1984.
Ed. J.R. Meyer : 80 - 98.
Pells, PnJ (1999). State of practice for socketed piles in rocks. Proc. 8th ANZ Conf.,
Hobart, 1, 307-327.
Pells & Turner (1979). Elastic solutions for the design and analysis of rock-socketed
piles. Can. Geotech. J.
Poulos, H.G. & Davis, E.H. (1980). Pile foundation analysis and design. Wiley.
Poulos, H.G. (2008). A practical design approach for piles with negative friction. ICE
Proc. Geotechnical Engineering, 6, 19-27.
Tomlinson, M. & Woodward, J. (2008). Pile design & construction practice (5 ed.).
Taylor & Francis, NY, USA.
Rowe, K. R. & Armitage, H.H. (1987). A design method for drilled piers in soft rock.
Can. Geotech. J., 24, 126- 142.
Semple, R. M. & Rigden, W. J. (1984). Shaft capacity of driven pipe piles in clay,
Symposium on Analysis and Design of Pile Foundations, American Society of Civil
Engineers, San Francisco, 59–79 (cited from Tomlinson & Woodward, 2008).
N.B. Not all references are listed here. Readers are suggested to refer to standard
texts listed above for all other references not listed here.
TABLE OF CONTENTS
Until relatively recently, the settlement of a single pile was calculated either analytically
based on many simplifying assumptions or on an empirical basis through correlations
with other pile tests in similar situations. Of course, it may be noted that in many
situations, the settlement was not calculated at all but assumed to be satisfactory if
the load did not exceed one third of the ultimate load.
Of the analytical methods, one of the more popular was to assume that for an end
bearing pile, the pile had its entire load concentrated on its base and it could be
analysed in the same way as a shallow foundation. For a friction pile, the same
assumption was made but the equivalent base was assumed to be founded at a depth,
which was two thirds of the length of the pile. Clearly, these assumptions oversimplify
a much more complex situation and cannot be expected to yield accurate predictions
of pile settlement.
There are many examples of the empirical approach including Meyerhof (1959), who,
for piles in clay, suggested that for loads not exceeding one third of the ultimate load,
the settlement is given by
db
30F
where db is the diameter of the base of the pile, and
F is a factor of safety
MR col
where col is the compression of the pile as a free standing column, and
MR is the movement ratio.
For short piles rigid piles of col < 8 mm, MR = 1. For long compressible piles for
col > 8 mm, MR = 0.5.
As it is nearly always one of these three methods which is employed to estimate the
settlement of a single pile, each of these three methods will be considered in detail
below.
The load transfer method was originally suggested by Seed & Reese (1957) with
further refinements by Coyle & Reese (1966). The method is basically an iterative
technique, which is used to calculate the load/displacement characteristics at the head
of the pile and, thereby allow the construction of a full load/settlement curve.
The method requires the pile to be divided into a number of segments, which are
assumed to be connected by springs as shown. A small displacement is then applied
to the pile base, and by calculating the forces and displacements for each segment
progressively up the pile shaft, the load at the head of the pile can be determined along
with its vertical displacement. This process is repeated for a number of base
displacements until a sufficient range of pile loads and displacements are obtained to
construct a complete load/settlement curve.
and that conventional elastic analyses applies. For example, for a rigid circular
base,
2DEz nbot
Fnbot
1 2
It is of course possible to assume that the pile base displacement is zero, which is
appropriate for piles with only small loads or for piles on rigid bearing strata.
3. For the first trial, assume that znmid the displacement at the midpoint of segment
n, is equal to znbot.
4. Using the
adopted t-z
curve, for the
n mid
value of znmid,
Shear stress ,
Shear displacement
Fntop Fnbot
Fnmid
2
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3B.5
Topic 3B: Settlelment of Piles
7. Calculate the elastic deformation of the bottom half of the segment, n, as follows
Fnmid Fnbot Ln
z nbot 2 2
AE
where AE is the pile stiffness.
9. Compare the value of znmid assumed after step (3) with the value calculated after
step (8). If they are outside a specified tolerance, then assume znmid to have the
value calculated after step (8) and repeat steps (4) to (9) until the tolerance is
satisfied. Normally, this can be achieved after a few iterations.
10. Proceed to the next segment up the pile shaft, segment n -1. Repeat steps (2) to
(10) with the following load and displacement values at the bottom of segment n -1
Fn 1 bot Fntop
Fntop Fnmid Ln
z ntop 2 2
AE
11. Continue in this manner up through the pile segments until the load and
displacement at the head of the pile are obtained.
12. Steps (1) to (11) are then repeated for different values of displacement until a
range of load and displacement values are obtained for the head of the pile to
enable a plot of load against settlement.
It is clear from this iterative process that the method is ideally suited for computers.
However, it would be fair to say that the method can only be regarded as approximate
and does suffer from three main disadvantages.
The first is that the method assumes that there is no significant interaction
between successive pile segments and the soils surrounding these segments.
This is not realistic. However, the model has nevertheless been found to be
reasonably accurate.
Finally, the method is only as accurate as the assumed t-z curve, which can be
difficult and costly to obtain. There would seem to be some limited success with
large laboratory shear box methods (Johnston et al., 1987).
Despite the above limitations, Reese, Houdan and Vijayvergiya (1969) have
suggested that for bored piles in clay, the load transfer curve can be generalised by
az a max 2
so so
A very useful form of the load transfer method has been developed by Randolph
(1986) and the program is called RATZ. This program is capable of many refined tasks
including the analysis of dynamic loading.
Probably the most successful forms of analysis for the settlement of piles foundations
is based on methods of elasticity in which Mindlin’s equations for the displacements
within a soil mass caused by loading within the mass are used. Several investigators
have used this approach but perhaps one of the most complete sets of solutions have
been developed by Poulos and Davis (1980). Vesic (1977) also presented a rigorous
method based on elastic analysis.
While there are arguments that pile behaviour is not necessarily totally elastic, the
solutions obtained do seem to capture the essential features of pile settlement and the
results obtained do seem to compare well with piles that have been instrumented. Pile
working loads are generally sufficiently below the pile ultimate capacity to ensure
essentially linear load/settlement response up to these loads. Therefore, the
application of elastic methods can be justified.
This technique involves the discretisation of the pile into a number of elements. It is
then necessary to obtain mathematical expressions for the vertical displacement of
Department of Civil Engineering, Monash University
Revised: 6/2020
CIV4249: Foundation Engineering 3B.7
Topic 3B: Settlelment of Piles
the pile and the soil at each element in terms of the unknown stresses on the pile. By
applying the compatibility conditions to the pile and the soil, the various equations can
be solved to yield the displacement at any given pile head load. Details of these
analytical methods are discussed in some detail in Poulos and Davis (1980).
One of the strong points of the various analyses that have been carried out is that
design charts for a wide range of pile conditions have been developed using
dimensionless parameters. These have made the analysis of pile settlements
comparatively easy for a wide range of commonly encountered conditions.
In the various analyses that follow, two variables are frequently used. The first is the
area ratio, RA, defined as
Ap
RA
As
where Ap is the area of the material of the pile section,
and As is the area bounded by the perimeter of the pile.
For a solid concrete pile, Ap = As, and it follows that RA = 1. However, for a steel pipe
pile, Ap < As, and therefore, the area ratio would be much less than one.
Ep
K RA
Es
where Ep is the Young’s modulus of the pile material,
and Es is the Young’s modulus of the soil.
It may sometimes be necessary to establish how much of the total pile load is
transferred to the base of the pile. While this is not directly connected to pile
settlement, it is given here as the analytical techniques used are the same as for the
settlement analysis. At working load, many piles carry very little of their load through
base resistance, even though much of their ultimate resistance may be derived though
end bearing. The longer and more flexible the pile, the higher will be the percentage
of load carried by shaft resistance at working load.
0.04 0.9
Values of K
0.8
10000
C
0.02 1000
0 10 20 30 40 50
0.7
L/d 100
0.6
0.5
0 0.1 0.2 0.3 0.4 0.5
s
The first of these charts also considers the effect of piles with enlarged bases. From
these figures, it is clear that the percentage of load reaching the pile base increases
considerably with enlarged bases.
20
Values of K Values of K
20 000 5000 20 000
5000
1000 1000
10
8
500 500
6
Cb
4
100
100
2
(a) L
d
= 75 (b) L = 50
d
1
1 10 100 1000 1 10 100 1000
Eb Eb
Es Es
10
10 000
8 1000 500
100
6
Values of K
4
20
Values of K Cb
20 000 1000
500
10 2
8 100 L
(d) = 10
d
6
Cb 1
1 10 100 1000
4 Eb
L
Es
(c) = 25
d
Values of K
2 5
500, 1 000, 5 000, 20
4 000
100
3
11 10 100 1000 Cb
Eb 2
Es
L
(e) =5
d
11 10 100 1000
Eb
Es
The values of Io, RK, RL and R are given in the following figures.
1.0 3
0.8
100 Values of L d
0.6
50
R
0.4 RK 25
2
Values of db /d
10
1
0.2
2
3
Io 5
2
0.1
0.08 1 1
10 100 1000 10 000
0.06 K
0.04 1.0
For L d = 100
I o = 0.0254
Values of L d
For 3 > d b d > 1 50
0.8
0.02 25
0 10 20 30 40 50 Rh 10
5
L 0.6
d
2
0.4 1
0.2
0
0 1 h 2
L 0.5 L 0
h
1 .0 0
K = 100
0 .9 0
500
R 2000
1000
0 .7 5
0 0 .1 0 .2 0 .3 0 .4 0 .5
s
PI P
E sd
Pile -
where, I IoRK RbR
Young's modulus Ep
d
In addition to the definitions
given above for Io, RK and R, L
Soil -
Rb is a correction factor for the Young's modulus Es
stiffness of the stiffer layer and is Poisson's ratio s
given in the following set of
charts.
Bearing Stratum -
Young's modulus Eb
Poisson's ratio b
1.0 1.0
100 100
0.8 500 0.8
1000 500
0.6 0.6
Rb > 20 000 5000 Rb > 20 000 1000
0.4 0.4
Values of K Values of K 5000
0.2 0.2
(a) L d = 75 (b) L d = 50
0 0
1 10 E b 100 1000 1 10 E b 100 1000
Es Es
1.0 1.0
100
0.8 0.8
0.6 0.6
Rb 500
0.4
Rb 100
0.4 1000
1000 > 10 000
0.2 Values of K 0.2 L
500
L
(c) d = 25 > 2000 (d) d = 10
0 0
1 10 E b 100 1000 1 10 E b 100 1000
Es Es
1.0
Values of K
0.8
0.6
Rb
0.4
(e) L d = 5 100
0.2
5000 500 1000
20 000
0
1 10 Eb 100 1000
Es
This figure shows that for very long compressible piles, the presence of a stiffer layer
at the pile base has little effect.
An alternative method of solution has been proposed by Poulos and Davis to include
the concept of the movement ratio, MR, where
settlement of pile
MR
elastic shortening of pile
PL
and MR
EA
1.0 L =5
d
10
25
0.8
50
100
0.6
MR
0.4
PL
= M
0.2 EP A P R
s = 0.5
s = 0
0
100 200 500 1000 2000 5000 10 000
K
For a pile bearing on a rigid stratum, this figure gives the relevant values of MR. This
method is although relatively simple but it does not give as general a solution as was
presented earlier.
Poulos and Davis (1980) provide some illustrations of how slip between the pile and
the surrounding soil can influence pile settlement when elastic analyses are employed.
For the purpose of illustration presented here, only the relatively simple case of a
floating pile in a cohesive soil will be considered. A slip factor, Ms, is defined so that
1.0 1.0
1
> 2000
2
0.8 0.8
3
100
0.6 0.6
K = oo Values of K
Ms Ms 10
vs = 0.5
0.4 0.4 L = 0.5
d
L ca s = 0.5
Curve d cu ca = 1
100 1 cu
0.2 1 0.2
2 25 1
3 25 0.4
0 0
2.25 2.0 1.75 1.5 1.25 1.0 6 5 4 3 2 1
Factor of Safety Factor of Safety
These two figures show the effect of pile slip for a constant adhesion along the shaft
of the pile in terms of factor of safety for a few variations in L/d and ca/cu, and for
constant L/d and ca/cu ratios. It would appear from this that provided the pile is
relatively slender (L/d > 25), the influence of L/d is not very pronounced. It would
appear that to working loads of about half ultimate loads, the pile response is primarily
elastic therefore providing some justification of the use of elastic methods. As the ratio
ca/cu decreases, the effect of slip becomes more pronounced. A further point, which
is evident, is that the lower the value of pile stiffness, K, the more significant is the
effect of soil slip.
Based on these observations, it would appear that long steel piles in soft clay are
unlikely to significantly be influenced by soil slip and therefore are ideal for analysis by
elastic methods. On the other hand, relatively short, large diameter concrete piles in
stiff clay will be influenced by side slip as the side shear resistance peaks long before
the base resistance (Whitaker and Cooke, 1966). Therefore, for these pile types,
purely elastic analyses to calculate settlements at working loads may not be
appropriate. Clearly, for these same pile types cast in soft rock where the pile stiffness
is even lower, the load/deformation curve departs linearity even earlier (Williams, et
al. 1980). For such piles, an alternative method of calculating settlement after the
development of significant slip is presented below.
As was pointed out in the previous section, slip on the pile-soil interface has little
influence on piles, which are relatively long, with a fairly stiff section and in weaker
soils. It follows that the load/settlement curves for these piles will be relatively linear to
For concrete piles with low L/d ratios in stiff clay however, the shaft interface usually
experiences full slip long before the base resistance reaches its ultimate value. The
following method gives a means of estimating the load/settlement curve.
The first step involves calculating the settlement experience by the shaft. The load on
the shaft is estimated from
Ps P (1 )
where is estimated from the previous methods for floating piles or for piles bearing
on a stiffer layer.
If it assumed that the settlement is linear up to the full capacity of the shaft, Psu, then
the settlement of the pile head with respect to Ps can be calculated from
where
The load carried by the base, Pb, is calculated from the total load as
Pb P
Again, by assuming that the displacement of the base alone up to failure is linear, the
settlement of the base is given by
I Pb
E sd b
Where, Es is again the Young’ modulus along the pile shaft. It should be noted that if
the modulus below the base is different from around the pile shaft, then this is allowed
for by use of the appropriate elastic model.
However, this equation does not include the compression of the pile shaft for loads
greater than Psu. For a shaft load of Psu, the load at the head of the pile is given by
.
Ps
and the proportion of this which reaches the pile base is given by
1
Therefore, the additional compression of the pile shaft for loads on the shaft greater
than Psu, is given by
P L
Pb su
1 AE
I Pb P L
Pb su
E s d (1 ) AE
Psu
Pys
1
I
ys Pys
E sd
Total
Pu
Base
Load
Pbu
Pys
Shaft
P su
ys u
Settlement
It is now possible to construct the load/settlement curve for the shaft as shown in this
figure, as well as the first leg of the full load/settlement curve.
It is now necessary to calculate the load/settlement curve for the base. This requires
calculation of the ultimate base load, then using
I Pbu P L
u Pbu su
E sd 1 AE
the settlement of the base, u, at which the base ultimate is reached can be calculated.
This then allows the base load/settlement curve to be constructed as shown above.
By adding the two curves together, the peak for the full load/ settlement curve is at Pu
= Pbu + Psu at a settlement of u.
Vesic (1977) proposed a semiempirical method based on elastic analysis where the
total elastic pile head settlement (s) has been divided into three components: (i) elastic
shortening of the pile (se1), (ii) settlement (se2) due to the working load at the pile point
(Qwp), and (iii) settelemnt (se3) due to the working load along the pile shaft (Qws).
Where, 𝑠
The settlement of the pile point (se2) due to the working load (Qwp) can be calculated
from the elastic equation presented for the shallow foundation case previously.
A semi-empirical method for calculating se2 was also proposed by Vesic (1977) as:
The elastic settlement of the pile (se3) due to the average working load along the shaft
was given by:
A semi-empirical method for calculating se3 was also proposed by Vesic (1977) as:
If single piles are placed in overconsolidated clay then the consolidation component
of the settlement could be small and often neglected. However, for group piles in clay,
the consolidation settlement of the group is to be estimated separately. This can be
done using the consildation settlement concepts discussed previously for shallow
footing.
In the foregoing sections, it is clear that the critical soil parameter required for the
estimation of pile settlement is the modulus of the soil. The accuracy with which this is
achieved will determine how accurate the predictions of settlement are. The available
methods for determining modulus are many and varied ranging from small laboratory
sized tests on samples obtained as undisturbed as possible from the field to large
scale in-situ tests at the site on which construction is to take place. Clearly, the choice
of which method should be used in a particular situation will be governed by a number
of variables, one of which is cost. Some general discussions on the relative methods
of a number of approaches are presented in Johnston (1983).
3𝐸′
𝐸
2 1 𝜇
It is generally considered that moduli obtained from small-scale tests in the laboratory
are not advisable for soils, which exhibit considerable disturbance on sampling and
soils that contain significant fabric can cause major scale effects. This basically means
that all granular soils and cohesive soils containing fissures, joints, lenses, inclusions,
etc. are not really suitable for small-scale laboratory tests. These types of soil may be
better to test in-situ. The best form of test is a pile-loading test, which is back-analysed
by elastic methods to obtain a modulus. This modulus can then be used to predict the
behaviour of piles of different sizes with a reasonable degree of accuracy. These tests
are however, not done without some cost.
Also, Young’s modulus of sands for bored piles is very difficult to estimate due to the
effect of installation method affecting the density of sand.
Poisson’s ratio of sands can vary between 0.25 and 0.35. An average value of 0.30
can be used when no test results are available.
3. Numerical Methods
Numerical techniques such as the finite element method and the boundary element
method have been developed for the analyses of pile behaviour by many investigators.
While these methods are relatively complex to use, they do provide facilities for the
realistic analyses of much more complex problems than can be handled by the load
transfer method or by methods based on elastic theory. However, as even a partial
description of the techniques that are available for piles would require significantly
more space than is available, it is suggested that more specific texts on the subject
are consulted.
Individual piles are surrounded by a zone of influence - stressed soil, which accepts
the loads transferred from the pile shaft and base. If the pile load is sufficiently small,
it can be assumed that there is no slip between pile and soil, and that the settlement
can be computed using elastic analysis methods.
As discussed with reference to the axial capacity of piles, when the piles are placed in
a group, the individual zones of influence overlap, and the same soil may be subjected
to stresses from many piles. The net result of this effect is that the capacity of an
individual pile in a group will reduce relative to its single pile capacity (There may be
other factors, such as densification of sands in the group, which may counteract this
effect, and even result in a capacity increase).
For the same reasons, when piles are in a group, the settlement effects of an individual
pile will tend to be additive, and the settlement of a pile group will therefore be greater.
The settlement of group piles depends on such factors as the number of piles in the
group, whether the piles carry their settlement by shaft friction or end bearing, and the
spacing/diameter ratio.
There have been a number of approaches proposed for determination of pile group
settlement, and these will be discussed hereafter. Empirical methods are perhaps the
least reliable, and the shallow footing analogy is overly simplistic. Elastic methods of
analysis are perhaps the most suitable method for the majority of engineering
problems. For more critical structures, finite element or boundary element analyses
may be more appropriate.
It is also important to note that the zone of influence of a pile group will extend
considerably further below the pile bases than the zone of influence of the individual
piles. If conditions below the pile group base improve, the net increase in pile group
settlement will be reduced. However, if a soft layer lies within the extended group zone
of influence, special consideration must be taken to include the increased settlements
due to compression of that layer.
Empirical approaches to pile group settlement analysis are based on the results of a
limited number of pile group load tests, either in field or laboratory situations. The test
data is particularly limited because of the difficulty and cost involved in test loading pile
groups.
empirical relationships such as the one suggested for driven piles in sand by Meyerhof
(1959):
s (5 s3 )
g p
(1 r1 ) 2
It is normally assumed that the pile group acts as a block, the dimensions of which are
defined by the perimeter of the group and a depth which is equal to either the depth of
the group for piles which are principally end bearing or to two thirds the depth of the
group for friction piles as shown below.
2 D
3
D
However, this method fails to take into account the way in which piles transfer their
load to the ground, the compressibility of the piles and the very important interaction
between the piles in a group. The fundamental assumption of block behaviour must
be questioned, because as discussed with respect to the axial capacity of piles in a
group, depending on the particular soil conditions and group geometry, the piles may
preferentially fail individually. In such a case, the block analogy would be
inappropriate. This method is at best approximate, and may give misleading
predictions.
The use of elastic analysis methods for pile group settlement has been aided
significantly by the advent of high-power computers since the 1960s. Although some
of these methods are relatively simple and can be performed manually, the majority
make use of computers, which can perform the numerical iterations necessary to
obtain a converged solution. Poulos and Davis (1980) have obtained solutions for a
vast range of boundary conditions, and have presented the results in dimensionless
chart form for general design use. The use of these charts will be discussed in the
following section.
The methods used for this form of analysis make use of Mindlin's equations for a
loaded and buried cylinder, and a loaded and buried base in an elastic halfspace.
However, for a pile group, it is necessary to consider the interaction between adjacent
piles and this can be done by integration of the Mindlin equations for adjacent piles.
By developing expressions for vertical displacements in terms of the unknown
stresses, and by applying the compatibility conditions, the displacements can be
obtained by solving the resultant equations. The results have initially been expressed
in terms of an interaction factor, , which can be defined as :
It should be noted that for all the solutions that are to be considered in this section, the
piles are considered as a free standing group with the pile cap located above the
principal founding layer and not contributing to the overall performance of the group.
The reason for this is simply that with a capped pile group, the interaction between the
piles and the cap acting as a raft is extremely complex and would probably be better
solved using numerical methods.
The following graph shows that as the spacing of the piles increases, so the interaction
factor decreases. Conversely, as the pile stiffness, K, increases so does the
interaction factor. Studies carried out on piles of different L/d ratio shows that
increasing values of L/d will cause the interaction factor to increase.
1.0
L
d = 25
0.8
s = 0.5
Values of K
0.6
F oo
1000
0.4
500
100
0.2
10
0
0 1 2 3 4 5
s 0.2 0.15 0.1 0.05 0
d d
s
𝛼 𝛼 𝑁
Where, 𝛼 is the interaction factor for floating piles (values of 𝛼 for L/d=10, 50 and
100 can be found in Poulos & Davies, 1980).
1.0 oo
Values of h
5 L
0.8
2.5
0.6
Nh
1.5
0.4
s
L
0.2 h 1.2
L = 25
d
K = oo
s = 0.5
0
0 s 5
d
0.2 d 0
s
The graph shows that as the bases of the piles get closer to the rigid layer, the
interaction factor decreases.
The settlement of a floating pile group in a semi-infinite soil formation can be calculated
from
g Rs p
The values of Rs are presented in Table 6.2 (p.121) (Poulos & Davies, 1980) as a
function of the L/d and s/d ratios, the number of piles in a group and the pile stiffness,
K. This table shows that the Rs factor increases as the s/d ratio decreases, but it
decreases as the number of piles decreases. As the pile stiffness increases, the value
of Rs also increases. The Rs values given in Table 6.2 are for square pile groups.
However, it should be noted that the shape of the pile group has little influence on the
individual values of Rs, and therefore, this table may be used for other similar shapes
of groups with the same number of piles.
For groups containing more than 16 piles, Poulos and Davis (1980) point out that Rs
varies approximately linearly with the square root of the number of piles in the group.
Therefore, for pile groups containing more than 25 piles, the Rs factor may be
reasonably estimated by selecting R16 and R25 for the same L/d , s/d and K, and
applying the following expression
R s (R 25 R 16 )( n 5) R 25
It may be of some interest to examine the effect of increasing the number of piles in a
group and the effect this has on settlement. To do this it may be better to consider an
alternative reduction factor which is defined as follows
First, consider the situation of increasing the number of piles so that the piles have the
same s/d ratio. Clearly, this means that the width and breadth of the pile group is
increasing. A solution to this problem is shown below for an L/d ratio of 25.
1.0
Single pile
0.9 L = 25
d
h = oo
Group reduction factor , R G 0.8 L
Group s = 0.5
0.7
2
0.6
3
2
0.5 2
0.4 5
3x2 2
3
0.3
42
0.2
52
0.1
0
0 1 2 3 4 5 d
s s
d 0.2 0.1 0
It can be seen that such an increase will reduce the settlement of the group with the
load being shared between more piles. However, it should be pointed out that the
settlement reduction seems to become increasingly less effective as the pile numbers
are increased. This is particularly true of piles of relatively close spacing. For example,
consider a pile group of s/d ratio of 3. For a 2 x 2 group, RG is about 0.71. If the
number of piles is increased to give a 3 x 3 group (i.e. increasing the number of piles
from 4 to 9), the value of RG decreases to 0.58. This is not a particularly large
reduction in settlement for over twice as many piles. If the group was further increased
to 5 x 5 piles (giving an increase in the number of piles to 25), the RG value would
decrease to 0.44. It follows that if the number of piles were increased from 4 to 25, the
total settlement would be a little over half the settlement for 4 piles. This seems a fairly
high price to pay. It should also be borne in mind that such an increase in piles for the
same s/d ratio would mean that the pile cap would have to increase dramatically,
which when considering the finite size of a building, may not be at all practical.
1.0
0.9
L d = 25
0.8 Design curve for h = oo
driven pile groups L
Group reduction factor , RG in sand s = 0.5
0.7 (Skempton,1953) K = oo
Points used by Skempton
0.6 to derive design curve
0.5
Group
0.4
2
2
2
0.3 3
2
4 2
0.2 5
0.1
0
0 10 20 30 40
Group breadth / d
This figure shows that for a given breadth to diameter ratio, there is little reduction in
settlement beyond a certain number of piles. It would appear then that unless the
original pile spacing is greater than about 6, there is little to be gained by increasing
the number of piles supporting a given pile cap, especially when there are over about
25 piles in the group. This general conclusion seems to hold for a very wide range of
different pile conditions.
Where a rigid layer exists below the pile bases, the group settlement will be reduced
to some extent. The correction factor for settlement for some particular cases is given
in the following figure. It will be seen that the settlement reduces as the depth to the
rigid layer decreases. It also decreases as the size of the pile group increases. More
details of this situation may be found in Poulos (1968).
1.0
22
32
0.8 42
52
2
2
0.6
h s = 0.5
32
s d = 2.5
42
0.4
52
L = 25
d
0.2 L = 50
d
0
1 h 2
L L
0.5 h 0
The settlement of such a group of piles can be calculated using the following equation
in conjunction with an appropriate table of values of Rs below.
g Rs p
It will be noted that whereas for friction piles, an increase in pile stiffness leads to an
increase in Rs, for end bearing piles, an increase in stiffness leads to a reduction in
Rs. Where the bearing layer is non-rigid, Mattes and Poulos (1974) give methods of
calculating group settlement.
It should be noted that where the number of piles exceeds 25, the relationship given
previously may be used:
R s (R 25 R 16 )( n 5) R 25
considered for projects where settlement is extremely critical to the structure, and
where a highly sophisticated analysis can be justified.
It should be noted that the accuracy of these methods is critically dependent on the
failure criterion adopted, and particularly the constitutive laws that govern the pile-
soil interface. The reader is referred to the numerous texts and papers, which deal
specifically with the numerical methods.
References
Poulos, H.G. & Davis, E.H. (1980). Pile foundation analysis and design. Wiley.
Tomlinson, M. & Woodward, J. (2008). Pile design & construction practice. Taylor &
Francis, NY, USA.
TABLE OF CONTENTS
1. Introduction ...................................................................................... 2
2. Design methods ............................................................................... 2
2.1 Brinch Hansen’s method .................................................................... 4
2.2 Broms method ................................................................................. 10
2.3 Elastic method ................................................................................. 12
2.4 p-y method for lateral capacity ......................................................... 14
3. Laterally loaded piles in a group .................................................. 18
1. Introduction
In the previous chapter, the axial load capacity of piles was discussed. However, piles
are used to support lateral loads and moments in addition to the axial load. In general,
lateral loads and moments are caused by wind, earthquake, wave and earth
pressures. Similar to axially loaded piles, piles subjected to lateral loads are to be
designed based on both the ultimate limit state and serviceability limit state criteria.
2. Design methods
The lateral load capacity of a pile is derived from the passive resistance of the soil
surrounding it. The distribution of the soil reaction depends on the stiffness of the pile,
the stiffness of the soil and the fixity of the pile. If a pile head is free and the length to
width ratio is below 10-12 then the pile tends to show a rotational failure. The lateral
movement of the pile causes to the mobilisation of passive resistance at the front of
the pile and soil deforms elastically. However, with larger pile head displacements, the
soil in the front fails plastically and load transfers at greater depth. On the contrary, a
long free-head pile fails due to the insufficient moment capacity of the section, as the
large embedment depth of the pile does not allow the lower part to move or rotate. In
the case of a short fixed-head (capped) pile, the pile may fail by translation whereas
the failure occurs at cap level for the long fixed-head piles due to the development of
high bending stress. These observations ultimately lead to the development of
appropriate design methods by considering the short or long pile type with or without
the head restraint.
This section will first demonstrate how a pile can be classified as a short or long pile.
Soil-pile interaction is considered using stiffness factors (R for OC clay and T for
sands, NC clay) for this classification purpose.
For soils with constant modulus with depth (i.e., stiff OC clay), the stiffness factor “R”
is defined as:
For soils with linear increase of modulus (i.e., soft NC, sands) the stiffness factor “T”
is defined as:
Effect of saturation on nh is to
be considered.
The above stiffness factors can be used to classify the free-head piles in a constant
or linearly increasing soil modulus case as per the table below (Tomlinson &
Woodward, 2008).
Rigid (free-head) L ≤ 2T L ≤ 2R
The lateral load calculation methods can be classified in to two broad categories:
Brinch Hansen’s (1961) method is suitable for short piles in uniform and layered soils.
It can be used to calculate the ultimate lateral resistance of the pile based on lateral
earth pressure theory and generally simple. Also, the method can be used for longer
semi-rigid piles to obtain a first approximation of the stiffness and embedment length
for a more rigorous analysis.
The resistance of the rigid pile is taken as the algebraic sum of the passive resistances
above and below the point of rotation “X” in the figure below (Tomlinson & Woodward,
2008). The location of the point of rotation is determined by trial-and-error solution.
For the calculation of passive resistance, the pile is divided into “n” number of
horizontal elements of length (L/n), where L is the length of the pile.
The unit passive resistance at a depth “z” below the ground surface can be calculated
using the following equation:
𝑝 𝑝 𝐾 𝑐𝐾
where, poz=effective overburden pressure, c=cohesion, Kqz &
Kcz=passive pressure coefficients for frictional and cohesive components at depth “z”.
To find depth of the point of rotation “X”, take M=0 at the point of horizontal load
application.
𝐿 𝐿
𝑀 𝑝 𝑒 𝑧 𝐵 𝑝 𝑒 𝑧 𝐵 0
𝑛 𝑛
By the process of trial and error, determine the depth “z” at which M=0.
If the pile head carries a moment then it can be replaced by a force and eccentricity
(M=H.e).
For a fixed-head pile the equivalent height (e1) of a force H acting on free head is
0.5(e+zf) where, zf=depth from ground surface to point of virtual fixity. Tomlinson &
Woodward (2008) provided the following guideline for the depth of point of virtual fixity
(zf).
Example 1
Piles are divided n=6 segments ie. L/n=6/6=1m and 1m width of pile
For cu=120kPa and u=0 (short-term), find Kc values from figure at “z”
The passive pressure variations with depth are tabulated below and the
distributions are shown in figure.
z (m) 0 1 2 3 4 5 6
z/B 0 1.1 2.2 3.3 4.4 5.5 6.6
Kc 2.2 5.5 6.2 6.7 7 7.2 7.3
cuKc 264 660 744 804 840 864 876
∑MH=(462x1x4.5)+(702x1x5.5)+(774x1x6.5)+(822x1x7.5)-[(852x1x8.5)+(870x1x9.5)]
= +1629 kN-m/m width of pile
Department of Civil Engineering, Monash University
Created: 6/2020
CIV4249: Foundation Engineering 3C8
Topic 3C: Laterally loaded piles
Try x=3.9m
To find out Hu, take moment around the point of rotation “X” as:
∑MX=Hux7.9=(462x3.4)+(702x2.4)+(774x1.4)+(820.2x0.45)+(838.2x0.1x0.05)+(852x
0.6)+(870x1x1.6).
Distributions of passive resistances and forces with depth (a) undrained-short term (b)
drained-long term (Tomlinson & Woodward, 2008)
z (m) 0 1 2 3 4 5 6
z/B 0 1.1 2.2 3.3 4.4 5.5 6.6
Kc 5.8 16 20 23 26 27 28
By algebraically summing the moment around the load application point (similar to the
undrained calculation) and varying the x-value,
Assume pile is fixed at an arbitrary depth below ground surface (zf values
were given previously)
Then calculate deflection as for simple cantilever with a free or fixed head
Note: The ultimate lateral capacity (Hu) for a long pile depends on the pile bending
capacity rather than soil yielding due to the very large passive resistance. This is a
structural design and one has to follow the structural design guidelines (perhaps
discussed in CIV2226) for determining the ultimate lateral load capacity by knowing
the ultimate moment capacity (Mu) as follows.
Broms (1965) introduced his method for calculating the ultimate lateral load of short
(rigid) piles based on simplified earth pressure distributions along the length of the
pile. The ultimate capacity can be calculated for a single layer cohesive or frictional
soil but cannot be used for c- soils or multi-layered soils. The pressure distribution
along the length of the in cohesive and frictional soils can be found in standard texts
and Das & Sivakugan (2019).
By knowing the pile head conditions (free or fixed), short (rigid) or long (flexible) piles,
eccentricity, and type of soils (clay or sand), the ultimate lateral load capacity can be
determined using the design charts provided below (taken from Das & Sivakugan,
2019). The deflection of the pile head xz (z=0) under working load can be estimated
using Broms (1965) charts available in standard texts (Das & Sivakugan, 2019).
Design chart for short piles in (i) sand and (ii) clay (Das & Sivakugan, 2019).
Design chart for long piles in (i) sand and (ii) clay (Das & Sivakugan, 2019).
In general, the allowable lateral load capacity of a pile depends on the acceptable pile
head deflection. Matlock & Reese (1960) proposed a method based on the subgrade
reaction, where the simpler Winkler soil model was considered.
The values of different coefficients used in the above solutions can be found in the in
Das & Sivakugan (2019). Readers are advised to refer this text for further details.
Step 1: Establish a set of p-y curves for the soil-pile system (separate
procedure for sand and clay, follow standard texts)
Step 2: Calculate characteristic lengths T (linearly increasing modulus, sands)
or R (constant modulus, stiff clay) for a given or estimated nh or k.
Step 3: Determine deflection with depth based on T or R using subgrade
reaction methods
Step 4: Determine p with depth from the p-y curves based on the calculated
deflections
Step 5: Find k=p/y then nh=k/x or k1=k and finally T=(EI/nh)0.2 or R=(EI/k)0.25.
Step 6: Check T or R value with Step 2. If they match then fine, otherwise go
for second trial with k or nh from Step 5. Continue until assumed and
calculated values agree. Use final R or T for y and M.
Establishment of p-y curves (sands)
Example 3
p (kN/m)
pu at 8m= 110 kN/m L= 20 m 2m
60
1m
0.5m
Construction of p-y curves p/pu=0.5(y/y50)^0.33 use this eqn to find p/pu ratio relevant to y/y50 40
y/y50 y(m) p/pu p(2m) p(4m) p(6m) p(1m) p(0.5m)
0 0.000 0.00 0 0 0 0 0 20
0.25 0.008 0.32 21 31 35 16 14
0.5 0.015 0.40 27 39 44 21 18 8y50
Select 0
Es, kPa
1 0.031 0.50 34 49 55 26 22 0.000 0.050 0.100 0.150 0.200 0.250 0.300 0.350
intervals 0 500 1000 1500 200
2 0.061 0.63 42 62 69 33 28 y (m) 0
4 0.122 0.79 53 77 87 41 35 trial 1
1
8 0.244 0.99 67 97 109 52 44 trial 2
10 0.305 0.99 67 97 109 52 44 2
1st TRIAL, assume nh= 500 kN/m^3 T= 3.9 m y=Ay(PT^3/EI)+By(MT^2)/EI General eng for deflection 3
Depth (m)
Range (nh) for soft NC clay=278-556 kN/m^3 Zmax=L/T 5.2 >4 Long pile 4
(Davisson & Prakash, 1963) y=yA+yB=(PT^3/EI)Ay+(MT^2/EI)By
5
TRIAL 1 From Table (refer books)-Reese & Matlock p-y curve kN/m/m
depth, m Z=L/T Ay By yA (mm) yB (mm) y=yA+yB depth, m y,m p (kN/m) Es=p/y 6
0 0 2.435 1.623 130 45 175 0 0.17 0 0 7
0.5 0.13 2.21 1.37 118 38 156 0.5 0.16 38 244 8 y = 0.0044x
1 0.26 2.03 1.215 108 33 142 1 0.14 41 289 y = 0.0054x
9
2 0.52 1.644 0.874 88 24 112 2 0.11 51 456
4 1.03 0.962 0.364 51 10 61 4 0.06 60 978 kN/m^3 m
6 1.55 0.46 0.81 25 22 47 6 0.05 62 1323 Trial nh T
8 2.07 0.142 -0.07 8 -2 6 8 0.01 10 1770 1 500 3.9
10 2.58 0.0335 -0.0795 2 -2 0 10 0.00 2 227 4.5
T= 4.5 m 3 200 4.7 close to 4.5
TRIAL 2 use nh = 227 kN/m^3 p-y curve kN/m/m
depth, m Z=L/T Ay By yA (mm) yB (mm) y=yA+yB depth, m y,m p (kn/m) Es=p/y
0 0 2.435 1.623 209 61 270 0 0.27 0 0
y (m) and M (kN-m)
0.5 0.11 2.26 1.43 194 54 248 0.5 0.25 44 178
0 500 1000 1500 2000 2500
1 0.22 2.07 1.263 177 48 225 1 0.23 50 222 0
2 0.44 1.7 0.935 146 35 181 2 0.18 59 326
4 0.88 1.08 0.44 93 17 109 4 0.11 75 687
6 1.32 0.64 0.16 55 6 61 6 0.06 65 1067
8 1.76 0.24 -0.03 21 -1 19 8 0.02 30 1543 Deflection (y)
10 2.21 0.0335 -0.0795 3 -3 0 10 0.00 5 Moment (M)
TRIAL FINAL use nh = 200 kN/m^3 T= 4.7 m kN-m p-y curve kN/m/m
depth, m Z=L/T Ay By yA (mm) yB (mm) y=yA+yB Am Bm M depth, m y,m p (kn/m) Es=p/y
Depth, m
(DRAFT release) Please check the formulas and output prior confirmation of the whole calculation.
References
Tomlinson, M. & Woodward, J. (2008). Pile design & construction practice. Taylor &
Francis, NY, USA.
Prakash, S. & Sharma, H.D (1990). Pile foundations in engineering practice. John
Wiley & Sons Inc.