Ashour 16 PH D
Ashour 16 PH D
Ashour 16 PH D
By
SAMIR ASHOUR
A Thesis submitted to
DOCTOR OF PHILOSOPHY
ABSTRACT
Soft clay soils in railway track can be problematic as, unless they are treated, they can result
in increased deformation of track. This will inheritably mean lower track speed. A number of
techniques including stone columns are available for improving strength of weak soils.
The use of stone columns, for improving both bearing capacity and settlement is well
rehearsed for static loading; little is understood about their response, when subjected to cyclic
loading. This study is focused on investigating the behaviour of stone columns when
A series of monotonic and cyclic loading conditions were undertaken on two laboratory
models undrained triaxial (diameter 100 mm, and 200 mm height) and large scale model
(diameter 300 mm, and 300 mm height). Tests were conducted on both soft soils (no column)
and soil/ stone column composite. All tests were performed on normally consolidated
specimens of soft clay (undrained shear strength of ≈ 12 kPa) and for treated soils they were
reinforced with 28 mm diameter stone columns. Three cyclic stress (50, 60 and 70 kPa) on
subgrade level and three loading frequencies (0.5, 1 and 3 Hz) simulating different train
speeds (35, 70 and 225 km/hr) were used to study the performance of both soft soils (with and
without column). The effect of both cyclic stresses and loading frequency on the permanent
deformation, soils stiffness and pore water pressure generation were investigated.
Generally, soft clay bed reinforced with stone column showed a significant improvement in
approximately 30% increase in failure load of soil with the stone column compared to that soil
only.
i
Samir Ashour Abstract
It was found that threshold dynamic stress of soil for cyclic loading increased from 50 kPa for
soil only to 60 kPa for soil with stone column. This is equivalent to CSR of 0.7.
Changes in frequencies from 0.5 to 3 Hz did not significantly influence the permanent strain
of reinforced soil, but these changes do affect the stiffness. Stiffness of the soil with the stone
column was about 25% higher at 3.0 Hz compared to that at 0.5 Hz.
Stone columns also helped reduce pore water pressure build up under cyclic loading by
providing a drainage path. This resulted in an increase in cyclic stress ratio from 0.6 to 0.7. It
also decreases the permanent deformation by about 70% in most tests, when compared to
ii
Samir Ashour Dedication
DEDICATION
This thesis work is dedicated to my parents, who have always loved me unconditionally and
whose good examples have taught me to work hard for the things that I aspire to achieve.
This work is also dedicated to my wife and children, who have been a constant source of
support and encouragement during the challenges of life. I am truly thankful for having you in
my life.
This thesis is dedicated to the sake of Allah, my Creator and my Master, My great teacher and
iii
Samir Ashour Acknowledgments
ACKNOWLEDGMENTS
I would like to acknowledge gratefully for the assistance received from supervisors Dr
Gurmel Ghataora and Professor Ian Jefferson for their constant supervision, guidance and
I would also express my sincere thanks to Lap technicians Mr Sebastian Ballard and Mr Jim
I am very much grateful to Libyan Ministry of Higher Education for granting the scholarship
Finally, I thank my family for their patience and understanding to make this piece work
possible.
iv
Samir Ashour Table of Contents
TABLE OF CONTENTS
ABSTRACT ............................................................................................................................................. i
DEDICATION ....................................................................................................................................... iii
ACKNOWLEDGMENTS...................................................................................................................... iv
TABLE OF CONTENTS ........................................................................................................................ v
LIST OF FIGURES.............................................................................................................................. xiv
LIST OF TABLES ............................................................................................................................... xiv
LIST OF ABBREVIATION ................................................................................................................. xv
CHAPTER 1............................................................................................................................................ 1
1 INTRODUCTION ........................................................................................................................... 1
1.1 Background ...................................................................................................................... 1
2 CHAPTER 2.................................................................................................................................... 6
2 LITERATURE REVIEW ................................................................................................................ 6
2.1 Introduction ..................................................................................................................... 6
v
Samir Ashour Table of Contents
3 CHAPTER 3.................................................................................................................................. 57
3. METHODOLOGY ........................................................................................................................ 57
3.1 Introduction ................................................................................................................... 57
3.3 Physical Properties of the Model Stone Column and Surrounding Soil Material ......... 59
vi
Samir Ashour Table of Contents
vii
Samir Ashour Table of Contents
viii
Samir Ashour List of Figures
LIST OF FIGURES
Figure 2.1 Deep vibrator details (Moseley and Kirsch, 2004) ................................................. 11
Figure 2.5 Unit cell concept (Barksdale and Bachus, 1983) .................................................... 21
Figure 2.6 Stone Column arrangements in (a) Triangular (b) Square (Balaam & Booker, 1981)
.................................................................................................................................................. 22
Figure 2.8 Settlement diagram for stone column in uniform soft clay (Greenwood, 1970) .... 27
Figure 2.10 Load – settlement curves for stone column (Balaam et al., 1977) ........................ 34
Figure 2.11 Vesic’s cylindrical cavity expansion factors (Barksdale and Bachus 1983)......... 37
Figure 2.12 Failure mechanism of single stone column (Barksdale and Bachus, 1983) .......... 40
Figure 2.13 Photographs of sand columns beneath circular footing at beginning, middle and
end of foundation loading process: (a) TS-01, 150 mm; (b) TS-02, 250 mm ......................... 42
Figure 2.14 Effect of undrained shear strength of the surrounding soil on the ultimate bearing
Figure 2.15 Effect area replacement ratio on the ultimate bearing capacity of the treated
ground ....................................................................................................................................... 46
Figure 2.17 Design chart for vibro replacement (Priebe; 1998) ............................................... 51
Figure 2.18 Cyclic loading sequences for railway embankment (Kempfert et al., 1990). ....... 53
Figure 3.1 Typical particle size distribution of Kaolin clay and crushed basalt ...................... 60
xiv
Samir Ashour List of Figures
Figure 3.2 Dry density-water content relationship for Kaolin clay .......................................... 61
Figure 3.3 Relationship between undrained shear strength and water content for Kaolin clay62
Figure 3.4 The undrained shear strength – liquidity index relationship ................................... 63
Figure 3.5 Typical stress-strain, direct shear box test (a) wet condition and (b) dry condition64
Figure 3.7 Consolidation curve square root of time step 50 to 100 kPa................................... 67
Figure 3.9 Cofficients of consolidation and compressibilty and stress relationship ................ 68
Figure 3.10 Model I consolidation process for 100 mm dia. Samples .................................... 74
Figure 3.11 Model II consolidation chamber for preparing 300 mm dia. Samples ................. 74
Figure 3.13 Dynamic measurement of subgrade stress (Yoo and Selig, 1979)........................ 79
Figure 3.14 Dynamic stresses mesurments at subgrade level (under 300mm of ballast) ......... 80
Figure 3.15 Typical settlement /√time relationship for 100 mm dia. triaxial sample .............. 83
Figure 3.16 Settlement and total vertical stress relationship for 100 mm dia. triaxial sample. 84
Figure 3.17 Void ratio/ log σv relationship for 100 mm dia. triaxial sample ........................... 84
Figure 3.18 Typical settlement / time relationship for 300 mm specimens ............................. 86
Figure 3.19 Settlement and total vertical stress relationship for 300 mm specimes ................ 86
Figure 3.20 Variation of excess pore water pressure at PPT1, PPT2, and PPT 3 positions
xv
Samir Ashour List of Figures
Figure 3.25 Variation in water content before testing for both Models (I and II) .................... 93
Figure 3.26 Shear strength determinations in test bed (pilot stage) ......................................... 94
Figure 4.2 Stress strain behaviour of soft soil (no column) and soil/stone column composite.
................................................................................................................................................ 101
Figure 4.3 Improvement ratio versus area replacement ratio ................................................. 101
Figure 4.4: Excess pore water pressure of soil (no column) and soil/ stone column composite.
................................................................................................................................................ 102
Figure 4.6 Typical stress-strain curves: (a) soil (no column); (b) soil/stone column composite
................................................................................................................................................ 106
Figure 4.7(b) Typical change in pore water presures-strain curves for soil/ stone column
Figure 4.8 Normalised excess pore water pressure – strain rate relasionship ........................ 108
Figure 4.9 Soil modulus of both soils only and soil/ stone column composit ........................ 110
Figure 4.10 The subgrade modulus of reaction (𝑘𝑠) for both soils (with no column and soil/
Figure 4.11 The stress-strain for different column densities .................................................. 114
Figure 4.12 Improvement ratio versus the relative density .................................................... 115
Figure 4.13 Change of pore pressure with column density .................................................... 115
Figure 4.14 normalised deviator stress against column density ............................................. 116
xvi
Samir Ashour List of Figures
Figure 4.18 Typical results for changes in pore water pressure for soil (no column) specimens
................................................................................................................................................ 126
Figure 4.19 Typical results for changes in pore water pressure for soil/stone column
Figure 4.20 Hand vane shear test results before and after testing .......................................... 127
Figure 4.21 Water content variation before and after testing ................................................. 128
Figure 5.1 Pemanent axial strain and number of cycles relationship: at 0.5Hz loading
Figure 5.2 Pemanent axial strain and number of cycles relationship: at 1 Hz loading frequency
................................................................................................................................................ 133
Figure 5.3 Pemanent axial strain and number of cycles relationship: at 3 Hz loading frequency
................................................................................................................................................ 134
Figure 5.4 Effect of loading application and cyclic stress leve at frequency of 1 Hz ............ 136
Figure 5.5 Effect of loading frequency and number of cycles under cyclic stress of 50 kPa. 137
Figure 5.6 Cyclic deformation during cyclic triaxial test: at 0.5 Hz loading frequency ........ 140
Figure 5.7 Cyclic deformation during cyclic triaxial test: at 1 Hz loading frequency ........... 140
Figure 5.8 Cyclic deformation during cyclic triaxial test: at 3 Hz loading frequency ........... 141
Figure 5.9 Variation of resilient modulus with number of cycles at loading frequncy of 0.5 Hz
................................................................................................................................................ 142
Figure 5.10 Variation of resilient modulus with number of cycles at loading frequncy of 1 Hz
................................................................................................................................................ 143
xvii
Samir Ashour List of Figures
Figure 5.11 Variation of resilient modulus with number of cycles at loading frequncy of 3 Hz
................................................................................................................................................ 143
Figure 5.12 Variation of resilient modulus with respect to dynamic stress and loading
Figure 5.13 Accumulation of excess pore water pressure at loading frequency of 0.5 Hz .... 145
Figure 5.14 Accumulation of excess pore water pressure at loading frequency of 1Hz ........ 145
Figure 5.15 Accumulation of excess pore water pressure at loading frequency of 3Hz ........ 146
Figure 5.16 Effect of cyclic stress on pore water pressure ..................................................... 147
Figure 5.17 Deformed shape for soil only specimes .............................................................. 149
Figure 5.18 Deformed shape for the soil/column specimen ................................................... 149
Figure 5.21 Pore water pressure measurment at the centre of the specime ............................ 153
Figure 5.22 Pore water pressure measurment at 50 mm from the centre of the specime ....... 154
xviii
Samir Ashour List of Tables
LIST OF TABLES
Table 2.2 Equivalent unit cell diameter (adopted from Balaam & Booker, (1981)) ................ 21
Table 2.4 Comparison between predictions of tank’s settlement (Ellouze et al., 2010) .......... 33
Table 2.5 Estimation of ultimate bearing capacity (Bergado et al., 1991) ............................... 36
Table 4.3 Measured and calculated vertical stress obtained from different studies ............... 103
Table 4.6 Deviance, BIC, R-Squared and Adjusted R2 values .............................................. 122
Table 5.1 Summary varibles investigated using cyclic triaxial test ....................................... 132
Table 5.2 Material parameters (after Li and Selig (1996)) ..................................................... 139
xiv
Samir Ashour List of Abbreviation
LIST OF ABBREVIATION
Cc Compression index
cv Coefficient of consolidation
dc Column diameter
De Equivalent diameter
f Frequency
Gs Specific gravity
h Depth of ballast
Ir Rigidity index
xv
Samir Ashour List of Abbreviation
L Column length
LI Liquidity index
Mr Resilient modulus
mv Coefficient of compressibility
N Number of cycles
p Mean stress
PC Subgrade pressure
q Deviator stress
w Water content
σv Vertical stress
xvi
Samir Ashour List of Abbreviation
εp Permanent strain
xvii
Samir Ashour Chapter 1: Introduction
CHAPTER 1
1 INTRODUCTION
1.1 Background
Soil improvement methods have received much attention over the last few decades, as a
result of shortage of good quality land for development, together with the need to develop
transport infrastructure in areas which were previously regarded as being unsuitable. These
often included soft soil deposits. In many instances, the proposed structures on improved/
stabilised soft soils are considered to be viable on economic grounds due to the high cost of
virgin land and the environmental need to site infrastructure away from developed areas.
Soft soil deposits are characterised as having low bearing capacity, and high compressibility
leading to potential instability and large settlement. Therefore, in order to reduce these
problems, either the structure or the underlying soils or both needs to be modified. Most often,
it is cheaper to undertake ground improvement than to modify the structure. Although in some
A range of ground improvement methods can be used. They can be categorised as:
biotechnical stabilization. Sometimes it may be more efficient to use more than one type of
treatment.
Stone columns are constructed by using the same equipment as vibro-compaction. They
are considered as a densification method and also can act as a vertical drain that can lead to
speedier the consolidation progress of the soil. It is also considered as a ground reinforcement
1
Samir Ashour Chapter 1: Introduction
technique due to the granular nature of the column material. Granular columns are installed
into the soft soil by means of compacting gravel or crushed rock into the cylindrical void
created by a vibrating poker. This technique has been used extensively over the last few
decades in many parts of the world. The design and construction techniques have been
improved to an extent that stone columns could be used in some cases as an alternative to the
traditional foundations such as piling (McCabe et al., 2009; Barksdale & Bachus, 1983;
Greenwood & Kirsch, 1984; Mitchell and Huber 1985; Priebe, 1995; Sivakumar et al., 2010).
Many full-scale, laboratory experimental and analytical investigations have been undertaken,
which have led to improvement of the stone column technique and methods of analysis,
providing better prediction for both bearing capacity and settlement (Bergado and Lam, 1987;
Greenwood, 1991). These studies demonstrate that the behaviour of the stone columns was
the area replacement ratio; column spacing; column and surrounding soil stiffness; the stress
ratio of both the column and surrounding soil; and the method of installation. Also they have
shown that columns can fail by bulging, bending, punching and shearing, depending on the
length of the column and the method of applying the loads on top of the column.
Despite of the fact that, the previous investigations have provided a wide understanding on
using stone columns to improve the properties of soft soils, these investigations were focused
the behaviour of stone columns under the application of static and monotonic loading and
there was very limited information on their behaviour when subjected to cyclic loading. The
key difference in this study and the previous ones is the mode of loading, where the applied
loads will vary with time. Thus this research will provide a better understanding of the
behaviour of stabilised stone columns under static and dynamic loading conditions. This
2
Samir Ashour Chapter 1: Introduction
includes evaluation of deformation under loading and mechanisms of stress transfer in both
There is increasing need for higher performance track substructure systems (ballast, sub-
ballast and subgrade) in order to provide constant and uniform support for the sleepers and the
rail allowing them to cope with increased demand for high speed trains. Without any
upgrades, tracks built over soft ground can be problematic, leading to increase in maintenance
Many researchers (Seed et al., 1955; Brown et al., 1975; Li and Selig, 1996; Miller et al.,
2000; Li-Zhong Wang et al., 2011) have studied the effect of cyclic loading on the
compressibility and strength characterisation of soft soil. They demonstrated that such soils
have high compressibility and low bearing capacity. Therefore, under cyclic stresses above
the threshold stress ratio (i.e. ratio between the dynamic stresses to the soils static strength,
which is typically between 0.4 and 0.65 depending on the soil type and conditions), it is
expected to have a high residual settlement resulting in a reduction in strength. Hence leading
In order to overcome this problem, it is generally better (i.e cost-effective) to improve the
ground rather than import materials of suitable quality to replace poor soils. This has
of vibro stone columns can be considered for improving subgrade soil (Raju, 2003; Fatahi et
al., 2012).
Vibro stone columns have been widely used to reinforce and stabilise railway tracks built on
soft soils to control the permanent deformation and the generation of excess pore pressure
3
Samir Ashour Chapter 1: Introduction
(Abdullah et al., 2009; Fatahi et al., 2012). However, studies on the dynamic performance of
stone columns are mostly limited to liquefaction mitigation potential in silty soils (Munfakh,
1984; Priebe, 1995; Munfakh, 2003; Rollins et al., 2009). Additionally, very limited research
has been carried out on the performance of stone columns subjected to repeated vertical
loading. Therefore, aspects such as identifying failure mechanisms and quantifying the
amount of settlement, especially at low replacement ratio (under 10%) require investigation.
For this study, vibro stone column were modelled in both triaxial tests and small scale tank
tests. Stone columns were constructed in soft soils (cu <15 kPa) and loaded to failure under
both monotonic and dynamic loading. The main objective of the testing program was to
examine the effect of loading rate, loading frequency, and cyclic stress ratio on the behaviour
of both soft soils and stone column reinforced soil. Thus, allowing for a better understanding
of the response of stone column under dynamic loading application can be achieved.
This research aim is to assess the behaviour of stone column subjected to cyclic loading as a
approached by developing the appropriate loading frequency and dynamic stresses on both
treated and untreated soft soils in order to investigate the stiffness, permanent deformation
and pore water pressures. In order to achieve the above aim, the following objectives were set:
ii. Examine the effect of loading strain rate during the monotonic conditions on the
iii. Examine the effect of relative density of the column material on column bearing
capacity;
4
Samir Ashour Chapter 1: Introduction
iv. Investigate the effect of loading frequency on the behaviour of stone column;
v. Investigate the effect of dynamic stress level on the behaviour of stone column.
Brief descriptions of the remaining chapters of this thesis are outlined below:
Chapter 2 provides a brief background of stone column foundations and reviews previous
The properties of the materials used; the design and manufacture of testing apparatus and
instrumentation used in this laboratory base study; and the specimen preparation and general
The monotonic testing programme and results of the effect of varying strain rate and column
density are described in Chapter 4. Results from current study are analysed and compared
with other research findings from both laboratory tests and field studies. In addition to this,
the failure mechanisms for a single stone column reinforced foundation in both models
Chapter 5 shows results and discussion of tests on stone columns, subjected to cyclic loading.
(i.e the effect of both frequency and cyclic stress level) in order to provide a further
5
Samir Ashour Chapter 2: Literature Review
2 CHAPTER 2
2 LITERATURE REVIEW
2.1 Introduction
Vibro stone column method has been widely used to improve soft soils over the last four
decades. This chapter reviews and the design concepts, installation methods, failure
Geotechnical engineers are often faced with the challenges of building on poor ground and
have to match the level of ground improvement with project requirements. Apart from
abandoning the project, Mitchell and Jardine, (2002) suggested four alternatives:
- bypass the area of poor ground laterally by relocating the facility, or vertically by
using piles or deep foundations;
- modify the natural condition of the poor ground to meet the project requirements.
The rising cost of land and increasing awareness of the impact that construction has on the
environment have been major contributors to the increasing use of ground improvement
techniques.
Ground improvement techniques generally aim to modify some of the soil characteristics in
increasing the density and shear strength in order to: improve bearing capacity; reduce soil
6
Samir Ashour Chapter 2: Literature Review
the rate of consolidation; and improve the homogeneity of soil. Therefore, ground
improvement can be defined as “the controlled alteration of the state, nature or mass
A wide range of ground treatments are available. Different classifications are proposed
depending on: equipment used; required period of treatment (temporary or permanent), or the
purpose of use (Munfakh, 1997; Van Impe et al., 1997; Mitchell and Jardine, 2002; Munfakh,
2003). Munfakh (1997) divides ground improvement methods into the following categories:
Chemical stabilization: permeation grouting; jet grouting; deep soil mixing; and
lime/ cement columns.
Vibro stone columns can be considered one of the most effective methods to improve
performance of the soft ground. stone columns have the advantages of accelerating
consolidation, increasing the bearing capacity of poor ground, and improving the slope
drainage path allowing the pore water pressure to dissipate. Therefore, they have been used in
a verity of applications, such as road and railway embankments, tanks and marine structures
7
Samir Ashour Chapter 2: Literature Review
In 1937 there was the first reported application of vibro compaction used to densify a loose
sand deposit under a building in Berlin by penetrating a vibrating poker to the ground
In the 1950’s this technique was adopted in the UK, however, its application was limited due
to the geological nature of the ground (soft soil; shear strength < 15 kPa). This led to the
development of vibro compaction technique to cover the treatment of fine grained and
cohesive soils by introducing a coarse granular backfill into a vertical void in soft ground
formed using a vibratory poker. This early method of stone column construction received
further improvement in Europe and the United State of America until it became a well-
Over the last three decades this technique has been widely used for both compaction of
cohesion-less soils and for reinforcing soft soils with granular columns. The method has been
shown to be an effective and reliable ground improvement method, especially, in cases, where
high sensitivity to settlement is not critical. It has proven to be a good alternative to the
traditional deep foundation methods such as piles (Serridge, 2006; McCabe et al., 2009).
Soil type and condition are the key factors that govern the method of installation of stone
columns to ensure the required quality of improvement (McCabe et al., 2009). The commonly
applicable installation techniques for soft soils (and the equipment used) are discussed below.
8
Samir Ashour Chapter 2: Literature Review
cu L dc
Installation method Commentary
(kPa) (m) (mm)
Stone columns are constructed first by creating a vertical void in soft ground using a vibratory
poker (Figure 2.1). As the poker is slowly withdrawn, stone is inserted in the void in typical
9
Samir Ashour Chapter 2: Literature Review
stages of 300 to 500 mm thickness and then compacted by the vibrating poker for 30 to 60
The equipment used to form the stone column typically comprises a hollow cylindrical steel
pipe with diameter ranging between 300 mm and 450 mm, and length ranges from 2.00 to
3.50 meter. The vibrator body unit consists of two main parts: the upper part known as the
extension or follower tube as shown in Figure 2.1. The main function of the follower tube is
to allow deep penetration into the ground (Watts, 2000; Moseley and Kirsch, 2004).
The lower section of the vibrating unit contains an eccentric weight which is powered by a
motor at the top of an axial shaft (typical motor power capacity ranges between 50 to 150 kW
and operating frequency ranging from 10 to 30Hz). This eccentric weight typically ranges
between 15 and 40 kN depending on the level of improvement and the condition of the
ground. There is a vibratory force of 150 to 700 kN can be transmitted from the vibrator
casing to the surrounding soil. (Greenwood and Kirsch, 1984; Raju and Sondermann, 2005).
In order to improve the installation process with a high level of quality control, the equipment
was modified by including useful features such as bottom feed delivery system, and
computerised monitoring system that could transfer data from site (Slocombe et al., 2000).
10
Samir Ashour Chapter 2: Literature Review
This technique has been used successfully in soft soils (fines content over 20 %; and
undrained shear strength of soil, cu, in range of 15 to 35 kPa), it can also be used for deep
ground treatments below water level (Moseley and Kirsch, 2004; McCabe et al., 2009).
In order to increase the bearing capacity and the stiffness of soft soil, a 10 to 35 % of the soil
may need to be replaced with uniformly graded stone, column diameters ranging between 0.6
11
Samir Ashour Chapter 2: Literature Review
and 1.2 meter and column depths between 5 and 15 meters. In this method it is preferable to
build columns that extend to end bearing stratum (Barksdale and Bachus, 1983).
The installation sequence of this technique is illustrated in Figure 2.2. The vibrator is inserted
into the soil with the assistance of jetted water. Both the flow of flushing water and the
vibration can reduce soil resistance and allow the vibrator to penetrate the soil under its own
weight. Once the required depth is achieved, granular backfill is placed from the bottom of the
hole to the top where gravel falls against a continuous upward flow of water. As the gravel is
approximately 0.50 meter intervals. A continuous water flow protects the bore from the
collapsing. As the vibrational energy is dissipated in radial waves through the backfill
materials into the surrounding soil and additional expansion of the column is prevented by the
passive resistance of the cohesive soil. (Greenwood and Kirsch, 1984; Watts, 2000).
However, due to environmental issues relating to disposal of flush arising, application of this
12
Samir Ashour Chapter 2: Literature Review
Dry top feed method is considered to be the simplest installation method, and most preferred
in the case of shallow to medium treatment depths of stable soils with undraind shear strength
larger than 30 kPa. This method is not applicable for treatment of soft clay soils due to the
fact that the borehole will collapse when the poker withdrawn due to lack of lateral support
In this method, (Figure 2.3), once the required depth is reached the poker is lifted and the
granular backfill of material introduced from the ground surface into the hole. The vibro-float
is reinserted and stone in the column is compacted. Typical geometries of formed stone
columns are usually 400 to 800 mm diameter and 10 to 15 meters in length (Munfakh, 1984;
McKelvey, 2002).
13
Samir Ashour Chapter 2: Literature Review
This method was introduced into the UK in the 1980’s and currently considered the most
commonly adopted installation method for stone column (McCabe et al., 2009). In terms of
application it is used in soils (cu ranges between 15 and 50 kPa), and unaffected by the
presence of ground water. During the backfill process the vibro-float remains inside the bore
providing more stability and preventing unwanted inclusion. The stone material is supplied
using a hopper to a pipe fixed to the side of the vibro-float shown in Figure 2.4. This can
allow the column to be constructed clean with the same back fill material. As with the top
feed method the stone is then compacted by repeated withdrawal and insertion of the poker
(McCabe et al., 2009). The average treatment depth achieved by this method is 15 meter
(Watts, 2000).
14
Samir Ashour Chapter 2: Literature Review
As a result of the variety of equipment and construction methodologies and due to the lack of
energy consumption of the vibrator, and air/water pressure (Moseley and Kirsch, 2004).
Continuous recording of parameters as a function of time can provide reliable data that can be
used to control stone column construction and help prevent problems. This includes the
possibility of having to alter the design of columns and ensures column uniformity over entire
length.
In addition to the above, cone penetration test (CPT) and standard penetration test (SPT) are
generally used to confirm the achievement of a satisfactory level of vibro compaction (Bell,
2004; McCabe et al., 2009). However, it should be noted that, whilst the data collected during
column construction are for quality control, the penetration tests are used to ensure the degree
In order to evaluate the compressibility of treated soil, it is better to conduct compaction tests
within at least one week of the compaction work to prevent the effect of ageing factor, which
influences the soil properties and causes an increase to its strength ranging from 50 to 100%.
This increase in strength occurs as a result of the reduction in water pressure, and also could
be due to the rearrangement of the physical and chemical bonding forces between the soil
Stone columns are primarily used to reduce settlement. Therefore, large-scale loading tests,
which can be carried out by loading a rigid plate on top of one or more columns are used to
assess the degree of improvement in terms of settlement. However, as this test is costly and
15
Samir Ashour Chapter 2: Literature Review
Recently, geophysical methods including continuous surface wave have been used to assess
settlement improvement (Madun et al., 2012). Although there is no large database to validate
the data collected by geophysical investigation, this method has an advantage of covering
large areas and could be used for both short and long term investigations, which could make it
cost effective (Redgers et al., 2008; Madun, 2012; Madun et al., 2012).
The key parameter in stone column design is the angle of shearing resistance, which is a
function of several parameters such as degree of compaction, grading and material strength
(McKelvey and Sivakumar, 2000; Jefferson et al., 2010). Most of the design methods
developed for predicting the bearing capacity of soils improved by stone column
reinforcement (Hughes and Withers, 1974; Greenwood and Kirsch, 1984) are related to the
undrained shear strength of surrounding soil and the internal angle of friction of the stone
materials; a reduction of 10o in the internal angle of friction can reduce the bearing capacity
In general, 45o can considered to be the maximum friction angle, while in the UK 40o
considered to be the most typical value used in the design (Serridge, 2006). Therefore, it is
important that the stone aggregate within the column have sufficient shear resistance whereby
the particles are strong enough to withstand local stress concentrations which occur during
construction. In addition aggregate used must be durable in the long-term (Watts, 2000).
16
Samir Ashour Chapter 2: Literature Review
There are three main sources of aggregate, which can be used in stone column construction.
Until recently natural or primary stone aggregate was the main source. However, increasing
demands for sustainable construction has led to the use of recycled aggregate as a second
source (e.g. construction demolitions and recycled railway ballast) and when possible the
secondary aggregates (i.e. industrial processed products such as steel slag, rock waste, and
Stone columns were installed in closely spaced groups beneath pad and strip foundations to
provide bearing pressures of 100 kN/m2 (associated with the main portal frame structure) and
on a square grid pattern beneath ground-bearing floor slab areas to provide bearing pressures
Serridge (2005), provides a detailed case history of a project in Coatbridge, Scotland, where
recycled (crushed concrete) aggregate was used successfully for stone column construction.
Columns with diameter of 600 mm and varied in lengths from 2.5 to 6 meters were installed
in a square pattern beneath pad and strip foundations to provide bearing pressure of 100
kN/m2 (for the main frame structure) and 30 kN/m2 under floor slab areas. The average
deformation modulus of 48 MN/m2 was achieved when assessed using plate load tests. He
suggested that, this value was comparable with typical results if natural aggregate was used in
McKelvey et al. (2002) examined the shear strength of different recycled aggregates (i.e.
quarry waste, crushed concrete, and building debris). The performance of these materials
were compared with the performance of crushed rock primary aggregate, in terms of the
influence of dry and wet conditions and the effect of the fines content. They found that the
primary aggregate performed better than the other recycled aggregate in all test conditions.
17
Samir Ashour Chapter 2: Literature Review
For example, in the dry condition, the angle of internal friction for the primary aggregate was
51o, whereas for the secondary aggregate was ranging from (46 o to 37o). These values were
decreased by about 6 % in the wet condition and by about 30 % when 20 % of Kaolin slurry
was add to the aggregate. This study indicated that the internal angle of friction of natural
aggregate had been influenced by all the tested conditions while the quarry waste showed
independency toward the slurry content and affected by the water content, where it showed
different trend in case of crushed concrete aggregate, where it was influenced by the slurry
This could indicate that, regardless of the stone source used, the achievable quality of the
stone column is mostly dependent on the quality of aggregate adopted. The ICE (1987) and
BSI (2005) suggested that in order to achieve a good interlock between the stone particles and
to allow an adequate level of drainage through the column, the grading of aggregates should
range between 20 to 75 mm, also the percentage of fine content should not exceed 10 %.
Stone column technique has been used successfully in many projects around the world for
consolidation process; increase the bearing capacity and to mitigate liquefaction) in various
types of soil from coarse gravel to fine sand, silt and clay (Woodward, 2004). The same
benefits have been found when this technique was also applied to improve the characterisation
of soft marine clay, non-engineering fills and layered soils (McKelvey, 2002; Raju and
18
Samir Ashour Chapter 2: Literature Review
2.5.1 Limitations
Although stone columns have been installed successfully in soils with undrained shear
strengths ranging between 10 kPa and 50 kPa (Hu, 1995; McKelvey, 2002), they are not
suitable for soils with very low undrained shear strength and in soils with thick layers of peat.
This is due to insufficient lateral support that these materials provide and excessive settlement
that could occur due to the loss of water content during consolidation (Greenwood and Kirsch,
1984). Barksdale and Bachus (1983) suggest that stone columns can be more effective when
used for soil stabilisation rather than structural foundations. Most recommended treatment
depths are in the range of 6 to 15 m although columns have been constructed to a depth
greater than 30 m (Watts, 2000). Floating columns are not recommended to be used in weak
soils and therefore, a competent end bearing for stone columns is generally specified (Killeen
As the column and the surrounding soils work together and share the stresses subjected to
them, the bearing response of the soil/ column system is influenced by the properties of both
materials and their interaction between each other. In order to theoretically solve this complex
problem a certain level of idealisation (unit cell idealisation) is used. Most existing theories
consider stone and clay as perfect elastic or elastic-plastic materials (Hughes and Withers,
1974; Barksdale and Bachus, 1983; Priebe, 1995). Therefore the design of stone columns
requires cognisance of parameters such as used grid pattern, area replacement ratio and stress
concentration factor (Sections 2.6.1.1 to 2.6.1.3). The typical design procedure suggested by
19
Samir Ashour Chapter 2: Literature Review
1. Predict the ultimate bearing capacity of the single column using the undrained shear
strength of the soil and the internal angle of friction,
2. Determine the column diameter and the allowable spacing between the column,
Balaam et al. (1977) stated that the unit cell idealises for the behaviour of single column and
the soil around it, and the behaviour of one unit cells within a group is the same. This concept
is used to determine the area within which vertical stresses are considered to be distributed
(Barksdale and Bachus, 1983). Stone columns can be installed in three main patterns: triangle,
square and rectangular patterns. The equilateral triangle pattern is considered the most usable
and efficient arrangement with regards to large areas and uniform rate of densification, while
the square and rectangular pattern is applied in condition of isolated spread footing (Watts,
Typical layouts of stone columns in different patterns are shown in Figure 2.5. It is
convenient to associate the tributary area of soil surrounding each stone column. The tributary
area can be closely approximated as an equivalent circle (unit cell) with equivalent diameter
The spacing (S) between the columns is a critical factor since it influences the degree of
cover a wide area of treated soils. Greenwood (1970) stated that, in practice, a narrow spacing
20
Samir Ashour Chapter 2: Literature Review
However, the spacing is generally ranges from 2 to 3 times the column diameter (Hughes and
Withers, 1974).
Depending on the stone column arrangements, spacing and the equivalent diameter can be
Table 2.2 Equivalent unit cell diameter (adopted from Balaam & Booker, (1981))
12 0.25
Triangular ( 2) 𝑆 ≅ 1.05 𝑆
𝜋
16 0.25
Square ( 2) 𝑆 ≅ 1.13 𝑆
𝜋
De
S
D
Plan view
De Stone column
Rigid
frictionless side
L
Unit cell
21
Samir Ashour Chapter 2: Literature Review
1⁄
12 4
De 𝐷𝑒 = ( ) 𝑆
𝜋2
=1.05 S
1⁄
16 4
𝐷𝑒 = ( ) 𝑆
De 𝜋2
=1.13 S
Figure 2.6 Stone Column arrangements in (a) Triangular (b) Square (Balaam & Booker, 1981)
Based on the concept of unit cell, it is possible to use the geometry of a mesh of stone
columns in order to determine the amount of soil replaced and the area replacement ratio As,
which is defined as the ratio between the area of each column Ac and soil area A (Barksdale
2
𝐴𝑐⁄ 𝐷𝑐
𝐴𝑠 = 𝐴 = [ ⁄𝐷𝑒 ] 2.1
Where: Dc is the compacted stone column diameter; and De is the equivalent unit cell
diameter.
The area replacement ratio, also can be presented in terms of the diameter and spacing of the
22
Samir Ashour Chapter 2: Literature Review
𝐷
𝐴𝑠 = 𝑐𝑔 〈 𝑆𝑐 〉2 2.2
In general, increasing the area replacement ratio leads to improved behaviour of the
composite soil (i.e. increased bearing capacity and reduced settlement). For improvement in
bearing capacity larger than 30 %, Wood et al. (2000) indicated that the area replacement ratio
A large portion of the stress concentration is transferred from the loaded foundation to the
stone column. The remaining stresses are transferred to the soil as it is weaker than the
column material. The vertical stress distribution within the unit cell (concentration factor) is
defined as the ratio of the stress on the stone column to that on the surrounding soil within a
unit cell (Han & Ye, 1991). The magnitude of the concentration factor generally ranges from
𝜎𝑠
𝑛= ⁄𝜎𝑐 2.3
Where: σs is stress in the stone column; and σc is the stress in the surrounding soil.
There are several factors that could influence the magnitude of stress concentration such as
type of foundation, length of the column and time of consolidation. Juran and Guermazi
(1988) stated that the magnitude of (n) increases with time of consolidation and decreases
with the length of the column. in addition, loading the soil/column system through a rigid
23
Samir Ashour Chapter 2: Literature Review
foundation would lead to larger stress concentration than if flexible a foundation was used
Ambily and Gandhi (2007) indicated that the modular ratio had a direct impact on the stress
concentration, they found that the n factor increased with any increase in modular ratio; and
The total vertical stresses σ in stone columns and the surrounding soil over the unit cell area
𝜎 = 𝜎𝑠 . 𝐴𝑠 + 𝜎𝑐 (1 − 𝐴𝑠 ) 2.4
Thus the stress concentration factor can be estimated (as follow) assuming equal vertical
displacement would occur and using the elastic theory as a function of the modular ratio of
the column and the soil (Barksdale and Bachus, 1983; Babu et al., 2013).
𝜎
𝜎𝑐 ≤ 𝜇𝑐 𝜎 = 𝜇𝑐 ( 𝑠⁄𝜇𝑠 )2 2.5
Where: 𝜇𝑐 𝑎𝑛𝑑 𝜇𝑠 are the stress ratios in the clay and the stone column respectively.
𝜇𝑠 = 𝑛⁄[1 2.6
+ (𝑛 − 1)𝐴𝑠 ]
24
Samir Ashour Chapter 2: Literature Review
One of the main purposes of using stone columns is to reduce the overall settlement of the
treated soil. This reduction is governed by the magnitude of applied load distribution that
In order to satisfy the unit cell concept condition (Figure 2.7), most of the design approaches
for predicting settlement of soil/column system assume an infinitely wide loaded area
reinforced with granular column with a constant diameter and spacing. The column is
assumed to be loaded via a rigid plate and ending in an undeformable bearing layer (Van
There are various methods for estimating the settlement of the stone column/ soil system.
25
Samir Ashour Chapter 2: Literature Review
1 0.5 + 𝑓(𝜇, 𝐴𝑠 )
= 1 + 𝐴𝑠 [ − 1]
𝑅 (𝑘𝐴 )𝑠 𝑓(𝜇, 𝐴𝑠 )
Priebe Priebe
1 − 𝜇2 (1 − 2𝜇)(1 − 𝐴𝑠 )
method (1976, 𝑓(𝜇, 𝐴𝑠 ) = [ 2
][ ] (𝑘𝐴 )𝑠 Overestimated results
1 − 𝜇 − 2𝜇 1 − 2𝜇 + 𝐴𝑠
(2.6.2.3) 1995)
∅
= tan2 (45 − )
2
Finite
Balaam et Its accuracy dependent on
element [𝐾𝐸 ]{∆𝜎 (𝑚−1) } = {{∆𝐹𝐸 } + {𝐾𝑐(𝑚) } {∆𝜎 (𝑚) } + {∆𝐹𝐷𝑁
(𝑚)
}}
al. (1977) the input parameters
method
Greenwood (1970) was one of the first to introduce empirical design curves to estimate the
settlement of stone column/ soil systems under widespread loading. The empirical curves
showed in Figure 2.8, represents the settlement reduction as a function of column spacing and
the undrained shear strength of the surrounding soil, which was limited between 20 and 40
kPa. The shaded zone in the curve represents the expected settlement reductions when
applying a wet process of construction. Greenwood (1970) suggested that these charts should
be used with caution within the indicated range. Although Balaam and Booker (1985)
indicated that these curves showed a close agreement with their finite element analysis
26
Samir Ashour Chapter 2: Literature Review
approach for settlement estimation, they suggested that at column spacing of over 2.5 m, the
Figure 2.8 Settlement diagram for stone column in uniform soft clay (Greenwood, 1970)
This method has been used (mainly in Japan) to estimate the settlement of sand compaction
piles (Barksdale and Bachus, 1983). Aboshi et al. (1979) presented this simple approach for
estimating the reduction in settlement of ground reinforced by stone columns based on the
unit cell assumptions and the one dimensional consolidation theory. This approach require
estimation of stress concentration factor n using previous work experience and the past results
27
Samir Ashour Chapter 2: Literature Review
The variation of the vertical stress in the clay (σc) due to the applied external stress can be
expressed by:
𝜎𝑐 = 𝜇𝑐 𝜎 2.8
Where: σ is the average external applied stress; and µc is the ratio of stresses in the clay (Eq.
2.7).
By applying the one dimensional consolidation theory the primary consolidation settlement
𝐶 𝜎0′ + 𝜎𝑐
𝑆𝑡 = (1+𝑒𝑐 ) . log10 ( ). 𝐻 2.9
0 𝜎0 ′
Where: St is the primary consolidation settlement over distance H of stone column treated
ground; H is the height of stone column; e0 is the initial void ratio; Cc is the compression
index; σ0 is the average initial effective stress; and σc is the changing in stress in the clay
layer.
Additionally, the following equation may be used to estimate the settlement of unreinforced
𝑆 = 𝑚𝑣 . 𝜎. 𝐻 2.10
Where, S and St, as defined previously, are the consolidation settlements of unreinforced and
28
Samir Ashour Chapter 2: Literature Review
The ratio of the settlement of the stone column (i.e. improved soil to the unimproved soils)
𝑆𝑡⁄ 1
𝑅= 𝑆 = [1 + (𝑛 − 1)𝐴𝑠 ] 2. 12
Where n is the stress concentration factor and As is the area replacement ratio.
The above equation showed that the settlement ratio is a function of the area replacement ratio
and the stress concentration factor, and the settlement of the treated soil decreases as those
equilibrium method would be sufficient for preliminary design if the stress concentration
Priebe method is considered one of the most widely adopted semi-empirical methods used to
obtain the improvement factor (estimating settlement reduction due to soil treatment with
stone column). This method has undergone a series of improvements and modifications to its
additional factors such as the effect of compressibility of the column material and
confinement from overburden, which considers the effect of unit weight of both column and
soil materials by adding the depth factor to the design calculation (Barksdale and Bachus,
The original solution proposed contained a number of simplifying assumptions that are listed
below:
29
Samir Ashour Chapter 2: Literature Review
long elastic hollow cylinder solution is assumed in order to find out the radial
deformation of the soil and;
According to these assumptions it can be concluded that the column is designed to only fail in
Based on these assumptions, Priebe (1995) presented in Figure 2.9 the ratio of settlement of
untreated to treated ground S/St as a function of the area replacement ratio As and internal
Therefore, the radial deformation of the column/soil system can be determined from the
solution of expanding cylindrical cavity in half elastic. This expression gives the radial
deformation of the column based on vertical stress acting on it and on the ground (Priebe,
1995). Establishing equality of vertical deformations in the soil and the column and by the
equilibrium condition of vertical stresses, the following improvement factor n0 expression can
be obtained:
𝐴𝑐 1⁄ +𝑓(𝜇 ,𝐴 /𝐴)
𝑠 𝑐
𝑛0 = 1 + . ( 𝐾 2 .𝑓(𝜇 − 1) 2.13
𝐴 𝑎𝑐 𝑠 ,𝐴𝑐 /𝐴)
(1−𝜇𝑠 ) . (1−𝐴𝑐⁄𝐴)
Where: 𝑓(𝜇𝑠 , 𝐴𝑐 /𝐴) = 𝐴 ; 2.14
1−2𝜇𝑠 + 𝑐⁄𝐴
30
Samir Ashour Chapter 2: Literature Review
µs is the Poisson’s ratio; and (𝐾𝑎𝑐 ) coefficient of active pressure of the column is
defined as:
∅𝑐⁄
(𝐾𝑎𝑐 ) = 𝑡𝑎𝑛2 (45𝑜 − 2) 2.15
Although this method has been widely used in ground improvement industry, it has
limitations in settlement prediction, such as not being suitable for estimating the behaviour of
Ellouze et al. (2010) reported that some of the assumptions in Preibe’s method, such as that
the unit weights of column material and initial soil are neglected, are questionable and not
well defined. Additionally, they highlight that the design procedure is not connected
mathematically and inconsistence. For example, in the first step of the design, Priebe (1995)
considers a cylindrical cavity subjected to lateral expansion during which zero vertical
second step, this solution is incorporated in the soil/column model for which there is a
distribution of vertical stress generating non-null vertical deformation, and, consequently, the
31
Samir Ashour Chapter 2: Literature Review
Using a case study of 54 m diameter tank built on reclaimed ground reinforced by stone
columns to ensure the stability of the tank at Zarzis terminal (Tunisia), Ellouze et al. (2010)
demonstrated that Priebe’s method was overestimated the settlement of reinforced soil by
stone columns comparing with other linear elastic models (Balaam and Booker (1981); Chow
(1996) Bouassida et al. (2003) and the French recommendations (Françaises, 2005)).
Working load of the tank was approximated as a quasi-uniform stress of 120 kPa. The
reinforcement was performed along an average depth of 7 meters with a columns diameter of
1.2 meter installed in a triangular pattern. The soil was improved on a circular area with a
The settlement of unreinforced soil was estimated to be about 230 mm underneath the centre-
line of the tank, whereas the settlement at the edge of tank was estimated to be about 60 mm.
32
Samir Ashour Chapter 2: Literature Review
The settlement predictions of the reinforced soil at the centre-line of the tank and its edge that
were obtained using the software columns (Bouassida et al., 2009) and Priebe’s method are
Table 2.4 Comparison between predictions of tank’s settlement (Ellouze et al., 2010)
Settlement of Settlement of
Settlement reduction Settlement reduction
reinforced soil at reinforced soil at
Method factor at the centre factor at the edge of
the centre line of the edge of the
line of the tank the tank
the tank, (mm) tank, (mm)
Recorded - 30 - 2.00
French recommendation
55 26 4.18 2.30
(Françaises, 2005)
The finite element method could provide an appropriate theoretically approach to model stone
column reinforced ground, where nonlinear material properties, the column/ soil interface and
the boundary conditions can all be sufficiently modelled (Gniel and Bouazza, 2007; Killeen,
2012; Babu et al., 2013; Killeen and McCabe, 2014; Mohanty and Samanta, 2015). Most of
the finite elements studies have utilized the axisymmetric unit cell model to analyse the
conditions of either a uniform load on a large group of stone columns or a single stone
column. Balaam et al. (1977) used finite element method to investigate the behaviour of stone
column, and developed a design curves for predicting settlement reduction. Figure 2.10
showed an example of these curves. Balaam et al. (1977) highlited that column diameter and
33
Samir Ashour Chapter 2: Literature Review
ratio of column length to depth of the soil layer had significant influence on settlement
q ̸ cu
δ ̸d
Figure 2.10 Load – settlement curves for stone column (Balaam et al., 1977)
Barksdale and Bachus (1983) developed a design curves for estimating settlement in both low
compressibility and compressible soils using a nonlinear finite element method combined
with the idealisation of unit cell to predict the primary consolidation settlement of column/soil
system. In this approach the soil modulus of elasticity was assumed to be constant with the
depth of the soil. In addition to this, different area replacement ratios (10, 20, 25, and 35 %) at
different L/D ratios (5, 10, and 20) were considered and each case has a different design chart.
Poorooshasb and Meyerhof (1997) introduced an elastic soil model to predict the settlement
reduction of raft foundation resting on end bearing stone column reinforced soft soils. This
model indicated that the design pattern, spacing and the degree of compaction of the column
material were the most viable parameters to control the level of settlement.
34
Samir Ashour Chapter 2: Literature Review
Generally, finite element solution claimed by authors (Ambily and Gandhi, 2007; Killeen,
2012; Killeen and McCabe, 2014) to provide a good agreement with the full scale and site
investigation results; however, their utility is governed by the accuracy of the input
parameters.
Stone columns are often constructed penetrating through soft soil to end bearing layers with a
critical length of 4 to 6 times the column diameter. Thus the most possible mode of failure
that might develop for an isolated stone column over a depth of 2 to 3 diameters from the
surface is bulging failure (Barksdale and Bachus, 1983; Bergado et al., 1991). Several
approaches have been developed to predict the ultimate capacity of an isolated stone column
surrounded by soft soils. Table 2.5 summaries the most common of these corresponding to the
These approaches were developed on the basis that the lateral confining stress supporting the
column is usually considered as the ultimate passive resistance which the surrounding soil can
mobilize as the column bulges. Since then the ultimate vertical stress(𝜎𝑣 ) can be predicted as
the coefficient of the passive pressure of the column (𝑘𝑝 ) times the lateral confining pressure
(𝜎𝑟 ).
𝜎𝑣 = 𝜎𝑟 ∗ 𝑘𝑝
1 sin ΄
= 𝜎𝑟 ( ) 2.16
1 sin ΄
35
Samir Ashour Chapter 2: Literature Review
1 + sin ∅𝑠
𝑞𝑢𝑙𝑡 = (𝛾𝑐 𝑧𝑘𝑝𝑐 + 2𝑐𝑜 𝑧√𝑘𝑝𝑐 ) Greenwood (1970)
1 − sin ∅𝑠
1 + sin ∅𝑠
𝑞𝑢𝑙𝑡 = (𝐹𝑐, 𝐶𝑜 + 𝐹𝑞′ 𝑄𝑜 ) Vesic (1972)
1 − sin ∅𝑠
Bulging
1 + sin ∅𝑠
𝑞𝑢𝑙𝑡 = (𝜎𝑟𝑜 + 4𝐶𝑢 ) Hughes and Withers (1974)
1 − sin ∅𝑠
3 2
General shear 𝑞𝑢𝑙𝑡 = 0.5𝛾𝑐 𝐵𝑡𝑎𝑛𝜓 + 2𝐶𝑜 𝑡𝑎𝑛𝜓 + 2 (1 − 𝑎𝑠 )𝐶𝑜 𝑡𝑎𝑛𝜓
Barksdale and Bachus (1983)
tan(𝜇𝑠 𝑎𝑠 tan ∅𝑠 )
𝜓 = 45𝑜 +
2
Greenwood (1970) noted that there was no exact mathematical method to estimate the bearing
capacity of cohesive soils treated by stone columns, because of the dilation that occurs within
the column and the resulting lateral stress to the surrounding soil which can be resisted by
passive pressure; (Greenwood, 1970) hypothesed that the column will behave as if it was in a
triaxial chamber, and the degree of improvement in the bearing capacity would be governed
by the lateral support from the surrounded clay to the column and the internal angle of
friction. In addition, Greenwood (1970) highlighted that the carrying capacity of the column
increases until either a local shear failure in the clay or end bearing failure at the bottom of the
column accords.
Hughes and Withers (1974) used a laboratory based model and radiography device, to study
the behaviour of both sand columns and surrounding clay by tracking the deformations
occurred within and outside the column. They concluded that the cylindrical cavity expansion
theory can be used to define the column behaviour. They proposed the following equation for
36
Samir Ashour Chapter 2: Literature Review
1+sin ∅′ ′
𝜎𝑣 = (𝜎𝑟𝑜 + 4𝑐) 2.16
1−sin ∅′
Vesic (1972), based on the theory of cylindrical cavity expansion, developed an expression
for estimating lateral resistance where the ultimate lateral resistance is:
σ3 = c Fc + q Fq 2.17
where: c is the cohesion of the soil, q is the mean stress at the equivalent failure depth and
Fc, Fq are cavity expansion factors which determine from Figure 2.11 using the rigidity index,
Ir,:
E
Ir = 2.18
2(1 v)(c q tan ΄ )
Where: E is the modulus of elasticity of the soil in which cavity expansion occur, c is the
cohesion of the soil, v is the Poisson’s ratio of the soil, ΄ is the effective angle of friction of
the soil and q is the mean stress within the zone of failure.
Figure 2.11 Vesic’s cylindrical cavity expansion factors (Barksdale and Bachus 1983).
37
Samir Ashour Chapter 2: Literature Review
Therefore assuming (𝜎3 = 𝜎𝑟 ) substituting (2.18) in (2.16) the ultimate stress that can be
1 sin ΄
𝜎𝑣 = [c Fc + q Fq] ( ) 2.19
1 sin ΄
The mean stress q used in the above equations should be taken as stress at the average depth
of the bulge taking into account the initial and the final stresses in the ground. Vesic (1972)
expressions can be used both for short and long term calculation.
Barksdale and Bachus (1983) developed a simple approach for estimating the bearing
capacity of a single column. (Equation 2.21) based on authors ‘past experience and utilising
𝑞𝑢𝑙𝑡 = 𝑐𝑢 ∗ 𝑁𝑐 2.20
Where, 𝑞𝑢𝑙𝑡 is the ultimate bearing capacity of the stone column; 𝑐𝑢 is the undrained shear
strength of the surrounding soil; and 𝑁𝑐 is the bearing capacity factor for the stone column
usually ranging between 18 and 22, and 5 for estimating the bearing cabacity of the untreated
soil. The value of 𝑁𝑐 is highly dependent on the compressibility of the soil surrounding the
column, where it is increased with the increase of soil stiffness (Barksdale and Bachus, 1983;
McKelvey, 2002).
To conclude, determination of load bearing capacity of a single stone column is complex and
there is no exact mathematical solution to predict it. This is due to the uncertainty of the
interaction behaviour between the stone column and the surrounding soils. However, the
relationship based on the laboratory model tests proposed by Hughes and Withers (1974) is
38
Samir Ashour Chapter 2: Literature Review
widely considered realistic for vibro stone columns design (Greenwood, 1991; Babu et al.,
Generally, stone columns are created using high friction granular material where the stiffness
of the column is dependent on the lateral support given by the surrounded soil. If this support
is not sufficient the column will fail (Barksdale and Bachus, 1983; Bergado et al., 1991).
The failure mechanism of a single column can be a function of two main factors in addition to
its length whether it is considered to be long or short (i.e. the length of the column is greater
𝐿
or shorter than the critical length(𝐷 ≈ 6)).
The first one is the method of column construction as either end bearing on a firm competent
The second factor is the load application on the column, which could be either directly to the
top of the column or through a rigid cap over an area larger than the diameter of the stone
column (Barksdale and Bachus, 1983; Sivakumar et al., 2004a). Figures 2.12 shows the
39
Samir Ashour Chapter 2: Literature Review
Rigid foundation
Side Friction
2-3D
End Bearing
D
a. Long stone column with rigid base, loaded b. Long stone column with rigid base c. Short column with rigid base d. Short floating column
to area larger than the column diameter Bulging failure Shear failure Punching failure
Bulging failure
Figure 2.12 Failure mechanism of single stone column (Barksdale and Bachus, 1983)
Bulging failure: for columns overlying on rigid base, having length over their critical
length and loaded directly over the column area, a bulging failure would occur over 2
to 3 diameter of the depth. However, if the load was applied over larger area, then
there would be an increase in the vertical and lateral stresses as in the surrounding soft
soil, which would then affect the response of the column to loading. Normally this
leads to smaller bulging and higher ultimate bearing capacity (Barksdale and Bachus,
1983).
𝐿
Shear and punching failure: columns shorter than the critical length (𝐷 ≈ 6) were
likely to fail in shear failure if they were end bearing on rigid base, or in punching
failure if they were floating columns (Barksdale and Bachus, 1983).
Bulging of the column is more noticeable in the upper portion of the column as indicated by
Greenwood (1970), Hughes and Withers (1974), Barksdale and Bachus (1983), Charles and
Watts (1983), Greenwood (1991), and Sivakumar et al. (2004a). Barksdale and Bachus (1983)
suggested that bulging would occur within 2 to 3 times the column diameter whereas Hughes
and Withers (1974) and Sivakumar et al. (2004a) observed that the column bulged at depth of
approximately four times the diameter of the column. Bae et al. (2002), suggested that the
column diameter could be the main parameter affecting the depth of this bulging zone, while
40
Samir Ashour Chapter 2: Literature Review
other factors such as soil strength and the depth ratio could be considered to have less
influence. on the other hand Sivakumar et al. (2007) noted that the degree of bulging is
The mode of failure of stone columns in group has been investigated by a number of
researchers (e.g. Barksdale and Bachus (1983); Hu (1995); McKelvey (2002); Sivakumar et
al. (2004a)). In general they indicated that stone columns constructed in groups showed
different failure behaviour than the isolated columns; in groups, each column can interact and
restrain the expansion of the neighbouring column leading to increase the bearing capacity.
Hu (1995) studied the behaviour of stone column groups and suggested that there are three
different modes of failures (bulging, shearing and lateral deflection) which depend on the
geometric configuration of the columns. Different geometry parameters (L/r0 = 2, 3.2 and 3.4;
rc = 5.5 and 8.75 mm; As = 24 and 30 %) were tested in group of 5 and 7 sand column
reinforced by a footing with a radius (r0) of 50 mm. Hu (1995) showed that the shear planes
through the columns would occur towards the edge of the footing, while bulging could be
seen deep directly under the foundation. The depth of the bulging was found to increase as the
area replacement ratio As increases. Additionally he observed that short columns tend to
Sivakumar et al. (2004a) noted that the central columns beneath the foundation bulge
uniformly, whereas the others at the edge they bulge away toward the surrounding soils.
A range of methods can be used to examine the deformation and failure patterns of stone
columns. For example, Hughes and Withers (1974) used an X ray technique to monitor the
deformation of isolated stone column. Whereas (Hu (1995), Wood et al. (2000), Ambily and
41
Samir Ashour Chapter 2: Literature Review
Gandhi (2007)) successfully formed a plaster-cast of the vacuumed holes after exhuming the
columns material.
McKelvey (2002) and Sivakumar et al. (2004a) used a transparent medium as can be seen in
Figure 2.13 (with shear strength properties similar to the soft clay) allowing constant
Figure 2.13 Photographs of sand columns beneath circular footing at beginning, middle and end
of foundation loading process: (a) TS-01, 150 mm; (b) TS-02, 250 mm
42
Samir Ashour Chapter 2: Literature Review
Sivakumar et al. (2004b) investigated the behaviour of single stone column under the
condition of undrained triaxial, they examined the mode of failure by splitting the specimen
Another method that can be used for this purpose by grouting the stone column by pouring the
concentrated cement slurry after testing and allowing it to set for some time depending on the
concentration of the slurry and then carefully removing the surrounding soils (Sivakumar et
al., 2010).
Hughes et al. (1976) defined the critical length of the column as the shortest length that allows
the column to carry the designed ultimate load regardless of the deformation that would occur.
Where, 𝑃𝑢𝑙𝑡 is the ultimate column load; 𝑐̅ 𝑎𝑛𝑑 𝑐 are respectively the average shaft cohesion
and the cohesion of the soil at the bottom of the critical length; 𝐴𝑐 is the column area with
diameter D; and (𝜋𝐷𝐿𝑐 ) representing the surface area of the side of the column with critical
length 𝐿𝑐 .
Hughes and Withers (1974) introduced the L/D ratio and they found that there is no
improvement in the bearing capacity of the column beyond (L/D = 4.1). Samadhiya et al.
(2008) studied the effect of the column length on the column bearing capacity, they identified
that the critical length was at 4.5D and found that any increase of the bearing capacity was
marginal after this length. Black et al. (2011) presented results obtained from a large physical
triaxial model, they investigated the influence of L/D ratio and the area replacement ratio on
the settlement improvement factor. Black et al. (2011) suggested that the L/D ratio greater
43
Samir Ashour Chapter 2: Literature Review
than 6 is required to provide the full limiting axial stress on the column. in addition they
found that an L/D ratio beyond 8 offers insignificant improvement in terms of settlement
control (i.e less than 5%). They also suggested that settlement can be controlled in short
columns (L/D < 6) using a relatively high area replacement ratio (> 30 %).
In conclusion, column length larger than the critical length may not provide further
improvement in the column bearing capacity, but it may be used to control the level of
settlement
As explained in previous sections (2.7.2 and 2.7.3), surrounding soils provide the lateral
confinement to stone columns, therefore when designing stone column it is important to have
an appropriate knowledge of the soils undrained shear strength. Figure 2.14 shows the
influence of this factor on the ultimate carrying capacity of the stone column according to
both an empirical design relation presented in Equation 2.16 and different laboratory
investigation results carried out by Hughes and Withers (1974), Ambily and Gandhi (2004),
Kim and Lee (2005), Ambily and Gandhi (2007), Zahmatkesh and Choobbasti (2010) , Black
et al. (2011), and Ali et al. (2014). Although the obtained results did not show a perfect fit
with the prediction lines, they show that the column bearing capacity increases with the
increase of the soil shear strength. However, the degree of improvement at similar soil
strength is dependent on other factors such as the area replacement ratio and the internal angle
of friction. For example, Ambily and Gandhi (2007) investigated the effect of the soil shear
strength cu= 7, 14, and 30 kPa) and the area replacement ratio (As = 5, 10, and 19 %) on the
behaviour of single stone column, they noted that the ultimate bearing capacity of the
44
Samir Ashour Chapter 2: Literature Review
improved soil increased in linearly from 175 kPa to 740 kPa when increasing the undrained
600
400
qult Prediction ∅′ 35o
200
0
0 5 10 15 20 25 30 35 40
Undrained shear strength, (kPa)
Figure 2.14 Effect of undrained shear strength of the surrounding soil on the ultimate bearing
Data from previous research studies Wood et al. (2000), Kim and Lee (2005), Ambily and
Gandhi (2007), Najjar et al. (2010), Zahmatkesh and Choobbasti (2010), Black et al. (2011),
were collected and plotted together in Figure 2.15 in order to observe the impact of the area
replacement factor on the bearing capacity of the improved ground. From Figure 2.15 it can
be seen that as spacing between the columns increases, the area replacement ratio decreases,
leading to a decrease in the axial capacity of the column which may lead to increase degree of
settlement. An area replacement ratio from 5 to 25 % could improve the column capacity by
45
Samir Ashour Chapter 2: Literature Review
occurs when the area replacement ratio between 0 and 20 % and beyond As of 20 % the
1000
Ultimate bearing capacity of the treated
cu = 30 kPa
800 cu=35 kPa
Kim and Lee (2005); cu=4kPa
Zahmatkesh and Choobbasti (2012); cu=5kPa
ground, (kPa)
cu = 5 kPa
200 cu = 4kPa
Figure 2.15 Effect area replacement ratio on the ultimate bearing capacity of the treated ground
The shear strength of the surrounding soil has a large impact on the improved ground bearing
capacity, where at As of 20 %, for example, the column capacity was increased from 170 to
For significant improvement in bearing capacity, Wood et al. (2000) indicated that the area
replacement ratio should be 25% or above, and Ambily and Gandhi (2007) suggested that the
Black et al. (2011) examined the influence of area replacement ratio and the column length on
the performance of a footing supported on stone column using a developed large triaxial cell.
The cell accommodated a consolidated specimen (cu =35 kPa) of 300 mm diameter by 400
mm height. Granular columns with different diameters of 25, 32, and 38 mm and various
lengths 125, 250, and 400 mm were tested. Results indicated that the settlement improvement
factor any increased with the increase in the area replacement ratio; however, the
46
Samir Ashour Chapter 2: Literature Review
improvement appeared to have a threshold of about 30–40 % area replacement ratio. For
similar to the individual columns at the same area replacement ratio. Therein settlement
improvement factors for the area replacement ratios of 28 and 40 % were 3.2 and 3.8
respectively. The settlement improvement factors for the corresponding single columns were
about 6.5–7.5 indicating that the performance of the group is not as good as that of the single
columns. Limited results for pressures recorded in the column and in the clay in the group
Bergado et al. (1987) conducted a full scale investigation on six footings by applying load on
single rammed aggregate piers (crushed gravel) with a diameter of 0.3 m and installed to a
layer of very soft clay with 6 meter thickness. The undrained shear strength ranged from 30 to
40 kPa for the upper clay layer, and from 15 to 25 kPa for the soft clay. Columns were loaded
in maintained stress increments with footings with diameters of 0.3, 0.45, 0.60, 0.75, 0.9, and
1.2 m, representing replacement ratios of 100, 44, 25, 16, 11, and 6 %, respectively. Results
indicated that the settlement improvement factor decreased from 8 to 1.5 as the area
Dynamic loading applications on the subsurface of soft soil layers induced by vibration
sources (e.g. earthquakes, traffic loads, or offshore waves) can cause a complex stress field, in
which these stresses vary with time (Gu et al., 2012). The stress strain behaviour of clay soils
47
Samir Ashour Chapter 2: Literature Review
under cyclic loading is dependent on factors such as stress level, drainage condition, type and
Undrained long-term cyclic loading for normally consolidated clays may lead to cyclic
failure. This failure can be identified by number of loading applications in which an arbitrary
predetermined double amplitude failure axial strain is reached (Andersen et al., 1980; Jianhua
Wang et al., 2006; Li-Zhong Wang et al., 2011). Yasuhara et al. (1992) considered residual
Loading frequency is also considered an important issue in cyclic loading. For example,
Ishihara (1996) classified the dynamic problems according to the time of loading as shown in
Figure 2.16. in other wards instances where a load application stays for more than tens of
seconds are considered as static whereas those with a shorter time of loading are considered as
dynamic loading. For example, a shaking wave during earthquakes involves 10 to 20 times
repetition of loads with different amplitudes, and the period of each impulse ranges between
0.1 and 3.0 seconds, which gives a time of loading of 0.02 to 1.0 seconds. On the other hand,
in case of repetitive loads induced by traffic, the soils in the subgrade underneath railways or
road embankments are subjected to a large number of load cycles during the life span. The
time of loading may be deemed on the order of 0.1 second to a few seconds.
48
Samir Ashour Chapter 2: Literature Review
The influence of loading frequency on soft clay behaviour is not well understood. Some
published studies show that both the accumulated pore pressure and shear strain induced by
cyclic loading increases with the reduction in loading frequency (Matsui et al., 1980; Procter
and Khaffaf, 1984; Wang et al., 1998). However, Ansal and Erken, (1989) and Hyde et al.,
(1993) reported different results, indicating that frequency has negligible influence on cyclic
Stone columns have been used to provide the efficient support for different type of
foundations in different soils. They have been used to support infrastructure projects, where
dynamic loading is imposed such as railway and road embankments. Stone columns also have
been used to mitigate liquefaction induced by earthquakes, by increasing the density of the
surrounding soil, allowing drainage, which controls the level of pore water pressure under the
foundation. Also the presents of the stone column will increase the total load carrying
49
Samir Ashour Chapter 2: Literature Review
capacity and at the same time the stress level on the surrounding soil will decrease (Munfakh,
1984; 1997; Van Impe et al., 1997; Raju, 2003; Adalier and Elgamal, 2004).
Liquefaction mitigation
Liquefaction occurs in loose to medium dense, saturated soils with fairly uniform grain size
distribution, covering the silty sandy range. In this soil condition, the dynamic forces imposed
by earthquake lead to a rearrangement of the grain structure to a denser state. Therefore, if the
drainage is insufficient, pore water pressure will increase and shear resistance of the soil will
Priebe (1998) presented a method for evaluating the potential impact of liquefaction with
stone columns. This method is similar to some extent to the design of vibro replacement in
In the dynamic condition, in order to be more realistic and to simplify the problem, Priebe
(1998) considers deformation of soil with the constant volume to calculate with Poisson’s
ratio µs = 0.5. In this method, the improvement factor n0 is determined using some
𝐴𝑐 1
𝑛0 = 1 + ∗[ 𝐴 − 1] 3.14
𝐴 𝐾𝑎𝑐 ∗(1− 𝑐⁄𝐴)
∅′𝑐⁄
𝐾𝑎𝑐 = 𝑡𝑎𝑛2 (45𝑜 − 2) 3.15
Where: A is the area within the compaction grid; Ac is cross section of stone column; and ∅′𝑐
50
Samir Ashour Chapter 2: Literature Review
Baez and Martin (1992) reviewed experimental data on the use of stone column for
valid only if the pore pressure ratio remains below 0.5. In addition, they noted that pore
pressure within the gravel drain is not constant, contrary to the assumption that they remain
essentially unchanged.
A case study presented by Rudolph et al. (2011) looked at the effectiveness of using rammed
aggregate piers (RAP) in mitigating the liquefaction potential in a site containing artificial
fills over liquefiable sandy clay soils; and a groundwater level was approximately 3.4 meters
below ground surface. RAP used had a 0.40 meter diameter, 8.5 meters depth and spaced at
2.1 meters centre to centre. The study focused on the densification of the soil surrounding the
RAPs. pre- and post-RAP CPTs measurements were used to calculate the Liquefaction
Potential Index (LPI) of the matrix soil and residual liquefaction and seismic settlement
51
Samir Ashour Chapter 2: Literature Review
Their results indicated that using rammed aggregate piers improved the soil properties of the
site and has been effective at mitigating liquefaction potential to acceptable levels, for
example, within the spread footing area, the reduction in LPI from the pre-RAP CPT versus
post-RAP CPT shows a reduction from very high risk (LPI of approximately 17.8) to low risk
Rollins et al. (2006) presented a case history in which prefabricated vertical drains and stone
columns were used in combination to provide treatment for a 4 meters thick layer of
liquefiable silty and sandy silt. They suggested that this method can only be suitable for the
soils with fines percentage lower than 20%, otherwise less effectiveness of this method will
be achieved. Also they reported that using stone column together with wick drain will
Embankments:
Building embankments on soft clay soils is usually a very challenging geotechnical task due
to the potential of bearing failure, excessive settlement, and local and global instability under
dynamic load (Mitchell and Jardine, 2002). Vibro stone column is one of the most commonly
adopted soil improvement method, it has been utilized to increase the bearing capacity and
Kempfert et al., (1999) investigated the bearing and deformation behaviours of a geotextile
encased sand column foundation at area ratios of 4; 12 and 16 % under static and cyclic using
both small scale models with scales of 1:6.5; 1:11.5 and 1:13 and full scale models with scale
of 1 : 1. The cyclic loading was divided into two stages as shown in Figure 2.18, the first
sequence was simulating an embankment with low height (i.e low preloading pressure (25
52
Samir Ashour Chapter 2: Literature Review
kPa) and high cyclic stresses expected on the soil (62.4 kPa)), whereas in the second cyclic
loading sequence a deep embankment was simulated by applying a low cyclic stresses with an
amplitude of 30.5 kPa and high preloading pressure 100 kPa). At the end of the cyclic loading
Results showed that at the same level of stresses the system developed a higher settlement
during the cyclic loading compared with the static; also they indicated that during the first
loading sequence the settlement was developed very quickly in comparison with the second
loading sequence. They showed that at a constant Area ratio of 12 %, 100 kN/m2 of loading
stress and a scale model of 1:6.5, using a geotextile coated sand columns can reduce
Figure 2.18 Cyclic loading sequences for railway embankment (Kempfert et al., 1990).
Kolekar et al. (2011) presented results from an experimental study investigating the behaviour
of stone columns under repeated loading. The cyclic test consists of two stages, in stage one
35 % of the static load failure was applied in a form of sinusoidal wave at a frequency of 0.1
Hz for a number of cycles of 500, whereas in the second stage the load was increased to about
53
Samir Ashour Chapter 2: Literature Review
100 % of the load at the static failure. The study was focused on the vertical stiffness ratio.
Kolekar et al. (2011) claimed that the stiffness ratio was increased dramatically (from 7 to 80
kPa/mm) by moving from stage one to stage two loading condition and suggest that the first
loading stage compacted the column making it more dense, however, they defined the
stiffness ratio as the ratio between the change in the vertical stress to the average vertical
settlement and as the deformation is much smaller in stage one compared with that in stage
two the material stiffness should be larger than that in the second stage.
Fatahi et al. (2012) employed a finite element model using PLAXIS to assess the relation between
the column position beneath the train track and overall settlement of the ballast rail formation. The
study also includes the effect of using geogrid reinforcement (where two layers of geogrid were placed
at the interface between the subgrade and sub-ballast; and sub-ballast and ballast). The model
geometry considers the typical track cross section with concrete sleepers as recommended on the NSW
rail network (i.e. layers depth from top to bottom including sleeper, ballast, sub-ballast, are 150, 300
and 150 mm). The subgrade layer depth was assumed to be 10 m and the gauge length of the track
was 1.4 m. Also train load was considered as 125 kN/m. Stone columns (1 m diameter) were
arranged in rectangular grid pattern at with a spacing of 1.5 m along the rail track.
Results indicated that the overall settlement was reduced by the presence of stone columns spaced
more closely at the centre of the track and not just under the rail (i.e. when the offset distance of the
columns from the track centreline increases, the maximum horizontal displacement of the improved
sub-grade also increases). Also indicated that the optimum column pattern will occur, when stone
columns just overlap in the cross section and located right in the centreline of the track.
In addition, Fatahi et al. (2012) showed that the use of stone columns (no geogrids reinforcement) was
less effective in reducing vertical settlements than geogrids only (no column), when distance of
columns from the centre line exceeds 1 m. Furthermore, the use of two stone columns with geogrids
54
Samir Ashour Chapter 2: Literature Review
has limited benefits over the use of just geogrids when the distance of columns from the track
Soft clay deposits usually have a low bearing capacity, low permeability, and high
compressibility. Thus it is important that these soils receive sufficient improvement before
construction activities commence. Whilst there are a range of ground improvement techniques
available, the use of stone columns is still regarded as one of the most popular and effective
increasing its shear strength, reducing excessive and differential settlement whilst speeding up
the consolidation progress by shortening horizontal drainage paths. Thus, it has been
successfully adopted for projects such as highway embankment, industrial and residential
structures.
Over the last three decades a considerable amount of research has been conducted to
investigate the behaviour of stone column in soft soils. This chapter reviewed different
aspects (i.e. critical column length, area replacement ratio, surrounding soil strength, internal
angle of friction of the column material, different column material, and column mode of
failure) that have been examined by authors to identify their influence on the behaviour of the
stone column. In addition, this chapter discussed design approaches and construction methods
that are commonly adopted. Almost all the previous study parameters mentioned above were
based on the assumptions that the external loading is applied monotonically with a constant
strain rate during the whole duration of the test. However, in many engineering projects, such
as tank filling and discharging, highway embankments, ocean banks, etc., the surcharge
loading is not applied instantaneously, but changes with time. Therefore, it must be
55
Samir Ashour Chapter 2: Literature Review
acknowledge that the exact behaviour of foundation on unimproved soft soil and on stone
high demands for high speed railway upgrade became a topic of discussion. Modern railway
infrastructure demands have increased not only to provide a high level of performance in
terms of settlements and stability of the railway track but also to reduce maintenance cost. In
areas where loose or soft cohesive deposits are found, ground improvement methods are often
required to ensure the required level of performance. Therefore, is it still viable to consider
Bulging or stone column expanding under applied load is important to understand the general
behaviour of vibro stone column foundation, whilst some researchers have studied this topic
under the application of static loading, the question here is would a column bulge in the same
way if the loading sequence was changed (i.e. cyclic loading)? In addition to this, what are the
implications for the threshold stress such improved soils can carry when compared to the
static stress failure? And what is the influence of the loading frequency on the overall
The current research was undertaken in an attempt to investigate the above questions and to
provide a better understanding to the behaviour of stone column under cyclic loading
application.
56
Samir Ashour Chapter 3: Methodology
3 CHAPTER 3
3. METHODOLOGY
3.1 Introduction
As noted in the previous chapter, laboratory investigations have played an important part in
trough full scale trials. Thus, mindful of the cost of full scale trials; and in order to further
advance the development of stone columns for cyclic loading a laboratory investigation was
undertaken. This study includes investigation of behaviour of stone columns under both static
Properties of material used together with both the description of apparatus used and procedure
are described in this chapter. It also includes a programme of the laboratory investigation
undertaken.
3.2 Materials
This laboratory scale investigation focused mainly on the fundamental mechanical behaviour
of stone column foundations subjected to cyclic loading. Twenty two laboratory based studies
(see Table 3.1) were reviewed in order to identify best practice and to identify suitable
materials for both the stone column construction and that for the surrounding soil. It was
found that over 50% of the studies were based on Kaolin clay bed because of its high rate of
consolidation compared to that of other fine soils derived from natural deposits. It could also
be used to produce high quality repeatable samples as Kaolin is a processed clay. In terms of
stone column construction, crushed rocks were found to give the highest angle of internal
57
Samir Ashour Chapter 3: Methodology
friction (over 40o) and were used by about 50 % of the researchers. Others used sand and river
gravel. Therefore a crushed basalt and kaolin clay were chosen to simulate the column and the
Surrounding
Column material φ' (o) cu (kPa) Reference
material
1 Leighton buzzard sand Kaolin clay 35 19 Hughes and Withers (1974)
2 River gravel Grey silty clay 38 22 Hughes et al. (1975)
3 Gravel Boulder clay - 4.4 Greenwood (1975)
4 Uniform gravel Boulder clay 47, 51, 53 30 Charles and Watts (1983)
5 Quartz Kaolin clay - 14.4-19.1 Barksdal and Bachus (1984)
6 Sand/gravel Bangkok marine clay 35, 43 20-45 Bergado and Lam (1987)
7 River sand Silt 38 - Juran and Guermazi (1988)
8 Loch Aline fine sand Kaolin clay 30 - Hu (1995).
9 River sand - 40.5 - Rajagopal et al. (1999)
10 Stone Ash/ clay 45 40 Watts et al. (2000)
11 - Kaolin clay 35 28 Sivakumar et al. (2004)
Crushed rock, Quarry Transparent clay-type
12 51, 46, 39 17.5-21.5, 32 McKelvey (2002)
waste, Crushed concrete material, Kaolin clay
13 Coarse sand Sand/clay 44 - Ayadat and Hanna (2005)
14 Sand/gravel CH 38, 49 12 Kim and Lee (2005)
15 Crushed basalt Kaolin clay 35 Black et al. (2006)
16 Crushed limestone CL clay 43 30 White et al. (2007)
17 stones CH 43 7, 14, 30 Ambily and Gandhi (2007)
18 Granite chips Lake clay 41 2.5 Murugesan and Rajagopal (2008)
19 Commercial sand Kaolin clay 35 5 Gniel and Bouazza (2009)
20 Quartz Kaolin clay 33 - Najjar et al. (2010)
21 Crushed basalt Kaolin clay 35 - Black et al. (2011)
22 Crushed basalt Kaolin clay 35 - Sivakumar et al. (2011)
23 Gravel Kaolin clay - 41 Cimentada et al. (2011)
24 Gravely sand Oxford clay - 16 Madun et al. (2012)
25 Aggregate Soft clay 41 5 Vekli et al. (2012)
26 Sand Clay 38-45 10.5-16.5 Hanna et al. (2013)
27 Hostun sand and Gravel Kaolin clay - - Frikha et al. (2014)
28 Stone chips Kaolin clay - 6-6.8 Ali et al. (2014)
29 Aggregate Pulverized clay 45 54, 15 Mohanty and Samanta (2015)
58
Samir Ashour Chapter 3: Methodology
Physical properties of Kaolin clay were determined in accordance with British Standard BS
1377 (1990) (The tests data available in Appendix A). These properties are summarised in
Table 3.2.
3.3 Physical Properties of the Model Stone Column and Surrounding Soil Material
Both the Liquid and Plastic limits for the Kaolin clay were determined according to the
procedures in Clauses 4.3 and 5.0 of BS1377: Part 2 (1990) (BSI, 1990). Liquid Limit of
Kaolin was 56%, and the Plastic Limit was 27%. These values are considered typical for
The particle size distribution of the Kaolin clay material was determined, using a hydrometer
in accordance with Clause 9.5 of BS1377: Part 2 (1990) (BSI, 1990). Results showed that
about of 50 % of the material was in the fine silt range and 50 % was clay (sub 0.002 mm).
The particle size distribution of the crushed rock was determined by using a dry sieve analysis
method in accordance with Clause 9.3 of BS1377: Part 2 (1990) (BSI, 1990), all the material
59
Samir Ashour Chapter 3: Methodology
used was between 1.18 and 2 mm. Typical grading curves for both the Kaolin clay and the
100
Crushed Basalt
Percentage of passing (%)
80 Kaolin clay
60
40
20
0
0.001 0.01 0.1 1 10
Particale size (mm)
Figure 3.1 Typical particle size distribution of Kaolin clay and crushed basalt
Specific gravity was determined in accordance with Clause 8.3 of BS1377: Part 2 (1990)
(BSI, 1990). The specific gravity of the Kaolin clay ranged from 2.61 to 2.65 with an average
of 2.63, and the values for crushed rock used to form the column, ranged from 2.70 to 2.73
with an average of 2.71. These values fall within the range typical of these material types.
The compaction test for Kaolin clay was conducted using a 2.5 kg rammer compaction falling
from 300mm (light compaction method). The test was carried according to BS 1377: Part 4
(BSI, 1990). The dry - moisture content relationship is given in Figure 3.2. The maximum dry
density was estimated as 1410 kg/m3 which occurred at an optimum moisture content of 29%.
60
Samir Ashour Chapter 3: Methodology
The maximum density of the crushed aggregate was determined in accordance with the
in a mould under a controlled water level using a vibrating hammer. This procedure was
repeated three times to ensure the results repeatability. The maximum density was found to be
The minimum density was determined by pouring the material from 0.5 m height into the
mould without disturbing the soil (BSI, 1990). The minimum density for the crushed basalt
1.50
1.40
Dry Density (Mg/m3)
1.35
1.30
0%
1.25
5%
1.20
10 %
1.15
1.10
10 20 30 40 50
Moisture Content (%)
In order to simulate very soft ground that has shear strength below 15 kPa, a series of quick
undrained triaxial tests and hand vane shear tests were conducted on compacted Kaolin clay,
at different water contents ranging from 30 % to 48 %. Shear strength of soil ranged between
61
Samir Ashour Chapter 3: Methodology
53 kPa and 5 kPa for moisture content of 33 % to 48 % respectively. Correlation between the
water content and undrained shear strength is shown in Figure 3.3. From these test result, it
can be concluded that soil shear strength of 15 kPa can be achieved with a water content of
42%.
60
40
y = 9569.1e-0.159x
30
20
10
0
30.00 35.00 40.00 45.00 50.00
Figure 3.3 Relationship between undrained shear strength and water content for Kaolin clay
Undrained shear strength can be related to liquidity index of soil (Skempton, 1957), as shown
𝐿𝐼 ≅ 1 − 0.192 ln 𝑐𝑢 3.1
Vardanega and Haigh (2014) examined a database of 641 fall cone tests in different soil types
in order to determine a statistical relationship between the undrained shear strength and the
liquidity index. Although they used a different method (fall cone tests) to carry out the
undrained shear strength of the soil, a considerable agreement between their statistical model
(Eq. 3.2) and the one presented in the current study was found.
62
Samir Ashour Chapter 3: Methodology
0.8
y=-0.226ln(cu) + 1.12
0.7 R2 = 0.949
0.6
Liquidity index, LI
0.3
Quick undrained triaxial
0.2 Vane shear test
(Varadanega and Haigh, 2014)
0.1
0
1 10 100
Undrained shear strength, cu, (kPa)
A series of direct shear tests were carried out on dry and wet samples of crushed basalt. In
order to simulate the condition where stone columns are installed in wet ground, the crushed
basalt was left to soak in water for 24 hr before testing. Each test was repeated three times to
ensure result repeatability. All tests were performed in a 100 mm x 100 mm direct shear box,
under vertical stresses of 25 kPa, 50 kPa, 100 kPa and 200 kPa. These normal stresses were
chosen to represent the typical application in the field. The material was placed into the shear
box in three layers; each one was compacted using a small wooden block. The average
density of the compacted material was 1540.6 kg/m3, with variation of about ± 20 kg/m3.
Shearing rate of 1 mm/min was used for all the tests, as suggested in the BS 1377: Part 7
(BSI, 1990), for the aggregate material. During the test, horizontal; vertical displacements and
the shear force measurements were recorded. Typical test results are presented in Figure 3.5.
63
Samir Ashour Chapter 3: Methodology
The general behaviour of the crushed basalt was followed a typical pattern for dense granular
The peak value of the shear stress was plotted against the corresponding normal stress in order
to determine the failure envelop (see Figure 3.6). for example at 100 kPa effective pressure,
under dry condition, the angle of internal friction was found to be 50o, whilst in wet condition
300 300
(a) Wet condition (b) Dry condition
250 200 kPa 250 200 kPa
Shear stress, (kPa)
100 100
50 kPa 50 kPa
50 50
25 kPa 25 kPa
0 0
0 5 10 0 5 10
Horizontal displacement, (mm) Horizontal displacement, (mm)
4.5 4.5
4 (b) Wet condition 4 (b) Dry condition 25 kPa
25 kPa
3 3 100 kPa
2.5 50 kPa 2.5
2 2
1.5 1.5
100 kP
1 200 kPa 1
0.5 0.5
0 0
0 5 10 15 0 5 10 15
-0.5 -0.5
Horizontal displacement, (mm) Horizontal displacement, (mm)
Figure 3.5 Typical stress-strain, direct shear box test (a) wet condition and (b) dry condition
64
Samir Ashour Chapter 3: Methodology
300
Wet condition
250
Dry condition
200
Shear stress; kPa
150
100
50
0
0 50 100 150 200 250
Normal stress, (kPa)
In order to determine the primary consolidation behaviour of Kaolin clay within tests, the
consolidation time, settlement and required load were needed. One dimensional consolidation
tests were carried out using the oedometer apparatus. A small quantity of slurry (1.5 * L.L
water content (Head, 1996; Sivakumar et al., 2004b)) was prepared then placed in the 75 mm
diameter consolidation ring. Considering the condition of the soil (slurry), the specimens were
firstly consolidated under the load of the apparatus hunger only for 24 hour, and then they
were monotonically loaded to pressures of 25 kPa, 50 kPa, 100 kPa, 200 kPa and 400 kPa. At
this stage the samples were unloaded to 50 kPa then reloaded to 800 kPa. Each pressure
maintained for a period of 24 h. During the increment loading period, the settlement of the
sample was carefully measured and the coefficient of consolidation, cv, was determined using
Equation 3.3:
𝑇𝑣 𝑑 2
𝑐𝑣 = 𝑡90
3.3
65
Samir Ashour Chapter 3: Methodology
The square root time method was used to determine t90 (Head, 1998). A typical result for the
The relationship between the void ratio and applied consolidation stresses is presented in
Figure 3.8. The void ratio, e, at the end of each increment period was calculated using
Equation 3.4:
∆e 1 + e0
= 3. 4
∆H H0
e0 is initial void ratio and can be found using the relationship between the water content and
Where, w is the water content, and Gs is the specific gravity of the clay.
The coefficient of compressibility, mv, at each pressure increment was determined using the
following equation:
1 𝑒 −𝑒
𝑚𝑣 = ( 0 1)
1+𝑒0 𝜎1 − 𝜎0
3.6
66
Samir Ashour Chapter 3: Methodology
In Equation 3.6, e0 and e1 are the initial and final void ratio respectively; and σ0, and σ1 are the
The coefficient of consolidation, cv, ranged from 7.8 to 27.2 m2/minute, whereas the
coefficient of compressibility, mv, ranged between 0.51 and 3.62 (m2/MN). Both of these
coefficient values are comparable with other typical values of normally consolidated clays
(Head, 1998). The compressibility and consolidation coefficients for each loading increment
0 5 10 15 20 25 30 35 40
0.00
0.10
Settlement, (mm)
0.20
0.30
0.40
0.50
0.60
0.70
Figure 3.7 Consolidation curve square root of time step 50 to 100 kPa
67
Samir Ashour Chapter 3: Methodology
1.3
Cc= 0.275,
1.2 Cr= 0.069
1.1
1.0
void ratio, e , (%)
0.9
0.8
0.7
10 100 1000
3 30
2 25
(m2/MN)
mv, (m2/MN)
mv,
cv, (mm2/minute)
cv, (mm2/minute)
1 20
0 15
10 100 1000
Δσ́, (kPa)
68
Samir Ashour Chapter 3: Methodology
laboratory scale test. However, it was not possible in physical model tests undertaken to
maintain a suitable scale factor that satisfies completely all the parameter that governs the
maintain the prototype stress level (Wood, 2004); however, this approach was unavailable for
Hu (1995) noted that parameters affecting response of soil-stone column systems and
influencing the load and settlement relationship can be divided into two categories of major
q ult G G γ d γ
⁄cu = f1 (S⁄D , L⁄D , L⁄d , c⁄G , As , ∅′ ) . f2 ( s⁄cu , s⁄γc , g⁄d , DB , c⁄cu ) 3.7
s c
Where:
S: the spacing between the column and represented by the area replacement ratio As
d
As = (mS)2 3.8
69
Samir Ashour Chapter 3: Methodology
According to this, some of the main aspects of the field situation such as L/D, L/d, As, and ϕ’
were maintained in the laboratory scale testing. It was considered that this could provide an
adequate response and would generate high enough quality data so that it can be compared
with other studies of similar conditions and would suggest expected behaviour from full scale
columns.
On the other hand, adopting the unit cell concept, and using the relationship between area of
the column to the area of the treated soil then refer it to the diameter ratio (N) as in equation
3.9.
2
A d
As = col⁄A = [ c⁄d ] = 1⁄ 2 3.9
e N
Where: Acol and A is the area of the column and the soft soil area, respectively. Typical value
of this ratio is 5 to 30%. This range results in approximate values of diameter ratio (N)
between 2 and 5. This shows an agreement with (Hughes and Withers, 1974) indicated that
the column can improve the strength of the surrounding soil that falls within a diameter of 2.5
times the diameter of the column. Also Bowels (1988) identified the failure zone extended
radially to a distance of 1.5 times the diameter of the foundation and over a depth of
To satisfy the above conditions and taking account of the facilities available, the maximum
size of 100 mm diameter foundation could be used. Therefore, N could be used as 3.5. This
leads to the diameter of the column of 28 mm. In terms of area replacement, these values
70
Samir Ashour Chapter 3: Methodology
correspond to replacement rates of 7%. With this, assumption the horizontal scale of the
The diameter of the stone column (28 mm) and the spacing (100 mm) were selected on the
The length of sample was determined by taking account of the typical value of length/ column
diameter ratio (L/d), which is generally between 6 and 20. This ratio may reduce to 4 in case
After selecting the diameter and length of stone column, it is important to select an
appropriate particle size of material to be used. Typical ratios of column diameter to average
particle diameter range from 20 to 40, thus granular materials of 2 mm particle sizes were
adopted. The stone column was subsequently tested in two sizes of clay beds as described
below.
Model I: A single 28 mm diameter stone column was installed in a clay specimen of 100 mm
diameter and 200 mm height and subjected to triaxial loading. This set simulated of constant
pressure surrounding the clay. This approach was first used by Hughes and Withers (1974)
Model II: The test is conducted on a single column of 28 mm diameter constructed in clay
bed of 300 mm diameter and 300 mm height. The column was loaded through a 70 mm
diameter steel plate which was 2.5 times the diameter of a single column. The lateral
dimension of the tank (300 mm) was chosen according to the minimum free distance between
71
Samir Ashour Chapter 3: Methodology
The values of the main parameters relating to columns used in the field and those used in this
study are summarised in Table 3.4. Comparing the typical stone column diameter to the
model diameter implies the scale factor of the model column was about 1/20.
Parameter As L d dg
L/d
Unit (%) (m) (mm) (mm)
The test programme involved application of both static and cyclic loading application, in two
models (3.4.2). Investigating the behaviour of soft soil /stone column system under cyclic
loading requires examining the influence of loading frequencies, in addition to the effect of
the cyclic stress ratios. Therefore, in order to fulfil this aim, the test apparatus and procedures
The clay bed for all tests was prepared by consolidating slurry with an initial moisture
content of 84 %. This allowed for a uniform and homogenous soil sample in which
In order to eliminate the boundary effect of the tank walls on the specimen, a tank
diameter of 300 mm was used (i.e. 4.25 times the model footing diameter (70 mm).
72
Samir Ashour Chapter 3: Methodology
The clay bed was prepared by placing kaolin slurry into five plastic tanks. Two tanks were
designed to prepare 100 mm triaxial specimens. The containers were 100 mm in diameter and
400 mm height and the set up was shown in Figure 3.10. Pressures were applied by adding
weights to a load hanger connected to a piston in top of the slurry in four stages giving
pressures equivalent to 12, 25, 50, and 100 kPa. The piston was perforated to allow water to
Three tanks measuring 300 mm in diameter with 550 mm depth were used to create the larger
specimens of 300 mm diameter by 300 mm height and the set up was shown in Figure 3.11.
For safety reasons the large quantity of weights that would be necessary in order to generate
the required pressure could not be used. Therefore, loading was applied via hydraulic jack
attached to a pump connected to air compressor. Proving rings were fitted between the
hydraulic jack and the perforated plate placed on top of the slurry sample to monitor and
control the amount of loading. Two dial gauges were placed on top of each sample to monitor
and record the settlement of the sample as consolidation progressed, also to recognise a
73
Samir Ashour Chapter 3: Methodology
Dial gage
Frame
Top
drainage
Perforated
Piston
Slurry
Filter
Bench
Bottom
drainage
Weight
Slurry
550 mm
PPT 3
300 mm
400 mm
Ø10mm
PPT1
PPT2 PPT3
50 mm
100 mm
150 mm
300 mm
Figure 3.11 Model II consolidation chamber for preparing 300 mm dia. Samples
74
Samir Ashour Chapter 3: Methodology
A 50 kN capacity load frame (Figure 3.12) was used to apply both static and cyclic loading.
Although it is specified that this machine can provide a range of cyclic loading frequency
between 10 Hz and 0.0001 Hz, several trials were conducted to check capacity of the load
frame. Loading frequencies adopted ranged from 0.5 Hz to 3 Hz. The minimum applied load
Triaxial
Soil specimen
chamber
Data logger
Pores disk
PWP
transeducer
Computer
Soil
Pores disk
Stone column
Filter papers
75
Samir Ashour Chapter 3: Methodology
3.6.2 Instrumentation
Load cell: A 10 kN capacity load cell was used for this study. It was calibrated using
SANS electromechanical universal testing machine (CMT 5000 series). A dynamic
test was run in this cell three times for both compression and tension, in order to
ascertain margins of error, which was found to be less than 1 %.
Data acquisition: A logging system supplied by Servocon Digital Control Ltd was
used. to control and monitor the load and pressure transducers. The signals were read
up to 0.001 accuracy, which was considered acceptable for this research. A logging
interval of 0.1 second was used. This leads to a very large number of data, especially
for the static loading test which lasted up to 10 hours at 0.03 mm/min. For cyclic
loading test of data can was recorded every 10 cycles. The large quantity of data was
considered important as it would help to better define the behaviour of specimen under
loading.
Several factors affect the behaviour of soft subgrade soil subjected to train induced cyclic
loading. these include load frequency, level of applied cyclic stress and stress condition, and
the physical condition of the soil (Brown et al., 1975; Li and Selig, 1998). These are discussed
below.
It has been suggested that loading frequencies and number of loading application have a small
influence on the deformation behaviour of the granular materials (Peacock and Bolton, 1968;
76
Samir Ashour Chapter 3: Methodology
Brown and Hyde, 1975). However, effect can be significant in case of soft cohesive soils
(Brown et al., 1975; Yasuhara et al., 1982; Zhou and Gong, 2001).
Table 3.5 shows the type of soil and test loading frequency applied by other researchers. This
shows that all the cyclic tests were conducted under triaxial conditions with loading
frequencies ranging between of 0.1 to 5 Hz to simulate the loading condition of traffic and
passing trains. The number of load applications ranged between 500 and 1,000,000 depending
on the material.
Table 3.5 Cyclic load testing conditions for previous laboratory investigation
Test condition
Sample
Material Sample σ3’ Frequency Number of Application Reference
size (mm)
preparation (kPa) f, (Hz) cycles
DXL
4 clay 39.1X80 OCR=(1-4) 110-240 0.01- 1 200-3000 - (Zhou and Gong, 2001)
5 Ballast 300X600 compaction 1-240 20 500000 High-speed (Lackenby et al., 2007)
trains
6 Kaolin clay 300X600 anisotropic - 5 3500 PVD (Indraratna et al., 2009)
K0 = 0.6 installation
beneath rail
tracks
7 Marine clay - CSR=1 41 0.01, 0.1, 1 1000 - (Li-Zhong Wang et al.,
2011)
8 Clay 50X100 - 100, 200 0.1, 0.5 400 Traffic loading (Gu et al., 2012)
9 Silt (ML) 76X152 - 100-250 1 100 Traffic loading (Ng et al., 2013)
(25% Kaolin (2 m × 2.2 - - 0.77 100 Subgrade (Sun et al., 2014)
+ 75% m × 2m) stabilisation
Kansas River
sand)
10 Nanjing fine 50X200 - 50-150 0.1 10000 Traffic loading (Cai et al., 2015)
sand
77
Samir Ashour Chapter 3: Methodology
Load frequency of a train, at the subgrade level is related to train speed (𝑣), and carriage
length (𝐿), Sekine (1996) developed a relationship to estimate the loading frequency at the
subgrade level:
𝑉
𝑓 = 3.6𝐿 3.10
In this study, a sinusoidal loading wave form was used to simulate train movement with
loading frequencies of 0.5, 1.0 and 3.0 Hz representing train speed of 35, 70 and 225 km/hr
Based on the depth of ballast and the axle load, there are several theoretical relationships for
determining the pressure at the subgrade level. Perhaps the most well-known is the Talbot
16.8𝑃𝑚
𝑃𝐶 = 3.11
ℎ1.25
where PC is subgrade pressure in psi, Pm is the applied stress on ballast in (psi) and h is the
Japanese National Railway also used similar formula to estimate the stresses at the subgrade
level.
50𝑃
𝑃𝐶 = 10+ℎ𝑚
1.35
3.12
where PC and Pm are as defined in equation 3.11 and h is the depth of ballast in cm.
78
Samir Ashour Chapter 3: Methodology
Field observation of dynamic stresses conducted by Yoo and Selig (1979) using a four axle
hopper cars each weighing 131 tons, with ballast layer thickness of 300 mm, showed that the
maximum stress on the subgrade level ranged between 80 and 110 kPa (Figure 3.13).
Compressive stress
Compressive stress
(kPa)
(psi)
Figure 3.13 Dynamic measurement of subgrade stress (Yoo and Selig, 1979)
However, there are other empirical equations where the speed of the train is also taken it
account together with axle load. Zhou, (1996) derived an empirical formula to predict the
where 𝜎𝑑 =dynamic vertical stress on the subgrade surface (kPa); 𝑃 is the axel load of the
train (kN); 𝑉 is the train speed (km/hr) and α is speed coefficient and equals to 0.005, 0.004
and 0.003 for ordinary grade railways speed, quasi-high speed and high speed railway
respectively.
Furthermore, field measurements in existing railways in China indicates that the average
magnitude of the vertical dynamic stress on the subgrade level for passenger trains is about 46
kPa (Chinese Academy of Railway Science, CARS., 2006). Therefore, for a 225 km/hr train
speed and 100 kN axle load the cyclic stress will be about 43 kPa.
79
Samir Ashour Chapter 3: Methodology
A large scale test at the University of Birmingham (unpublished work) included application of
125 kN wheel load applied at 2 Hz for 2 million cycles on rail across three sleepers (moor
information include in Appendix G) indicated that the stresses subgrade level were
approximately 80 kPa (Figure 3.14) showing good agreement with Yoo and Selig (1979)
study.
100
PT1
PT2
Stresses at the subgrade level, (kPa)
PT3
80
60
40
20
0
0 50 100 150
Time, (minute)
Figure 3.14 Dynamic stresses mesurments at subgrade level (under 300mm of ballast)
Based on the above, cyclic stresses of 50 kPa, 60 kPa and 70 kPa were used in this study.
These values are equivalent to a cyclic stress ratio (CSR) of 0.6, 0.7 and 0.8 (CSR = the cyclic
deviator stress (qcyclic)/ the static deviator stress of reinforced soil at failure (qfailure)).
80
Samir Ashour Chapter 3: Methodology
In order to create homogenous bed of soft clay, Kaolin clay was mixed with water content of
1.5 times the Liquid Limit (i.e. water content of 84 %). A Hobart A-120 mixer with capacity
of 12 litres was used for mixing 5 kg of Kaolin with 4.2 litre of de-aired water for each
required batch. Soil and water was mixed a constant speed of 60 rpm for 60 minutes. At the
middle of the mixing process, in order to check the quality of mixing, the mixer was stopped
and the slurry was checked by hand and freed from any created masses (lumps). Once fully
mixed the slurry mixture was stored in sealed plastic containers for 24 hours before using it to
insure saturation. Each batch was just enough to prepare two samples for triaxial testing. For
the large sample five more batches were required to prepare one sample.
3.8.1.2 Consolidation
approach by many researchers to create homogeneous clay bed (Juran and Guermazi, 1988;
Sivakumar et al., 2004b; Maakaroun et al., 2009). However, whilst this technique can provide
satisfactory results in short consolidation chambers (< 200 mm) such as the Rowe
consolidation cell (Rowe and Barden, 1966), some difficulties have been reported by a
number of researcher (Anderson et al., 1991; McKelvey, 2002; Ahmadi and Robertson, 2004;
Black et al., 2011) when using this method for preparing larger depth samples. This because
friction occurs between the soil and the walls of the consolidation chamber which could lead
clay particles to arrange themselves into an internal concave structure. In addition, there might
be a loss in consolidation pressure with the depth of the samples. (Valls-Marquez, 2009;
81
Samir Ashour Chapter 3: Methodology
Black et al., 2011). To mitigate these problem, as proposed by the previous researchers inner
walls were carefully lubricated with silicon grease which also aids sample extrusion after
consolidation. A hydraulic pump was adopted in order to keep the consolidation pressure
Approximately 4000 g of the prepared slurry was poured into each consolidation chamber to
an initial height of 294 mm in three stages. The cell was vibrated for 30 seconds to eliminate
trapped air.
According to BS 1377-5:1990, under the fully saturation conditions the initial void ratio (e0)
and the initial height (H0) were calculated using equations 3.2 and 3.3, they were 2.209 and
294 mm respectively. Before and after placing the soil, a 100 mm diameter saturated porous
plastic filter (Vyon sheet) with 1.2 mm thickness was placed on top of the slurry surface to
Every load step was held for a minimum of 24 hours and a maximum of 72 hours to ensure
that primary consolidation was complete. Typical consolidation plots are shown in Figures
The void ratio, e, at the end of each increment period was calculated using equation 3.2. The
relationship between void ratio and applied consolidation stresses is presented in Figure 3.17.
After completion of consolidation, the sample was either directly extruded from the
consolidation chamber (no column installed) or the stone column was first installed and then
extruded. A 100 mm diameter Piston extruder was used to remove the samples in the same
direction as consolidation loading. A steel ruler was used to trim the surface of the specimen
to ensure the flatness and height of 200 mm. In each test samples were taken to measure water
82
Samir Ashour Chapter 3: Methodology
content. Filter papers were fitted around the sample and then a rubber membrane was placed
Square root of time √Time, (minute) Square root of time √Time, (minute)
0 10 20 30 40 50 0 10 20 30 40 50
0 0
Settelement, (mm)
12 kPa 25 kPa
Settelement, (mm)
5 5
10
10
15
15
Square root of time √Time, (minute) Square root of time √Time, (minute)
0 10 20 30 40 0 10 20 30 40
0 0
100 kPa
Settlement, (mm)
50 kPa
Settelement, (mm)
5 5
10
10
15
15
Figure 3.15 Typical settlement /√time relationship for 100 mm dia. triaxial sample
83
Samir Ashour Chapter 3: Methodology
0 20 40 60 80 100 120
0
20
Settlement (mm)
40
60
80
100
Figure 3.16 Settlement and total vertical stress relationship for 100 mm dia. triaxial sample
1.8
1.6
Void ratio (e)
1.4
1.2
1.0
1 10 100
log σv (kPa)
Figure 3.17 Void ratio/ log σv relationship for 100 mm dia. triaxial sample
84
Samir Ashour Chapter 3: Methodology
A similar procedure was followed when preparing larger specimens, with the exception of
drainage which was allowed only from the top in order to measure pore water pressure from
the bottom of the sample at different positions (PPT1, PPT2, and PPT3) as shown in Figure
approximately 45 kg of slurry to be poured into each tank giving an initial slurry height of
448 mm. A saturated plastic porous filter was placed on top surface of the slurry. A perforated
25 mm thickness plate was placed at the top of the filter, and used as a consolidation loading
plate. A similar loading sequence was used then adopted to give pressures of 6, 12, 25, 50,
and 100 kPa. After the completion of the consolidation stage, the specimens surface was
flattened using metal ruler with 298 mm length, samples were taken for determining water
content. Stone column was installed in the clay bed (as explained in the following section
3.8.2), and the tank was moved to the loading frame for testing.
Figure 3.18 shows typical consolidation curves for each loading sequence. Accumulated
85
Samir Ashour Chapter 3: Methodology
Settlement, (mm)
12 kpa 25 kpa
5
5
10
10 15
Settlement, (mm)
Settlement, (mm)
50 kpa 0
100 kpa
5
5
10 10
15 15
20
Settlement (mm)
40
60
80
100
120
Figure 3.19 Settlement and total vertical stress relationship for 300 mm specimes
86
Samir Ashour Chapter 3: Methodology
In the case of 300 mm diameter specimens, settlement during the consolidation stage was
measured together with change in pore water pressure. Load was only increased when pore
pressure had reached an equilibrium value as there was a minimal increase in settlement
beyond this value as can be observed in Figures 3.18. Figure 3.20 shows typical response of
pore water pressure during loading stages during consolidation of slurry. It was difficult to
take measurements in the initial loading stage (7 kPa) where the data logger was showing un-
stabilised readings, this was probably as a result of trapped air in the system. During the
subsequent loading stages, pore water pressure increased to approximately the increment in
50
100 kPa
40
Pore water pressure (kPa)
PPT1
30
12.5 kPa 50 kPa PPT2
PPT3
20 25 kPa
10
0
0 5000 10000 15000 20000
Time (minute)
Figure 3.20 Variation of excess pore water pressure at PPT1, PPT2, and PPT 3 positions during
consolidation
87
Samir Ashour Chapter 3: Methodology
As stated in the literature review (section 2.4) stone columns can be constructed in the field
either by replacement or displacement method. Wood et al. (2000) and Sivakumar et al.
(2004b) investigated various methods for installing model stone columns in laboratory. They
used methods such as forming the column by freezing the stone material in a mould of the
same size of the required columns, another method was by displacing the soil using a rod with
the same size of the column and then filling the formed cavity by the column material, and
also the most common used method which was done by replacing a part of the soil and
compacting the stones in the cavity. A reduction in the density of the column upon thawing
was found in case of preformed frozen columns, while it was difficult to implement the
displacement technique in a small scale model due to the generated suction that may occur
In this study the replacement method was adopted because it proved to be repeatable. Thus a
28 mm diameter cavity at the centre of the sample was created by inserting a pipe with this
diameter to the clay bed and removing the soil included in the tube (Figures 3.21 and 3.22).
After creating the cavity, the granular aggregate was introduced in four stages and compacted
using a 1.5 Kg metal rod free falling through a fixed distance of 50 mm. 15 bowls were
applied to each layer. Aggregate between 1.18 mm and 2.1 mm in the size was used to form
The average bulk density of the column was calculated in each test. It was 1900 ± 45 kg/m3. It
is suggested that some of the variation in density may have been due to small increases in
88
Samir Ashour Chapter 3: Methodology
89
Samir Ashour Chapter 3: Methodology
The prepared sample was located in the triaxial cell, and saturated using increments of cell
pressure and back pressure. A Skempton 𝐵̅ parameter of over 0.95 was achieved. The sample
was then consolidated at an effective confining pressure of 100 kPa with a back pressure of
300 kPa.
i) Static loading: where the specimens were sheared at a constant axial strain rate of 0.03 mm/
minute.
ii) Cyclic loading: load was applied at different cyclic stress ratios (CSR) 0.6, 0.7, and 0.8,
each at three frequencies 0.5, 1 and 3 Hz. The cyclic stress was applied in sine wave form as
In both loading conditions, pore water pressure was measured at the bottom of the specimen.
qmax
q0=20 kPa
Time (s)
90
Samir Ashour Chapter 3: Methodology
Following installation of the column, the specimen was saturated slowly from bottom to top
using piezometer. This process took about 4 hours, in order to ensure that the stone column
was saturated and the air trapped between the particles was forced out, so as not to affect the
quality of pore water pressure measurement beneath the column. Then a sand blanket of 3 mm
thickness was placed covering the entire top surface of both soil and column to ensure that the
stress imposed from the surcharge load will be uniformly distributed to cover the entire area.
The footing was then located at the centre using a guide plate (see Figure 3.24). the sample
was then reconsolidated for 48 hours using a slightly higher pressure (120 kPa) than that used
to create samples (100 kPa), the extra pressure simulated the surcharge pressure on top of the
footing. Following this, in similar manner to the triaxial model, both static and cyclic loading
Loading ram
Plastic plate
Inflatable rubber ring
Loading foundation
3 mm Sand blanket
Clay bed
300 mm
Stone column
300 mm
400 mm
91
Samir Ashour Chapter 3: Methodology
Water content in the clay bed can provide a good indication of its strength and uniformity.
This was carried before testing (Figure 3.25) in order to ensure the quality and the general
A plastic pipe (similar to the one used to install the column) with 28 mm diameter and 45 mm
height was used to take the soil samples for water content test. Soil within the pipe was
loosened and a vacuum was used to lift out the soil from within. This technique can provide
samples along the depth of the clay bed and can be applied in different locations.
As observed in Figure 3.25 the water contents through the clay bed were consistent and
ranged between about 43 and 42.5 %. This range of water content was very close to that
estimated using the relationship of undrained strength and liquidity index (3.3.5). At about 43
Results showed that in most of the tests the variation in water content was comparatively
larger at the middle height of the specimen (± 0.5 %); this could be related to length of the
drainage path.
92
Samir Ashour Chapter 3: Methodology
50
Model I
100
Model II
150
200
Figure 3.25 Variation in water content before testing for both Models (I and II)
In order to investigate the quality of the prepared clay bed, a hand vane was used to determine
the variation in shear strength by testing clay at a range of positions near where samples were
taken for water content determinations. Hand vane tests were conducted in accordance with
BS 1377:1990. A 12.7 mm hand vane was used to limit clay disturbance to minimum. For the
100 mm diameter specimens, due to size of the specimen it was difficult to perform the hand
vane test at more than one location; however, three pilot tests in three specimens were carried
at different depth. In the large sample, the undrained shear strength was determined at four
(50, 100, 150, and 200 mm) depth and four horizontal positions (80 mm around the centre).
This procedure was carried out in three pilot tests only. After that the shear vane test was taken
in each sample before testing only in the centre and after testing in various locations to collect
as much data as possible to observe the changes that occurred before and after testing.
93
Samir Ashour Chapter 3: Methodology
Figure 3.26 shows the undrained shear strength values obtained from these tests including
estimated values using Equation 3.2 (3.3.5). The average undrained shear strength was 11.5
kPa with a variation of ± 0.5 kPa. This result was less than the estimated values with a
The results from both water content and undrained shear strength demonstrated excellent
50
Depth of clay bed, (mm)
100
150
200
Measured undrained shear strength (Hand vane)
Estimated undrained shear strength using Eq.3.1
250
Investigating the shape of the column after testing can provide a good indication of the failure
mechanisms and the interaction between the soil and the stone column. As it is not easy to
monitor changes during the test, observation at the end of the test is the only available guide
about the behaviour of the materials. Several techniques were investigated to assess condition
and shape of the column at the end of the test. These included grouting the column material,
94
Samir Ashour Chapter 3: Methodology
freezing, using both wax and resin to hold the aggregate material together and section slicing
method.
The simplest was slicing the specimen in half vertically through a selected section. This
method appeared to be a very effective, quick and showed the shape of the column. It was
95
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
4 CHAPTER 4
4.1 Introduction
This chapter presents the results of both 100 mm diameter triaxial specimens and foundation
Sections 3.6. The general behaviour of both soft clay (no column) and clay/stone column
composite system are discussed together with the influence of loading rate, and column stone
density on the bearing capacity. The discussion includes results from previous studies with
similar testing conditions in order to validate the current study work. The test programme
TRI-C-01 0.003
TRI-C-02 - 0 0.03 -
TRI-C-03 0.3
TRI-C.C-04 0.003
TRI-C.C-08 1805
TRI-C.C-09 1640
0.03
TRI-C.C-10 1500
LM-C-01 - 0
- 70 mm Footing
LM-C-02 - 0
Large Model 100 12 100 0.03
LM-C.C-03 16 1
1805 70 mm Footing
LM-C.C-04 16 1
96
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
representing soil/ stone column composite. These two types of specimens were tested under
the same testing conditions. After extruding the specimen from the consolidation chamber,
they were installed in a triaxial cell for loading. Each sample was first saturated and then
isotopically consolidated under effective confining pressure of 100 kPa, with back pressure of
300 kPa. At this stage drainage was allowed from the bottom of the specimen, as such a
volume change device was connected, allowing volume changes to be recorded during
consolidation. After consolidation, the entire top of the sample was loaded under undrained
condition.
4.2.1 Saturation
The specimen was saturated by increasing the pore water pressure so that any air in the void
space can be eliminated. This allows for reasonable, reliable and repeatable readings of pore
water pressure changes during loading stage (Head, 1998). As described on BS1377: Part
8(1990), there are several established methods for saturation of samples, the one applied in
this study is saturating samples by using the application of back pressures, so that air in the
system is forced into solution. To achieve over 95 % of saturation, the process were carried
out by raising the confining pressure and back pressure in alternate increments as shown in
Table 4.2.
As the tested material is characterized by its low permeability and the size of the sample is
considered large, therefore the time required to reach 95 % of saturation was quite long
between 3 and 5 hours. The stone column reinforced samples took a shorter time to saturate
97
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
0 - -1.2 - - 0 - -2.5 - -
50 40 38 - - 50 40 39.6 - -
4.2.2 Consolidation
After the saturation stage was complete the cell and back pressures were adjusted to allow the
sample to be isotopically consolidated under effective confining pressure of 100 kPa and back
pressure of 300 kPa. By allowing the water to drain from the sample to the back pressure
system, the pore water pressure was gradually decreased until it reached approximately the
value of the back pressure. 95% dissipation of pore water pressure was achieved after 24
hours.
The degree of consolidation was determined by expressing volume change at a given time as a
function of the final volume change. As expected consolidation took place quicker when stone
columns were installed compared to soil only (without stone column), where 90%
consolidation was achieved seven times faster, see results presented in terms of degree of
consolidation versus the square root of time in Figure 4.1. This result showed similar
98
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
behaviour to previous research by Sivakumar et al. (2004b) under similar test conditions.
consolidation between this current study and the one conducted by Sivakumar et al. (2004b),
this difference could be due to using the filter papers around specimens used in this current
study as there was no indication that they used such method in their study. In addition to this,
in case of reinforced specimens, the size of the column material particles used and the
100
90
80
Degree of consolidation, (%)
70
60
50
40
30
The unit cell approach was used to analyse the load – deformation data, where the column and
the surrounding soil were assumed to act as a single element with a homogeneous distribution
of stresses and strains as proposed by Balaam et al. (1977). Whilst this assumes that the
sample is homogeneous, this approach has proved satisfactory for many practical design
99
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
applications (Lee and Pande, 1998; Maakaroun et al., 2009; Shahu et al., 2000; Sivakumar et
The influence of the inclusion of stone column in soft clay bed was observed in all tests
curried out. A typical stress – strain relationship is presented in Figure 4.2. The results from
reinforced clay tests are compared with tests carried out on soft clay specimens only. It was
observed that the reinforced samples developed higher bearing capacity than the soil only
samples. The deviator stress at failure for soil samples was approximately 60 kPa, while with
the presence of full length stone column this increased to about 90 kPa. This is an increase of
especially; when many previous researches have reported that the area over 15 % is required
to show significant improvement in the bearing capacity. This significant improvement can be
attributed to higher internal angle of friction of the material used in this study for building the
Figure 4.3 shows results of tests conducted by Sivakumar et al. (2004b); Black et al. (2006);
Andreou et al. (2008); Najjar et al. (2010) at similar testing condition (undrained triaxial
condition at the same effective confining pressure) are shown together with results determined
in the current study. In terms of the degree of improvement of bearing capacity, results show
that a higher angle of internal friction would provide a greater improvement ratio at a lower
area replacement ratio. For example, using material with an angle of friction at 48o in this
study (comparing with 35o at the other studies) enhanced the degree of improvement by about
100
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Despite of the limited number of studies on stone columns using undrained triaxial conditions,
it is possible to obtain a relationship between the improvement ratio and the area replacement
100
80
Deviator stress, q, (kPa)
60
40 Reinforced soil
Soil/column (current
(current study)
study)
Reinforced
Soil/columnsoil (Sivakumar
(Sivakumar et al,
et al, 2004)
2004)
Soilonly
Soil only(Sivakumar,
(Sivakumar 2004)
et al, 2004)
20
Soil only (current study)
0
0 5 10 15
Vertical Strain, (%)
Figure 4.2 Stress strain behaviour of soft soil (no column) and soil/stone column composite.
0.4
0.35
Improvement ratio
0.3
Y = 4.0034(As) - 10.831(As)2
0.25
R² = 0.7338
101
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
In all tests, pore water pressures were measured, and expected application of stress on low
permeability soils induces excess pore water pressure. Through the process of consolidation,
pore water is reduced resulting in a gradual increase in effective stress. The speed of
consolidation amongst other factors depends on drainage path. Installation of stone column
Figure 4.4 shows that the pore water pressure in all cases increased tending towards their
highest value as the deviator stresses reaches their maximum value at axial strain of about 5
%. There was a reduction in pore water pressure at about 15 % in case of the reinforced soils.
This reduction during the undrained loading stage is possibly due to the dilation of the stone
column.
80
Increase in pore water pressure, (kPa)
60
40
Reinforced soil
Soil/column (current
(current study)
study)
20 Soil only (current study)
Reinforced soil
Soil/column (Sivakumar
(Sivakumar et al, 2004)
et al,2004)
Soil only (Sivakumar, 2004)
(Sivakumar et al, 2004)
0
0 5 10 15
Vertical Strain, (%)
Figure 4.4: Excess pore water pressure of soil (no column) and soil/ stone column composite.
The vertical deviator stresses at failure for both unreinforced and reinforced specimens under
undrained conditions are presented in Table 4.3, for tests under similar condition to this study.
The improvement occurred in terms of bearing capacity, which ranged between 27 and 56 %,
102
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
and was mainly dependent on the area ratio and the internal angle of friction of the column
materials.
The estimation of vertical stress at failure of a single column is included in Table 4.3. using
1+sin ∅′
σv = ( ) (4c + σr′0 ) 4.2
1−sin ∅′
In the above equation σr′0 was considered firstly equivalent to 2𝑐𝑢 as suggested by Hughes
and Withers (1974), and secondly was considered to be equivalent to the effective confining
stress in the test (100 kPa). These calculated values are presented in Table 4.3 as
σv(calculated) (1) and (2) respectively. Whilst both assumption provide an over estimation of
values for all tests, it shows that the first assumption (σr′0 =2𝑐𝑢 ) would provide a closer
prediction (within an average of 25 %) than considering the amount of the effective stress.
Table 4.3 Measured and calculated vertical stress obtained from different studies
After each test, the specimens were removed from the cell and split vertically along their
diameter to examine the mode of failure. It was expected that the columns would bulge at a
depth of 2 to 3 times the column diameter as suggested by Hughes and Withers, (1974). In all
cases the column bulging was found to be relatively uniform with the depth of the sample as
shown in Figure 4.5. This probably was due to the constant confining pressure along the depth
103
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
of the specimen (triaxial test condition), whereas in the field condition the confining stress
200
150
100
50
0
-50 0 50 100 150
Strain rate has been recognised as having a significant influence on the undrained shear
strength of saturated soils (Kimura and Saitoh, 1983; Nakase and Kamei, 1986). Skempton
and Bishop (1954) summarised that the strength of cohesive soil are generally sensitive to the
change in strain rate and any reduction in the strain rate would decrease the strength of the
soil, whereas the granular soils are generally independent of this factor.
The strain rate for shearing saturated soils under undrained condition can be determined using
𝜀𝑓 𝐿
𝑅𝑎𝑡𝑒 𝑜𝑓 𝑠𝑡𝑟𝑎𝑖𝑛 (𝑚𝑚⁄𝑚𝑖𝑛𝑢𝑡𝑒) = 100𝑡 4.3
𝑓
104
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Where: 𝜀𝑓 is the strain at failure and typically ranging from 15 to 20 % for normally
consolidated clay; and 𝐿 is the specimen length (mm); and finally 𝑡𝑓 (minute) is the time to
failure which is calculated from 𝑡100 which is the time required to ensure at least 95 %
The axial strain rate for consolidated undrained triaxial test with pore water pressure
measurement was carried out using the above expression and it was found to be 0.03 mm/
minute. This study included investigation of the influence of the strain rate on the undrained
shear strength of both soft soils and stone column reinforced soft soil.
Three strain rate values were applied to cover the range recommended for cohesive soils (0.3,
0.03, and 0.003 mm/minutes). Each test was repeated at least two times to ensure repeatability
in results.
Figure 4.6 (a) and (b) show the typical stress-strain curves of both treated and untreated soils.
Considering the load capacity of the both soils to be at 5 % of the specimen strain (equivalent
to 10 % of the foundation width) as suggested by Bowles (1996). Results show that there is
about increase of about 10 % change in deviator stress associated with increasing the
magnitude of rate of strain from 0.003 to 0.3 mm/minute in unreinforced specimens. In case
of reinforced specimens the trend was in reverse order where the deviator stress showed an
Conversely, if failure zone was considering to be at 8 % of the deformation (i.e. 58% of stone
column diameter (Hughes and Withers, 1974)), the deviator stress in the reinforced
specimens was almost fall toward same point, this could be due to the drainage path provided
by the presents of the stone column which results in reduction in the effect of pore water
105
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
pressures. Whereas in case of unreinforced soils the deviator stress was continued to increase
This suggests that with or without stone column results are not very sensitive to rate of strain.
This may be due to the fact that the soil examined is kaolin clay which has a high percentage
of silt content (about 50%). In addition to this it is worth noting that the soil with stone
80
60
Deviator stress, (kPa)
40
0.003 mm/min
20 0.03 mm/min
0.3 mm/min
0
0 2 4 6 8 10 12 14
Vertical strain, (%)
100
80
Deviator stress, (kPa)
60
40
0.003 mm/min
20 0.03 mm/min
0.3 mm/min
0
0 2 4 6 8 10 12 14
Vertical strain, (%)
Figure 4.6 Typical stress-strain curves: (a) soil (no column); (b) soil/stone column composite
106
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
The typical observed relationships of the change in pore water pressures and the axial strains
for both soil only and stone column/soil specimens for strain rates of 0.3 to 0.003 mm/minute
are given in Figure 4.7 (a) and (b). Results show that there is an approximate 20 % increase in
pore water pressure for reduction in speed of second order of magnitude for clay only
specimens. For the case of stone column reinforced soil, change in pore water pressure was
about 8 %. This was most likely due to shorter drainage path. Thus, the effect of rate of
testing stone column reinforced soil can be discounted when assessing their behaviour in case
Figure 4.8 shows the change in pore water pressure- strain rate at the peak deviator stress
level. The pore water pressure was normalised in terms of the effective confining pressure
(𝜎3′ ). Findings tend to confirm that reduction in pore water pressure at tested strain rate as
observed in this study. It also showed the difference in pore water pressure between treated
60
Change in pore water pressure, (kPa)
40
20 0.003 mm/min
0.03 mm/min
0.3 mm/min
0
0 2 4 6 8 10 12 14
Vertical strain, (%)
Figure 4.7(a) Typical change in pore water pressures-strain curves for soil only
107
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
60
0.003 mm/min
20
0.03 mm/min
0.3 mm/min
0
0 2 4 6 8 10 12 14
Vertical strain, (%)
Figure 4.7(b) Typical change in pore water presures-strain curves for soil/ stone column
composit
0.6
0.5
Δu/σ'₃
0.4
Soil only
soil/stone column
SivaKumar et al. (2004)
SivaKumar et al. (2004)
Andreou at al. (2008)
Andreou at al. (2008)
Maakaroun et al. (2009)
0.3
0.001 0.01 0.1 1
Strain rate, (mm/minutes)
Figure 4.8 Normalised excess pore water pressure – strain rate relasionship
108
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Young’s modulus of soil (E) is the proportion between the changes in deviator stress to
∆𝜎
𝐸= ∆𝜀
4.4
Secant modulus, which is defined as the slope of the secant drawn from the origin to a
particular point on the stress-strain curve (Schanz et al., 1999), can provide a convenient
measure of soil stiffness. This modulus is an important parameter for the design of projects
such as roads and railways, where it is used to predict settlement (Selig and Waters, 1994).
During this study the effect of rate of strain and inclusion of stone column on soil stiffness
were investigated by tracking the secant modulus behaviour of both reinforced and
unreinforced soils at different stress between 50 and 100% of the soils shear strength.
Figure 4.9 shows the relationship between the secant modulus of both soil only and soil
reinforced by stone column and the strain at different strain rates. Other data from different
studies are also included in this graph. Generally the deformation modulus of soil showed a
independency on the tested strain range for both treated and untreated soils specimens. In
addition to this, there is a trend of reduction in deformation modulus with increase in strain to
about 4 %. After that a very small reduction was noted for both soil and soil with column.
In case of soil/stone column composite, the secant modulus was increased by about 30%
compared to soil (no column). This improvement was higher between 1 and 3 % strain. For
strains higher than about 4 % the secant modulus decreases and the degree of improvement in
modulus decreased by about 17%. These trends showed a very close agreement with results
109
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
conducted by Sivakumar et al. (2004b) and Andreou et al. (2008) in terms of percentage of
7
Esoil/col _0.005 mm/min_SivaKumar et al. (2004)
6 Esoil _0.03 mm/min_Andreou at al. (2008)
Esoil/col _0.03 mm/min_Andreou at al. (2008)
5
4
3
2
1
0
0 1 2 3 4 5 6
Axial Strain, (%)
Figure 4.9 Soil modulus of both soils only and soil/ stone column composit
In order to relate soil stiffness to road and railway applications results were used for
evaluation of subgrade reaction, which is used mainly in the design of foundation under wheel
and concentrated loads. The modulus of subgrade reaction (𝑘𝑠 ) can be defined as the relation
𝑞
𝑘𝑠 = 𝛿
4.5
Bowles (1996) provided an empirical equation to estimate the modulus of subgrade reaction
Where: Factors (40, 50, 83, and 160) are dependent on the settlement(∆𝐻) at 25, 20, 12 and 6
mm respectively; and
110
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
𝑞𝑢𝑙𝑡
𝑞𝑎 = ⁄𝑆𝐹 4.7
Therefore, the subgrade modulus of reaction for both soil (with no column) and soil/ stone
Department of Transport (as for pavement foundations the modulus should be limited to small
strains in the elastic range) and at 6 mm and 10 mm as suggested by Bowles (1996) for
foundation design. Results illustrated in Figure 4.10 (a, b, and c) show that the strain rate has
no significant effect on modulus of subgrade reaction for both soils (with and with no
column) for each level of modulus deformation. However, the installation of the stone column
increased 𝑘𝑠 by about 30 % compared with the soil only specimens; this improvement was
observed in all deformation levels. It was also observed that the subgrade modulus of reaction
A comparison between test results and the calculated modulus of subgrade reaction using
equation (4.3) can be seen in Table 4.4 which shows good agreement as the empirical
111
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
35 Soil only
25
Ks, (MN/m2/m)
20
15
10
0
0.003 0.03 0.3
(a)
Strain rate, (mm/minutes)
35
at 6 mm deformation Thousands
30
25
Ks, (MN/m2/m)
20
15
10
0
0.003 0.03 0.3
(b) Strain rate, (mm/minutes)
35
at 10 mm deformationThousands
30
25
Ks, (MN/m2/m)
20
15
10
0
0.003 0.03 0.3
(c) Strain rate, (mm/minutes)
Figure 4.10 The subgrade modulus of reaction (𝑘𝑠 ) for both soils (with no column and soil/ stone
10 mm
112
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Most of the widely adopted design methods used for predicting bearing capacity of stone
columns are mainly dependent on two factors, shear strength of the surrounding soil and the
internal angle of friction of the column material (as discussed in Section 2.6). In the specific
case of constant surrounding soil properties, the friction angle of a fill material (i.e. granular
soil) plays a crucial role in estimating stresses that the soil/ stone column composite can carry.
This friction angle is dependent on many factors such as the magnitude of stresses; relative
density; and on the size and shape of soil particles (Bolton, 1991).
The impact of density on the shear strength of granular material has been widely investigated,
however, very limited information on the influence of this factor on the stone column
application is available. This might be due to that the final diameter of the column is only
Generally, for a given material, the internal angle of friction decreases with the reduction in
In this study, columns were constructed with three different densities 1805, 1642, and 1500
kg/m3 representing relative densities of 80, 42, and 2% respectively. All tests were conducted
at the same effective confining pressure of 100 kPa and same strain rate of 0.03 mm/minute.
The stress-strain results are shown in Figure 4.11. These results show that the degree of
Firstly, for columns made at the loose density (1500 kg/m3) (columns were installed by
pouring the aggregate from 500 mm height with no compaction was applied (simulating a
condition close to the minimum dry density of 1490 kg/m3, (maximum void ratio))) there was
about 25 % improvement.
113
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
At the second test where the density was increased by applying compaction to the column (the
material was filled in two layers and 10 tamps were applied each time by 1.5 kg rod) the
At the final test the column material was compacted to achieve the a density close to the
maximum density, the deviator stress reach its maximum at about 90 kPa which is about 31 %
The variation of bearing improvement ratio qt/qunt versus the relative density Dr is shown in
Figure 4.12. This indicates that the relative density may have a significant influence in
increasing the strength of soil/column system and accordingly the bearing capacity of the
ground.
100
80
Deviator stress, (kPa)
60
40
Soil only
1805kg/m3
1642kg/m3
20
1500kg/m3
0
0 2 4 6 8 10 12 14
Vertical strain, (%)
114
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
1.8
1.4
1.2
1
0 10 20 30 40 50 60 70 80 90
The pore pressure behaviour showed independency from the change in column density as can
be seen in Figure 4.13 which shows the measurements of pore water pressures against the
axial strain.
60
Change in pore water pressure, (kPa)
40
1805 kg/m3
20
1642 kg/m3
1500 kg/m3
Soil only
0
0 2 4 6 8 10 12 14
Vertical strain, (%)
115
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Results of this study together with these of other researchers are shown in Figure 4.14 in
terms of maximum deviator stress to effective confining pressure ratio (q/𝜎3′ ). Results show a
𝑞
𝜎3′
= 0.056 (𝛾𝑑 ) + 0.0972 4.8
y = 0.056x - 0.0972
R² = 0.9836
0.8
y = 0.0475x - 0.0124
R² = 0.5157
q/σ'3
0.6
Current study
Juran et al. (1988)
0.4 SivaKumar et al. (2004)
Andreou at al. (2008)
Najjar et al. (2009)
0.2
13 15 17 19
Column density, kN/m3
To conclude any reduction in relative density of the column leads to a reduction in the bearing
capacity of about 14 % (midpoint of upper and lower bounds). This could be due to a
reduction in the peak internal angle of friction of the column material and this might be
related to the increase in the dilation angle of the column particle. Therefore, it might be ideal
when designing stone columns to use the critical angle of friction which is function of both
116
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
The specimens were tested under a confined pressure of 120 kPa which was maintained
during the whole test procedure, aiming to simulate the surcharge on the ground in the field.
The model foundation was loaded at a constant strain rate of 0.03 mm/min, because It was
concluded from the previous tests (i.e. triaxial model) that this rate was sufficient to allow
excess pore water pressures to dissipate. The test was stopped when the footing displacement
Three pore water pressures transducers placed in different locations at the base of the
The main aim of these large tests was to investigate the condition within which the lateral
pressure distributed along the depth of the column (typically increase with the increase in
depth (Bowles, 1996)). Thus the mode of failure is investigated in more detail than that in the
triaxial condition.
Typical load-displacement relationship for both the reinforced and unreinforced soils are
presented in Figure 4.15. The maximum load carrying capacity of soil (no column) specimen
was found to be about 0.7 kN (equivalent to bearing pressure of 200 kPa). Penetration of 35
mm was required to achieve this. This bearing capacity was close to the predicted value of
117
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Where cu is the undrained shear strength of the soft soil (12 kPa), and Nc, Nq are the bearing
capacity factor (value of 5.14 and 1 respectively for zero internal angle of friction (φ = 0)) and
This value was approximately doubled when the 28 mm diameter column was installed
beneath the foundation of 70 mm diameter. Bearing capacity of the stone column reinforced
500
400
Foundation stress, (kPa)
300
200
100
Soil/column
Soil only
0
0 2 4 6 8 10 12 14
Axial strain, (%)
A reasonable agreement (within 13 %) between measured and predicted values was achieved,
which using the expression proposed by (Hughes and Withers, 1974) (Equation 4.10). This
was not surprising due to the fact that the empirical method considered only the influence of
strength of the surrounding soil and the degree of the internal friction, but not the effect of
118
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
area replacement ratio and the column density as discussed in chapter two (sections 2.7.1.2-
2.7.1.4).
1+sin ∅′
σv = (1−sin ∅′ ) (4𝑐𝑢 + σr′0 ) 4.10
Where ∅′ is the friction angle of the stone column material; σr′0 is the effective in situ lateral
stress and assumed to be 2𝑐𝑢 ; and 𝑐𝑢 is the undrained shear strength of the soil and equals to
12 kPa. The predicted value was 485 kPa assuming a factor of safety of 1 to allow the
comparison.
In addition to this, the degree of improvement was comparable with other researcher’s results
as can be noted in Table 4.5. The load bearing capacity ratios for before and after
Results included laboratory, field, and numerical studies on the improvement provided by
stone column are summarised in Table 4.5. Most of methods used to predict the bearing
capacity of a stone column considered the soil shear strength and the angle of shearing
resistance of the column material and ignoring the influence of column diameter (area
replacement ratio).
This database was used to developed a simple linear model to predict the ultimate bearing
capacity of reinforced soil as a function of cu, As, and ∅. (more information available in
Appendix D)
Figure 4.16 shows a scatter plot matrix between all variables. The bottom row of this scatter
plot matrix gives the scatter plots of qult against each of the other three input variables.
Additionally, the univariate relationship between our outcome variable qult and the input
119
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
120
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
φ qult
cu Cu
φ
φ
As AS
qult qult
Cu AS
Regression analysis is a powerful technique that can be used to address various research
questions. In this report, this was used to check how q levels are affected by As; cu and ∅. In
Multiple linear regressions were used to process the data presented database, in order to
determine the best-fitting line through the data points (this line is sometimes referred to as the
regression line). Multiple means we have more than one input variable (also known as
predictor); hence, we are trying to fit a plane or hyper-plane rather than a line (i.e. Combining
all the input variables by conducting a coefficient multiplied by each variable, and then
summing these variables). The idea was to use a linear combination of input variables to
121
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
model their relationship with an output variable (qult). After using R's lm() function, the model
looks as follows:
By examining equation 4.11, we observe that when fixing ∅ and AS, an increase, or decrease,
of cu by one unit, causes an increase, or decrease, in qult by 11.313 units. Similarly, when
fixing cu and AS, an increase, or decrease, of ∅ by one unit, causes an increase, or decrease, in
𝑞𝑢𝑙𝑡 by 14.629 units. Also, when fixing cu and ∅, an increase, or decrease, of AS by one unit,
Table 4.6 shows the values of four diagnostics that were used to examine the goodness of fit
of a model. Namely, these were: Deviance, BIC, R2, Adjusted R2 and the Residuals (these
factors were defined in Appendix D). As shown in this table the regression model has a
coefficient of determination R2 of 0.44, in other words 44% of the variation in the data can be
These data is also presented in Figure 4.17. The ultimate bearing capacity of the reinforced
soil measurements were normalised by dividing it by the soils bearing capacity before
reinforcement and plotted against the area replacement in terms of internal angle of friction of
the column material creating a linear relationships between the bearing capacity ratio and the
area replacement ratio at different internal angle of friction as presented in equations 4.12-
4.15:
122
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Φ=(42-43o), q ratio = 4* As + 2
As, (%)
Settlement estimation:
Various design approaches used for estimating settlement of a stone column – soil system, as
explained in the literature (section 2.6.2). They are based on the unit cell assumption. Priebe’s
method for settlement prediction is one of these approaches and considered to be the most
123
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
commonly used in practice. In order to compare test results, from this study firstly, one
dimensional settlement for the soil only was calculated using the following formula:
𝑆0 = 𝑚𝑣 ∆𝜎𝑧 𝐻 4.16
Where mv is the coefficient of volume compressibility of clay; Δσz is the vertical stress at
depth z due to the applied pressure at the surface and equals to (Ig P0);( Ig is the Influence
facto for the increase in vertical stress, P0 is the imposed load) and H is the thickness of the
soil layer.
Secondly, the settlement improvement factor was determined using the area replacement and
illustrated in chapter two (section 2.6.2.3). The internal angle of shearing resistance ranged
between 35o and 45o, thus a value of 45o was used to predict the improvement factor.
For example, based on the test conducted in this study, settlement estimation was based on
Equation 4.16, where Ig = 0.64644, and P0 = 100 kPa, leading to Δσz = 64.64 kPa, also using
mv as 9.803* 10-4 m2/kN (Section (3.3.6)). Settlement was estimated to be about 18 mm for
the clay only; and then the settlement of the reinforced soil can be estimated by dividing the
estimated settlement of the unreinforced soil by the improvement factor (n0) of 2.8 (obtained
from Priebe’s charts, Figure 2.9) resulting in a value of 6.4 mm, which was about six times
higher than the measured one. Although Priebe’s method demonstrated an overestimation to
the settlement of the reinforced soil, which confirmed the findings of Bouassida et al. (2003)
and Ellouze et al. (2010), the measured improvement factor of 2.1 showed reasonable
agreement with similar studies (Ambily and Gandhi, 2007; Black et al., 2011).
124
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Pore water pressure transducers PPT1,PPT2, and PPT3 were fitted at the centre of base, 50
mm from the centre and 100 mm from the centre respectively. Figures 4.18 and 4.19 show
results during the foundation loading on the soil (only) and soil/ stone column composite
specimens. When soils (no column) were tested, pore water pressure was seen to increase
only at PPT1 (the midpoint) when loading the foundation to about 3 kPa, while in the other
two (PPT2 and PPT3) transducers there were no observed change during the test.
On the other hand, when the column was installed, the PPT3 (close to the tank wall) did not
show any change in pore water pressure measurements during the loading. There was gradual
increase by about 5 kPa in pore water pressure at PPT2 (50 mm from the column); and
beneath the column at the centre of the specimen (PPT1), the pore water pressure build-up
started dissipating. This confirms that the column carried larger part of the load compared to
the surrounding soil. Also when the pore water pressure stopped building up under the
column and still increasing beneath the soil could indicate that the failure of the column had
occurred.
125
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
4
PPT1
Pore water pressure, (kPa) PPT2
0
0 5 10 15
Axial strain, (%)
Figure 4.18 Typical results for changes in pore water pressure for soil (no column) specimens
30
PPT1
PPT2
Pore water pressure, (kPa)
20
10
0
0 5 10 15
Axial strain, (%)
Figure 4.19 Typical results for changes in pore water pressure for soil/stone column specimens
126
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
In order to investigate the quality of the prepared clay bed, a series of shear vane tests were
carried out at different positions to determine the undrained shear strength of the clay
surrounding the column. Strength measurements are shown in Figure 4.20. The undrained
shear strength was determined at (50, 100, 150, and 200 mm) depth and at four horizontal
positions (around the midway between the column and the tank wall). In addition to this, water
content measurements were made at the same levels as shown in Figure 4.21. Comparing both
before and after testing for the shear vane results and water content, it can be observed that at
the first 50 mm thickness of the soil there was an increase in the soil strength of about 12%
Before testing
After testing
50
Depth of clay bed (kPa)
100
150
200
Figure 4.20 Hand vane shear test results before and after testing
127
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
Before testing
Depth of clay bed (mm)
50 After testing
100
150
200
After each test, the specimens were removed from the cell and split vertically along their
diameter to examine the mode of failure. A typical example is shown in Figure 4.22. In most
cases column bulge, took place at approximately at 120 mm below the surface, which is about
1.7 times the column diameter. This was very close to the expected value of 2 to 3 times the
128
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
300
250
200
150
100
50
0
0 50 100 150 200 250 300
order to investigate the influence of strain rate and column material density on the behaviour
Results showed that there was a considerable improvement in load carrying capacity and
stiffness of soft clay when stone column was installed. For example, replacing a small part of
soft soil (7 %) with granular materials that have a comparably a high angle of shear resistance
(48o) leads to an increase in bearing capacity and soil stiffness of approximately 30 % and
reduces the settlement by about 80 % compared with the settlement of unreinforced soils.
Although results from this study showed reasonable agreement with empirical estimations and
other published works in terms of load bearing capacity, observed settlement behaviour
129
Samir Ashour Chapter 4: Monotonic Loading Results and Discussion
showed a conservative figure not least when compared with the settlement estimated using
Priebe’s method.
Strain rate did not appear to have a significant influence on the level of improvements, where
results show that there is a decrease of less than 10 % in deviator stress for second order of
There was a strong indication suggests that the relative density of the column material may
have a significant influence on column bearing capacity. For example, increasing the column
density by about 25 % leads to an increase in soil/ stone column composite capacity of about
17 %.
130
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
5. CHAPTER 5
5.1 Introduction
This chapter presents results of both the undrained cyclic triaxial tests (Model I) and cyclic
foundation tests on large scale model (Model II) carried out on both soft clay and soft clay
/stone column composite. Dynamic stress levels ranging between 50 to 70 kPa (equivalent to
cyclic stress ratio (CSR) from 0.6 to 0.8) and frequencies ranging between 0.5 to 3Hz were
used, as described in Sections 3.7.1 and 3.7.2. All tests were conducted at an effective
confining pressure of 100 kPa. Comparison is made between the behaviour of the soft clay
Change in pore water pressure and deformation with number of cycles were examined. The
latter was used to identify threshold value of stress ratio. The variables examined during this
131
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
1
TRI-C-24 0.6 0.2 50 0.6
Typical results of cyclic undrained triaxial tests are shown in Figures 5.1-5.3, which shows
permanent axial strain against number of cycles. Results show the influence of cyclic stress
levels (50, 60 and 70 kPa) at different loading frequency (0.5, 1 and 3 Hz) on the deformation
pattern for soft soil specimens (no column) and soft soil/ stone column specimens. All
specimens were subject to a number of cycles of 10000 cycles, apart from those whom failed
132
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
16 CSR=0.6_Soil
qcyclic=50kPa_Soilonly
only
14 CSR=0.7_soil only
qcyclic=60 kPa_Soil only
CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
10 qcyclic=70 kPa_Soil/column
CSR=0.8_Soil/column
0
10 100 1000 10000
Number of cycles
Figure 5.1 Pemanent axial strain and number of cycles relationship: at 0.5Hz loading frequency
16 qcyclic=50 kPa_Soil
qcyclic=50kPa_Soil only only
14 qcyclic=60 kPa_soil
qcyclic=60 kPa_Soil onlyonly
Permanent axial strain, (%)
qcyclic=50 kPa_Soil/column
qcyclic=50 kPa_Soil/column
12
qcyclic=60 kPa_Soil/column
qcyclic=60 kPa_Soil/column
10 qcyclic=70 kPa_Soil/column
qcyclic=70 kPa_Soil/column
0
10 100 1000 10000
Number of cycles
Figure 5.2 Pemanent axial strain and number of cycles relationship: at 1 Hz loading frequency
133
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
16
CSR=0.6_Soil only
qcyclic=50kPa_Soil only
14
qcyclic=60 kPa_Soil
CSR=0.7_soil onlyonly
Permanent axial strain, (%) 12 qcyclic=50 kPa_Soil/column
CSR=0.6_Soil/column
qcyclic=60 kPa_Soil/column
CSR=0.7_Soil/column
10
qcyclic=70 kPa_Soil/column
CSR=0.8_Soil/column
8
0
10 100 1000 10000
Number of cycles
Figure 5.3 Pemanent axial strain and number of cycles relationship: at 3 Hz loading frequency
It can be noted that, for both soil (no column) and soil/ stone column composite, the progress
of cumulative strain was highly influenced by increases in cyclic deviator stress applied hence
CSR. However, in the range of the loading frequencies examined, soils (no column)
specimens were more sensitive to the change of frequency than the soil/ stone column
composite.
Results in Figures 5.1 to 5.3 showed that the permanent axial strains for soils specimens (no
column) were larger at lower frequency (0.5 Hz) and decreases with the increase in loading
frequency, for instance, when specimens were subjected to a cyclic stress of 50 kPa, the
permanent axial strains after 10000 cycle were 6.3, 4.4 and 6.3 % for frequencies 0.5, 1 and 3
permanent axial strain was very small (at 50 kPa cyclic stress they were 1.9, 1.8 and 1.6 % for
134
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
In addition, the effect of the cyclic stress on the development of permanent strain was
significant for both types of specimens. For example, when tests were performed at loading
frequency of 1 Hz, results (Figure 5.2) showed that soil specimens (no column) at an applied
cyclic stress of 50 kPa stabilised after 900 cycles at a strain level of 4.4 %, when the cyclic
stress was increased to 60 kPa a shear failure occurred at less than 100 cycles of load
applications. On the other hand, when the stone column was installed and specimens tested
under a 50 kPa cyclic stress application, the progress of the permanent strain stabilised after
1000 cycle at a value of 1.8 %. This value increases to reach 5.6 % by increasing the cyclic
stress to 60 kPa, and when the test was performed at a cyclic stress of 70 kPa the specimen
was failed in a shear failure after 30 cycles. Similar behaviour was shown when specimens
tested under the condition of other loading frequencies (0.5 and 3 Hz).
These results indicate that there can be a threshold stresses of 50 kPa and 60 kPa for soil (no
column) specimens and soil/ stone column composite specimens respectively above which the
total permanent strain increased rapidly with repeat load application; and failure can occur
during the first few cycles. Normalising these dynamic stresses in terms of static vertical
stress at failure for both types of specimens (with and without column) leads to have a CSR of
about 0.75. From literature this value varies from 0.4 to 0.8 and often is assumed to be 0.5.
(Seed et al., 1955; Heath et al., 1972; Brown et al., 1975; Brown, 1996; Andersen, 2004; Jiang
et al., 2010; Cai et al., 2013) This variation in the value of the threshold stress could be due
Figure 5.4 shows the influence of cyclic deviator stress on the permanent strain of both soil
only and soil/column composite specimens at load applications of 100 and 1000 cycle and
loading frequency of 1 Hz. Results showed that below the threshold stress level in soil (no
135
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
column) specimens were more sensitive to number of loading application than the reinforced
samples. For example, at 50 kPa cyclic stress, the permanent strain increased from 1.7 to 4.4
% by increasing number of cycles from 100 to 1000 cycle. On the other hand, the influence of
number of cycles were much smaller in case of soil/column specimens, where difference in
permanent strain between 100 and 1000 cycles at similar stress level was about 0.9 %. This is
probably due to the fact that the column materials were compacted due to load application,
increasing the density of the column leading to greater resistance with lower settlement.
Soil/column_100 cycle
Soil/column_1000 cycle
Permanent strain, (%)
10
15 Failure
20
25
Figure 5.4 Effect of loading application and cyclic stress leve at frequency of 1 Hz
Frequency of cyclic load is believed to have an important influence on the dynamic properties
of the saturated soft clay (Andersen et al., 1980), however, its effect is still attracting debate.
Jian-hua Wang et al. (1998), Zhou and Gong (2001), and Jiang et al. (2010) indicated that
136
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
both accumulative permanent strain and pore water pressure induced by cyclic loading will
increase with reduction in loading frequency. Whereas other published work indicated a
different observation, where they reported that the change in frequency has little or no
influence on the strength and deformation properties of soil (Ansal and Erken, 1989; Hyde et
al., 1993).
In this current study, an increase in permanent deformation for the soil only specimens was
observed with reduction in frequency of loading, as shown in Figure 5.5, for both soil
specimens (with and without column) at load applications of 100, 1000 and 10000 cycles. All
the test were undertaken at the same cyclic stress level (50 kPa) which was less than the
threshold stress. For stone column reinforced specimens, the influence of loading frequency
on the permanent stain was less pronounced compared with that for unreinforced specimens.
Also the effect of number of cycles was negligible after 1000 cycle.
Frequency, (Hz)
0 0.5 1 1.5 2 2.5 3 3.5
0
2
Permanent strain, (%)
Figure 5.5 Effect of loading frequency and number of cycles under cyclic stress of 50 kPa.
137
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
There are several models to predict permanent strain in soils subjected to cyclic loading
conditions. The most commonly used is the power model which takes into account a number
of factors such as number of loading cycles and both static and dynamic stresses (Li and
𝜀𝑝 = 𝐴𝑁 𝑏 5.1
Where 𝜀𝑝 is the permanent strain (cumulative strain) in %; N is the number of cycles; and A
and b are material parameters (dependent in the soil properties and stress level). Li and Selig
(1996), quantified the coefficient A by relating it to the physical state of the soil and the
𝜎
𝐴 = 𝑎( 𝑑⁄𝜎𝑠 )𝑚 5.2
Where a and m are the correlation parameters; 𝜎𝑑 𝑎𝑛𝑑 𝜎𝑠 are the deviator stress and static
strength respectively. Therefore, parameter (A) is highly dependent on the applied deviator
stress and the maximum deviator stress that the soil can take at static test condition.
Li and Selig (1996), recommend values for the above parameters for a number of soil types
based upon back calculation from a number of different studies (Table 5.2). Values for the
soft soils (no column) and the soil/ stone column composite parameters calculated in this
study through regression analysis of cyclic triaxial results are included (see Appendix F for
138
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Soil classification
Kaolin/stone
Model parameters Kaolin clay
ML MH CL CH column
(test condition)
(test condition)
As explained by Li and Selig (1996), that exponent b reflects the accumulative strain rate
under cyclic loading (i.e. that the higher the value of b the more the soil affected by number of
cycles). Whereas exponents a and m are related to the deformation occurred at the first cycle
and its relation to the cyclic stress ratio applied and also to the degree of softening of the
As can be seen in Table 5.2 the parameters value of the used soil (no column) is comparable
with the CL values, the reason for this is that the kaolin clay considered as a clay soil with
high plasticity. However, for soil/ stone column composite the exponent b showed a smaller
value than other type of soils; this might indicate that the stone column reinforced soils are
likely to be less sensitive to number of loading application. Whereas a value showed a higher
level (23.39) than the other soils, which indicates that, the soil/ column system could have a
139
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Figures 5.6 to 5.8 show the cyclic deformation behaviour for soil (no column) and soil/ stone
column composite. The influence of loading frequency; dynamic stress level and number of
0.9
qCSR=0.6_Soil
cyclic=50kPa_Soilonly
only
0.8 qcyclic=60 kPa_Soil only
Series4
Cyclic deformation, (mm)
Figure 5.6 Cyclic deformation during cyclic triaxial test: at 0.5 Hz loading frequency
qcyclic=50 kPa_Soil/column
CSR=0.6_Soil/column
0.5 qcyclic=60 kPa_Soil/column
CSR=0.7_Soil/column
q =70 kPa_Soil/column
cyclic
0.4 CSR=0.8_Soil/column
0.3
0.2
0.1
10 100 1000 10000
Number of cycles
Figure 5.7 Cyclic deformation during cyclic triaxial test: at 1 Hz loading frequency
140
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
0.7
CSR=0.6_Soil only
qcyclic=50kPa_Soil only
0.6
Elastic deformation, (mm)
CSR=0.7_soil onlyonly
qcyclic=60 kPa_Soil
CSR=0.6_Soil/column
qcyclic=50 kPa_Soil/column
0.5
CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
0.4 CSR=0.8_Soil/column
qcyclic=70 kPa_Soil/column
0.3
0.2
0.1
10 100 1000 10000
Number of cycles
Figure 5.8 Cyclic deformation during cyclic triaxial test: at 3 Hz loading frequency
In all specimens tested, the cyclic deformation decreases as the number of cycles increases.
This was due to hardening of soil specimens. This phenomenon occurs only at stresses below
the threshold level (i.e. 50 kPa for soil (no column) and 60 kPa for soil/ stone column
composite) above these stresses the cyclic strain sharply increases due to the shear failure
occurred.
Figures 5.6 to 5.8 have also indicated that the cyclic strain in both soils (with and without
stone column) are slightly influenced by the change in loading frequency as there was a small
For instance, when the specimens subjected to 50 kPa cyclic stress, the values of cyclic
deformation for soil (no column) specimens after 10000 cycle were 0.24, 0.22 and 0.18 mm
for frequencies 0.5, 1 and 3 Hz respectively. Whereas in the soil/ stone column specimens the
cyclic deformations were 0.23, 0.22, and 0.17 mm for frequencies 0.5, 1 and 3 Hz
respectively. This would indicate that there is a threshold elastic strain of 0.2 % for soil (no
141
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Soil resilient modulus (Mr) can be defined as the ratio of the dynamic deviator stress (𝑞𝑑 ) to
the elastic (recoverable) strain (𝜀𝑟 ) (Brown et al., 1975; Selig and Waters, 1994; Brown,
1996).
𝑞𝑑
𝑀𝑟 = 5.3
𝜀𝑟
This parameter is considered to be one of the fundamental factors of understanding the elastic
behaviour of soil foundations subjected to dynamic stresses. For instance, in railway track
analysis, resilient modulus is important for determining stresses and settlement in the track
Typical results of resilient modulus determination of both specimens type (soil (no column)
and soil/ stone column composite) under dynamic stresses of 50 and 60 kPa at different
loading frequencies (i.e. 0.5, 1 and 3.0 Hz) are shown in Figure 5.9 to 5.11 respectively.
45
40
Resilient modulus, (MPa)
35
30
CSR=0.6_Soil onlyonly
qcyclic=50kPa_Soil
25 CSR=0.6_Soil/column
qcyclic=50 kPa_Soil/column
CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
20
10 100 1000 10000
Number of cycles
Figure 5.9 Variation of resilient modulus with number of cycles at loading frequncy of 0.5 Hz
142
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
50
45
Resilient modulus, (MPa)
40
35
30
CSR=0.6_Soil only
qcyclic=50kPa_Soil only
25 CSR=0.6_Soil/column
qcyclic=50 kPa_Soil/column
CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
20
10 100 1000 10000
Number of cycles
Figure 5.10 Variation of resilient modulus with number of cycles at loading frequncy of 1 Hz
60
55
Resilient modulus, (MPa)
50
45
40
35
CSR=0.6_Soil only
qcyclic=50kPa_Soil only
30
CSR=0.6_Soil/column
qcyclic=50 kPa_Soil/column
25 CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
20
10 100 1000 10000
Number of cycles
Figure 5.11 Variation of resilient modulus with number of cycles at loading frequncy of 3 Hz
Results indicated that the resilient modulus increased with frequency in both specimen types
(20 % for soil/ stone column specimens and about 12 % for soil (no column) specimens). For
example, during the application of 50 kPa dynamic stress, the increase was from 40.6 MPa to
143
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
49.4 MPa in case of soil/ stone column composite specimens, while it was from 33.5 MPa to
In addition to this, Figure 5.12 shows that any increase in the dynamic deviator stress has a
negative impact on resilient modulus, which decreases as dynamic stress increases. For soil/
stone column specimens, when the dynamic deviator stress was increased from 50 to 60 kPa,
there was a general reduction in the amount of resilient modulus this reduction influenced by
the change of frequency as well (i.e. the higher the frequency the smaller the reduction). For
example, there was a reduction percentage in resilient modulus of 6.5, 7.6 and 13.3% when
Comparing the resilient modulus results with the secant modulus obtained from the static
triaxial test in the previous chapter (Section 4.2.4.3), which was about 9.5 and 8 MPa for
soil/column specimens and soil only specimens respectively; it was found that the resilient
modulus could reach up to four times the secant modulus value when the cyclic deviator stress
was 50 kPa.
55 Soil/column_3 Hz
Soil/column_1 Hz
50
Soil/column_0.5 Hz
Resilient modulus, (MPa)
45 Soil only_3Hz
Soil only_1 Hz
40 Soil only_0.5 Hz
35
30
25
20
15
40 50 60 70 80
Dynamic deviator stress, (kPa)
Figure 5.12 Variation of resilient modulus with respect to dynamic stress and loading frequency
144
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Under continued application of cyclic loading, excess pore water pressure gradually built up
over each cycle. Figures 5.13 to 5.15 show changes of pore water pressure corresponding to
the cyclic deviator stress at the three loading frequencies for both soil (no column) and soil/
90 CSR=0.6_Soil onlyonly
qcyclic=50kPa_Soil
CSR=0.7_Soil onlyonly
qcyclic=60 kPa_Soil
80
Increase in pore water pressure, kPa
CSR=0.6_Soil/column
qcyclic=50 kPa_Soil/column
70
CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
60 CSR=0.8_Soil/column
qcyclic=70 kPa_Soil/column
50
40
30
20
10
10 100 1000 10000
Number of cycles
Figure 5.13 Accumulation of excess pore water pressure at loading frequency of 0.5 Hz
90 CSR=0.6_Soil onlyonly
qcyclic=50kPa_Soil
Increase in pore water pressure, kPa
80 CSR=0.7_Soil onlyonly
qcyclic=60 kPa_Soil
CSR=0.6_Soil/column
q =50 kPa_Soil/column
70 cyclic
CSR=0.7_Soil/column
q =60 kPa_Soil/column
cyclic
60 CSR=0.8_Soil/column
q =70 kPa_Soil/column
cyclic
50
40
30
20
10
10 100 1000 10000
Number of cycles
Figure 5.14 Accumulation of excess pore water pressure at loading frequency of 1Hz
145
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
100 CSR=0.6_Soil
qcyclic=50kPa_Soilonly
only
80 CSR=0.6_Soil/column
qcyclic=50 kPa_Soil/column
CSR=0.7_Soil/column
qcyclic=60 kPa_Soil/column
70
CSR=0.8_Soil/column
qcyclic=70 kPa_Soil/column
60
50
40
30
20
10
10 100 1000 10000
Number of cycles
Figure 5.15 Accumulation of excess pore water pressure at loading frequency of 3Hz
As expected, the excess of pore water pressure increased with cyclic deviator stress at all
loading frequencies examined. Results show that in the case of most specimens, the pore
water pressure increased significantly at the initial stage of loading application (i.e. number of
cycle less than 150) then there was a gradual increase up to 1000 cycle before it started to
stabilise.
Figures 5.13 to 5.15 Also indicated that, the frequencies of loading did not show a noticeable
influence on the excess of pore water pressure in the soil /column specimens. Tests carried out
at 50 kPa cyclic stress the pore water pressures after 10000 cycles were 38, 37, and37 kPa for
frequencies 0.5, 1 and 3 Hz. At higher cyclic stress level (70 kPa) the specimen failed before
the pore water pressure could stabilise. Before failure it had increased rapidly to over than 75
kPa.
In contrast, for soil only samples results (Figures 5.13-5.15) showed that at a cyclic stress of
60 kPa, the excess pore water pressure increases rapidly to reach 70 kPa in few cycles leading
to a rapid failure in the sample, however, there was a small effect on the soil only specimens,
146
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
where the build-up of pore water pressure decreased with increase in the loading frequency
especially at the early stage of loading application (in the application of 50 kPa cyclic stress
and after 10000 loading cycle, pore water pressures were 51, 45.4 and 42.5 kPa for
Generally, in most tests conducted below the threshold stress levels pore water pressure
stabilises after 2000 cycles of loading applications. This behaviour is comparable with the
general behaviour of soft clay under cyclic loading condition (Miller et al., 2000; Indraratna et
al., 2009).
Figure 5.16, indicates that there was a critical excess of pore pressure and found at around 60
kPa for the soil specimens (no column) and about 70 kPa for the soil/column specimens ( pore
water pressure measurements at failure). These pore water pressure measures can be
normalised in terms of effective confining pressure (𝜎3′ = 100 𝑘𝑃𝑎) leading to ratios of 0.6
for soil (no column) and 0.7 for soil/ stone column composite. This also indicates that the
failure occurred when the pore water pressure approached the level of the dynamic stress.
80
70
Pore water pressure, (kPa)
60
50
Soil/column_0.5 Hz
Soil/column_1 Hz
40 Soil/column_3 Hz
Soil only_0.5 Hz
Soil only_1Hz
Soil only_3 Hz
30
40 50 60 70 80
Cyclic deviator stress, (kPa)
147
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Observation indicates that both soil specimens (no column) and soil/ stone column composite
specimens fail in shear failure mode when subjected to cyclic stress above the threshold cyclic
stress 50 kPa for soil and 60 kPa for soil/column specimens; this was irrespective of
frequency of loading. Also below these stresses both specimens did not fail at the end of the
Figures 5.17 and 5.18 show the mode of failure for both reinforced and unreinforced
specimens. Soil specimens (no column) specimens they failed within a range from 60 to 100
cycles of loading application (Figure 5.17) a slip surface was clear along weakest plan.
A combined mode of failure was observed in case of reinforced specimens (bulging and shear
failure (Figure 5.18)). Due to the large strain happened bulging failure took place first and
then by the increase in loading cycle shear failure accrued. The shear surface for the column
was at approximately the bottom 50 mm resulting in L/D ratio of 5.4 which is close to the
148
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
12
10
Sample height, (cm)
0
0 10 20
Horizontal measurment, (cm)
-5 5 15 25 35
20
15
Sample height, (cm)
10
149
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
The procedure for undertaking these tests is described in Sections 3.6.7 and 4.3. As
summarized in Table 5.3, the testing program included a series of cyclic zone load tests (Load
applied to column and portion of the surrounding soil in case of soil/column test). These tests
were conducted for both treated and untreated soils under two loading frequencies (1 and 3
Hz). Loading was applied in a sinusoidal wave form using two dynamic stress levels of 165
and 235 kPa, which were equivalent to cyclic stress ratios of 0.5 and 0.7 respectively.
The range of stress applied was based on the results of triaxial tests; particularly findings
relating to the threshold value, which was 70 % of the reinforced soil strength at failure. Thus
cyclic stress ratios of 0.5 and 0.7 were chosen for large scale test to evaluate this observation.
Stress at Loading
Pmax Pmin σmax σmin
Soil type Test failure CSR frequency
(kN) (kN) ( kPa) ( kPa)
(kPa) (f)
410
LM-C.C-09 1.05 0.2 280 45 0.7 1
Investigating the behaviour of permanent deformation of the reinforced soil on large scale
specimens indicated that the permanent stain was independent from the change in frequency
of loading; also the permanent deformation is highly dependent on the amount of the dynamic
150
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
stress applied and its relation to monotonic stress at failure. This confirms to the results
Figure 5.19 shows the effect of number of cycles; loading frequency, and cyclic stress ratio on
the measured permanent strain for both soil (no column) and soil/ column composite
specimens. It was observed when testing soil only specimens (no column) at a dynamic
loading stress of 165 kPa and loading frequency of 1 and 3 Hz, that the soil failed after 20
cycles of loading application; this was the cases at both loading frequencies. This was
expected as the loading stress was larger than the threshold ratio of the soil obtained in
Section 5.2.2. (As discussed in the methodology chapter this load was chosen due to the
limitation of the apparatus used in this study). On the other hand, soils/ stone columns
composite showed increased resistance to load, and the permanent deformation did not go
beyond 2 % at both loading frequencies. This increased to about 6 % when the dynamic stress
12
A Soil/column_CSR = 0.5; f=1 Hz
Soil/column _CSR = 0.5; f = 3 Hz
10 Soil/column_ CSR = 0.7; f = 1 Hz
Soil only_CSR = 0.5; f = 1 Hz
Soil only_CSR = 0.5; f = 3 Hz
Permanent strain, (%)
0
1 2001 4001 6001 8001
Number of cycles
151
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Predicted values of permanent strain obtained from Equation 5.2 are also shown in Figure
5.19. Material parameters a, b, and m, determined from the trixial test and explained in
During the application of dynamic stress at 165 kPa, permanent strain was predicted to be
increasing with the increase of number of applications reaching approximetly 1.6 % at 1500
cycles and then maintaining this level. Comparing this with the laboratory results, at the same
stress level the permanent strain increased to 1.8% at 1500 cycles and then continued at that
level to the end of the test (10000 cycles). A similar trend was observed when the stress level
was increased to 235 kPa. The estimated permanent strain was supposed to reach a maximum
of 4.6 % and stabilize after 1500 cycles, however, test results showed that the development of
permanent strain reached a maximum of 5.3 % and stabilised at this level after 4000 cycles.
Although the predicted permanent strain was underestimated compared with the measured
values, the difference seems to be considerable where it was within a range of 15 %. This
variation could be due to the difference in test conditions, as the parameters used in the
prediction were calculated from the triaxial test and for a limited CSR.
Figure 5.20 shows the deformation pattern after the cyclic loading application of 10000 cycles
at a cyclic stress of 165 kPa and loading frequency of 3 Hz. Although it was expected that the
columns will bulge at a depth of 2 to 3 times the column diameter as suggested by Hughes
and Withers, (1974), it appears that the column bulging under dynamic loading was different
from that of monotonic loading. In other words it was observed that the column bulge in two
sections (double bulging). These bulges were at 20 and 60 mm depth from the top, this can
152
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Double Bulging
300
280
260
240
Sample depth, (mm)
220
200
180
160
140
120
100
80
60
40
20
0
0 50 100 150 200 250 300
Sample diameter, (mm)
During cyclic loading pore water pressures were measured in two positions, at the centre
(PPT1) of the specimen base and at a distance of 50 mm from the centre (PPT2). The pore
water pressure measurements are illustrated in Figures 5.21 and 5.22; both show the pore
pressure changes that occurred under peak cyclic loading on the foundation plate.
12
Soil/column_CSR = 0.5; f=1 Hz
Soil/column _CSR = 0.5; f = 3 Hz C PPT1
10 Soil/column_ CSR = 0.7; f = 1 Hz
Change in pore water pressure, (kPa)
0
10 100 1000 10000
Number of cycles
Figure 5.21 Pore water pressure measurment at the centre of the specime
153
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
4
Soil/column_CSR = 0.5; f=1 Hz
Soil/column _CSR = 0.5; f = 3 Hz PPT2
Change in pore water pressure, (kPa)
Soil/column_ CSR = 0.7; f = 1 Hz
Soil only_CSR = 0.5; f = 1 Hz
3
Soil only_CSR = 0.5; f = 3 Hz
0
10 100 1000 10000
Number of cycles
Figure 5.22 Pore water pressure measurment at 50 mm from the centre of the specime
As was the case for triaxial tests, under continues application of cyclic loading, excess pore
water pressure was initially built up rapidly before it start to decrease (Figures 5.13 to 5.15).
Figure 5.21 shows the change in pore water pressure under the stone column at CSR of 0.5
was increased from about 1.5 kPa at the initial stage to about 8.7 kPa during the loading
application of 1000 cycle, after this the pore water pressure start to decrease gradually.
Similar behaviour was observed using CSR of 0.7 with larger magnitude of pore water
pressure (about 10 kPa). However, this behaviour was different at the other position (PPT2);
the change in pore water pressure was almost steady during the first 2000 cycle then started to
rapidly increase.
Soil stiffness (K) is defined as the ratio of the dynamic applied load (𝑃𝑑 ) on the foundation
154
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Stiffness of the treated and untreated soil tested at CSR of 0.5 and 0.7 at loading frequencies
of 1 and 3.0 Hz are shown in Figure 5.23 . Results show that after applying 10000 cycles the
stiffness for soil/ stone column composite increased from 1770 to 2050 kN/m as frequency
Soil specimens (no column) at both frequencies 1 and 3 Hz showed low stiffness (745 and
960 kN/m respectively) before failure occurred after 100 and 330 cycles respectively. This
was as a result of the high deformation occurred as shown before in Figure 5.19 (Points A and
B).
In addition, it was noticeable that the stiffness at soil/ stone column composite, at CSR of 0.5,
was decreased gradually (from1700 to 1630 kN/m at frequency of 3 Hz and from 1540 to
1450 kN/m at frequency of 1 Hz) during the first 100 cycles then increased by 20 % before it
stabilised after 1500 cycles. This is likely due to the stiffening effect of stone column and its
provision as a drainage path leading to reduction in pore water pressure (Figure 5.21 points C
and D). As pore water pressures start decreasing from 8.7 to around 6 kPa after 1000 cycles
This behaviour was different when the dynamic stresses was increased (CSR = 0.7), where
the stiffness of the reinforced soil was almost constant during the first 100 cycles at 1320
kN/m then it started to decrease slowly to reach 1100 kN/m at the end of the test (10000
155
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
2500
1500
1000
Soil/column_CSR = 0.5; f=1 Hz
Soil/column _CSR = 0.5; f = 3 Hz
Soil/column_ CSR = 0.7; f = 1 Hz
500
Soil only_CSR = 0.5; f = 1 Hz
Soil only_CSR = 0.5; f = 3 Hz
Failure
0
10 100 1000 10000
Number of cycles
The applied dynamic stress through the rail is transmitted through the ballast layer to the
subgrade level, played the main role in controlling the amount of permanent deformation of
the subgrade. Therefore, as suggested by Heath et al. (1972); Brown (1996) and Frost et al.
(2004) identifying the threshold stress would help for determining the thickness of unbound
layer.
From looking at literature and monotonic triaxal test results presented here it can be seen that
the bearing capacity of soft clay subgrade soils can be improved by (15 to 40 %) if the soil
was replaced by granular material having an angle of internal friction between 35 and 48 o
formed in columns with an area replacement ratio between 5 and 15 %. However, in this
current study using an angular crushed aggregate with an internal angle of friction and area
replacement ratio of 7 % have led to improved the bearing capacity of the soil by 30 % under
156
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
This improvement in the bearing capacity under the monotonic loading showed an increase in
the threshold cyclic stress ( under the cyclic loading condition) of about 15 % and 70 %
reduction in settlement; also showed an increase by 25 % in the resilient modulus of the soil.
The threshold relation between reinforced and reinforced soils can be related to limiting the
allowable stress applied on the soft subgrade via the relationship between the internal angle of
friction, area replacement ration and the undrained shear strength of the soil (Section 4.3.1).
This limit then can be applied in the design procedure of, for instance, railway tracks.
For example, this improvement in the soil properties could help to reduce the thickness of the
ballast layer. As one of the important performing functions of this layer is to reduce stresses
transferred through the sleepers from the train to the acceptable stress levels that the subgrade
soils can carry. Thus for a 100 kN axel wheel load using Talbot’s equation (Section 3.7.2), the
ballast depth (h) can be reduced by 15 % from 500 mm to 425mm. This can save about 225m3
of ballast per one kilometre of railway track which can cost around 2000 Pounds sterling per
Another impact is that, by simulating train movement with loading frequencies of 0.5, 1.0 and
3.0 Hz representing train speed of 35, 70 and 225 km/hr, stone columns reinforced subgrades
showed stability and less pronounced influence toward both the change in loading frequency
and long term loading application. This can increase the period of maintenance and provide
A comparison between the behaviour of both reinforced and unreinforced soils under the
157
Samir Ashour Chapter 5: Cyclic Loading Results and Discussion
Results indicated that the presence of stone column enhanced the threshold cyclic stress of the
soft soil by about 15 % and reduced the amount of permanent strain by more than double
when testing the soils below their threshold cyclic stress. Stone columns also help reduce pore
water pressure by providing a drainage path. This might be one of the reasons that allow the
Changes in frequencies had no significant influence on the total strain of reinforced soil, but
these changes do affect the stiffness. Stiffness of the soil with the stone column was about
25% higher at 3.0 Hz compared to that at 0.5 Hz. Resilient modulus of the reinforced soil
158
Samir Ashour Chapter 6: Conclusion
5. CHAPTER 6
6 CONCLUSION
6.1 General
The vibro stone column technique is commonly used to improve bearing capacity and reduce
settlement of soft soils. They have been used worldwide to provide economic, flexible
(applicable to various types of week soils) and sustainable solution. More recently, stone
columns have been used in the UK and Europe to improve soft subgrade soils subjected to
Performance and failure mechanisms of both isolated and grouped stone columns have been
investigated under the application of monotonic loading but not under cyclic loading. This
This research focused on the load/deformation behaviour of single stone column subjected to
monotonic and cyclic loading in both small scale models (triaxial specimen 100 mm diameter)
and large scale models (300 mm diameter * 300 mm height). In both cases tests were
kPa) reinforced with 28 mm diameter stone columns. Effects of loading frequencies (0.5-3Hz)
at a range of cyclic stress ratios (0.5-0.8) were examined. In addition to the influence of rate
of strain on shear strength on both soils (no column) and soil/ stone column composite in the
monotonic loading conditions (section 6.2) and cyclic loading conditions (Section 6.3). The
159
Samir Ashour Chapter 6: Conclusion
6.2.1 General
Replacing 7 % of soft clay with granular materials, which had an internal of friction angle of
48o, led to increase the bearing capacity and stiffness of the soil by a factor of 1.3. Compared
to predicted values and previously reported results reasonable agreement was found in terms
of bearing capacity, However, measured value for settlement were about 83 % lower than
Both soils (no column) and soil/ stone column composite specimens showed a small
influence from rate of strain; for soils (no column) there was an increase of about 10
% deviator stress associated with increasing the magnitude of rate of strain from 0.003
to 0.3 mm/minute. However, in the case of reinforced specimens the trend was
reversed order where the deviator stress showed an increase of 8 % with the reduction
The presence of stone column helped reduce build-up of pore water pressure and there
was limited influence from rate of strain on pore pressure as compared to soils (no
column). The latter showed pore water pressure reduction with increase in strain rate.
In case of soil/stone column composite, the secant modulus increased by about 30%
compared to soil (no column). This improvement was higher between 1 and 3 % strain
and for strains higher than about 4 % the secant modulus decreases and the degree of
160
Samir Ashour Chapter 6: Conclusion
The secant modulus of soil was seen to be independent of strain across the range
examined (0.3 to 0.003) for both soils (no column) and soil/ stone column composite
specimens.
In model I (Triaxial tests), the column bulging was found to be relatively uniform with
the depth of the sample, this was probably due to the constant confining pressure
along the depth. However, in model II (Large scale tests), bulging took place at
approximately 120 mm below the surface, which was about 1.7 times the column
diameter. This was very close to the expected value of 2 to 3 times the column
diameter
capacity of the soil/ stone column composite by about 14 %, this could be due to a
reduction in the peak internal angle of friction of the column material and this might
A change in aggregate density of the column has showed negligible effect on pore
pressure measurement.
Below the threshold stress level, reinforced soils were less sensitive to the number of
loading application than soils (no column). This is likely due to compaction of the
column materials that occurs during load application, leading to increase in density
161
Samir Ashour Chapter 6: Conclusion
Presence of a stone column enhanced the threshold cyclic stress of the soil by about 15
% and reduced the amount of permanent strain by half relative to soil without column.
Stone columns helped reduce pore water pressure by providing a drainage path. This
may be one of the reasons that allows the reinforced soil to support higher cyclic stress
Both reinforced and unreinforced soils failed in shear mode when subjected cyclic
stress above the threshold cyclic stress. This failure mode was not affected by
frequency of loading and below threshold stress both specimens did not fail
Large scale tests showed that the column bulging that forms under dynamic loading
was different from that of monotonic loading. It was observed that the column bulge
in two sections (double bulging) under cyclic loading. These bulges were at a depth of
0.7 and 2 times the column diameter below the surface of the samples.
influence by the change in frequency of loading. However, in the case of soil (no
those at higher frequency, about 55 % increase when frequency was changed from 3 to
0.5 Hz.
Below the threshold stress, the resilient modulus of both reinforced and unreinforced
soils increase with the increase of loading frequency. The value for soils/ stone
162
Samir Ashour Chapter 6: Conclusion
reinforced soil subjected to cyclic loading there are still important areas requiring further
Future work also needed to be carried out to investigate the effectiveness of stone
columns under the condition of fluctuating ground water table associated with the
It is anticipated that fines from surrounding soil may ingress into the stone column.
Investigating the rate of movement of fine and their effect on durability of the stone
column.
Effect of embankment height on the general behaviour of the soil/ stone column
composite soils.
Study the influence of stone columns grid pattern and spacing on deformation of track
Relative density of the column has an influence on the bearing capacity of the soil/
stone composite soils. Therefore, it might be ideal to investigate the effect dilatancy of
163
Samir Ashour Reference
REFERENCE
Aboshi, H., Ichimoto, E., Harada, K., et al. (1979) "The composer: A method to improve the
characteristics of soft clays by inclusion of large diameter sand columns". Proc., Int. Conf.
on Soil Reinforcement., E.N.P.C.,. Paris.
Adalier, K. and Elgamal, A. (2004) Mitigation of liquefaction and associated ground
deformations by stone columns. Engineering Geology, 72: (3): 275-291.
Ahmadi, M. and Robertson, P. (2004) Calibration chamber size and boundary effect for CPT
q (c) measurements. Geotechnical and Geophysical Site Characterization. vols, 1: 829-
833.
Ali, K., Shahu, J.T. and Sharma, K.G. (2014) Model tests on single and groups of stone
columns with different geosynthetic reinforcement arrangement. Geosynthetics
International, 21: 103-118(115).
Ambily, A. and Gandhi, S. (2004) Experimental and theoretical evaluation of stone column in
soft clay. ICGGE-2004, 201-206.
Ambily, A. and Gandhi, S. (2007) Behavior of Stone Columns Based on Experimental and
FEM Analysis. Journal of Geotechnical and Geoenvironmental Engineering, 133: 405-
415.
Andersen, K.H. (2004) "Cyclic clay data for foundation design of structures subjected to wave
loading". Proceedings of the International Conference on Cyclic Behaviour of Soils and
Liquefaction Phenomena, CBS04, Bochum, Germany.
Andersen, K.H., Rosenbrand, W.F., Brown, S.F., et al. (1980) Cyclic and static laboratory
tests on Drammen clay. Journal of the Geotechnical Engineering Division, 106: (5): 499-
529.
Anderson, W., Pyrah, I. and Fryer, S. (1991) A clay calibration chamber for testing field
devices. Geotechnical Testing Journal, 14: (4).
Andreou, P., Frank, R., Frikha, W., et al. (2008) Experimental study on sand and gravel
columns in clay. Proceedings of the ICE - Ground Improvement, 161: 189-198.
Ansal, A. and Erken, A. (1989) Undrained Behavior of Clay Under Cyclic Shear Stresses.
Journal of Geotechnical Engineering, 115: 968-983.
Babu, M.D., Nayak, S. and Shivashankar, R. (2013) A critical review of construction, analysis
and behaviour of stone columns. Geotechnical and Geological Engineering, 31: (1): 1-22.
Bae, W.-S., Shin, B.-W., An, B.-C., et al. (2002) "Behaviors of Foundation System Improved
with Stone Columns". The Twelfth International Offshore and Polar Engineering
Conference. Kitakyushu, Japan.
Balaam, N.P. and Booker, J.R. (1981) Analysis of rigid rafts supported by granular piles.
International Journal for Numerical and Analytical Methods in Geomechanics, 5: (4):
379-403.
Balaam, N.P., Brown, P.T., Brown, H.G., et al. (1977) Settlement Analysis of Soft Clays
Reinforced with Granular Piles. University of Sydney, School of Civil Engineering, Civil
Engineering Laboratories.
164
Samir Ashour Reference
Barksdale, R.D. and Bachus, R.C. (1983) "Design and construction of stone column ". Federal
Highway Administration, RD-83/026.
Bell, A. (2004) The development and importance of construction technique in deep vibratory
ground improvement. Ground and Soil Improvement, 103-111.
Bergado, D., Alfaro, M. and Chai, J. (1991) The granular pile: its present state and future
prospects for improvement of soft Bangkok clay. Geotechnical Engineering, 22: 143-175.
Bergado, D., Huat, S.H. and Kalvade, S. (1987) "Improvement of soft Bangkok clay using
granular piles in subsiding environment". In: Proceedings, 5th international geotechnical
seminar on case histories in soft clay. Singapore.
Bergado, D. and Lam, F.L. (1987) Full scale load test of granular piles with different densities
and different proportions of gravel and sand on soft Bangkok clay. Soil Found, 27: (1): 86-
93.
Black, J., Sivakumar, V. and Bell, A. (2011) The settlement performance of stone column
foundations. Geotechnique, 61: 909-922.
Black, J., Sivakumar, V., Madhav, M., et al. (2006) An Improved Experimental Test Set-up to
Study the Performance of Granular Columns. Geotechnical Testing Journal, 29: 14195.
Bolton, M.D. (1991) A guide to soil mechanics. Universities Press.
Bouassida, M., Guetif, Z., De Buhan, P., et al. (2003) Variational approach for settlement
estimation of a foundation on soil reinforced by columns. Revue Française de
Géotechnique, 102: (1): 21-29.
Bowles, J.E. (1996) "Foundation analysis and design". Fifth ed., McGraw-Hill, New York.
Brough, M., Ghataora, G., Stirling, A., et al. (2003) Investigation of railway track subgrade. I:
In-situ assessment. Proceedings of the ICE-Transport, 156: (3): 145-154.
Brown, S.F. (1996) Soil mechanics in pavement engineering. Geotechnique, 46: 383-426.
Brown, S.F. and Hyde, A. (1975) Significance of cyclic confining stress in repeated-load
triaxial testing of granular material. Transportation Research Record, (537).
Brown, S.F., Lashine, A.K.F. and Hyde, A.F.L. (1975) Repeated load triaxial testing of a silty
clay. 25: 95 - 114.
BSI, B.S.I. (2005) "Execution of special geotechnical works - Ground treatment by deep
vibration". London, BS EN 14731:2005, BSI.
Cai, Y., Gu, C., Wang, J., et al. (2013) One-Way Cyclic Triaxial Behavior of Saturated Clay:
Comparison between Constant and Variable Confining Pressure. Journal of Geotechnical
and Geoenvironmental Engineering, 139: 797-809.
Cai, Y., Sun, Q., Guo, L., et al. (2015) Permanent deformation characteristics of saturated
sand under cyclic loading. Canadian Geotechnical Journal, 52: (6): 1-13.
Charles, J. and Watts, K. (1983) "Compressibility of soft clay reinforced with granular
columns". Proceedings of the 8th European conference on soil mechanics and foundation
Engineering, Helsinki.
Chow, Y. (1996) Settlement Analysis of Sand Compaction Pile. Soils and Foundations, 36:
(1): 111-113.
165
Samir Ashour Reference
Ellouze, S., Bouassida, M., Hazzar, L., et al. (2010) On settlement of stone column foundation
by Priebe's method. Proceedings of the ICE-Ground Improvement, 163: (2): 101-107.
Fatahi, B., Khabbaz, H. and Le, T.M. (2012) Improvement of rail track subgrade using stone
columns combined with geosynthetics. Advances in Transportation Geotechnics Ii, 202-
206.
Françaises, N. (2005) Recommendations for the design, calculation, installation and control of
stone column foundations under buildings and structures sensitive to settlement. Revue
Française de Géotechnique, 111: (2): 3-16.
Frikha, W., Bouassida, M. and Canou, J. (2014) Parametric Study of a Clayey Specimen
Reinforced by a Granular Column. International Journal of Geomechanics.
Frost, M.W., Fleming, P.R. and Rogers, C.D. (2004) Cyclic triaxial tests on clay subgrades
for analytical pavement design. Journal of transportation Engineering, 130: (3): 378-386.
Gniel, J. and Bouazza, A. (2007) "Methods used for the design of conventional and
geosynthetic reinforced stone columns". 10th Australia New Zealand Conf. on
Geomechanics (Australian Geomechanics Society and New Zealand Geomechanics
Society 21/10/2007 to 24/10/2008). Geomechanics Society.
Gniel, J. and Bouazza, A. (2009) Improvement of soft soils using geogrid encased stone
columns. Geotextiles and Geomembranes, 27: (3): 167-175.
Greenwood, D. (1970) "Mechanical improvement of soils below ground surface".
Proceedings, Ground Engineering Conference, Institution of Civil Engineers London,
UK.
Greenwood, D. (1991) Load tests on stone columns. Deep Foundation Improvements:
Design, Construction, and Testing, 1089: 148.
Greenwood, D. and Kirsch, K. (1984) "Specialist ground treatment by vibratory and dynamic
methods".
Gu, C., Wang, J., Cai, Y., et al. (2012) Undrained cyclic triaxial behavior of saturated clays
under variable confining pressure. Soil Dynamics and Earthquake Engineering, 40: 118-
128.
Han, J. (2015) Recent research and development of ground column technologies.
Proceedings of the Institution of Civil Engineers - Ground Improvement, 168: (4): 246-
264.
Hanna, A., Etezad, M. and Ayadat, T. (2013) Mode of failure of a group of stone columns in
soft soil. International Journal of Geomechanics.
Head, K.H. (1996) Manual of Soil Laboratory Testing, Permeability, Shear Strength and
Compressibility Tests. Wiley.
Head, K.H. (1998) Manual of Soil Laboratory Testing, Effective Stress Tests. Wiley.
Heath, D., Waters, J., Shenton, M., et al. (1972) "Design of conventional rail track
foundations". ICE Proceedings. Thomas Telford.
Hu, W. (1995) Physical modelling of group behaviour of stone column foundations.
University of Glasgow.
166
Samir Ashour Reference
Hughes, J.M.O. and Withers, N.J. (1974) Reinforcing of soft cohesive soils with stone
columns. GROUND ENGINEERING, 7.
Hughes, J.M.O., Withers, N.J. and Greenwood, D.A. (1976) A field trial of the reinforcing
effect of a stone column in soil. Ground treatment by deep compaction.
Hyde, A., Yasuhara, K. and Hirao, K. (1993) Stability Criteria for Marine Clay under One‐
Way Cyclic Loading. Journal of Geotechnical Engineering, 119: 1771-1789.
ICE (1987) "Specification for Ground Treatment". Thomas Telford.
Indraratna, B., Attya, A. and Rujikiatkamjorn, C. (2009) Experimental Investigation on
Effectiveness of a Vertical Drain under Cyclic Loads. Journal of Geotechnical and
Geoenvironmental Engineering, 135: 835-839.
Jefferson, I., Gaterell, M., Thomas, A.M., et al. (2010) Emissions assessment related to vibro
stone columns. Proceedings of the ICE-Ground Improvement, 163: (1): 71-77.
Jiang, M., Cai, Z., Cao, P., et al. (2010) "Effect of Cyclic Loading Frequency on Dynamic
Properties of Marine Clay". Soil Dynamics and Earthquake Engineering. American Society
of Civil Engineers 240-245.
Juran, I. and Guermazi, A. (1988) Settlement Response of Soft Soils Reinforced by
Compacted Sand Columns. Journal of Geotechnical Engineering, 114: 930-943.
Killeen, M. (2012) Numerical modelling of small groups of stone columns. PhD, National
University of Ireland, Galway.
Killeen, M. and McCabe, B.A. (2014) Settlement performance of pad footings on soft clay
supported by stone columns: A numerical study. Soils and Foundations, 54: (4): 760-776.
Kim, B.-I. and Lee, S.-H. (2005) Comparison of bearing capacity characteristics of sand and
gravel compaction pile treated ground. KSCE Journal of Civil Engineering, 9: (3): 197-203.
Kimura, T. and Saitoh, K. (1983) The influence of strain rate on pore pressures in
consolidated undrained triaxial tests on cohesive soils. Soils and Foundations, 23: (1): 80-90.
Kolekar, Y.A., Mir, O.S. and Dasaka, S.M. (2011) " Behaviour of stone column reinforced
marine clay under static and cyclic loading". Proceedings of Indian Geotechnical
Conference. December 15-17, Kochi.
Lackenby, J., Indraratna, B., McDowell, G., et al. (2007) Effect of confining pressure on
ballast degradation and deformation under cyclic triaxial loading. Geotechnique, 57: (6):
527-536.
Li, D. and Selig, E. (1996) Cumulative Plastic Deformation for Fine-Grained Subgrade Soils.
Journal of Geotechnical Engineering, 122: 1006-1013.
Loh, R.B. and Nikraz, H.R. (2012) Compressibility and Strength of Clay Subgrade of
Railroad Foundation in Highly Saturated Condition. Journal of GeoEngineering, 7: (2): 043-
051.
Maakaroun, T., Najjar, S. and Sadek, S. (2009) "Effect of Sand Columns on the Load
Response of Soft Clays". Contemporary Topics in Ground Modification, Problem Soils,
and Geo-Support. American Society of Civil Engineers 217-224.
Madun, A. (2012) Seismic evaluation of vibro-stone column. PhD, University of
Birmingham.
167
Samir Ashour Reference
Madun, A., Jefferson, I., Foo, K., et al. (2012) Characterization and quality control of stone
columns using surface wave testing. Canadian Geotechnical Journal, 49: (12): 1357-1368.
McCabe, B.A., Nimmons, G.J. and Egan, D. (2009) A review of field performance of stone
columns in soft soils. Proceedings of the Institution of Civil Engineers-Geotechnical
engineering, 162: 323–334.
McKelvey, D. (2002) The performance of vibro stone column reinforced foundations in
deep soft ground. PhD, Queen's University of Belfast.
McKelvey, D. and Sivakumar, V. (2000) "A review of the performance of vibro stone column
foundations". Proceedings of the 3rd International Conference on Ground Improvement
Techniques, Singapore.
McKelvey, D., Sivakumar, V., Bell, A., et al. (2002) Shear strength of recycled construction
materials intended for use in vibro ground improvement. Proceedings of the ICE-Ground
Improvement, 6: (2): 59-68.
Miller, G.A., Teh, S.Y., Li, D., et al. (2000) Cyclic Shear Strength of Soft Railroad Subgrade.
Journal of Geotechnical and Geoenvironmental Engineering, 126: 139-147.
Mitchell, J.M. and Jardine, F.M. (2002) A Guide to Ground Treatment. Construction
Industry Research & Information Association.
Mohanty, P. and Samanta, M. (2015) Experimental and Numerical Studies on Response of the
Stone Column in Layered Soil. International Journal of Geosynthetics and Ground
Engineering, 1: (3): 1-14.
Moseley, M.P. and Kirsch, K. (2004) Ground improvement. London; New York: Spon
Press.
Munfakh, G.A. (1984) "Soil Reinforcement by stone columns–varied case applications".
Munfakh, G.A. (1997) Ground improvement engineering–the state of the US practice: part 1.
Methods. Ground Improvement, 1: 193–214.
Munfakh, G.A. (2003) Ground improvement in transportation projects: from old visions to
innovative applications. Proceedings of the ICE - Ground Improvement, 7: 47-60.
Najjar, S. (2013) A State-of-the-Art Review of Stone/Sand-Column Reinforced Clay Systems.
Geotechnical and Geological Engineering, 31: 355-386.
Najjar, S., Sadek, S. and Maakaroun, T. (2010) Effect of sand columns on the undrained load
response of soft clays. Journal of Geotechnical and Geoenvironmental Engineering.
Nakase, A. and Kamei, T. (1986) Influence of strain rate on undrained shear characteristics of
K0-consolidated cohesive soils. Japanese Socity of Soil Mechanics and Foundation
Engineering, 26: (1): 85-95.
Ng, C., Zhou, C., Yuan, Q., et al. (2013) Resilient modulus of unsaturated subgrade soil:
experimental and theoretical investigations. Canadian Geotechnical Journal, 50: (2): 223-
232.
Peacock, W. and Bolton, S. (1968) Sand liquefaction under cyclic loading simple shear
conditions. Journal of Soil Mechanics & Foundations Div.
Poorooshasb, H.B. and Meyerhof, G.G. (1997) Analysis of behavior of stone columns and
lime columns. Computers and Geotechnics, 20: 47–70.
168
Samir Ashour Reference
169
Samir Ashour Reference
Sivakumar, V., Boyd, J.L., Black, J.A., et al. (2010) Effects of granular columns in compacted
fills. Proceedings of the ICE - Geotechnical Engineering, 163: 189-196.
Sivakumar, V., Glynn, D., Black, J., et al. (2007) "A laboratory model study of the
performance of vibrated stone columns in soft clay". Proceedings of the 14th European
Conference on Soil Mechanics and Geotechnical Engineering, Madrid, Spain.
Sivakumar, V., McKelvey, D., Graham, J., et al. (2004b) Triaxial tests on model sand
columns in clay. Canadian Geotechnical Journal, 41: 299-312.
Skempton, A. (1957) Discussion: Further data on the c/p ratio in normally consolidated clays.
Proceedings of the Institution of Civil Engineers, 7: 305-307.
Slocombe, B., Bell, A. and Baez, J. (2000) The densification of granular soils using vibro
methods. Geotechnique, 50: (6): 11.
Sun, X., Han, J., Wayne, M.H., et al. (2014) "Quantifying the Benefit of Triaxial Geogrid in
Stabilizing Granular Bases over Soft Subgrade under Cyclic Loading at Different Intensities".
Geo-Congress 2014 Technical Papers at Geo-characterization and Modeling for
Sustainability. ASCE.
Unnikrishnan, N., Rajagopal, K. and Krishnaswamy, N.R. (2002) Behaviour of reinforced
clay under monotonic and cyclic loading. Geotextiles and Geomembranes, 20: 117–133.
Valls-Marquez, M. (2009) Evaluating the capabilities of some constitutive models in
reproducing the experimental behaviour of stiff clay subjected to tunnelling stress paths.
The University of Birmingham.
Van Impe, W.F., De Cock, F., Van der Cruyssen, J.P., et al. (1997) Soil improvement
experiences in Belgium: part II. Vibrocompaction and stone columns. Proceedings of the
ICE-Ground Improvement, 1: 157–168.
Vardanega, P. and Haigh, S. (2014) The undrained strength–liquidity index relationship.
Canadian Geotechnical Journal, 51: (9): 1073-1086.
Vekli, M., Aytekin, M., Ikizler, S.B., et al. (2012) Experimental and numerical investigation
of slope stabilization by stone columns. Natural hazards, 64: (1): 797-820.
Wang, J.-h., Liu, Y.-f. and Xing, Y. (1998) Estimation of undrained bearing capacity for
offshore soft foundations with cyclic load. China Ocean Engineering, 12: (2): 213-222.
Wang, J., Li, C. and Shi, M. (2006) Experimental study on cyclic mechanics characteristic of
saturated soft clay strata. Trans. Tianjin Univ., 12: (2).
Wang, L.-Z., Li, L.-L. and Dan, H.-B. (2011) Undrained behavior of natural marine clay
under cyclic loading. Ocean Engineering, 38: 1792-1805.
Watts, F. (2000) Specifying vibro stone columns. British Research Establishment (BRE)
Report, Watford, UK.
Wood, D.M., Hu, W. and Nash, D.F.T. (2000) Group effects in stone column foundations:
model tests. Geotechnique, 50: 689–698.
Woodward, J. (2004) An Introduction to Geotechnical Processes. CRC Press.
Yasuhara, Hirao, K. and Hyde, A.F.L. (1992) Effects of cyclic loading on undrained strength
and compressibility of clay. Soils and Foundations, 32: (1): 100-116.
170
Samir Ashour Reference
Yasuhara, Yamanouchi, T. and Hirao, K. (1982) Cyclic strength and deformation of normally
consolidated clay. Soils and Foundations.
Yoo, T. and Selig, E. (1979) Field observations of ballast and subgrade deformations in track.
Transportation Research Record, (733).
Zahmatkesh, A. and Choobbasti, A.J. (2010) Settlement evaluation of soft clay reinforced
with stone columns using the equivalent secant modulus. Arabian Journal of Geosciences.
Zhou, J. and Gong, X. (2001) Strain degradation of saturated clay under cyclic loading.
Canadian Geotechnical Journal, 38: 208-212.
171
Samir Ashour Appendix A
APPENDIX A
Material Properties
100.00
90.00 Test (1)
Percentage of passing (%)
172
Samir Ashour Appendix A
100
90
80
Percentage passing, %
70
60
50
40
30
20 Hydrometer test
10 Supplied by manufacturer WBB
0
0.00 0.01 0.10
Particle size, mm
173
Samir Ashour Appendix A
Soil Index
Liquid Limit
dial gage reading 14.00 13.80 14.10 19.40 19.00 18.90 19.90 20.50 20.30 23.20 22.90 23.50
Average
13.97 19.10 20.23 23.20
Penetration
Container # p6 pla3 y11 LLL m45 120 E 342y w17 8A mc51 mc8
Mass of container
3.15 3.99 3.19 7.84 5.49 4.41 3.34 3.39 4.72 4.73 4.76 5.40
g
Mass of container +
16.31 19.24 23.67 21.59 26.45 20.52 17.11 18.34 25.51 20.64 28.46 26.06
Wet Soil (g)
Mass of container +
11.75 13.95 16.56 16.70 19.00 14.79 12.17 12.98 18.06 14.83 19.80 18.51
Dry Soil (g)
Mass of Moisture
4.56 5.29 7.11 4.89 7.45 5.73 4.94 5.36 7.45 5.81 8.66 7.55
(g)
Moisture content
53.02 53.11 53.18 55.19 55.14 55.20 55.95 55.89 55.85 57.52 57.58 57.59
(%)
Ave. moisture
53.10 55.18 55.89 57.56
content (%)
174
Samir Ashour Appendix A
25.00
20.00
15.00
10.00
50.00 52.00 54.00 56.00 58.00 60.00
Water content, (%)
Plastic Limit
175
Samir Ashour Appendix A
Moisture content
Container No.: f28t31 P12 66 S10 B22-8 212ms2 m9t2 T25 M7T
Wet soil & container, g 8.43 8.64 7.76 14.67 12.19 21.01 11.97 17.55 13.51
Dry soil & container, g 7.41 7.72 6.98 12.20 10.66 17.33 9.84 14.43 11.41
Container,( g) 3.17 3.92 3.67 3.39 5.32 4.22 3.36 4.70 4.78
Dry soil , (g) 4.24 3.80 3.31 8.81 5.34 13.11 6.48 9.73 6.63
Moisture loss, (g) 1.02 0.92 0.78 2.47 1.53 3.68 2.13 3.12 2.10
Moisture content,( %) 24.06 24.21 23.56 28.04 28.65 28.07 32.87 32.07 31.67
Average moisture, (%) 23.94 28.25 32.20
Dry density, (Mg/m3) - 1.41 1.38
Test NO.: 4 5 6
Cylinder & soil, A g 6738 6677 6593
Cylinder, B g 4978.5 4978.5 4979.3
Wet soil, A-B g 1759.5 1698.5 1613.7
Volume, cm3 0.998 0.998 0.998
Wet density, ρ Mg/m3 1.76 1.70 1.62
Vane shear strength 23 14 -
176
Samir Ashour Appendix A
Moisture content
Container No.: f28t 39s CL17 SA01 SA02 SA03 m90u f28t31 39s
Wet soil & container, g 18.65 8.65 21.11 12.59 12.14 18.65 19.50 11.17 13.18
Dry soil & container, g 14.96 7.19 17.22 9.93 9.65 14.22 14.46 8.33 9.64
Container, g 4.77 3.18 6.54 3.29 3.50 3.19 5.32 3.18 3.18
Dry soil , g 10.19 4.01 10.68 6.64 6.15 11.03 9.14 5.15 6.46
Moisture loss, g 3.69 1.46 3.89 2.66 2.49 4.43 5.04 2.84 3.54
Moisture content, % 36.21 36.41 36.42 40.06 40.49 40.16 55.14 55.15 54.80
Average moisture, % 36.35 40.24 55.03
Dry density ρd Mg/m3 1.29 1.21 1.04
1.90
1.80
Dry Density (Mg/m3)
1.70
1.60 5%
10%
1.50 0%
Dry Density
1.40 Wet Density
1.30
1.20
1.10
1.00
10 20 30 40 50
Moisture Content %
A B C Average A B C Average
177
Samir Ashour Appendix A
80
50
40
30
20
10
0
30.00 35.00 40.00 45.00 50.00
Water content; %
178
Samir Ashour Appendix A
250
Wet condition 200kPa
200 Test 1_200kPa
Test 2_200kPa
Shear stress, (kPa)
Test 1_100kPa
150 100kPa
Test 2_100kPa
Test 1_50kPa
100 Test 2_50kPa
50kPa Test 1_25kPa
Test 2_25kPa
50
25kPa
0
0 2 4 6 8 10
Horisontal displacement, (mm)
179
Samir Ashour Appendix A
300
Dry condition
250 200 kPa
Test 1_ 25 kPa
200
Shear stress, (kPa)
300
Wet condition
250
Dry condition
200
Shear stress; kPa
150
100
50
0
0 50 100 150 200 250
Normal stress, (kPa)
180
Samir Ashour Appendix A
Test No.: 1 2 3 4
Mass of bottle (g) m1 48.16 46.88 49.48 48.57
Mass of bottle & soil (g) m2 49.97 49.51 51.59 50.74
Mass of bottle water & soil (g) m3 104.05 108.38 108.49 105.51
Mass of bottle full water (g) m4 102.93 106.75 107.18 104.17
Specific gravity 2.62 2.64 2.64 2.61
average 2.63
Test No.: 1 2 3 4
Mass of bottle (g) m1 48.16 46.88 49.48 48.57
Mass of bottle & soil (g) m2 60.98 61.03 59.62 58.94
Mass of bottle water & soil (g) m3 110.97 115.72 113.60 110.67
Mass of bottle full water (g) m4 102.90 106.78 107.19 104.15
Specific gravity 2.70 2.72 2.72 2.69
average 2.71
181
Samir Ashour Appendix A
182
Samir Ashour Appendix A
Test (1)
Test (2)
20.00
Test (3)
Settlement (mm)
40.00
60.00
183
Samir Ashour Appendix A
Test (1)
Test (2)
Test (3)
Settlement (mm)
10.00
20.00
5.00
10.00
184
Samir Ashour Appendix A
10.00
15.00
10.00
15.00
185
Samir Ashour Appendix B
APPENDIX B
1. Test Procedure
186
Samir Ashour Appendix B
187
Samir Ashour Appendix B
Sample consolidation
188
Samir Ashour Appendix B
189
Samir Ashour Appendix C
APPENDIX C
stress Change
Load Pk Stress corr. Disp. Strain PWP PWP σ3
Posn Pk (kN) (kPa) (kPa) (mm) (%) (kPa) (kPa) (kPa)
-23.11 0.00 0.00 0.00 0.00 0.00 294.77 0 400.58
-23.05 0.04 4.46 4.45 0.38 0.18 295.15 0.54 400.67
-22.99 0.03 3.69 3.68 0.45 0.22 295.32 0.72 400.66
-22.92 0.02 2.67 2.67 0.52 0.25 295.54 0.93 400.69
-22.86 0.02 3.06 3.05 0.58 0.28 296.03 1.42 400.64
-22.79 0.09 10.95 10.92 0.65 0.31 297.80 3.20 400.77
-22.73 0.10 12.73 12.69 0.71 0.34 300.62 6.02 400.79
-22.67 0.16 19.74 19.66 0.77 0.37 302.96 8.36 400.53
-22.60 0.19 24.70 24.60 0.84 0.40 304.98 10.38 400.60
-22.54 0.19 24.70 24.59 0.91 0.44 306.82 12.22 400.56
-22.47 0.21 26.87 26.74 0.97 0.47 308.27 13.67 400.61
-22.40 0.24 31.07 30.91 1.04 0.50 309.80 15.20 400.85
-22.34 0.22 27.88 27.74 1.10 0.53 311.01 16.41 400.84
-22.27 0.25 32.21 32.03 1.17 0.56 312.15 17.54 400.49
-22.21 0.24 31.07 30.88 1.23 0.59 313.54 18.94 400.26
-22.14 0.29 37.31 37.07 1.30 0.62 314.67 20.07 400.63
-22.08 0.28 36.16 35.92 1.36 0.66 315.88 21.27 400.40
-22.01 0.32 40.74 40.46 1.43 0.69 316.77 22.17 400.81
-21.95 0.31 39.22 38.93 1.49 0.72 317.77 23.17 400.63
-21.88 0.30 38.45 38.16 1.56 0.75 318.45 23.84 400.57
-21.82 0.30 37.94 37.65 1.63 0.78 318.70 24.10 400.52
-21.75 0.34 43.29 42.94 1.69 0.81 319.69 25.09 400.40
-21.69 0.36 46.09 45.70 1.76 0.84 320.74 26.14 400.65
-21.62 0.38 47.75 47.33 1.82 0.87 321.64 27.03 400.67
-21.56 0.36 45.84 45.42 1.89 0.91 322.10 27.50 400.66
-21.49 0.41 52.71 52.22 1.95 0.94 322.87 28.27 400.82
190
Samir Ashour Appendix C
191
Samir Ashour Appendix C
192
Samir Ashour Appendix C
193
Samir Ashour Appendix C
194
Samir Ashour Appendix C
195
Samir Ashour Appendix C
196
Samir Ashour Appendix C
197
Samir Ashour Appendix C
198
Samir Ashour Appendix C
199
Samir Ashour Appendix C
200
Samir Ashour Appendix C
201
Samir Ashour Appendix C
202
Samir Ashour Appendix C
203
Samir Ashour Appendix C
100
90
80
Deviator stress, (kPa)
70
60
50
(a) Stone column test results
40
30 Test 5-a (clay/column_0.03mm/min)
Test 5-b(clay/column_0.03mm/min)
20
Test 5-c (clay/column_0.03mm/min)
10 Test 5-d (clay/column_0.03mm/min)
0
0 2 4 6 8 10 12 14 16 18
Axial strain, (%)
70.00
60.00
Deviator stress, q, (kPa)
50.00
40.00
(b) Soil only test results
30.00
Test 1_a (soil only_0.03mm/min)
20.00 Test 1_b (Soil only_0.03mm/min)
Test 1_c (Soil only_0.03mm/min)
10.00
Test 1_d (Soil only_0.03mm/min)
0.00
0.00 2.00 4.00 6.00 8.00 10.00 12.00 14.00 16.00
Vertical strain, (%)
204
Samir Ashour Appendix C
100
80
Deviator stress, q, (kPa)
60
40
Reinforced soil (current study)
0
0 5 10 15
Vertical Strain, (%)
205
Samir Ashour Appendix C
206
Samir Ashour Appendix C
200 0.56 0.19 71.73 24.82 0.37 -19.95 -20.18 2.53 2.29 0.237 1547.82 1.18 1.07 352.07 44.07 354.47 46.47
210 0.57 0.20 73.06 25.59 0.37 -19.87 -20.11 2.60 2.36 0.235 1573.00 1.21 1.10 352.31 44.31 355.00 47.00
220 0.57 0.20 73.18 25.85 0.37 -19.83 -20.07 2.64 2.40 0.258 1581.87 1.23 1.12 352.34 44.34 354.85 46.85
230 0.58 0.20 73.30 25.51 0.38 -19.77 -20.03 2.70 2.45 0.246 1454.92 1.26 1.14 352.79 44.79 355.06 47.06
240 0.56 0.20 71.85 26.02 0.36 -19.72 -19.97 2.75 2.51 0.236 1463.21 1.29 1.17 352.94 44.94 355.38 47.38
250 0.58 0.20 73.30 25.76 0.37 -19.67 -19.91 2.80 2.56 0.239 1582.03 1.31 1.20 353.17 45.17 355.62 47.62
260 0.58 0.20 73.42 26.02 0.37 -19.64 -19.87 2.84 2.60 0.241 1557.74 1.33 1.21 353.56 45.56 355.87 47.87
270 0.58 0.21 73.54 26.27 0.37 -19.59 -19.83 2.89 2.64 0.227 1540.41 1.35 1.24 353.75 45.75 355.92 47.92
280 0.57 0.21 73.06 26.36 0.37 -19.55 -19.78 2.92 2.69 0.247 1615.73 1.36 1.26 353.61 45.61 356.14 48.14
290 0.57 0.20 72.09 26.02 0.36 -19.51 -19.76 2.96 2.72 0.223 1464.98 1.38 1.27 353.85 45.85 356.14 48.14
300 0.57 0.20 72.70 26.02 0.37 -19.52 -19.74 2.96 2.73 0.24 1643.95 1.38 1.28 353.92 45.92 356.15 48.15
310 0.58 0.21 73.54 26.53 0.37 -19.47 -19.71 3.00 2.76 0.247 1538.46 1.40 1.29 354.04 46.04 356.55 48.55
320 0.58 0.21 73.42 26.10 0.37 -19.44 -19.69 3.03 2.78 0.227 1504.57 1.42 1.30 353.99 45.99 356.42 48.42
330 0.58 0.20 74.03 26.02 0.38 -19.41 -19.64 3.06 2.84 0.247 1661.01 1.43 1.33 354.28 46.28 356.56 48.56
340 0.57 0.21 73.06 26.45 0.37 -19.39 -19.64 3.08 2.84 0.223 1482.19 1.44 1.32 354.29 46.29 356.64 48.64
350 0.58 0.21 73.42 26.19 0.37 -19.37 -19.60 3.10 2.88 0.237 1663.50 1.45 1.34 354.69 46.69 356.79 48.79
360 0.58 0.20 73.42 26.02 0.37 -19.34 -19.58 3.13 2.89 0.226 1570.89 1.46 1.35 354.57 46.57 357.07 49.07
370 0.58 0.21 73.42 26.19 0.37 -19.32 -19.54 3.16 2.93 0.255 1641.42 1.48 1.37 354.77 46.77 357.04 49.04
380 0.57 0.21 73.06 26.62 0.36 -19.28 -19.54 3.19 2.93 0.236 1430.43 1.49 1.37 354.77 46.77 357.13 49.13
390 0.58 0.21 73.91 26.62 0.37 -19.27 -19.51 3.20 2.97 0.235 1573.77 1.50 1.39 354.66 46.66 356.98 48.98
400 0.58 0.21 74.03 26.79 0.37 -19.25 -19.49 3.22 2.98 0.249 1578.81 1.50 1.39 354.61 46.61 357.19 49.19
410 0.57 0.21 73.06 26.27 0.37 -19.22 -19.47 3.25 3.01 0.258 1475.66 1.52 1.40 354.65 46.65 357.16 49.16
420 0.58 0.20 73.66 25.85 0.38 -19.19 -19.45 3.28 3.02 0.226 1455.58 1.53 1.41 354.73 46.73 357.27 49.27
430 0.58 0.21 73.78 26.27 0.37 -19.20 -19.43 3.27 3.04 0.248 1651.06 1.53 1.42 354.98 46.98 357.25 49.25
440 0.58 0.21 74.03 26.27 0.38 -19.17 -19.41 3.31 3.06 0.24 1512.26 1.55 1.43 355.18 47.18 357.28 49.28
450 0.57 0.21 72.94 26.10 0.37 -19.15 -19.39 3.32 3.08 0.236 1532.63 1.55 1.44 355.21 47.21 357.56 49.56
460 0.58 0.21 73.42 26.53 0.37 -19.14 -19.38 3.33 3.10 0.233 1560.51 1.56 1.45 355.36 47.36 357.57 49.57
207
Samir Ashour Appendix C
470 0.58 0.21 73.66 26.36 0.37 -19.11 -19.35 3.36 3.13 0.245 1594.51 1.57 1.46 355.26 47.26 357.70 49.70
480 0.58 0.20 73.91 26.02 0.38 -19.10 -19.34 3.37 3.13 0.241 1535.10 1.58 1.46 355.21 47.21 357.69 49.69
490 0.57 0.21 73.18 26.10 0.37 -19.08 -19.32 3.39 3.15 0.23 1534.15 1.59 1.47 355.32 47.32 357.61 49.61
500 0.58 0.20 73.91 25.93 0.38 -19.07 -19.30 3.41 3.18 0.226 1638.13 1.59 1.48 355.32 47.32 357.64 49.64
510 0.58 0.21 73.42 26.53 0.37 -19.07 -19.29 3.41 3.18 0.249 1629.56 1.59 1.49 355.35 47.35 357.66 49.66
520 0.58 0.21 73.66 26.10 0.37 -19.04 -19.29 3.43 3.19 0.229 1500.12 1.60 1.49 355.50 47.50 357.82 49.82
530 0.58 0.21 73.30 26.36 0.37 -19.03 -19.26 3.45 3.22 0.259 1609.91 1.61 1.50 355.75 47.75 357.89 49.89
540 0.58 0.21 73.66 26.27 0.37 -19.01 -19.27 3.46 3.20 0.231 1437.03 1.62 1.50 355.71 47.71 357.85 49.85
550 0.58 0.21 73.30 26.45 0.37 -19.01 -19.24 3.46 3.23 0.25 1593.07 1.62 1.51 355.57 47.57 357.97 49.97
560 0.58 0.20 73.54 25.93 0.37 -18.99 -19.24 3.49 3.24 0.235 1495.68 1.63 1.51 355.69 47.69 358.12 50.12
570 0.58 0.21 73.78 26.27 0.37 -18.98 -19.21 3.49 3.26 0.234 1587.83 1.63 1.52 355.65 47.65 357.90 49.90
580 0.58 0.21 73.54 26.45 0.37 -18.98 -19.21 3.50 3.26 0.238 1580.77 1.63 1.53 355.59 47.59 357.92 49.92
590 0.58 0.21 73.66 26.27 0.37 -18.96 -19.20 3.51 3.27 0.227 1563.82 1.64 1.53 355.71 47.71 358.08 50.08
600 0.58 0.21 73.91 26.10 0.38 -18.95 -19.18 3.52 3.29 0.234 1653.88 1.64 1.54 355.70 47.70 357.87 49.87
610 0.57 0.21 72.94 26.53 0.36 -18.95 -19.18 3.53 3.29 0.236 1557.61 1.65 1.54 355.83 47.83 357.95 49.95
620 0.57 0.21 73.18 26.36 0.37 -18.92 -19.16 3.55 3.31 0.237 1558.14 1.66 1.55 355.82 47.82 358.05 50.05
630 0.58 0.21 74.03 26.53 0.37 -18.92 -19.16 3.55 3.31 0.247 1573.97 1.66 1.55 355.98 47.98 358.16 50.16
640 0.58 0.21 73.66 26.27 0.37 -18.91 -19.16 3.56 3.31 0.237 1506.84 1.66 1.55 355.86 47.86 358.24 50.24
650 0.58 0.21 74.03 26.27 0.38 -18.91 -19.15 3.56 3.33 0.221 1582.45 1.66 1.55 355.99 47.99 358.32 50.32
660 0.58 0.21 73.30 26.19 0.37 -18.89 -19.12 3.58 3.36 0.24 1674.25 1.67 1.57 355.94 47.94 358.47 50.47
670 0.59 0.21 74.51 26.62 0.38 -18.88 -19.12 3.59 3.35 0.235 1567.33 1.68 1.56 355.87 47.87 358.21 50.21
680 0.58 0.21 74.39 26.96 0.37 -18.89 -19.12 3.59 3.35 0.241 1585.23 1.68 1.57 355.96 47.96 358.28 50.28
690 0.58 0.21 73.91 26.70 0.37 -18.87 -19.11 3.61 3.36 0.242 1538.34 1.68 1.57 356.02 48.02 358.16 50.16
700 0.58 0.21 73.91 27.04 0.37 -18.86 -19.10 3.61 3.37 0.247 1520.91 1.69 1.57 356.02 48.02 358.39 50.39
710 0.58 0.21 73.54 26.27 0.37 -18.85 -19.09 3.63 3.38 0.237 1503.00 1.70 1.58 356.08 48.08 358.41 50.41
720 0.59 0.21 74.75 26.45 0.38 -18.84 -19.08 3.63 3.40 0.243 1600.84 1.70 1.59 356.04 48.04 358.56 50.56
730 0.58 0.21 74.15 26.27 0.38 -18.84 -19.08 3.64 3.40 0.225 1547.28 1.70 1.59 356.05 48.05 358.40 50.40
208
Samir Ashour Appendix C
740 0.58 0.21 74.15 26.27 0.38 -18.83 -19.06 3.64 3.41 0.218 1671.07 1.70 1.60 356.14 48.14 358.47 50.47
750 0.58 0.21 74.15 26.19 0.38 -18.83 -19.05 3.64 3.43 0.255 1727.80 1.70 1.60 356.19 48.19 358.62 50.62
760 0.58 0.21 73.54 26.53 0.37 -18.81 -19.07 3.66 3.41 0.236 1447.96 1.71 1.59 356.23 48.23 358.41 50.41
770 0.58 0.21 73.78 26.19 0.37 -18.82 -19.06 3.65 3.42 0.232 1583.94 1.71 1.60 356.09 48.09 358.35 50.35
780 0.58 0.21 73.42 26.36 0.37 -18.80 -19.03 3.67 3.44 0.237 1593.19 1.72 1.61 356.08 48.08 358.40 50.40
790 0.57 0.21 73.18 26.19 0.37 -18.80 -19.03 3.68 3.44 0.229 1557.22 1.72 1.61 356.07 48.07 358.48 50.48
800 0.58 0.21 73.42 26.19 0.37 -18.79 -19.02 3.68 3.45 0.245 1619.91 1.72 1.61 356.19 48.19 358.55 50.55
810 0.57 0.21 73.18 26.79 0.36 -18.79 -19.03 3.69 3.44 0.245 1487.22 1.72 1.61 356.23 48.23 358.82 50.82
820 0.58 0.21 74.15 26.87 0.37 -18.77 -19.02 3.70 3.46 0.247 1515.51 1.73 1.61 356.31 48.31 358.96 50.96
830 0.58 0.21 74.15 26.19 0.38 -18.76 -19.01 3.71 3.46 0.234 1524.94 1.73 1.62 356.34 48.34 358.67 50.67
840 0.58 0.21 73.91 27.04 0.37 -18.76 -19.00 3.71 3.48 0.23 1572.91 1.73 1.62 356.32 48.32 358.68 50.68
850 0.58 0.21 73.78 26.19 0.37 -18.76 -18.99 3.71 3.48 0.247 1625.26 1.73 1.63 356.41 48.41 358.66 50.66
860 0.58 0.21 73.54 26.45 0.37 -18.75 -19.00 3.72 3.48 0.223 1497.57 1.74 1.62 356.30 48.30 358.61 50.61
870 0.58 0.21 73.66 27.04 0.37 -18.75 -18.98 3.72 3.50 0.242 1641.97 1.74 1.63 356.23 48.23 358.56 50.56
880 0.57 0.21 73.06 26.53 0.37 -18.74 -18.98 3.74 3.49 0.255 1510.04 1.75 1.63 356.35 48.35 358.69 50.69
890 0.58 0.21 74.27 26.53 0.37 -18.72 -18.98 3.75 3.50 0.23 1470.31 1.75 1.63 356.31 48.31 358.87 50.87
900 0.58 0.21 73.54 26.79 0.37 -18.73 -18.96 3.74 3.51 0.232 1596.61 1.75 1.64 356.24 48.24 358.65 50.65
910 0.58 0.20 73.66 26.02 0.37 -18.72 -18.96 3.75 3.52 0.235 1612.93 1.75 1.64 356.37 48.37 358.79 50.79
920 0.58 0.21 74.27 26.36 0.38 -18.73 -18.96 3.75 3.51 0.236 1601.15 1.75 1.64 356.51 48.51 358.87 50.87
930 0.57 0.21 72.94 26.36 0.37 -18.72 -18.95 3.76 3.52 0.238 1550.08 1.76 1.65 356.73 48.73 358.99 50.99
940 0.58 0.21 74.15 26.27 0.38 -18.71 -18.95 3.76 3.52 0.251 1579.79 1.76 1.65 356.37 48.37 358.91 50.91
950 0.58 0.21 73.42 26.10 0.37 -18.70 -18.95 3.77 3.52 0.232 1480.60 1.76 1.65 356.57 48.57 358.70 50.70
960 0.58 0.21 74.03 26.62 0.37 -18.69 -18.92 3.79 3.55 0.246 1605.00 1.77 1.66 356.30 48.30 358.65 50.65
970 0.58 0.21 73.78 26.45 0.37 -18.70 -18.95 3.77 3.53 0.245 1511.38 1.76 1.65 356.43 48.43 358.72 50.72
980 0.58 0.21 74.15 26.36 0.38 -18.69 -18.94 3.78 3.54 0.221 1531.92 1.77 1.65 356.41 48.41 358.76 50.76
990 0.58 0.21 73.78 26.10 0.37 -18.70 -18.92 3.78 3.56 0.231 1694.48 1.76 1.66 356.53 48.53 358.66 50.66
1000 0.58 0.21 74.15 26.62 0.37 -18.69 -18.92 3.79 3.56 0.241 1616.06 1.77 1.66 356.65 48.65 359.02 51.02
209
Samir Ashour Appendix D
APPENDIX D
1.1 Introduction
Regression analysis is a powerful technique that can be used to address various research
questions. In this report, we are going to use it to check how qult levels are affected by cu
and φ.
In particular, the type of regression we are going to use is ”Multiple Linear Regression”.
Linear regression is the process of finding the best-fitting straight line through data points
(this line is sometimes referred to as the regression line). Multiple means we have more
than one input variable (also known as predictor), hence, we are trying to fit a plane or
hyper-plane rather than a line. The input variables in our case are cu, phi and as. Linear
means that we are trying to find a combination of the input variables such that each
variable is multiplied by a coefficient and then we sum the products. The idea is to use this
linear combination of input variables to model their relationship with an output variable (in
for statistical computing and graphics. R makes it easy to manipulate data and perform
210
Samir Ashour Appendix D
R-squared
This value calculates the percentage of variation of the output explained by the input
variables in the model. This means the higher the value of R-squared the better the model.
R-squared adjusted
This value is similar to R-squared but it accounts for the number of input variables in the
Residuals
A residual is the difference between the actual value and predicted value for each point (or
record) in the data. Histograms are often used to check the distribution of residuals. Also,
they are plotted against each input variable. If a model fits well, the residuals will be small
and will be no pattern of their distribution around zero (i.e. they should be evenly spread
around zero).
Deviance
The deviance is a statistic that is used to determine the quality of fit for a model. It is a
measure of how much better a model with more parameters fits the data. It is used to
compare nested models. A nested model is a model which is a subset of another model.
For example, if we have two models, the first describes the relationship between one input
variable x and an output variable y and the second describes the relationship between two
input variables x and z and an output variable y, then the first model is nested within the
second.
211
Samir Ashour Appendix D
When there are a limited number of models, the BIC is often used for model selection.
Often, the model with the smallest BIC is preferred. The BIC is a penalised version of the
P-values
propabilities (i.e. their values always lie between zero to 1) and they show how likely
certain situations are. P-values which are close to zero (usually <= 0.05) are more likely to
occur if the study has shown something positive. It is said that the result is significant if
the p-value is close to zero. On the other hand, the result is said to be non-significant if the
Confidence Intervals
It is known that in statistics we use the sample data at hand to draw inferences about the
entire population (i.e. all the data) and make an estimate of the value(s) we are trying to
precision. This measure of precision depends on the sample size and normally takes the
form of a 95% Confidence Interval or a standard error value (the former is calculated from
the latter). The 95% confidence interval gives the range of poulation parameters that the
sample leads us to believe are possible. The 95% confidence interval is presented as a
range of two values (a, b) and is interpreted as: we can be 95% confident that the
result/effect we are trying to measure what will happen by an overage of at least a and
maybe as much as b
212
Samir Ashour Appendix D
Data base
213
Samir Ashour Appendix D
Table 1: Statistical Summaries of the Variables with 95% Confidence Interval for the
Mean values
Property cu ∅ As qult
Valid 31 31 31 31
Missing - - - -
By examining the table, we observe that we have data for all the points (i.e. no missing
values). Also, we notice that there is no big deference between the mean and median
values of each variable, hence outliers are unlikely (outliers usually have a big influence
on the deference between the mean and the median). Another values that can analyse from
the table are the minimum and maximum values for each input variable. They appear to be
within possible ranges for all of the three input variables. Finally, as standard deviation is
an indicator of how spread out the data is, we can check the validity of our data by going 2
standard deviations on each side of the mean for the outcome variable Q. We notice that
more than 95% of all values of this variable lie within that range.
214
Samir Ashour Appendix D
close look at these histograms indicate some "bunching" around multiples of 5 for the cu
input variable (Figure 1a) and around multiples of 10 for the AS variable (Figure 1c).
plot matrix gives the scatter plots of qult against each of the other three input variables.
Additionally, we show the univariate relationship between our outcome variable qult and
the input variables in Table 2. From this table, we can notice that the correlation between
cu and qult and between φ and qult is significant at the p < 0.05 level.
215
Samir Ashour Appendix D
Pearson Spearman
Information also are provided about the correlations between our input variables in table 3.
From this table we can observe that none of the correlations between the input variables is
216
Samir Ashour Appendix D
As cu ∅
As 1 -0.2596521 0.2605386
(0.1584) (0.1569)
cu -0.2596521 1 -0.2176558
(0.1584) (0.2395)
∅ 0.2605386 -0.2176558 1
(0.1569) (0.2395)
As cu ∅
As 1 -0.3944457 0.204013
(0.0281) (0.271)
cu -0.3944457 1 -0.2161602
(0.0281) (0.2428)
∅ 0.204013 -0.2161602 1
(0.271) (0.2428)
would seem reasonable to model CONFOUNDER FIRST and then enter cu into the model
and see whether it is associated with qult. We can do this in the opposite way by entering cu
into the model first and then follow it by φ and AS in turn to see what impact each of them
217
Samir Ashour Appendix D
In this section, a simple linear regression model is built using just cu as input and qult as
output. After using R's lm() function, the model looks as follows:
Now after building the model that describes the relationship between cu and qult, let us plot
the input variable cu against the residuals. As Figure 4 shows, the residual values do not
seem to have a particular pattern and they are randomly scattered around zero.
218
Samir Ashour Appendix D
In this section, we are going to build a simple linear regression model using just AS as
input and Q as output. After using R's lm() function, the model looks as follows:
219
Samir Ashour Appendix D
Now after building the model that describes the relationship between AS and qult, let us plot
the input variable AS against the residuals. As Figure 6 shows, the residual values do not
seem to have a particular pattern and they are randomly scattered around zero (as we stated
in Section 4.3).
220
Samir Ashour Appendix D
In this section, a simple linear regression model was built using just φ as input and qult as
output. After using R's lm() function, the model looks as follows:
221
Samir Ashour Appendix D
Now after building the model that describes the relationship between φ and qult, let us plot
the input variable PHI against the residuals. As Figure 8 shows, the residual values do not
seem to have a particular pattern and they are randomly scattered around zero.
222
Samir Ashour Appendix D
In this section, we are going to build a simple linear regression model using CU and AS as
inputs and Q as output. After using R's lm() function, the model looks as follows:
by one unit, causes an increase, or decrease, in qult by 10.194 units. Similarly, when fixing
4.174 units.
Now after building the model that describes the relationship between cu, AS and qult, let us
plot the input variables cu and AS against the residuals. As Figure 9 shows, the residual
values do not seem to have a particular pattern and they are randomly scattered around
zero.
223
Samir Ashour Appendix D
In this section, we are going to build a simple linear regression model using cu and φ as
inputs and qult as output. After using R's lm() function, the model looks as follows:
one unit, causes an increase, or decrease, in qult by 10.29 units. Similarly, when fixing CU,
units. Now after building the model that describes the relationship between cu, PHI and
qult, let us plot the input variables cu and φ against the residuals. As Figure 10 shows, the
residual values do not seem to have a particular pattern and they are randomly scattered
around zero.
224
Samir Ashour Appendix D
In this section, we are going to build a simple linear regression model using AS and φ as
inputs and qult as output. After using R's lm() function, the model looks as follows:
225
Samir Ashour Appendix D
by one unit, causes an increase, or decrease, in qult by 1.536 units. Similarly, when fixing
10.8 units. Now after building the model that describes the relationship between AS, φ and
qult, let us plot the input variables AS and φ against the residuals. As Figure 11 shows, the
residual values do not seem to have a particular pattern and they are randomly scattered
around zero.
226
Samir Ashour Appendix D
1.7.7 Modelling the Relationship between cu, AS, φ and qult (Model 7)
In this section, we are going to build a simple linear regression model using CU, AS and
PHI as inputs and Q as output (recall this is the purpose of this study). After using R's lm()
By examining Equation 7, we observe that when fixing φ and AS, an increase, or decrease,
of cu by one unit, causes an increase, or decrease, in qult by 11.313 units. Similarly, when
fixing cu and AS, an increase, or decrease, of φ by one unit, causes an increase, or decrease,
in qult by 14.629 units. Also, when fixing cu and φ, an increase, or decrease, of AS by one
227
Samir Ashour Appendix D
Now after building the model that describes the relationship between cu, AS, φ and qult, let
us plot the input variables cu, AS and φ against the residuals. As Figure 12 shows, the
residual values do not seem to have a particular pattern and they are randomly scattered
around zero.
228
Samir Ashour Appendix D
Five diagnostics that can be used to examine the goodness of fit of a model were
mentioned In Section 1.3. Namely, these were: Deviance, BIC, R-Squared, Adjusted R-
Squared and the Residuals. Table 4 shows the values of the remaining four diagnostics for
Table 4: Deviance, BIC, R-Squared and Adjusted R-Squared of the Seven models
229
Samir Ashour Appendix D
As we mentioned before, lower BIC indicates a better fitting model. By analysing table 4
we observe that Model 5 has the lowest BIC with Model 7 in second. Howerver, when we
examine the value of R-Squared, we realise that Model 7 has the highest R-Squared
amongst all the models (see Section 4.1). This gives us confidence that from amongst the
seven models that we created using various combinations of the input variables, Model 7 is
the best model that describes the relationship between the input variables CU, AS and PHI
230
Samir Ashour Appendix F
APPENDIX F
Iteration Historyb
Residual Parameter
Iteration Sum of
a
Number Squares a b m
1.0 71288.417 .000 .000 .000
1.1 29769.136 1.807 .000 .000
2.0 29769.136 1.807 .000 .000
2.1 23243.121 3.354 .202 3.196
3.0 23243.121 3.354 .202 3.196
3.1 4830.984 3.625 .145 2.993
4.0 4830.984 3.625 .145 2.993
4.1 4333.676 4.377 .091 2.411
5.0 4333.676 4.377 .091 2.411
5.1 4316.993 4.332 .097 2.477
6.0 4316.993 4.332 .097 2.477
6.1 4316.966 4.348 .096 2.468
7.0 4316.966 4.348 .096 2.468
7.1 4316.965 4.346 .096 2.470
8.0 4316.965 4.346 .096 2.470
8.1 4316.965 4.346 .096 2.469
Derivatives are calculated numerically.
a. Major iteration number is displayed to the left of the decimal,
and minor iteration number is to the right of the decimal.
b. Run stopped after 16 model evaluations and 8 derivative
evaluations because the relative reduction between successive
residual sums of squares is at most SSCON = 1.000E-8.
231
Samir Ashour Appendix F
Parameter Estimates
95% Confidence Interval
Std. Lower Upper
Parameter Estimate Error Bound Bound
a 4.346 .197 3.960 4.733
b .096 .005 .085 .107
m 2.469 .170 2.136 2.803
Correlations of Parameter
Estimates
a b m
a 1.000 -.400 .550
b -.400 1.000 .540
m .550 .540 1.000
ANOVAa
Sum of Mean
Source Squares df Squares
Regression 66971.451 3 22323.817
Residual 4316.965 3058 1.412
Uncorrected
71288.417 3061
Total
Corrected Total 4932.510 3060
Dependent variable: VAR00005
a. R squared = 1 - (Residual Sum of Squares) / (Corrected
Sum of Squares) = .125.
232
Samir Ashour Appendix F
Iteration Historyb
Residual Parameter
Iteration Sum of
a
Number Squares a b m
1.0 105048.595 .000 .000 .000
1.1 59039.952 1.285 .000 .000
2.0 59039.952 1.285 .000 .000
2.1 14058.800 3.055 .145 2.207
3.0 14058.800 3.055 .145 2.207
3.1 15227.522 6.360 .098 4.206
3.2 9802.091 4.171 .146 3.403
4.0 9802.091 4.171 .146 3.403
4.1 10647.892 7.604 .105 4.544
4.2 8249.791 5.485 .145 4.075
5.0 8249.791 5.485 .145 4.075
5.1 7672.376 8.405 .108 4.564
6.0 7672.376 8.405 .108 4.564
6.1 6139.380 11.429 .082 4.654
7.0 6139.380 11.429 .082 4.654
7.1 6122.613 17.495 .034 4.895
8.0 6122.613 17.495 .034 4.895
8.1 5059.439 20.596 .030 4.991
9.0 5059.439 20.596 .030 4.991
9.1 5023.065 23.229 .019 5.066
10.0 5023.065 23.229 .019 5.066
10.1 5017.609 23.381 .019 5.082
11.0 5017.609 23.381 .019 5.082
11.1 5017.607 23.393 .019 5.080
12.0 5017.607 23.393 .019 5.080
12.1 5017.607 23.392 .019 5.081
Derivatives are calculated numerically.
a. Major iteration number is displayed to the left of the decimal,
and minor iteration number is to the right of the decimal.
b. Run stopped after 26 model evaluations and 12 derivative
evaluations because the relative reduction between successive
residual sums of squares is at most SSCON = 1.000E-8.
233
Samir Ashour Appendix F
Parameter Estimates
95% Confidence Interval
Std. Lower Upper
Parameter Estimate Error Bound Bound
a 23.392 .445 22.519 24.264
b .019 .002 .015 .024
m 5.081 .039 5.004 5.157
Correlations of Parameter
Estimates
a b m
a 1.000 -.751 .335
b -.751 1.000 .355
m .335 .355 1.000
ANOVAa
Sum of Mean
Source Squares df Squares
Regression 100030.989 3 33343.663
Residual 5017.607 6050 .829
Uncorrected
105048.595 6053
Total
Corrected Total 26624.392 6052
Dependent variable: permanent deformation
a. R squared = 1 - (Residual Sum of Squares) / (Corrected
Sum of Squares) = .812.
234
Samir Ashour Appendix G
APPENDIX G
application of 125 kN wheel load on rail across three sleepers as shown in Figure 1a. The
cross section of the setup is shown in Figure 1b. Load was applied at 2 Hz for 2 million
cycles. Transducer positioned at top of subgrade under the central sleeper (where load was
applied) in order to monitor the stresses. Figure 2 showed that the dynamic stresses at the
subgrade level were approximately 80 kPa showing good agreement with Yoo and Selig
(1979) study.
235
Samir Ashour Appendix G
Rail Pressure
215 mm Load ram Transducer
Rail conneted to three sleepers Sleeper
Pressure Transducer
Half sleeper PT
Ballast 2
600mm
X
300mm
X
1800mm
PT
570mm
Subgrade material
1
600mm
225mm
100mm
Building sand
730mm
PT
3
3000.00mm
Load
application
100
PT1
PT2
Stresses at the subgrade level, (kPa)
PT3
80
60
40
20
0
0 50 100 150
Time, (minute)
236
Samir Ashour Appendix G
Based on the above, cyclic stresses of 50 kPa, 60 kPa and 70 kPa were used in this study.
These values are equivalent to a cyclic stress ratio (CSR) of 0.6, 0.7 and 0.8 (CSR = the
cyclic deviator stress (qcyclic)/ the static deviator stress of reinforced soil at failure (qfailure)).
Surcharge Pressure
Surcharge pressure, equivalent to the weight of the track component (see Figure 3) was
applied as the lowest pressure level during cyclic loading. this could be ranged between 10
and 25 kPa depending on the rail and sleeper types and the depth of the used ballast layer
(Selig and Waters, 1994; Brough et al., 2003). Using ballast thickness of 300 mm with a
density of 19 kN/m3, rail weight of 60 kg/m, and a concrete sleeper of 300 kg with typical
to centre, the surcharge pressure should be equivalent to 9.4 kPa. However, due to
SSleeper
lSleeper
Figure 3 Structure of a rail track system used to calculate the surcharge pressure.
237