Title: Numerical Modelling of Small Groups of Stone Columns
Title: Numerical Modelling of Small Groups of Stone Columns
Title: Numerical Modelling of Small Groups of Stone Columns
Publication 2012-09-21
Date
Downloaded 2019-09-25T01:34:53Z
Some rights reserved. For more information, please see the item record link above.
NUMERICAL MODELLING OF SMALL
GROUPS OF STONE COLUMNS
by
Micheál Killeen B.E.
June 2012
I hereby declare that this thesis has not been submitted in whole or in part to any other
University as an exercise for a degree. I further declare that, except where reference is
given in the text, it is entirely my own work.
I would like to express my sincere appreciation and deep gratitude to my supervisor, Dr.
Bryan McCabe, who devoted so much time and effort to this research project. I am extremely
grateful for all his support, encouragement and advice throughout my studies. Most of all, I
appreciate his friendship.
I would also like to extend a special thanks to my academic mentors, Dr. Annette Harte and
Prof. Padraic O' Donoghue, for their contributions throughout the duration of this research.
I would like to acknowledge Keller Foundations (UK) for their generous financial support
provided to this research project. I am especially grateful to Dr. Alan Bell and Dr. Derek
Egan for their highly practical advice and valuable contributions to many technical problems.
This thesis on the behaviour of small groups of stone columns was conducted in collaboration
with research partners. I would like to thank Dr. V. Sivakumar (Queen's University, Belfast),
Dr. Jonathan Black (University of Sheffield, UK) and also Dr. Minna Karstunen and Dr.
Daniela Kamrat-Pietraszewska (Strathclyde University, Glasgow) for their assistance to this
research project. I would also like to thank Dr. Jorge Castro and Dr. Bostjan Pulko for their
help in patiently answering my many questions. I would also like to acknowledge the advice
from the PLAXIS support team.
I greatly appreciate the financial scholarship provided by the Irish Research Council for
Science, Engineering and Technology (IRCSET) to support this research project. I would also
like to express my sincere gratitude to the Geotechnical Society of Ireland for funding my trip
to Egypt to participate in the International Young Geotechnical Engineers conference in
2009. I thoroughly enjoyed this experience and it has encouraged me to pursue future studies
in geotechnical research.
I am indebted to the staff in the Department of Civil Engineering at NUI Galway who
through their dedication and enthusiasm helped to instil a passion for engineering in me
during my undergraduate studies.
I am extremely grateful to my fellow postgraduate students, all my friends and my relatives
(including those who have passed away) for their support, advice and warm friendship they
have provided to me.
Finally, my last and most important acknowledgement goes to my parents and sisters. I am
forever indebted to you for the never-ending support provided to me and all the sacrifices
made.
Míle buíochas
Abstract
Vibro stone columns, installed using the vibro replacement technique, are a cost-
effective form of ground improvement for enhancing the bearing capacity and settlement
characteristics of various soil types. Large groups of stone columns (such as used to
support embankments) are conventionally modelled using the unit cell concept, which is
based on an infinite grid of columns supporting an infinitely wide load area. Therefore
all columns are equally confined on all sides and are subject to a constant increment of
vertical stress with depth. The behaviour of small groups of stone columns supporting
small area footings is quite complex as peripheral columns are subject to a loss of lateral
confinement and the increment of vertical stress decays sharply with depth. This research
is the first comprehensive three-dimensional numerical study of the factors affecting
both the mechanisms of load transfer from columns to soil and the settlement
performance of small column groups at working loads.
PLAXIS 3D Foundation is used for this research in conjunction with the advanced
elasto-plastic Hardening Soil model, which is adopted for the parent soil and stone
backfill. The soft soil modelled is that at Bothkennar, Scotland, the former UK
geotechnical test bed. The influence of key design parameters such as area ratio, column
length, stiffness, strength and installation effects upon the settlement performance and
deformational behaviour of small groups of stone columns was investigated through a
total of 45 numerical analyses.
The modelling has shown that the area ratio and column length have a significant
influence on the settlement performance of stone columns. Moreover, they appear to be
inter-dependent as the effect of column length becomes more pronounced at low area
ratios. The influence of column confinement (i.e. increasing number of columns) was
also found to have a positive influence on the settlement performance of small groups of
stone columns. It was also shown that the influence of key design parameters upon the
settlement performance of stone columns is dependent upon the mode of deformation.
New parameters called compression and punching ratios were defined to help identify
three distinct mechanisms referred to as "punching", "block failure" and "bulging". The
occurrence of these mechanisms was verified by analysing the distribution of total shear
strain within columns and the surrounding soil and also examining the variation of stress
and strain along the length of columns. It was found that area ratio and column length,
rather than the number of columns, dictates the load transfer mechanism for small groups
of stone columns.
A more in-depth analysis of the deformational behaviour reveals that some combination
of punching and bulging occurs simultaneously, with one particular mode of deformation
more influential for a given area ratio and column length. This is consistent with the
finding that settlement improvement factors increase with column length for all
configurations of columns and suggests that a unique critical length, as proposed by
previous laboratory studies, does not exist for small groups of stone columns. The
presence of a stiff crust, an important feature of soft soil stratigraphy not captured in the
laboratory tests, was shown to have a significant influence upon the deformational
behaviour of columns. The observation of a critical length from laboratory studies is
shown in part to be due to the absence of a stiff crust (i.e. homogeneous soil samples) as
columns are more susceptible to bulging in the upper layers and thus cannot transfer the
applied load to their base.
The stress concentration ratios at the ground surface were also examined and it was
found that they are related to the mode of deformation. Moreover, it was shown that
stress concentration ratios vary considerably with column position and as such, do not
uniquely reflect the settlement performance of stone columns. Instead, the stress
concentration ratios with depth were noted and it was observed that they are constant
with depth in the yielded sections of the column and decrease towards unity at the base
of floating columns thereafter.
The numerical output in this thesis has been developed into a simplified design method
which allows the settlement of a column group to be related to that of a unit cell with
knowledge of the footing to column length ratio and the column length to layer thickness
ratio (and thereby caters for floating column groups).
Table of Contents
List of Figures i
List of Tables viii
Nomenclature ix
1 Introduction
1.1 Background 1
1.2 Aims and objectives 2
1.3 Methodology 2
1.4 Outline of thesis 3
2 Literature review
2.1 Background 5
2.1.1 Introduction 5
2.1.2 Construction processes 6
2.1.3 Applications of vibro replacement 8
2.1.4 Limitations of vibro replacement 10
2.2 General deformational behaviour of stone columns 11
2.2.1 Behaviour of single columns 11
2.2.2 Behaviour of column groups 13
2.2.3 Summary of general deformational behaviour of stone columns 18
2.3 Settlement performance of stone columns 18
2.3.1 Case studies 19
2.3.2 Laboratory studies 20
2.3.3 Numerical studies 22
2.3.4 Summary of settlement performance of stone columns 28
2.4 Column installation effects 29
2.4.1 Laboratory investigations 29
2.4.2 Observed field measurements 30
2.4.3 Simulation of installation effects in numerical models 32
2.4.4 Summary of column installation effects 35
2.5 Design methods 36
2.5.1 Unit cell concept 36
2.5.2 Cylindrical Cavity Expansion Theory (CCET) 36
2.5.3 Ultimate bearing capacity 38
2.5.4 Magnitude of settlement 40
2.5.5 Rate of settlement 46
2.5.6 Summary of design methods 48
2.6 Summary of literature review 49
References 224
List of Figures i
List of Figures
2 Literature review
Figure 2.1 Range of soil types treatable by vibro compaction and vibro replacement
techniques (courtesy of Keller Ground Engineering) 6
Figure 2.2 Depth vibrator and principle of vibro compaction (Sondermann & Wehr, 2004) 7
Figure 2.3 Dry bottom feed system using ‘vibrocat’ 8
Figure 2.4 Area ratios for (i) square grids, (ii) triangular girds and (iii) pad footings 9
Figure 2.5 Development of stress concentration ratios in stiffer stone columns 9
Figure 2.6 (a) Vertical displacement within the column against depth and (b) Radial
displacement at the edge of the column / initial column radius against depth
(Hughes & Withers, 1974) 12
Figure 2.7 Typical test setup examined by Narasimha Rao et al. (1992) 12
Figure 2.8 Photographs of deformed sand columns exhumed at the end of footing
penetration (Muir Wood et al., 2000) 14
Figure 2.9 Photographs of sand columns beneath strip footings at beginning, middle and end
of foundation loading process: (a) L/d = 6; (b) L/d = 10 (McKelvey et al., 2004) 16
Figure 2.10 Deformed group of (a) short columns and (b) long slender columns with
additional column outside the footing, upper part, u = 10 mm, clay (white) and
column (grey) (Wehr, 2004) 17
Figure 2.11 Settlement improvement factor against area replacement ratio for sites with
widespread loading (McCabe et al., 2009) 19
Figure 2.12 Predicted versus measured settlement improvement factors for all widespread
loadings and footings (McCabe et al., 2009) 20
Figure 2.13 Comparison of Ks for isolated and group formation (Black, 2006) 21
Figure 2.14 Illustration of block failure in group columns (Black, 2006) 21
Figure 2.15 Settlement variation with (a) EP/ES ratio and (b) degree of penetration of piles
(Balaam et al., 1977) 23
Figure 2.16 Measurement of stress concentration below an embankment (Kirsch &
Sondermann, 2003) 25
Figure 2.17 Settlement correction factor versus size ratio (a) case of layered soil, (b) case of
10.8 m thick soft clay layer and (c) case of 30 m thick soft clay layer (Elshazly et
al., 2008a) 27
Figure 2.18 Soil displacement due to pile installation (Yu, 2000) 29
Figure 2.19 Factor of restraint measured during the installation of stone columns (Kirsch,
2006) 31
Figure 2.20 Development of ground stiffness during column installation (Kirsch, 2006) 31
Figure 2.21 Zones of influence for square and triangular arrangements of stone columns 36
Figure 2.22 Equilibrium of soil element (Powrie, 2004) 37
List of Figures ii
Figure 2.23 Relationship between allowable vertical stress on stone column and undrained
shear strength (Thorburn & MacVicar, 1968) 39
Figure 2.24 Settlement diagram for stone columns in uniform soft clay (Greenwood, 1970) 40
Figure 2.25 Design curves for basic settlement improvement factor n0 with Poisson’s ratio υ =
⅓ (Priebe, 1995) 42
Figure 2.26 Settlement of single footings on groups of columns (Priebe, 1995) 43
Figure 2.27 Boundary conditions for solutions A & B proposed by Balaam & Booker (1981) 44
Figure 2.28 Dissipation of excess pore water pressure (Han & Ye, 2001) 47
Figure 2.29 Stress concentration factor. Influence of radial deformation and plastic strain
(Castro & Sagaseta, 2009) 48
Figure 4.4 Comparison of field load test data with numerical simulations for columns spaced
on (i) 1.2 m × 1.5 m, (ii) 1.75 m × 1.75 m and (iii) 2.1 m × 2.1 m 73
Figure 4.5 Bothkennar test site location (Nash et al., 1992a) 75
Figure 4.6 Geotechnical profile at Bothkennar test site (Nash et al., 1992a) 75
Figure 4.7 (a) Yield stress and (b) yield stress ratio from incremental load consolidation tests
(Nash et al., 1992a) 76
Figure 4.8 Variation of (a) lateral total stress (b) K0 with depth (Nash et al., 1992a) 77
Figure 4.9 Variation of (a) compression index Cc and (b) initial voids ratio e0 with depth
(Nash et al., 1992b) 78
Figure 4.10 Variation of load with time (Jardine et al., 1995) 80
Figure 4.11 Load-displacement behaviour for pad footings (Jardine et al., 1995) 81
Figure 4.12 Characteristic stiffness-strain behaviour of soil with typical strain ranges for
laboratory tests and structures (after Atkinson & Sallfors, 1991) 82
Figure 4.13 Settlement profile of Bothkennar test site for wide area loading, calculated using
one-dimensional compression theory and PLAXIS 3D Foundation 85
Figure 4.14 Comparison of settlement improvement factors for an infinite grid of end-bearing
columns with field data collated by McCabe et al. (2009) 86
Figure 4.15 Variation of normalised oedometric moduli with vertical effective stress for
Bothkennar test site 89
Figure 4.16 Distribution of minor principal stress σ3 in end-bearing stone column and
surrounding soil (A/Ac = 3.5; L = 13.9 m; Infinite grid of columns) 90
Figure 4.17 Comparison of settlement improvement factors for an infinite grid of end-bearing
columns with analytical design methods 91
Figure 5.10 Normalised settlement improvement factors for small groups of end-bearing
columns 102
Figure 5.11 Normalised settlement for small groups of end-bearing columns 102
Figure 5.12 Variation of settlement improvement factors with (i) column length and (ii)
volume of stone for various arrangements of stone columns 103
Figure 5.13 Influence of column position upon settlement improvement factors 104
Figure 5.14 Influence of column compressibility upon settlement improvement factors 105
Figure 5.15 Distribution of plastic points (Mohr-Coulomb) within a cross-section of (i) 4, (ii)
5 and (iii) 9 columns beneath a 3 m square footing (L = 6 m) 105
Figure 5.16 Influence of column strength upon settlement improvement factors 106
Figure 5.17 Influence of coefficient of lateral earth pressure upon settlement improvement
factors 107
Figure 5.18 Influence of the lower Carse clay and the stiff crust upon (i) settlement and (ii)
settlement improvement factors for (a) 4, (b) 5 and (c) 9 columns beneath a 3 m
square footing 109
Figure 5.19 Displacement at surface of footing, base of columns and soil surrounding base of
columns 110
Figure 5.20 Variation of punching ratios with column length for an infinite grid of stone
columns, with usoil measured in (i) a 1 m square zone and (ii) a square zone,
whose width is equal to the column spacing 111
Figure 5.21 Variation of (i) punching and (ii) compression ratios with column length for an
infinite grid of columns 112
Figure 5.22 Deformation of an infinite grid of floating stone columns at an area ratio A/AC =
14.1 for column lengths of (i) L = 1 m, (ii) L = 2 m and (iii) L = 3 m 113
Figure 5.23 Deformation of an infinite grid of 3 m long columns at area ratios of (i) 3.5, (ii)
8.0 and (iii) 14.1 113
Figure 5.24 Variation of (i) punching and (ii) compression ratios with column length for
single columns 115
Figure 5.25 Variation of (i) punching and (ii) compression ratios with column length for 2×2
groups of columns 115
Figure 5.26 Variation of (i) punching and (ii) compression ratios with column length for a 3×3
group of columns spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1 117
Figure 5.27 Variation of (i) punching and (ii) compression ratios with column length for a 4×4
group of columns spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1 118
Figure 5.28 Deformed shapes from PLAXIS 3D Foundation for a 3×3 group of columns at
area ratios of (i) 3.5, (ii) 8.0 and (iii) 14.1 for column lengths (a) L = 3 m and (b)
L=6m 119
Figure 5.29 Influence of column confinement upon (i) punching and (ii) compression ratios
for various lengths and configurations of columns spaced at area ratios of (a) 3.5,
(b) 8.0 and (c) 14.1 121
Figure 5.30 Variation of (i) punching and (ii) compression ratios with column length for
various configurations of stone columns beneath a 3 m square footing 123
List of Figures v
Figure 5.31 Variation of (i) punching and (ii) compression ratios with column length for
various column spacings beneath a 3 m square footing 124
Figure 5.32 Influence of column compressibility upon the variation of (i) punching and (ii)
compression ratios with column length for a group of (a) 4, (b) 5 and (c) 9
columns beneath a 3 m square footing 125
Figure 5.33 Influence of column strength upon the variation of (i) punching and (ii)
compression ratios with column length for a group of (a) 4, (b) 5 and (c) 9
columns beneath a 3 m square footing 126
Figure 5.34 Influence of column installation effects upon the variation of (i) punching and (ii)
compression ratios with column length for a group of (a) 4, (b) 5 and (c) 9
columns beneath a 3 m square footing 129
Figure 5.35 Influence of stiff crust upon the variation of (i) punching and (ii) compression
ratios with column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3
m square footing 130
Figure 5.36 Determination of stress concentration ratios for (i) field/laboratory studies and (ii)
numerical/analytical studies 131
Figure 5.37 Stress concentration ratio determined by (i) averaging stress over entire soil area
and (ii) averaging stress within the zone of influence for each column 132
Figure 5.38 Distribution of vertical effective stress (σ'vert) immediately beneath a 4.5 m square
footing, which is (i) unreinforced (L = 0 m) and (ii) reinforced with end-bearing
stone columns (L = 13.9 m) 132
Figure 5.39 Variation of stress concentration ratio with column length for an infinite grid of
columns 133
Figure 5.40 Variation of stress concentration ratios with column length for (i) single columns
and (ii) 2×2 groups of stone columns 134
Figure 5.41 Variation of stress concentration ratio with column length for a 3×3 group of
columns, spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1 136
Figure 5.42 Variation of stress concentration ratio with column length for a 4×4 group of
columns, spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1 137
Figure 5.43 Variation of stress concentration ratio for small groups and infinite grids of
columns spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1 140
Figure 5.44 Vertical stress beneath an unreinforced footing for a (i) single column, (ii) 2×2,
(iii) 3×3 and (iv) 4×4 groups of columns at an area ratio of 3.5 141
Figure 5.45 Variation of stress concentration ratio with column length for a group of 4, 5, 7
and 9 columns beneath a 3 m square footing 143
Figure 5.46 Variation of stress concentration ratio with column length for various columns
positions beneath a 3 m square footing 145
Figure 5.47 Vertical effective stress beneath a 3 m square footing for columns spaced at (i) 1.0
m, (ii) 1.5 m and (iii) 2.0 m 145
Figure 5.48 Influence of column compressibility upon the variation of stress concentration
ratios with column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3
m square footing 147
List of Figures vi
Figure 5.49 Influence of column strength upon the variation of stress concentration ratios with
column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square
footing 149
Figure 5.50 Influence of column installation effects upon the variation of stress concentration
ratios with column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3
m square footing 150
Figure 5.51 Influence of stiff crust upon the variation of stress concentration ratios with
column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square
footing 152
Figure 6.15 Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration
ratio with depth for 13.9 m long columns, within a 4×4 group for A/AC of (a) 3.5,
(b) 8.0 and (c) 14.1 186
List of Tables
2 Literature review
Table 2.1 Increase in K0 to account for stone column installation effects by various authors 33
Table 2.2 Final solution for stresses and strains in column and soil (Balaam & Booker,
1981) 44
Nomenclature
Unless otherwise stated, the following abbreviations and symbols are used in this thesis.
However, when referring to specific publications, the original notation has been used.
Abbreviation:
CCET Cylindrical Cavity Expansion Theory
FEA Finite Element Analysis
FEM Finite Element Method
Symbols:
A Tributary area of soil per column in a large grid
AC Area of the stone column
A/AC Area ratio
B Width of square footing
CC, CS Compression and swelling indices, respectively
Ck Hydraulic change index (= ∆e/∆log(k))
D Diameter of circular footing
E Young's modulus
E50 Secant Young's modulus at 50% deviatoric stress
Eur Young's modulus for unload-reload
Ecol/Esoil Modular ratio
G Shear modulus
H Thickness of soil deposit
IL Liquidity index
IP Plasticity index
K0 Coefficient of lateral earth pressure at rest
KA, KP Coefficient of active and passive earth pressure, respectively
L Column length
N Diameter ratio (= de/dc)
Rf Failure ratio (Hardening Soil model)
Rinter Strength reduction factor (Hardening Soil model)
a Column radius
b Unit cell radius
c Cohesion
cu Undrained shear strength
Nomenclature x
Chapter 1
Introduction
1.1 Background
Vibro replacement is a popular form of ground improvement which is used to enhance
the settlement and bearing capacity characteristics of soft soils. Vertical columns of
compact stone are formed in the ground using either the top or bottom feed systems. The
vibro replacement process typically involves replacing 10–35% of the in situ soil with
crushed gravel. The high stiffness properties of the crushed gravel reduce the overall and
differential settlements of treated soil. The consolidation time is also reduced due to the
high permeability of the crushed gravel. Stone columns can be used to provide support
for a variety of loading scenarios ranging from small footings (i.e. pad/strip footings) to
wide area loadings (i.e. embankments and large floor slabs).
The design of foundations on soft soils is usually governed by settlement criteria rather
than bearing capacity due to their high compressibility. The majority of analytical design
methods developed to date contain many simplifying assumptions, such as the unit cell
concept, which assumes an infinite grid of columns supporting an infinitely wide load
area. Therefore, they do not account for the loss of lateral confinement associated with
groups of columns supporting small footings. In addition, the reduction in vertical stress
with depth beneath small footings is much sharper than that beneath wide area loadings
and, therefore, offers the possibility of partial depth treatment. While some correction
factors exist to account for the loss of lateral confinement, current design methods do not
consider the loss of lateral confinement for small groups of floating stone columns.
While the bearing capacity of small groups of stone columns has been well researched, a
dearth of information exists regarding the settlement performance of small groups of
columns. This is highlighted by McCabe et al. (2009) whose settlement database of over
20 case studies comprised only three case studies related to small groups of columns.
Similarly, the majority of numerical studies conducted are axisymmetric analyses on
Introduction 2
large groups of columns. The lack of information regarding the settlement performance
of small groups of stone columns was identified by Black (2006) who conducted some
high quality laboratory research; however, it is difficult to extrapolate the findings due to
scale effects associated with laboratory tests and also as some of the area ratios
considered are at the high end of typical values used in practice.
1.3 Methodology
A series of three-dimensional finite element analyses are proposed to investigate the
settlement performance and deformational behaviour of small groups of stone columns in soft
soil. PLAXIS 3D Foundation is a FE program which is specially developed for geotechnical
engineering and is adopted for this thesis. The stress-strain behaviour of the soil and stone
columns is simulated using advanced constitutive models. The soil profile adopted is that
of the well characterised Bothkennar test site, which consists of soft uniform clay overlain by
a stiff crust. The adopted soil profile and material parameters are validated and the analysis is
extended to study the influence of key design parameters such as area ratio, column length,
stiffness, strength, installation effects and the presence of a stiff crust.
Introduction 3
Two parameters are introduced to study the deformational behaviour of stone columns,
column compression and column punching ratios. These ratios are determined from the
displacement at the top and bottom of columns and are a simple method to investigate the
influence of various design parameters upon the mode of deformation of columns.
A more detailed analysis of load-transfer mechanisms for small groups of stone columns is
conducted by examining the distribution of total shear strains within stone columns and the
surrounding soil. This analysis is conducted for a select number of column lengths, which are
specifically chosen to display each mode of deformation. The effect of column spacing and
confinement upon column-column interaction and column-soil interaction for each mode of
deformation is identified.
The distribution of stress and strain along the column is also examined. These parameters
play a key role in the settlement performance of stone columns and allow a link between the
settlement performance and the deformational behaviour of columns to be established.
A background to the finite element analysis is presented in Chapter 3. The material models
used in PLAXIS 3D Foundation are described and the results of preliminary studies on mesh
sensitivity and boundary effects are presented. Other issues such as simulating the long term
behaviour of soft soil, modelling the column-soil interface and incorporating column
installation effects into the finite element model are discussed.
parameters are validated by simulating a load test reported by Jardine et al. (1995) at the test
site. The settlement performance of an infinite grid of columns at the Bothkennar test site is
also simulated with PLAXIS 3D Foundation and compared with one-dimensional
compression theory, a settlement improvement database and analytical design methods.
The results from the primary analysis of this thesis regarding the settlement performance,
deformational behaviour and stress concentration ratios for small groups of columns are
presented and interpreted in Chapter 5. The interaction between column confinement, length
and area ratio and their influence upon the settlement performance of columns is examined.
Three modes of modes deformation are identified, which are referred to as "punching",
"block failure" and "bulging", at typical working loads levels. Finally, the relationship
between the stress concentration ratio at the surface and the mode of deformation is
examined.
The modes of deformation identified in Chapter 5 are verified by examining the distribution
of total shear strain, normal strain and stress concentration ratios within columns for specific
column lengths in Chapter 6. The distribution of shear strains within columns and the
surrounding soil allows the column-column and column-soil interactions to be investigated.
The characteristic behaviour of individual columns (i.e. corner, edge and centre columns) is
examined though the distribution of vertical and horizontal strains along the length of
columns. The variation of stress concentration ratio with depth is also examined and its
relationship with the mode of deformation is highlighted.
A simplified design method which accounts for the loss of lateral confinement and reduction
of vertical stress with depth for groups of columns supporting small loaded areas is developed
in Chapter 7. The findings from PLAXIS 3D Foundation are discussed in relation to previous
research and the existence of a critical length and the effect of the stiff crust are highlighted.
The influence of the mode of deformation upon the settlement performance of stone columns
is also discussed. Finally, the stress concentration ratios at the surface and at depth are
compared with a field and laboratory measurements and also analytical design methods.
Finally, the conclusions of the research are presented in Chapter 8. A better understanding of
small groups of stone columns is gained and the findings can be used by practitioners to
relate the settlement of small groups and infinite grids of columns.
5
Chapter 2
Literature review
2.1 Background
2.1.1 Introduction
Deep vibratory techniques such as vibro compaction and vibro replacement are forms of
ground improvement which are used to enhance the bearing capacity and settlement
characteristics of weak soils. The vibro compaction technique was developed in 1936 by the
Keller group and one of its first applications was to densify non-cohesive river-borne
granular soils (McCabe et al., 2007). The technique consists of lowering a vibrating poker
into the ground, which imparts horizontal forces to the surrounding soil and causes the soil
particles to re-arrange into a denser state. However, it was found that vibro compaction
reaches its technical and economic limits in saturated sands with high silt contents, as fine
particles attenuate the horizontal forces imparted by the vibrating poker (Sondermann &
Wehr, 2004).
The limitation of vibro compaction in cohesive soils (i.e. soils composed of a high proportion
of fine particles) was overcome in 1956 by the development of the vibro replacement
technique. As with the vibro compaction technique, a vibrating poker penetrates the soil to
form a borehole. The resulting borehole is then backfilled in stages with coarse aggregate,
which is compacted by re-lowering the poker. This process results in stone columns which
are tightly inter-locked with the surrounding soil. Stone columns can readily be formed up to
15 m and typically replace 10–35% of the in situ soil (McKelvey et al., 2004). Figure 2.1
illustrates how the vibro replacement technique extends the range of soils types treatable by
deep vibratory techniques. Vibro stone columns offer a cost-effective alternative to traditional
piled solutions when supporting moderately loaded residential, commercial and industrial
structures.
Other variations of vibro stone columns include geotextile-encased columns and vibro
concrete columns. Geotextile-encased stone columns are used in very weak soils where
Literature Review 6
insufficient lateral support is provided by the surrounding soil and excessive bulging of stone
columns occurs. The geotextile material develops tensile forces which constrain the column
and reduce bulging.
Figure 2.1 - Range of soil types treatable by vibro compaction and vibro replacement techniques (courtesy of
Keller Ground Engineering)
Figure 2.2 - Depth vibrator and principle of vibro compaction (Sondermann & Wehr, 2004)
crushed stone is attached to the vibrator and bends inwards at the tip. Therefore columns are
formed without removing the poker, which allows the poker to case the borehole thereby
preventing collapse. Once at the design depth the poker is raised by 100–500 mm, depending
on soil conditions, and crushed stone flows into the cavity through the supply tube. The poker
is then raised and lowered a number of times until satisfactory compaction of the stone
backfill is achieved. This cycle is repeated until a column of compact stone is formed up to
the ground surface.
1.5
= = =
# ×
(i) (ii) (iii)
Figure 2.4 - Area ratios for (i) square grids, (ii) triangular girds and (iii) pad footings
Vibro replacement primarily enhances the settlement performance of treated soils due to the
high stiffness of granular columns. If equal strain (i.e. uniform deformation) is assumed
beneath the base of footings, it can be shown from equilibrium that a larger proportion of the
applied load is taken by the stiffer stone columns (Figure 2.5). Therefore, stress concentration
ratios develop in stone columns. Field experience indicates that stress concentration ratios
typically range from 2.5–5.0 (Barksdale & Bachus, 1983). The stress concentrations which
occur in columns reduce the vertical stress on the surrounding soil which, consequently,
reduces the settlement.
A comprehensive review of over twenty case studies by McCabe et al. (2009) highlights the
effectiveness of vibro stone columns in enhancing the settlement performance of soft soils. In
addition to reducing the magnitude of settlement, vibro stone columns also homogenise
treated soil deposits and thus reduce differential settlement. Therefore, it is not necessary to
span ground beams between compaction points, as is required with piled solutions, and
shallow footings can be used to support the foundations.
Literature Review 10
Vibro stone columns also act as vertical drains due to the high permeability of the stone
backfill. This allows excess pore pressure to dissipate rapidly which decreases consolidation
time. However, the drainage capacity of stone columns may be reduced slightly by the
creation of a smear zone around the vibrating poker upon penetration.
Slope stability
The vibro replacement technique can be used to increase the stability of embankment side
slopes if stone columns are installed deeper than the rotational failure surface. Stone columns
have a higher shear strength than the surrounding soil and therefore increase the shearing
resistance along the failure surface. Moreover, stress concentrations ensure that a larger
proportion of the applied load is carried by the stronger column material. The effectiveness of
stone columns at increasing the stability of two full scale test embankments, founded in deep
deposits of very soft cohesive soil, was investigated by Munfakh et al. (1984). The authors
report that in addition to reducing the settlement and consolidation time, stone columns
increased the load-carrying capacity of in situ soil by approximately 50%.
Mitigation of liquefaction
Liquefaction is initiated when excess pore pressures generated from seismic motions increase
to equal the inter-granular stresses. As a result, soil loses all its shear strength and large
deformations follow. The vibro replacement technique mitigates against liquefaction by
densifying in situ soils and acting as vertical drains. The initial poker penetration densifies
granular soils which increases inter-granular forces between soil particles. Stone columns
also act as vertical drains which allows for a rapid dissipation of excess pore pressure. Vibro
stone columns were employed to support a wastewater treatment plant in Santa Barbara,
USA, which was located in a seismically active zone affected by liquefaction (Mitchell &
Huber, 1985).
columns should not be installed in soils which contain layers of peat or decomposable
organics greater than 1d–2d (d = column diameter) in thickness (Barksdale & Bachus, 1983).
McCabe et al. (2009) state that considerable evidence exists to suggest that stone columns
can be formed successfully in soils with undrained shear strengths (cu) well below 15–20 kPa,
which is often quoted as the lower practical limit of the system. A small scale laboratory
study conducted by Wehr (2006a) indicates that columns can be installed in soil with cu ≥ 4
kPa.
The vibro replacement technique is also not appropriate in high sensitivity soils. The shear
strength of high sensitivity soils dramatically reduces once the grain structure is disturbed.
Therefore, the vibratory forces imparted from the poker during column installation may
remould the soil and, thus, significantly reduce the soil strength.
(a) (b)
Figure 2.6 - (a) Vertical displacement within the column against depth and (b) Radial displacement at the edge
of the column / initial column radius against depth (Hughes & Withers, 1974)
Figure 2.7 - Typical test setup examined by Narasimha Rao et al. (1992)
Literature Review 13
The findings from the laboratory tests were verified by a plane strain Finite Element Analysis
(FEA). Four modes of deformation were observed; bulging, shearing, punching and bending.
The occurrence of column bulging was found to be dependent on the lateral restraint provided
by neighbouring columns (A in Figure 2.8); the magnitude of bulging increases and the depth
of bulging becomes shallower with increasing area ratio (A/AC). Shearing of columns
occurred when low lateral restraint is provided to columns and is most noticeable beneath the
edge of footings near the ground surface (B in Figure 2.8). The zone of influence for the
footing may be defined by a cone of undeforming soil, located directly beneath the footing
and bound by the shear planes occurring in the soil. The depth of the zone of influence is
dependent on the mobilised angle of friction of the composite soil mass, which increases at
low A/AC. Punching of columns is evident in Figure 2.8(a) and occurs when columns are not
long enough to transfer the load to depth; this mode of deformation becomes more
pronounced for closely-spaced columns, i.e. low A/AC. Bending of stone columns occurs
when columns are subject to lateral loads and can be seen in Figure 2.8(c).
The findings from the laboratory tests also indicate that columns located at mid-radius (i.e.
between central and edge columns) are the most heavily loaded, followed by central columns
and then edge columns. The authors attribute the reduced load carrying ability of edge
columns to their decreased lateral restraint. However, the reason for the reduced performance
of central columns is less clear, as the stiffness of the soil increases towards the footing
centre. The authors also suggest that a critical length may exist; however, the dominant strain
in the column is dependent on the footing diameter (D) and not the column diameter. Wehr
(2004) simulated the laboratory tests with a plane strain FEA and suggests that a critical
column length exists at 1.5D, where a marked change in column behaviour occurs. Central
and edge columns shorter than this length exhibit punching and buckling failure, respectively.
Literature Review 14
The findings from the laboratory tests also indicate that columns located at mid-radius (i.e.
between central and edge columns) are the most heavily loaded, followed by central columns
and then edge columns. The authors attribute the reduced load carrying ability of edge
columns to their decreased lateral restraint. However, the reason for the reduced performance
of central columns is less clear, as the stiffness of the soil increases towards the footing
centre. The authors also suggest that a critical length may exist; however, the dominant strain
in the column is dependent on the footing diameter (D) and not the column diameter. Wehr
(2004) simulated the laboratory tests with a plane strain FEA and suggests that a critical
column length exists at 1.5D, where a marked change in column behaviour occurs. Central
and edge columns shorter than this length exhibit punching and buckling failure, respectively.
deeper. The authors conclude that the ultimate capacity of stone columns is influenced by
column spacing and diameter rather than column length and the granular mat.
(a) (b)
Figure 2.9 - Photographs of sand columns beneath strip footings at beginning, middle and end of foundation
loading process: (a) L/d = 6; (b) L/d = 10 (McKelvey et al., 2004)
Wehr (2004)
Wehr (2004) examined the behaviour of a single stone column and a group of stone columns
using a plane strain FEA, which employed an elasto-plastic constitutive law in a Cosserat
continuum. The results from the FEA for the single column and the group of columns are
compared with experimental studies by Brauns/Witt (1978) and Hu (1995), respectively. A
wedge of undeforming stone, which induces radial and vertical displacement, is observed
from the experimental tests on a single stone column. This failure mechanism, along with a
secondary failure mechanism in the form of a shear plane along the column-soil interface, is
identified from the FEA. A wedge shaped zone of undeforming soil, which is bounded by
shear planes, is also observed from the laboratory tests on a group of stone columns (Figure
2.10). The shear zones intersect beneath the centre of the footing and cause excessive bulging
Literature Review 17
to occur in central columns. In addition to inclined shear planes, the FEA identifies extra
shear zones along the interface of internal and external stone columns. The author suggests
that the shear zones between the internal columns may unite if the spacing of the stone
columns is reduced and the lateral deformation is limited. The unification of shear zones is
also dependent upon the relative movement of the column and soil. Another interesting
observation from the FEA is that the thickness of the inclined shear planes increases with
column length and that the ‘bulging’ observed by Hu (1995) in outer columns, is actually a
large shear zone. The influence of an edge row of columns, as investigated by Hu (1995), is
quantified using the FEM and its influence on bearing capacity was found to be negligible
(peak force increased by 7%).
(a) (b)
Figure 2.10 - Deformed group of (a) short columns and (b) long slender columns with additional column
outside the footing, upper part, u = 10 mm, clay (white) and column (grey) (Wehr, 2004)
Wehr (2006b)
Wehr (2006b) extended the previous analysis (Wehr, 2004) to study the effects of footing
flexibility. A wedge shaped zone of undeforming soil was not observed for a flexible footing
and columns tended to bulge rather than buckle. However failure mechanisms in the form of
a broad vertical shear zone at the edge of the footing and several approximately wedge
shaped and parallel shear zones were observed. The number of shear zones increased with the
flexibility of the footing. The shear zones also extend to a limited depth and are dependent on
the movement of the columns relative to the soil. Wehr (2006b) concluded that a flexible
footing has a better load-carrying capacity than a rigid footing due to the formation of wider
shear zones near the footing edge and the additional shear zones between columns.
Literature Review 18
It appears that column interaction influences the deformational behaviour of stone columns as
groups of stone columns exhibit bulging, punching, shearing and bending. The occurrence of
each mode of deformation was found to be dependent upon the geometrical properties of
columns and the footing, i.e. column spacing, length, position beneath footing and footing
diameter. The critical length for groups of columns was found to be dependent on the column
diameter and the footing diameter, with L = 1.5D and L = 6d suggested by Muir Wood et al.
(2000) and McKelvey et al. (2004), respectively.
The stress concentration ratio which develops in stone columns appears to be related to the
mode of deformation (McKelvey et al. 2004). Stress concentration ratios are low for columns
which are punching into the underlying soil and increase with column length to a maximum
when column bulging is occurring. It was also observed that the position of columns beneath
footings influences the magnitude of stress concentration ratio as edge columns carry a lower
proportion of the applied load at failure due to the loss of lateral confinement.
Figure 2.11 - Settlement improvement factor against area replacement ratio for sites with widespread loading
(McCabe et al., 2009)
The authors also examine the effectiveness of different construction techniques by comparing
measured and predicted settlement improvement factors in Figure 2.12. It appears that the
construction technique has a significant influence on settlement performance of stone
columns and also that the bottom feed system is the preferred choice for installing vibro stone
columns in soft soils. Some of the data points in Figure 2.12 suggest that vibro stone columns
are behaving worst than predicted, which may be due to a number of factors such as
workmanship and uncertainties in the measured data. The authors stress the importance of
workmanship in installing stone columns and state that the disturbance caused by imparting
excessive energy in the installation sequence may offset the benefit of a stone column.
Uncertainties also exist with the measured data, such as the methods used to predict
settlement and the ‘as-built’ dimensions of columns. On consideration of these facts and the
Literature Review 20
excellent performance of the bottom feed system, Figure 2.12 indicates that the assumed
angle of internal friction of 40° for the stone backfill is somewhat conservative.
Figure 2.12 - Predicted versus measured settlement improvement factors for all widespread loadings and
footings (McCabe et al., 2009)
Based on the initial slope of the 'load versus settlement' curve for both area ratios, it appears
that the group of stone columns are under-performing as the stiffness for groups is half of that
for single columns (Figure 2.13). This is attributed to a ‘block failure’ that manifests itself in
column groups, where all the columns act together and punch into the underlying soil, and
also to the reduced shoulder length between the edge of columns and the edge of the footing
(Figure 2.14).
The deformation patterns of columns were found to be dependent on both the column length
and arrangement. It was observed that both single columns and column groups that were 125
mm long, with a length to diameter ratio (L/d) ranging from 3 to 5, all punched into the
underlying soil. As the length was increased to 250 mm, it was observed that the single
column (L/d = 7–10) bulged whereas the group of columns (L/d = 11–14) continued to punch
into the underlying soil. However, column groups act as a 'block' and re-defining the L/d
ratios on the basis that the diameter for the group is more appropriate (i.e. d = 60 mm) leads
to revised L/d ratios of 4 and 6. Punching is precluded for end-bearing columns and it was
Literature Review 21
Figure 2.13 - Comparison of Ks for isolated and group formation (Black, 2006)
observed that all 400 mm long columns bulged. Black (2006) suggests that the critical L/d
ratio which defines the cross-over from punching to bulging failure is approximately 8 for
these drained analyses.
Another interesting finding from the tests is the relationship between the stress concentration
ratio (σcol/σsoil) and column length. No increase of vertical stress was observed in short
columns, as columns failing by end-bearing cannot sustain any extra vertical stress. However,
stress concentration ratios are observed to increase with column length, which reflects the
Literature Review 22
higher ultimate capacity of columns which are bulging. It was also observed for column
groups that the pressure recorded beneath the centre of the footings was higher than in the
columns. The increased pressure in the centre of the footing imparts a lateral force on the
surrounding columns, forcing them to bulge outwards and away from neighbouring columns.
This may explain the reduced settlement performance of groups of columns compared to
single columns, which are positioned beneath the centre of the footing.
The long term settlement performance of stone columns is estimated using the FEM and
accordingly both the column and soil are modelled as drained materials. The authors
examined the settlement performance of stone columns using an elastic and elasto-plastic
analysis. While the elastic analysis was found to slightly under-predict the settlement of the
column-soil unit at working loads (maximum discrepancy = 6%), it was deemed sufficiently
accurate and was adopted for the subsequent analyses. It is shown in Figure 2.15(a) that the
ratio of the column and soil moduli has a minor influence on the settlement performance of
stone columns. It is also noted by the authors that columns shorter than L/h = ¼ (h =
thickness of soil deposit) or columns spaced wider than de/d = 5 (A/AC = 25) yield negligible
settlement improvement factors (n). It can be seen in Figure 2.15(b) that the influence of
column length is more pronounced for closely-spaced columns. It is also interesting to note
the identical settlement performance of floating columns spaced at de/d = 2 and end-bearing
columns spaced at de/d = 5, which occupy 12.5 % and 4 % of the soil volume, respectively.
This emphasises the importance of column length upon the settlement performance of large
arrays of columns.
The rate of settlement is determined from a finite difference solution, which uses diffusion
theory to approximate consolidation theory developed by Biot (1941). The diffusion theory is
Literature Review 23
(a) (b)
Figure 2.15 - Settlement variation with (a) EP/ES ratio and (b) degree of penetration of piles (Balaam et al.,
1977)
The influence of the column spacing was examined by conducting analyses over a range of
area ratios A/AC = 3–10. As expected, settlement improvement factors increase (n =
1.2→2.0) at lower A/AC. However, no change in differential settlement was observed. Stone
Literature Review 24
columns act as vertical drains due to their high permeability and it was found that the rate of
consolidation increased at lower A/AC. The horizontal displacement at the column-soil
interface was also found to decrease with lower A/AC, which may be attributed to the
confining effects of neighbouring columns.
The influence of column compressibility was also investigated by varying the modular ratio
(λcol/λsoil) from 10 to 100, for columns spaced at A/AC = 5.3. It was observed that both
settlement improvement factors (n = 1.3→3.8) and differential settlement increase with
increasing column stiffness. The rate of consolidation is also influenced by the column
stiffness as the time taken to reach consolidation reduces for higher column stiffness, i.e. time
taken to reach δ = 0.9δmax (δ and δmax = average settlement and final settlement on the ground
surface, respectively) was 16 weeks and 10 weeks for λcol/λsoil = 20 and 100, respectively.
The stress concentration ratio was observed by Domingues et al. (2007b) to increase linearly
from 3.9 to 14.0 with increasing modular ratio.
study. The authors present the results of the development of stress concentration as a function
of applied loading for an embankment site in Kuala Lumpur. The curve presented in Figure
2.16 indicates that the stress concentration increases with applied load; this curve is specific
to columns installed on a square grid (A/AC = 4) to 14 m, in a 16 m deposit.
Figure 2.16 - Measurement of stress concentration below an embankment (Kirsch & Sondermann, 2003)
Elshazly et al. (2008a) develop modification factors to relate the settlement of finite groups
of columns (s) with the settlement of an infinite grid of columns (suc). The column length (L)
and area ratio (A/AC = 3.4) were both constant, while the number of columns beneath
footings was varied to yield normalised footing widths (B/L) from 0.5 to 4.7. Two soil
profiles were adopted for the investigation; the first soil profile is a layered estuarine deposit,
which is described in detail by Mitchell & Huber (1985), and the second soil profile is
developed from typical parameters for soft soil. The ratio of the Young’s modulus of the
stone column to the soil ranges from 1.3 to 2.6 for the first soil profile and is 8.5 for the
second soil profile. The installation effects associated with stone column installation are
incorporated by increasing the coefficient of lateral earth pressure (K0) of the soil
immediately surrounding the stone columns to 1.5 and 1.2 for the first and second soil
profiles, respectively. The loads are applied through a stone distribution blanket and the long
term settlements are determined using an axisymmetric FEA.
The authors identify two counteracting effects which are associated with the unit cell model.
The first effect concerns the distribution of vertical stress beneath pad footings and wide area
loadings. Vertical stress reduces sharply with depth beneath small footings and is negligible
beyond a depth of z = 2B. This is in contrast to the unit cell which stresses the soil profile to
the full depth and therefore increases settlement. The second effect is related to the loss of
lateral confinement that is experienced by exterior columns beneath small footings, which
increases the settlement. However, full lateral confinement is provided on all sides of
columns within an infinite grid, which results in enhanced lateral support and less settlement.
It can be seen in Figure 2.17(a) for a layered deposit that the settlement of small footings
(B/L < 2) is less than the settlement of a unit cell. Therefore it appears that the beneficial
Literature Review 27
effect of stress reduction with depth overrides the loss of lateral confinement for small
footings. However, as footing width increases (B/L > 2) it can be seen that the settlement of
small footings exceeds the settlement of a unit cell. This may be attributed to an increase in
the vertical stress in the soil with increasing footing width.
It can be seen in Figures 2.17(b) and 2.17(c) for the soft soil deposit that the settlement of
finite groups is always less than that of a unit cell. The increased footing settlement may be
attributed to the low lateral support provided by the soft soil deposit.
(a)
(b)
(c)
Figure 2.17 - Settlement correction factor versus size ratio (a) case of layered soil, (b) case of 10.8 m thick soft
clay layer and (c) case of 30 m thick soft clay layer (Elshazly et al., 2008a)
Literature Review 28
Kirsch (2008)
Kirsch (2008) conducted a FEA on two footings, which are supported by 25 columns.
Installation effects, which are based on a previous study by Kirsch (2006), are incorporated
into the FE model. The stress concentration ratio (σcol/σsoil) was measured in a sand layer
above the stone columns and was found to increase with load from 1.4 to 1.6 for central
columns and from 1.4 to 2.8 for corner columns. The measured stress-concentration ratios
increased upon reloading. The FEA is compared with field data and despite over-predicting
the ultimate capacity of columns it predicts the settlement behaviour quite well at working
loads. The influence of increasing the column length and the angle of internal friction is
found to be more pronounced at low A/AC. Another interesting finding from the study is that
column stiffness has a negligible influence on the settlement behaviour of the footing. Kirsch
(2008) suggests that columns undergo plastic deformations at relatively low loads due to
stress concentrations and are therefore insensitive to elastic stiffness parameters.
The settlement performance of small groups of stone columns was also reviewed. It appears
that the mode of deformation changes with the number of columns, as closely-spaced
columns act together and punch as a 'block' into the underlying soil. A comparison of small
groups of columns and infinite grids of columns highlights the importance of the lateral
support to columns and the distribution of vertical stress beneath pad footings in determining
Literature Review 29
the settlement of small groups of columns. The influence of column length and the angle of
internal frictionn for the stone columns w
were found to be more pronounced at low A/AC.
Stress concentration ratios were found to vary with applied load and also with position
beneath pad footings, as corner columns carry a higher proportion of the vertical load than
centre columns.
beyond the soil particles. The ultimate displacement pattern of the soil is similar to that of
cylindrical cavity expansion.
The use of a video camera to observe soil displacement surrounding the installation of a flat
bottomed penetrometer was adopted by Gill & Lehane (2001). This non-intrusive method
determines the soil displacement by tracking the movement of dark beads against a light
background, which consisted of a transparent artificial material. This artificial material has
properties similar to that of a lightly over-consolidated natural clay and consists of silica
particles and paraffin. The findings show that the soil is displaced both outwards and
downwards as the penetrometer tip advances towards the beads. Once the tip passes the
beads, the soil displacement becomes radial and upwards. It was also observed that negligible
displacement occurred once the penetrometer tip was more than five penetrometer diameters
deeper than the level of the beads. The results compare well with experimental data from
other authors, which reinforces the assumption by Baligh (1985) that the kinematics of deep
penetration are such that the displacement paths are relatively insensitive to the properties of
the soil.
Figure 2.19 - Factor of restraint measured during the installation of stone columns (Kirsch, 2006)
Figure 2.20 - Development of ground stiffness during column installation (Kirsch, 2006)
Castro (2007) recorded the development and dissipation of excess pore water pressure which
resulted from the installation of seven stone columns in a normally consolidated clay. The
pore water pressure rose upon poker penetration and peaked as the poker tip passes the level
of the piezometers. Despite significant heave occurring at the ground surface, it was deemed
reasonable to assume that plane strain conditions exist and cylindrical cavity expansion
theory may be used to simulate the poker installation. The increase in excess pore water
pressure recorded during the initial poker penetration agreed favourably with theoretical
values, which were predicted using an analytical formula derived by Randolph et al. (1979b).
An additional assumption was required that the undrained shear strength of the surrounding
soil iss reduced, due to the increase in excess pore water pressure, which is a common
phenomenon associated with driven pile installation. However, the agreement between field
and theoretical values breaks down with the installation of subsequent columns, as the
assumed boundary conditions of plane strain are violated. The dissipation of excess pore
water pressure occurred very quickly (within 15 minutes), which is 100 times less than the
Literature Review 32
theoretical dissipation time, predicted using a finite difference technique based upon theory
developed by Soderburg (1962). This is attributed to fractures in the clays, which result from
the high pressures due to column installation and act as drainage channels.
Gäb et al. (2007) measured the pore water pressure and settlement due to the construction of
a large embankment. The embankment was founded on 11 m of loose-medium dense
lacustrine sand which was underlain by 50 m of clayey silt. The stone columns were installed
to 14.5 m and on a triangular grid at A/AC = 7.7. In total 37 columns were installed in 4 rings,
with construction starting at the outer ring and progressively moving towards the centre. The
excess pore water pressure increased as the construction sequence moved closer to the
piezometers. This increase in excess pore water pressure was a maximum at 12 m, near the
column base and was even noticeable at 20 m depth. The excess pore water pressure
dissipated very quickly in the sand (< 1 day) and very slowly in the clay. Slight heave was
also measured at the site.
Egan et al. (2008) collate data relating to the installation of stone columns from various sites
including that of Castro (2007). The authors observe that heave may occur during the
installation of stone columns and tentatively suggest a relationship between column density
and the amount of heave. The authors suggest that heave is also a function of column size,
spacing and construction method. The footing arrangement also influences the amount of
heave as small groups and strips of stone columns do not generate as much heave as large
grids.
incorporate any change to soil properties due to column installation. The authors then
attribute any discrepancy in the load-settlement curve to a change in the stress state of the
surrounding soil and adjust K0 until the load-settlement curves match. Assuming full
confidence in the adopted soil parameters, an axisymmetric FEA indicates that K0 = 1.5
yields the best match between the actual and predicted load-settlement curves.
Elshazly et al. (2008b) extended to the previous analysis to investigate the effect of column
spacing upon the post-installation stress state of the soil. A range of column spacings were
analysed, which corresponded to A/AC ranging from 2.5 to 4.8. The back-calculated values
for the K0 ranged from 0.7 to 2.0, depending on the level of confidence in the soil parameters.
The best conservative estimates of K0 are 1.7, 1.2 and 0.85 for A/AC = 2.5, 3.7 and 4.8,
respectively. This indicates that installation stresses reduce at high A/AC (i.e. widely-spaced
columns).
Table 2.1 - Increase in K0 to account for stone column installation effects by various authors
Reference A/AC Coefficient of lateral earth pressure, K
Balaam & Booker (1977) 4 – 100 1.00
Barksdale & Bachus (1983) 4 – 10 0.75
Mitchell & Huber (1985) 2.0 – 4.9 1.00
Elshazly et al. (2006) 3.4 1.50
Domingues et al. (2007a) 3.3 – 10.0 0.70
Elshazly et al. (2008b) 2.0 – 4.9 0.85-1.70
Both authors conducted numerical investigations into the influence of this technique upon the
stiffness development of the surrounding soil and zone of influence of the expanding column.
The increase in soil modulus is estimated from the change in mean effective stress, using an
experimental relationship developed by Biarez et al. (1998). It is postulated by Brinkgreve &
Broere (2006) that soil stiffness is directly proportional to the mean effective stress for soft
soils. The undrained cavity expansion induces large excess pore water pressures in the
surrounding soil, which are magnified at low A/AC. Following consolidation of the soil, the
mean effective stress and hence the soil stiffness are increased. Debats et al. (2003) report a
30% and 40% increase in soil stiffness within a cylindrical zone of radius 2d in the
surrounding soil for column spacings of 6 m and 10 m, respectively. However, these findings
are specific to the Mohr-Coulomb model and even higher increases in soil stiffness in a
greater zone of influence are observed when the Hardening Soil model is adopted. Guetif et
al. (2007) also observe that a stress concentration occurs in the stone column and that the
lateral coefficient of earth pressure is increased above unity in the surrounding soil.
Kirsch (2006) investigated the influence of both the individual and the global column
installation effects upon the settlement improvement factor for a 7.2 m square footing,
supported by 25 columns. The results indicate that increasing the radial expansion of the
individual column to 8% results in a 45% increase in the settlement improvement factor. An
examination of the global installation effects suggests that increasing the stiffness of this zone
(2d–5d from the centre-line of the outer row columns) by three times its initial value yields a
25% increase in the settlement improvement factor. However, it is observed that the
Literature Review 35
additional improvement is highly dependent upon the loading stage, and that at high loads the
additional improvement gained drops once stone columns yield and plastic deformation
becomes dominant. Therefore, the installation effects do not influence the ultimate behaviour
of stone columns but they do play a positive role in enhancing the settlement behaviour.
Field experience indicates that horizontal stress and soil stiffness increase within a zone
which ranges from 4d–8d from the centre-line of columns. However, increases in horizontal
stress and soil stiffness are offset by the effects of remoulding and dynamic excitation within
4d. It was also found that an increase in excess pore water pressure due to initial poker
penetration agrees favourably with the analytical formula developed by Randolph et al.
(1979b). Column installation effects tend to improve settlement behaviour of stone columns.
However, the additional improvement drops once stone columns yield and plastic
deformation becomes dominant. Therefore, the installation effects do not influence the
ultimate behaviour of stone columns.
Figure 2.21 - Zones of influence for square and triangular arrangements of stone columns
The radial (∆εr) and tangential (∆εθ) strain are defined as:
Δ = −
=−
(2.1)
Δ = − =−
(2.2)
Literature Review 37
Assuming the soil can be idealised as a uniform, isotropic and linear elastic material which
obeys Hooke’s laws, then the principal total stress and strain increments in the radial (r),
tangential (θ) and vertical (z) directions are related as follows:
%& &
where = '%'% ; $ =
'%
Considering equilibrium of the element in Figure 2.20(b) and 2.20(c), the equation for radial
equilibrium is as follows:
* + * = * (2.8)
:; <:; :=
+ = (2.9)
<
Assuming plane strain conditions (i.e. ∆εz = 0) and combining the stress-strain equations (2.3,
2.4 and 2.5) and equation of equilibrium (2.10):
C:; C:= C:; '
=
= ) Δ + Δ + Δ" + 2$Δ , = )2$Δ − 2$Δ , (2.11)
'
) Δ + Δ + 2$Δ , = )2$Δ − 2$Δ ,
(2.12)
'
) + 2$Δ + Δ , = )2$Δ − 2$Δ ,
(2.13)
'
D + 2$ E− F+ Δ E− FG = D2$ E− F − 2$ E− FG
(2.14)
H '
− + 2$ H − + 2$ + + 2$ H = 0 (2.16)
H '
H
+ − H = 0 (2.17)
Powrie (2004) states that equation 2.17 is a differential equation in terms of a single variable
y (y = radial movement). Its general solution, which may be verified by substitution, is:
J
I = + (2.18)
Figure 2.23 - Relationship between allowable vertical stress on stone column and undrained shear strength
(Thorburn & MacVicar, 1968)
Assuming columns are in a critical state, the ultimate vertical stress (σv') is related to the
limiting radial stress by:
'QRSTU
O ′ = E'QRSTUF ′ (2.20)
If the vertical forces exceed the shear resistance along the sides of the column and the
ultimate bearing capacity at the base, then column punching will occur.
Figure 2.24 - Settlement diagram for stone columns in uniform soft clay (Greenwood, 1970)
= V 1V + 1 − 1V (2.22)
The stress concentration ratio is estimated based on field experience. Aboshi et al. (1979)
present stress concentration ratios measured at several construction sites which range from
Literature Review 41
1.6–11.5. However, Barksdale & Bachus (1983) suggest that stress concentration ratios
typically range from 2.5–5.0. The proposed design method is based on an infinite grid of
columns and the settlement of an untreated S and treated S' soil deposit is:
Y = 0Z [\ (2.25)
Y′ = 0Z X [\ (2.26)
Priebe (1995)
Priebe (1995) developed a semi-empirical design method to estimate the settlement of an
infinite grid of end-bearing stone columns. The columns are assumed to be in an active state
and bulge uniformly along their length. The surrounding soil is idealised as an isotropic
elastic material, with an increased coefficient of lateral earth pressure (K0 = 1) to account for
the effects of column installation. The relationship between the applied vertical stress and the
expansion of the column is determined using CCET, the fundamental equations of which are
outlined in 2.5.2. The columns are assumed to be incompressible and a basic settlement
improvement factor (n0), which is defined as the ratio of settlements for an untreated footing
to a treated footing, is developed:
1⁄2 + X_ , ⁄
.M = 1 + ^ − 1b
aJ, X_ , ⁄
(2.28)
It is clear from Figure 2.25 that the settlement improvement factors are significantly
influenced by both column spacing (A/AC) and the strength of the stone backfill (φc).
Literature Review 42
Figure 2.25 - Design curves for basic settlement improvement factor n0 with Poisson’s ratio υ = ⅓ (Priebe,
1995)
Priebe (1995) modifies the basic settlement improvement factor (n0) to account for column
compressibility and the effect of overburden stress. The assumption of column
incompressibility yields infinite settlement improvement factors for an area ratio A/AC = 1.
This is clearly unrealistic as settlement improvement factors should at best be equal to the
ratio of the compression moduli of the column and soil. The effect of column compressibility
is accounted for by increasing the area ratio A/AC by an amount ∆A/AC, which is dependent
on the ratio of the compression moduli of the column and soil. The modification for column
compressibility results in lower settlement improvement factors (n1).
The basic settlement improvement factor also neglects the effect of overburden stress. The
pressure difference at the column-soil interface is assumed to be constant with depth and does
not account for the unit weight of the column and the surrounding soil. Overburden stress
increases with depth and reduces the pressure difference at the column-soil interface.
Therefore, consideration of the overburden stress reduces column bulging and yields higher
settlement improvement factors. A depth factor (fd) is introduced to account for the effect of
overburden stress. The depth factor is defined as the ratio of the original pressure difference
to the ‘new’ pressure difference (which incorporates overburden stress). The settlement of
stone columns is directly related to the depth factor and modified settlement improvement
factors are calculated as n2 = fd×n1.
Priebe (1995) develops design charts relating the settlement of a group of columns beneath
pad footings (s) to the settlement of an infinite grid of columns (s∞). The design curves
account for the stress distribution beneath pad footings and assume a reduced bearing
capacity for the outer columns. It appears from Figure 2.26 that the settlement ratio (s/s∞) is
Literature Review 43
independent of footing area (A). However Priebe states that the footing area must be
calculated from the area ratio A/AC, which compensates for footing size as larger footing
areas have higher A/AC ratios and hence yield lower settlement improvement factors. Priebe
also states that this compensation provides acceptable results for area ratios A/AC up to 10.
The curves in Figure 2.26 indicate that the settlement ratio (s/s∞) reduces rapidly with depth,
especially for smaller groups of columns. This results from a decay of vertical stress with
depth beneath pad footings and hence the influence of the depth factor (fd) is reduced for pad
footings. Therefore, Priebe suggests dividing the subsoil into layers and calculating the
settlement of each layer individually to avoid over-estimating settlements of pad footings.
The settlement is calculated using the following formula:
f
∆ = )⁄h L ?L − ⁄h i ?i ,
gV .
(2.29)
where dL and dU are the lower and upper bound depths of the layer considered.
substratum are insensitive to a smooth or rough boundary condition being assumed at the
substratum. Therefore, it is possible to assume a smooth substratum and hence to find an
exact analytical solution. This solution is based on CCET, the fundamental equations of
which are outlined in 2.5.2.
As a first approximation (solution A), the column is assumed to be laterally restrained. The
authors report that this assumption yields reasonable results, however a stress discontinuity
(∆σr) occurs at the column-soil interface. As the column is much stiffer than the surrounding
soil, it will attract more load and hence will develop a higher radial stress than the
surrounding soil. This stress discontinuity displays itself as column bulging in reality and it is
necessary to account for this. Therefore a second solution (B) is developed for zero vertical
movement of the raft and a laterally expanding column, which imparts a radial stress equal
and opposite to ∆σr at the column-soil interface (Figure 2.27).
Figure 2.27 - Boundary conditions for solutions A & B proposed by Balaam & Booker (1981)
The final solution is found by super-imposing solutions A & B (Table 2.2). The relationship
between strain and the average applied stress qA is determined by integrating the vertical
stresses across the soil surface.
Table 2.2 - Final solution for stresses and strains in column and soil (Balaam & Booker, 1981)
Region 1 Region 2
Stone column Clay
"
zH {| H H }
I y Dy G
|H zH
zH ~ |H
) ' − 2 ' + $' y, D + |H zH E + $ + $ H FG
zH ~ |H
) ' − 2 ' + $' y, D + |H zH E + $ − $ H FG
~zH
" ) ' + 2$' − 2 ' y, D + 2$ + 2 |H zH G
Literature Review 45
& &
Lame’s parameters: = '%'% ; $ =
'%
] = /) ' + 2$' ⁄ + + 2$ 1 − ⁄ − 2 ⁄ ' − y, (2.34)
satisfy the Mohr-Coulomb yield criterion and have a non-associated flow rule. Appropriate
material parameters and geometrical dimensions were selected to comprehensively test the
validity of the assumptions. A good agreement exists between both methods and, therefore,
the solution proposed by Balaam & Booker (1985) is an efficient and accurate way to
calculate the reduction in settlement due to stone columns.
The effectiveness of stone columns in increasing the rate of consolidation is shown Figure
2.28, where the proportion of excess pore pressure dissipation due to the stiffness of the stone
column is shown. It can also be seen that vertical stress transfers from the soil to the stone
columns as consolidation proceeds. The proposed solution was compared with a numerical
study by Balaam & Booker (1981) and a reasonable agreement was observed. However, the
proposed solution tends to under-estimate the rate of consolidation initially, while a reverse in
this trend is observed for rates of consolidation greater than 40 %. The may be explained as
the proposed solution does not account for lateral expansion of columns. Initially, the pore
water pressure in the soil takes the entire load which forces the column to compress. This
mechanism acts as a relief for the soil which increases the rate of consolidation. However, as
consolidation occurs, stress concentration increases in columns which cause columns to
expand into surrounding soil. This increases the excess pore water pressure in the soil and
reduces the rate of consolidation.
Figure 2.28 - Dissipation of excess pore water pressure (Han & Ye, 2001)
at any stage in the loading history. A comparison with Han & Ye (2001) of the development
of stress concentration ratio with time is shown in Figure 2.29.
Figure 2.29 - Stress concentration factor. Influence of radial deformation and plastic strains (Castro & Sagaseta,
2009)
Empirical and analytical design methods to determine the ultimate bearing capacity of
columns were presented in the previous section. These design methods indicate that the
ultimate bearing capacity of stone columns is highly dependent upon the passive resistance of
the surrounding soil, especially in the upper sections of columns (i.e. near the ground
surface).
A series of design methods, ranging from empirical to analytical, are presented in the
previous section to determine the magnitude of settlement for stone columns. Priebe (1995)
develops a semi-empirical design method which contains some simplifying assumptions.
Literature Review 49
Columns are assumed to be in an active state and uniform bulging occurs along the length of
columns. Modification factors are introduced to account for column compressibility and
overburden stresses. Balaam & Booker (1981) propose a more rigorous theoretical solution
than Priebe (1995), but model the stone columns and surrounding soil as linear elastic
materials and, therefore, fail to account for column yielding. This can significantly over-
estimate the effectiveness of stone columns at reducing the settlement of foundation. Balaam
& Booker (1985) and Pulko & Majes (2005) extend this solution to account for column
yielding through an interaction analysis and analytical design method, respectively. The
behaviour of stone columns is idealised as elastic-rigid plastic and column yielding is
captured using dilation theory developed by Rowe (1962).
Stone columns are an effective method to increase the rate of consolidation as their high
permeability allows them to act as vertical drains and also as stress concentrations, which
develop in stone columns, reduce the vertical stress on the surrounding soil. Han & Ye (2001)
develop a simplified solution to compute the rate of consolidation for stone column
reinforced foundations. This solution is extended by Castro & Sagaseta (2009) who model the
stone column as an elasto-plastic dilatant material and account for lateral expansion. This
design method can be used to determine the depth and time of yielding, which allows for an
accurate determination of stresses and strains occurring in columns at any stage in the loading
history.
The deformational behaviour of stone columns is well understood and it is established that
single columns or large groups of columns exhibit either a punching or bulging mode of
deformation. However, the behaviour of small groups of stone columns is more complex and
Muir Wood et al. (2000) demonstrate that small groups of columns can fail by punching,
bulging, bending and shearing. This may be attributed to the loss of lateral confinement for
Literature Review 50
peripheral columns and also the reduction in vertical stress with depth beneath small loaded
areas.
Stone columns can be adopted to treat soft clayey soils, which are typically characterised by
poor strength, stiffness and drainage properties. Consequently, foundations on these soils
undergo large displacements at relatively low loads, and settlement is usually the governing
criterion in foundation design. While the deformational behaviour of small groups of columns
is well understood, there exists a dearth of information regarding the settlement performance
of small groups of columns. This is highlighted by McCabe et al. (2009), who collated a
settlement database from field records, and observed that only three out of 20 case studies
relate to small groups of stone columns. Analytical design methods contain many simplifying
assumptions, such as the unit cell concept which is only applicable to large groups of stone
columns. Elshazly et al. (2008) conduct an axisymmetric FEA to examine the settlement
performance of finite groups of columns; however, the majority of numerical studies are
based on large groups of columns supporting wide area loadings. The lack of information
regarding the settlement performance of small groups of stone columns was identified by
Black (2006), who conducted a series of high quality laboratory tests. However, it is difficult
to extrapolate the findings due to scale effects and also as some of the area ratios considered
are at the high end of typical values used in practice.
It was shown in the laboratory tests that column interaction plays an important role in the
deformational behaviour of small groups of stone columns. Therefore, it is necessary to
model columns at discrete locations, rather than approximating their presence by cylindrical
rings, as is the case with axisymmetric analyses. Furthermore, the decay of vertical stress
with depth will have an important role on the settlement performance of stone columns, and
must be captured when analysing the settlement performance of small groups of stone
columns. It is proposed in this thesis to conduct a three-dimensional FEA using PLAXIS 3D
Foundation to examine the influence of various design parameters, such as area ratio, column
length and the number of columns, upon the settlement performance and deformational
behaviour of small groups of stone columns.
51
Chapter 3
Background and preliminary checks for
Finite Element Analysis
3.1 Introduction
The Finite Element Method (FEM) is a powerful tool used for analysing complex engineering
problems. The FEM is an approximate technique which uses the principle of virtual work to
estimate the distribution of stresses and strains throughout a continuum. PLAXIS 3D
Foundation is a three-dimensional FE program which is specifically tailored for geotechnical
applications. This program is ideal for capturing the complex behaviour of small groups of
stone columns and was adopted for the subsequent FEA. The behaviour of the soil and stone
is simulated with advanced constitutive models, which are described in the following chapter.
In addition, as the FEM is an approximate technique it is necessary to carry out a number of
preliminary checks, such as mesh sensitivity and distance to the boundary, to ensure accurate
numerical analyses.
Mohr-Coulomb model
This is a first order model which idealises soil as an elastic-perfectly plastic material. The
behaviour of soil before failure is approximated by Hooke’s law of elasticity. The failure of
soil is based upon the Mohr-Coulomb failure criterion which is defined by two parameters,
angle of internal friction (φ) and cohesion (c). This failure criterion is an extension of
Coulomb’s friction theory and its yield surfaces in principal stress space are shown in Figure
3.1. The model does not generate irreversible strains below the yield surfaces (Mar, 2002).
However, irreversible plastic strains resulting from shearing are captured using a non-
associated flow rule, which is defined by an angle of dilation (ψ).
Figure 3.1 - Mohr-Coulomb yield surface in principal stress space (c = 0 kPa) (Brinkgreve & Broere, 2006)
A constant Young’s modulus defines the behaviour of soil before failure and it is important to
choose an appropriate value which reflects the stress path and stress level experienced by the
soil. A drawback to the Mohr-Coulomb model is that it fails to accurately capture the
stiffness response of soils due to the simplistic assumption of linear elasticity before failure
and is only used as a first approximation of soil behaviour (Brinkgreve & Broere, 2006).
While the Mohr-Coulomb model captures the failure behaviour of soil quite well in drained
conditions, it over-estimates the undrained shear strength, as continuous dilation is assumed
once yielding occurs (Potts & Zdravkovic, 2001).
dilatancy and introduces a yield cap (Schanz et al., 1999). It accounts for both shear and
volumetric hardening, thus capturing irreversible strains caused by deviatoric and
compression loadings, respectively.
The basis for the Hardening Soil model is that the relationship between the deviatoric stress
(q) and the vertical strain (ε1) for a primary triaxial loading may be approximated by a
hyperbola. Therefore, the yield curve for a standard drained triaxial test, which is shown in
Figure 3.2, can be defined by:
6z 6
' = 6 < 6
2M 6z − 6
(3.1)
The secant Young’s modulus is used instead of the initial modulus (Ei) as this is easier to
define from triaxial tests. The Hardening Soil model also captures the stress path dependency
of soil stiffness using an unload-reload modulus (Eur). The ultimate deviatoric stress (qf) is
derived from the Mohr-Coulomb failure criterion and related to the asymptotic value of shear
strength (qa) by a failure ratio (Rf):
6-. 6z
6 = f + . 8 =
3 − -. 6 (3.2)
Figure 3.2 - Hyperbolic stress-strain relation in primary loading for a standard drained triaxial test (Schanz et al.,
1999)
Background and preliminary checks for FEA 54
Soils tend to exhibit a stress level dependency which results in an increasing stiffness with
increasing confining pressure. The stress dependency of soil stiffness is captured by the
Hardening Soil model using a power law:
=
(3.3)
where Eref is the reference stiffness corresponding to the reference confining pressure (σref).
The confining pressure is taken as the minor or major principal stress for triaxial or
oedometer tests, respectively. The stress dependency is determined by the parameter, m,
which ranges from 0.5 for Norwegian sands and silts (Janbu, 1963) to 1.0 for soft clays
(Brinkgreve & Broere, 2006).
Once the ultimate deviatoric stress is reached, the Mohr-Coulomb failure criterion is satisfied
and perfectly plastic yielding occurs. In addition to capturing plastic shear strains, which
result from decreasing stiffness, the Hardening Soil model also accounts for the volumetric
strains due to dilatancy. Plastic volumetric strains are determined from the stress-dilatancy
theory developed by Rowe (1962). A cap yield surface is introduced to account for
volumetric strains resulting from isotropic loadings. The yield surfaces for the Hardening Soil
model in principal stress space are shown in Figure 3.3.
Figure 3.3 - Hardening Soil yield surface in principal stress space (c = 0 kPa) (Brinkgreve & Broere, 2006)
The Hardening Soil model does not include anisotropic strength and stiffness or time
dependent behaviour of soil. Also, it does not account for the strain dependent nature of soil
stiffness and therefore its application to dynamic loadings is limited.
Background and preliminary checks for FEA 55
(i) (ii)
Figure 3.4 - Distribution of (i) nodes and (ii) stress points within 15 node wedge elements
The influence of the number of elements upon the accuracy of the FEM was investigated by
conducting a mesh sensitivity analysis. In the subsequent parametric studies, both the number
of columns and column spacing are varied beneath pad footings which results in various
footing sizes. Mesh sensitivity analyses are conducted for six different footing sizes where
the accuracy of medium and fine meshes are compared against very fine (vf) meshes.
The vertical displacement (uy) and mean effective stress (p') were measured at three points A,
B and C which are 0, 1 and 2 m below the centre of footings, respectively (Figure 3.5). When
conducting mesh sensitivity analyses it is prudent to examine stresses in the zone of interest,
as the distribution of stress within elements is derived from lower order equations than the
displacement. Therefore, stress converges slower than displacement with increasing mesh
density and, consequently, the distribution of stress within an element will not be as accurate
as the displacement. Vertical displacements were also examined for the mesh sensitivity
analysis as the settlement performance of stone columns is the primary focus of this thesis.
The soil profile adopted for the subsequent FEA and the mesh sensitivity analysis is that of
the well characterised Bothkennar test site, which is described in detail in Section 4.4.
Background and preliminary checks for FEA 56
The accuracy of medium and fine meshes is determined by comparing the normalised error
for vertical displacement (uy) and mean effective stress (p') against very fine meshes:
,
Normalised error for vertical displacement, uy = ,
× 100
Normalised error for mean effective stress, p' =
× 100
10.5 m
14.5 m
Figure 3.6 - Continuous displacement contours along inter-element boundaries beneath a 3 m square footing
In the subsequent FEA, the settlement performance and deformational behaviour of various
configurations of columns is examined at 50 kPa, which is a typical working load for the
Background and preliminary checks for FEA 57
Bothkennar test site. However, most of the mesh sensitivity analyses were conducted at 25
kPa in order to save computational time.
A flowchart outlining the steps taken to ensure accurate numerical analyses is shown in
Figure 3.7.
The number of elements used for each mesh and the normalised error for the vertical
displacement and mean effective stress is shown in Table 3.1. It is clear that the values for
medium and fine meshes are converging towards the very fine mesh and also that the
normalised error between the fine and very fine meshes for the vertical displacement and
mean effective stress is quite low (maximum error is no greater than 2%).
Table 3.1 - Vertical displacement and mean effective stress at points A, B and C
Ftg width, B No. of elements No. of elements Vert displ, uy (mm) Normalised Error (%) Mean eff stress, p' (kPa) Normalised Error (%)
Mesh
(m) in top work plane in 3D mesh A B C A B C A B C A B C
Medium 302 7852 11.6 10.8 9.2 0.9 2.4 8.1 21.5 24.6 25.5 3.1 0.3 0.9
Fine 302 12684 11.7 10.5 8.5 0.5 0.4 0.3 20.8 24.5 25.3 0.4 0.1 0.0
2
Very fine 700 29400 11.7 10.6 8.5 20.9 24.5 25.3
Medium 356 11036 48.3 45.4 38.4 1.5 1.5 0.9 26.5 29.7 29.2 2.2 0.5 0.1
3 Fine 356 16376 48.7 45.7 38.6 0.8 0.8 0.5 25.8 29.7 29.2 0.5 0.3 0.1
Very fine 678 31188 49.1 46.1 38.7 25.9 29.8 29.3
Medium 536 16616 59.4 55.7 43.5 1.1 1.1 0.9 25.7 30.4 24.4 0.4 0.4 0.0
4 Fine 536 23048 59.9 56.2 43.8 0.3 0.3 0.3 25.8 30.5 24.4 0.1 0.1 0.0
Very fine 822 35346 60.1 56.3 43.9 25.9 30.5 24.4
Medium 480 20640 85.4 82.1 71.3 1.8 1.8 1.6 25.9 30.7 32.6 0.8 0.6 0.5
4.5 Fine 1124 48332 86.8 83.5 72.4 0.2 0.2 0.2 26.1 30.9 32.8 0.1 0.1 0.1
Very fine 1370 58910 87.0 83.6 72.5 26.1 30.9 32.8
Medium 392 10976 122.8 118.8 105.6 2.7 2.6 2.5 273.2 32.3 36.2 25.0 1.0 2.0
6 Fine 738 31734 125.5 121.3 107.8 0.6 0.6 0.5 370.9 32.6 35.4 1.7 0.2 0.2
Very fine 1476 63468 126.2 122.0 108.4 364.5 32.6 35.5
Medium 738 20664 170.2 165.1 149.2 1.3 1.2 1.2 30.7 34.5 37.6 6.0 0.4 0.4
8 Fine 738 31734 171.5 166.2 150.2 0.5 0.5 0.5 28.8 34.6 37.7 0.2 0.2 0.2
Very fine 1532 58216 172.4 167.1 151.0 28.9 34.6 37.8
Table 3.2 - Influence of the distance to boundary upon the settlement of a 3 m square footing
Distance to boundary Vertical displacement, uy (mm) Normalised Error (%) Mean effective stress, p' (kPa) Normalised Error (%)
(B = footing width) Point A Point B Point C Point A Point B Point C Point A Point B Point C Point A Point B Point C
2B 61.9 58.4 50.2 26.2 26.8 29.3 26.8 30.6 29.9 0.5 2.8 2.0
4B 53.2 50.0 42.1 8.4 8.5 8.5 26.2 30.1 29.5 1.7 1.0 0.6
6B 49.4 46.3 38.9 0.6 0.6 0.3 25.9 29.8 29.3 2.9 0.1 0.1
8B 49.1 46.1 38.7 0.0 0.0 0.2 25.9 29.8 29.3 2.9 0.1 0.2
10B 48.9 45.9 38.7 0.3 0.3 0.2 26.6 29.8 29.3 0.1 0.1 0.1
12B 49.1 46.1 38.8 0.0 0.0 0.0 26.7 29.8 29.3 0.0 0.0 0.0
Background and preliminary checks for Finite Element Analysis 59
positioning the boundary closer than 6B to the centre of the footing induces more settlement
in the pad footing. The distance to the boundary from the centre of the footing is
conservatively chosen at 8B for all the subsequent FEA.
Method (i) is a closer simulation of reality than method (ii), however it is far more time
consuming. In method (i), the soil is specified to behave in an undrained manner during the
application of footing loads. Initial undrained settlements are computed and long term
settlements are then determined by conducting a consolidation analysis. The soil can be
defined using either total or effective strength parameters for this method. However, it is not
possible to capture the increase in soil shear strength with consolidation or the stress
dependency of soil stiffness when using total strength parameters and, consequently, effective
strength parameters are used in all the subsequent FEA. The initial undrained response of the
soil is simulated by adding a very large bulk modulus to the pore water stiffness matrix,
which ensures that all of the stresses generated from the footing loading will be taken by the
pore water. The consolidation analysis is then conducted by removing the extra bulk
modulus. The stresses in the pore water pressure are then transferred to the soil matrix and
the resulting settlements for primary consolidation are computed. The behaviour of the soil
for method (ii) is again defined using effective strength and stiffness parameters. However, it
is not possible to separate the settlements from the initial undrained response and primary
consolidation for the second method.
The consistency of methods (i) and (ii) is compared in Figure 3.8, which shows the variation
of settlement improvement factors with column length for three arrangements of stone
columns beneath a 3 m square footing. For the purposes of comparing the consistency of both
methods it suffices to examine their relative performance; a detailed analysis of the results is
provided in chapter 5. It appears that method (ii) predicts slightly higher settlement
Background and preliminary checks for Finite Element Analysis 60
improvement factors than method (i) for footings supported by 9 columns, while the opposite
is observed for the footing supported by 5 columns. However, both methods predict a similar
variation of settlement improvement factor with column length for all configurations of
columns (maximum normalised error is no greater than 7%). To save computational time,
method (ii) is therefore adopted for all the subsequent FEA, unless otherwise stated.
4 Method (ii)
10
12
14
Figure 3.8 - Influence of analysis type upon the settlement improvement factors for various groups of columns
The element behaviour is modelled as elastic-plastic, with the Coulomb criterion adopted to
distinguish between elastic and plastic behaviour. The loss of strength at the interface is
modelled with a strength reduction factor (Rinter), which relates the interface strength to the
soil strength through friction angle (φ) and cohesion (c):
Figure 3.9 - Interface elements for a group of 5 columns beneath a 3 m square pad footing
In reality elements have zero thickness, but are assigned a virtual thickness in order to
determine element stiffness. Gap and slip displacements are calculated from the oedometric
(Eoed,i) and shear (Gi) moduli, respectively. The moduli are related by the expression:
1 − R
,R = 2$R
1 − 2R
(3.6)
where υi = 0.45
Guetif et al. (2007) adopt rigid interface elements (i.e. Rinter = 1) on the basis that stone
columns are tightly interlocked with the surrounding soil and a perfect bond exists along the
column-soil interface. However Gäb et al. (2008), Elshazly et al. (2008a), Domingues et al.
(2007a) and many other authors model a perfect bond along the column-soil interface by
omitting interface elements.
The influence of interface elements also appears to become more pronounced with an
increasing number of columns. It will be shown in chapter 5 that closely-spaced columns
punch into the underlying soil, whereas columns at higher A/AC tend to bend and bulge.
Therefore, as A/AC decreases, a larger proportion of the applied load is transferred along the
side of columns and along the interface elements, which induces more displacement.
It is clear from Figure 3.10 that a similar variation of settlement improvement factors with
column length was observed for columns modelled with and without interface elements. As
columns are tightly interlocked with the surrounding soil it was deemed appropriate to model
the column-soil interface in the subsequent FEA by omitting interface elements, in keeping
with the prior work of several other authors.
4
No interface
10
12
14
Figure 3.10 - Influence of interface elements upon settlement improvement factors for various groups of
columns
The change in the stress state of the soil due to column installation may be predicted by an
undrained CCE formulation developed by Gibson & Anderson (1961). As the poker expands
a hole from an initial diameter of zero to the final diameter of a constructed stone column, the
lateral strain is effectively infinite. Consequently, Egan et al. (2008) state that it is necessary
to assume cavity expansion pressure reaches its limit value (plim) in order to accurately
simulate column installation effects.
The relationship between cavity pressure and lateral expansion is derived from CCET and is
shown in Figure 3.11, where the cavity pressure (p) is normalised by the limit pressure and
the current borehole radius (a) is normalised by the initial radius (a0). It can be seen that
considerable lateral expansion is required to reach the limit pressure.
1.0
Normalised Cavity Pressure (p/plim )
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
a = 1.15a0 a = 1.4a0
0.1
0.0
1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0
Normalised Expansion (a/a0 )
Figure 3.11 - Variation of cavity pressure with radius (undrained CCE) (McCabe et al., 2008)
McCabe et al. (2008) simulated the lateral expansion of a 600 mm diameter column in the
Bothkennar soil profile using PLAXIS 3D Foundation. The Mohr Coulomb and Hardening
Soil models were adopted to simulate the behaviour of the stone column and soil,
respectively. Various degrees of lateral expansion were applied (a/a0 = 1.03, 1.06, 1.1, 1.33,
1.67) and the excess pore pressure and radial stress in the soil were measured at mid-depth
along a 5 m long column.
The variation of excess pore pressure and radial stress with normalised radial distance from
the centre of the stone column is shown in Figures 3.12(a) and 3.12(b), respectively. The
excess pore pressure and radial stress is also compared with theoretical curves, developed by
Randolph et al. (1979b) and Gibson & Anderson (1961), respectively. Two curves are
Background and preliminary checks for Finite Element Analysis 64
presented for radial stress in Figure 3.12(b); the first curve corresponds to limit pressure in
the expanding cavity (i.e. field conditions) and the second curve corresponds to the actual
lateral expansion applied in the numerical analysis.
200 200
100 100
50 50
0 0
-50 -50
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Normalised radial distance (r/a) Normalised radial distance (r/a)
(a-i) (b-i)
200
100 100
50 50
0 0
-50 -50
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
(a-ii) (b-ii)
200 200
Change in radial total stress (kPa)
PLAXIS 3-D
150 PLAXIS 3-D 150 CCE
Randolph et al (1979) CCE (plim)
Excess pwp (kPa)
100 100
50 50
0 0
-50 -50
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Normalised radial distance (r/a) Normalised radial distance (r/a)
(a-iii) (b-iii)
Figure 3.12 - Variation of (a) Excess pore pressure (pwp) and (b) total radial stress with normalised radial
distance for lateral expansions (a/a0) of (i) 1.03, (ii) 1.1 and (iii) 1.33 (McCabe et al., 2008)
It can be seen in Figure 3.12 that PLAXIS 3D Foundation data tends towards the plim curve as
lateral expansion increases. It can also be seen that PLAXIS 3D Foundation captures the
variation of radial stress with radial distance from the stone column quite well. However, a
large scatter for excess pore pressure and radial stress is evident, especially for a < r < 2a.
The level of scatter increases with lateral expansion and would not be acceptable at the level
required to reach plim. Furthermore, the lateral extent of the scatter precludes this technique
Background and preliminary checks for Finite Element Analysis 65
for simulating the installation effects of groups of stone columns. The results suggest that
PLAXIS 3D Foundation is currently incapable of simulating the column installation process
by applying large strain cavity expansions. An approximate technique to increase the stress
state of the soil by increasing the coefficient of horizontal stress (K0) is recommended by the
authors and this is investigated in chapter 5.
66
Chapter 4
Development and validation of soft soil
profile
4.1 Introduction
The FEM is a powerful investigative tool which has many advantages over laboratory
experiments and large scale field tests. However, the FEM is still an approximate technique
which idealises real-life situations into a set of continuum components and adopts
constitutive models to simulate soil behaviour. As a result, it is necessary to validate the
output of the FEM to ensure that the real-life situation is accurately modelled.
The main aim of this thesis is to examine the behaviour of small groups of stone columns in
soft clay. A model of the soft soil profile at the well-characterised Bothkennar test site is
developed (using the Hardening Soil model in PLAXIS) in this chapter. The Bothkennar test
site consists of a soft uniform clay, which is overlain by a stiff crust and is representative of
many sites where the applicability of stone columns is of growing interest. The adopted soil
profile and material parameters are validated by back-analysing a load test on an unreinforced
rigid pad footing at the Bothkennar test site, described by Jardine et al. (1995).
theory, a settlement improvement database (collated by McCabe et al., 2009) and a selection
of analytical design methods.
While the majority of the parameters for the Hardening Soil model can be readily converted
from the Duncan & Chang (1970) material model, the latter model does not account for stress
path dependency or formulate a cap yield surface. Therefore, it is necessary to define unload-
reload (Eur) and oedometric (Eoed) moduli for the soil layers. However, the settlement of the
field load tests are insensitive to Eur and Eoed, as the field load tests are monotonically loaded
and do not develop large isotropic stresses (Elshazly et al., 2008a). Therefore, practicable
values of Eur = 5E50 and Eoed = E50 are selected for the Hardening Soil model. A Poisson’s
Development and validation of soft soil profile 69
ratio υ = 0.2 is adopted for both the estuarine and marine deposits, which yields reasonable
lateral deformations and does not under-estimate final settlements.
Elshazly et al. (2008b) state that it is not appropriate to use high permeability values for the
stone column and cohesionless layers as the column is infiltrated with silt and clay particles
resulting from column installation and also as a large proportion of fine particles exist in the
cohesionless layers. Accordingly, a vertical permeability kz = 10-5 m/s and 10-6 m/s is adopted
for the stone columns and cohesionless soil layers, respectively. The vertical permeability for
cohesive layers is within the range 10-8 < kz < 10-7 m/s, as suggested by Lambe & Whitman
(1979) for low liquid limit silts and clays. Elshazly et al. (2008b) adopt a ratio of horizontal
to vertical permeability kx/kz = 2 and state that this ratio does not affect the settlement of the
field load tests, as stone columns act as vertical drains and therefore reduce the influence of
horizontal permeability.
A summary of the material parameters for the Santa Barbara wastewater treatment plant
quoted by Elshazly et al. (2008b) are shown in Table 4.1. It should be noted that the failure
ratio (Rf) for the cohesive marine deposit is different to the original value of Rf = 0.84 quoted
by Mitchell & Huber (1985) and is possibly an error.
Table 4.1 - Material parameters for Santa Barbara wastewater treatment plant (Elshazly et al., 2008b)
Parameter Estuarine Estuarine Marine Marine Gravel
cohesive cohesionless cohesive cohesionless
Dry unit weight, γd (kN/m3) 15 15 17 17 18.6
Saturated unit weight, γsat (kN/m3) 19 19 20 20 21.6
Cohesion, c' (kPa) 0 0 0 0 0
Angle of internal friction, φ' (°) 34 38 34 37 41
Poisson’s ratio, υ 0.2 0.2 0.2 0.2 0.2
Secant Young’s modulus, E50ref (kPa) 8500 17000 8700 12600 29200
Power for stress dependency, m 0.65 0.69 0.65 0.90 0.59
Reference pressure, pref (kPa) 100 100 100 100 100
Failure ratio, Rf 0.87 0.69 0.67* 0.67 0.86
Coefficient of vert. permeability, kz (m/s) 1 × 10-8 1 × 10-6 1 × 10-8 1 × 10-6 1 × 10-5
Coefficient of horiz. permeability, kx (m/s) 2 × 10-8 2 × 10-6 2 × 10-8 2 × 10-6 2 × 10-5
* Mitchell & Huber (1985) adopt Rf = 0.84
was applied in 45 kN increments once the settlement rate dropped to 0.25 mm/hr and was
increased to a maximum of 350–400 kN, with the final load maintained for 6 hrs after the
settlement rate slowed to 0.25 mm/hr. The criterion for increasing the vertical load at 0.25
mm/hr was arbitrarily chosen to reduce the length of the load tests and additional settlement
would have occurred if a longer time increment was adopted (Mitchell & Huber, 1985).
Therefore, the recorded load-settlement curves are not specific to fully drained conditions.
Figure 4.1 - Soil profile for Santa Barbara waste water treatment plant (columns spaced at 1.75 m x 1.75 m)
Development and validation of soft soil profile 71
Figure 4.2 - Conversion of surrounding stone columns into a cylindrical ring for axisymmetric FEA
The behaviour of the soft clay is modelled using methods (i) and (ii), as outlined in Section
3.4.1. The loading process is accounted for in method (i) by applying the load in 45 kN
increments and allowing the soil to consolidate until the settlement rate dropped below 0.25
mm/hr. The extra settlement due to the pause periods can be seen as vertical drops in the
load-settlement curves shown in Figure 4.4. In contrast all the layers are assumed to behave
in a drained manner for method (ii) and do not develop excess pore pressures.
Development and validation of soft soil profile 72
Stone columns:
- d = 1.07 m
- L = 11.1 m
CL
19.2 m
LEGEND
Fill
Concrete footing
Stone column
Estuarine cohesionless
Estuarine cohesive
Marine cohesionless
Marine cohesive
4.3.1 Comparison of PLAXIS 3D Foundation with field records and axisymmetric FEA
The observed load-settlement curves from the 28 field load tests and numerical simulations
are shown in Figure 4.4. A large degree of variation is observed in the field load tests, which
Mitchell & Huber (1985) attribute to the non-uniformity of soil properties in the estuarine
deposits and the variable depth to older marine deposits.
It can be seen in Figure 4.4 that all the numerical analyses over-predict the load-settlement
response from the field tests and this becomes more pronounced as column spacing
decreases. The distribution of fine particles throughout the soil profile is quite varied and
Mitchell & Huber (1985) state that soil layers defined as cohesionless may have contained
enough fine particles to prevent free drainage during the short duration of the load tests. This
would explain the stiffer response of the observed field measurements.
A comparison of the two axisymmetric analyses reveals that Mitchell & Huber (1985) tend to
predict less settlement than Elshazly et al. (2008b), especially for closely spaced columns.
This may be explained as Mitchell & Huber (1985) position the cylindrical rings of stone
columns at radii equal to the orthogonal distance between columns, while Elshazly et al.
(2008b) position the cylindrical rings of stone columns at radii equal to average distance
between orthogonal and diagonal columns (Figure 4.2). In addition, both analyses are based
on different material models, as Mitchell & Huber (1985) adopt the Duncan & Chang (1970)
Development and validation of soft soil profile 73
2 2
2
4 4
4
Settlement (mm)
Settlement (mm)
Settlement (mm)
6 6
6
8 8
8
10 10
Field data Field data Field data
Mitchell & Huber (1985) Mitchell & Huber (1985) Mitchell & Huber (1985)
10
12 Elshazly et al. (2008) 12 Elshazly et al. (2008) Elshazly et al. (2008)
PLAXIS 3D: Method (i) PLAXIS 3D: Method (i) PLAXIS 3D: Method (i)
PLAXIS 3D: Method (ii) PLAXIS 3D: Method (ii) PLAXIS 3D: Method (ii)
14 14 12
(i) (ii) (iii)
Figure 4.4 - Comparison of field load test data with numerical simulations for columns spaced on (i) 1.2 m × 1.5 m, (ii) 1.75 m × 1.75 m and (iii) 2.1 m × 2.1 m
74
hyperbolic model and Elshazly et al. (2008b) adopt the more advanced Hardening-Soil
model. It can be seen that Mitchell & Huber (1985) tend to predict less settlement as the
applied load increases. This may be attributed to the limitations of the hyperbolic model,
which is based on the theory of elasticity.
The presence of the surrounding columns and the material model are two important factors to
consider when modelling the load-deformation behaviour of stone columns. PLAXIS 3D
Foundation appears to capture the load-settlement response of the columns closet to the field
data (Figure 4.4). This is most noticeable for closely spaced columns and reflects the accurate
modelling of the surrounding columns at discrete locations, rather than approximating their
presence by cylindrical rings. This highlights the advantage of PLAXIS 3D Foundation over
the axisymmetric analyses.
It can also be seen in Figure 4.4 that methods (i) and (ii) predict very similar load-
displacement curves. This indicates that the excess pore pressures generated in the cohesive
layers from the undrained loading in method (i) are almost fully dissipated by the time the
next load increment is applied. This may be due to the short drainage length to the
cohesionless layers and to the presence of the stone columns which reduce the radial drainage
length. The close agreement between both methods gives confidence to the use of method (ii)
for the subsequent FEA.
Figure 4.6 - Geotechnical profile at Bothkennar test site (Nash et al., 1992a)
Hight et al. (1992) suggest that post depositional processes such as erosion, changes in
groundwater levels and bonding occurred at the Bothkennar test site. According to the
Development and validation of soft soil profile 76
geological history a maximum drop of 15 kPa in effective overburden pressure could have
occurred at the Bothkennar test site. This implies that over-consolidation ratios should be
high in the upper layers of the Carse clay and reduce with depth. The over-consolidation ratio
determined on the basis of this stress history is compared with data from incremental load
tests by Nash et al. (1992a) in Figure 4.7, where yield stress ratio is equivalent to over-
consolidation ratio. While the profile agrees favourably in the upper layers, it under predicts
the over-consolidation ratio below 4 m where a constant value of 1.55 better represents the
actual stress state. It is suggested that the soil below 4 m may be influenced by ageing, which
accounts for the observed over-consolidation ratio being larger than that due to the stress
history. The profile adopted for all subsequent FEA is shown in Figure 4.7 where the over-
consolidation ratio based on a 15 kPa drop in effective overburden pressure defines the upper
2.5 m of the soil profile and a constant over-consolidation ratio of 1.5 defines the soil below
this level.
Adopted OCR
profile for FEA
(a) (b)
Figure 4.7 - (a) Yield stress and (b) yield stress ratio from incremental load consolidation tests (Nash et al.,
1992a)
The profiles of the lateral total stress measured in situ and the coefficient of lateral earth
pressure (K0) with depth are presented by Nash et al. (1992a) in Figures 4.8(a) and 4.8(b),
respectively. The coefficient of lateral earth pressure was determined from a self-boring
pressuremeter, spade cells and dilatometer tests. The coefficient of lateral earth pressure is
high in the upper layers, which Hight et al. (1992) state is consistent with recent groundwater
fluctuations, and decreases to 0.6–0.9 in the lower Carse layers. The distribution of horizontal
stress with depth adopted for the FEA is shown in Figure 4.8.
Development and validation of soft soil profile 77
Adopted K0
profile for FEA
(a) (b)
Figure 4.8 - Variation of (a) lateral total stress (b) K0 with depth (Nash et al., 1992a)
Depth (m)
8 8
10 CC / (1+e0) CC 10
12 12
14 14
16 16
18 18
20 20
(a) (b)
Figure 4.9 - Variation of (a) compression index Cc and (b) initial voids ratio e0 with depth (Nash et al., 1992b)
The behaviour of the samples prior to yielding was not routinely measured by Nash et al.
(1992b). However, Allman & Atkinson (1992) also investigated the stiffness behaviour of
reconstituted Carse clay and found the slopes of the normal compression (λ) and swelling (κ)
lines to be 0.181 and 0.025, respectively. The ratio between these indices (λ/κ = 7.2) was
used to determine the coefficients of swelling (CS) from the one-dimensional data presented
by Nash et al. (1992b), i.e. λ/κ = CS/CC = 7.2. The one-dimensional stiffness parameters are
converted into three-dimensional parameters for the Hardening Soil model using the
following expressions (Brinkgreve & Broere, 2006):
It is expected that the influence of the stiff crust will, to a certain extent, mask the behaviour
of stone columns in the upper sections, i.e. prevent bulging. In addition, the stiff crust will
absorb a higher proportion of the applied load and thus minimise the influence of elevated
stress levels beneath the edge of rigid footings. However, the stiff crust is an important
feature of soft soil stratigraphy; for example, the soil profile at Kinnegar, Belfast consists of
soft Belfast ‘sleech’ overlain by relatively stiff layers of fill and silty sand (McCabe, 2002).
Therefore, the presence of the stiff crust will yield a more realistic analysis of the settlement
performance and deformational behaviour of small groups of stone columns.
A summary of the parameters developed for the Bothkennar test site is outlined in Table 4.2.
Table 4.2 - Summary of adopted material parameters for the Bothkennar test site
Soil Parameter Crust Upper Carse Lower Carse
clay clay
Depth (m) 0.0 – 1.5 1.5 – 2.5 2.5 – 14.5
Bulk unit weight, γ (kN/m3) 18.0 16.5 16.5
Over-consolidation ratio - - 1.5
Pre-overburden stress (kPa) 15 15 -
Coefficient of lateral earth pressure, K0 1.5 1.0 0.75
Effective cohesion, c' (kPa) 3 1 1
Angle of internal friction, φ (°) 34 34 34
Initial voids ratio, e0 1.0 1.2 2.0
Compression index, CC 0.07 0.25 1.12
Swelling index, CS 0.01 0.03 0.16
Reference pressure, pref (kPa) 13 20 30
Vert. coefficient of permeability, kvert (m/day) 6.9 × 10-5 6.9 × 10-5 6.9 × 10-5
Horz. coefficient of permeability, khorz (m/day) 1.0 × 10-4 1.0 × 10-4 1.0 × 10-4
Development and validation of soft soil profile 80
80
Test A
40
0
0 20 40 60 80 100 120 140
Time (hrs)
Figure 4.10 - Variation of load with time (Jardine et al., 1995)
The footings were instrumented with pneumatic piezometers, spade cells, inclinometers and
magnetic extensometers. The surface settlement was measured with a precise level, giving a
nominal resolution of 0.1 mm. In total ten targets were set into the concrete pads and 16
targets were augured into the surrounding soil.
The load test on Pad B was not simulated with PLAXIS 3D Foundation as this load test
contained an unload-reload loop and the recorded footing displacement includes both primary
and secondary settlement (Jardine et al., 1995). The aim of this thesis is to examine the long
term primary settlement response of small groups of stone columns supporting pad footings.
In addition, the Hardening Soil model does not account for secondary settlement and is
therefore unable to accurately simulate the load test on Pad B.
140
120
Jardine et al. (1995)
Applied pressure (kPa)
Plaxis 3D Foundation
100
80
60
40
20
0
0 20 40 60 80 100 120 140 160 180 200
Vertical displacement (mm)
Figure 4.11 - Load-displacement behaviour for pad footings (Jardine et al., 1995)
Figure 4.12 - Characteristic stiffness-strain behaviour of soil with typical strain ranges for laboratory tests and
structures (after Atkinson & Sallfors, 1991)
settlements, and not just the settlement improvement factors, are in the correct order of
magnitude.
The settlement performance of an infinite grid of stone columns is then analysed over a
typical range of A/AC values presented in the settlement database by McCabe et al. (2009).
The unit cell concept was adopted to simplify the analysis of an infinite grid of columns to
one column and its surrounding zone of influence (see Section 2.5.1). Settlement
improvement factors determined from PLAXIS 3D Foundation are compared with a
settlement database of field tests (collated by McCabe et al., 2009) and a selection of current
analytical design methods.
The settlement performance of stone columns at the Bothkennar test site was determined at
50 kPa. The ultimate bearing capacity recorded by Jardine et al. (1995) for an unreinforced
footing at the Bothkennar test site was 138 kPa (Figure 4.11). Applying typical factors of
safety (2.5–3.0) yields allowable bearing pressures in the range of 46–55 kPa. Therefore, 50
kPa is deemed a typical working load for the Bothkennar test site.
The angle of internal friction adopted for stone columns φ = 45° is chosen on the basis of
findings by McCabe et al. (2009), who conducted a review of field tests on stone columns
and found that φ = 40° is conservative for columns formed using the bottom feed system. A
nominal value for cohesion c = 1 kPa is adopted to ensure numerical stability. The angle of
dilatancy is determined from the empirical relationship ψ = φ – 30°, developed by Bolton
(1986). A positive angle of dilatancy ensures that columns dilate when subject to high stress
ratios.
Development and validation of soft soil profile 84
Barksdale & Bachus (1983) report that Young's moduli (E) for stone columns, which are
back-calculated from measured settlements and recommended by other authors, typically
range from 30–58 MPa. However, these values are specific to columns formed using the top
feed system and may under-estimate Young’s moduli for stone columns following the
assertion by McCabe et al. (2009) that columns formed using the dry bottom feed system
yield a better settlement performance. Consequently, a higher Young’s modulus E50 = 70
MPa is adopted for the stone columns. A similar value was also adopted by Gäb et al. (2008),
who conducted a numerical investigation of stone columns supporting an embankment in
loose-medium compacted sand and weak clayey silt. The unload-reload Young’s modulus Eur
= 210 MPa is also similar to Gäb et al. (2008) and its selection is based on the relationship
(Eur = 3E50) proposed by Brinkgreve & Broere (2006).
V ′,M + ∆′
_ = \ ^ 5 b
1 + 9M ′,M
(4.3)
where σ'y, 0 and σ'y, max = in situ and maximum vertical effective stress, respectively
The Bothkennar soil profile is divided into seven layers and the settlement is calculated on
the basis of the vertical effective stress at the centre of each layer (Table 4.4). The settlement
for each layer and the variation of settlement with depth is shown in Figures 4.13(i) and
4.13(ii), respectively. It can be seen in Figure 4.13(i) that a slight variation exists between
PLAXIS 3D Foundation and one-dimensional compression theory for the settlement of each
Development and validation of soft soil profile 85
layer with depth, which may be attributed to a difference in the definition of the cap yield
surface. However, it is clear from Figure 4.13(ii) that PLAXIS 3D Foundation predicts the
settlement response of an untreated infinitely wide footing at the Bothkennar test site quite
well (normalised error for the total settlement is less than 3%).
2 2
4 4
Depth (m)
Depth (m)
6 6
8 8
10 10
Plaxis 3D Foundation Plaxis 3D Foundation
12 12
1D compression theory 1D compression theory
14 14
(i) (ii)
Figure 4.13 - Settlement profile of Bothkennar test site for wide area loading, calculated using one-dimensional
compression theory and PLAXIS 3D Foundation
However, McCabe et al. (2009) suggest that settlement improvement factors for some of the
case studies may have been under
under-estimated due to the following reasons:
(i) The time at which settlement was measured may not be consistent for all tests and
untreated settlement may be under
under-estimated,
estimated, as primary consolidation was not complete
when the measurements were taken.
(ii) Not all of the field data are specific to end
end-bearing
bearing columns and, therefore, the measured
settlement includes the settlement of the soil deposit underneath the base of columns.
The high n values predicted by PLAXIS 3D Foundation may also be related to the high
compressibility of the lower Carse clay layer relative to the stone columns, i.e. high modular
ratio (Ecol/Esoil). Field experience suggests that Ecol/Esoil typically ranges from 10–50
10 (Castro
& Sagaseta, 2009). However, columns at the Bothkennar test site are formed in a soft clay
with an average Ecol/Esoil = 100 (see Figure 6.7) for the lower Carse clay layer. This is
significantly higher than the upper limit suggested by Castro & Sagaseta (2009), but is not
unrealistic as a finite element study by Barksdale & Bachus (1983) indica
indicates
tes Ecol/Esoil up to
100 for columns formed in soft cohesive soil.
11
Venmans (1998) Cooper and Rose (1999)
10 Munfakh et al (1983) Greenwood (1991)
Settlement Improvement Factor, n
A finite element study by Domingues et al. (2007a) shows that high Ecol/Esoil leads to an
increase in n values. The authors observed that increasing Ecol/Esoil from 10 to 100 and from
40 to 100 yielded an increase of n values in the range 80
80–200%. These increases
creases are applied
to Priebe’s (1995) basic design (n0) and are shown in Figure 4.14.. It can be seen that PLAXIS
3D Foundation agrees quite well with field data when adjusted for a high Ecol/Esoil. This gives
Development and validation of soft soil profile 87
confidence to the selection of the material parameters adopted for the stone backfill (see
Section 4.6.1) and also the ability of PLAXIS 3D Foundation to predict the settlement
performance of large groups of stone columns in soft clay.
The coefficients of lateral earth pressure and material parameters adopted for the stone
backfill in the analytical design methods are outlined in Table 4.5.
The parameter 'm' for the Bothkennar soil profile can be determined from one-dimensional
stiffness data presented by Nash et al. (1992b). The values presented by Nash et al. (1992b)
are converted into oedometric moduli (Eoed) using equation 4.1. The oedometric moduli are
then normalised by the effective overburden stress (σ'y,0), at the level which the samples were
taken from the soil, and plotted against σ'y,0 in Figure 4.15. In this case the major principal
stress is adopted for the confining pressure as the samples are subject to oedometric tests. It
can be seen for the lower Carse clay that Eoed/σ'y,0 is constant with an increasing σ'y,0, which
suggests that Eoed is directly proportional to σ'y,0. This indicates that m = 1 is appropriate for
the lower Carse clay layer, which is consistent with Brinkgreve & Broere (2006) for soft
soils. Insufficient data is available to determine the parameter 'm' for the crust and the upper
Development and validation of soft soil profile 89
Carse clay and similar values to the lower Carse clay were adopted. A lower value of 'm' is
more appropriate for the granular material and m = 0.3 is used for stone columns, which is
similar to Gäb et al. (2008).
40
60 Lower
Carse Clay
80
100
120
140
Figure 4.15 - Variation of normalised oedometric moduli with vertical effective stress for Bothkennar test site
While the parameter 'm' is now determined for the Bothkennar soil profile and stone columns,
it is still necessary to define the confining pressure for stone columns in the field condition. A
FEA conducted by Balaam & Booker (1985) indicates that columns are in a triaxial state of
stress during loading. Therefore, the minor principal stress (σ3) defines the confining pressure
for the stone column. The surrounding soil also carries a proportion of the applied vertical
load and, consequently, is compressed in vertical direction. As with stone columns, the minor
principal stress is used to define the confining pressure for the surrounding soil.
The variation of minor principal stress with depth for stone columns and the surrounding soil
is presented in Figure 4.16 for an infinite grid of end-bearing columns (A/AC = 3.5). It can be
seen that the average minor stress is slightly higher in the stone column, which may be
explained as this material is stiffer than the surrounding soil and therefore takes more load.
The average minor stress for each layer is tabulated in Table 4.6 and the resulting Young’s
moduli are determined from equation 3.3.
Development and validation of soft soil profile 90
Depth (m)
6
10
12
14
Figure 4.16 - Distribution of minor principal stress σ3 in end-bearing stone column and surrounding soil
(A/Ac = 3.5; L = 13.9 m; Infinite grid of columns)
Table 4.6 - Average minor principal stresses (σ'3) and Young’s moduli (E) for Bothkennar soil profile
Soil Stone backfill
Depth σ'ref Eref (= E50) σ'3 E σ'ref Eref (= E50) σ'3 E
(m) (kPa) (kPa) (kPa) (kPa) (kPa) (kPa) (kPa) (kPa)
0.6 – 1.5 13 1068 28 2304 100 70000 35 51215
1.5 – 2.5 20 506 28 703 100 70000 32 49848
2.5 – 14.5 30 231 53 410 100 70000 57 59120
The design method proposed by Balaam & Booker (1981) significantly over-estimates the
settlement performance of stone columns, which is due to the simplified assumption of linear
elasticity for the stone backfill (Castro & Sagaseta, 2009). This design method is shown as
it forms the basis of Pulko & Majes (2005) and, also, as it demonstrates that linear
elasticity is not a valid assumption for stone columns.
Development and validation of soft soil profile 91
19
Plaxis
PLAXIS3D 3D
Foundation
Foundation
17 Priebe:
Priebe:n0
n0 (1995)
11
7
Priebe: n2
5 K0, lower Carse clay = 1.00
K0, lower Carse clay = 0.75
3
1
2 4 6 8 10 12 14 16
Area ratio, A/Ac
Figure 4.17 - Comparison of settlement improvement factors for an infinite grid of end-bearing columns with
analytical design methods
The design method proposed by Pulko & Majes (2005) is an extension of Balaam & Booker
(1981) which accounts for column yielding and, hence, predicts more realistic n values.
While Pulko & Majes (2005) predicts a similar settlement performance to PLAXIS 3D
Foundation for A/AC > 8.0, it tends to over-predict n values at low A/AC. This may be
attributed to the assumption of rigid-plastic behaviour for stone columns by Pulko & Majes
(2005), in contrast to elasto-plastic behaviour (i.e. Hardening Soil model) adopted in PLAXIS
3D Foundation. The main difference between these models is that the rigid-plastic model
does not account for elastic strains developed when columns are in a plastic state. Therefore,
the rigid-plastic model predicts less settlement, and thus higher n values, than the elasto-
plastic model. The extent of plasticity becomes more pronounced for an infinite grid of
columns at low A/AC, as columns are better confined and can carry a larger proportion of the
applied load. Consequently, the discrepancy between the rigid-plastic model and the elasto-
plastic model is most pronounced at low A/AC.
The basic (n0) and modified (n2) settlement improvement factors determined from Priebe
(1995) are also shown in Figure 4.17. It appears that the basic design curve (n0) under-
estimates the settlement performance of stone columns. Priebe (1995) modifies the basic
design curve to account for column compressibility and the effect of overburden stress, which
leads to a significant increase in n values. The modified n2 values agree favourably with the
numerical analysis, although the settlement performance is slightly over-estimated. This over-
estimation may be attributed to the assumption of K0 = 1 for the lower Carse clay layer,
Development and validation of soft soil profile 92
which is higher than that adopted in Pulko & Majes (2005) and PLAXIS 3D Foundation.
Settlement improvement factors were re-calculated for Priebe (1995) n2 with K0 = 0.75 for
the lower Carse clay layer. It can be seen that the re-calculated n values are very close to
Pulko & Majes (2005) and PLAXIS 3D Foundation.
A soil profile and set of material parameters were developed for the well characterised
Bothkennar test site, which consists of a soft uniform clay overlain by a stuff crust. This
soil profile is representative of many sites where the applicability of stone columns is of
growing interest and, therefore, forms the basis for the subsequent FEA. The adopted
soil profile and choice of material parameters were validated by simulating a field load
test, which was conducted on an unreinforced (i.e. no stone columns) rigid pad footing at
the Bothkennar test site.
Chapter 5
Results of FEA: Settlement performance,
deformational behaviour and stress
concentration ratios
5.1 Background
The main focus of this research is to develop a better understanding of the behaviour of small
groups of stone columns in soft soils. A series of FEA were conducted using PLAXIS 3D
Foundation to examine the influence of key design parameters upon the settlement
performance and deformational behaviour of small groups of stone columns. The key design
parameters and considerations investigated are as follows:
• Column length (L)
• Area ratio (A/AC)
• Column confinement
• Column arrangement
• Column position relative to footing edge
• Column compressibility
• Column strength
• Column installation effects
• The presence of a stiff crust
The design of foundations on soft soils is usually governed by settlement rather than bearing
capacity criteria, due to their high compressibility (Priebe, 1976). Therefore, the settlement
performance of stone columns at working load levels is of the upmost importance. The
degree of confinement provided to columns is dependent upon the number of columns
and the area ratio. The influence of column confinement is investigated by analysing the
settlement performance of different configurations of columns (i.e. different group sizes)
over a range of area ratios. The settlement performance of each configuration of columns
is examined for various column lengths, which allows the relationship between column
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 94
length, area ratio and column confinement to be determined. Finally, the influence of the
key design parameters and practical considerations outlined above is investigated for a
select number of column configurations, which were specifically chosen to cover a wide
range of area ratios.
The modes of deformation for small groups of stone columns are examined with the aid of
new parameters called compression and punching ratios. In previous studies the modes of
deformation have usually been determined at ultimate conditions; however, this study
examines the modes of deformation at typical working loads. This allows a direct link to be
made between the deformational behaviour and settlement performance of stone columns.
The stress concentration ratio is an important parameter which compares the vertical stress in
columns to the vertical stress in the surrounding soil. For this research the vertical stress used
to determine this parameter is measured beneath the base of footings (i.e. at the top of stone
columns). Stress concentration ratios are used in analytical design methods to determine the
settlement performance and the potential of stone columns to stabilise slopes; the variation of
this ratio with the modes of deformation and the key design parameters is investigated.
and infinite grids of columns (see Figure 5.2). Single columns and infinite grids provide a
useful frame of reference for the other configurations of columns as these two cases
correspond to zero and full column confinement, respectively. Columns are positioned on a
square grid and spaced at 1.0, 1.5 and 2.0 m, which corresponds to A/AC of 3.5, 8.0 and 14.1,
respectively. This is a typical range of A/AC for small loaded areas, as can be seen in the
settlement database developed by McCabe et al. (2009); columns at A/AC < 3.5 would not be
practicable or economic to construct and columns at A/AC > 14.1 would typically be used to
support wide area loadings such as embankments, which can tolerate larger settlements. The
edge of the footing is located at a distance of half the column spacing (0.5×s) from the
centreline of the outer row of columns; the influence of the footing overhang is examined in a
subsequent parametric study.
Column arrangement
The column arrangements adopted to support a 3 m square pad footing are shown in Figure
5.3. The arrangements of 4, 5 and 9 columns are carefully chosen to allow the benefit of
individual columns (i.e. centre, edge and corner) to be assessed. A comparison of the 4 and 5
column groups allows the benefit of an extra central column to be assessed, while a
comparison of the 5 and 9 column groups allows the benefit of four edge columns to be
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 96
assessed. A group of 7 columns may also be adopted to support small pad footings in practice
and this configuration was compared with a group of 9 columns.
Column position
The effect of footing overhang and column position was examined for a 3 m square pad
footing by positioning a 2×2 group of columns progressively closer to the footing edge
(Figure 5.4). It is well known from elastic theory that the distribution of vertical stress
beneath rigid footings is non-uniform and high vertical stress develops beneath the edge of
footings. The column spacings investigated were 1.0, 1.5 and 2.0 m.
Figure 5.4 - Column configurations to examine influence of column position relative to the footing edge
Figure 5.5 - Column configurations for parametric study of column compressibility and strength
Figure 5.6 - Various profiles adopted to investigate the effect of the lower Carse clay (Profile 2) and the stiff
crust (Profile 3)
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 98
4 4
6 6
8 8
10 10
A/Ac = 3.5
12 A/Ac = 8.0 12
A/Ac = 14.1
14 14
(i) (ii)
Figure 5.7 - Variation of (i) footing settlement and (ii) settlement improvement factor with column length for an
infinite grid of columns
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 99
It can be seen in Figure 5.8(a-i) that settlement improvement factors for columns at low area
ratios (i.e. A/AC = 3.5) increase continuously with column length. A jump in settlement
improvement factors is observed when moving from floating to end-bearing columns (i.e. L =
13→13.9 m) at low area ratios. The improved settlement performance may be explained as
the proportion of applied load transferred to the base of columns is supported by a rigid
stratum in the case of end-bearing columns, rather than a compressible stratum in the case of
floating columns. The jump in settlement improvement factors suggests that closely-spaced
columns transfer a significant proportion of the applied load to their base. The magnitude of
the jump increases for larger groups of columns, which indicates that the proportion of
applied load transferred to the base of columns increases for larger groups of columns. This
reflects the enhanced lateral confinement associated with larger groups of columns.
The variation of settlement improvement factors with column length for various groups of
columns at A/AC of 8.0 and 14.1 is shown in Figures 5.8(a-ii) and 5.8(a-iii), respectively.
Similar to columns at low area ratios, a continuous increase in settlement improvement
factors is observed with increasing column length for all columns. However, the rate of
increase in settlement improvement factors reduces for column lengths longer than L = 5 m.
This suggests that long widely-spaced columns transfer a smaller proportion of the load to the
base of columns and may be deforming closer to the ground surface.
It is clearly shown in Figure 5.8(a) that the enhanced levels of lateral confinement associated
with larger groups of columns yields higher settlement improvement factors. However, it
appears that the variation of settlement improvement factors with column length is dictated
by area ratio rather than the number of columns, as a continuous increase in settlement
improvement factors is observed with column length at low area ratios (i.e. A/AC = 3.5),
while the increase in settlement improvement factors with column length tends to reduce
beyond L = 5 m at higher area ratios (i.e. A/AC = 8.0 and 14.1).
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 100
The influence of the footing width (B) upon the settlement performance of small groups of
columns is shown in Figure 5.8(b). For the column configurations analysed, it appears that
A/AC governs the variation of settlement improvement factors with column length. The
classical Boussinesq (1885) solution for vertical stress distribution beneath an unreinforced
square footing suggests that the stress increment applied to the footing is no longer perceived
at L/B = 2 and so improving the settlement characteristics of the soil below this level, by
installing stone columns, should not enhance the settlement performance. However, it appears
that this significant depth (L/B = 2) does not influence settlement improvement factors of
columns at low area ratios, as a continuous increase in settlement improvement factors is
observed with increasing column length beyond this level.
Settlement improvement factors for small groups of end-bearing stone columns (n) are
normalised by settlement improvement factors for an infinite grid of end-bearing stone
columns (nuc) and presented in Figure 5.10. The n/nuc ratio allows the influence of column
confinement upon small groups of columns to be assessed, as an infinite grid of columns has
full confinement and is an ideal benchmark. An increase in n/nuc is observed with both an
increasing number of columns and an increasing area ratio. It can be seen in Figure 5.10 that
n/nuc increases with area ratio, from 0.19 to 0.72 for single columns and from 0.45 to 0.82 for
a 4×4 group of columns. The increase in n/nuc for larger groups of columns reflects the
increased level of lateral confinement. However, the increase in n/nuc with increasing area
ratio may be explained as the influence of lateral confinement upon settlement improvement
factors reduces at high area ratios. Consequently, the difference between settlement
improvement factors for small groups (n) and infinite grids (nuc) of stone columns reduces at
high area ratios, which yields an increase in n/nuc.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 101
6 6 6
8 8 8
10 10 10
12 12 12
14 14 14
(a-i) (a-ii) (a-iii)
4 4 4
6 6 6
8 8 8
The relationship between the settlement of small groups of columns and infinite grids (s/suc)
is shown in Figure 5.11. Priebe (1995) also developed a design chart (Figure 2.24) which
relates the settlement of small groups of columns to infinite grids. However, Priebe (1995)
does not directly account for the influence of footing area and relies on A/AC to compensate
for a decrease in s/suc with an increasing footing area. It can be seen in Figure 5.11 that both
the area ratio and the number of columns have a significant influence on the settlement
performance of small groups of columns. For similar area ratios, it can be seen that s/suc
increases for a larger group of columns. This may be attributed to an increase in footing size
for larger groups of columns. It can also be seen that the increase in s/suc for larger groups of
columns is more pronounced at high area ratios (A/AC = 14.1). This highlights the positive
confining effects of closely-spaced columns.
11
1 col
Settlement improvement
2×2
9
3×3
4×4
factor, n
7 Unit Cell
Increasing
confinement
5
1
3 5 7 9 11 13 15
Area ratio, A/Ac
Figure 5.9 - Settlement improvement factors for small groups of end-bearing columns
1.0
Normalised n value, n/nuc
0.8
0.6
0.4
1 col
2×2
0.2
3×3
4×4
0.0
3 5 7 9 11 13 15
Area ratio, A/AC
Figure 5.10 - Normalised settlement improvement factors for small groups of end-bearing columns
0.5
Settlement ratio, s/suc
0.4
0.3
0.2
A/Ac=3.5
0.1 A/Ac=8.0
A/Ac=14.1
0.0
0 5 10 15 20
No. of columns
Figure 5.11 - Normalised settlement for small groups of end-bearing columns
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 103
The arrangement of columns can be tailored to achieve a specific settlement performance, i.e.
long columns at high area ratios or short columns at low area ratios. The volume of stone
required for each configuration of columns is shown in Figure 5.12(ii). It appears that
installing an extra centre column in a 2×2 group is only beneficial for long columns. It also
appears that short columns at low area ratios are the most efficient configuration to achieve
high settlement improvement factors. However, it must be noted that the volume of stone
does not entirely reflect the overall cost of construction. Other factors such as the time taken
to form stone columns must also be taken in account, i.e. a large number of short columns
will take longer to construct than a small number of long columns due to the time taken to
remove the poker. The longer length of time taken to construct a larger number of columns
translates itself as an additional cost and must also be taken into consideration.
4 7 cols 7 cols
9 cols 9 cols
15
6
20
8
25
10
30
4 s = 2.0 m
10
12
14
Figure 5.13 - Influence of column position upon settlement improvement factors
reduces for larger groups of columns. This may be explained as the average applied load
carried per column reduces for larger groups. This suggests that larger groups of columns are
in more of an elastic state and hence are more sensitive to changes in elastic stiffness moduli.
Settlement improvement factor, n
1.0 2.0 3.0 4.0
0
4 cols - s = 2.0 m
5 cols - s = 1.0 m
2 9 cols - s = 1.0 m
E50, col = 30 MPa
E50, col = 50 MPa
10
12
14
Figure 5.14 - Influence of column compressibility upon settlement improvement factors
Table 5.1 - Influence of column compressibility upon settlement improvement factors for end-bearing columns
Settlement improvement factors % reduction
No. of cols Column compressibility, E50 (MPa) Column compressibility, E50 (MPa)
30 50 70 30 50 70
4 1.57 1.64 1.69 -7 -3 -
5 1.76 1.89 1.94 -9 -2 -
9 2.60 3.23 3.68 -29 -12 -
Crust
4
φcol = 50°
10
12
14
Figure 5.16 - Influence of column strength upon settlement improvement factors
Table 5.2 - Influence of column strength upon end-bearing settlement improvement factors
Settlement improvement factors % reduction / increase
No. of cols Angle of internal friction, φ (°) Angle of internal friction, φ (°)
40 45 50 40 45 50
4 1.43 1.69 2.09 -15 - 24
5 1.57 1.94 2.49 -19 - 29
9 2.46 3.68 4.32 -33 - 17
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 107
4 K0 = 1.00
K0 = 1.25
6
10
12
14
Figure 5.17 - Influence of coefficient of lateral earth pressure upon settlement improvement factors
Table 5.3 - Influence of coefficient of lateral earth pressure upon end-bearing settlement improvement factors
Settlement improvement factors % reduction / increase
No. of cols Coefficient of lateral earth pressure, K0 Coefficient of lateral earth pressure, K0
0.75 1.00 1.25 0.75 1.00 1.25
4 1.69 1.88 2.10 - 11 24
5 1.94 2.24 2.49 - 16 29
9 3.68 4.00 3.91 - 9 6
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 108
Settlement performance
Profile 2 consists of a stiff crust which is underlain by upper Carse clay. This soil profile is
stiffer at depth than Profile 1, as it does not contain the lower Carse clay layer. This is
reflected in Figure 5.18(i) as the untreated (L = 0 m) and treated settlements are lower than
Profile 1. Profile 3 does not contain the stiff crust and consists entirely of upper Carse clay.
While it is softer than Profile 1 near the surface, it is stiffer at depth and this is again reflected
in Figure 5.18(i) as the treated settlements are lower than Profile 1. It is interesting to note
that the untreated settlements (L = 0 m) for Profiles 1 and 3 are similar, which suggests that
the increase in settlement due to the absence of the stiff crust is offset by the increased
stiffness at depth, as the lower Carse clay is replaced by upper Carse clay.
It can also be seen in Figure 5.18(ii) that settlement improvement factors increase
significantly with column length up to L = 3 m in Profile 3. However, the rate of increase in
settlement improvement factors with column length slows down thereafter. This is most
pronounced for smaller groups of columns (i.e. higher A/AC) and suggests that a critical
length may exist for columns formed in Profile 3. This profile is similar to that used in
laboratory studies described in the literature review (see Sections 2.2 and 2.3) and may
explain why authors have postulated critical lengths in the past. However, it can be seen that
critical lengths are less well defined for more realistic soil profiles. Therefore, this highlights
the benefits of the FEM which accurately captures the influence of the stiff crust.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 109
2 2
Column length, L (m)
6 6
8 8
10 10
4 cols - 2.00 m 4 cols - 2.00 m
Profile 1 Profile 1
12 12
Profile 2 Profile 2
Profile 3 Profile 3
14 14
(a-i) (a-ii)
2 2
Column length, L (m)
4 4
6 6
8 8
10 10
5 cols - 1.00 m 5 cols - 1.00 m
Profile 1 Profile 1
12 12
Profile 2 Profile 2
Profile 3 Profile 3
14 14
(b-i) (b-ii)
2 2
Column length, L (m)
4 4
6 6
8 8
10 10
9 cols - 1.00 m 9 cols - 1.00 m
Profile 1 Profile 1
12 12
Profile 2 Profile 2
Profile 3 Profile 3
14 14
(c-i) (c-ii)
Figure 5.18 - Influence of the lower Carse clay and the stiff crust upon (i) settlement and (ii) settlement
improvement factors for (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 110
Punching ratio
Depending on the adopted configuration of columns, a proportion of the applied load may be
transferred to the base of the columns which results in columns punching vertically into the
underlying soil. In order to compare the punching of columns beneath different footing sizes,
the displacement at the base of columns is normalised by the displacement of the footing:
7_
− 7Q
R
[7.ℎ-.5 18- =
7
RS¡ (5.1)
The displacement of the soil surrounding the base of columns (usoil) is determined by
averaging soil displacements within a zone surrounding each column. The size of the zone
adopted to determine usoil will influence the magnitude of the punching ratio, as soil
displacement decreases with distance from columns. Two zones, shown in Figure 5.19, were
compared in order to find a consistent method of calculating the punching ratio. The first
zone (i) is a 1 m square zone surrounding the column and the second zone (ii) is also a square
zone, whose width is equal to the column spacing. The size of the second zone will always be
greater than or equal to the first zone as column spacing ranges from 1.0–2.0 m in this study.
Figure 5.19 - Displacement at surface of footing, base of columns and soil surrounding base of columns
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 111
The variation of punching ratio - determined by methods (i) and (ii) - with length for an
infinite grid of columns is shown in Figure 5.20. The variation of punching ratio with length
for A/AC = 3.5 is not shown as this corresponds to a column spacing of 1.0 m; the zones for
methods (i) and (ii) coincide at this spacing and no difference will be observed. It can be seen
in Figure 5.20 that methods (i) and (ii) predict a similar variation of punching ratio with
length. It is also clear that the punching ratio is always higher for method (ii). This may be
explained as the displacement of the surrounding soil is averaged over a larger area for the
method (ii). This leads to a lower average displacement of the soil and, consequently, higher
punching ratios, as soil displacement decreases with distance from the column. The punching
ratio is determined using the method (i) for the subsequent analyses, as the zone in which the
soil displacement is averaged is the same size for all column spacings and thus allows for a
consistent analysis.
2
Column length, L (m)
Unit cell
10
(i) - A/Ac = 8.0
(i) - A/Ac = 14.1
12 (ii) - A/Ac = 8.0
(ii) - A/Ac = 14.1
14
Figure 5.20 - Variation of punching ratios with length for an infinite grid of columns, with usoil measured in (i)
a 1 m square zone and (ii) a square zone, whose width is equal to the column spacing
Compression ratio
The compression ratio is an important indicator of the mode of deformation and is defined as
the proportion of surface settlement transferred to the column base. Low compression ratios
suggest that columns are transferring the load to depth, rather than deforming along their
length. One limitation of the compression ratio is that is does not differentiate between axial
and radial deformation. However, axial deformation due to column compressibility results in
a steady increase in compression ratios, while radial deformation due to column bulging can
be identified by a rapid increase in compression ratios.
7
RS¡ − 7_
0f9-. 18- =
7
RS¡ (5.2)
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 112
6 6
8 8
Unit cell
10 10
A/Ac = 3.5
12 A/Ac = 8.0 12
A/Ac = 14.1
14 14
(i) (ii)
Figure 5.21 - Variation of (i) punching and (ii) compression ratios with column length for an infinite grid of
columns
It can be seen in Figure 5.21(i) that the punching ratio increases with column length up to a
maximum at L = 3 m. This increase may be attributed to a reduction in the stiffness of soil
layers in the upper 3 m of the Bothkennar test site, i.e. the crust is the stiffest layer, followed
by the upper Carse clay and then the lower Carse clay (see CC values in Table 4.2).
Therefore, as column length increases, the columns bases are founded in more compressible
soil strata which results in higher punching ratios (Figure 5.22). It is interesting to note that
the maximum punching ratio occurs at L = 3 m, despite the fact that the base of 2 m long
columns coincides with the weakest part of the soil profile, i.e. the top of the lower Carse clay
layer. It seems intuitive that the maximum punching ratio would occur when the column
bases are formed in the weakest part of the soil profile. However, it can be seen in Figure
5.22(ii) that the base of 2 m long columns extends just beyond the bottom of the upper Carse
clay layer, which reduces the displacement of the soil at the top of the lower Carse clay layer,
thus, reducing the punching ratio.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 113
It is also evident in Figure 5.21(i) that the punching ratio increases with increasing area ratio
and this is most pronounced at L = 3 m. The deformation of an infinite grid of 3 m long
columns at various area ratios is shown in Figure 5.23. It can be seen in Figure 5.23(i) that
columns at low area ratios act with the surrounding soil and punch as a single entity into the
underlying soil. This mode of deformation is referred to as 'block' failure and is a form of
column punching. It occurs for large groups of closely-spaced columns and is characterised
by low punching and low compression ratios. High punching ratios usually denote punching
failure; however, punching ratios are low for 'block' failure as columns deform with the
surrounding soil. A similar mode of deformation is observed by Black (2006) for a group of
three columns at low area ratios (A/AC = 2.5 and 3.6). It can be seen in Figures 5.23(ii) and
5.23(iii) that column interaction decreases at higher area ratios, as columns tend to act more
individually and punch into the underlying soil.
While it can be seen in Figure 5.21 that columns at all area ratios exhibit similar variations of
punching and compression ratios with column length up to L = 3 m, a marked difference is
observed with area ratio thereafter:
• Columns at low area ratios continue to exhibit relatively low and constant punching
ratios for column lengths L > 3 m. Compression ratios also remain relatively low, which
indicates that columns do not deform significantly along their length. This suggests that
'block' failure remains the dominant mode of deformation with increasing length for an
infinite grid of columns at low area ratios. A slight increase in the punching ratio is
observed for column lengths L > 11 m.
• Columns at high area ratios exhibit a decrease in punching ratio and increase in
compression ratio for column lengths L > 3 m. The increase in compression ratios
indicates that columns are deforming along their length and may be bulging closer to the
ground surface. Columns which are bulging cannot transfer the applied load to the base
of columns which consequently results in low punching ratios.
It can be seen in Figures 5.24(i) and 5.25(i) that an increase in area ratio leads to a decrease in
punching ratio for 1 m long single columns and 2×2 groups of columns, respectively. This
may be explained as the soil surrounding at the base of 1 m long columns is significantly
influenced by the footing overhead. An increase in area ratio for groups of columns
corresponds to an increase in footing width, and wider footings induce more displacement in
the soil at z = 1 m (z = depth beneath footing). As a consequence, this reduces the differential
displacement between the base of columns and the surrounding soil, which results in lower
punching ratios.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 115
6 6
8 8
1 column
10 10
A/Ac = 3.5
12 A/Ac = 8.0 12
A/Ac = 14.1
14 14
(i) (ii)
Figure 5.24 - Variation of (i) punching and (ii) compression ratios with column length for single columns
4 4
A/Ac = 14.1
6 6
8 8
2×2 cols
10 10
A/Ac = 3.5
12 A/Ac = 8.0 12
A/Ac = 14.1
14 14
(i) (ii)
Figure 5.25 - Variation of (i) punching and (ii) compression ratios with column length for 2×2 groups of
columns
However, it can also be seen that an increase in area ratio leads to an increase in punching
ratios for L = 2, 3 and 4 m. It is established that punching is the dominant mode of
deformation for short stone columns, which implies that columns develop shear stress and
end-bearing pressure along the side and at the base of columns, respectively. An increase in
area ratio corresponds to an increased footing width for groups of columns. This implies that
the total load taken per column also increases with area ratio. Columns of similar length
which are 'punching' have the capacity to carry similar shear force along their sides.
Therefore, an increase in the total load taken per column increases the portion of the load
transferred to the base of columns, which consequently increases the punching ratio.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 116
A reduction in the punching ratio is coupled with an increase in compression ratios for single
columns and 2×2 groups which are longer than L = 3 m. This suggests that columns are
deforming along their length and that bulging is the dominant mode of deformation. It is
interesting to note that the reduction in punching ratios and increase in compression ratios
becomes more pronounced with increasing area ratio. This is consistent with the findings
from an infinite grid of columns and indicates that column bulging becomes more
pronounced with increasing area ratio.
The deformed shape of a 3×3 group of columns at L = 3 m and 6m is shown in Figure 5.28. It
can be seen in Figure 5.28(i) that closely-spaced columns are acting as a 'block' with the
surrounding soil. This is consistent with relatively low punching and compression ratios,
evident in Figures 5.26 and 5.27 for 3×3 and 4×4 groups, respectively. It appears that 'block'
failure remains the dominant mode of deformation with increasing column length for closely-
spaced columns. In contrast, it can be seen that columns at higher area ratios which are longer
than L = 3 m deform significantly with increasing column length. A reduction in punching
ratios and increase in compression ratios suggests that columns are deforming along their
length and do not transfer the applied load to the base of columns (see Figure 5.28). It can
also be seen that external columns tend to bend outwards and away from central columns.
This mechanism was also observed by McKelvey et al. (2004) in small scale laboratory tests.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 117
8 8
3×3 cols
10 A/Ac = 3.5 10
Centre
12 Edge 12
Corner
14 14
(a-i) (a-ii)
4 4 Edge
Corner
6 6
8 8
3×3 cols
10 A/Ac = 8.0 10
Centre
12 Edge 12
Corner
14 14
(b-i) (b-ii)
4 4 Edge
Corner
6 6
8 8
3×3 cols
10 A/Ac = 14.1 10
Centre
12 Edge 12
Corner
14 14
(c-i) (c-ii)
Figure 5.26 - Variation of (i) punching and (ii) compression ratios with column length for a 3×3 group of
columns spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 118
8 8
4×4 cols
10 A/Ac = 3.5 10
Centre
12 Edge 12
Corner
14 14
(a-i) (a-ii)
4 4 Edge
Corner
6 6
8 8
4×4 cols
10 A/Ac = 8.0 10
Centre
12 Edge 12
Corner
14 14
(b-i) (b-ii)
4 4 Edge
Corner
6 6
8 8
4×4 cols
10 A/Ac = 14.1 10
Centre
12 Edge 12
Corner
14 14
(c-i) (c-ii)
Figure 5.27 - Variation of (i) column and (ii) column ratios with column length for a 4×4 group of columns
spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 119
Figure 5.28 - Deformed shapes from PLAXIS 3D Foundation for 3×3 group of columns at area ratios of (i) 3.5, (ii) 8.0 and (iii) 14.1 for column lengths (a) L = 3 m and (b) L
=6m
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 120
The variation of punching ratios for individual columns (i.e. centre, edge and corner columns)
within a 3×3 and 4×4 group can also be seen in Figures 5.26 and 5.27, respectively. For
columns shorter than L = 3 m, it appears that punching ratios are lowest for centre columns,
followed by edge and corner columns. This trend is also evident for all lengths of closely-
spaced columns (i.e. A/AC = 3.5), which suggests that this variation of punching ratios with
column length is associated with a punching mode of deformation. Columns which are
punching impart shear stresses to the surrounding soil which tends to drag the surrounding
soil downwards. Punching ratios are lowest for centre columns as the soil surrounding the
base of the column is dragged downwards on all sides by edge and corner columns. However,
it can be seen for widely-spaced columns (i.e. A/AC = 8.0 and 14.1) which are longer than L
= 3 m that punching ratios for centre columns exceeds punching ratios for corner and edge
columns. It is evident in Figure 5.28 that long widely-spaced edge columns tend to bulge and
bend outwards, which reduces punching at the base of columns. Centre columns are better
confined and can transfer more of the applied load to the base.
A reduction in the punching ratio with an increasing number of columns is also observed in
Figures 5.29(b-i) and 5.29(c-i) for columns at higher area ratios. It is shown in Figure
5.28(iii) that column bending occurs in the outer row of widely-spaced columns. Bending
reduces the ability of columns to transfer the applied load to depth and hence results in lower
punching ratios. For a constant area ratio, footing width increases with an increasing number
of columns. Columns become increasingly slender as footing width increases and,
consequently, are more susceptible to bending. Therefore, punching ratios reduce with an
increasing number of columns at high area ratios.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 121
8 8
A/AC = 3.5
10 1 col 10
2×2
3×3 - average
12 12
4×4 - average
Unit Cell
14 14
(a-i) (a-ii)
3×3 - average
4 4
4×4 - average
Unit Cell
6 6
8 8
A/AC = 8.0
10 1 col 10
2×2
3×3 - average
12 12
4×4 - average
Unit Cell
14 14
(b-i) (b-ii)
4 4 3×3 - average
4×4 - average
Unit Cell
6 6
8 8
A/AC = 14.1
10 1 col 10
2×2
3×3 - average
12 12
4×4 - average
Unit Cell
14 14
(c-i) (c-ii)
Figure 5.29 - Influence of column confinement upon (i) punching and (ii) compression ratios for various
lengths and configurations of columns spaced at area ratios of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 122
It can be seen in Figure 5.29(i) that unit cell conditions yield the lowest punching ratios for
columns at low area ratios. However, this is not the case for high area ratios and it can be
seen that punching ratios for the unit cell are larger than those for groups of columns longer
than 7 m and 5 m in Figures 5.29(b-i) and 5.29(c-i), respectively. The external columns in
groups experience a loss of lateral confinement and bend outwards, which reduces the
punching ratio. An infinite grid of columns cannot bend due to symmetry and punching ratios
are higher, despite the increased level of column confinement.
The relationship between column confinement and compression ratios is shown in Figure
5.29(ii). It can be seen for closely-spaced columns that the average compression ratios of a
group reduces as the number of columns increases. This is consistent with findings from the
settlement performance, where the jump in settlement improvement factors observed for end-
bearing columns with an increasing number of columns. This indicates that larger groups of
columns transfer higher proportions of the applied load to the base, which reduces the
deformation along the length of columns.
While all configurations of columns exhibit a similar variation of punching and compression
ratios with increasing length for short columns, a marked change is observed with increasing
length thereafter. The groups of 4 and 5 columns exhibit a significant reduction in punching
ratio and a simultaneous increase in compression ratio. This indicates that these column
groups are not transferring the applied load to the base of columns and are deforming closer
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 123
8 8
B=3m
10 4 cols 10
5 cols - average
12 7 cols - average 12
9 cols - average
14 14
(i) (ii)
Figure 5.30 - Variation of (i) punching and (ii) compression ratios with column length for various
configurations of stone columns beneath a 3 m square footing
to the ground surface. This is consistent with previous findings for long widely-spaced
columns. In contrast to this behaviour, relatively low compression ratios are observed for the
groups of 7 and 9 columns. This suggests that punching failure remains the dominant mode of
deformation with increasing column length.
8 8
4 cols - B = 3 m
10 s = 1.0 m 10
s = 1.5 m
12 s = 2.0 m 12
s = 1.0 m (B = 2 m)
14 14
(i) (ii)
Figure 5.31 - Variation of (i) punching and (ii) compression ratios with column length for various column
spacings beneath a 3 m square footing
The deformational behaviour of a 2×2 group spaced at 1.0 m and supporting a 2 m and a 3 m
square pad footing is also compared in Figure 5.31. It appears that the mode of deformation is
different for the group of columns beneath the 2 m square footing (i.e. A/AC = 3.5). Increased
punching ratios and lower compression ratios indicate that punching remains the dominant
mode of deformation for all column lengths. Therefore, it is postulated that area ratio is more
critical than column spacing in governing the deformational behaviour of columns.
8 8
4 cols - s = 2.0 m
10 10
E_50 = 30 MPa
12 E_50 = 50 MPa 12
E_50 = 70 MPa
14 14
(a-i) (a-ii)
4 4
E_50 = 70 MPa
6 6
8 8
5 cols - s = 1.0 m
10 10
E_50 = 30 MPa
12 E_50 = 50 MPa 12
E_50 = 70 MPa
14 14
(b-i) (b-ii)
E_50 = 50 MPa
4 4
E_50 = 70 MPa
6 6
8 8
9 cols - s = 1.0 m
10 10
E_50 = 30 MPa
12 E_50 = 50 MPa 12
E_50 = 70 MPa
14 14
(c-i) (c-ii)
Figure 5.32 - Influence of column compressibility upon the variation of (i) punching and (ii) compression ratios
with column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 126
8 8
4 cols - s = 2.0 m
10 10
φ = 40°
12 φ = 45° 12
φ = 50°
14 14
(a-i) (a-ii)
4 4
φ = 50°
6 6
8 8
5 cols - s = 1.0 m
10 10
φ = 40°
12 φ = 45° 12
φ = 50°
14 14
(b-i) (b-ii)
4 4
φ = 50°
6 6
8 8
9 cols - s = 1.0 m
10 10
φ = 40°
12 φ = 45° 12
φ = 50°
14 14
(c-i) (c-ii)
Figure 5.33 - Influence of column strength upon the variation of (i) punching and (ii) compression ratios with
column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 127
can be seen that compression ratios decrease with an increasing angle of friction for the stone
column. This results in a corresponding increase in punching ratios. The changes in punching
and compression ratios are more pronounced for fewer columns and only negligible changes
are observed when increasing the angle of internal friction from 45° to 50° for the group of 9
columns.
Punching ratio
It can be seen in Figure 5.35(i) that punching ratios increase with column length up to a
maximum at L = 3 m for Profile 1. A similar variation is observed for Profile 2, however
punching ratios reach a maximum at L = 2 m and are lower than Profile 1 for column lengths
L > 3 m. This reflects the increased stiffness at depth for Profile 2, as the lower Carse clay in
the standard Bothkennar soil profile (Profile 1) is replaced with the stiffer upper Carse clay.
Profile 3 consists entirely of upper Carse clay and, therefore, allows for the influence of the
stiff crust to be assessed. It can be seen that punching ratios are largest near the surface for
this soil profile and decrease with increasing column length thereafter. Punching ratios for
Profile 3 are higher than Profiles 1 and 2 for short columns, but this trend reverses for
columns longer than L = 3 m.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 128
Compression ratio
The variation of compression ratios with column length is shown in Figure 5.35(ii). It can be
seen that compression ratios are lowest for Profile 1, which indicates that columns formed in
Profiles 2 and 3 are more susceptible to a bulging mode of deformation. This may be
explained by a number of reasons:
(i) The lower Carse clay in Profile 1 is replaced by the stiffer upper Carse clay in Profiles 2
and 3. Therefore, the base of floating columns in Profile 2 and 3 is formed in a more
competent soil layer. The enhanced support at the base of columns reduces punching and
forces columns to compress more along their length.
(ii) The absence of the stiff crust (i.e. Profile 3) leads to the highest compression ratios. The
stiff crust prevents bulging in the upper sections of columns, where overburden stresses
are lowest.
It appears that the modes of deformation for columns formed in a homogenous soil profile, as
used in the laboratory tests described in the literature review (see Sections 2.2 and 2.3), are
very different to columns formed in a more realistic soil profile, as modelled in PLAXIS 3D
Foundation. This is discussed in more detail in Section 7.3.3.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 129
8 8
4 cols - s = 2.0 m
10 10
K0_soil = 0.75
12 K0_soil = 1.00 12
K0_soil = 1.25
14 14
(a-i) (a-ii)
K0_soil = 1.00
4 4
K0_soil = 1.25
6 6
8 8
5 cols - s = 1.0 m
10 10
K0_soil = 0.75
12 K0_soil = 1.00 12
K0_soil = 1.25
14 14
(b-i) (b-ii)
K0_soil = 1.00
4 4
K0_soil = 1.25
6 6
8 8
9 cols - s = 1.0 m
10 10
K0_soil = 0.75
12 K0_soil = 1.00 12
K0_soil = 1.25
14 14
(c-i) (c-ii)
Figure 5.34 - Influence of column installation effects upon the variation of (i) punching and (ii) compression
ratios with column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square
footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 130
2 2
Column length, L (m)
6 6
8 8
2 2
Column length, L (m)
4 4
6 6
8 8
2 2
Column length, L (m)
4 4
6 6
8 8
Figure 5.35 - Influence of stiff crust upon the variation of (i) punching and (ii) compression ratios with column
length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 131
(i) (ii)
Figure 5.36 - Determination of stress concentration ratios for (i) field/laboratory studies and (ii)
numerical/analytical studies
Stress concentration ratios from PLAXIS 3D Foundation are defined by two methods, which
are illustrated in Figure 5.37:
(i) The first method compares the average vertical stress in each column (σ'col) to the
average vertical stress in the entire soil area beneath the footing (σ'soil, average).
(ii) The second method compares the average vertical stress in each column (σ'col) to the
average vertical stress in the soil within a square zone of influence surrounding each
column (σ'soil); the width of the zone of influence is equal to the column spacing.
The first method allows the stress levels in corner, edge and centre columns to be directly
compared with each other, as the average stress over the entire area of soil is used. The
second method is a measure of the work each column is doing relative to the stress in the
surrounding zone of influence and, therefore, indicates the efficiency of each column.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 132
(i) (ii)
Figure 5.37 - Stress concentration ratio determined by (i) averaging stress over entire soil area and (ii)
averaging stress within the zone of influence for each column
Stress concentration ratios determined from both methods will be identical for single
columns, 2×2 groups and infinite grids as due to symmetry, the area of soil beneath footings
corresponds to the zone of influence for each column. However, this is not the case for
groups of columns larger than 2×2 configurations as the stress in the zone of influence varies
for different columns. By way of example, the distribution of vertical effective stress
immediately beneath a 4.5 m square concrete footing, which is (i) unreinforced (L = 0 m) and
(ii) reinforced with a 3×3 group of end-bearing (L = 13.9 m) stone columns, is shown in
Figure 5.38. It can be seen that high vertical stress develops beneath the edges of the rigid
unreinforced footing. The high stress levels beneath the edges of the footing reduce for the
reinforced case as columns absorb a larger proportion of the load and, therefore, reduce the
stress on the surrounding soil.
(i)
(ii)
Figure 5.38 - Distribution of vertical effective stress (σ'vert) immediately beneath a 4.5 m square footing, which
is (i) unreinforced (L = 0 m) and (ii) reinforced with end-bearing stone columns (L = 13.9 m)
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 133
A slight reduction in the stress concentration ratio is observed at high area ratios. Column
confinement reduces at high area ratios, which consequently reduces the stiffness of stone
columns and hence yields lower stress concentration ratios.
2
Column length, L (m)
10 Unit cell
A/Ac = 3.5
12 A/Ac = 8.0
A/Ac = 14.1
14
Figure 5.39 - Variation of stress concentration ratio with column length for an infinite grid of columns
shear stress along the sides and end-bearing pressure at the base of columns. The load-
carrying capacity of columns increases with increasing column length as shear stress can
develop over a larger area. As a result, stress concentration ratios increase in concert.
2 2
Column length, L (m)
6 6
8 8
Stress concentration ratios continue to increase with column length for closely-spaced
columns (i.e. A/AC = 3.5), which is consistent with previous findings that closely-spaced
columns punch into the underlying soil. While large increases in stress concentration ratios
are initially observed with increasing column length, the rate of increase in stress
concentration ratios with column length reduces significantly for column lengths L > 3 m.
This may attributed to large percentage increases in column length for short columns (100%
for L = 1 → 2 m) compared to longer columns (10% for L = 10 → 11 m).
It can also be seen in Figures 5.40(i) and 5.40(ii) that stress concentration ratios for columns
at A/AC = 8.0 and 14.1 reach maximum values at L ≈ 4 m. Stress concentration ratios remain
relatively constant with increasing column length thereafter for columns at A/AC = 8.0, which
indicates that the mode of deformation has changed from punching to bulging. Column
bulging occurs at the point of lowest lateral support, which occurs near the top of the lower
Carse clay (z ≈ 3 m). Increasing the column length for this mode of deformation does not
enhance the load-carrying capacity of columns and, therefore, does not increase stress
concentration ratios. A slight reduction in stress concentration ratios for columns at A/AC =
14.1 is observed with increasing column length in the range 4–7 m and 4–5 m for single
columns and 2×2 groups, respectively. A combination of punching and bulging modes of
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 135
deformation may be occurring over this range of lengths. However, bulging becomes more
prevalent as column length increases and stress concentration ratios decrease slightly towards
constant values.
A reduction in stress concentration ratios is observed for higher area ratios at L = 1 m. This is
consistent with findings from the deformational behaviour of single columns and 2×2 groups,
as shown in Figures 5.24(i) and 5.25(i), respectively, where a reduction in column punching
is observed for higher area ratios at L = 1 m. The footing size increases at higher area ratios,
which stresses the soil deeper and induces more displacement in the soil surrounding the base
of columns at z = 1 m. Therefore, the end-bearing resistance at the base of columns, and
consequently the load-carrying capacity of stone columns, reduces with increasing area ratio.
However, it is interesting to note that stress concentration ratios increase at a faster rate at
higher area ratios for L = 2, 3 and 4 m. The average load carried per column increases with
footing area. Therefore, a larger force is carried along the sides and at the base of columns.
This results in high punching ratios, as shown in Section 5.4.3, and also forces columns to
carry a larger proportion of the applied load which results in higher stress concentration
ratios.
2 2
Column length, L (m)
6 6
8 8
3×3 cols 3×3 cols
A/Ac = 3.5 A/Ac = 3.5
10 10
Centre Centre
12 Edge 12 Edge
Corner Corner
14 14
(a-i) (a-ii)
2 2
Column length, L (m)
4 4
6 6
8 8
3×3 cols 3×3 cols
A/Ac = 8.0 A/Ac = 8.0
10 10
Centre Centre
12 Edge 12 Edge
Corner Corner
14 14
(b-i) (b-ii)
2 2
Column length, L (m)
4 4
6 6
8 8
3×3 cols 3×3 cols
A/Ac = 14.1 A/Ac = 14.1
10 10
Centre Centre
12 Edge 12 Edge
Corner Corner
14 14
(c-i) (c-ii)
Figure 5.41 - Variation of stress concentration ratio with column length for a 3x3 group of columns, spaced at
area ratios of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 137
2 2
Column length, L (m)
6 6
8 8
4×4 cols 4×4 cols
A/Ac = 3.5 A/Ac = 3.5
10 10
Centre Centre
12 Edge 12 Edge
Corner Corner
14 14
(a-i) (a-ii)
2 2
Column length, L (m)
4 4
6 6
8 8
4×4 cols 4×4 cols
10 A/Ac = 8.0 10 A/Ac = 8.0
Centre Centre
12 Edge 12 Edge
Corner Corner
14 14
(b-i) (b-ii)
2 2
Column length, L (m)
4 4
6 6
8 8
4×4 cols 4×4 cols
A/Ac = 14.1 A/Ac = 14.1
10 10
Centre Centre
12 Edge 12 Edge
Corner Corner
14 14
(c-i) (c-ii)
Figure 5.42 - Variation of stress concentration ratio with column length for a 4x4 group of columns, spaced at
area ratios of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 138
concentration ratios between individual columns becomes more pronounced for lower area
ratios. The external row of columns are positioned a distance 0.5×s (s = column spacing)
from the edge of footings. As column spacing decreases for lower area ratios, columns are
positioned closer to the edge of footings and hence closer to the zone of elevated stress levels.
Therefore, the variation in stress concentration ratios is most significant at low area ratios.
Variation of stress concentration ratio with column length and area ratio
The variation of stress concentration ratios with column length for 3×3 and 4×4 groups of
columns at A/AC = 3.5 is shown in Figures 5.41(a-i) and 5.42(a-i), respectively. An increase
in stress concentration ratios with increasing column length is observed for all columns up to
L = 4 m and L = 3 m for 3×3 and 4×4 groups of columns, respectively. A slight reduction in
stress concentration ratios with increasing column length is subsequently observed for corner
columns up to L ≈ 7 m, while stress concentration ratios in edge and centre columns continue
to increase or remain constant. As outlined previously, corner columns carry a greater
proportion of the applied load due to the position of columns relative to the footing edge.
However, corner columns are more susceptible to bulging and bending as column length
increases. Therefore, corner columns cannot sustain the high levels of vertical stress and
transfer the applied load to edge and centre columns. Consequently, stress concentration
ratios continue to increase with column length for edge and centre columns.
The transfer of vertical stress from corner columns to edge and centre columns is also evident
in the variation of stress concentration ratios, measured using the second method (σ'col/ σ'soil).
It can be seen in Figures 5.41(a-ii) and 5.42(a-ii) for 3×3 and 4×4 groups of columns,
respectively, that corner and edge columns initially carry a larger proportion of the applied
load relative to the surrounding soil. However, as column length increases a larger proportion
of the load is carried by centre columns relative to the surrounding soil. This reflects the
higher levels of lateral confinement associated with central columns.
The variation of stress concentration ratios with column length for columns spaced at high
area ratios (i.e. A/AC = 8.0 and 14.1) is shown in Figures 5.41(b-i) and 5.41(c-i) for 3×3
groups of columns and 5.42(b-i) and 5.42(c-i) for 4×4 groups of columns. As with closely-
spaced columns, it can be seen that stress concentration ratios increase with column length up
to L = 3 m. This is again consistent with a punching mode of deformation and is in keeping
with previous findings for short widely-spaced columns. A slight reduction in stress
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 139
concentration ratios with increasing column length is subsequently observed for corner
columns up to L ≈ 5 m, while stress concentration ratios continue to increase for centre
columns at A/AC = 8.0. This indicates that the heavily loaded corner columns transfer the
applied vertical load to centre columns. No change in stress concentration ratio with
increasing column length is observed thereafter, which suggests that bulging is the dominant
mode of deformation. While no transfer of stress is evident for columns at an area ratio A/AC
= 14.1, it can be seen that stress concentration ratios for all columns (i.e. corner, edge and
centre) reduce slightly with increasing column length up to L ≈ 6 m and are constant
thereafter. The absence of stress transfer from corner to centre columns may be attributed to
the relatively uniform stress levels in all columns. The slight reduction in stress concentration
ratios with is also evident from single and 2×2 groups of columns and suggests that both
punching and bulging modes of deformation are occurring. Bulging becomes more prevalent
as column length increases and no change in stress concentration ratios is observed with
increasing column length.
The positive influence of column confinement can be seen in Figures 5.41(b-ii) and 5.41(c-ii)
for 3×3 groups of columns and 5.42(b-ii) and 5.42(c-ii) for 4×4 groups of columns. Stress
concentration ratios - measured using the second method (i.e. the vertical stress in columns
relative to the local stress level) - are highest in the centre columns. This reflects the
increased levels of lateral confinement associated with centre columns. This is in contrast to
edge and corner columns which bulge and bend outwards away from central columns.
8 A/AC = 3.5
1 col
10 2×2
3×3 - average
12 4×4 - average
Unit Cell
14
(a)
2
Column length, L (m)
8 A/AC = 8.0
1 col
10 2×2
3×3 - average
12 4×4 - average
Unit Cell
14
(b)
2
Column length, L (m)
8 A/AC = 14.1
1 col
10 2×2
3×3 - average
12 4×4 - average
Unit Cell
14
(c)
Figure 5.43 - Variation of stress concentration ratio for small groups and infinite grids of columns spaced at
area ratios of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 141
the potential to absorb a larger proportion of the applied load. Consequently, the average
stress concentration ratios for the group of columns increases.
A similar trend is observed for short columns at higher area ratios (i.e. A/AC = 8.0 and 14.1)
in Figures 4.43(b) and 4.43(c). It was established previously that punching is the dominant
mode of deformation for this configuration of columns. An increase in stress concentration
ratios with an increasing number of columns is again due to the position of columns relative
to the edge of the footing. However, once column length exceeds L = 3 m, the mode of
deformation changes to bulging and no difference in stress concentration ratios is observed.
Columns which are bulging fail at the same location (i.e. the weakest part of the soil profile).
Therefore, the ultimate load is the same for all columns and is not dependent upon the
position of columns relative to the edge of footings.
The stress concentration ratios for a unit cell are relatively uniform as columns are located
within an infinite grid and therefore are not subject to the elevated stress levels near the edge
of rigid footings.
columns. The lower stress concentrations ratios for the group of 5 columns are observed due
to the location of the central column in a region of lower stress levels (see Figure 5.44).
Therefore, this column takes less of the applied load and the average stress concentration
ratio for the group is reduced.
B=3m
10
4 cols
5 cols - average
12
7 cols - average
9 cols - average
14
Figure 5.45 - Variation of stress concentration ratio with column length for a group of 4, 5, 7 and 9 columns
beneath a 3 m square footing
The variation of stress concentration ratios with length for the groups of 7 and 9 columns is
quite different to groups of 4 and 5 columns. Stress concentration ratios for groups of 7
columns and 9 columns increase at a slower rate with column length up to L = 3 m than the
smaller groups of columns. Furthermore, stress concentration ratios for the groups of 7 and 9
columns continue to increase with length longer than L = 3 m, which suggests that punching
remains the dominant mode of deformation with increasing column length. However, the rate
of increase in stress concentration ratio reduces with increasing column length and stress
concentration ratios are eventually constant for the group of 7 columns longer than L = 8 m.
The magnitude of stress concentration ratios is higher for the group of 9 columns than the
group of 7 columns. This is due to the position of columns relative to the footing edge. The
outer ring of columns for a 9 column group are positioned closer to the footing edge and
therefore closer to the elevated stress levels which develop in this region.
concentration ratios increase with column length up to L = 3 m, which indicates that all
column configurations are punching into the underlying soil. However, no increase in the
stress concentration ratio is observed with increasing column length thereafter, which
indicates that the mode of deformation has changed to bulging. A slight decrease in stress
concentration ratios with column length is observed in the range 3–6 m for columns spaced at
2.0 m, which indicates that a transition from punching to bulging occurs for this range of
column lengths.
It can be seen in Figure 5.46 that the position of columns relative to the footing edge has a
significant influence upon stress concentration ratios. Stress concentration ratios increase as
columns are positioned closer to the footing edge; an increase of c.50% is observed at L = 3
m. This is due to elevated stress levels which occur beneath the edge of rigid footings, as
shown in Figure 5.47. Positioning columns closer to this zone allows columns to absorb more
load and hence develop higher stress concentration ratios. However, it was shown in Section
5.3.5 that the position of columns relative to the footing edge has a relatively minor influence
upon the settlement performance of stone columns. Therefore, stress concentration ratios
measured at the surface do not define the settlement behaviour of a group of columns. This is
an interesting finding and suggests that the variation of stress concentration ratios with depth,
which is not routinely measured in laboratory or field tests, may play an important role in
determining the settlement performance of stone columns. The variation of stress
concentration ratio with depth is investigated as part of this thesis in Section 6.3.
The influence of footing overhang upon stress concentration ratios is also investigated by
comparing a 2×2 group of columns spaced at 1.0 m beneath a 2 m and a 3 m square footing.
It can be seen in Figure 5.46 that stress concentration ratios for the 3 m square footing
increase with column length to a maximum at L = 3 m and are constant thereafter. In contrast,
a continuous increase in stress concentration ratios with increasing column length is observed
for the 2 m square footing. Furthermore, it can be seen that stress concentration ratios are
higher for the 2×2 group of columns beneath a 2 m footing. This is again due to the position
of the columns close to the footing edge.
It appears that the position of columns and the mode of deformation both play important roles
in the magnitude of stress concentration ratios.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 145
8
4 cols - B = 3 m
10 s = 1.0 m
s = 1.5 m
12 s = 2.0 m
s = 1.0 m (B = 2 m)
14
Figure 5.46 - Variation of stress concentration ratio with column length for various columns positions beneath a
3 m square footing
(i)
(ii)
(iii)
Figure 5.47 - Vertical effective stress beneath a 3 m square footing for columns spaced at (i) 1.0 m, (ii) 1.5 m
and (iii) 2.0 m
stiffness is controlled by the secant Young’s modulus (E50) and it can be seen that an increase
in E50 leads to higher stress concentration ratios.
It can be seen for columns shorter than L = 3 m that the influence of columns stiffness upon
stress concentration ratios becomes more pronounced as column length increases up to L = 3
m. However, the variation in stress concentration ratios with E50 does not change with
increasing column length thereafter for groups of 4 and 5 columns, as shown in Figure
5.48(a) and 5.48(b), respectively. This may be related to the bulging mode of deformation
which becomes more prevalent for groups of 4 and 5 columns which are longer than L = 3 m.
It was shown previously that the load-carrying capacity, and hence the stress concentration
ratio, of columns which are bulging does not increase with column length.
The variation of stress concentration ratios with column length for a group of 9 columns is
shown in Figure 5.48(c). As with the groups of 4 and 5 columns, the influence of column
stiffness becomes more pronounced with increasing column length up to L = 3 m. Stress
concentration ratios continue to increase at a slower rate with increasing column length
thereafter for E50 = 70 MPa. However, no increase in stress concentration ratios is observed
with column length for E50 = 30 MPa, which suggests that the mode of deformation has
changed from punching to bulging.
A linear relationship appears to exist between column stiffness and stress concentration ratios
as shown by stress concentration ratios for end-bearing columns in Table 5.4. Reductions in
Young’s moduli of 29% (E50 = 70 MPa→50 MPa) and 57% (E50 = 70 MPa→30 MPa) yield
average reductions in stress concentration ratios of 12% and 31%, respectively.
Table 5.4 - Influence of column compressibility upon stress concentration ratios for end-bearing columns
Stress concentration ratio, σ'col / σ'soil % reduction
No. of cols Column compressibility, E50 (MPa) Column compressibility, E50 (MPa)
30 50 70 30 50 70
4 4.6 5.7 6.4 -28 -11 -
5 4.3 5.4 6.1 -30 -11 -
9 4.1 5.3 6.3 -35 -15 -
4 cols - s = 2.0 m
10
E_50 = 30 MPa
12 E_50 = 50 MPa
E_50 = 70 MPa
14
(a)
2
Column length, L (m)
5 cols - s = 1.0 m
10
E_50 = 30 MPa
12 E_50 = 50 MPa
E_50 = 70 MPa
14
(b)
2
Column length, L (m)
9 cols - s = 1.0 m
10
E_50 = 30 MPa
12 E_50 = 50 MPa
E_50 = 70 MPa
14
(c)
Figure 5.48 - Influence of column compressibility upon the variation of stress concentration ratios with column
length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 148
variation of stress concentration ratio with column length for groups of 4 and 5 columns is
shown in Figures 5.49(a) and 5.49(b), respectively. Negligible changes in stress concentration
ratios are observed with increasing column strength for short columns. However, it can be
seen that increasing the strength of columns which are longer than L = 3 m leads to
significant increase in stress concentrations ratios. This suggests that the influence of column
strength upon stress concentration ratios is dependent on the mode of deformation. Column
strength is controlled by the angle of internal friction and this parameter is only mobilised
when columns are in a plastic state. Columns which are punching do not mobilise the column
strength as stresses are concentrated along the column-soil interface and at the base of
columns. Therefore, stress concentration ratios do not increase with column strength for
columns which exhibit punching.
The variation of stress concentration ratio with column length for a group of 9 columns is
shown in Figure 5.49(c). As with the groups of 4 and 5 columns, no variation of stress
concentration ratios is observed with column strength for short columns. However, a decrease
in column strength, from φ = 45° to 40°, leads to a reduction in stress concentration ratios.
Stress concentration ratios are relatively constant with increasing column length for columns
longer than L = 3 m and defined by φ = 40°. This suggests that this configuration of columns
is bulging. In contrast, a slight increase in stress concentration ratios with column length is
observed for columns longer than L = 3 m and defined by φ = 45°. This is consistent with a
punching mode of deformation. Furthermore, no variation of stress concentration ratios is
observed when increasing column strength from φ = 45° to 50°. Therefore, the stronger group
of columns also exhibits punching and, consequently, stress concentration ratios are not
influenced by a change in the column strength.
8
4 cols - s = 2.0 m
10
φ = 40°
12 φ = 45°
φ = 50°
14
(a)
2
Column length, L (m)
8
5 cols - s = 1.0 m
10
φ = 40°
12 φ = 45°
φ = 50°
14
(b)
2
Column length, L (m)
8
9 cols - s = 1.0 m
10
φ = 40°
12 φ = 45°
φ = 50°
14
(c)
Figure 5.49 - Influence of column strength upon the variation of stress concentration ratios with column length
for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 150
8
4 cols - s = 2.0 m
10
K0_soil = 0.75
12 K0_soil = 1.00
K0_soil = 1.25
14
(a)
2
Column length, L (m)
8
5 cols - s = 1.0 m
10
K0_soil = 0.75
12 K0_soil = 1.00
K0_soil = 1.25
14
(b)
2
Column length, L (m)
8
9 cols - s = 1.0 m
10
K0_soil = 0.75
12 K0_soil = 1.00
K0_soil = 1.25
14
(c)
Figure 5.50 - Influence of column installation effects upon the variation of stress concentration ratios with
column length for a group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 151
A further increase in stress concentration ratios is observed for columns formed in Profile 3.
This is due to the absence of the stiff crust at the ground surface. The soil at the ground
surface in Profile 3 consists of upper Carse clay, which yields a higher modular ratio
(Ecol/Esoil) than Profile 1. This results in higher levels of vertical stress in the stone columns,
which consequently leads to higher stress concentration ratios.
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 152
10 4 cols - s = 2.0 m
Profile 1
12 Profile 2
Profile 3
14
(a)
2
Column length, L (m)
10 5 cols - s = 1.0 m
Profile 1
12 Profile 2
Profile 3
14
(b)
2
Column length, L (m)
10 9 cols - s = 1.0 m
Profile 1
12 Profile 2
Profile 3
14
(c)
Figure 5.51 - Influence of stiff crust upon the variation of stress concentration ratios with column length for a
group of (a) 4, (b) 5 and (c) 9 columns beneath a 3 m square footing
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 153
The influence of these parameters is related to the degree of plasticity within columns,
and hence the mode of deformation.
• The coefficient of lateral earth pressure has a positive influence on the settlement
performance of stone columns and this is most pronounced at high area ratios.
• The presence of the stiff crust has a significant influence upon the settlement
performance of stone columns. It was found that the improvement in settlement
improvement factors tails of significantly for homogenous soils samples, i.e. without a
stiff crust. This profile is similar to that adopted in laboratory studies (see Sections 2.2
and 2.3) and may explain why authors have postulated critical lengths in the past.
However, critical lengths are less well defined for more realistic soil profiles.
Three modes of deformation were identified from the numerical analysis: (i) punching; (ii)
‘block’ failure and (iii) bulging. The punching mode of deformation is observed for all short
columns (i.e. L < 3 m), as an increase in punching ratios is coupled with relatively low
compression ratios. Punching remains the primary mode of deformation for long single
columns at low area ratios, as relatively high punching and low compression ratios are
observed with increasing column length. ‘Block’ failure is characterised by low punching
and compression ratios and this mode of deformation is observed for large groups of columns
at low area ratios. ‘Block’ failure is effectively an extension of punching, as columns act as a
single entity with the surrounding soil and punch uniformly, thus yielding low punching
ratios. The bulging mode of deformation is observed for long columns at high area ratios, as a
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 155
Stress concentration ratios do not increase with column length for an infinite grid of columns.
This is due to symmetry which ensures that the behaviour of an infinite grid of columns is
Results of FEA: Settlement performance, deformational behaviour and stress concentration ratios 156
identical and, consequently, no shear stress develops along the side of the columns near the
ground surface. Therefore, stress concentration ratios are only dependent upon the modular
ratio (Ecol/Esoil). In contrast, stress concentration ratios for a group of stone columns are
dependent upon the mode of deformation. Columns which are punching develop low stress
concentration ratios, which increase with column length. However, no increase in stress
concentration ratios is observed with increasing column length for columns which exhibit a
bulging mode of deformation. Bulging occurs at the point of lowest lateral resistance and
increasing the column length does not enhance the load-carrying capacity of stone columns.
Chapter 6
Results of FEA: Distribution of total shear
strains and characteristic column
behaviour
6.1 Background
The characteristic behaviour of individual columns is examined in the following chapter. The
load-transfer mechanisms are identified for various column configurations at three specific
column lengths: 3, 8 and 13.9 m. The column lengths are chosen on the basis of previous
findings from the deformational behaviour of columns (see Section 5.4) to examine specific
modes of deformation; 3 m long columns are chosen to investigate punching; 8 m long
columns are chosen to investigate a combination of punching-bulging failure and 13.9 m end-
bearing columns are chosen to investigate bulging failure, as punching is precluded.
It is established that the deformational behaviour of columns is highly dependent upon the
support provided by the surrounding soil. The type of stress imparted from columns to the
surrounding soil depends upon the dominant mode of deformation, i.e. columns which are
bulging primarily impart radial stress, while columns which are punching impart a
combination of shear and end-bearing stress. The response of the surrounding soil to each of
these stresses plays a key role in determining the settlement performance of stone columns.
The distribution of shear strains within columns and in the surrounding soil is therefore very
important. This is investigated in the following sections and allows the modes of deformation
identified in Section 5.4, which were determined by compression and punching ratios, to be
verified.
The mode of deformation for 8 m long floating columns is shown in Figure 6.1(b). It can be
seen in Figure 6.1(b-i) that the majority of shear strain continues to develop beneath the base
of columns at low area ratios. It appears that an infinite grid of columns at low area ratios do
not develop shear stress along the sides of columns and transfer the applied load to the base
of columns, which is again consistent with ‘block’ failure. A change in the mode of
deformation is observed for 8 m long columns at higher area ratios in Figures 6.1(b-ii) and
6.1(b-iii). While a fluctuation in shear strain is observed within columns, it is clear that
bulging develops at the top of the lower Carse clay layer (i.e. weakest part of the soil profile).
Columns at higher area ratios experience a loss of lateral confinement which results in
column bulging. The loss of lateral confinement becomes more pronounced with increasing
area ratio and the magnitude of shear strain increases in concert.
No punching is observed in Figure 6.1(c) as the end-bearing columns (L = 13.9 m) are resting
on a rigid stratum. It can be seen that bulging always occurs at the top of the lower Carse clay
layer for all area ratios. The magnitude of column bulging increases at high area ratios, which
reflects the loss of lateral confinement and, also, the higher force taken per column with
increasing area ratio.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 159
The distribution of total shear strain for 8 m long floating columns is shown in Figure 6.2(b).
It can be seen in Figure 6.2(b-i) that shear strain develops along the side and at the base of
columns at low area ratios, which indicates that punching is the dominant mode of
deformation. It can also be seen that the magnitude of shear strain at the base of columns is
considerably reduced compared to 3 m long columns. This may be explained as longer
columns can generate a larger shear force along the side of columns and also a higher base
resistance is available as columns are formed in a deeper, more competent soil stratum. It can
be seen in Figures 6.2(b-ii) and 6.2(b-iii) that the mode of deformation changes from
punching to bulging with increasing area ratio for 8 m long columns. While some shear strain
develops along the side and at the base of columns, the majority of shear strain develops
within columns near the top of the lower Carse clay layer.
The distribution of shear strains for single end-bearing columns is shown in Figure 6.2(c) and
indicates that no bulging occurs in columns at low area ratios. It appears that columns at low
area ratios develop shear strains along the sides of columns in upper sections. In contrast, the
majority of shear strain develops within columns at higher area ratios, which suggests that
bulging is the primary mode of deformation.
Localised shear zones also develop beneath the edges of pad footings for columns at higher
area ratios. These shear zones develop as a results of the large displacement discontinuity
between the pad footing and the surrounding soil; however, they do not appear to extend
beyond the stiff crust.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 161
Bulging at top of
Shear strain develops ‘lower Carse clay’
along sides and at layer
base of column
As column length increases, Figure 6.3(b-i) illustrates how the magnitude of shear strain
beneath the base of columns reduces and a more uniform distribution of shear strain develops
along the side and at the base of columns. It can be also seen that no shear strain develops in
the central zone of soil bounded by the columns which again suggests that 'block' failure
remains the dominant mode of deformation with increasing column length for 2×2 groups at
low area ratios. In contrast, the mode of deformation appears to change from punching to
bulging with increasing column length for 2×2 groups of columns at higher area ratios.
Column bulging is evident in Figures 6.3(b-ii) and 6.3(b-iii) as the majority of the shear strain
develops within columns, near the top of the lower Carse clay layer. Shear zones also develop
from the corners of pad footings to the point of maximum bulging in columns at high area
ratios. It also appears that shear zones bridge between the points of maximum bulging in
columns and create a zone of undeforming soil between the base of the footing and the points
of column bulging. This zone of undeforming soil was also observed by Muir Wood et al.
(2000) and Wehr (2004) in small scale laboratory tests and numerical studies, respectively.
This zone of soil moves downwards and displaces the surrounding soil, causing the columns
to bulge and bend outwards.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 163
As A/AC ↑
Uniform strain at Localised failure
base of columns
⇒ ‘block’ failure
at column bases
Development of
shear zone
Low γ in soil
bounded by outer
columns
Development of
Shear strain shear zone
develops in soil
The distribution of shear strains for end-bearing columns is shown in Figure 6.3(c). While no
bulging occurs in columns at low area ratios, it appears from Figure 6.3(c-i) that shear strain
develops in the soil bounded by the columns. This causes columns to bend outwards and
away from the centre of the footing. Shear strains are also observed along the side of
columns, which indicates that the columns may still be transferring the applied load through
side friction. Figures 6.3(c-ii) and 6.3(c-iii) indicate that bulging failure continues to occur for
end-bearing columns at high area ratios. The development of a shear zone, which extends
from the corners of pad footings to the points of maximum bulging, and bridges the soil
between columns is again evident.
The punching mode of failure continues to occur with increasing length for columns at low
area ratios, as shown in Figure 6.4(b-i) for 8 m long columns. While some shear strain
develops within the upper sections of external column, which may result in column bending,
the majority of shear strain develops along the outer side of external columns and at the base
of columns. Therefore, it appears that ‘block’ failure remains the dominant mode of
deformation with increasing column length for a group of columns at low area ratios. In
contrast, it can be seen in Figures 6.4(b-ii) and 6.4(b-iii) that columns at higher area ratios
tend to bulge with increasing column length. Bulging occurs at the top of the lower Carse
clay layer. The shear zone which extends from the corners of pad footings to the points of
maximum bulging in external columns becomes more defined with increasing column length.
The shear zone bridges between the points of bulging in the central and outer columns, which
creates a zone of undeforming soil.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 165
Development of
shear zone
Low γ in soil
bounded by outer
columns
The mode of deformation changes from punching to bulging once columns are end-bearing
on a rigid stratum. In contrast to single columns and groups of 2×2 columns, bulging occurs
in 3×3 groups of end-bearing columns at low area ratios (Figure 6.4(c-i)). It can be seen in
Figures 6.4(c-ii) and 6.4(c-iii) that the bulging mode of deformation continues to occur with
increasing column length for 3×3 groups of columns at high area ratios. Furthermore, shear
strain in external columns is more acute, which highlights the positive confining effect
provided to the central column. It also appears that larger deformation occurs at a slightly
shallower depth in external columns compared to central columns. Initially all columns tend
to bulge at the same depth, i.e. at the top of the lower Carse clay. However, the bulging of
external and central columns mobilises the passive resistance of the soil which is bounded by
the central and exterior columns. This increases the lateral resistance at this depth, which
enhances the confinement provided to the central column, and hence forces bulging deeper.
The magnitude of bulging is also reduced, as bulging occurs in a deeper and stiffer soil
stratum. Bulging in external columns continues to occur at the top of the lower Carse clay
layer, as columns bulge and bend outwards towards the unconfined sides.
It can be seen in Figure 6.5(b-i) that punching remains the dominant mode of deformation
with increasing column length for columns at low area ratios. Shear strains develop beneath
the base of columns and along the outer side of external columns. It is also clear that low
levels of shear strain develops in central columns and in the soil bounded by external
columns, which provides further evidence to suggest that columns at low area ratios are
deforming as a ‘block’. However, column bulging becomes more pronounced with increasing
Results of FEA: Distribution of total shear strains and characteristic column behaviour 167
Development of
shear zone
area ratio, as shown in Figures 6.5(b-ii) and 6.5(b-iii). Bulging occurs at the top of the lower
Carse clay layer in external columns and occurs slightly deeper in central columns. The
development of a shear zone, which extends from the corner of pad footings to the point of
maximum bulging, is again observed.
The distribution of shear strain for end-bearing columns is shown in Figure 6.5(c). It can be
seen that bulging occurs for all area ratios as column punching is precluded. It appears that
the magnitude and location of bulging is similar for central and external columns at low area
ratios. However, as the area ratio increases it can be seen that the magnitude of shear strain
increases and bulging occurs at shallower depths in external columns compared to central
columns. The development of shear zones from the corner of the pad footing is again evident
and they join together beneath the centre of the footing to form a zone of undeforming soil.
horizontal strain increases with area ratio, which reflects the loss of lateral confinement and
increase in the total load taken per column with increasing area ratio. In addition, it is
interesting to note that no horizontal strain develops in the soil below the base of columns.
This indicates that the unreinforced soil is subject to one-dimensional boundary conditions,
imposed by the wide area loading.
The distribution of vertical strain with depth for an infinite grid of 8 m long columns is
shown in Figure 6.6(b-i). It can be seen for columns at low area ratios (i.e. A/AC = 3.5) that
the majority of vertical strain develops beneath the base of columns. In contrast, significant
vertical strain develops within columns at higher area ratios. It is interesting to note that
vertical strain within columns increases with area ratio, while the vertical strain beneath the
base of columns remains unchanged. This may be explained as the wide area loading stresses
the soil profile to the full depth and, therefore, induces similar vertical strain in the soil
beneath the base of columns. As with 3 m long columns, it appears that the distribution of
vertical strain within 8 m long columns is related to the horizontal strain (Figure 6.6(b-ii)). It
appears that both vertical and horizontal strain within columns increases with area ratio and
that the maximum horizontal strain occurs at the same depth (z ≈ 2.5 m), close to the top of
the lower Carse clay layer (i.e. the weakest point in the soil profile). A significant amount of
horizontal strain occurs within columns at high area ratios (i.e. A/AC = 8.0 and 14.1) which is
consistent with column bulging.
The distribution of vertical strain for an infinite grid of end-bearing stone columns is shown
in Figure 6.6(c-i). Vertical strain increases with depth to a maximum at z ≈ 2.5 m and
decreases thereafter, eventually becoming negligible at z = 8 m. The horizontal strain is
shown in Figure 6.6(c-ii) and is similar to the distribution of vertical strain with depth. This
indicates that vertical strains are primarily developed from horizontal strains.
4
Depth, z (m)
2 σ'col/σ'soil ↑ at
transition to
softer layers
4
Depth, z (m)
4
Depth, z (m)
Figure 6.6 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for an infinite grid of columns at lengths L of (a) 3 m, (b) 8 m and (c) 13.9 m
Results of FEA: Distribution of total shear strains and characteristic column behaviour 171
Depth, z (m)
6
10
A/Ac = 3.5
12 A/Ac = 8.0
A/Ac = 14.1
14
Figure 6.7 - Distribution of modular ratio Ecol/Esoil with depth for an infinite grid of end-bearing stone columns
at various area ratios
with depth according to the power law outlined in equation 3.3. The rate of increase in
stiffness with overburden pressure is controlled by the parameter, ‘m’, which is defined as 0.3
and 1.0 for the column and soil, respectively. The higher ‘m’ value results in a faster increase
in Esoil compared to Ecol with depth. Therefore, Ecol/Esoil reduces with depth in the crust,
which leads to lower stress concentration ratios.
An increase in stress concentration ratios is observed at the transition between the crust and
the softer upper Carse clay layer. This may be explained as Ecol/Esoil increases sharply for this
layer and columns, therefore, carry a larger proportion of the applied load. A slight reduction
in stress concentration ratios with depth is again observed throughout the upper Carse clay
layer, followed by a sharp increase at the transition to the softer lower Carse clay layer. Stress
concentration ratios continue to reduce with depth in the lower Carse clay layer and revert to
unity (i.e. σcol/σsoil = 1) at the base of floating columns, as no extra vertical load can be
sustained in the columns at this level.
It can be seen in Figure 6.6(c-iii) that stress concentration ratios are related to the modular
ratio for an infinite grid of columns at A/AC = 3.5, as stress concentration ratios reduce by
c.55% in the lower Carse clay layer (σcol/σsoil = 5.0 to 2.3) which coincides with a similar
drop in modular ratios for the same layer (Ecol/Esoil = 160 to 68), as shown in Figure 6.7.
While a similar distribution of stress concentration ratio with depth is observed for all short
columns (Figure 6.6(a-iii)), a marked change is observed with increasing area ratio for 8 m
Results of FEA: Distribution of total shear strains and characteristic column behaviour 172
and 13.9 m long columns, as shown in Figures 6.6(b-iii) and 6.6(c-iii), respectively. Despite a
large degree of scatter, it appears that stress concentration ratios reduce with increasing
horizontal strain. It can be seen in Figure 6.6(ii) that horizontal strain increases with
increasing area ratio. This reflects the loss of lateral confinement and, also, the higher load
carried per column with increasing area ratio. It can be seen that sections of columns which
develop high horizontal strain cannot sustain the vertical stress and, consequently, yield low
stress concentration ratios. However, it is interesting to note that in sections of columns
which develop low horizontal strain, stress concentration ratios increase for higher area ratios
(Figure 6.6(c-iii)).
The distribution of vertical and horizontal strain along 8 m long columns is shown in Figures
6.8(b-i) and 6.8(b-ii), respectively. A marked changed in the distribution of vertical and
horizontal strain within columns is observed with increasing area ratio. It can be seen in
Figure 6.8(b-i) for low area ratios that the majority of vertical strain develops beneath the
base of columns, which suggests that punching remains the dominant mode of deformation.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 173
4
Depth, z (m)
2 σ'col/σ'soil ↓ in
section which
4 is bulging
Depth, z (m)
2 σ'col/σ'soil ↓ in
section which
4 is bulging
Depth, z (m)
Figure 6.8 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a single column at lengths L of (a) 3 m, (b) 8 m and (c) 13.9 m
Results of FEA: Distribution of total shear strains and characteristic column behaviour 174
It is interesting to note that the magnitude of vertical strain at the base of columns is reduced
compared to 3 m long columns. This may be explained as the base of 8 m long columns is
founded in a deeper, more competent soil stratum. It can also be seen in Figure 6.8(b-i) for
columns at high area ratios that vertical strain develops within columns and beneath the base
of columns. This suggests that punching and bulging modes of deformation may be occurring
simultaneously. However, it can be seen that the majority of the vertical strain develops
within columns rather than at the base, which suggests that bulging is the dominant mode of
deformation. The development of high horizontal strains within columns is also evident in
Figure 6.8(b-ii), which is consistent with a bulging mode of deformation. It can be seen in
Figure 6.8(b-i) that the magnitude of vertical strain and the extent of the column which is
strained both increase at high area ratios. This is due to the loss of lateral confinement and
increase in the total load taken per column. Maximum bulging occurs at a depth between 2.0–
2.5 m, which coincides with the top of the lower Carse clay layer and is consistent with the
depth of maximum bulging for an infinite grid of columns.
The distributions of vertical and horizontal strain within end-bearing columns are shown in
Figures 6.8(c-i) and 6.8(c-ii), respectively. It can be seen that no horizontal strain develops in
single columns at low area ratios which indicates that columns are not bulging. However, a
small amount of vertical strain is observed near the top of columns at low area ratios. This
indicates that columns are transferring the applied load through shear stress along the sides of
columns, which is similar to a punching mode of deformation. The bulging mode of
deformation remains prevalent for columns at higher area ratios, as a significant degree of
vertical and horizontal strain is observed in the upper sections of columns. Similar to 8 m
long floating columns, the magnitude of strain and the extent of the column which is strained
both increase with increasing area ratio.
stress concentration ratio, reduces with depth. Stress concentration ratios approach unity at
the base of columns as no extra vertical stress in columns can be sustained at this point.
The distribution of stress concentration ratios with depth for 8 m and 13.9 m long columns is
shown Figures 6.8(b-iii) and 6.8(c-iii), respectively. It can be seen that stress concentration
ratios increase with area ratio for sections of columns which are subject to low horizontal
strains. This may be attributable to an increased load that is taken per column with increasing
area ratio and is consistent with the findings from infinite grids of columns. However, it
appears that the magnitude of stress concentration ratios is significantly reduced in the upper
sections of columns, where column bulging occurs. Columns in a state of plasticity cannot
sustain high vertical stress and, consequently, a significant reduction in vertical stress, and
hence stress concentration ratios, occurs. The reduction of stress concentration ratios appears
to be most pronounced at the top of the bulged section of columns.
4
Depth, z (m)
2 σ'col/σ'soil ↓ in
section which
4 is bulging
Depth, z (m)
2 σ'col/σ'soil ↓ in
section which
4 is bulging
Depth, z (m)
Figure 6.9 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 2×2 group of columns at lengths L of (a) 3 m, (b) 8 m and (c) 13.9 m
Results of FEA: Distribution of total shear strains and characteristic column behaviour 177
The distribution of vertical and horizontal strain with depth for 2×2 groups of 8 m long
columns is shown in Figures 6.9(b-i) and 6.9(b-ii), respectively. It can be seen in Figure
6.9(b-i) for low area ratios that the majority of vertical strain develops beneath the base of
columns. This is again consistent with a punching mode of deformation. In contrast, it can be
seen that the majority of vertical strain develops within columns at high area ratios, which
suggests that bulging is the dominant mode of deformation. It can be seen in Figure 6.9(b-ii)
that columns at high area ratios develop large horizontal strains, which is consistent with a
bulging mode of deformation.
The distribution of vertical and horizontal strain with depth for 2×2 groups of end-bearing
columns is shown in Figures 6.9(c-i) and 6.9(c-ii), respectively. It can be seen that the
vertical and horizontal strain in columns at low area ratios is quite low, which suggests that
columns do not fail by bulging. This is quite interesting as punching failure is precluded for
end-bearing columns. It appears that closely-spaced columns develop shear stress along the
side of the columns to resist the applied load. A significant increase in horizontal strain is
observed with an increasing area ratio, which indicates that columns at high area ratios are
bulging.
A change in the distribution of stress concentration ratio at higher area ratios (A/AC = 8.0 and
14.1) is observed as column length increases. It can be seen in Figures 6.9(b-iii) and 6.9(c-iii)
that stress concentration ratios reduce with depth in the stiff crust (i.e. z < 0.9). This again
may be explained as columns carry the applied load through shear stress along the sides of
the columns in this section. However, the magnitude of stress concentration ratios at this
depth is lower than stress concentration ratios for closely-spaced columns (i.e. A/AC = 3.5),
which is a direct result of the bulged column in the Carse clay. Columns which are bulging
Results of FEA: Distribution of total shear strains and characteristic column behaviour 178
are in a contained state of plasticity which reduces their load-carrying capacity and, therefore,
the vertical stress in bulged sections of columns. The vertical stress increases with depth in
bulged sections of columns due to overburden stress, which increases column confinement
and hence allows columns to carry a larger proportion of the load. Stress concentration ratios
decrease thereafter as vertical stress is carried through shear stress and also as the modular
ratio reduces with increasing overburden stress.
The distribution of vertical and horizontal strain with depth for 3×3 groups of 8 m long
columns is shown in Figures 6.11(i) and 6.11(ii), respectively. It can be seen in Figure 6.11(a-
i) for 8 m long columns at low area ratios that the majority of vertical strain develops beneath
the base of columns, which suggests that punching remains the dominant mode of
deformation. A change in the mode of deformation is observed with increasing area ratio for
8 m long columns. A significant proportion of horizontal strain is observed within columns at
area ratios of 8.0 and 14.1 in Figures 6.11(b-ii) and 6.11(b-iii), respectively. While a certain
amount of vertical strain develops at the base of columns, it is clear that column bulging is
the dominant mode of deformation. The variation in strain levels with column position is only
observed for columns at low area ratios, where the positive effects of column confinement are
again evident as vertical and horizontal strains are lowest in central columns.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 179
Stress concentation,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -5 -10 -15 -20 0 3 5 8 10 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 3.5 A/Ac = 3.5 A/Ac = 3.5
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(a-i) (a-ii) (a-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 -40 0 5 10 15 20 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 8.0 A/Ac = 8.0 A/Ac = 8.0
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(b-i) (b-ii) (b-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Radial strain, εr (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 -40 0 5 10 15 20 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 14.1 A/Ac = 14.1 A/Ac = 14.1
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(c-i) (c-ii) (c-iii)
Figure 6.10 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 3 m long column, within a 3×3 group for A/Ac of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Distribution of total shear strains and characteristic column behaviour 180
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -2 -4 -6 -8 -10 0 1 2 3 4 5 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 3.5 A/Ac = 3.5 A/Ac = 3.5
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(a-i) (a-ii) (a-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -5 -10 -15 -20 0 5 10 15 20 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 8.0 A/Ac = 8.0 A/Ac = 8.0
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(b-i) (b-ii) (b-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 -40 0 10 20 30 40 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 14.1 A/Ac = 14.1 A/Ac = 14.1
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(c-i) (c-ii) (c-iii)
Figure 6.11 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 8 m long column, within a 3×3 group for A/AC of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Distribution of total shear strains and characteristic column behaviour 181
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -1 -2 -3 -4 -5 0 1 2 3 4 5 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 3.5 A/Ac = 3.5 A/Ac = 3.5
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(a-i) (a-ii) (a-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -5 -10 -15 -20 0 5 10 15 20 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 8.0 A/Ac = 8.0 A/Ac = 8.0
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(b-i) (b-ii) (b-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 -40 0 10 20 30 40 1 3 5 7 9
0
4
Depth, z (m)
8
3×3 cols 3×3 cols 3×3 cols
10 A/Ac = 14.1 A/Ac = 14.1 A/Ac = 14.1
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(c-i) (c-ii) (c-iii)
Figure 6.12 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 13.9 m long column, within a 3×3 group for A/AC of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Distribution of total shear strains and characteristic column behaviour 182
The distribution of vertical and horizontal strain for 13.9 m long columns is shown in Figures
6.12(i) and 6.12(ii), respectively. Punching failure is precluded for end-bearing columns and
bulging appears to be the dominant mode of deformation for all columns. A significant
increase in the magnitude of vertical and horizontal strain is observed with increasing area
ratio, which may be attributed to the loss of lateral confinement and increase in the total load
taken by columns with increasing area ratio. No variation in the levels of vertical and
horizontal strain is observed with column position for end-bearing columns.
The distribution of stress concentration ratios for 8 m long columns is shown in Figure
6.11(iii). It can be seen for closely-spaced columns in Figure 6.11(a-iii) that stress
concentration ratios reduce steadily with depth, which is consistent with a punching mode of
deformation. However, a change in the variation of stress concentration ratios with depth is
observed for higher area ratios (Figures 6.11(b-iii) and 6.11(c-iii)). It appears that stress
concentration ratios for columns at high area ratios are related to horizontal strain. A
significant proportion of horizontal strain occurs in long columns at high area ratios which
reduces the load-carrying capacity of columns and hence results in lower stress concentration
ratios. This is similar to findings for single columns and 2×2 groups. The largest drop in
stress concentration ratios is evident at the top of the bulged section of columns. It can be
seen that stress concentration ratios increase with depth to the bottom of the bulged section,
due to overburden stress, and then decrease to unity at the base of floating columns.
The variation of stress concentration ratios with column position is also evident in Figure
6.11(iii). Vertical stress at the surface is highest for corner columns, followed by edge and
Results of FEA: Distribution of total shear strains and characteristic column behaviour 183
centre columns. This is again attributable to the position of columns relative to the footing
edge. No change in stress concentration ratio with column position is observed with depth
thereafter for closely-spaced columns. However, it can be seen for columns at higher area
ratios that stress concentration ratios are highest for centre columns, followed by edge and
corner columns over the bulged section of columns. This highlights the positive effects of
column confinement provided to centre columns.
A similar variation of stress concentration ratios for 8 m long columns is observed for end-
bearing columns in Figure 6.12(iii). It can be seen for closely-spaced columns in Figure
6.12(a-iii) that stress concentration ratios reduce steadily with depth, which is consistent with
a punching mode of deformation. The distribution of stress concentration ratios for columns
at higher area ratios is shown in Figures 6.12(b-iii) and 6.12(c-iii). Similar to single columns
and 2×2 groups, it can be seen that the magnitude of stress concentration ratios reduces with
increasing levels of horizontal strain.
The distribution of vertical and horizontal strain for 8 m long columns is shown in Figures
6.14(i) and 6.14(ii), respectively. It can be seen in Figure 6.14(a-i) that punching remains the
dominant mode of deformation for closely-spaced columns, as the majority of vertical strain
occurs beneath the base of columns. However, a change in the mode of deformation is
observed with increasing area ratio. It can be seen in Figure 6.14(ii) that considerable
horizontal strain occurs in columns, which indicates that column bulging is occurring.
Results of FEA: Distribution of total shear strains and characteristic column behaviour 184
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 -40 0 5 10 15 20 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 3.5 A/Ac = 3.5 A/Ac = 3.5
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(a-i) (a-ii) (a-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -20 -40 -60 -80 0 10 20 30 40 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 8.0 A/Ac = 8.0 A/Ac = 8.0
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(b-i) (b-ii) (b-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -30 -60 -90 -120 0 10 20 30 40 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 14.1 A/Ac = 14.1 A/Ac = 14.1
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(c-i) (c-ii) (c-iii)
Figure 6.13 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 3 m long column, within a 4×4 group for A/AC of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Distribution of total shear strains and characteristic column behaviour 185
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -2 -4 -6 -8 0 1 2 3 4 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 3.5 A/Ac = 3.5 A/Ac = 3.5
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(a-i) (a-ii) (a-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 0 5 10 15 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 8.0 A/Ac = 8.0 A/Ac = 8.0
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(b-i) (b-ii) (b-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -15 -30 -45 -60 0 10 20 30 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 14.1 A/Ac = 14.1 A/Ac = 14.1
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(c-i) (c-ii) (c-iii)
Figure 6.14 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 8 m long column, within a 4×4 group for A/AC of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Distribution of total shear strains and characteristic column behaviour 186
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -2 -4 -6 -8 0 1 2 3 4 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 3.5 A/Ac = 3.5 A/Ac = 3.5
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(a-i) (a-ii) (a-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -10 -20 -30 0 5 10 15 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 8.0 A/Ac = 8.0 A/Ac = 8.0
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(b-i) (b-ii) (b-iii)
Stress concentration,
Vertical strain, εy (× 10-3) Horizontal strain, εh (× 10-3) σ'col / σ'soil, average
0 -20 -40 -60 -80 0 10 20 30 40 1 3 5 7 9
0
4
Depth, z (m)
8
4×4 cols 4×4 cols 4×4 cols
10 A/Ac = 14.1 A/Ac = 14.1 A/Ac = 14.1
Centre Centre centre
12 Edge Edge edge
Corner Corner corner
14
(c-i) (c-ii) (c-iii)
Figure 6.15 - Distribution of (i) vertical strain, (ii) horizontal strain and (iii) stress concentration ratio with depth
for a 13.9 m long column, within a 4×4 group for A/AC of (a) 3.5, (b) 8.0 and (c) 14.1
Results of FEA: Distribution of total shear strains and characteristic column behaviour 187
While it appears that column bulging is the dominant mode of deformation, it can still be
seen that some of the applied load is transferred to the base of columns as vertical strains are
developed in this region. The variation of strain level with column position is clearly visible
for closely-spaced columns. The positive influence of column confinement is evident as
vertical and horizontal strains are lowest for centre columns. The variation of strain level with
column position is less well defined as area ratio increases. It appears that centre columns
exhibit slightly lower vertical and horizontal strain.
The distribution of vertical and horizontal strain for end-bearing columns is shown in Figures
6.15(i) and 6.15(ii), respectively. It can be seen that bulging is the dominant mode of
deformation for all area ratios. A significant increase in vertical and horizontal strain is also
evident with increasing area ratio. The variation of strain levels with column position can also
be seen for columns at higher area ratios. It can be seen in Figure 6.15(b-i) that the magnitude
of vertical strain is highest and maximum strain occurs at the deepest level for centre columns
spaced at A/AC = 8.0. The maximum vertical strain also occurs deeper for centre columns
which are spaced at A/AC = 14.1. However, the maximum vertical and horizontal strain is
highest for corner columns.
The distribution of stress concentration ratios for 8 m long columns is shown in Figure
6.14(iii). The distribution of stress concentration ratio with depth for closely-spaced columns
again reduces steadily with depth. This indicates that columns are punching into the
underlying soil and is consistent with previous findings from single columns, 2×2 groups and
3×3 groups. However, a marked change in the distribution of stress concentration ratios is
observed with increasing area ratio. It appears that stress concentration ratios are related to
the horizontal strains. The development of horizontal strain in columns reduces the load-
carrying capacity of columns, which results in lower stress concentration ratios. It is also
Results of FEA: Distribution of total shear strains and characteristic column behaviour 188
interesting to note that while stress concentration ratios are lowest for centre columns at the
surface, a reverse in this trend is observed in lower sections of columns which are bulging.
The distribution of stress concentration ratios for end-bearing columns is shown in Figure
6.15(iii). It can be seen that stress concentration ratios follow a similar trend to 3×3 groups of
columns and reduce steadily with depth for closely-spaced columns. However, it can be seen
that the magnitude of stress concentration ratios reduces for higher area ratios, as vertical and
horizontal strain increase. It can be seen that stress concentration ratios increase for centre
columns in sections of columns which are bulging.
within columns. This mode of deformation occurs for long columns at high area ratios. In
addition to the three modes of deformation, the development of shear planes was also
observed. The shear planes extend from the corner of pad footings to the point of maximum
bulging and create a zone of undeforming soil. This is consistent with previous findings by
Muir Wood et al. (2000); however, it is important to note that the shear planes are a
secondary mode of deformation, but may become more pronounced at higher load levels.
6.4.2 Distribution of vertical strain, horizontal strain and stress concentration ratios
The distribution of vertical and horizontal strain is consistent with the shear strains. It was
also observed that punching or bulging modes of deformation are not clearly defined and
some combination of both modes of deformation occurs for all configurations of columns.
The distribution of stress concentration ratios is significantly influenced by the vertical and
horizontal strain within columns. Sections of columns which exhibit large vertical strains are
in a contained state of plasticity and cannot sustain high vertical loads. Consequently, stress
concentration ratios reduce in these regions.
190
Chapter 7
Development of simplified design method
and comparison of findings with previous
research
7.1 Introduction
The results from the numerical analysis regarding the settlement performance, deformational
behaviour and stress concentration ratios for small groups of stone columns are placed in
context to previous research in this chapter.
The loss of lateral confinement for small groups of end-bearing stone columns has been
studied previously, and the findings of these studies are compared with those obtained using
PLAXIS 3D Foundation. The previous research is extended through investigation of the
influence of column confinement upon the settlement performance of small groups of floating
stone columns. A simplified design method is proposed which allows the settlement of small
groups and infinite grids of columns to be related with knowledge of the footing-width to
column-length ratio and column-length to soil-thickness ratio.
Three modes of deformation referred to as "punching", "block failure" and "bulging", were
identified in the previous chapters for small groups of stone columns. The development of
these modes of deformation is compared with small scale laboratory studies, conducted by
Muir Wood et al. (2000) and Black (2006). The importance of capturing the influence of the
stiff crust, a salient feature of soft soil profiles, is highlighted and the existence of a critical
length, proposed by previous laboratory research, is investigated.
The influence of the mode of deformation upon the settlement performance of stone columns
is also investigated in the following chapter. The modes of deformations were identified
using compression and punching ratios and the variation of these parameters with settlement
improvement factors for different configurations of columns is presented and analysed.
Development of simplified design method and comparison of findings with previous research 191
Stress concentration ratios at the surface of an infinite grid of columns are compared with a
database, which consists of field and laboratory measurements, and also analytical design
methods. The distribution of stress concentration ratios with depth in not well understood and
this is examined and compared with analytical design methods.
• The design curves are based on numerous calculations which account for the reduction in
vertical stress with depth beneath small loaded areas and also consider a lower bearing
capacity for the outer ring of columns in the column group.
• No direct reference is made to the footing area (A) in Figure 2.24 and it appears that s/suc
is independent of the area ratio (A/AC).
• However, for a given number of columns, an increase in footing area leads to larger
settlements (s) due to the greater depth of stress. The corresponding increase in A/AC
leads to an increase in suc. Priebe (1995) claims that the increase in s and suc
"compensate" such that s/suc ratios are acceptable for A/AC < 10.
The variation of s/suc with B/L for the two soil profiles is shown in Figure 7.1:
• It can be seen that s/suc values increase somewhat asymptotically with B/L for profile (i).
For profile (ii) a rapid increase in s/suc values is observed at low B/L followed by a
gradual decrease with increasing B/L.
• It can also be seen that s/suc are much higher for profile (ii). Elshazly et al. (2008a)
attribute this to the high modular ratio (Ecol/Esoil). The Young's modulus of the stone
column (Ecol) is similar for both profiles and a higher Ecol/Esoil indicates that the soil
stiffness is lower for profile (ii). This suggests that the lateral support provided by the
surrounding soil has a significant influence on s/suc.
Development of simplified design method and comparison of findings with previous research 193
The design curve developed by Priebe (1995) is adapted to match the numerical analysis by
Elshazly et al. (2008a) (L/d = 10.8/1.06 = 10.2) and is shown in Figure 7.2. As mentioned
previously, Priebe (1995) does not directly relate s/suc to the footing area. However, the
footing width (B) can be determined by assuming the columns are positioned on a square grid
and the distance from the centre-line of the outer row of columns to the edge of the footing is
half the column spacing.
3.0
- Profile (ii); L = 30 m
- Ecol/Esoil = 8.5
2.5
Settlement ratio, s/suc
1.0
0.5
0.0
0 1 2 3 4 5
Normalised footing width, B/L
Figure 7.1- Variation of settlement ratio with normalised footing width for a layered soil and a soft clay, as
determined by Elshazly et al. (2008a)
3.0
- Profile (ii); L = 30 m * K0 refers to soil surrounding the footing
- K0* = 0.55 (not immediately beside the columns)
2.5
Settlement ratio, s/suc
1.0
Priebe (1995)
0.5 - K0 = 1.0
0.0
0 1 2 3 4 5
Normalised footing width, B/L
Figure 7.2 - Comparison of settlement ratios determined by Elshazly et al. (2008a) with Priebe (1995)
While the variation of s/suc with B/L for profile (i) is quite similar to Priebe (1995), it can be
seen in Figure 7.2 that Priebe (1995) predicts lower s/suc. This may be related to the lateral
support provided by the surrounding soil:
Development of simplified design method and comparison of findings with previous research 194
• Elshazly et al. (2008a) account for column installation effects by increasing K0 in the soil
immediately surrounding stone columns. However, the soil surrounding the footing (i.e.
outside the outer row of columns) is modelled as a normally consolidated soil and the
lateral stress is defined by K0 = 1 – sin(φ') (Jaky, 1944). This results in K0 = 0.38–0.44
and 0.55 for profiles (i) and (ii), respectively.
• Priebe (1995) accounts for column installation effects by adopting K0 = 1.0 for the
surrounding soil. A higher lateral stress state in the surrounding soil provides enhanced
support to the columns and, therefore, results in lower s/suc.
Another possible explanation of the high s/suc values determined by Elshazly et al. (2008a)
relative to Priebe (1995) is the flexible loading conditions. Elshazly et al. (2008a) determine
s/suc from the maximum settlement, which occurs near the edge of footings (Figure 7.3). This
is higher than the settlement of rigid footings, assumed by Priebe (1995), and is another
contributing factor to the higher s/suc values.
Figure 7.3 - Deformed mesh of a foundation with B/L = 3 (case of 10.8 m-thick soft clay layer with a 19×19
column arrangement) (Elshazly et al., 2008a)
7.2.4 Comparison of s/suc for end-bearing columns from PLAXIS 3D Foundation &
Priebe (1995)
The design curve developed by Priebe (1995) in Figure 2.24 is adapted to match the
numerical analyses conducted using PLAXIS 3D Foundation (i.e. L/d = 13.9/0.6 = 23.2) and
the variation of s/suc with the number of columns is shown in Figure 7.4. In order to remain
consistent with the assumptions of Priebe (1995), only s/suc values for end-bearing columns
from PLAXIS 3D Foundation are shown in Figure 7.4.
Development of simplified design method and comparison of findings with previous research 195
As expected, it can be seen in Figure 7.4 that an increase in the number of columns leads to
higher s/suc values. An increase in the number of columns results in higher levels of column
confinement and also higher levels of vertical stress with depth (due to larger footing areas).
Therefore, the boundary conditions for groups of columns approach unit cell conditions with
an increasing number of columns and s/suc values tend towards unity.
It can also be seen in Figure 7.4 that s/suc values generally increase with A/AC and this
becomes more significant for larger groups of columns. The design curve developed by
Priebe (1995) agrees favourably with PLAXIS 3D Foundation for A/AC = 3.5, but diverges at
higher area ratios. This confirms that Priebe’s (1995) method is only acceptable for column
groups having low area ratios and is un-conservative for more widely-spaced columns.
0.5
0.4
Settlement ratio, s/suc
0.3
Priebe (1995)
A/Ac = 3.5
0.2
PLAXIS 3D Foundation
0.1 A/Ac=3.5
A/Ac=8.0
A/Ac=14.1
0.0
0 5 10 15 20
Number of columns
Figure 7.4 - Comparison of settlement ratios for groups of end-bearing stone columns from PLAXIS 3D
Foundation with Priebe (1995)
Influence of L/H
For a given B/L, Figure 7.5 shows that s/suc reduces with decreasing L/H. A decrease in L/H
corresponds to a thicker layer of soil beneath the base of columns. The settlement of this
layer is much higher for the infinite grid of columns relative to small groups, due to the
constant stress with depth. Therefore, a decrease in L/H leads to much higher suc and
consequently lower s/suc.
1.0
PLAXIS 3D Foundation
A/Ac = 3.5
A/Ac = 8.0
0.8 A/Ac = 14.1
Settlement ratio, s/suc
Priebe (1995)
A/Ac = 3.5
0.6
L/H = 1.00
0.4
L/H = 0.58
L/H = 0.22
0.2
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
Normalised footing width, B/L
Figure 7.5 - Variation of settlement ratios with normalised footing width (B/L) and column length (L/H)
Influence of A/AC
It appears that A/AC only influences s/suc for end-bearing stone columns (L/H = 1.0). It was
shown in Section 5.3.1 that the settlement performance of an infinite grid of stone columns
improves significantly at low A/AC. This is attributed to the increased levels of lateral
confinement and it can be seen in Figure 5.7(ii) that the influence of lateral confinement
becomes most pronounced for long columns (L > 10 m). Therefore, a significant reduction in
suc occurs for long, end-bearing columns at low A/AC, which leads to higher s/suc.
In contrast to end-bearing columns, it appears that s/suc does not vary with A/AC for floating
stone columns (L/H = 0.22 and 0.58). For a given L/H and B/L, it can be seen in Figure 7.5
that s/suc values are similar for different A/AC. The influence of A/AC upon the settlement
performance of stone columns may be hidden in s/suc. The settlement of floating stone
columns is comprised of settlement from two sections: (i) over the length and (ii) beneath the
base of columns. The settlement of the soil beneath the base of columns is quite high for an
infinite grid of columns, due to the constant stress with depth. This results in high suc which
yields low s/suc and hence hides the influence of A/AC upon s/suc.
Development of simplified design method and comparison of findings with previous research 197
Influence of B/L
A somewhat linear relationship appears to exist between s/suc and B/L for floating stone
columns. However, care must be exercised in assuming a linear relationship for all B/L
values as it can be seen from the distribution of s/suc with B/L for Priebe (1995) that s/suc
approaches unity somewhat asymptotically at high B/L ratios.
The variation of s/suc with B/L for a wide range of column lengths is shown in Figure 7.6. It
can be seen that the influence of A/AC upon s/suc becomes more pronounced as L/H
increases. This may be explained as the settlement of soil beneath the base of columns
reduces at high L/H values. However, it must be noted that for the majority of column
lengths, A/AC does not influence s/suc and a near-linear relationship exists between s/suc and
B/L for a given L/H.
1.0
PLAXIS 3D Foundation
L/H = 0.22
Normalised settlement, s/suc
L/H = 0.36
0.8
L/H = 0.58
Priebe (1995) L/H = 0.79
A/Ac = 3.5 L/H = 1.0
0.6
s/suc = 0.76(B/L)
R2 = 0.78 s/suc = 0.51(B/L)
R2 = 0.97
0.4 s/suc = 0.31(B/L)
R2 = 0.99 s/suc = 0.17(B/L)
R2 = 0.97 s/suc = 0.1(B/L)
R2 = 0.95
0.2
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5
Normalised footing width, B/L
Figure 7.6 - Influence of footing width (B/L) and column length (L/H) upon settlement ratios for groups of
columns
The lines of best fit between s/suc and B/L are shown in Figure 7.6, along with their
corresponding R2 values. R2 is the coefficient of determination and is sometimes referred
Development of simplified design method and comparison of findings with previous research 198
It can be seen that R2 = 1 for this line which indicates that the quadratic equation predicts
the relationship between α and L/H perfectly.
1.0
0.8
α = 0.61(L/H)2 + 0.1(L/H) + 0.06
α = (s/suc) / (B/L)
R2 = 1
0.6
0.4
0.2
0.0
0.0 0.2 0.4 0.6 0.8 1.0 1.2
Normalised column length, L/H
Figure 7.7 - Relationship between α (which relates the s/suc and B/L) with L/H
/_ = £. />
The above equation (7.3) is an interesting and useful finding which allows designers to relate
the settlement of small groups and infinite grids of columns by knowing B/L and L/H.
Development of simplified design method and comparison of findings with previous research 199
It should be noted that equation 7.3 is specific to the benchmark parameters, developed in
Sections 4.4 and 4.6.1. It was shown previously that design parameters such as modular ratios
(Ecol/Esoil) and the column strength (φ) have a significant influence upon the settlement
performance of stone columns. However, it is difficult to predict their influence upon s/suc
values, as the relative influence of Ecol/Esoil and φ upon s and suc may cancel each other.
A short summary of the main characteristics of each mode of deformation is given below:
(i) Punching
- Typically occurs for short columns and small groups of long columns at low A/AC.
- Characterised by low compression ratios and high punching ratios.
- Shear stress and end-bearing pressure develop along the sides and at the base of columns.
(ii) 'Block' failure
- Typically occurs for large groups of columns at low A/AC.
- Characterised by low compression ratios and low punching ratios.
- Columns act as a single entity and displace as a 'block' into the underlying soil, which
ensures low differential settlement at the base of columns and consequently low
punching ratios.
- Shear stress develops along the outer sides of external columns and end-bearing pressure
at the base of columns.
(iii) Bulging
- Typically occurs for long columns at high A/AC.
- Characterised by high compression ratios and low punching ratios.
- Columns bulge at the point of lowest lateral resistance, which generally occurs near the
ground surface where overburden stresses are lowest.
- Bulging precludes the transfer of the applied load to depth which yields low punching
ratios.
• Black (2006) observed that the displacement at the base of single columns becomes
negligible for L/d > 8 and postulates that the mode of deformation changes from
punching to bulging at this point.
• 'Block' failure was observed for small groups of columns whereby the columns acted as a
single entity. The mode of deformation also changes from punching to bulging at L/d =
8, however d is defined as the outer diameter of the group in this instance.
• The findings by Black (2006) are consistent with McKelvey et al. (2004) who conducted
a small scale laboratory study on small groups of stone columns and suggested that a
critical column length exists and lies in the range L/d = 6–10.
5.0
10.0
15.0
While the influence of the stiff crust upon the deformational behaviour of single columns was
not investigated exclusively, it can be seen in Figure 7.9 that compression ratios increase
dramatically with the absence of the stiff crust for a group of 4 columns beneath a 3 m square
footing (i.e. A/AC = 8.0) and agree more closely with Black (2006).
The stiff crust tends to confine columns near the surface, where columns are subject to the
lowest overburden stress and hence are most susceptible to bulging. Therefore, the presence
of the stiff crust enhances the load-carrying capacity of columns and allows a greater
proportion of the load to be transferred to the base of columns. This highlights the advantages
of the FEM, which is capable of capturing the presence of the stiff crust and, therefore,
capable of modelling more realistic soil profiles.
5.0
10.0
Absense of crust
increases column
compression
15.0
However, it appears that the transition between punching and bulging is more gradual at
typical working load levels, as determined by PLAXIS 3D Foundation. It can be seen in
Figures 7.8 and 7.9 that compression ratios are less than 1.0 for all floating columns. This
Development of simplified design method and comparison of findings with previous research 203
indicates that displacement occurs at the base of all floating columns - even in the case of
long columns at high area ratios (A/AC = 8.0 and 14.1), where bulging is established as the
dominant mode of deformation. This finding is quite interesting and suggests that the mode
of deformation at working load levels is a combination of punching and bulging, with one
particular mode more influential for a given A/AC and L.
It was shown previously in Figure 5.15 that the extent of plasticity within columns increases
for fewer supporting columns beneath a 3 m square footing, i.e. as the average load taken per
column increases. Column plasticity reduces the ability of columns to transfer the applied
load to the base of columns and, consequently, results in higher compression ratios. This may
explain the high compression ratios observed by Black (2006) as the deformed shape of
columns is determined at ultimate conditions, which ensured that columns were in an
advanced state of plasticity. In contrast, columns at low working load levels are in more of an
elastic state and can transfer a larger proportion of the load to the base of columns.
Black (2006) used a basalt aggregate to form the stone columns. This was subject to direct
shear and triaxial compression tests which indicated angles of internal friction (φ') of 46° and
43°, respectively. These values are very similar to those adopted in PLAXIS 3D Foundation
(φ' = 45°) and therefore would not contribute towards the discrepancy in the results.
Compression ratios from PLAXIS 3D Foundation are also presented in Figure 7.10. As with
Black (2006), the laboratory tests described by Muir Wood et al. (2000) are conducted in a
homogeneous soil sample and for consistency, the data from PLAXIS 3D Foundation is
specific to soil profiles with and without a stiff crust. The column configurations from
PLAXIS 3D Foundation correspond to groups of 4, 5 and 9 columns beneath a 3 m square
footing, which correspond to A/AC = 8.0, 6.4 and 3.5, respectively. In the interest of
consistency, only compression ratios for 3 m long columns (i.e. L/B = 1.0) from PLAXIS 3D
Foundation are presented in Figure 7.10. The position of columns for PLAXIS 3D
Foundation is defined as the normalised radial distance of columns from the footing centre
(r/r0), where r0 is the equivalent footing radius (r0 = B/√π).
1.0
Figure 7.10 - Influence of column location upon compression ratios measured by Muir Wood et al. (2000) and
determined PLAXIS 3D Foundation
It can be seen in Figure 7.10 that the findings from PLAXIS 3D Foundation are consistent
with Muir Wood et al. (2000), as higher compression ratios are observed with increasing
A/AC and for columns positioned towards the footing edge.
While the absence of a stiff crust leads to an increase in compression ratios from PLAXIS
3D Foundation, it can be seen that Muir Wood et al. (2000) observed higher compression
ratios for all A/AC. This may be explained as follows:
• The angle of shearing resistance (φ' ≈ 30°) for the sand columns in the laboratory tests is
much lower than PLAXIS 3D Foundation (φ' = 45°). It was shown in Section 5.4.7 that
decreasing φ from 50° to 40° had a significant influence on compression ratios, with an
increase of compression ratios in the order of 0.05–0.15 observed for 3 m long columns.
Development of simplified design method and comparison of findings with previous research 206
• The deformed shape of columns was determined at ultimate conditions. As with Black
(2006), loading was terminated at high strain levels which ensured the columns were in
an advanced state of plasticity and hence yield high compression ratios.
However, these findings are in contrast to McKelvey et al. (2004) and Black (2006) who
suggest that a critical length may exist in the range L/d = 6–10, and is not dependent on the
footing width. In light of these more recent findings, the data presented by Muir Wood et al.
(2000) was re-examined in this thesis plotted in Figure 7.11.
Muir Wood et al. (2000) observed punching and bulging for different configurations of
columns and these are plotted against L/D in Figure 7.11(i). The mode of deformation
changes from punching to bulging in the range L/D = 1.0–1.5, which prompted the authors to
postulate that the mode of deformation is controlled by the footing diameter. However, the
observed modes of deformation are plotted against L/d in Figure 7.11(ii) and it is clear that,
with the exception of two tests out of total of 14, that the mode of deformation changes from
punching to bulging in the range L/d = 9.1–13.6. This is consistent with more recent findings.
The results from this thesis indicate the area ratio and column length, rather than the footing
width, dictate the mode of deformation. It was established in Chapter 5 that different
configurations of area ratio and column length yield either a punching-type (i.e. punching of
‘block’ failure) or a bulging mode of deformation. For example, it was shown in Figure
6.3(b-i) that a 2×2 group of 8 m long columns at A/AC = 3.5 (i.e. B = 2 m; L/B = 4) exhibit a
punching mode of deformation, while a 4×4 group of 8 m long columns at A/AC = 14.1 (i.e.
B = 8 m; L/B = 1) exhibit a bulging mode of deformation (Figure 6.5(b-iii)).
Development of simplified design method and comparison of findings with previous research 207
0.0 0
D=100 mm D=100 mm
0.2 2
d d
0.4 4
0.6 d = 11 mm Punching d = 11 mm
6 d = 17.5 mm
Punching d = 17.5 mm
0.8
8
1.0
10 Possible
1.2 Postulated critical length
critical length 12 at L/d = 9.1
1.4
1.6 14 Bulging
1.8 Bulging 16
L/D L/d
(i) (ii)
Figure 7.11 - Observations of punching and bulging for groups of columns, reported by Muir Wood et al.
(2000), plotted against (i) L/D and (ii) L/d
1.0 0.20
1 col 1 col
A/Ac = 3.5 A/Ac = 3.5
0.8 A/Ac = 8.0 A/Ac = 8.0
Compression ratio
0.15 L=3m
Punching ratio
11 m A/Ac = 14.1 A/Ac = 14.1
5m
0.6
8m 0.10
8m
0.4
5m 11 m
0.05
0.2
L=3m
0.0 0.00
1.0 1.5 2.0 2.5 3.0 3.5 1.0 1.5 2.0 2.5 3.0 3.5
Settlement improvement factor, n Settlement improvement factor, n
(a-i) (a-ii)
1.0 0.20
2×2 2×2
A/Ac = 3.5 A/Ac = 3.5
0.8 A/Ac = 8.0 A/Ac = 8.0
Compression ratio
0.15
Punching ratio
A/Ac = 14.1 A/Ac = 14.1
0.6 11 m L=3m
0.10 5m
0.4 8m
8m
11 m
0.05
0.2
5m
L=3m
0.0 0.00
1.0 1.5 2.0 2.5 3.0 3.5 1.0 1.5 2.0 2.5 3.0 3.5
Settlement improvement factor, n Settlement improvement factor, n
(b-i) (b-ii)
1.0 0.20
3×3 3×3
A/Ac = 3.5 A/Ac = 3.5
0.8 A/Ac = 8.0 A/Ac = 8.0
Compression ratio
0.15
Punching ratio
1.0 0.20
4×4 4×4
A/Ac = 3.5 A/Ac = 3.5
0.8 A/Ac = 8.0 A/Ac = 8.0
Compression ratio
0.15
Punching ratio
Figure 7.12 - Variation of (i) compression ratio and (ii) punching ratio with settlement improvement factor for
(a) 1 column, (b) 2×2, (c) 3×3 and (d) 4×4 groups
Development of simplified design method and comparison of findings with previous research 209
between area ratio and column length is an inter-dependent one, as different configurations of
columns result in different modes of deformation. Each mode of deformation has its own
characteristics which makes its advantageous in different soil profiles.
It was shown in Figure 7.12 that the settlement performance of columns which exhibit
punching, and especially 'block' failure, are most sensitive to changes in compression ratios.
This is consistent with previous findings that the influence of column stiffness upon the
settlement performance of small groups of columns is most pronounced for these modes of
deformation. It was also observed that columns which exhibit punching and 'block' failure are
most sensitive to changes in punching ratios. This reflects the load-transfer mechanism, as a
large proportion of the applied load is transferred to the base of columns. Therefore these
modes of deformation are not highly dependent upon the properties of the surrounding soil
and are most efficient is soil profiles with weak upper layers.
In contrast, it can be seen that influence of compression and punching ratios upon n values
diminishes for the bulging mode of deformation. Column bulging is highly dependent upon
the lateral support provided by the surrounding soil and generally occurs near the ground
surface, where overburden stresses are lowest. The bulging mode of deformation makes
efficient use of the high stiffness properties of the stiff crust. Columns are confined in the
upper layers which tends to push bulging deeper and hence increase the load-carrying
capacity of columns.
Charles & Watts (1983) are also shown. It can be seen in Figure 7.13 that stress concentration
ratios predicted by PLAXIS 3D Foundation agree quite well with field and laboratory data.
15
Plaxis 3D Foundation
7
Barksdale & Bachus (1983)
5
1
2 4 6 8 10 12 14 16
Area ratio, A/AC
Figure 7.13 - Comparison of stress concentration ratios for an infinite grid of end
end-bearing
bearing columns with field
data
It can also be seen in Figure 7.14 that Priebe (1995) generally predicts
dicts lower stress
concentration ratios than Pulko & Majes (2005),
(2005) especially at low area ratios. This is due to
the different assumptions inherent in each design method. Priebe (1995) assumes columns are
in an active state throughout loading, while Pulko & Majes (2005) account for the elastic
behaviour of columns at low load levels. Stress concentration ratios
ratio are much higher for
columns in an elastic state and consequently Pulko & Majes (2005) predict higher stress
Development of simplified design method and comparison of findings with previous research 212
concentration ratios than Priebe (1995). This is most pronounced at low area ratios as
columns are better confined and hence carry a larger proportion of the applied load when in
an elastic state.
Perhaps the most important finding from Figure 7.14 is that both Priebe (1995) and Pulko &
Majes (2005) predict higher surface stress concentrations ratios than PLAXIS 3D Foundation.
This is quite interesting given that both of these design methods predict similar settlement
improvement factors to PLAXIS 3D Foundation (see Section 4.6.4). Considering the equation
of equilibrium developed by Aboshi et al. (1979) in equation (2.27), it seems intuitive that
higher stress concentration ratios at the surface would yield higher settlement improvement
factors. However, this is not the case and it appears that the distribution of stress
concentration ratios with depth, rather than at the surface, may have more of an influence
upon the settlement performance of stone columns and thus merits further investigation.
15
Plaxis 3D Foundation
Stress concentration ratio, σcol/σsoil
13 Priebe (1995)
Pulko & Majes (2005)
11 Balaam & Booker (1981)
1
2 4 6 8 10 12 14 16
Area ratio, A/Ac
Figure 7.14 - Comparison of stress concentration ratios at the surface for an infinite grid of end-bearing
columns with analytical design methods
corresponds to the top of the lower Carse clay layer (i.e. the weakest part of the soil profile).
Stress concentration ratios are lower in the overlying soil, even though these soil layers are in
an elastic state. This may be explained as an infinite grid of columns do not develop shear
stress along their sides and hence the maximum vertical load which columns can carry at the
surface is limited by the vertical stress where bulging occurs. The extent of plasticity within
columns at higher area ratios (i.e. A/AC = 8.0 and 14.1) is more significant and despite the
large scatter, a constant distribution of stress concentration ratios with depth is observed in
Figure 7.15 for yielded sections of columns. It can also be seen that stress concentration
ratios reduce with depth in the lower sections of columns which are in an elastic state.
2 2
4 4
Depth, z (m)
Depth, z (m)
6 6
8 8
PLAXIS 3D Foundation PLAXIS 3D Foundation
A/Ac = 3.5 A/Ac = 3.5
10 A/Ac = 8.0 10 A/Ac = 8.0
A/Ac = 14.1 A/Ac = 14.1
Priebe (1995) Pulko & Majes (2005)
12 A/Ac = 3.5 12 A/Ac = 3.5
A/Ac = 8.0 A/Ac = 8.0
A/Ac = 14.1 A/Ac = 14.1
14 14
(i) (ii)
Figure 7.15 - Comparison of stress concentration ratios with depth from PLAXIS 3D Foundation with (i) Priebe
(1995) and (ii) Pulko & Majes (2005) for an infinite grid of end-bearing columns
It can be seen in Figure 7.15(i) that Priebe (1995) predicts lower stress concentration ratios
than PLAXIS 3D Foundation within yielded sections of columns. This is most pronounced
for closely-spaced columns (A/AC = 3.5) at z = 1.9 m and may be explained as Priebe (1995)
fails to account for stress concentration ratios developed when columns are in an elastic state
(i.e. at low load levels). However, stress concentration ratios for closely-spaced columns
reduce at a faster rate with depth for PLAXIS 3D Foundation, as the lower sections of
columns are in an elastic state. It can also be seen in Figure 7.15(i) that stress concentration
ratios at higher area ratios (i.e. A/C = 8.0 and 14.1) are relatively constant with depth in
yielded sections of columns compared to Priebe (1995). This is again attributed to the
assumption by Priebe (1995) that columns are in an active state throughout loading.
It is interesting to note that Priebe (1995) generally under-predicts stress concentration ratios
but still yields similar settlement improvement factors to PLAXIS 3D Foundation. This may
explained as stress concentration ratios determined from Priebe (1995) do not contain the
modification factor for overburden stress, which is taken into account when determining
settlement improvement factors (n2).
It also appears that Pulko & Majes (2005) slightly over-estimate stress concentration ratios
compared to PLAXIS 3D Foundation and this may be due to the high modular ratios adopted
for the analytical design methods. The soil profile adopted for the analytical design methods
consists of three layers, each defined by a single Young's modulus. In contrast, PLAXIS 3D
Development of simplified design method and comparison of findings with previous research 215
Foundation accounts for the stress dependency of soil stiffness and captures the increase of
soil stiffness with depth. The modular ratios adopted for each layer in Pulko & Majes (2005)
are higher than PLAXIS 3D Foundation which increases the magnitude of stress
concentration ratios and also extent of plasticity within columns.
The influence of higher stress concentration ratios upon the settlement performance of stone
columns is evident in Figure 4.16, as Pulko & Majes (2005) predict higher settlement
improvement factors than PLAXIS 3D Foundation at low area ratios (A/AC = 3.5). This may
be explained as the majority of columns are in an elastic state at low area ratios and are
therefore more susceptible to differences in modular ratios. However, influence of the higher
stress concentration ratios is negligible at higher area ratios, as Pulko & Majes (2005) predict
similar settlement improvement factors to PLAXIS 3D Foundation (Figure 4.16). The extent
of plasticity within columns is more pronounced at higher area ratio and hence the influence
of modular ratios upon the settlement performance of columns reduces.
The distribution of stress concentration ratios with depth clearly indicates when columns are
in an elastic and plastic stress state. This is quite useful for stability analyses where a
rotational failure surface may intersect columns at different depths. However, it difficult to
form a direct relationship between stress concentration ratio and settlement improvement
factors for the following reasons:
(i) Settlement improvement factors are based on the increment of vertical stress (∆σcol and
∆σsoil) with depth rather than the total vertical stress(σcol/σsoil), which is used to define
stress concentration ratios, i.e. σcol/σsoil == (∆σcol + γ'col.z)/(∆σsoil + γ'soil.z). Furthermore,
the increment of vertical stress reduces with depth beneath small loaded area.
(ii) The distribution of stress concentration ratios only shows the stress state at the end of
loading and does not capture the proportion of the load applied to columns in an elastic
and plastic state. The majority of footing settlement is comprised from plastic strain and
it is necessary to predict the yield load.
Chapter 8
Conclusions
8.1 Introduction
A series of three-dimensional FEA were undertaken to examine the influence of key design
parameters upon the settlement performance and deformational behaviour of small groups of
stone columns. PLAXIS 3D Foundation is adopted for this research in conjunction with the
advanced elastic-plastic Hardening Soil model. A set of material parameters and soil profile
were developed for the well characterised Bothkennar test site, which formed the basis of this
research. The key design parameters investigated were area ratio, column length,
confinement, stiffness, strength, installation effects and the presence of a stiff crust. The main
conclusions are summarised in this chapter.
solution than Priebe (1995), and was shown to predict a similar settlement response to
PLAXIS 3D Foundation.
columns and the surrounding soil. The main characteristics of each mode of deformation are
outlined below:
(i) Punching:
- Characterised by high punching and low compression ratios.
- Typically occurs for short columns at all area ratios and small groups of columns at low
area ratios.
- High shear strains develop along the sides and beneath the base of columns which
exhibit punching.
(ii) 'Block' failure:
- Characterised by low punching and high compression ratios.
- Most prominent in large groups of columns at low area ratios.
- Columns act as a single entity and punch uniformly into the underlying soil, thus
resulting in low punching ratios. 'Block' failure is an extension of the punching mode of
deformation, as both modes of deformation develop similar load transfer mechanisms.
- High shear strains develop along the external sides of outer columns and beneath the
base of columns. No shear strains develop in the central zone of soil bounded by
columns, which confirms that columns punch as a 'block' with the surrounding soil.
(iii) Bulging:
- Characterised by low punching and high compression ratios.
- Prevalent in long columns at high area ratios.
- This mode of deformation is highly dependent upon the lateral support provided by the
surrounding soil and as such, bulging generally occurs near the ground surface where
overburden stresses are lowest.
- High shear strains develop within columns near the top of the lower Carse clay layer,
i.e. the weakest part of the soil profile.
In addition to the three modes of deformation, the development of shear planes which extend
from the corners of pad footings to the point of column bulging beneath the centre of footings
were also observed. It should be noted that the shear planes are a secondary mode of
deformation; however they may become more prominent at higher load levels, as observed by
Muir Wood et al. (2000).
The deformational behaviour of small groups of stone columns is governed by the area ratio
and column length. While the number of columns enhances the settlement performance of
Conclusions 219
stone columns for a given area ratio, it does not fundamentally change the mode of
deformation. This is an important findings which was used in the development of a simplified
design method.
The increase in settlement improvement factors with column length is most pronounced for
columns at low area ratios, for which punching is the dominant mode of deformation. Shear
stress develops along the sides of columns and increasing the column length directly
improves the settlement performance of stone columns. A certain degree of punching is
evident for all the configurations of columns analysed, which is consistent with the finding
that settlement improvement factors increase with column length for all configurations of
columns. Therefore, the existence of a unique critical column length, proposed from previous
laboratory studies, is not supported by this numerical modelling.
The presence of the stiff crust also contributes towards a combination of punching and
bulging modes of deformation. Previous research from laboratory tests observed well-defined
crossovers between punching and bulging modes of deformation at column lengths in the
range L/d = 6–10. However, this has been shown to be in part due to the absence of the stiff
crust in the laboratory tests (i.e. homogeneous soil samples). Columns in a homogeneous soil
Conclusions 220
sample are more susceptible to bulging near the ground surface, which reduces the ability of
columns to transfer the applied load to the base of columns and hence reduces punching
potential.
Stress concentration ratios from PLAXIS 3D Foundation are in the typical range (2.5–5.0)
proposed by Barksdale & Bachus (1983).
In contrast to an infinite grid of columns, stress concentration ratios for small groups of stone
columns are dependent on the mode of deformation. An increase in stress concentration ratios
is observed with increasing column length for columns which are punching. However, no
increase in stress concentration ratios is observed with increasing column length for bulging.
Settlement ratios for end-bearing columns increase with area ratio and were found to agree
favourably with Priebe’s (1995) s/suc values at low area ratios. Previous research is only
focussed on end-bearing columns at low area ratios. However, this thesis examines the
influence of both area ratio and floating columns (i.e. partial-depth treatment) upon s/suc for
the first time.
It was observed that area ratio only influences s/suc values for end-bearing stone columns.
Moreover, for a given normalised column length (L/H), a near-linear relationship exists
between s/suc and the normalised footing width (B/L). Therefore, it appears that the footing
width, and hence the stress distribution in the underlying soil, governs s/suc values. A
simplified design method was developed which allows the settlement of a small groups of
columns to be quickly estimated from B/L, L/H and suc. It is recommended that Pulko &
Majes (2005) be used as an estimate of suc.
settlement at various depths, using extensometers, for unreinforced and reinforced footings.
There is also further potential to extend the design method to other soil profiles:
- It was shown that the modular ratio and stress state in the surrounding soil both play an
important role in the settlement performance of stone columns. It would be valuable to
examine their influence upon settlement ratios for floating stone columns.
- Although the range of B/L and L/H investigated in this thesis covers a practicable range,
investigating the settlement performance of small groups of stone columns in a soil
profile of varying thickness would allow an even greater range of B/L and L/H to be
examined.
This thesis has focussed on the long term settlement performance of small groups of columns.
In addition to improving the settlement performance, stone columns also act as vertical drains
and enhance the drainage capacity of soils. It is well established that the soil carries a large
proportion of the vertical stress at the start of loading due to the high undrained stiffness.
However, vertical stress is transferred from the soil to the stone columns as consolidation
proceeds. The dissipation of excess pore pressure for small groups of stone columns has not
been investigated to date and it would also be useful to examine how the position of columns
relative to the footing edge affects the stress transfer during consolidation.
The installation effects associated with vibro stone columns have been touched on in this
thesis. Installing stone columns imposes a complex stress regime on the surrounding soil; the
displacing effects of the poker increase the soil stiffness, while the horizontal vibrations from
the poker can lead to remoulding and consequently a loss of soil strength. More field
measurements are essential to examine the influence of different soil types and column
configurations upon column installation effects.
224
References
Aboshi, H., Ichimoto, I., Enoki, M. and Harada, K. (1979). The compozer - a method to improve
characteristics of soft clays by inclusion of large diameter sand columns. Proc. Int. Conf. on Soil
Reinforcement: Reinforced Earth and Other Techniques, Paris, Vol. 1, 211–216
Allman, M. A. and Atkinson, J. H. (1992). Mechanical properties of reconstituted Bothkennar soil.
Géotechnique, Vol. 42, No. 2, 289–301
Ambily, A. P. and Gandhi, S. R. (2007). Behaviour of stone columns based on experimental and FEM
analysis. Journal of Geotechnical and Geoenvironmental Engineering, Vol. 133, No. 4, 405–415
Andreou, P. and Papadopoulos, V. (2006). Modelling stone columns in soft clay. Proc. 6th European
Conf. on Numerical Methods in Geotechnical Engineering, Graz, 777–780
Atkinson, J. H. and Sallfors, G. (1991). Experimental determination of soil properties. Proc. 10th
European Conf. on Soil Mechanics and Foundation Engineering, Florence, Vol. 3, 915–956
Barksdale, R. D. and Bachus, R. C. (1983) Design and construction of stone columns, Report no.
FHWA/RD-83/026, National Information Service, Springfield, Virginia, USA
Balaam, N. P. and Booker, J. R. (1981). Analysis of rigid rafts supported by granular piles. Int.
Journal for Numerical and Analytical Methods in Geomechanics, Vol. 5, 379–403
Balaam, N. P. and Booker, J. R. (1985). Effect of stone column yield on settlement of rigid
foundations in stabilised clay. Int. Journal for Numerical and Analytical Methods in
Geomechanics, Vol. 9, 331–351
Balaam, N. P., Poulos, H. G. and Brown, P. T. (1977). Settlement analysis of soft clays reinforced
with granular piles. Proc. 5th Southeast Asian Conf. on Soil Engineering, Bangkok, Vol. 1, 81–91
Baligh, M. M. (1985). Strain path method. Journal of the Geotechnical Engineering Division, ASCE,
Vol. 111, No. 9, 1108–1136
Biarez, J., Gambin, M., Gomes-Corriea, A., Falvigny, E. and Branque, D. (1998). Using
pressuremeter to obtain parameters of elastoplastic models for sands. Proc. 1st Int. Conf. on Site
Characterization, Atlanta, Georgia, USA, 747–752
Biot, M. A. (1941). General theory of three-dimensional consolidation. Journal of Applied Sciences,
Vol. 12, 155–164
Black, J. A. (2006). The settlement performance of a footing supported on soft clay reinforced with
vibrated stone columns. PhD thesis, Queen's University of Belfast
Bolton, M. D. (1986). The strength and dilatancy of sands. Géotechnique, Vol. 36, No. 1, 65–78.
Boussinesq, J. (1885). Applications des Potentiels à L’étude de L’équilibre et du Mouvement des
Solides Élastiques, Gauthier-Villars, Paris
References 225
Brauns, J. (1978). Die anfangstraglast von schottersäulen im bindigen untergrund. Die Bautechnik,
263–270
Brinkgreve, R. B. J and Broere, W. (2006). PLAXIS 3D Foundation Manual Version 2. PLAXIS BV
Castro, J. and Sagaseta, C. (2009). Consolidation around stone columns. Influence of column
deformation. Int. Journal for Numerical and Analytical Methods in Geomechanics, Vol. 33, Issue
7, 851–877
Castro, J. (2007). Pore pressures during stone column installation. Proc. 18th European Young
Geotechnical Engineers' Conference, Ancona
Charles, J. A. and Watts, K. S. (1983). Compressibility of soft clay reinforced with granular columns.
Proc. 8th European Conf. on Soil Mechanics and Foundation Engineering, Helsinki, Vol. 1,
347–352
Cunze, G. (1985). Ein beitrag zur abschätzung des porenwasserüberdruckes beim rammen von
verdrängungspfählen in bindigen böden. Dissertation University of Hannover
Debats, J. M., Guetif, Z. and Bouassida, M. (2003). Soft soil improvement due to vibro-compacted
columns installation. Proc. Int. Workshop on Geotechnics of Soft Soils - Theory and Practice,
Noordwijkerhout, Netherlands, 551–556
Domingues, T. S., Borges, J. L. and Cardoso, A. S. (2007a). Parametric study of stone columns in
embankments on soft soils by finite element method. Applications of Computational Mechanics
in Geotechnical Engineering V: Proc. 5th Int. Workshop, Guimaraes, 281–291
Domingues, T. S., Borges, J. L. and Cardoso, A. S. (2007b). Stone columns in embankments on soft
soils. Analysis of the effect of the gravel deformability. 14th European Conf. on Soil Mechanics
and Geotechnical Engineering, Madrid, Vol. 3, 1445–1450
Duncan, J. M. and Chang, C. Y. (1970). Nonlinear analysis of stress and strain in soil. ASCE Journal
of the Soil Mechanics and Foundation Division, Vol. 96, No. 5, 1629–1653
Egan, D., Scott, W. and McCabe, B. A. (2008). Installation effects of vibro replacement stone
columns in soft clay. Proc. 2nd Int. Workshop on the Geotechnics of Soft Soils, Glasgow, 23–30
Elshazly, H. A., Hafez, D. and Mosaad, M. (2006). Back calculating vibro-installation stresses in
stone columns reinforced ground. Proc. of the ICE - Ground Improvement, Vol. 26, Issue 2, 47–
53
Elshazly, H. A., Hafez, D. H. and Mossaad, M. E. (2008a). Reliability of conventional settlement
evaluation for circular foundations on stone columns. Journal of Geotechnical and Geological
Engineering, Vol. 26, No. 3, 323–334
Elshazly, H., Elkasabgy, M. and Elleboudy, A. (2008b). Effect of inter-column spacing on soil
stresses due to vibro-installed stone columns: Interesting findings. Journal of Geotechnical and
Geological Engineering, Vol. 26, No. 2, 225–236
Gäb, M., Schweiger, H. F., Thurner, R. and Adam, D. (2007). Field trial to investigate the
performance of floating stone columns. Proc. 14th European Conf. on Soil Mechanics and
Geotechnical Engineering, Madrid, 1311–1316
Gäb, M., Schweiger, H. F., Kamrat-Pietraszewska, D. and Karstunen, M. (2008). Numerical analysis
of a floating stone column foundation using different constitutive models. Proc. 2nd Int.
Workshop on Geotechnics of Soft Soils, Glasgow, 137–142
References 226
Institution of Civil Engineers. (1992). Bothkennar soft clay test site: characterization and lessons
learned. Géotechnique, Vol. 42, No. 2, 161–378.
Gibson, R. E. and Anderson, W. F. (1961). In-situ measurements of soil properties with the
pressuremeter. Civil Engineering and Public Works Review, Vol. 56, No. 658, 6l5–6l8
Gill, D.R. and Lehane, B.M. (2001). An optical technique for investigating soil displacement patterns.
Geotechnical Testing Journal, Vol. 24, 324–329
Goughnour, R. R. and Bayuk, A. A. (1979). Analysis of stone column-soil matrix interaction under
vertical load. Proc. Int. Conf. on Soil Mechanics Reinforcement, Paris, 271–277
Greenwood, D. A. (1970). Mechanical improvement of soils below ground surface. Geotechnical
Engineering, Proc. of Institution of Civil Engineers, London, 9–20
Greenwood, D. A. (1991). Load tests on stone columns. In Deep Foundation Improvements: Design,
Construction and Testing, Esrig M. I. and Bachus R. C. (eds.), American Society for Testing and
Materials, Philadelphia, ASTM STP 1089, 141–171
Guetif, Z., Bouassida, M. and Debats, J.M. (2007). Improved soft clay characteristics due to stone
column installation. Computers and Geotechnics, Vol. 34, Issue 2, 104–111
Han, J. and Ye, S. L. (2001). Simplified method for consolidation rate of stone column reinforced
foundations. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, Vol. 127, No.
7, 597–603
Hight, D. W., Bond, A. J., and Legge, J. D. (1992). Characterization of the Bothkennar clay: an
overview. Géotechnique, Vol. 42, No. 2, 303–347
Hu, W. (1995). Physical modelling of group behaviour of stone column foundations. PhD thesis,
University of Glasgow
Hughes, J. M. and Withers, N. J. (1974). Reinforcing of soft cohesive soils with stone columns.
Ground Engineering, Vol. 7, No. 3, 42–49
Jaky, J. (1944). The coefficient of earth pressure at rest. Journal of the Union of Hungarian Engineers
and Architects, 355–358 (in Hungarian).
Janbu, J. (1963). Soil compressibility as determined by oedometer and triaxial tests. Proc. European
Conf. on Soil Mechanics and Foundation Engineering, Wiesbaden, Vol. 1, 19–25
Jardine, R. J., Lehane, B. M., Smith, P. R. and Gildea, P. A. (1995). Vertical loading experiments on
rigid pad foundations at Bothkennar. Géotechnique, Vol. 45, No. 4, 573–597
Kirsch, F. and Sondermann, W. (2003) Field measurements and numerical analysis of the stress
distribution below stone column supported embankments and their stability. Int. Workshop on
Geotechnics of Soft Soils - Theory and Practice, Noordwijkerhout, Netherlands, 595–600
Kirsch, F. (2006). Vibro stone column installation and its effect on the ground improvement. Int.
Conf. on Numerical Simulation of Construction Processes in Geotechnical Engineering for
Urban Environment, Bochum, 115–124
Kirsch, F. (2008). Evaluation of ground improvement by groups of vibro stone columns using field
measurements and numerical analysis. Proc. of the 2nd Int. Workshop on the Geotechnics of Soft
Soils, Glasgow, 241–248
Lambe, T. W. and Whitman, R. V. (1979). Soil mechanics (SI Version). John Wiley and Sons, Inc.
References 227
Leroueil, S., Lerat, P., Hight, D. W. and Powell, J. J. M. (1992). Hydraulic conductivity of a recent
estuarine silty clay at Bothkennar. Géotechnique, Vol. 42, No. 2, 275–288
Mar, A. (2002). How to undertake finite element based geotechnical analysis. NAFEMS Publication
McCabe, B. A. (2002). Experimental investigations of driven pile group behaviour in Belfast Soft
Clay. PhD Thesis, Department of Civil Engineering, University of Dublin, Trinity College
McCabe, B. A., McNeill, J. A. and Black, J. A. (2007). Ground improvement using the vibro-stone
column technique. Transactions of the Institution of Engineers of Ireland, Galway
McCabe, B. A., Killeen, M. M. and Egan, D. (2008). Challenges faced in 3-d finite element modeling
of stone column construction. Proc. Joint Symposium on Concrete Research and Bridge and
Infrastructure Research in Ireland, NUI Galway, 393–400
McCabe, B. A., Nimmons, G. J. and Egan, D. (2009). A review of field performance of stone columns
in soft soils. Proc. Institution of Civil Engineers - Geotechnical Engineering, Vol. 162, No. 6,
323–334
McKelvey, D., V. Sivakumar, Bell, A.L. and Graham, J. (2004). Modelling vibrated stone columns in
soft clay. Proc. Institution of Civil Engineers - Geotechnical Engineering, Vol. 157, No. 3, 137–
149
Mitchell, J.K. and Huber, T. R. (1985). Performance of a stone column foundation. Journal
Geotechnical Engineering, ASCE, Vol. 11, No. 2, 205–223
Muir Wood, D., Hu, W. and Nash, D. F. T. (2000). Group effects in stone column foundations: model
tests. Geotechnique 50(6): 689–698.
Munfakh, G. A., Sarkar, S. K. and Castelli, R. J. (1984). Performance of a test embankment founded
on stone columns. Proc. Int. Conf. on Advances in Piling and Ground Testing, London, 259–265
Narasimha Rao, S., Madhiyan, M. and Prasad, Y. V. S. N. (1992). Influence of bearing area on the
behaviour of stone columns. Indian Geotechnical Conf., Calcutta, 235–237
Nash, D. F. T., Powell, J. J. M. and Lloyd, I. M. (1992a). Initial investigations of the soft clay test site
at Bothkennar. Géotechnique, Vol. 42, No. 2, 163–181
Nash, D. F. T., Sills, G. C. and Davison, L. R. (1992b). One-dimensional consolidation testing of the
soft clay from Bothkennar. Géotechnique, Vol. 42, No. 2, 241–256
Paul, M. A., Peacock, J. D. and Wood, B. F. (1992). The engineering geology of the Carse clay of the
national soft clay research site, Bothkennar. Géotechnique, Vol. 42, No. 2, 183–198
Potts, D.M. and Zdravkovic, L. (1999). Finite element analysis in geotechnical engineering: theory.
Thomas Telford
Potts, D.M. and Zdravkovic, L. (2001). Finite element analysis in geotechnical engineering:
application. Thomas Telford
Powrie, W. (2004). Soil mechanics: concepts and applications. Spon Press
Priebe, H.J. (1976). Evaluation of the settlement reduction of a foundation improved by vibro-
replacement. Bautechnik, No. 5, 160–162 (in German)
Priebe, H. J. (1995). The design of vibro replacement. Ground Engineering, Vol. 28, No. 10, 31–37
References 228
Pulko, B. and Majes, B. (2005). Simple and accurate prediction of settlements of stone column
reinforced soil. 16th International Conference on Soil Mechanics and Geotechnical Engineering,
Osaka, Vol. 3, 1401–1404
Randolph, M. F., Carter, J. P. and Wroth, C. P. (1979a). Driven piles in clay - the effects of
installation and subsequent consolidation. Géotechnique, Vol. 29, No. 4, 361–393
Randolph, M. F., Steenfelt, J. S. and Wroth C. P. (1979b). The effect of pile type on design
parameters for driven piles. Proc. 7th European Conf. on Soil Mechanics and Foundation
Engineering, London, 2, 107–114
Rowe, P. W. (1962). The stress-dilatancy relation for static equilibrium of an assembly of particles in
contact. Proc. Royal Society A, Vol. 269, No. 1339, 500–527
Schanz, T., (1998). Zur modellierung des mechanischen verhaltens von reibungsmaterialen.
Habilitation, Stuttgart Universität
Schanz, T., Vermeer, P.A., Bonnier, P.G., (1999). The hardening-soil model: Formulation and
verification. Beyond 2000 in Computational Geotechnics - 10 Years of PLAXIS International, R.
B. J. Brinkgreve (ed.), A. A. Balkema, Rotterdam, Netherlands, 281–290
Soderberg, L. O. (1962). Consolidation theory applied to foundation pile time effects. Géotechnique,
Vol. 12, No. 3, 217–225
Sondermann, W. and Wehr, W. (2004). Deep vibro techniques, Ground Improvement, 2nd Edition,
edited by M.P. Moseley and K. Kirsch, 57–92, Spon Press
Terzaghi, K. (1925). Erdbaumechanik auf bodenphysikalischer grundlage, F. Deuticke, Vienna
Thornburn, S. and MacVicar, R. S. L. (1968). Soil stabilization employing surface and depth
vibrators. The Structural Engineer, Vol. 46, No. 10, 309–316
Wehr, J. and Herle, I. (2006). Exercise on calculation of stone columns - Priebe method and FEM.
Proc. 6th European Conf. on Numerical Methods in Geotechnical Engineering, Graz, 773–776.
Wehr, W. (2004). Stone columns - single columns and group behavior. Proc. 5th Int. Conf. on Ground
Improvement Techniques, Malaysia, 329–340
Wehr, J. (2006a). The undrained cohesion of the soil as a criterion for the column installation with a
depth vibrator. TRANSVIB 2006, Gonin, Holeyman and Rocher-Lacoste (ed.), Editions du LCPC.
Paris
Wehr, W. C. S. (2006b). Stone columns - group behaviour and influence of footing flexibility. Proc.
6th European Conf. on Numerical Methods in Geotechnical Engineering, Graz, 767–772
Witt, K. J. (1978). Versagensmechanismus einzeln belasteter schottersäulen im bindigin untergrund
bei plötzlicher belastung. B. Sc. thesis, University of Karlsruhe
Woo-Seok Bae, Bang-Woong Shin, Byung-Chul An and Joo-Sub Kim (2002). Behaviours of
foundation system improved with stone columns. Proc. 12th Int. Offshore and Polar Engineering
Conf., Kitakyushu, Japan, 675–678
Yu, H. (2000). Cavity expansion methods in geomechanics. Kluwer Academic Publishers BV,
Dordrecht, NL